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

книги из ГПНТБ / Процессы переноса в близи поверхности раздела океан - атмосфера

..pdf
Скачиваний:
9
Добавлен:
21.10.2023
Размер:
9.6 Mб
Скачать

8.3. Интегральный перенос каплями тепла и пара при шторме

Вертикальный перенос тепла (Р + ЕЕЕ)к и массы Ек всеми кап­ лями определяется выражениями:

( Р + ^ Р К) = ] [P«(r) + ^ ( r ) ] - g - d r .

(8.16)

rmln

(8.17)

где у (г) — функция распределения по размерам вертикального по­ тока капель (см~2/с) на нулевом уровне 2= 0. Входящие в (8.16) и (8.17) величины испарения и теплоотдачи одной капли Pi (г) и (г) вычисляются для данных условий приводного слоя воздуха так, как это было сделано выше (см. табл. 8.4). Если предположить допу­ стимость замены реального потока брызг монодисперсным потоком

со средним эффективным радиусом капель г, то выражения (8.16) и (8.17) упрощаются:

 

P *= P i{r)j .

(8.18)

 

Рк= E i {r)j.

(8.19)

Значения / (см_2/с)

при различных скоростях ветра

(для г —

= 0,05 см) определяются формулой (4.3).

различ­

Относительные

величины {Р+ЕЕ Е)К/(Р+ ЕЕ Е)т при

ных скоростях ветра приведены на рис. 8.1. Перенос тепла каплями

(Р + Е Е Е ) к вычислен здесь с использованием

(8.18) и (8.19) и дан­

ных табл. 8.4 при г = 0,05 см, а турбулентный

поток

{ Р +

ЕЕ Е ) т —

с использованием модели теплообмена (см. главу 6)

и зависимости

(7.28)

для

коэффициента

сопротивления

морской

 

поверхности.

Видно,

что

отношение (Р + ЕЕ Е ) К/ ( Р + Е ? Е ) Т

быстро

возрастает

при усилении скорости ветра и уже

при

«ю= 204-25

м/с прибли­

жается к единице.

 

 

 

 

 

 

 

Принятое нами значение

г = 0,05

см, как уже отмечалось, обе­

спечивает получение незавышенных

оценок

переноса

импульса,

тепла и пара брызгами. Представляет интерес оценка эффекта из­

менения

принятой величины г.

Воспользовавшись данными

табл. 8.4

и введя (при сохранении

данной водности, г/см3, привод­

ного слоя) приближенное условие

(8.20)

201

получаем:

(Р + ^ Е ) к

» 3 ,6 ,

( 8.21

г = 0,01 см

(P + J3TE) к

 

 

г = 0,05 см

 

 

(Я + ^ £ ) к ;Г = 0,10 см

^ 0 ,4 .

(8.22)

(Р + ^ £ ) к

 

 

г = 0,05 см

При более точном, чем по (8.20), определении / (гг) отношение потоков в (8.21) получилось бы несколько больше, чем 3,6, а в (8.22) — меньше, чем 0,4 (т. е. если бы были учтены различии

(Р*ХЕ)к

СЕ

 

СЕЮ

10

20

 

 

30 им/с

 

 

 

 

 

 

Рис. 8.1.

Отношение

вертикального

 

 

 

 

 

 

переноса

тепла брызгами

(P + J3?E

 

 

 

 

 

 

к турбулентному потоку

(P + JS?E) т

 

 

 

 

 

 

в зависимости от скорости ветра.

Рис. 8.2. Отношение коэффициента

 

 

 

 

 

в вертикальной

скорости капель

испарения СЕ(ит) при различных

разного размера

на уровне изме­

скоростях ветра к значению СЕ1о при

скорости ветра

и10= 10 м/с.

рения водности),

 

однако,порядок

 

1 — расчет по

модели,

учиты ваю щ ей испа­

величин не изменился бы.

рение

с

бры зг;

2 — эксперим ентальны е

Величины отношений (8.21) и

значения

по данны м

О куда и Х аям и (1959).

(8.22) показывают, что при всей

капель

по

размерам реальную

значимости учета

распределения

оценку переноса тепла и массы брызгами

можно

получить, введя

средний эффективный радиус капель. Удовлетворительное согласие рассчитанной нами зависимости коэффициента испарения

СЕ—

Ет+ Ек

(8.23)

^z) uz

 

 

от скорости ветра с экспериментами Окуда и Хаями (1959) и Эстлунда (1970) подтверждает этот вывод (рис. 8.2). Окуда и Хаями (1959), измеряя профили скорости ветра и влажности воздуха в аэрогидродинамическом канале, обнаружили, что коэффициент испа-

202

рения

 

*Дг2 da

du

 

СЕ

dz

dz

(8.24)

(aw az) uz

практически постоянный при скоростях ветра, меньших 15 м/с, на­ чинает быстро расти при больших скоростях.

Оценки испарения с поверхности океана при ураганной скоро­ сти ветра Эстлунд (1970) получил косвенным путем, использовав данные измерений влажности, концентрации пара тритиевой воды (НТО) и скорости ветра в тропическом урагане. Уравнения баланса массы воздуха, обычной Н20 и тритиевой НТО воды для слоя, где радиальная составляющая скорости ветра направлена к центру ура­ гана, составляют систему:

™ = 4 ? ^ Н ')’

(8-25)

w<lH = - l R ( v H ' y ) + E ~ P f ’

(8-26)

™ ^ = - ^ ( ъ Н ' д * ) + Е * - Р * г

(8.27)

Здесь w — вертикальная скорость воздуха; qH и -q'jj — отношения

смеси обычной и тритиевой воды на верхней границе слоя, толщина которого v — средняя для слоя Н' радиальная составляющая

скорости ветра (относительно центра

урагана);

q и q*— отноше­

ния смеси, осредненные по слою от 0

до Я ';

£

и Е* — скорости

испарения обычной и тритиевой воды; Pf и

Р* — интенсивности

выпадения осадков, образовавшихся в слое 0—Н' (см/с). Решение системы уравнений приводит к выражению для скорости испарения в кольцевой зоне поверхности океана, заключенной между окруж­ ностями с радиусами Ri и R2 ,(центром является центр урагана):

db

 

1

(8.28)

Здесь v и q — величины, осредненные по радиусу от Ri до /?2; 0-=

= q*/q\ индекс w относится к поверхности

океана, индекс b —•

к нижней границе облачного слоя.

 

Вычисленные Эстлундом по формуле (8.28) величины испарения

вместе с данными о перепаде влажности

в подоблачном слое и

о скорости ветра вблизи нижней границы облаков позволили нам

найти значения коэффициента испарения для

подоблачного слоя:

СЕН'

Е

(8.29)

Q'b) Мb

 

 

203

При этом переход от приведенных Эстлундом величин отношения смеси q (г/г) к плотности водяного пара (г/см3)

217 • 10~6

рд

Т623 + q

на высоте 1200 м (нижняя граница облаков) был выполнен с ис­ пользованием средних данных о вертикальном распределении тем­ пературы воздуха Т и давления р в ураганах (Риль, 1963; Данн, Миллер, 1964). Переход к коэффициенту испарения приводного слоя

СЕ— (а. «г) «Z

ЦЬ

(8.30)

«Z «г ’

где уровень z не превышает примерно 50 м, требует знания отно­ шений (awаь)/(аю ciz) и Ub/uz, Изменения скорости ветра с вы­ сотой в урагане весьма малы; вычисленные по данным Пальмена и Риля (1957) отношения скорости ветра на уровне 1000 мб (мюоо) к скорости на различных уровнях (ир) приведены в табл. 8.5 (для

двух значений R — расстояния

от центра урагана).

Независимо

от выбора R для уровня

р = 880 мб

(z~1200

м) можно принять

Uiooo/weeo = 0,95.

 

 

 

 

 

 

 

 

 

Таблица 8.5

р м б ............................... . . .

1000

900

800

700

«1000/ “р при:

., . .

1,00

0,95

1,00

1,00

R — 60 м. миль . . .

R — 120 м. миль . . .

., . .

1,00

0,93

0,96

1,01

Так как в урагане «Фейс», исследованном Эстлундом (1970), давление у поверхности океана составляло 1004—1005 мб, можно полагать, что давление 1000 мб соответствует высоте 40—50 м. Предположив, что вертикальный профиль скорости на высотах от 10 до 50 м логарифмический, и приняв, что «юоо/мь= 0,95, находим

ц10—

0,95ий

(8.31)

 

Здесь х = 0,4 — постоянная Кармана; Сию— коэффициент сопротив­

ления морской поверхности, оценки которого для больших скоро­ стей ветра приведены в главе 5.

Среднее отношение перепада влажности в приводном слое к пе­ репаду в слое поверхность воды—уровень 880 мб, по Пальмену и Рилю (1957), можно принять равным 0,60.

Величины коэффициента испарения для уровня 10 м, найденные по данным Эстлунда (1970) с использованием (8.30) и (8.31) и эм­ пирических соотношений, приведены на рис. 8.3. Они неплохо со­ гласуются с зависимостью, полученной для интервала скоростей

204

ветра 15—25 м/с в результате анализа роли брызг в переносе тепла и влаги, и проэкстраполированной линейно до скорости 30—35 м/с.

С£ Ю3

Рис. 8.3. Зависимость коэффициента испарения Се от скорости ветра на высоте 10 м цщ.

1 — расчет с использованием модели, учитывающей испа­ рение с брызг; 2 — экстраполяция с учетом значений СЕ,

вычисленных по данным Эстлунда (1970) (3).

Такое согласие значений коэффициента испарения СЕ, полученных на совершенно различной физической основе, подтверждает бли­ зость наших оценок к реальным величинам.

8.4.Роль штормов в тепло- и влагообмене океана

иатмосферы за большие интервалы времени

Относительный вклад штормов в макровзаимодействие океана и атмосферы определяется как рассмотренным выше усилением энерго- и массообмена при штормах, так и повторяемостью послед­ них. Некоторые оценки этого вклада были сделаны Решетовой (1972), использовавшей гипотезы о линейной зависимости коэффи­ циента сопротивления Си от скорости ветра и о равенстве коэффи­ циентов теплообмена и сопротивления морской поверхности Св = Си-

Оценка климатологических значений коэффициента Се с учетом рассмотренной выше специфики штормового взаимодействия оке­ ана и атмосферы сделана Ариель и др. (1972) . В этой работе рас­ смотрена повторяемость скорости ветра, превышающей 17 м/с, для различных районов Мирового океана; выше отмечалось, что именно при скорости ветра 17 м/с происходит «включение» штормовых ме­ ханизмов и начинается быстрый рост коэффициентов обмена. Сред­ няя для всего Мирового океана повторяемость штормовых

205

скоростей оценивалась по ежемесячным картам Американского мор­ ского климатического атласа мира (1969); среднегодовая повторяе-

Рис. 8.4. Годовой ход вероятности скоростей ветра, превышающих 17 м/с.

1 — средний для всего Мирового океана, по данным Аме-

риканского морского климатического атласа мира

(1969);

2 — средний для океанов

северного полушария,

широта

больше 30° с., по данным Американского морского атласа

(1969); 3 — вычисленный

для той

же области (широта

больше 30° с.) по формуле

(8.32) с

использованием сред­

них значений дисперсий по данным Крауса и Моррисона

(1966).

мость штормов (и ^ 1 7 м/с) по этим данным всего 3,2% (см. кри­ вую 1 на рис. 8.4). Эта величина явно занижена, так как на картах проведены изолинии с повторяемостью штормов больше 5%, а мень­

 

 

 

шие повторяемости

не учитываются.

 

 

 

Кроме того, пути транспортных су­

 

 

 

дов, по наблюдениям которых соста­

 

 

 

влены

карты, прокладывают

так,

 

 

 

чтобы

по

возможности

избежать

 

 

 

встречи со штормами.

Наблюдения

 

 

 

кораблей

погоды

позволяют

более

 

 

 

достоверно и детально оценить функ­

 

 

 

цию

вероятности

распределения

 

 

 

скоростей ветра,

по крайней

мере

 

 

 

для северной части

Атлантического

Рис. 8.5. Повторяемость скоро­

океана. На рис. 8.5 приведены функ­

ции распределения вероятности зна­

стей

ветра,

превышающих

чений скорости

ветра,

превышаю­

17 м/с,

по данным кораблей

погоды А (1) и Е (2) за зим­

щих 17 м/с; они

построены по дан­

ние месяцы

1958—1960 гг.

ным наблюдений кораблей

погоды

 

 

 

в зимние

месяцы

(декабрь—фев­

раль)

за три года (1958—1960). Используя эти кривые, можно вы­

числить интегральную повторяемость штормовых скоростей ветра Q

(« ^ 1

7 м/с). Средняя по всем кораблям погоды Q (м ^17) оказа-

206

лась равной 17% (табл. 8.6). Если по Американскому морскому кли­ матическому атласу мира (1969) вычислить Q (ы^Э= 17) для района расположения кораблей погоды в Атлантике (широта больше 30°)» то среднее значение для трех зимних месяцев оказывается равным 9,2% (см. кривую 2 на рис. 8.4). Как и предполагалось, карты Аме­ риканского морского атласа занижают повторяемость штормов.

 

 

 

 

 

 

Таблица 8.6

 

а (и ^1 7 )

 

 

 

 

 

Корабль

по натур­

 

 

и м/с

 

^ 0

К Л И М

погоды

по форму­

“шт м/с

и 0 М / С

 

 

ным дан­

ле (8.32).

 

 

 

с 0

 

ным

 

 

 

 

А

0,21

0

,2 6

12,4

2 0 ,2

10,4

1,35

в

0 ,1 6

0

,1 9

11,9

2 1 ,2

10,1

1,28

с

0 ,2 0

0,21

11,9

2 1 ,3

9 ,5

1 , 35.

D

0 ,2 0

0,21

12,0

20,7

9 ,9

1,35

I

0 ,1 7

0

,2 6

12,7

19,9

11,2

1 , 26,

J

0 ,2 0

0

,2 9

12,8

2 0 ,2

11,0

1,32

К

0 ,1 5

0

,1 0

10,3

19,9

8 ,6

1 , 29'

м

0 ,1 9

0

,1 0

10,5

19,8

8 ,3

1 , 36,

Е

0 ,1 0

0

,0 6

9 ,8

19,7

8 ,7

1,20

Среднее

0,17

0 ,1 9

11,6

2 0 ,3

9 ,9

 

Q

Интересно сравнить

интегральную повторяемость

штормов;

17), полученную

по наблюдениям, с величиной,

которую

можно вычислить по средней скорости ветра и и дисперсии скоро­

сти а2и в предположении

максвелловского закона

распределения

модуля скорости:

 

 

 

 

Q(«KP) =

j q (и) d u = e “кР/"2 .

(8.32)

2и

“кр

 

 

 

Здесь q (ы )= -^ -е ~ “2/“2— функция распределения модуля скорости

и2

ветра; цКр=17 м/с; значения и и о2 для всех месяцев года и всех

кораблей погоды приведены в работе Крауса и Моррисона (1966). Сопоставление реальной повторяемости штормов и вычисленной по формуле (8.32) показывает, что разность между этими величи­ нами зависит от средней скорости ветра (см. табл. 8.6). Значения Q(«kP= 17), средние по всем кораблям погоды, иллюстрирует кри­ вая 3 на рис. 8.4; можно полагать, что она дает незавышенную по­ вторяемость штормов для летних месяцев.

Значения климатологического коэффициента

теплообмена

С9кл1ш оценивались на основе простого равенства

 

^ в к л и м ^ ^ ® [1 Q (« к р )] И о Д 0с А ( « 0 ,

 

+ Q ( ^кр) Cg (^l'jt) ^штД0шт.

(8.33)

207

Здесь и и Д0— средние скорость ветра и перепад температуры вода—воздух за весь рассматриваемый интервал времени; и0, Л0о

и Де0 (перепад влажности в миллиметрах) — средние значения, по­ лученные при исключении из ряда наблюдений случаев со штормом;

ышт и Д0шт — средние значения тех же величин при штормах.

Значения «о, и, ишт приведены в табл.8.6. Вопрос о соотношении перепадов температуры (и влажности) при штормовых и нештор­ мовых условиях весьма важен. Можно было полагать, что усиление турбулентного обмена в приводном слое воздуха во время шторма

приводит к уменьшению перепадов Д0ШТ и ДеШт по сравнению

t Д0 и Ае; в этом случае вклад штормов в тепло- и влагообмен оке­ ана и атмосферы уменьшился бы. Для проверки этого предположе­

ния были рассчитаны Д0ШТ и Аештдля всех кораблей погоды за зим­

ние месяцы 1958—1960 гг. и сравнены с Д0 и Ае за тот же период при всех скоростях ветра. Как видно из табл. 8.7, перепады темпе-

Таблица 8.7

Корабль

погоды

А

В

С

D

Г

J

К

м

Е

Среднее

£1 1

О

 

о

+3,19

2,50

1,37

3,09

3,16

1,80

0,85

3,14

1,09

2,24

Дб °с

Д<?шт мб

Ае Мб

+2.89

+3,07

+2,28

3,27

2,22

2,07

1,35

2,50

1,90

3,55

3,47

4,87

2,85

3,49

2,81

1,81

3,14

2,56

1,07

3,09

2,55

3,92

3,31

2,64

1,35

5,75

4,58

2,37

3,33

2,91

ратуры и влажности при штормах практически не отличаются от средних перепадов. Это обстоятельство позволило использовать

формулу

(8.33)

для

получения

грубых

оценок

отношения

Секлим/Се («,

Д0,

Ае)

(см. табл.

8.6).

При этом

величина

Св(и,

Д0,

Ае)

определялась по номограммам (Бортковский,

1971),

а для

коэффициента теплообмена при шторме было принято

посто­

янное значение Се (мшт) = 3,0 • 10~3, которое, очевидно, не является завышенным. Приведенные значения (см. табл. 8.6) отношения

Секлим/Се (и, Д0, Ае) иллюстрируют вклад штормов в теплообмен океана с атмосферой, составляющий для Северной Атлантики в зимние месяцы около 30%: от 20% для корабля погоды Е до 36% для корабля М. Тепло- и влагообмен в умеренных широтах в лет­ ние месяцы значительно меньше, чем в зимние. Поэтому, несмотря

208

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

Однако формула (8.33) дает лишь грубые оценки по двум при­ чинам: во-первых, коэффициент теплообмена при шторме сильно зависит от скорости ветра (Бортковский, 1972) и вклад штормов в теплообмен нужно считать с учетом дифференциальной функции распределения штормовых скоростей ветра; во-вторых, следует бо­ лее точно, с учетом флуктуаций скорости ветра и перепадов темпе­ ратуры и влажности,-вычислять потоки при нештормовых условиях.

/9 Л 1

I

I

I______ I______1_______ I______ I____I_______ I_____ I

2

3

4

5

6

7

8

9

10

11 й м /с

Рис. 8.6. Эмпирическая зависимость безразмерной дисперсии модуля скорости ветра от средней скорости.

1 — по данным Американского морского климатического атласа мира (1969); 2 — по данным Крауса и Моррисона (1966). Стрелки показывают средние квад­ ратичные отклонения данных.

Следовательно, эффективное значение коэффициента испарения

для данного

интервала скорости ветра А и определяется

выраже­

нием

Ц+Ди

 

 

j q (и) СЕ (и) du

 

 

---------------•

(8-34)

 

| q (и) du

 

 

U

 

где q (и) — функция распределения вероятности скорости ветра. Распределение вероятности среднечасовых значений модуля скоро­ сти ветра над океаном достаточно хорошо описывается формулой Максвелла (8.32); однако анализ натурных данных показал, что

значения отношения «2/н = (о2/ы ) +1 отклоняются от величины

4/я, соответствующей максвелловскому распределению. Это откло­ нение закономерно связано со средней скоростью ветра в данном районе в данный сезон (рис. 8.6). Эмпирическая зависимость

14 Заказ № 154

209

аи/и от и, приведенная на рисунке, означает,

что при

8 м/с

функция распределения становится более узкой,

чем максвеллов­

ская, а при ы ^ 8 м/с — белее широкой. Сопоставление реальных

повторяемостей штормового ветра различной силы,

рассчитанных

по данным Американского морского климатического

атласа мира

(1955—1959), с результатами расчета по кривой

распределения

Максвелла, при котором значения аи/и определялись по кривой на рис. 8.6, показало практически полное их совпадение. Поэтому для

всех районов Мирового океана повторяемость

скорости можно

определить по формуле Максвелла

 

q { u ) = ^ e - ^ ^ \

(8.35)

если величина и2 взята из наблюдений. Для удобства вычислений эффективного значения коэффициента испарения С'Е формула

(8.34) записывается в виде

2

Q (“»'■ “ г + Ди) С Е (“г + -у -)

Q '

____________.___________

Е

2 Q

+ Ди)

 

i

 

где Q («,, Ui + Au) — вероятность нахождения скорости ветра в ин­ тервале А и; CE{^U i+ ^ - j — коэффициент для середины интервала

u .i + Au/2. Из формулы (8.35) следует

и2

_ (и + Дц)2

 

Q(u, и + А и )= е * - е

* .

(8.36)

По формуле (8.36) вычислены распределения вероятности для еди­ ничных интервалов скорости штормового ветра (17—18, 18—19 и т. д. до 30 м/с) при различных значениях средней за месяц или сезон скорости ветра. Поскольку в Американском морском клима­

тическом атласе мира и в других пособиях приведена

сила ветра

в баллах, эффективные значения С'

вычислены для

интервалов

скорости 17—21 м/с (8 баллов), 21—24 м/с (9 баллов)

и 24—30 м/с

(10—11 баллов)

(табл. 8.8)

в зависимости от средней скорости и.

 

 

 

 

Таблица 8.8

Интервал

 

Средняя скорость ветра (м/с)

 

 

 

 

 

скорости (м/с)

6

8

10

12

 

17-21

1,43 • Ю-з

1,49 ■Ю-з

1,50 • Ю-з

1,50 ■Ю-з

21—24

1,80

2 , 0 2

2,04

2,06

24-30

3,14

3,86

4.11

 

 

 

 

210

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