Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Мельникова И.И. Динамическая метеорология учеб. пособие для океанологов

.pdf
Скачиваний:
36
Добавлен:
23.10.2023
Размер:
5.62 Mб
Скачать

«Л- Я - К = 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*

 

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

 

 

 

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

107
С\ — с р ,

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

Соседние файлы в папке книги из ГПНТБ