CC BY
20
Вестник РУДН Серия Математика. Информатика. Физика. № 2 (1). 2010. С. 5-10
Математика
УДК 517.955.8
Асимптотики решений двух задач динамики экспоненциально стратифицированной жидкости
А. В. Глушко, Е. Н. Свиридова
Кафедра уравнений в частных производных и теории вероятностей Воронежский государственный университет Университетская пл., д.1, г. Воронеж, 394006, Россия
Работа посвящена изучению компонент решений двух задач динамики, описывающих малые колебания экспоненциально стратифицированной и равномерно вращающейся жидкости в декартовой системе координат (хг,х2,х3), жёстко связанной с вращающейся жидкостью. Жидкость стратифицирована вдоль оси Охз, совпадающей с осью вращения: ро(х3) = Аехр(—2^ж3), где ¡3 > 0 — параметр стратификации. Построены асимптотики при £ ^ компонент решения обеих задач.
Ключевые слова: дифференциальные уравнения, системы, переменные коэффициенты, существование решения, асимптотическое поведение.
1. Введение
В настоящей работе рассматриваются две задачи, описывающие двумерное движение стратифицированной жидкости, т.е. такое, при котором компоненты решения не зависят от одной из пространственных переменных, х\ или жо (для определённости — от ж0).
В первой задаче движение жидкости изучается в рамках модели, принадлежащей А.Г. Свешникову, Ю.Д. Плетнер, Л.В. Перовой. В работе [1] этих авторов была доказана стабилизация решения при большом времени. Нами, при несколько других требованиях на функцию из граничных условий, выписаны точные асимптотики компонент решения. Рассматриваемая система уравнений имеет вид
V о.
~т - аУо +
1 др
ро(хз) дх\
0,
Ж + =0.
д2Уз дг2
+ + ро(хз) дЪ ' 0.
( др ^
\дхз )
(1)
т + т = 0
дх1 дх3 '
где х1 € К, х3 > 0, £ > 0; а = (0, 0, а) — вектор Кориолиса; = 2/Зд — квадрат частоты Вейсяля-Брента, р — динамическое давление, ро(хз) — плотность жидкости в невозмущённом состоянии. С помощью замены
а = рр-1(хз) = А~1ре0^Х3
Статья поступила в редакцию 13 октября 2009 г.
приходим к системе уравнении
( дУх да п
—^ - аУ2 + — = 0, дъ дх\
§ + =0,
Ж + ^ + д (II - 2*) =0,
(3)
т + т =0
дхх дхз
Дополним систему (3) следующими начальными и граничными условиями:
дУз
Ук(х, 0) = 0, к = 1, 2, 3, (х, 0) = 0, а(х, 0) = 0,
дф (4) уз (хl,^, ¿)и=о = - дхх~(х1, г).
Замечание. В силу условия соленоидальности, задаваемого четвёртым уравнением системы, с компонентами Ух, Уз вектора скорости У = {У1, У2,Уз} можно связать функцию тока Ф = Ф(хх,хз, ¿) следующим соотношением: {Ух, Уз} = д Ф дФ 1
——, - —— >. Этим соотношением обусловлено граничное условие. дхз дхх )
В рассматриваемой модели А.Г. Свешников, Ю.Д. Плетнер, Л.В. Перова предполагали, что учёт переменного коэффициента плотности ро (хз) в уравнении неразрывности не является необходимым, обосновывая это экспериментальными соображениями. Мы решили посмотреть, как с математической точки зрения этот эффект учитывается не только в уравнениях Эйлера, но и в уравнении неразрывности. С этой целью изучена вторая задача.
Во второй задаче четвёртое уравнение системы (1) рассматривается в общем виде
ё1у (ро (х3) У) = 0 (5)
или
1Г- + Р - 2"У' = 0-
дхх дхз
Полученную таким образом систему преобразуем с помощью замены (2), как и в первой модели. Начальные условия из (3) для новой системы сохраняются, а граничное условие имеет вид
У <х-, *)!—о = - ^^ ■ (6)
Для этой модели также выписаны точные асимптотики компонент решения.
Обобщённое решение обеих задач строится с помощью преобразований Фурье и Лапласа. С помощью них же доказывается существование решения. Получение асимптотических представлений компонент решения задач основано на изучении свойств фазовых функций в интегралах, образующихся при вычислении обратных преобразований Лапласа и Фурье, и применении метода стационарной фазы (см. [2]).
Приведём краткий обзор доказанных результатов для обеих моделей.
Предполагается, что функция ф (хх,£) удовлетворяет следующим условиям
Условие 1. Функция ф (хх,£) равномерно по (хх,£) € ограничена:
1ф (хх^)| < с.
Условие 2. В ¿1 (К х (0; то)) существуют следующие производные функции ф(жь ■£):
( д0\п Хп (Х1, ¿) = - ф (Х1, V), п = 2, 3,4.
Условие 3. Имеет место оценка
(1 + Ы)
+С
/ (1 +Ы) (
д0
1 - эщ) ф (Х1,
д,х1 < то, п = 2,4.
Условие 4. Функция ф (х1, Ь) финитна: ф (х1, Ь) = 0 при Ь > N.
Замечание. Из условий 1 и 4 следует, что существует ё > 0 такое, что справедлива оценка |ф (ж1, t)| < сехр(—5¿), (ж1, ¿) € М+.
Введём нормы
1СЮ С
I I (1+^2) ШХ1, ¿)]| (1 +
1/0
. 0 -с
+ С
Н^ & *)М|£2(®2х(0;с)) = I ! \в ,
0 К2
+С
{{{я (х, ¿)>»^2(К2Х(0;С)) = //\ ?~013хзд (х, г) \ 0 ЛхЛг.
(7)
Через 1|-Уь2(Е3) , 1|-Н^2сж+)> ^•^Ь2(^2Р(К)Х(0,С)) будем обозначать стандартные
нормы в соответствующих пространствах ¿0(К3), Ь0(М+), Ь0(Ж^К) х (0, то)), где (М) — пространство Соболева-Слободецкого с индексом р на К.
2. Задача 1. «Модель Свешникова»
В работе [3] доказано существование решения задачи, получены оценки для норм компонент решения и их производных, входящих в систему, в пространствах ¿2(К+), ¿2("^''(К) х (0, то)), выполнена проверка начальных и граничных условий. Приведём результаты для компонент вектора скорости.
< С ^ , ¿)|Ь2(Кх(0;с)) +
дф
дх1
Ь2(Кх(0;с))
+ {{Ф>> —1/0 + {{Ф>>
1/0
|МVo(Xl,Xз, ¿)М|Ь2(К+Х(0;с))
<
< С'
д ,, тф(х1,
+
Ь2(Кх(0;с))
д0ф
дЪдхл^
Ь2(Кх(0;с))
+ {{Ф>> —1/0 + {{Ф>>
r3(Xl,Xз, ¿)М|Ь2(К+х(0;с)) < ||ф(жЬ ¿)М|Ь
2(Кх(0;с))
+ с{{ф>>
1/0
1/0
В работе [4] исследовано асимптотическое поведение компонент решения при £ ^ +то. Доказана следующая теорема
Теорема 1. Пусть выполнены условия 1, 2, 3, 4- Тогда при £ ^ для компонент решения задачи (3), (4) справедливы следующие асимптотические представления
У (х, ¿) = а (х3) £
-1.5
д2 \4
1 - ЯГТ2 ф (у 1, т)йтёг/1+
о
д 12
+ (1 + |х1|)(х3-1 + х-2)О( Г2)
( д2 \4 (1 + - о^) ф (т)^т
У2 (х, ¿) = а/ ехр (/хз)
31
д
-аехр(//хз) 1
ф (з 1, 0) ехр ^-хз з21 + /2
ф (в 1, 0) ехр ( -хз\/+ /
+
+ Ь(хз)Ь
1.5
о
д2 4
1 - ЯГГ2 ) ф (т)^Т<1у 1 +
д 12
+ (1 + |х11) (х-1 + х-2) О (Г2)
(1 + - ф (
СЮ сю 2 4
Уз (х, ¿) = к (хз) Г55 I I (х1 - г/1) - ф (г/1, т)ётёу 1 +
-ю 0
+ (1 + |х1|)х-зО( Гз)
д
— (х, = Ро (хз) с (хз) г-1.5
-ю о
+ (1 + |х11) (х-1 + х-2)О(Г2)
(1 + - ф (
д у\.
1 - ф (т)ётё2/1+ (1 + И)(1 - ф (
д д
— (х, ^ = Ро (хз) ехР (//хз)
ф (з 1, 0) ехр ( -хзлЫ +/2
+
сю сю ^ 4
+ ё (хз) ро (хз) Г 2 I I (х1 - 2/1)^1 - ф (У1, г)ёгё?/1+
-ю о
+ (1 + |х1|)х-зО( Гз)
(1 + - ф (т)ёгёт
2
Здесь О (t-2) и О (t-3) есть функции переменных t > 0, х3 > 0, x1 G R, равномерно not > to > 0, х3 > 0, xi G R, удовлетворяющие соответственно оценкам: |О(Г2) | ^ct-2, |О(Г3) | ^ ct-3; коэффициенты k (х3),а (х3), b(x3), с(х3), d(x3) равномерно по х3 : 0 < е < х3 < N < то ограничены и выписаны в явном виде в работе [4].
3. Задача 2. «Наша модель»
Аналогично первой задаче проведено полное исследование всех компонент решения системы и оценка их норм в шкале пространств Соболева-Слободецкого. Приведём результаты для компонент вектора скорости V = (У, V2, V3). Справедливы следующие оценки
, х3, ^))))l2(r+x (0;ю))
<
||ф(хь t)\\L
2(Rx(0;c))
+
дф
дх1
l2(rx(0;c))
+ «ф))-1/2 + т)
' 1/2
, х3, ^))))L2(R+x (0;^))
<
< С'
д
тф(х1,t]
+
L2(Rx(0;^))
д2ф
д д х1
L2(Rx(0;^))
+ «ф))-1/2 + т)
(((V3(x1,x3, t))))L2(Rl х(0;ю)) < ||ф(хъ
2(Kx(0;m))
+ с{{ф))
1/2
1/2 .
Теорема 2. Пусть выполнены условия 1, 3, 4, условие 2 при п = 2 . Тогда при t ^ +то для компонент вектора скорости решения второй задачи справедливы следующие асимптотические представления
сю сю 2
Ух (х, ¿) = Ш1 (х) Г3/2 I ! - ^2) ф (уи т)ётёу 1+
-с 0
/ сю сю 2 N
+ (1 + |х 11) х-20 (Г2) П | (1 + И)(\ - ф (г, т)^ёт
I а . / \
m1 (х)
л/2 exp (2/3x3) n3/2x3po (0) ' V а2 - ш2
тг , . аexp (2/х3) х2 — х3
2 (х, t)= /^Г 3) •—1-232 (X1, 0) +
Po (0)^ (х1 + х|)
СЮ СЮ 2
+ m2 (х) t-3/2 f J - дд/f) Ф (У1, T)drdy 1+
-С o
+ (1 + |х 11) х-2О (t-2)if j (1 + N)(l - д^) Ф (z, r)dzdr
o -С
\/2exp (2/3x3) I а / \
^3/2x!PO (0) V O2-^! ' °°Vto - ~4J ;
m2 ( х)
д2
2
Уз (х, t) = m3 (х) Г1 J J - д^) ф (у 1, r)drdyi +
-œ 0
2
+ (1 + |х 11 ) х-20 (t-3/2)lj J (1 + N)(l - ф (z, r)dzdr
0
2 ехр (2(Зх3) [ Л ео8(Лж1)
ТО3 (х) = - п2р0 (0) ео8 ('Шо)] Ц+Л2? ёЛ-
о
Здесь О (¿-2) и О (¿-3/2) есть функции переменных Ь > 0, хз > 0, х1 € К, равномерно по £ > ¿0 > 0, х3 > 0, х1 € К, удовлетворяющие соответственно оценкам: |О (¿-2)| < сГ2, |О (¿-3/2)| ^сГ3/2.
4. Заключение
Итак, нами построены точные асимптотики компонент решения двух задач динамики, описывающих малые колебания экспоненциально стратифицированной и равномерно вращающейся жидкости.
Литература
1. Перова Л. В., Плетнер Ю. Д., Свешников А. Г. О колебаниях в стратифицированной и вращающейся жидкости, возбуждаемой плоской, бегущей по дну волной // Журнал вычислительной математики и математической физики. — 2000. — № 1. — С. 136-143.
2. Федорюк М. В. Метод перевала. — М.: Наука, 1977.
3. Глушко А. В., Свиридова Е. Н. Оценка поведения при £ ^ решения задачи о малых колебаниях вращающейся стратифицированной жидкости в полупространстве // Труды математического факультета ВГУ. — 2007. — № 11. — С. 35-48.
4. Свиридова Е. Н. Асимптотика при £ ^ компонент решения задачи о малых колебаниях вращающейся стратифицированной жидкости в полупространстве. Часть 1. // Вестник ВГУ, серия «Физика, математика». — 2009. — № 1. — С. 150158.
UDC 517.955.8
The Asymptotic Behavior of Solutions of Two Problems of Dynamic of Exponentially Stratified Fluid
A. V. Glushko, E. N. Sviridova
Department of Partial Differential Equations and Probability Theory Voronezh State University University square, 1, Voronezh, Russia 394006
This article is devoted to studying the components of solutions of two problems of dynamics. The problems describe the small fluctuations of the exponentially stratified and uniformly rotating fluid in the Cartesian system of coordinates (xi,x2,%3) rigidly connected with the rotating fluid. The fluid is stratified along the axis Ox3 coinciding with the axis of rotation: Po (x3) = A exp(-2^x3), where /3 > 0 is a parameter of stratification. The asymptotics of the components of solution are constructed for both problems.
Key words and phrases: differential equations, systems of differential equations, variable coefficients, existence of solution, asymptotic behavior.
