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

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

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

V. ПРИЗЕМНЫЙ СЛОЙ АТМОСФЕРЫ

Приземным слоем называется сравнительно тонкий слой у поверхности земли (толщиной 10—100 м), в котором наблюда­ ются, как правило, максимальные градиенты метеорологических элементов (по мере приближения к поверхности земли градиен­ ты возрастают как \/z). Важной особенностью приземного слоя, в значительной мере оправдывающей выделение в рамках по­ граничного слоя, является постоянство по высоте турбулентных потоков.

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

§ I. Формулировка задачи о строении приземного слоя

Проинтегрируем от г0 до z уравнения движения для стацио­ нарного, горизонтально-однородного пограничного слоя атмо­ сферы (4.4.1) для случая, когда ось х совпадает с направлением приземного ветра

к

du

,

du

Z

2u)z (v—G sin

a) dz — 0,

 

dz

к

-гг-

-f I

 

 

 

dz

*0

 

 

 

 

 

 

 

 

 

 

 

(5.1.1)

 

dv

 

dv

 

 

 

 

 

к

 

— j 2и>г (u —G cos a) dz =

0.

 

dz

 

dz

 

 

 

 

*0

zo

 

 

 

 

 

 

 

 

 

 

du

 

dv

 

Тогда на

основании (4.4.8)

k dz

-V

а k

dz

= 0 , так как

ось X направлена по приземному ветру и совпадает с направле­ нием касательного напряжения у поверхности земли т (рис. 30). Действительно, если « — угол между направлением приземного ветра с и осью х, а а' — угол между направлением касательного напряжения т и осью х, то

80

tg а

 

V

 

и

 

-»0

 

г - * 0

tg*'

_

~і/

 

 

dv'dz du dz

Z->

dvtdz du.dz

2-+0

2-*0

 

 

z->0

 

 

>■ Г

Перепишем систему (5.1.1)

в виде

 

Рис. 30. Соотношение между

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

касательного

напряжения

 

.

du

I

 

 

 

 

 

 

k- -у-

I

,

, 2шг

г

 

 

 

 

dz

г

 

 

 

 

Т>*2

— *

i

J (■о—G-sina) dz,

 

 

 

dv

 

 

 

 

 

(5.1.2)

 

 

 

 

 

 

 

 

 

 

dz

г

2o>z

г

 

 

 

 

 

v j

~

V»2 I' {u—G' cos а) dz.

 

Введем средние для слоя от z0 до z величины и, ѵ и а и пе­

рейдем в правой части (5.1.2)

к

безразмерным

переменным:

un = u!G,

vn—vjQ.

z„= 2 ,7г. Тогда

 

 

 

,

du

 

 

 

 

 

 

 

к

dz

 

1

2ш7■h -ü

-

 

 

 

v J

 

V*

(T'n- s in a )-z„;

 

 

 

 

 

 

 

 

 

dv

 

 

 

 

 

 

 

dz

2щ -h-G (un- cos а)-zn.

 

 

 

Vf.

 

 

 

 

v*‘

 

 

 

Так как по определению (см. раздел II § 6) все безразмер­ ные величины имеют порядок единицы и, кроме того, cos а = 0(1),

а s in « < < l (по наблюдениям а = 1 0 —15°), то

,

du

 

 

к- -і~

 

 

 

dz

2

uv а-А 0 ( 1),

 

 

,

dv

 

(5.1.3)

 

 

Т = - - Й ^ [ 0 ( 1 ) - 0 ( 1 ) Ь

6

S1

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

практически постоянны^

по высоте, т. е. высоту

приземного'

слоя; она должна удовлетворять условию

 

 

2i»g-Q‘h

 

 

•X1, 2

 

 

Если 2иг —1C“4, G = 103 см/сек, а ;;2=103 см2/сек2, тогда

 

//

100 и/.

(5.1.4)

Из второго уравнения

(5.1.3)

следует, что для Іг < . 100 м

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

и ііи

•>

,

* • 77; =

 

 

(°л -0)

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

Проинтегрируем теперь от го до г уравнение притока тепла

(4.1.8), тогда

d » \

Л ä O \ _

/? (г)

R (0)

d z \

1

d z l ~

^

где R (о) —R(Zo).

Согласно определению турбулентного потока тепла

о _

.

аѳ

!

,

‘ о

Р'сp-jr.T-h •

,

I

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

можно записать

R(z)

R(0) ;

■k-

 

 

dz

R ( 0 )

ОlPCp

Po ( R n - О.

82

п о с к о л ь к у на

основании уравнения теплового баланса

 

 

-

-

 

I

= 0(1), а в пределах

/і<

100 м R „~> 1,0, то (R„ — I)

.1

п

 

 

 

 

 

(5.1.6)

Уравнение

(5.1.6)

показывает, что в пределах

приземного

слоя можно пренебречь радиационным притоком тепла (впро­ чем, как и фазовым). Следует обратить внимание на то, что для условий суши действительно Pq~ R (о ) , так как R(o) = Р + LE + В, а затраты тепла на испарение LE и теплообмен поверхности с нижележащими слоями В малы. Для морских условий турбу­ лентный поток тепла Pq может быть значительно меньше радиа­ ционного баланса поверхности и тогда лучистый приток тепла будет играть важную роль в формировании распределения тем­ пературы в приземном слое.

По аналогии с (5.1.6) можно сразу записать уравнение для определения удельной влажности

(5.1.7)

k

где Е0— скорость испарения; a q— — ' ; kq — коэффициент тур­

булентности для потока влаги. Уравнение (5.1.7) предполагает, что в приземном слое можно пренебречь притоком влаги за счет фазовых переходов.

Система уравнений (51.5), (5.1.6), (5.1.7) содержит четыре неизвестные величины: //, Ѳ, q, k. Для замыкания системы ис­ пользуем уравнение баланса энергии турбулентности (2.4.20), соотношения, полученные па основании теории подобия (2.4.17, 2.4.18), и выражение для пути перемешивания (2.4.19"'). Для малых 2 эти уравнения примут следующий вид:

(5.1.8)

k

k = i V ь-

(5.1.9)

(5.1.10)

( 5 .1 .П )

83

§ 2. Замыкание системы уравнений для приземного слоя атмо­ сферы на основании теории подобия

Качественный

анализ системы

(5.1.8) —(5.1.11)

показывает,

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

турбулентности

должен

однозначно

опреде-

ляться следующими параметрами

[12]: с*.

СГ

р

 

 

е

Z, ,

— — и г, т

 

 

"

7

Ѵ'СР

'

1

 

А = /

 

 

 

 

(5.2.1)

Раскроем функциональную зависимость

(5.2.1)

с

помощью

П-теоремы, для чего выпишем формулы размерности всех вхо­

дящих в

(5.2.1) величин: [k]— L2- 7_1, [тт*]— L ■

= І Х

 

Г

X 7’-2-

рс„- = 1 - гр -Г ^, \z]^L .

 

Гак как в данном случае ищется связь пяти физических ве­ личин, из которых три имеют независимую размерность, то функциональную зависимость (5.2.1) можно раскрыть только с точностью до функции от безразмерного аргумента. Нетрудно показать, что из величин, входящих под знак'функцпи в (5.2.1), единственным образом с точностью до числового множителя составляется безразмерная комбинация zfL, где L — имеет раз­ мерность длины и называется масштабом длины Монина—Обу­ хова. Действительно, на основании П-теоремы

g

J \ \

k

• Z'1 — F

Приравняв показатели степени при независимых размерно­ стях, определим а и <х2, ßi и ß2:

7 '

I 2 CCi + ß i ,

2~CC2~hß2

Т

I —1 ==—ОЫ,

—3 = —<х2

«i = l, ß - 1, а 2 = 3, ß2 = —1.

В таком случае

k = v * - z - F ,

\ g Т-Р, -л:р

или, обозначив через

84

, , ё '

&

(5.2.2)

 

' T

f,Cp

 

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

А—г ѵ 2- ^ -f-)•

(5.2.3)

В (5.2.2) постоянная Кармана у. введена для удобства неко­

торых последующих выкладок, а знак выбран

так, чтобы было

L> 0 при устойчивой стратификации.

Покажем, что zjL можно рассматривать как критерий устой­ чивости; для этого установим связь между z/L и используемым обычно для описания устойчивости числом Ричардсона. Число Ричардсона характеризует отношение энергии турбулентности, возникающей за счет сил плавучести, к продукции энергии тур­ булентности за счет среднего движения. В соответствии с этим" градиентная форма числа Ричардсона имеет следующий вид:

Если предположить, что профили температуры и скорости ветра подобны до самой подстилающей поверхности, то из (5.2.4) нетрудно получить выражение для приближенного (ин­ тегрального) числа Ричардсона

tw

Ш>= Г '

где иа — скорость ветра на г —а. При z 6 м, Г = 300°, ссх =1,0, получаем,что

 

 

 

Rb—0,2-

L - L

 

 

(5.2.5)

 

 

 

 

и /

 

 

 

Если ввести так называемые профильные коэффициенты

JL

du

_

1

dt

^ _

1

dg

иа

' d\n z'

т~

ta — tw

‘ сПпз’

ч~~ qaqw " d \ n z

и учесть,

что согласно наблюдениям до высот

4—8 м и страти-

_ фикации от безразличной до умеренной

(над морем)

ГЫ= ГТ — -

= Г д —0,1,

то для

выбранных величин z и Т получим

 

85

(5.2.6)

Наряду с градиентной формой числа Ричардсона часто ис­ пользуется число Ричардсона, выраженное через потоки

(5.2.7)

dz

Если выразить поток количества движения т и градиент ско­ рости через соотношения

дс _ г»*2

то нетрудно установить связь между Rf и г/L

(5.2.8)

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

Хотя соотношение (5.2.3) позволяет замкнуть систему урав­ нений для приземного слоя, им трудно воспользоваться, если не раскрыть вид функции F(z/L).

§ 3. Частные выражения для профилей коэффициента турбу­ лентности и метеорологических элементов

Для некоторых предельных случаев стратификации удается получить-явный вид функции F(z/L).

Безразличная стратификация

Так как при этом равен нулю поток тепла, то L—ooyi z/L = 0. Измерения в аэродинамических трубах и в приземном слое по­ казывают, что F(0) » к » 0,40 и при г > г0

86

k = -/.-vä.-z.

(5.3.1)

Из (5.2.8) следуех, что при стратификации, стремящейся к безразличной, z j L Rf.

Выражения для профилей метеорологических элементов мож­ но получить, если проинтегрировать уравнения (5.1.5), (5.1.7) от zQд о 2

и =

~ ln z,z0

 

(5.3.2)

 

ln z ‘za

 

(5.3.3)

и — lit - 1In Zt za

;

 

E

.

,

(5.3.4)

q~q0= --------------In z;z0

 

PV'-'1'*

 

 

,

\

In z/z0

(5.3.5)

<7 «.=(9=

9.)

 

 

где c/o— удельная влажность

на уровне г0; иi,

q\, <?2 — скорость

ветра и удельная влажность-.на 2 = 2 | и 2= 22. Для морских усло­ вий можно считать, что c/o = qw —- насыщающей удельной влаж­ ности при температуре поверхности i w.

Итак, в случае безразличной стратификации ■коэффициент турбулентности является линейной функцией высоты, а профили метеорологических элементов — логарифмической функцией вы­ соты (рис. 31, 32).

Продифференцировав (5.3.3) и (5.3.5) по z и использовав (5.3.1), можно получить выражения для определения турбулент­ ных потоков

а,

— У ..Х )

. р . . - --------1

 

 

 

\n z Lz0

(5.3.6)

Г

 

 

—а

Я\ —Яі

Е

• р ■ у.•V * • --- ---- —

 

4

111

Z ; Zj

Стратификация, близкая к безразличной

Б этом случае мал турбулентный поток тепла (или мало г)., а F(z/L) как функцию от малого аргумента можно разложить в ряд Маклореиа

87

f

Рис. 31. Профиль скорости

 

Рис.

32.

Профиль

удельной

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

 

влажности

в

приземном

слое

безразличной

стратификации

 

для

безразличной стратифика­

 

 

 

 

 

 

 

 

ции

 

 

 

f (Z;d = f (о ) - ~

F ' ( 0 ) + ( ~ y .

 

+ . . .

 

Ограничимся первыми двумя членами ряда,

так как

(z/L)2 уже

величина второго порядка малости. Тогда

 

 

 

 

 

 

F(z/L)=F(0).

Р'

2

 

 

 

(5.3.7)

 

F'(Ö)

L

 

 

 

где ß=

 

 

 

 

 

 

 

 

 

F(0) (величина ß, по современным наблюдениям,

нахо­

дится в пределах от 3 до 10)

и (5.2.3) запишется в виде

 

 

 

 

: У. • V.,, •Z •

1 -ß

2

 

 

 

(5.3.8)

Из

(5.3.8)

следует,

что

при

устойчивой

стратификации

(F>0) коэффициент турбулентности меньше, а при неустойчи­ вой стратификации (L < 0) — больше, чем при безразличной.

Проинтегрировав (5.1.5) — (5.1.7) от го до г, с учетом (5.3.8) получим выражения для профилей метеорологических элементов

33

du

VJ

 

V#

_

“D*

 

dz

k

7.

2 ( 1 — 3 2

L )

7.2

 

откуда

 

 

 

 

 

 

 

 

 

и =

(ln

— +^-7-

(5.3.9)

 

 

 

 

 

L

 

(считая, что

<

L Y

 

 

 

 

 

 

Ѳ — ö 0-j- Г* (ln — + ? y ) )

(5.3.10)

 

 

 

 

 

 

 

 

 

+

(/*

(ln — -Т.З-Г-),

 

 

 

 

 

+>f\

 

L, i

 

где

Г*: ^ -У *ят ?*=

Итак, в случае стратификации, близкой к безразличной, полу­ чаются линейно-логарифмические профили метеорологических элементов (рис. 33, 34). Очевидно, что при устойчивой стратифи-

89

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