ТОЧНЫЕ РЕШЕНИЯ В ТЕОРИИ НЕЛИНЕЙНЫХ ДЛИННЫХ ВОЛН ЖИДКОСТИ С ВИХРЯМИ И ДАВЛЕНИЕМ

Описание распространения волн в завихренной жидкости является актуальной проблемой математического моделирования вообще и математической теории волн в частности

Процесс распространения длинных волн (по сравнению с толщиной завихренного слоя жидкости) описывается системой уравнений Бенни

Здесь t — время, x — пространственная координата, h(t,x) — неизвестная высота свободной поверхности однородной тяжёлой жидкости, (g ≠ 0), u(t,x,y) — горизонтальная компонента скорости жидкости. Индекс внизу означает производную по соответствующей переменной. В рассматриваемой системе координат ускорение свободного падения g=1 плотность жидкости также равна 1.

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

Введем бесконечную систему моментов

Бенни было показано, что система уравнений (1) порождает бесконечную систему уравнений для моментов

Пусть существует функция распределения (плотность вероятности) f(t,x,v) которую далее пока будем записывать как f(v) рассматривая переменные t,x как параметры. Для f(v) определим характеристическую функцию как

(4)

Здесь:

Тогда корректно следующее выражение функции f(v), через моменты функции f:

Обозначение

Справедлива следующая теорема:

Если функция распределения определена с помощью формулы

Доказательство. Оно состоит в прямой проверке. Именно:

Рассматривая написанные выражения как обобщенные функции, вычислим их действия на основные функции из пространства Шварца vnη(v) где η(v) — некоторая произвольная бесконечно дифференцируемая функция, равная 1 в некоторой окрестности нуля и равная нулю вне большой окрестности нуля. Опуская для кратности в записи η(v), будем писать vn вместо vnη(v). Применим обобщенную функцию

Аналогично

Если моменты цепочки позволяют получить непрерывную быстро убывающую на бесконечности функцию f(v), то цепочка уравнений может быть получена прямым вычислением интегралов от уравнения (7). Уравнение (7) является уравнением с самосогласованным полем, аналогичное тем уравнениям, которые были построены для плазмы А.А. Власовым.

Пусть горизонтальная компонента скорости жидкости u(t,x) не зависит от переменной y.

Тогда

Характеристическая функция в этом случае принимает вид:

Тогда функция распределения будет иметь вид

При этом функции h(t,x), u(t,x) удовлетворяют классической системе уравнений мелкой воды

Рассмотрим теперь случай линейной зависимости функции скорости по переменной y таким образом u(t,x,y)=u0(t,x)+ u1(t,x)y. Заданный функциональный вид решения с необходимостью приводит к тому, что u1(t,x)=c=const и возникает уравнения мелкой воды в завихренной жидкости

В этом случае ряд

Аналогично рассматриваются произвольные полиномиальные аппроксимации горизонтальной скорости по переменной

Примером точного решения системы (12) являются обобщение автомодельного решения уравнений мелкой воды

Постоянная величина c находится из начальных условий.

Уравнение (7) имеет «автомодельную» переменную η = x/t. Фактически подобная ситуация наблюдается, если задача не содержит конечный горизонтальный масштаб. Классическая проблема указанного типа — задача о разрушении плотины. В стандартной теории мелкой воды она решается с помощью метода характеристик

Роль частного решения — «разрушения плотины» играет автомодельное решение f(x/t,v). Для автомодельного решения получается уравнение от двух независимых переменных

Следуя методам, развитым для кинетической теории плазмы

Следовательно,

Откуда следует, что можно заменить функцию h на μ, и получить уравнения

Решая методом характеристик уравнения (15), можно найти общий интеграл характеристической системы

где R(μ) является произвольной гладкой функцией. Таким образом, при R=0 решение (16) соответствует классическому автомодельному решению уравнения.

Таким образом, доказана следующая теорема:

Если

Смысл получаемых решений выясняется из уравнений (13), это решения типа простых волн. Следовательно, решения уравнений (15) получается из простых волн с помощью функционального преобразования (16), так что практически во всех значимых случаях, решение будет иметь особенности, типа бора (гидравлический скачок = ударная волна при газодинамической аналогии).

Указанные решения обобщают классические решения в теории мелкой воды и в газовой динамике. Этим самым определяется обширный класс решений в жидкости с распределенными вихрями, которые можно назвать неоклассическими

Другой обширный класс решений появляется, когда наличие вертикального перемешивания в жидкости становится определяющим. Например, при разрушении плотины «недалеко» от плотины. Но и в этом случае в качестве нулевого приближения берется неоклассическое решение, а малая вихревая дисперсия может быть учтена методами возмущений, что предполагает, конечно, численное моделирование

Плоскопараллельное движение жидкости в модели Бенни без источников и стоков точно соответствует бесстолкновительной динамике нейтрального газа, т.е. рассматривается следующее уравнение:

Простота уравнения (17) и возможность его общего решения полезны для качественного анализа задачи, причем с обобщением выводов на нелинейную область. Важно сравнивать на классической задаче динамику плоскопараллельных движений в модели Бенни и классической нелинейной теории мелкой воды.

Рассмотрим задачу о распаде плотины в рамках чисто инерционного подхода, описывающегося уравнением (17). Общее решение уравнения (17) имеет вид

Здесь F — произвольная функция. Необходимо отметить, что в отличие от обычно рассматриваемого случая здесь одна степень свободы в пространстве скоростей v. При t=0 имеем

Здесь m,k,T — константы задачи. Искомая функция распределения примет следующую форму

Отсюда находятся плотность числа частиц газа (высота свободной поверхности жидкости) и средняя макроскопическая скорость газа

Где

Изложенной модельной задаче соответствует классическая задача о разрушении плотины в рамках гиперболической теории мелкой воды. Начальные условия

Решение находится с помощью теории простых волн, а именно

при

Свободная поверхность между фронтами имеет параболическую форму.

При сравнении полученных решений достаточно сопоставить нулевой момент решения уравнения (17), поскольку именно высота свободной поверхности является наблюдаемой величиной. Во-первых, автомодельные переменные одинаковы. В точке x=0 при всех t из теории мелкой воды имеем

Отсюда можно сделать вывод об удовлетворительном при ограниченном времени наблюдения качественном анализе опрокидывания свободной поверхности в рамках приближения длинных волн на основе уравнения (17). С учетом неоднородной глубины можно анализировать на его основе такие явления, как цунами.

Система уравнений Бенни с давлением имеет следующий вид:

Здесь P(t,x) — заданное внешнее давление. С помощью уравнения (17) можно строить точные решения уравнения (24), используя внешнее давление, заранее не известное, как условие согласования нулевых моментов. А именно, пусть выполнено условие

Тогда уравнение (17) переходит в (24). Общее решение уравнения (17) имеет вид (18), то есть

Здесь F — произвольная начальная гладкая функция. Основным результатом, на котором основывается эта конструкция, является одинаковость по форме нулевого моментного уравнения для уравнений (17), (24)

Тогда мы можем выбрать момент A0(t,x) общим для уравнений (17), (24), так что Px(t,x) определяется однозначно для всех моментных уравнений (3).

Рассмотрим, например задачу о распространении газа, сосредоточенного в полупространстве x<0 при t=0 в вакууме. Гидродинамическая задача, которой ей соответствует – задача о распаде плотины, так что вода находится при x<0. Тогда мы имеем решение (19), (20), (21). Решение (19) – (21) автомодельно, так как в задаче нет характерного размера. Тогда это соответствует нелинейной задаче (24) с внешним давлением

Аналогичные результаты можно получить, положив

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

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