одного итерационного процесса, а обычный путь - трех. После решения
уравнения Пуассона (23) методом установления [ 5 ] — переходим к ре-
шению уравнения (22) на новом временном 
интервале. В случае сжимае-
мой жидкости внутри только что описанного цикла 
между первым и вто-
рым его этапом 
решается система уравнений (2) с целью нехождения
гидростатического давления, плотности воды, её минерализеции и теп-
ловых характеристик. Граничные условия для поля "вихря" рассчитнва-
ются по известным методикам 
| 4 ] 
-
Теперь следует сказать несколько слов об определении уровня сво-
бодной поверхности, значение котброго,Как правило, входит в формулу
для вычисления гидростатического давления 
[ 4 ] 
.
Поле уровней свободной поверхности обычно вычисляют, решая следу-
ющее уравнение 
[ 4 ] ;
ЁЁ” 
+ аш (а - #) Йа т 
` 
(24)
тде #, , %9 
— кооішинатн уровня` свободной поверхности и дна соответ-
ственно, о» — вектор горизонтального перемещения слоя воды 
#, -#9 .
Но уравнение (24) можно привести к виду более удобному для расче-
та и требующему меньшей оперативной памяти ЭВМ.
Для этого уравнение несжимаемости (8) проинтегрируем от уровня
свободной поверхё'ости ВодЫ В івпцоеме ло его дна
_]%щ_-[%д+}_д‘;”_а‚ъ 
(25)
Внося интегЁзЁл под знак пЁ›пи:‹ввшюй и п%имв!-шя теорему о среднем
значении интеграла, получим 
)
агаа й,. 
(26)
— \/ 
= аш (#„т 9) И5 -Цр Е
#= #
Складывая уравнения (25) и (24), получим новую Форму записи урав-
нения (24): 
:
а#л -
ЭЁ 
а-
Следует только помнить, что \\ - это та вертикальная составляю-
щая скорости, которая получается при расчете из уравнения неслимае-
мости. Уравнение (27) позволяет, в целях повышения устойчивости, поль-
зоваться методом сдвинутых сеток [ 1 ].
С целью количественной проверки по предложенным методем счета бы-
ли рассчитаны поля течений в водоеме, анелогом которого является де-
тально изученная физическая модель [ 2] 
„ Гидродинемическая задача
решелась неявно, методом прогонок, РеЗУЛЬТВТН расчетов совпадаели ка-
чественно с опытными данными уже на первых шегах по времени. Количес-
твенно же соответствие рассчитанных и измеренных результатов наступа-
ло через один час на ЭВМ ВС-ТО22 ( 
рис. 1). Ошибка не превышела
(:2'7)
50