Файл: Океанография и морская метеорология учебник..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 175

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

теорологііческих элемента, чтобы придать составленной системе замкнутый вид. Таким подобием является одно­ значная зависимость между T u t . На практике такой способ оказывается нерациональным. Более приемлемым является метод, основанный на установлении порядка величины метеорологического элемента. Под порядком величины понимается значение этого элемента с обеспе­ ченностью не менее 99%.

Порядок величин метеорологических элементов опре­ деляется с округлением до ближайшей степени числа 1 0 и обозначается 0 (...). Давление воздуха, например,

имеет порядок ()(Я) =

1 0 2, так оно

меняется

в пределах

от 950 до

1050 мбар

(95—105 сбар).

Для

определения

порядков

производных, входящих

в

составленную си­

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

рядок

проекций силы

Кориолиса

 

будет

 

 

0

(2ш sin <рѵ) 0 (2

(о sin щ ) ж

 

1 , 2

1 0 - 4 • 7 • 1 0 ° =

 

 

8 - 10- 4 s=: 10~ 3

(для

ф =

45°),

 

а диссипативных сил

 

 

 

 

 

 

 

 

0 ( д Т Д г ) ^ 0 ( - ^ ^ ) ä

5.10°.1,2-10-6ä 10-5

 

 

 

 

 

 

 

 

 

(при

р — 5 м2 /с).

Определим теперь порядок членов, входящих в пер­

вое уравнение составленной системы:

 

 

 

 

ди

ди

ди

ди

 

 

1

дР

 

~дГ +

и дх + ду_ + W dz_~

 

о

дх

 

10~4

- 4

ю

- 4

іо

- 6

 

“T ö ^

 

То

 

 

 

 

 

 

 

+

2 (о cos <pw2 ( 0 sin <pv +

pi

d2u

(9.39)

 

 

1 0 15

 

~'l0

 

 

 

p

d z 2

 

 

 

 

- 3

 

 

 

io - ß

 

Анализ установленных порядков показывает, что в первом уравнении (9.38) имеются члены, порядок кото­ рых значительно выше остальных, следовательно, они не оказывают заметного влияния на результат решения и ими можно пренебречь.

419


Таким образом, произведя определение порядков всех членов системы п оставляя только соизмеримые по порядкам члены, получим замкнутую систему уравнений, состоящую из шести уравнений, в которых содержится шесть неизвестных. Однако для интегрирования этой си­ стемы необходимо задать вертикальную составляющую ветра w, которая в метеорологии непосредственно не из­ меряется, способов ее вычислить нет. Более удобной формой выражения вертикальной составляющей ветра является замена аналогом, который образуется при за­ мене вертикальной оси г через геопотенциальные высо­ ты главных изобарических поверхностей. Это достигает­ ся переводом составленной системы в изобарическую систему координат. Тогда в этой системе независимыми переменными будут: хр, ур, Р, tp. Так как при таком пе­ реходе оси X, у и время t не меняются, то

х р = х) у р = у;

tp = t.

(9.40)

Ось же 2 связана с новой

осью

Р уравнением

статики

атмосферы

 

 

 

- £ =

- №

 

(9.41)

Тогда функции координат н времени в новой системе координат будут и, v, t, z = H{xp, ур, Р, tp), а функция

dP

дР .

дР

,

дР .

дР

x = l t r =

- W + u - W

+ v lj7 + W -d r является анало­

гом вертикальной

скорости

ветра.

Применяя правило

дифференцирования сложных функций многих перемен­

ных,

окончательно

получим

 

 

 

 

ди

 

и

ди

 

ди^

 

ди

 

dH

2 со sin «рп;

~дГ

 

дх

 

V ■

 

 

HF

-S - W

 

 

 

ду

 

 

 

 

дѵ

 

,

дѵ .

дѵ

 

 

дѵ

 

дН

-f 2

ш sin срн;

и г

-4- ^

— "4~

V ,

 

 

~дР = - g -

ду

 

 

дх

1

ду

 

 

 

 

 

 

 

 

 

 

 

Л - . р .

dH

 

 

(9.42)

 

 

 

 

 

 

 

 

R

дР

 

 

 

 

 

 

 

 

ди

“Ь

 

дх

'■0 ;

 

 

 

 

 

 

 

дх

 

дР

 

 

 

 

 

дТ

 

дТ

+

V

дТ

RT (Та — 3)

 

 

 

 

dt

 

дх

ду

 

Pg

 

 

Так выглядят основные уравнения гидро- и термоди­ намики в изобарической системе координат, пригодные

420


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

Уравнение вихря скорости и его анализ. В целях лик­ видации неравносильности членов, входящих в систему уравнений гидро- и термодинамики, на Них накладывают операцию вихря скорости.

В векторном исчислении вихрем поля называется век­ тор, определенный в каждой точке поля, который яв­ ляется объемной производной этого поля, взятой с об­

ратным знаком. Вихрь обозначается

rot V или П.

 

В прямоугольных координатах вихрь скорости дает

проекции

 

 

 

 

 

(Iw

дѵ

п

ди

дѵ

(9.43)

ду

дг ’

dz

дх

 

Для наших целей используется только вертикальная составляющая вихря скорости, так как другие состав­ ляющие на несколько порядков меньше основных членов уравнений, т. е. Q*= Qy = 0; QZ= Q. Для получения рав­ носильных уравнений продифференцируем первое урав­ нение системы по у, а второе — по х. После этого выч­ тем из второго уравнения первое. Тогда

д дѵ

I-

it

д

дѵ

 

ди дѵ

д дѵ

дѵ

дѵ

дх ~Ж

дх

дх

+

дх Их

+ Ѵ дх ду +

Их

ду

®

+ и

Ф

дх

+

ду

Их

+

 

® @

ду

 

\ (9.44)

ду

dt

Иу'

Ü ду

ду

+ Ну

+...J

д

ди

 

д

ди

 

ди

ди

 

 

д

ди

дѵ

ди

 

 

Рассмотрим разности: ф ; ф;...

 

 

 

 

 

 

^

д

дѵ

 

д

ди ___

д

/

дѵ

ди \ ___

ÖQ

(9.45)

 

 

дх

dt

 

ду

dt

~~

dt

\

дх

 

ду J

 

dt ’

 

 

 

 

 

 

421


Выражение в скобках есть вихрь скорости (9.43), по­ этому разность этих членов равна производной от вихря скорости по времени.

 

дѵ

д

ди

 

д

f

дѵ

ди \

 

dQ

; (9.46)

 

дх

^ ду

дх

 

 

 

дх

ду J

11 ~дх

ди

дѵ

ди

ди

__

ди {

дѵ

ди ^ __q

ди

 

^

 

 

‘1 *

 

 

 

 

~ д у ) ~

 

1 х ’ ' У Л /'

д

дѵ

д

ди

_

д

/

дѵ

ди X

__

дО

(9.48)

~дх' ~ д Р ~ Х ~ду '~дР

~"~~дР\~дх ~'~ду }~~~~дР

 

ИТ . Д .

Приводя подобные члены, получим уравнение в виде

дО .

dQ

дО

до .

дх

дѵ

дх

ди __

~дГ + и

дх

ду

~ дР

дх

дР

ду

дР

 

 

ди

. дѵ \

 

dl

dl

 

-<s+o(£+£)-

 

 

 

(9.49)

 

 

дх

ду

11

дх

® ду

 

где / — параметр

Кориолиса; / = 2ш sirup.

 

 

Произведя оценку порядка членов, входящих в урав­ нение (9.49), и оставляя члены с соизмеримыми поряд­ ками, получим уравнение вихря скорости в сокращенном виде:

до дО dt 4~ it дх

(9.50)

или окончательно

до . д (0 + 1) . „ д ( 0 + 1)

, / ди

дѵ \

ИГ + и ~~Тк----- *-ѵ — By—

------ 1 { 7 7 + ~д7)-

В полученном уравнении изменение вихря скорости по времени связано с изменением параметра силы Ко­ риолиса 2 ш sin ср = / и горизонтального переноса (адвек­

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

422