Постановка задачи модели
Необходимо отметить, что при исследовании шельфовой циркуляции, наиболее перспективными являются трехмерные математические модели, используемые для воспроизведения основных океанографических полей с учетом специфических процессов, обусловленных особенностями шельфовых зон морей.
Так была сформулирована задача о расчете поля течений в мелководной шельфовой области моря. При этом модель должна позволять исследовать синоптическую и сезонную изменчивость гидрологических полей, учитывать возможность неодносвязности зоны моделирования, особенности конфигурации береговой черты и донной топографии.
Исходя из поставленной задачи, с учетом анализа предшествующих наработок в этой области, выбор наиболее соответствующего типа математической модели представляется очевидным.
Процессы образования, трансформации и динамики водных масс в районах материкового шельфа играют важную роль, как в системе крупномасштабной океанической циркуляции, так и при формировании гидрологического режима отдельных морей или их составных частей. Учитывая нетрадиционно малые пространственные масштабы шельфовых областей, физические процессы формирования гидродинамических полей и их отдельных элементов с достаточной степенью точности могут быть описаны системой полных уравнений движения Рейнольдса для горизонтальных компонент скорости в сферической системе координат {λ,φ,z}, дополненной трехмерными уравнениями турбулентной диффузии тепла и соли. Математическая формулировка модели может быть представлена в виде системы уравнений (4.2.1-4.2.5) в приближении гидростатики:
(4.2.1)
(4.2.2)
(4.2.3)
(4.2.4)
. (4.2.5)
где u, v, w - компоненты скорости; R -радиус Земли; ωz - угловая скорость вращения Земли, РZ - гидростатическое давление на уровне z; Рa - атмосферное давление; KZ и KL - соответственно коэффициенты вертикальной и горизонтальной турбулентной вязкости; Ñ2 - плоский оператор Лапласа; ρ - плотность морской воды; z - превышение уровня над невозмущенной поверхностью; T-температура воды; S- соленость воды. Оператор
В качестве граничных условий на твердом контуре ставится условие непротекания (4.2.6) :
n=0, (4.2.6)
где vn - нормальная к берегу составляющая вектора скорости.
На жидком контуре ставится условие свободного протекания:
, (4.2.7)
где n - нормаль к жидкой границе.
На свободной поверхности задается касательное напряжение трения:
(4.2.8)
, (4.2.9)
где Wλ и Wφ - компоненты скорости ветра; ρa - плотность воздуха.
Величина коэффициента сопротивления са бралась постоянной и равной 0.0012.
В качестве условий на дне (при z=H) принимается квадратичная относительно придонной скорости течения аппроксимация напряжения трения Рейнольдса:
(4.2.10)
(4.2.11)
Коэффициент придонного сопротивления сv считался постоянным и равнялся 0.0025.
Для вертикальной скорости на поверхности используется свободное кинематическое условие (4.2.12),
, (4.2.12)
Условие для вертикальной скорости на дне сформулируем в виде:
(4.2.13)
В дополнение к разработанной модели расчета горизонтальных компонент вектора скорости течения в систему также включено уравнение для расчета вертикальной составляющей движения (w), базирующееся на уравнении сохранения массы в полном, непреобразованном виде.
Еще статьи по теме
Принципы социоприродной этики
Современный исторический момент характеризуется
противоречивостью, мозаичностью и разнообразием социальных форм жизни. Угрозой
сегодняшнему и будущему человечеству являются глобальные процессы деструкции
социального, человеческого, природн ...
Экологические проблемы Приморского края
Вопросы
улучшения экологической ситуации, устойчивого развития, уменьшения потребления
природных ресурсов на единицу валового продукта, роста энергоэффективности
экономики сегодня сформулированы как приоритетные задачи государства.
Не
...