Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

MATHEMATICAL MODELING OF SMALL INTERNAL MOTIONS IN NONVISCOUS INCOMPRESSIBLE STRATIFIED LIQUID

Sviridova E.N. 1
1 Military Educational and Scientific Center of the Air Force «N.E. Zhukovskiy and Yu.A. Gagarin Air Force Academy»
The two-dimensional small internal motions of the nonviscous incompressible stratified liquid are investigated. The liquid rotates uniformly around the vertical axis. Mathematical model of this physical process is a closed system of partial differential equations, consisting of equations of hydrodynamics, in which the Coriolis force and gravitational forces are taken into account. It is assumed that the liquid density function changes in the direction of the gravitational field according to the exponential law. This leads to the conclusion that the Brunt-Vis frequency is constant for the problem under consideration. Oscillations that take into account the Coriolis force and gravitational forces are called gravity-gyroscopic waves. Equations are presented in a coordinate system rotating together with a fluid. The angular velocity of rotation is constant. The system of equations is completed by necessary initial and boundary conditions. The dominant terms of the asymptotic form of the velocity vector components are obtained. The influence of various physical parameters on the asymptotic form of the velocity vector components are analyzed. Graphs of the dominant terms of the asymptotic form of the velocity vector components against time at varying values of the different physical parameters are given.
stratified rotating liquid
fluid dynamic equations
asymptotic behavior
dominant terms of the asymptotic forms

Работа посвящена исследованию линейной задачи нестационарной теории внутренних волн, распространяющихся в стратифицированной жидкости.

Задачи о распространении внутренних волн и колебаний в стратифицированных жидкостях представляют собой большой интерес, а исследования в этой области имеют высокое практическое значение, так как позволяют прогнозировать и понимать природу обширного круга явлений, связанных с проблемами океана, физики атмосферы и геофизики.

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

Будем рассматривать стратификацию жидкости, вызванную силой тяжести. Распределение частиц, солей и взвешенных суспензий, создаваемое в жидкости этой силой, будет сопровождаться неоднородностью плотности жидкости в направлении гравитационного поля.

Постановка задачи

Рассмотрим движение невязкой несжимаемой стратифицированной жидкости, находящейся в состоянии равномерного вращения вокруг вертикальной оси. Введем декартову систему координат (x1, x2, x3), вращающуюся вместе с жидкостью, причем ось Ox3 направим вдоль оси вращения. Относительно этой системы невозмущенная жидкость является покоящейся. Малые движения этой жидкости в поле сил тяжести в отсутствие других внешних сил принято описывать следующей системой уравнений [1]:

svir02.wmf (1.1)

svir03.wmf (1.2)

svir04.wmf (1.3)

svir05.wmf (1.4)

где svir06.wmf – вектор скорости частиц жидкости; svir07.wmf – вектор Кориолиса (α – удвоенная частота вращения); ρ1(x, t) – изменение плотности, вызванное движениями жидкости; р – динамическое давление; e3 – орт оси Ox3; g – ускорение свободного падения. Вектор ge3 представляет из себя вектор ускорения свободного падения, включающий только гравитационную составляющую, тогда как составляющей, вызванной наличием центробежного ускорения, можно пренебречь [1]. Через svir08.wmf обозначен квадрат частоты Вейсяля – Брента, определяемый формулой

svir09.wmf. (2)

Частота Вейсяля – Брента ω0(x3) является важнейшей характеристикой динамических свойств стратифицированной жидкости. Предполагается, что распределение плотности жидкости ρ0(x3) удовлетворяет естественному физическому условию устойчивости стратификации: svir10.wmf. Стационарная плотность жидкости считается зависящей лишь от переменной x3 [1].

В системе (1) уравнение (1.1) представляет собой закон сохранения импульса; уравнения (1.2) и (1.3) следуют из уравнения неразрывности сплошной среды и предположения о несжимаемости жидкости. Предположение о несжимаемости приводит к уравнению

svir11.wmf,

а уравнение неразрывности, имеющее в линейном приближении вид

svir12.wmf,

будучи переписано в виде

svir13.wmf,

с учетом (1.3) приводит к уравнению

svir14.wmf,

которое с учетом представления (2) для частоты Вейсяля – Брента может быть переписано в виде (1.2).

Колебания, описываемые системой (1), которые учитывают силу Кориолиса и гравитационные силы, называются гравитационно-гироскопическими волнами [1–3].

Если рассмотреть явное выражение для частоты Вейсяля – Брента через стационарное распределение плотности (2), то легко видеть, что эта частота будет постоянной в случае, когда функция ρ0(x3) является экспоненциальной, то есть svir15.wmf, где A, β – положительные константы. Практическая важность рассмотрения этого случая обусловлена тем, что распределение плотности по экспоненциальному закону является больцмановским распределением в однородном поле сил тяжести и тем самым оказывается одним из наиболее часто встречающихся в физике распределений.

Далее будем изучать систему (1) с учетом экспоненциальной стратификации жидкости и, как следствие, постоянства частоты Вейсяля – Брента:

svir16.wmf. (3)

Будем рассматривать двумерные движения вращающейся стратифицированной жидкости, то есть такие движения, которые описываются функциями V и p, не зависящими от одной из пространственных переменных x1 или x2, например, от переменной x2 [1, 4]. Система (1) для таких движений с учетом (3) примет вид

svir17.wmf (4)

Здесь мы уже исключили функцию ρ1.

Из физических и геометрических соображений ясно, что нетривиальные двумерные движения могут происходить лишь в областях типа цилиндров, то есть в областях вида svir18.wmf. Здесь svir19.wmf – вращающаяся система координат.

С помощью замены svir20.wmf приходим к системе уравнений

svir21.wmf (5)

Система (4) дополняется следующими начальными и граничными условиями:

svir22.wmf (6)

Замечание 1. Граничное условие обусловлено существованием функции тока svir23.wmf (см., например, [1]), которую можно связать с компонентами вектора скорости частиц следующим соотношением: svir24.wmf.

Предполагается, что функция ψ(x1, t) удовлетворяет следующим условиям:

Условие 1. При некотором δ > 0 равномерно по svir25.wmf справедлива оценка svir26.wmf, svir27.wmf.

Условие 2. В svir28.wmf существуют следующие производные функции ψ(x1, t): svir29.wmf, svir30.wmf.

Условие 3. Имеет место оценка svir31.wmf.

Условие 4. Функция svir32.wmf финитна: svir33.wmf при svir34.wmf.

В работах [5, 6] доказано существование решения задачи (5)–(6) в пространствах интегрируемых функций С.Л. Соболева и выписано два члена асимптотических разложений компонент решения задачи при svir35.wmf.

Явные представления для главных членов асимптотик компонент вектора скорости

В данной статье проводится исследование зависимости главных членов асимптотик компонент вектора скорости svir36.wmf от времени t при разных значениях физических параметров, входящих в систему (5). Непосредственно из работ [5, 6] следует справедливость теоремы о представлениях главных членов асимптотик компонент вектора скорости задачи (5)–(6):

Теорема 1. Пусть выполнены условия 1, 2, 3, 4. Тогда при svir37.wmf главные члены асимптотик компонент вектора скорости решения задачи (5)–(6) с точностью до постоянного множителя определяются следующими соотношениями:

svir38.wmf;

svir39.wmf

svir40.wmf

svir41.wmf;

svir42.wmf,

где svir43.wmf – частичное преобразование Фурье по переменной x1; svir44.wmf – обратное частичное преобразование Фурье.

Заметим, что компонента svir45.wmf имеет стационарный, не зависящий от времени главный член асимптотики, поэтому для нее приведены два члена асимптотического разложения.

Анализ влияния различных физических параметров на асимптотики компонент вектора скорости

На рис. 1–3 представлены зависимости главных членов асимптотик компонент V1, V3 вектора скорости от времени при различных значениях физических параметров.

svirid1.tif

а) б)

Рис. 1. Зависимости главных членов асимптотик (а) компоненты V1, (б) компоненты V3 от времени при различных значениях частоты вращения жидкости α при фиксированных значениях параметра стратификации β и переменной x3

svirid2.tif

а) б)

Рис. 2. Зависимости главных членов асимптотик (а) компоненты V1, (б) компоненты V3 от времени при различных значениях параметра стратификации β при фиксированных значениях частоты вращения жидкости α и переменной x3

svirid3.tif

а) б)

Рис. 3. Зависимости главных членов асимптотик (а) компоненты V1, (б) компоненты V3 от времени при различных значениях переменной x3 при фиксированных значениях частоты вращения жидкости α и параметра стратификации β

Выводы

Проведенный анализ позволяет сделать следующие выводы:

1) в рамках рассмотренной модели во всех случаях наблюдается естественное затухание колебаний жидкости с течением времени;

2) увеличение параметра частоты вращения жидкости α приводит к росту значений V1(t), V3(t) при фиксированном времени;

3) увеличение параметра стратификации β приводит к уменьшению значения V1(t) и к увеличению значения V3(t) при фиксированном времени;

4) увеличение параметра x3 приводит к уменьшению значений V1(t), V3(t) при фиксированном времени.