Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МАЛЫХ ВНУТРЕННИХ ДВИЖЕНИЙ НЕВЯЗКОЙ НЕСЖИМАЕМОЙ СТРАТИФИЦИРОВАННОЙ ЖИДКОСТИ

Свиридова Е.Н. 1
1 Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина»
Исследуются двумерные малые внутренние движения невязкой несжимаемой стратифицированной жидкости, которая находится в состоянии равномерного вращения вокруг вертикальной оси. Математическая модель данного физического процесса представляет из себя замкнутую систему уравнений в частных производных, состоящую из уравнений гидродинамики, в которых учтены сила Кориолиса и гравитационные силы. Предполагается, что функция плотности жидкости изменяется в направлении гравитационного поля по экспоненциальному закону. Это приводит к заключению о постоянстве частоты Вейсяля – Брента для рассматриваемой задачи. Колебания, которые учитывают силу Кориолиса и гравитационные силы, называются гравитационно-гироскопическими волнами. Уравнения записаны в системе координат, вращающейся вместе с жидкостью с постоянной угловой скоростью. Система уравнений дополнена необходимыми начальными и граничными условиями. Выписаны явные представления для главных членов асимптотик компонент вектора скорости частиц жидкости. Исследовано влияние различных физических параметров (параметр стратификации жидкости, угловая частота вращения жидкости) на асимптотики этих компонент. Приведены графики зависимости главных членов асимптотик компонент вектора скорости частиц жидкости от времени при различных значениях физических параметров, входящих в исследуемую систему.
стратифицированная вращающаяся жидкость
системы уравнений динамики жидкости
асимптотическое поведение
главные члены асимптотик
1. Габов С.А. Линейные задачи теории нестационарных внутренних волн / С.А. Габов, А.Г. Свешников. – М.: Наука. Гл. ред. физ.-мат. лит., 1990. – 344 с.
2. Ильюшин А.А. Механика сплошной среды / А.А. Ильюшин. – М.: Ленанд, 2014. – 320 c.
3. Pedlosky J. Waves in the ocean and atmosphere: introduction to wave dynamics / J. Pedlosky. – Berlin–Heidelberg: Springer, 2010. – 260 p.
4. Перова Л.В. О колебаниях в стратифицированной и вращающейся жидкости, возбуждаемой плоской, бегущей по дну волной / Л.В. Перова, Ю.Д. Плетнер, А.Г. Свешников // Журнал вычислительной математики и математической физики. – 2000. – № 1. – С. 136–143.
5. Свиридова Е.Н. Асимптотика при t > компонент решения задачи о малых колебаниях вращающейся стратифицированной жидкости в полупространстве. Часть 1 / Е.Н. Свиридова // Вестник Воронежского государственного университета. Серия: Физика. Математика. – 2009. – №1. – С. 150–158.
6. Свиридова Е.Н. Асимптотика при t > компонент решения задачи о малых колебаниях вращающейся стратифицированной жидкости в полупространстве. Часть 2 / Е.Н. Свиридова // Вестник Воронежского государственного университета. Серия: Физика. Математика. – 2009. – № 2. – С. 101–111.

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

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

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

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

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

Рассмотрим движение невязкой несжимаемой стратифицированной жидкости, находящейся в состоянии равномерного вращения вокруг вертикальной оси. Введем декартову систему координат (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) при фиксированном времени.


Библиографическая ссылка

Свиридова Е.Н. МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ МАЛЫХ ВНУТРЕННИХ ДВИЖЕНИЙ НЕВЯЗКОЙ НЕСЖИМАЕМОЙ СТРАТИФИЦИРОВАННОЙ ЖИДКОСТИ // Современные наукоемкие технологии. – 2018. – № 7. – С. 88-92;
URL: https://top-technologies.ru/ru/article/view?id=37084 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674