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

книги из ГПНТБ / Кузьмин, А. Д. Физика планеты Венера

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

230 Гл. V. ОПТИЧЕСКИЕ СВОЙСТВА И ТЕПЛОВОЙ РЕЖИМ

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

Малкус [579], рассмотревший трехмерную модель, при­ шел к выводу о возможности возникновения вблизи вер­ шины тропосферы слабых ячеек гадлеевского типа и пере­ носа за счет температурных контрастов между экватором и полюсом энергии и углового момента к высоким широ­ там. Сэтой моделью можно, в принципе, совместить экспе­ риментальные данные о пространственных движениях, полученные «Маринером-10», если дополнительно пред­ полагать большое влияние подсолнечной области высокого давления на структуру циркуляции. Возникающие в этом случае меридиональные градиенты давления должны приводить к ускорению с ростом широты потоков, орто­ гональных изобарическим поверхностям, с образованием спиральных движений и струйных течений. Сток энергии будет происходить в полярных вихрях вследствие схожде­ ния векторов скоростей движений у полюсов при низких давлениях [581].

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

Вертикальный тепло- и массообмен на вращающейся планете описывается нестационарным уравнением теплового

*) Состояние проблемы общей

циркуляции

в

атмосферах пла­

чет подробно рассмотрено в работах

Минца [110,

429],

в моногра­

фии Голицына

[46]. Применительно

к

Венере режимы

циркуля­

ции рассматриваются также

в монографии Бурангулова,

Бйлитип-

кевича и др. [554]. С учетом

этого

мы

уделим

основное

внимание

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

в

рамках моделей ячейковой

и проникающей

конвекции (см. §

V.5).

Их следует

при этом рас­

сматривать скорее как механизмы,

которые м о г у т

иметь место на

Венере, а не как отражающие

реальную

ситуацию в венерианской

атмосфере,

 

 

 

 

 

 

 

 

 

§ V.2. ОБ ЭНЕРГЕТИЧЕСКОМ БАЛАНСЕ В АТМОСФЕРЕ 231

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

p C p ( ' S “ + 7 g r a d 7 ’) = = ? + g®-

(v -9)

Выравнивание поля температуры в горизонтальном на­ правлении при этом не учитывается (в (V.8) дц = 0). Тепловой поток q = div (A, grad Т) может быть определен из решения уравнения лучистого переноса. В (V.4) можно положить

д® (Л <

Л») =

0

и д® =

/ (h) HaTtdShoo,

 

 

 

 

Q

где / (h) — доля

солнечной энергии, поглощаемой ниже

уровня h [462].

член

в

левой

части (V.9) формально

Конвективный

представляется выражением для свободной конвекции:

Як — рСрУДГ — рСркн + Taj i (V.10)

где А Т — разность между установившейся температурой

вконвективной атмосфере при скорости V и температурой

сградиентом, равным сухоадиабатическому. уа опреде­ ляется формулой (IV. 12), а коэффициент турбулентной диффузии [313]

 

 

 

 

дТ

ч*

 

 

 

кн — оL2

dh + Та

 

дТ .

. А

 

 

 

прн -аГ +

Т а< 0 .

Здесь

L — масштаб

движений,

о — числовая константа

(о ~ 1,4).

Очевидно,

что кн = 0

в случае,

когда выпол­

няется

критерий устойчивости *):

 

 

 

 

 

■ § • + * ■ > о-

<V 1 I >

*)

Более

общим является

критерий устойчивости Пристли

дТ

Та >

Т

cicaa2v

i выражающий тот факт, что помимо соот-

ИГ '+

g

^

 

 

 

° дТ

 

 

 

ношения между Та и

устойчивость зависит от характерных раз-

232 гл. V. ОПТИЧЕСКИЕ СВОЙСТВА И ТЕПЛОВОЙ РЕЖИМ

Эти соотношения выражают условие переноса тепла за счет свободной конвекции в атмосфере вверх. В качестве характерного масштаба для глубокой атмосферы обычно принимают L ~ Н. Заметим, что соотношение (V.10) соот­ ветствует диффузионному уравнению в форме Фикка и описывает перенос тепла вихрями, размеры которых малы по сравнению с L.

Если вместо Т ввести потенциальную температуру 0 = Т (Pq/P)ri^cp, то (У.10) перепишется в виде

где за счет эмпирической константы ^названной Дирдорфом [264] «противоградиентом», обеспечивается при неиз­ менной энтропии выравнивание потенциальной температу­ ры при возникновении отклонения от гидростатического равновесия. Эта константа, по существу, выражает факт незначительного отличия dTIdh от уа (малость АТ7) в не­ устойчивой конвективной области, где, согласно Гирашу

dT/dh Т-

и Гуди [313], выполняется условие -------------< tjl.

Та К сожалению, подобный формализм неприменим к гори­

зонтальному тепло- и массообмену и к переносу момента количества движения. Это следует из изучения метео­ рологических процессов в земной атмосфере, где турбу­ лентный перенос связан с такими крупномасштабными возмущениями, как циклоны и антициклоны, и удовле­ творительной параметризации динамической структуры горизонтальных полей не найдено [325].

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

меров элементарного конвективного объема /0 и интенсивности тур­ булентного теплообмена с окружающей атмосферой (а, у — коэффи­ циенты температуропроводности и вязкости, с1жс2 — константы, за­ висящие от формы частиц). Отсюда следует, что с уменьшением /о растут потребные для конвекции вертикальные температурные градиенты из-за быстрой диссипации энергии всплывного движения малых объемов в процессах турбулентного перемешивания. Тео­ рия Пристли применима вне пограничного слоя (подробнее см. [174, 263, 466]).

§ V.2. ОБ ЭНЕРГЕТИЧЕСКОМ БАЛАНСЕ В АТМОСФЕРЕ

233

ренных параметров атмосферы Венеры разница между температурами на дневной и ночной сторонах должна быть крайне мала, что отмечалось уже в [117]. Действительно, запас тепла в столбе газа единичного поперечного се­ чения

' 9 =

UСр(Т) dT = А к

Т (h) п (.h)dh

 

о

 

составляет tj

4,2-1011 эрг-см 2 ~

104ккал-см 2. В то же

время потеря тепла за венерианскую ночь (tH = 58,4 зем­

ных суток): аТе tH~ 1012 эрг-смГ2 ~ 24 ккал-см~2, т. е.

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

температуры в

нижней атмосфере Венеры при

q% —

= q<? (1 — А )/4

TI2

1°К*).

не должны превышать -=—— ~

 

i 9 lg%

 

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

mCpT

(V.12 )

~атГ

характеризующее временную шкалу выравнивания темпе­ ратурных неоднородностей по отношению к возмущениям путем излучения [573, 574], и время радиационного теп­ лообмена

Тр.т. = CpT/q (Я)

(V.13)

—характерная временная шкала для переноса вверх энер­ гии, обеспечивающей компенсацию радиационных потерь на верхней границе атмосферы [398]. Здесь тп — масса единичного столба атмосферы, q (Я) — скорость охлаж­ дения на единицу массы с уровня порядка Я, т. е. соот­ ветствующего ~ V3m.

*) Очевидно, аналогичная оценка среднего суточного изме­ нения температуры в атмосфере дается формулой А Т =

g<sT\t

= . -----— , где Рь — давление у поверхности, t — время в сутках.

P S C p

234 Гл. V. ОПТИЧЕСКИЕ СВОЙСТВА И ТЕПЛОВОЙ РЕЖИМ

Если для углекислой атмосферы принять эффективность излучения е ~ 0,15, то

1 р . Т . ~

тСрт —. 2-t.

(Y.13')

т . р »

0,5зГр

Значения тр.т, и тт.Р( для Венеры приведены в табл. 17. Для сравнения там же показаны соответствующие оценки для Земли и Марса согласно Лиови и Поллаку [3981.

 

 

 

 

 

Т а б л и ц а

17

Характерные температуры и времена тепловой релаксации

 

П ланета

т С Г 1

Т е, °К

т, ° к

Т Т . р , с е к

Тр т

сек

 

ЬР\_мг-град

 

 

 

 

 

 

Венера

10е

230

~500

Ю Н )

2■1010

Земля

107

250

~

Г .

107

2-107

Марс

1,3 -105

200

~

Т '

3-105

6-10°

Интересно, что если бы Венера перестала получать солнечную энергию, ее остывание за счет потери тепла излучением до теплосодержания земной атмосферы

(г© ~ 24,5 ккал-см~2), соответствующего ночной потере тепла на Венере, при условии сохранения всех остальных параметров неизменными, произошло бы приблизительно за 65 земных лет. На самом деле эта оценка является очень грубой вследствие неучета изменения непрозрачно­ сти атмосферы при остывании, режима планетарной цир­ куляции, а также нарушения равновесных условий форми­ рования газовой оболочки Венеры, прежде всего равновесия между карбонатами и силикатами (см. § IV.8).

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

§ V.3. ОПТИЧЕСКИЕ СВОЙСТВА АТМОСФЕРЫ ВЕНЕРЫ 235

нению с временем высвечивания запасенной энергии, рас­ четную схему можно строить, исходя из того, что падаю­ щая на дневную сторону солнечная энергия на широте ср поглощается на миделе усеченного конуса высотой Лф и осредняется по всей его боковой поверхности на этой широте (что справедливо для вращающейся планеты). Ин­ тегрирование левой части соотношения (V.7) в такой расчетной схеме по Q дает следующий закон равномер­ ного распределения солнечной энергии по поверхности планеты на уровне hx (см. [5]):

<?9 (! — А)

(V.14)

q = - 2— ----- cos ф,

что приводит к Те, отличающейся от (IV.1). Формально это связано с неучетом вклада дц в энергетический баланс

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

§V.3. Оптические свойства атмосферы Венеры

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

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

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

236 Гл. V. ОПТИЧЕСКИЕ СВОЙСТВА И ТЕПЛОВОЙ РЕЖИМ

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

лекулярного рассеяния erf, селективного поглощения отдельными компонентами газовой фазы хх, поглощения

ах и рассеяния of на аэрозоле. Соответственно

Ах = 4 + хх + ах + Gf .

(V.15)

Коэффициент экстинкции и входящие в него слагаемые определяют оптическую толщину:

h

 

xk(h) = ^/«х (h) dh

(V.16)

о

 

и альбедо однократного рассеяния, характеризующее по­ глощательные свойства рассеивающей среды:

“ ох -= —

+

----,

пт лп\

 

(V.17)

которые вместе с индикатрисой рассеяния (ср) представ­ ляют собой основные оптические характеристики атмо­ сферы.

Рассеяние на компонентах газовой фазы наиболее су­ щественно в коротковолновой области спектра вне полос поглощения, а в длинноволновой преобладает аэрозоль­ ное рассеяние. Подробное изложение теории рассеяния света на малых частицах читатель может найти в монографияхван деХюлста[21], Соболева [146,149], Шифрина [173], Розенберга [138] (см. также [109, 137, 531]), а на круп­ ных частицах — в работах Ирвина [353, 357]. Из теории следует, что в силу некогерентности рассеяния газовой фа­ зой и частицами аэрозоля можно допустить аддитивность соответствующих компонент матрицы рассеяния атмосфе-

ры Dik (ц, ц0) для падающего и рассеянного под углами arccos ц0 и arccos ц потоков лучистой энергии. Интегри­ рование по углам соотношения

ХцШ/РШ

§ V.3. ОПТИЧЕСКИЕ СВОЙСТВА АТМОСФЕРЫ ВЕНЕРЫ 237

с учетом Dik = ([х, ц0) позволяет выразить эффективную

нормированную матрицу рассеяния для мутной среды в виде (см., например [21,138])

/а- (й) М-о)

° g f ik 0х- М + ° afik (Р> Р»)

(V.19)

+ <з“

 

 

Конкретный вид матрицы рассеяния (число независи­ мых и отличных от нуля компонент) определяется свой­ ствами рассеивающей среды (составом, размерами, фор­ мой частиц) и зависит от длины волны.

Молекулярное рассея­ ние, происходящее по за­ кону Рэлея, описывается коэффициентом рассеяния 0й и индикатрисой рас­ сеяния х (ф), которые вы­ ражаются хороню извест­ ными формулами:

ае

32я3 — 1)2

(Y.20)

 

 

3NX*

 

 

 

 

 

 

Жр(ф)=4 1

cos2 q>), (V.21)

 

 

где

т — средний

показа­

Рис. 79. Изменение оптической толщи­

ны для рэлеевского рассеяния

в ат­

тель

преломления

среды,

мосфере Венеры в зависимости

от X.

N — число

рассеивающих

 

 

молекул в единице объема,

 

 

Ф =

arccos р—угол рассеяния. Из формулы (V.20) видно,

что роль молекулярного рассеяния особенно велика в переносе видимой солнечной радиации. Результаты расче­ тов тр в интервале XX от 0,2 до 1,2 мкм для атмосферы Венеры (согласно приведенной в Приложении основной модели) показаны на рис. 79 и 80. На рис. 79 приведена зависимость тр (в нормировке на единицу давления) от X, а рис. 80 дает семейство кривых хр (X) по высоте. Посколь­ ку рэлеевское рассеяние прямо пропорционально массе рассеивающего газа, на Венере оно оказывается прибли­ зительно на два порядка интенсивнее земного. Вблизи максимума кривой видности (Х0 ~ 0,5 мкм) тр ~ 30, что

238 Гл. V. о п т и ч е с к и е Св о й с т в а и т е п л о в о й р е ж и м

примерно соответствует оптической толщине земных об­ лаков. Заметим, что спектральная зависимость ag выражается удобным отношением, очевидным из (V.20):

of = a?0(^-')4.

(V.20')

Рассчитанные значения

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

рассеяния сГх0(К)

и оптической толщины

тр\0 в зависимости от высоты в

атмосфере Венеры приведены в табл. 18.

Ь,км

Рис. 80. Семейство кривых (К) для атмосферы Венеры в зависимости от высоты.

При наличии взвешенной в газе аэрозольной состав­ ляющей оптические свойства атмосферы дополнительно зависят от концентрации, высотного распределения, рас­ пределения по размерам, природы и микроструктуры час­ тиц. В этом случае все 16 компонент матрицы рассея-

ния fih (р, р0) могут быть, в принципе, различными. При характерных размерах частиц г К удобное приближе­ ние дается теорией Ми, позволяющей рассчитывать сга в зависимости от показателя премломления т и пара­ метра р = 2я/7А,. С увеличением содержания аэрозоля и ростом этого параметра возрастает дисимметрия фазо­ вой функции (вытянутость индикатрисы рассеяния в

§ V.3. ОПТИЧЕСКИЕ СВОЙСТВА АТМОСФЕРЫ

ВЕНЕРЫ 239

Т а б л и ц а

1 8

Коэффициент рэлеевского рассеяния

(Л)

 

(А/0 = 0,5 мнм) и соответствующая оптическая толщина

Тр°

в зависимости от высоты (согласно [37])

 

h ,

 

 

h,

аХ0.1(Г

 

н ,

° v i o ?

к м

с м —1

Р

к м

с м ~ 1

Р

к м

с м ~ 1

• а о

0

165

2 8 , 4

24

44

5 , 6

48

6 , 2

0 , 5

3

144

2 3 , 8

27

37

4 , 4

51

4 , 9

0 , 3

6

124

1 9 ,7

30

29

3 , 4

54

3 , 6

0 , 2

9

106

1 6 ,3

33

24

2 , 6

57

2 , 5

0 , 1

12

90

1 3 ,4

36

19

1 , 9

60

1 , 6

0 , 0 3

15

76

1 0 ,9

39

16

1 , 4

70

0 , 2 8

0 ,0 1

18

64

8 , 8

42

11

1 , 0

80

0 , 0 4

0 ,0 0 3

21

53

7 , 0

45

8 , 7

0 , 7

 

 

 

переднюю полусферу) за счет определяющего влияния

аэрозольной матрицы

рассеяния

(р, р0); одновременно

возрастают

и ах по сравнению с o'? и хх.

 

С учетом полидисперсности и неоднородности состава

аэрозоля, второе

слагаемое в (V.18) можно

представить

в виде [138]:

 

 

 

 

 

 

 

 

оо

 

 

£>ift(P,Po)

= 2

\ я г/(г)

 

(V.22)

vо

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

чить выражение для эффективного сечения рассеяния:

оо

° а = 2

$ п У ^

3?У d r ■

(V -23)

v

о

 

 

Плотность распределения по

размерам частиц

сорта у,

пу (г) = dNv (г)/dr дается так называемым F-распределе- нием п (г) = n0e~irru, либо зависимостью типа Юнге [368, 568] при различных значениях параметра степенного рщ>

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