Приземный слой

Прежде чем рассматривать спектры турбулентных пульсаций, необходимо дать некоторое описание приземного слоя. (дополнить)

При теоретическом анализе процессов в приземном слое атмосферы мы будем исходить из общепринятой схемы потока над безграничной шероховатой поверхностью, свойства которой предполагаются достаточно однородными по горизонтали. Осредненные характеристики течения в этой схеме зависят только от вертикальной координаты . Наиболее важными характеристиками являются потоки импульса, тепла и влаги.

Поток импульса можно трактовать как напряжение турбулентного трения. Вместо турбулентного трения

, (13)

где и – пульсации горизонтальной и вертикальной компонент ветра, – плотность воздуха, а черта сверху означает осреднение, удобно рассматривать динамическую скорость

(14)

В пределах приземного слоя и турбулентный поток тепла можно считать практически не зависящими от высоты .

Условие постоянства потоков импульса и тепла (в пределах данного допуска) может служить определением самого понятия приземного слоя. Обухов в своей работе (ссылка) дает примерную оценку высоты приземного слоя на основе анализа возможных изменений . Исходить следует из осредненных уравнений гидромеханики в поле силы Кориолиса. Соответствующее уравнение для координаты (направление скорости ветра у Земли) в квазистационарном случае имеет следующий вид:

,(15)

где – градиент давления, – параметр Кориолиса, ­– компонента осредненной скорости ветра по направлению оси .

Проинтегрируем обе части уравнения по высоте в пределах слоя толщиной и оценим правую часть:

(16)

Отбрасывание члена приводит к усилению неравенства, так как сила Кориолиса частично компенсирует силу градиента давления. Вводя динамическую скорость и скорость геострофического ветра , можно записать полученное неравенство в следующей форме:

(17)

Определим так, чтобы относительное изменение в слое толщиной не превосходило допуска , т.е.

(18)

В силу неравенства (17) для выполнения (18) достаточно, чтобы

.(19)

Отношение скорости трения к скорости геострофического ветра можно оценить величиной порядка 0,05:

,

откуда следует, что

.

При и получаем (м):

При допуске получаем искомую оценку высоты приземного слоя: . В пределах этого слоя можно считать практически постоянным и пренебречь действием силы Кориолиса (вращением ветра с высотой). Полученная оценка достаточно хорошо согласуется с данными наблюдений. (ссылки на наблюдения)

(краткое описание видов стратификации)

СЛУЧАЙ БЕЗРАЗЛИЧНОЙ СТРАТИФИКАЦИИ.

Начнем с рассмотрения самого простого (но относительно редкого) случая безразличной (нейтральной) термической стратификации. В этом случае турбулентный поток тепла равен нулю и турбулентный режим приземного слоя определяется единственным размерным параметром – вертикальным турбулентным потоком импульса (т.е. напряжением трения) .

При условиях безразличной стратификации процессы турбулентного перемешивания в приземном слое могут быть описаны схемой логарифмического пограничного слоя. Соответствующие закономерности детально изучены в экспериментальной аэрогидродинамике и широко используются в метеорологии.

Вывод логарифмического закона распределения ветра осуществляется на основе гипотез подобия. Предположим, что для значений , где – высота травостоя (характерный масштаб микронеоднородностей подстилающей поверхности), статистические характеристики для относительных движений в потоке инвариантны по отношению к преобразованиям подобия , , , . При этих преобразованиях полупространство переходит само в себя, а уравнения движения остаются неизменными. Это обстоятельство является теоретическим основанием для принятой гипотезы подобия. Заметим также, что естественный масштаб скорости остается инвариантным по отношению к указанным преобразованиям. Рассмотрим стационарный режим и составим отношение разности осредненных скоростей на двух уровнях , к динамической скорости . Соответствующая безразмерная величина является функцией и и в силу предположения о самоподобии потока может зависеть только от отношения :

.(20)

Определим вид функции . Очевидно, что для трех высот

(21)

и вместе с тем

.(22)

Отсюда следует, что функция удовлетворяет функциональному уравнению

, (23)

, .

Естественным решением этого функционального уравнения является логарифмическая функция . Полагая , получаем:

, (24)

где – известная постоянная Кармана. Эта формула впервые была выведена еще около 1930 г. Л.Прандлем и Т. Карманом. По эмпирическим данным .

РАСПИСАТЬ ТОЧНОСТЬ ЗНАЧЕНИЙ æ

Уравнение (24) можно записать в обычной дифференциальной форме, рассматривая бесконечно близкие значения и :

.(25)

Уравнения (24) и (25) не содержат характеристик подстилающей поверхности и могут быть отнесены к любой подстилающей поверхности, если выполняется условие . Вместе с тем формула (25) определяет только изменения средней скорости ветра с высотой. Определение абсолютного значения требует учета свойств подстилающей поверхности.

Предположим теперь, что наблюдения за скоростью ветра проводятся на фиксированной высоте над некоторой определенной подстилающей поверхностью. Допустим, что мы располагаем возможностью проводить независимые измерения величины турбулентного трения и, следовательно, можем в каждом отдельном случае определять . Значение может быть определено, например, из термоанемометрических наблюдений за пульсациями и или суммарно на основании измерения напряжения трения у поверхности Земли. Последний метод практически используется при изучении турбулентного движения в трубах. В условиях атмосферы попытка применения динамического метода измерения была сделана Шеппардом(ссылка).

Сравнение рядов наблюдений за и позволит нам определить зависимость между этими величинами. Опыт аэродинамики учит, что при больших числах Рейнолдса и «шероховатой» поверхности зависимость от имеет квадратичный характер, и, следовательно,

, (26)

где – безразмерный коэффициент, зависящий от свойств подстилающей поверхности. При фиксированной высоте «коэффициент трения» может служить объективной характеристикой свойств подстилающей поверхности, описывающей ее динамическое воздействие на поток. Однако недостаток использования заключается в том, что мы связаны с определенным выбором высоты наблюдения. Зависимость от высоты наблюдения легко установить, подставив в формулу(26). Для любых двух высот мы будем иметь

. (27)

Из уравнения (27) следует, в частности, что убывает с высотой. Потенциируя и объединяя величины, содержащие соответственно и , получаем:

,(28)

т.е. величину, не зависящую от высоты. Таким образом, величина , имеющая размерность длины, определяется только свойствами подстилающей поверхности; она носит название «динамической шероховатости». Выразим коэффициент трения через :

. (29)

Тогда на основании формулы(26) можно получить искомое распределение скорости ветра:

. (30)

Преимущество приведенного выше способа введения понятия шероховатости подстилающей поверхности заключается в том, что он опирается исключительно на свойства потока на достаточно больших высотах, где с достаточным основанием можно пользоваться универсальными законами развитой турбулентности. В большинстве случаев, однако, мы не располагаем средствами для прямых измерений и, следовательно, , и в связи с этим при практическом определении динамической шероховатости приходится использовать характеристики профиля ветра, определяемого непосредственно из наблюдений.

СЛУЧАЙ ТЕРМИЧЕСКИ СТРАТИФИЦИРОВАННОГО ПРИЗЕМНОГО СЛОЯ.

Одной из важнейших характеристик турбулентного режима в приземном слое атмосферы является вертикальный турбулентный поток тепла

, (31)

где – удельная теплоемкость воздуха при постоянном давлении, – плотность, и – пульсации вертикальной составляющей скорости ветра и температуры соответственно, обусловленные прохождением через данную точку турбулентных «элементов», а черта сверху означает осреднение. Величина представляет собой среднее количество тепла, переносимое турбулентными пульсациями через единичную площадку в единицу времени. С достаточным основанием можно считать, что в приземном слое при стационарных условиях турбулентный поток тепла практически не зависит от высоты (если не рассматривать потоки лучистой энергии). Вместо величины можно пользоваться также «потоком температуры»

. (32)

Величина турбулентного потока тепла может быть определена экспериментально на основе безынерционных измерений пульсаций температуры и вертикальной составляющей скорости ветра . Современная техника открывает принципиальную возможность таких измерений. Вместе с тем на практике все же приходится пользоваться косвенными методами определения , опирающимися на более простые градиентные измерения. Для правильной интерпретации этих измерений необходимо исследовать связь характеристик турбулентности и с распределениями средней скорости ветра и температуры. При решении этой задачи воспользуемся методами подобия и постараемся установить систему минимального числа параметров, описывающих турбулентный режим в температурно-неоднородной среде.

Неоднородности турбулентного поля, носящие систематический характер (изменение температуры с высотой), оказывают определенное влияние на общий режим турбулентности (действие архимедовых сил). При условии, что температурные пульсации малы по сравнению со средней температурой слоя , уравнения динамики температурно неоднородной среды могут быть записаны в следующей форме:

, (33)

,

,

,

В этой системе и - отклонения от стандартных значений. Упрощения, сделанные при выводе системы уравнений, сводятся к пренебрежению силой Кориолиса и радиационным притоком тепла, а также к линеаризации относительно стандартного статистического распределения давления и температуры. Последнее означает, что пренебрегается изменениями плотности за счет изменений давления, и допускается, что отклонения плотности и температуры от стандартных значений пропорциональны.(ссылка на Ланд-Лиф 14гл 5). Эти упрощения, принятые в теории конвекции, позволяют описать архимедову силу слагаемым . Таким образом, в уравнения входит размерная константа , которую в дальнейшем должны учесть при составлении критериев подобия.

Заметим, что линеаризация уравнений относительно колебаний скорости недопустима, так как при этом была бы потеряна турбулентность. Кроме того, в уравнениях опущены члены, содержащие вязкость и теплопроводность. В условиях развитого турбулентного режима эти члены следует учитывать только при исследовании очень тонких деталей микроструктуры ветрового и температурного поля. Вертикальный перенос импульса и тепла обусловлен неоднородностями некоторого «среднего масштаба», для которых непосредственное влияние вязкости и теплопроводности достаточно мало. Естественно предположить, что изменения средней скорости и температуры с высотой могут быть выражены через координату , параметр и «внешние параметры» и , причем соответствующие уравнения допускают запись в безразмерной форме и не содержат других размерных констант. Принятая нами гипотеза подобия находится в согласии с уравнениями (33) и эквивалентна предположению, что система уравнений (33) вместе с уравнениями

, (34)

, являющимися аналогом граничных условий, однозначно определяют статистические характеристики турбулентного режима. Таким образом, три параметра: , и –можно считать исчерпывающими характеристиками турбулентности приземного слоя атмосферы. Из этих параметров можно однозначно (с точностью до числового множителя) построить масштаб длины и масштаб температуры , которые будут иметь следующий вид:

, (35) (масштаб обухова)