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

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

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

р=

'

k-ii- ат

(6 .4 .2 7 )

— ѳ і)

r.-x

 

ф

 

Р ^ - с А ^ о - ^ ^ - у

ukZj

(6.4.28)

~ p - A - Y L .

Напомним, что рассмотренное выше решение задачи о транс­ формации было получено для &= const и n = const. Оценки вы­ соты внутреннего пограничного слоя показывают, что для рас­ стояний около нескольких километров влияние подстилающей поверхности сказывается в основном в пределах приземного слоя. Поскольку именно в приземном слое скорость ветра и ко­ эффициент турбулентности заметно меняются с высотой, то же­

лательно отказаться от условия к = const

и м = const и задать

закон изменения их с высотой. Допустим, что

 

 

и

Z \г

 

k ~ k x

■ Л - £

(6.4.29)

 

 

 

 

 

 

 

 

 

 

 

 

н )

 

 

Тогда (6.4.10)

примет вид

 

 

 

 

 

 

ul -zi1—

 

0 &_

=V dz Z

дЬ

 

 

 

дх

dz'

 

Обозначим

постоянный

и

положительный коэффициент

U\Z\ 1— 2s'/kl через а2. В таком случае

 

 

 

 

a--ze

 

д

1

dft

(6.4.30)

 

дх

:ат—

2

dz’

 

 

 

dz

 

 

Полученное уравнение можно свести к обыкновенному диф­ ференциальному уравнению, если, как и раньше, ввести новую переменную g = £(х; z ) Вид этой зависимости найдем на основа­ нии теории подобия. Для всех переменных величин в (6.4.30) введем характерные масштабы

z=H -z',

-L • х ' ,

» '=

ft— ft,

 

(6.4.31)

f t 0- f t

1

 

 

 

 

и перепишем уравнение в безразмерном виде

a*-tf1+2s

dti'

О

(6.4.32)

L

' Z ' d x ' ~ ^ ' d z ' Z

~dzr’

Так как масштабы H и L можно выбрать какими угодно, то без потери общности задачи возмем их такими, чтобы

9

129

 

<f- ■H l + 2 г

 

 

 

L

 

 

 

T. e. L = a2-Д 1 ; 2г.

примет вид

'г"'л

С учетом этого (6.4.32)

 

 

_д_

,1-

дѴ

z Т х ' ~ У'т' dz'

z

'

dz''

Запишем граничные условия для d':

 

 

 

-О,

1 ,

г

 

х'—Ч

*'■--()

 

z'==оо

z ' X J

а-'; 0

 

х '> 0

Если решить (6.4.33) с учетом (6.4.34), то

d '= /(*', z')

или

д — ді = (до — ді) /

яг-/У

(6.4.33)

(6.4.34)

(6.4.35)

С другой стороны, если в (6.4.30) перейти к новой перемен­ ной д' и воспользоваться граничными условиями (6.4.34), то ре­ шение будет выглядеть так:

д' = /(*, г)

или

d - d i = ( d o - d ,) / ( J f , 2 , Й2).

(6.4,36)

Из (6.4.36) видно, что масштаб Н не входит в число парамет­ ров, определяющих' д, и значит в (6.4.35) под знаком функции Н должно сократиться, т. е.

12 г

»--»і=(»о—»і)*/( а-■/ г

2 і ’ н

2 s

(6.4.37)

Итак, при решении задачи целесообразно ввести новую пере­ менную

 

X

(6.4.38)

а

Tf+Ys’

At 1

 

которая позволит свести исходное уравнение к обыкновенному дифференциальному уравнению. Дальнейшее решение при гра-

130

ничных условиях (6.4.11) — (6.4.13) производится так же, как и в предыдущем случае, поэтому запишем сразу окончательное выражение для іУ

(6.4.39)

* 0 - .

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

нет 100%. Так как г = е/£

— упругость

водяного пара;

Е — насыщающая упругость

при

температуре воздуха), т. е.

относительная влажность зависит

от фактического влагосодер-

жания

и температуры воздуха, то ясно, что г

может достигнуть

1 0 0 %

как за счет повышения

влажности, так и за счет пониже­

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

Полученные ранее формулы для определения турбулентного потока тепла и затрат тепла на испарение могут быть использо­ ваны при прогнозе сроков замерзания водоемов. Для этого не­ обходимо знать теплозапас слоя воды в начальный момент и скорость уменьшения этого запаса из-за потерь тепла через по­ верхность.

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

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

\

-131

§ 5. Трансформация динамических характеристик воздушного потока при изменении шероховатости подстилающей поверхности

Ограничимся рассмотрением случая нейтральной стратифи­ кации и не очень больших расстояний от границы раздела, при которых можно считать, что внутренний пограничный слой не выходит за рамки приземного слоя. Будем считать, что ось .ѵ совпадает с направлением ветра, при д'<() расположена поверх­ ность с z0 = z'о, а при х > 0 — поверхность с z0 = z"0 (рис. .48), ось у направлена вдоль границы раздела.

z

Рис. 48. Качественная картина трансформации воздушного потока при резком изменении шеро­ ховатости

Для указанных условий система уравнений (6.4.1) -- (6.4.9) примет вид

 

ди

,

ди

д

ди

 

 

(6.5.1)

 

ох

 

 

dz

dz

dz

 

 

 

 

ди

,

dw

n-

 

 

 

(6.5.2)

 

'

дх

1

dz

 

 

 

 

 

 

 

 

 

 

db

,

 

db

 

, (du\•

 

b2

 

(6.5.3)

 

+

 

=

-

CT

 

 

 

k =

H

‘4 .1 * . ^

 

 

 

(6.5.4)*

 

 

 

 

 

d-u dz1

 

 

 

Если считать, что при х <0

существует

установившийся по­

ток воздуха, т. е. можно

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

модель горизонтально­

однородного приземного

слоя,

а

при л: >

0

на

 

верхней границе

у * Выражение (6.5.4) получается из (6.4.7.) —(6.4.9).

132

внутреннего пограничного слоя выполняется условие склейки скоростей и непрерывности потоков количества движения, гра­

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

 

'

-

 

х — 0.

u = —*~\nz z0', /;- с

1;а-г’Ѵ ;

(6.5.5)

-У• -У.*

^

 

 

 

 

 

X

20",

и — и ' = 0,

 

 

 

 

 

0

 

 

 

 

 

(6.5.6)

 

 

да

 

 

 

 

 

 

£ • dz

W;

 

 

 

 

 

;

. ,

,

,

ди

(6.5.7)

z — li,

и ——ln h z0 ,

#•-—= г ѵ ';

X

0

y-

 

 

d z

 

 

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

,

и

w„—

w \

к

Введем

оезразмерные величины:

ип—— ,

— , Ая= т - ,

 

X

__ Z и ^

IIА

 

сСд

К.>

 

 

 

 

 

хп ~ Т ’ 2"“ 77Ья~ Т 1;

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

гд;;;-гх(.щ/ ;ѵ )-

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

иновой подстилающей поверхностью.

Сучетом указанных выше соотношений получаем следую­ щие выражения для масштабов:

- Ц Н :

2 .V

ѵ'л •»0 >

 

У.

 

 

 

 

ѵл

V ■

(6.5.8)

 

 

Б таком случае систему уравнений и граничных условий за­ пишем через безразмерные величины:

133

 

 

 

 

 

ди.

\ п дх„ ^

п dz J

 

dzn

п

d ztt'

дип

, dw„

-0 :

 

 

дхп

dzn

 

 

 

db„

 

 

 

b 2

L . J U ™ »+ w ™ . n \ k (дЦ.)

1n .

 

 

 

 

 

k :

K = - V K

dUn'dZn

 

 

дЧп dzn2

 

 

 

 

Грашічные условия:

(6.5.9)

(6.5.10)

(6.5.11)

(6.5.12)

л-„-0 , un--\nzn, /?„— 1 ;

(6.5.13)

ди

z« ”/W’ un — wn —Z7/-/7, kn-~— (6.5.14) Xrt>U

 

zn=h,v

ип=\п h,„

kn

: 1 .

 

(6.5.15)

где

-*„>0

 

 

ozn

 

 

 

 

V * — V*

 

 

m ■

P =

:

=a.

(6.5.16)

 

v J

 

 

 

 

 

 

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

мл= « яо+а«лі» w n= w n()4 - x z i’n i,

bn= b n^-'J.bnV

(6.5.17)

 

Подставим (6.5.17) в систему

(6.5.9) — (6.5.15), упростим ее

(отбрасывая члены, содержащие

а 2, а 3 . . ., раскладывая УЬ„

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

Решение системы для а0 —1

- Для^ нулевого приближения система уравнений и граничных условий имеет следующий вид:

д и ди0

T z d t e ^ 0 ’

(€.5.18)

134

д а п

д w 0

О х

‘ d z

(6.5.19)

(6.5.20)

 

 

- _ | У .

()< dz .

 

(6.5.21)

 

 

1 0

(T - u - o d z 2 ’

 

 

z —

т .

«0—-Т£'0—0, k

o = m

- p

 

 

\/ о

 

 

 

 

 

 

(6.5.22)

 

 

d

t i 0

 

 

 

 

 

 

 

 

 

 

 

 

 

U 0 -

-

 

 

 

 

 

 

 

«„—In //.

 

 

 

 

(6.5.23)

X

0

 

 

 

 

 

 

 

 

 

получим* k0' Öi-C

 

 

Проинтегрировав

(6.5.18)

по

z,

const, или

с учетом (6.5.22)

дип= р2, откуда на основании (6.5.20)

 

Перепишем (6.5.18)

в виде

* Ь0= ра.

 

 

(6.5.24)

 

 

 

 

 

 

 

(>k0'dzdu.0Jâz + k0~ , ° —О

 

 

 

 

 

 

 

ÖZЛ

 

 

 

и используем (6.5.21),

тогда

 

 

 

 

 

 

 

 

dkv dz~ (

Ь„~р.

 

 

 

Проинтегрировав полученное выражение с учетом (6.5.22),

получим

 

 

ко~рг.

 

 

 

(6.5.25)

 

 

 

 

 

 

Подставив k0 в соотношение k0

ÖU

— р 2 , найдем

с

учетом

~

(6.5.22) выражение для профиля ветра

 

 

 

 

 

Uo — p\nz!m.

 

 

(6.5.26)

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

(6.5.19)

и граничного

условия

(6.5.22) определим

 

^dp

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.5.27)

■tt,o= —'^ (s -ln z .'/H - s + m).

 

Наконец, используем граничное условие(6.5.23) для опреде­ ления р и dpjdx

135

 

 

_

ln h

t

(6.5

.28)

 

P

ln h 'm

 

 

 

dp _

1

dh

 

ln

(6.5.29)

dx

h

dx

ixPhm

 

 

Из (6.5.19) — (6.5.23)

видно,

что

полученное решение сис­

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

Решение системы для и1 = а

Для этого случая система уравнений и граничных условий

примет такой вид:

 

 

 

 

 

 

 

 

 

 

, ди0__ L

k .дио

 

flux

(6.5.30)

 

^ w°lFz -

 

dz

 

dz

 

 

 

 

 

f)u1

flWy

 

 

 

(6.5.31)

 

 

 

dx п

 

 

 

 

 

 

а -іи dx ^ 0flz }

(du0y

 

ЯиV

Ally

 

1flz j

 

 

flz

 

 

 

 

2bnb

 

b \

 

 

 

(6.5.32)

 

 

 

J L I i J _ f c . ! » •

 

 

 

 

 

 

Kkq

'

v Кk*'о

 

 

 

 

d 7u 0

f l 2U y

, Г ZOilI

 

 

 

djiq

(6.5.33)

«Г - ^ 2

+ * 0 - Q

J

-

 

2

1

bn

flz

 

 

 

 

 

 

Граничные условия:

 

 

 

 

 

 

 

 

z— m, u l= w l= 0 , ky = 0 ,

 

 

 

 

 

x>

0

 

 

 

 

 

(6.5.34)

*• l ? + ^ = 0:

z—h

ux0 , k0 dux-j~k

Alin _

1—P5

(6.5.35)

x > 0

öz

dz

a

 

 

 

Подставив в (6.5.30) — (6.5.35) выражение для k0, b0, uQ, Wo, после некоторых преобразований получим решения для k u u\, bi, которые учитывают влияние горизонтальной и вертикальной адвекции

136

 

 

 

 

 

(6.5.36)

bl==p(ix lF (z)—a -z -lnzlm^

 

(6.5.37)

4 - 7Г 4 фі(г) - ^ Ф 2( 2 ) - — •

h—tn

X

А dx \

1V

w

2

 

X [АФ1(А )-/я Ф2 (А)]}

 

 

(6.5.38)

ГДС

F(z)—z • ln2 zjrn —(Зг-j-m) ln 2'm -f4(z—m).

cl)j (z)=ln2 z'm (6 а) ln z m -[- ^ 1 2 ----

jj- а

Ф2 (г)—ln* z/'m-j-6,0 ln 2 Im -j—М 2---- ^

(6.5.39)

a

ф і ( А ) = ф і (2) I - Фа(А )= ф г(2)

I •

 

z=h

z-^h

 

Определим теперь h — h(x), для чего используем граничное условие (6.5.35), которое с учетом (6.5.28) — (6.5.29) примет вид

1

dh

т

p(u\— 1

n

/ ,

ln2Л

\

, „

h

d x '

ln m

ln Am

/

 

,

(6.0.40)

ln2 Am

'

а

ln А

\

ln* Am

где

 

 

 

F(h)—F(z) I .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 —h

 

 

 

 

Вернемся к обозначению безразмерных переменных значком

п и введем новую безразмерную переменную х п

 

 

 

 

 

Сп

X />

Х п < Х - Х п 'Ѵ-

 

 

(6.5.41)

Тогда (6,5.40) примет вид

ЁЬ5=Я(А„), dx„

где

/?(Ля)=х**Ля. ln* A„/m-f ln h„lm-hnln hn ln hn’F(hn)

137

 

 

7

dhn

 

 

(6.5.42)

 

 

* п == i

 

 

 

 

 

•' R{hny

 

 

 

Таким образом, для определения

h„--~ hn(xn) нужно рассчи­

тать интеграл для разных т и //„.

 

 

 

 

 

Окончательные формулы для определения искомых величин

кп, u'n, Ьп имеют следующий вид:

 

 

 

 

1

ln h„

,

1

dhn

1

hnm

X

"

•lnz« w -

-472

 

nl

ln й„//н

 

 

 

Xh„ ' ln2 hn m

 

dx„

X j*„ ®1 («„)“

“n

Zn ~

К

1

!

e

/4

e

 

 

X [АяФЛА,,)—/пФ*(Ля)] ;

 

 

(6,5.43)

Ая=А-л0 яАв1

 

In2//,

1

dhn

1

ln/?/-ln//

ІП4Л„/Л

у.2

dx„

hn

ln3 h„ mX

 

 

 

 

 

 

 

 

 

 

 

X

\F{zn) - a - z nAn B„ m)\

 

 

(6.5.44)

 

i

 

^ri Ап

1

dhn

z n.

ln m

V

 

llX/x i

іпал>л

X*

' ~

 

ln2 hn m A

 

 

 

 

 

dxn

 

 

 

 

 

- - 4я Фз(0

+ ~

( 2 1 n x ^ : ' - 6 ,0 0 )+

4z,,

. hn

 

hn—m

 

 

X

[ln Ф,(АЙ)—.»та Фа (Ая)]|

 

 

(6.5.45)

где

 

 

 

/

 

 

 

 

 

 

фз(гл)= 1п®zn /?/—(4—a) ln zn"m-\- (б -

 

 

 

wn =* wnl.-faw«!

J_ yHbxnrr, 'Zn\nzn:m-rm-~

 

 

 

 

A„ ’

dx„

1n2z„m

 

 

Ä

СІгГ)

(z^

 

 

 

 

 

 

T

■ rfx2“ X

"2 X i n 2 zn 'm - { 7 - a ) ln z,, /н--(7~д) /2 -

— mzn(ln2 zn' m \ 4 ln z„///i+ 8 -

 

12 —

| - « ) x

 

 

 

 

 

 

X (2 „ - m ) ~ ^ 6 ---- 1- a ) m2 — aß (ân)■

 

3

,

(6.5.46)

 

~2

 

 

 

 

 

2m

 

 

 

138

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