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

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

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

меняются вдоль направления движения, то справедливо соотно­ шение

j = N V ,

(4.14)

где N — концентрация частиц (см-3). Это соотношение легко обоб­ щается на рассматриваемый нами случай, когда через выбранное сечение проходят два равных по величине и противоположно на­ правленных потока (все капли, взлетающие с поверхности воды, за­ тем падают обратно): '

 

j = N

ПГ+11 w.\

(4.15)

 

I W+\ + \W-\

 

 

 

Здесь

\ W+ \ — модуль вертикальной составляющей скорости капель

при движении вверх; \W-\ — то же при движении

вниз; в общем

случае

\W+ \ ф \ W- \.

 

 

Так как концентрация частиц N определяется практически по

данным измерений в конечном объеме, в пределах

этого объема,

имеющего высоту Аг, изменение скорости —-—Аг, предполагав- dz

шееся выше нулевым, может оказаться существенным. Обобщение формулы (4.15) дает

j

М'{г2- г х)

(4.16)

dz

dz

 

 

I W+ (г) I

I W- (г) |

 

Концентрацию брызг (см-3) можно вычислить по измеренной

Монагэном (1968) водности, предполагая, что

радиус капель г =

= 0,05 см. Концентрация,

вычисленная таким

образом, относится

к уровню 13 см над водой,

а объем,

в котором проводились измере­

ния, ограничен уровнями 22=14,5 см и 2i = 11,5 см. Значения верти­ кальной скорости в области Az= z2zx находились при решении системы уравнений движений (4.5). Полученное этим способом зна­ чение потока капель при «=15 м/с, /i5= 0,19 см- 2>с-1, весьма бли­ зко к приведенному выше. Таким образом, оба способа дают практи­ чески одно и то же число капель, взлетающих с единицы поверхно­ сти моря за 1 с при скорости, начиная с которой зависимость j(u) делается экспоненциальной.

4.4. Результаты расчетов

Порядок решения системы уравнений движения (4.5) был сле­ дующим:

1) вычисление числа

2r \ f w l + {Un- u f .

2) определение коэффициента сопротивления CD(Re„) по гра­ фику рис. 4.9;

9*

131

3) вычисление составляющих ускорения капли

 

{■§-),=тt

 

W

'■)/ ^ + № - я ) а.

 

d W

3 р

W , Y w l + ( U , ~ u f + g - ,

 

dt

 

4

 

 

 

4)

вычисление

составляющих

скорости следующего, (п+1)-го

шага

 

 

 

 

 

 

 

 

и . +,= и . + ы ( - § . \ ,

 

 

 

Wn+l= W n+ A t (t ddtW

и ординаты капли

 

W„,, + IF

 

 

 

 

 

 

Zn+ i^zn+ A t - ^ \ -----"

При

расчетах

принимались следующие значения параметров: г =

= 5*10-2 см,

р= 1,2 -10—3 г/см3, р№= 1,0 г/см3, g=980 см/с2, шаг по

времени Д^=0,01 с,

начальное

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

0= 300 см/с.

 

 

 

 

Таким образом были получены в последовательные моменты времени значения U, W, z, Re. Расчет заканчивался, когда оказыва­ лось, что 2П= 0, т. е. когда капля «падала в воду». Вычисления ве­ лись при значениях скорости ветра и от 10 до 25 м/с. Для вычисле­

ния переноса импульса нужно знать

приращение горизонтальной

составляющей количества движения

капли за время ее полета; оно

определяется при принятом нами условии (Н0= 0) выражением

xt= M U f (W 0,

и),

(4.17)

где Uf — горизонтальная составляющая

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

в момент

падения в воду. Значения Uf/u в зависимости от скорости ветра при­ ведены на рис. 4.10. Зависимость Uflu от и довольно слабая: при изменении скорости ветра от 10 до 25 м/с Uflu меняется от 0,75 до 0,60.

Полный перенос количества движения от атмосферы к воде всеми каплями в принятой нами модели описывается формулой

)

(4.18)

= ^ i = M U f { W a, u )j

1

 

или с учетом (4.3)

 

-zK=M(Jf (W0, и)у,5е°'40(“_15).

(4.18')

Подстановка в (4.18) M = kh,nr3 при г = 0,05 см,

Uf (W0, и)

с рис. 4.10 и величины j15= 0,2 см- 2-1 позволяет вычислить верти­ кальный поток количества движения, переносимого брызгами при скорости ветра и от 10 до 25 м/с (табл. 4.1).

132

и м / с .......................

тк г/ (СМ • С2) . . .

тт г/(см - с2) . . .

^kI^T .......................

 

 

 

 

Таблица 4.1

. .

10

15

20

25

. .

0,02

0,12

1,14

10,5

. .

2,2

6,0

12,5

22,3

. .

0,01

0,02

0,09

0,47

Величины турбулентного потока количества движения тт, при­ веденные в табл. 4.1, вычислены по формуле (4.12).

Сравнение величин тгк и тт показывает, что при скорости 25 м/с перенос импульса каплями становится сравнимым с турбулентным

потоком. Величинами пото­

Си-го3

 

 

ков импульса, приведенны­

 

 

ми в табл.

4.1,

можно

вос­

 

 

 

пользоваться

для получения

 

 

 

зависимости

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

 

 

 

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

морской

по­

 

 

 

верхности

Си

от

скорости

 

 

 

ветра и:

тт +

 

 

 

 

 

 

 

С„

 

(4.19)

 

 

 

 

 

ри 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Кривая Си (и), приведен­

 

 

 

ная на рис. 4.11, удовлетво­

 

 

 

рительно согласуется с не­

 

 

 

многочисленными

экспери­

 

 

 

ментальными и другими рас­

 

 

 

четными данными, получен­

 

 

 

ными при

штормовых усло­

 

 

 

виях. Согласие это довольно

 

 

 

сильно

зависит

от

выбран­

 

 

 

ной нами

зависимости

для

 

 

 

тт;

кроме

 

того,

следует

Рис. 4.11. Коэффициент сопротивления мор­

учесть,

что оценки Миллера

ской поверхности при штормовых и ураган­

(1964)

также получены при

ных

ветрах.

ряде существенных

допуще­

1 — по данным Миллера

(1964); 2 — по данным By

ний. Что касается

данных

(1969а); 3 — расчет по модели

переноса импульса

брызгами с использованием

выражения (4.12);

By (1969а) и Маковой (1968),

4 — по Кузнецову

(1970).

то

в них отчетливо

прояв­

 

 

 

ляется резкое увеличение коэффициента сопротивления при «ю~ =«15-4-17 м/с. При этой же скорости начинается интенсивное обра­ зование брызг на поверхности воды (Тоба, 1961; Окуда и Хаями, 1959). Следовательно, несмотря на количественные расхождения наших результатов при больших скоростях с данными By и Мако­ вой, последние качественно подтверждают роль переноса импульса брызгами. Не останавливаясь подробно на причинах количествен­ ных расхождений, отметим, что By (1969а) представляет осредненные величины; однако отмечалось (Блинов, 1950; Монагэн, 1968), что при одной и той же скорости ветра интенсивность образования брызг может сильно меняться в зависимости от особенностей вол­

нения. Очень высокие значения Си,

полученные Маковой

(до 12,9X

Х10“3 при ц10 = 23,8 м/с), видимо,

определяются тем,

что при

133

глубине моря в точке измерений 22 м высота волн 3%-ной обеспе­ ченности достигала 8 м; когда глубина стала значительно меньше длины волны, равной при высоте волн 8 м примерно 60 м, должно было происходить резкое возрастание крутизны волн и обрушение гребней.

Упрощения, использованные при решении задачи о переносе им­ пульса брызгами, позволяют рассматривать результаты как до­ вольно грубую оценку. Качественно можно оценить роль использо­ ванных упрощающих предположений. Предположение о вертикаль­ ности вылета капли из воды (t/o = 0), видимо, не приводит к большим погрешностям, так как зависимость Ujju слаба, а введение и 0фО эквивалентно некоторому, относительно небольшому изменению скорости. Специально проведенные вычисления при разных значе­ ниях Wo также показали, что варьирование этой величины не при­ водит к существенным изменениям значений Uf.

Введение профиля скорости ветра представляется весьма суще­ ственным. Однако на малых высотах полета капли профиль и (г) должен сильно меняться над различными точками волны; данные об изменении «(£) при ?< 1 м (£ — расстояние до волновой поверх­ ности), полученные при шторме, нам пока неизвестны. Если же для простоты считать, что профиль скорости остается логарифмическим вплоть до самой поверхности, т. е.

— In- у- ,

(4.20)

где иг— скорость на уровне 210 м, «ю — на уровне

10 м, и* =

= «юУ(1,00 + 0,07и1о) 10_3— динамическая скорость,

определяемая

с помощью (4.12), то величина Uf/ul0, вычисленная

по описанной

выше схеме, оказывается равной 0,39 (при ию=10 м/с) и 0,34 (при ы10 = 25 м/с). Таким образом, весьма приблизительный учет зависи­ мости скорости ветра от высоты уменьшает величину тк примерно на 40%; однако порядок величины сохраняется. Если неучет про­ филя «(г) несколько завышает результат, то неучет спектра брызг

и расчет при г = 0,05 см его занижают. Действительно, если учесть наличие в приводном слое более мелких брызг, то, вероятно, для них величины Uf/u окажутся ближе к единице, чем полученные для

/ =0,05 см.

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

Ч а с т ь II

ПРОЦЕССЫ ТЕПЛО-И МАССОПЕРЕНОСА

Г л а в а 5

ОСОБЕННОСТИ РАСПРЕДЕЛЕНИЯ ТЕМПЕРАТУРЫ В ПРИПОВЕРХНОСТНОМ СЛОЕ ВОДЫ

Во многих обобщающих работах по вопросам физики привод­ ного слоя, начиная, очевидно, с работ Монтгомери (1947), термиче­ ские особенности поверхностного слоя моря рассматриваются как один из важнейших аспектов проблемы термического взаимодейст­ вия воды и воздуха. Качественный физический анализ процессов теплообмена на границах раздела, выполненный, например, Тимо­ феевым (1963), Лайхтманом (1961), указывает на то, что в обычных условиях в верхнем слое воды, непосредственно участвующем в те­ плообмене вода—воздух, имеются отрицательные температурные градиенты, т. е. температура у поверхности ниже, чем на некоторой глубине. Причиной этого является то обстоятельство, что, за исклю­ чением редких аномальных условий, водная поверхность отдает те­ пло в атмосферу путем и лучистой, и турбулентной теплопередачи при малом поглощении солнечной радиации поверхностным слоем, и для выполнения условий баланса тепла необходимо наличие ком­ пенсирующего эти потери потока из нижележащих слоев к поверх­ ности. Этот эффект наблюдается на фоне ослабления процессов об­ мена в воде по мере приближения к поверхности, что создает доста­ точно большие для инструментального измерения величины обратных температурных градиентов. Ослабление процессов обмена в воде вблизи поверхности приводит к формированию своеобразного ре­ жима передачи не только тепла, но и момента количества движения и соли. Однако наибольший интерес вызывают особенности тепло­ передачи, в связи с чем в последнее время получил распространение термин «холодная пленка». Под пленкой обычно подразумевается слой молекулярной теплопроводности. Однако использование та­ кого определения не всегда удобно. Например, разность величин температуры поверхности, определенной инфракрасным радиомет­ ром и по методике обычных судовых наблюдений, не будет равна перепаду температуры на границе пленки, как видно из рис. 5.1.

На рис. 5.1 I — область молекулярной теплопроводности, II — область молекулярной и турбулентной теплопроводности, III — об­ ласть турбулентной теплопроводности, Н — граница исчезновения

135

температурных градиентов, Д0о-ьб— перепад температур на грани­

цах пленки, 0№— температура на некоторой глубине, принимае­ мая в качестве 0о при гидрологических измерениях. Видно, что об­ щий перепад температуры Л0о-н отличен от перепада температур на границах пленки A0o-fte на величину А0л6-н.

Рис. 5.1. Схема температурного профиля в верхнем слое воды.

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

5.1.Условия образования холодной пленки

Вобмене теплом вода—атмосфера за счет переноса длинно­ волновой радиации участвует слой воды толщиной в доли милли­ метра. Относительно малое поглощение коротковолновой радиации

вэтом слое приводит к тому, что в первом приближении уравнение теплового баланса может быть записано в виде

 

- Е ьъ = & Е + Р + В ,

 

(5.1)

где Еэф— эффективное

излучение; £ ? Е — затраты

тепла на

испа­

рение (скрытый поток

тепла); Р — турбулентный

поток

тепла

вода—атмосфера; В — поток тепла в верхнем слое воды. Так как толщина рассматриваемого слоя воды достаточно мала, то изме-

136

нением теплосодержания слоя можно пренебречь I

cppwH -----=

'

d t

= 0). Этот член становится сравнимым с В при Я, равном 10 м, для слоев воды толщиной в несколько сантиметров он на 2—3 порядка меньше В. Тогда

В --

дв

 

 

 

 

дН я-*о

 

Полагая линейной зависимость 0 (Я)

и полагая /Cn,= const в урав­

нении теплопроводности для этого слоя

 

д

К„

ав

О,

 

дН

 

дН

 

 

можно записать величину В в виде

 

 

 

 

 

 

(5.2)

где A = p wcPKw — среднее значение коэффициента

теплопроводно­

сти в слое 0—Я.

 

 

 

 

Тогда уравнение теплового баланса можно записать в виде

( в о - в » ) = - 4 ( £'»Ф+ -2’Е + Р ) -

(5-3)

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

Для иллюстрации влияния разных форм теплообмена на форми­

рование (0о— 0,о)

приведем в табл. 5.1

отношение величин слагае­

мых потока тепла

к общему потоку

тепла (в процентах) для

оз. Севан и Северной Атлантики (Тимофеев, Малевский—Малевич, 1967). Из таблицы видно, что роль лучистого теплообмена вода— атмосфера сравнима по величине с турбулентными составляющими, а в условиях озера может иногда иметь преобладающее значение.

 

 

 

 

Таблица 5.1

Месяц

I

II ш I V

V VI V I I V I I I I X

X XI XI I

оз е р сI С е в а н

£эф

.......................................

44

45

60

72

64

57

43

41

37

38

42

42

З '

Е ...................................

25

28

29

24

29

37

58

50

51

42

39

33

Р ..........................................

 

31

27

11 .

4

7

6

9

9

12

20

19

25

 

 

С в е р н а я

Ат1л а н т и к а

 

 

 

 

 

£Эф .........................

• 16

18

18

20

32

36

32

28

24

22

17

16

&

Е ...................................

50

51

64

59

55

54

58

60

66

58

50

52

Р ..........................................

 

34

31

18

21

13

10

10

12

10

20

33

32

137

Рассмотрим теперь вопрос о теплопроводности поверхностного слоя. По аналогии с приземным (или приводным) слоем атмо­ сферы, для которого справедливо условие

Нш К (г)—*- v, z ^ 0

где v — молекулярная вязкость воздуха, можно предположить на­ личие следующего условия для верхнего слоя воды:

lim Л (Н ) -*■ Лм,

 

 

я

— о

 

 

где Лм — коэффициент молекулярной теплопроводности воды.

На основании такого предположения Сондерс

(1967)

опреде­

ляет перепад температур на границах пленки как

 

 

 

 

 

(5.4)

где rs— безразмерная константа;

В — суммарный

поток

тепла

вода—атмосфера. При этом толщина вязкого подслоя определяется как

Сондерс указывает, что его наблюдения дали величины rs в пре­ делах 5—10 со средним значением rs = 7 при величинах (0о— 0Ю) =

= -(0,24-0,3)° С.

(1972) привели величины г3, полученные в ла­

Полсон и Паркер

бораторных опытах,

причем определенные ими

величины rs зна­

чительно больше, чем у Сондерса (rs= 1 5 ± l).

Отметим попутно,

что указанные авторы ставили своей целью выяснить влияние раз­ ных механизмов теплообмена на формирование холодной пленки, однако условия этих опытов не позволяют полностью отождествить рассмотренный теплообмен вода—-воздух с теплообменом в есте­ ственных условиях, так как единственным механизмом, вызываю­ щим охлаждение поверхности, в них являлся процесс испарения

(Р/ J?E = 0,014-0,07, ЕЭф/,£’Е = 0,004-0,03). В то же время в есте­ ственных условиях структура полного потока тепла вода—атмо­ сфера существенно иная (см. табл. 5.1).

Экспериментальное подтверждение существования вязкого под­ слоя в воде вблизи поверхности можно найти в работах Мак-Али-

стера

(1964), Мак-Алистера и Мак-Лейша (1969),

Мак-Алистера

и др.

(1971), где приводятся данные термического

зондирования

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

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

138

ловое излучение воды, относящееся к двум различным глубинам. Соответствующая методика пересчета дает возможность получить информацию о средних температурах этих слоев. Мак-Алистер ис­ пользовал радиометр с пропусканием в участках спектра 2,0—2,4 и 3,5—4,0 мкм, что позволяет определять средние температуры слоев воды толщиной 0,06 и 0,50 мм. Проведя в естественных усло­ виях измерения градиентов температуры в воде и рассчитав при этих же условиях теплообмен вода—атмосфера, Мак-Алистер при­ шел к выводу, что условие неразрывности потоков выполняется при наличии молекулярного теплообмена в рассматриваемом слое.

 

 

 

 

Таблица 5.2

№ экспери­

 

В '

В

В'

д 8 0- А 8

 

 

мента

кал/(см2*мин)

в

 

 

 

 

1

0 , 1 5

0 , 2 7

0 , 2 6

1 , 0

2

0 , 1 2

0 , 2 2

0 , 2 2

1 , 0

3

0 , 1 9

0 , 3 4

0 , 2 2

1 . 5

Среднее

0 , 1 5

0 , 2 8

0 , 2 3

1 , 2

Результаты Мак-Алистера (1964) приведены в табл.

5.2. Здесь

A0o-/j6— перепад температур в слое 0,06—0,50 мм;

В ' — поток

тепла в этом слое, обусловленный чисто молекулярной

теплопро­

водностью; В — полный поток тепла

вода—атмосфера. На основа­

нии этих результатов Мак-Алистер

и др. (1971) провели ряд на­

блюдений над океаном, используя полученные данные для опреде­ ления полного потока тепла вода—атмосфера и основываясь на предположении о молекулярном теплообмене в пленке. Мак-Али­ стер и Мак-Лейш (1969) оценили влияние разных видов теплооб­ мена в воде на формирование холодной пленки. Их расчеты пока­ зали, что влияние радиационного теплообмена в воде несущест­ венно для формирования температурного профиля в верхнем слое на глубинах, превышающих 20 мкм. Определяющую роль в форми­ ровании температурного режима верхнего миллиметрового слоя воды играет молекулярная теплопроводность. При этом важно, что, согласно полученным результатам, наличие ветра скоростью до 10 м/с не приводит к замене режима молекулярной теплопроводно­ сти соответствующим турбулентным режимом, во всяком случае, в пределах верхнего полумиллиметрового слоя. Это свидетельствует о возможности сохранения холодной пленки на поверхности океана при развитии волнения. Некоторые данные о влиянии волнения на градиенты температуры будут приведены ниже.

Можно упомянуть еще некоторые экспериментальные работы, показывающие определяющую роль молекулярного теплообмена вблизи поверхности воды. Таковы, например, работы Болла (1954), Ролля (1952а), Андреева и др. (1969).

139

Наличие слоя молекулярной теплопроводности в поверхностном слое воды позволяет определить общий поток тепла на границе раз­ дела вода—атмосфера. Однако реализация такого метода для ши­ рокого использования крайне затруднительна, так как существуют лишь единичные экземпляры двухканальных радиометров. Кроме работ Мак-Алистера, известна лишь работа Красавцева и Сивкова (1972). Поэтому представляет интерес рассмотрение интенсивно­ сти теплопередачи в слоях воды несколько большей вертикальной протяженности. С одной стороны, при этом упрощается техника измерений (необходимо измерение температуры поверхности с по­ мощью инфракрасного радиометра и измерение температуры воды на некоторой глубине). С другой стороны, как показано выше (см. рис. 5.1), перепад температур на границах пленки еще не дает воз­ можности непосредственной увязки величин 0Юи во.

Такие данные для верхнего двухсантиметрового слоя воды при­ ведены в работе Малевского-Малевича (1969). В этой работе ис­ пользованы данные Мак-Алистера (1964), Альтберга и Попова (1934), Тимофеева и Малевского-Малевича (1967).

Альтберг и Попов (1934) 1 провели измерения на р. Неве при помощи миниатюрных термосопротивлений. Результаты их измере­ ний приведены в табл. 5.3.

 

 

 

Таблица 5.3

Глубина

 

Ьщ ° с

 

 

 

 

(мм)

1 5 / X I

 

1 7 / X I

 

 

0

0,38

 

-0 ,2 0

2

.0,85

 

-0 ,1 1

4

1,00

от

-0 ,0 9

6

1,00

—0,09 до 0,00

8

от

0,00 до 0,02

10

от

0,00 до 0,02

Шуляковский (1960) обработал эти данные с целью расчета по­ токов тепла и определения коэффициента теплопроводности для глубин 2 и 4 мм. Тимофеевым и Малевским-Малевичем (1967) при­ ведены данные о коэффициенте теплопроводности для глубин 1 и 2 см (по наблюдениям на Валдайском озере). Результаты ряда ав­ торов приведены на рис. 5.2. Из рисунка следует, что изменение коэффициента теплопроводности с глубиной в верхнем слое воды происходит по закону, близкому к линейному. Это обстоятельство в значительной степени упрощает расчеты потока тепла в воде, так как величина А/Я (или Я/А) может быть принята постоянной в пределах верхнего слоя.

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

1 Очевидно, это одна из первых работ с подробными измерениями про филя температуры в поверхностном слое воды в естественных условиях.

140

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