теорологііческих элемента, чтобы придать составленной системе замкнутый вид. Таким подобием является одно значная зависимость между 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 <pw— 2 ( 0 sin <pv + |
pi |
d2u |
(9.39) |
|
|
1 0 —15 |
|
~'l0 |
|
|
|
p |
d z 2 |
|
|
|
|
- 3 |
|
|
|
io - ß |
|
Анализ установленных порядков показывает, что в первом уравнении (9.38) имеются члены, порядок кото рых значительно выше остальных, следовательно, они не оказывают заметного влияния на результат решения и ими можно пренебречь.
Таким образом, произведя определение порядков всех членов системы п оставляя только соизмеримые по порядкам члены, получим замкнутую систему уравнений, состоящую из шести уравнений, в которых содержится шесть неизвестных. Однако для интегрирования этой си стемы необходимо задать вертикальную составляющую ветра 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 |
|
|
Так выглядят основные уравнения гидро- и термоди намики в изобарической системе координат, пригодные
для интегрирования. Однако при решении этой системы обнаруживается неравносильность отдельных членов, входящих в уравнения, которая проявляется в том, что окончательные значения прогностических величин полу чаются в виде малых долей при оперировании большими числами. Это обстоятельство приводит к тому, что при малейшей ошибке в исходных данных получается боль шая ошибка в предвычисленных. Следовательно, состав ленная система в таком виде не решает поставленной задачи предвычисления значений метеорологических полей.
Уравнение вихря скорости и его анализ. В целях лик видации неравносильности членов, входящих в систему уравнений гидро- и термодинамики, на Них накладывают операцию вихря скорости.
В векторном исчислении вихрем поля называется век тор, определенный в каждой точке поля, который яв ляется объемной производной этого поля, взятой с об
|
ратным знаком. Вихрь обозначается |
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 ’ |
|
|
|
|
|
|
Выражение в скобках есть вихрь скорости (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. Связь же с давлением находится в скрытом виде, а так как нас интересует именно изменение давления, без чего невозможно построение прогностической карты, то выра зим входные параметры через высоты стандартных изо барических поверхностей. Для этого будем считать, что на уровне изобарической поверхности действительный ветер близок к геострофическому ветру. Тогда выраже-