книги из ГПНТБ / Процессы переноса в близи поверхности раздела океан - атмосфера
..pdf8.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