номер шага по времени. Разностное уравнение 
для \/]і‘- аналогично
уравнению (2.36). Схема устойчива при А% < Ах/ \5 ЕН, )
Если необходимо рассчитать станионарное поле скорости, задачу
можно несколько упростить, вводя функцию полного потока { такую,
что
У, = —д/дг, 
У, = дудх, 
(2.40)
н 
н
тде У, = 
{ У @2 
М, = 3; \„б2. Изолинией функции $ является ли-
ния тока для полного потока У = (М,, У,), а разность значений { в
двух точках определяет расход воды между этими точками . Уравнение
для функции полных потоков имеет вид /93/
_( МУ ду , _до ( Ку 99 оы
дх б дрдх) б оааа 
н° 
240
здесь Ку— коэффициент вертикального турбулентного обмена.
В модели Фельзенбаума К_у является функцией скорости ветра и
глубины
Ку= ЭУН/р 
(242)
где у = 0,065 кг/м3 — постоянная.
В монографии /93/ рассмотрены возможности численной реализации
модели на основе метода конечных разностей. В работе /77/ приведен
усовершенствованный вариант модели, позволяющий учесть трение на
дне (у Фельзенбаума на дне задается условие прилипания). Численное
решение дифференциального уравнения для функции полных потоков
осуществляется с использованием метода конечных элементов. Послед-
ний позволяет учесть реальную конфигурацию водоема точнее, чем ме-
тод конечных разностей. Кроме того, появляется возможность увели-
чить разрешающую способность в наиболее важных участках водоема
без увеличения общего числа узлов.
Как правило, для резлизации моделей требуются быстродействую-
щие ЭВМ. Модели реализуются в виде программ, время расчета зави-
сит от быстродействия ЭВМ, сложности модели, размеров области ре-
шения. Например, модель, основанная на разностной схеме (2.36)—
—(2.37) 
и реализованная в виде ФОРТРАН — программ на ЭВМ
ЕС-1045, требует для расчета на сутки с временным шагом 45 с для
пространственной сетки размером 40 х 30 около 16 мин машинного
времени. Выбор модели для конкретных водоемов должен осущест-
вляться с учетом особенностей данного водоема (глубины, наличия
стратификации и т. д.) и поставленной задачи.
59