книги из ГПНТБ / Мельникова И.И. Динамическая метеорология учеб. пособие для океанологов
.pdf«Л- Я - К = 4 л + ^ |
arc tg у + ІП |
- К ; (0.5.5) |
|
К = у і \ |
(5.5.6) |
|
1 —V4 |
(5.5.7) |
|
In- |
|
|
|
|
В выписанных выше |
соотношениях у — вспомогательная |
функция, связанная с г лформулой (5.5.3): при устойчивой стра
тификации (L>0) у< 1, при |
неустойчивой |
(Е<0) у > 1, |
при |
||
гриближении к безразличной |
стратификации |
оо) |
у -*• 1; |
||
й л= |
I ат• d Ѳп и <:/„= I OLq-dqn— приведенные |
значения |
темпера |
||
туры |
и ѵдедьноп влажности, |
учитывающие |
зависимость |
а т и |
zq от стратификации z„=z/L. Современные на_блюдения показы вают, что а т= а 9. На рис. 35 приводится зависимость « т от гп.
«г
Рис. 35’. Зависимость от стратификации обратно го турбулентного числа Праидтля. 1— РаупдХилл; 2 — Австралия
Если воспользоваться зависимостью а т = а г (*„), то после некоторых преобразований и численного интегрирования мож
но получить соотношения, связывающие Ѳл |
с Ѳл (аналогично |
|
<7„c q п), и установить |
зависимость Ѳл и q n |
от у (или z n). Таб |
лицы для определения |
ип, Ѳ„ (или q„) как |
функции z„ или у |
приводятся в [11]. Для перехода от безразмерных величин к раз мерным нужно найти величины турбулентных потоков, от кото рых зависят масштабы в (5.5.1).
100
§ 6. Определение турбулентных потоков на основании градиент ных или стандартных измерений
Рассмотренная в предыдущем параграфе нелинейная модель приземного слоя позволяет определить на основании градиент ных наблюдений L, v#, Po и Е.
Ф(і) iffФЩ
Рис. 36.' График для определения масштаба L по измерениям скорости ветра и температуры на высотах 0,5 и 2,0 м
Определение L
Воспользуемся соотношениями (5.5.1). Тогда
|
А Ѳ ( г 2, д ) |
\ в п _ |
Т„ |
|
|
ДW“ (z4, zn) |
Awn“ |
g'L |
|
Если считать, |
что Г0/Я ~ 28, то |
|
|
|
АѲ (*2, г,) ^ |
28 |
_ АѲя(з2/і , |
з ,,7.) |
|
A«2(z4. z3) |
Z, |
1 и п 2 ( г А 1 . |
z , J L ) |
|
Для фиксированных 2 ,- можно построить зависимость функ- % ции Ф от L, которая позволит определять L как функцию от
101
ДѲ/Дм2. В качестве примера на рис.. 36 показана зависимость L ■от ДѲ/Дм2 При 22 = 24 = 2,0 М, 2 і = 2 3 = 0,5 м.
Определение турбулентного потока тепла Р0 и затрат тепла на испарение LE
С учетом (5.5.1) можно записать, что
Ли-ДѲ = -------- |
|
0-,- д«„-дѳя, |
|
|
|
|
|
|||||
откуда |
|
|
Р• сп• У-"' |
|
|
|
|
|
|
|
|
|
|
|
|
Дм-ДѲ |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
(5.6.2) |
||||
|
|
|
- и •'<■ •с |
|
|
|
|
|
|
|||
Для определения турбулентного потока тепла по |
(5.6.2) нуж |
|||||||||||
но сперва на основании ДѲ и Дм найти L (допустим из рис. 36), |
||||||||||||
затем рассчитать zni |
и получить Дип и ДѲ„ . Более целесообраз |
|||||||||||
|
|
|
|
но построить номограммы, по |
||||||||
|
|
|
|
зволяющие |
определить |
Р0 как |
||||||
|
|
|
|
функцию ДѲ и Дм. Допустим, |
||||||||
|
|
|
|
что 22 = |
21 |
и |
Zi = |
2 3; зададим |
||||
|
|
|
|
для выбранных величин P0j- |
||||||||
|
|
|
|
ряд значений |
|
и.-;і ; найдем |
L h |
|||||
|
|
|
|
|
2 j |
L j, |
{lln2) j , |
(Un-j)j, |
|
n>\t |
||
|
|
|
|
(Ѳ„,)Б |
определим |
(ДиД) |
и |
|||||
|
|
|
|
(ДѲ„)/ |
и с помощью (5.6.2) пе |
|||||||
|
|
|
|
рейдем к Дм; и ДѲ,-. Построив |
||||||||
|
|
|
|
изолинии Я0у |
|
в системе |
коор |
|||||
|
|
|
|
динат Дм и ДѲ, получим номо |
||||||||
|
|
|
|
грамму |
для |
определения |
Р 0. |
|||||
|
|
|
|
В качестве примера на рис. 37 |
||||||||
|
|
|
|
показана |
зависимость |
/ДДДщ |
||||||
|
|
|
|
ДѲ) |
для z2 — 2.1 = 2,0 м, |
2 ) = |
||||||
|
|
|
|
= 2 а |
= 0,5 м. |
|
|
|
|
|
||
|
|
|
|
Нетрудно показать, что изо |
||||||||
|
|
|
|
линии Ро можно рассматривать |
||||||||
|
|
|
|
как изолинии LE, если счи- |
||||||||
|
|
|
|
тать ДѲ - *ч -і |
, |
Действитель- |
||||||
|
|
|
|
но, так |
как |
|
|
|
|
|
||
2 |
3 |
ли г* 'CFK |
LE* |
|
-U-E-L |
\u - \ q |
(5.6.3) |
|||||
Рис. 37. Номограмма для определе |
|
Лм„- Д(/„ |
||||||||||
|
|
|
|
|
|
|
|
|
||||
ния турбулетного потока |
тепла (Ро), |
|
|
|
|
|
|
|
|
|
||
и потока скрытого тепла |
(LE) |
по |
|
|
|
|
|
|
|
|
|
|
наблюдениям на высотах |
0 5 и 2,0 м |
|
|
|
|
|
|
|
|
|
102
и А</я —\Ѳ, |
(если- а т = я 7), то, приравнивая (5.6.2) |
и (5.6.3), |
получаем, что при условии PÜ— LE |
|
|
По найденной величине затрат тепла на испарение (потоку |
||
скрытого тепла) легко определить скорость испарения |
Е. |
|
|
Определение динамической скорости ѵ% |
|
Из (5.5.1) |
следует, что |
|
V
X• Ли
(5.6.4)
Дк„
Для определения динамической скорости по (5.6.4) нужно сперва найти І(Аи, ЛѲ), а потом Аи„. Более удобно, однако,
рассчитать номограмму |
для |
определения |
ѵ.А; (Ак, АѲ). При |
z2 = z4 и 2 \ = 2ъ зададим |
для |
выбранных величин АѲ , ряд вели |
|
чин L ;; найдем z2/L(-, Zi/L,, |
(и„2) , (ипѴп) |
, (ѲЯІ),-. (Ѳпі){, (Ѳя1)<, |
рассчитаем Лкя)г и (АѲ„),- и с помощью (5.6.1) перейдем к (Ак),-, после чего определим по (5.6.4) ѵ^:і. Результаты расчета пред ставлены на рис. 38, на котором в виде примера приводится за висимость ü* (Ам, АѲ) для z2 = zi = 2,0 м, Zi —г3 = 0,5 м.
Рис. 38. Номограмма для определения динамической скорости (i'*j по наблюдениям на высотах 0,5 и 2,0 м
Поскольку над водной поверхностью градиентные наблюде
ния выполняются редко (только |
в специальных экспедициях), |
то желательно уметь определять |
характеристики приземного |
слоя на основании стандартных гидрометеорологических наблю
дений. Предположим, что скорость ветра |
температура |
Ѳ„ и |
влажность воздуха ца измеряются на одной |
высоте г = о |
(в об |
щем случае высоты измерения и. Ѳ и q могут быть различными). В качестве второго уровня будем рассматривать уровень шеро ховатости Zo, при этом возникают по крайней мере два вопроса: как определить го и чему равны скорость ветра, температура и
.влажность при г = г0. Из простых физических соображений ясно,
103
что шероховатость водной поверхности должна зависеть от ло кальной скорости ветра и характеристик, определяющих разви тие волн. В общем случае следует считать, что
Zo = f(v*, g, Т, X, 1, р, р), |
(5.6.5) |
где g — ускорение силы тяжести; Т — поверхностное натяжение; А' — разгон ветра; t — продолжительность действия ветра; р —
плотность воздуха; р — плотность воды.
Используя теорию размерности, зависимость (5.6.5) можно
свести к безразмерному виду |
|
|
|
|
h l L —P T-g |
X- |
L s |
? |
(5.6.6) |
Так как р/р можно считать константой (значительно меньшей |
||||
■единицы), то Р(р/р) также |
стремится |
к |
постоянной |
величине. |
Из полученного соотношения можно оценить условия, при кото рые допустимо пренебречь влиянием отдельных факторов. Так как функция от значительно меньшего единицы безразмерного комплекса стремится к постоянной, то при выполнении условия
-р
Т «
g
можно не учитывать зависимость Zo от поверхностного натяже ния (она войдет в константу). Считая, что Г = 70 г/сек2 и ис пользуя известный из наблюдений факт, что ѵ... ~ 0,03 • ѵ, где о — скорость ветра на уровне судовых наблюдений, можно показать, что полученное неравенство выполняется для скоростей ветра около 10 м/сек и выше.
С другой стороны, можно считать, что функция от второго и третьего аргумента в (5.6.6) стремится к постоянной, если
x-g |
»1 |
и ~ > : |
V* |
|
V« |
откуда следует, что |
|
|
|
g' ’ |
t » L > |
|
g |
Это соответствует: А > > 1 см, t > > 0,03 сек. Очевидно, что шероховатость устанавливается быстро и для малых разгонов ветра. Следует оговорить, что го может зависеть от X и /, если происходит передача энергии от низкочастотной части к высоко
104
частотной части |
спектра |
(т. |
е. за |
счет |
неустойчивости |
крупных |
||||
волн образуются мелкие волны). |
|
|
|
|
||||||
Итак, при выполнении рассмотренных выше условий (5.6.6) |
||||||||||
можно записать в виде |
|
|
|
|
|
|
|
|||
|
|
|
|
|
z0= m |
|
|
|
(5.6.7) |
|
Для проверки справедливости |
(5.6.7) |
и определения |
величи |
|||||||
ны т запишем |
логарифмическую |
формулу для профиля ветра |
||||||||
|
|
|
J L |
— 2,3] / Г g |
. jK?.£o_ |
(5.6.8) |
||||
|
|
|
Г г |
|
* |
V |
т |
у |
г 20 ‘ |
|
В таком случае, |
если |
(5.6.7) |
выполняется, то в системе коорди- |
|||||||
нат |
'1 |
V- |
l g |
2 |
2 q |
результаты измерении должны рас- |
||||
У — и ( z , |
л |
— — — |
- ггг- |
V z zn
полагаться на прямой линии, тангенс угла наклона которой поз-
£q~ to
Рис. 39. Номограмма для определения турбулентных потоков и динамической скорости по стандартным" наблюдениям па
поверхности моря и |
на высоте |
14 м. Сплошные линии — |
Р и LE Е кал/слРмин-, |
пунктирные линии — ѵ* в см/сек; q— |
|
в |
г/кг; иа— в |
м/сек |
105
воЛяет определить т. Обработка наблюдений показала, что, не смотря на заметный разброс то'йек, они в среднем могут быть аппроксимированы прямой линией с /» = 0,075 [11]. Современные оценки т находятся в пределах 0,02—0,08.
Рассмотрим теперь второй из возникающих вопросов: чему равны скорость ветра, температура и влажность при z = z0. Скорость ветра на уровне параметра шероховатости должна быть равна скорости поверхностного течения. Так как последняя (за исключением сильных постоянных течений^ составляет всего
около 2% |
от скорости ветра, то |
можно считать, что |
при z —Zq |
|||||
// = 0 |
(т. е. |
Аи = ил). |
Обработка |
большого |
количества |
наблюде |
||
ний |
показывает, |
что |
для морских условий температура |
воздуха |
||||
и влажность на |
z = z0 практически равны |
температуре |
поверх |
|||||
ности и соответствующей ей насыщающей |
влажности. |
|
|
Итак, специфика расчета характеристик приземного слоя по стандартным гидрометеорологическим наблюдениям состоит в том, что іпіжнпй уровень (параметр шероховатости) не фикси рован, а находится в процессе расчета на основании (5.6.7). В качестве примера на рис. 39 приводятся номограммы для рас чета Р, LE и у* при z2 = z+ = 14 м.
106
I
VI. ПРОЦЕССЫ ВБЛИЗИ ПОВЕРХНОСТИ РАЗДЕЛА ВОЗДУХ-ВОДА
§ 1. Постановка задачи о динамическом взаимодействии погра ничных слоев океана и атмосферы
Воспользуемся физической аналогией процессов в погранич-
.ных слоях атмосферы и моря и на основании нелинейной теории пограничного слоя атмосферы сформулируем задачу о динами ческом взаимодействии пограничных слоев атмосферы и моря. Выпишем замкнутую систему уравнений, определяющих дина мические процессы в пограничных слоях [11]:
( 6. 1. 1)
/ |
(6.1.2) |
ЯД*,-); |
(6.1.3) |
k( = /, У Ь{; |
(6.1.4) |
\ |
(6.1.5) |
|
( 6. 1.6) |
Здесь І=1 соответствует атмосфере, і = 2 — морю: ugX, |
v gX— |
компоненты геострофического ветра; us2, vg2 — компоненты геострофического течения; Ѳі = Г, ©2= рг; осі = а т, а 2 = а р;
с4 = |
с’„ — удельная |
теплоемкость |
воды; Роі — турбулентный |
по |
|||||
ток |
тепла; |
Р02 — турбулентный поток плотности в море. |
|
|
|||||
Чтобы |
проинтегрировать (6.1.1) — (6.1.6) и определить |
вер |
|||||||
тикальное |
распределение |
г',-, |
k{, bt сформулируем следующие |
||||||
граничные условия: |
|
|
|
|
|
|
|||
|
|
2j->oo, Ui—Ugi, |
i'i |
-т'й/, bi~+0; |
|
(6.1.7) |
|||
|
|
zi~*zoir «.— «У, |
Vi-*V0- |
|
|
( 6 . 1. 8) |
|||
|
|
|
сіил |
— -- |
. |
du2 |
|
|
|
|
|
Pi |
äzl |
k, р2 |
—т- - |
|
|
|
|
|
|
|
|
U*•> |
|
|
|
||
|
|
|
->0 |
|
|
г»-*0 |
1 |
(6.1.9) |
|
|
|
|
d%\ |
|
. |
dv2 |
|||
|
|
|
— - |
|
|||||
|
|
*iPi |
|
|
|
||||
|
|
~dzy |
к*Рз |
dz. |
|
|
|
||
|
|
|
-0 |
|
|
-0 |
1 |
|
|
Ось 2 \ направлена вверх, ось z2 — вниз, начало координат на не возмущенной поверхности.
Условия (6.1.7) отражают известный факт, что за пределами пограничных слоев ветер н течения стремятся к геострофическим значениям, а турбулентность затухает. По мере приближения к поверхности раздела скорость ветра и течения стремятся к ско рости поверхностного течения и0, Ѵо (условие склейки скоро стей), а кинетическая энергия, как можно показать из анализа уравнения баланса энергии турбулентности, стремится к с"1;2г,2:і:і-. Наконец, можно считать, что у поверхности раздела должно вы полняться условие непрерывности потоков количества движения. Для"определения параметра шероховатости можно использовать
(5.6.7).
Zo;=p/M- |
(6.1.10) |
Если ввести новые переменные величины |
|
Ui=Ui-ut, V,—Vt-v,). ugl=ug{—a0, ve~ v g{—v 0 |
(6.1.11) |
и переписать для них систему уравнений и граничных условия, то задача сведется к рассмотренной ранее задаче о строении по граничного слоя атмосферы над неподвижной поверхностью (раздел IV, § 4). Таким образом, система уравнений для каж дого пограничного слоя может быть решена отдельно и искомые безразмерные величины выражены через некоторые универсаль ные функции anV kni, bni, зависящие для пограничного слоя атмосферы от
108
\Н = Ы - е |
—ахг„1 |
у9 |
|
pQl ?Ср |
||
|
, [t0r |
|
||||
|
" |
‘ ~T |
‘ г ^ Ѵ |
|||
|
|
|
|
|||
* . = - ! |
/ - |
l ^ . |
|
|
(6.1. 12) |
|
|
|
|
||||
nl |
l ; |
м |
2шг * |
|
|
|
для пограничного слоя моря от |
|
|
|
|||
|
|
~ а і~па |
|
g_ |
. Л 12. Pa |
|
И^Н’-ог^ |
г Ро2 = |
|
Рг |
‘ |
||
__ ^ |
f |
* л1 <> |
|
|
(6.1.13) |
|
|
|
|
||||
"л2“ |
Т 2’ |
|
|
|
|
|
При раздельном решении системы для атмосферы и моря ос тались неиспользованными условие склейки скоростей и непре рывности потоков количества движения у поверхности раздела. Именно эти граничные условия позволяют «склеить» два неза висимых решения, определить и0, .ѵ0 и связь между динамиче скими скоростями в атмосфере и море и получить выражение для геострофического коэффициента трения % и угла а между вектором касательного напряжения и вектором геострофическо го ветра (ось Ох направлена вдоль вектора,касательного напря жения на поверхности раздела)
xG “
1
(6.1.14)
|
|
•Д0 |
Здесь G — скорость геострофическрго ветра; |
||
*0п~ |
тп• 2шг• G х=«-х; |
m2u)2 • G 28 -X |
(считается, что Zo= Zq, т. е. равны размерные параметры шеро ховатости, со стороны воздуха и воды).
109