CC BY
40
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
УДК 532.5.013.4
А. А. Гаврильева, Ю. Г. Губарев
УСТОЙЧИВОСТЬ УСТАНОВИВШИХСЯ ПЛОСКО-ПАРАЛЛЕЛЬНЫХ СДВИГОВЫХ ТЕЧЕНИЙ ИДЕАЛЬНОЙ СТРАТИФИЦИРОВАННОЙ ЖИДКОСТИ В ПОЛЕ СИЛЫ ТЯЖЕСТИ
Исследована задача линейной устойчивости стационарных плоско-параллельных сдвиговых течений невязкой непрерывно стратифицированной жидкости, находящейся в поле силы тяжести между двумя неподвижными непроницаемыми твердыми параллельными бесконечными пластинами.
Прямым методом Ляпунова показано, что данные течения абсолютно неустойчивы относительно малых плоских возмущений (на полубесконечных интервалах времени). А именно построена априорная оценка снизу, которая свидетельствует об экспоненциальном во времени нарастании рассматриваемых возмущений. А также прямым методом Ляпунова получены достаточные условия практической (на конечных временных промежутках) неустойчивости данных течений по отношению к малым плоским возмущениям.
Ключевые слова: идеальная жидкость, стратифицированная жидкость, установившиеся течения, прямой метод Ляпунова, априорная оценка, условия неустойчивости.
A. A. Gavrilyeva, Yu. G. Gubarev
Stability of steady plane-parallel shear flow of ideal stratified fluid in a gravity field
The problem of linear stability of stationary plane-parallel shear flows of an inviscid continuously stratified fluid in a gravitational field between two fixed impermeable solid parallel infinite plates is studied.
The Lyapunov direct method shows that the data flow is absolutely unstable with respect to small plane perturbations (in the semi-infinite time intervals). Namely, a priori estimate below has been constructed, which shows an exponential increase in time under consideration of the perturbations. As well as sufficient conditions for the practical (on finite time intervals) instability of these flows with respect to small perturbations of the flat has been found by the direct Lyapunov’s method.
Key words: ideal fluid, stratified fluid, steady flow, the direct Lyapunov method, a priori estimate, the conditions of instability.
ГАВРИЛЬЕВА Анна Андреевна - ведущий инженер отдела материаловедения Института физико-технических проблем Севера им. В. П. Ларионова СО РАН.
E-mail: gav-ann@yandex.ru
ГУБАРЕВ Юрий Геннадьевич - к. ф.-м. н., доцент, с. н. с лаборатории вихревых движений жидкости и газа Института гидродинамики им. М. А. Лаврентьева СО РАН.
E-mail: Yu.G.Gubarev@mail.ru
Проблема устойчивости стационарных плоскопараллельных сдвиговых течений идеальной непрерывно стратифицированной жидкости в поле силы тяжести по отношению к малым плоским возмущениям является одной из фундаментальных задач гидродинамики, физики атмосферы, океанологии и других родственных областей науки [1], которая не утратила своей актуальности и в настоящее время.
Основные результаты изучения устойчивости установившихся плоско-параллельных сдвиговых течений невязкой непрерывно стратифицированной жидкости, находящейся в поле силы тяжести относительно малых плоских возмущений,
получены спектральным методом с использованием интегральных соотношений [2-4]. К сожалению, эти результаты справедливы лишь для возмущений в форме нормальных волн [1, 4]. Кроме того, следствием применения интегральных соотношений служит то, что реально исследуются не все возможные возмущения в виде нормальных волн, а только удовлетворяющие дополнительному требованию к дифференциальному оператору Тейлора-Гольдштейна [2-4].
Наконец, как вытекает из выполненных в работе [5] энергетических рассуждений, устойчивы лишь состояния покоя, у которых плотность рассматриваемой жидкости возрастает в направлении нижней пластины, причем по отношению к такому классу малых возмущений, что для него малые возмущения плотности на границе равны нулю. Наложение же сколь угодно слабого сдвига скорости дестабилизирует данные состояния покоя. Более того, за исключением описанных здесь состояний покоя, отыскать условия устойчивости каких-либо иных стационарных течений изучаемой жидкости относительно малых возмущений типа [5] с помощью энергетических соображений нельзя.
Все это дает возможность предположить, что имеется абсолютная неустойчивость установившихся плоско-параллельных сдвиговых течений идеальной непрерывно стратифицированной жидкости в поле силы тяжести по отношению к малым плоским возмущениям. Данное предположение может быть проверено прямым методом А. М. Ляпунова [6-8] -самым мощным аналитическим методом современной математической теории гидродинамической
устойчивости. Однако главная трудность при использовании этого метода заключается в том, что до сих пор неизвестны регулярные способы построения функционалов А. М. Ляпунова. Поэтому нами развивается новая аналитическая методика [6], позволяющая или демонстрировать теоретическую абсолютную неустойчивость, или получать теоретические и практические достаточные условия
неустойчивости тех либо других исследуемых
течений жидкости, газа или плазмы относительно малых возмущений. Суть новой методики состоит в алгоритмическом конструировании функционалов А. М. Ляпунова, которым, с одной стороны, характерен монотонный рост со временем в силу соответствующих начально-краевых задач для малых возмущений, а с другой - такая форма зависимости от малых возмущений, что из нарастания данных
функционалов во времени следует рост малых возмущений со временем в согласии с определением неустойчивости. При этом возникающие условия
нарастания во времени функционалов А. М. Ляпунова как раз и будут достаточными условиями неустойчивости рассматриваемых течений по отношению к
малым возмущениям. Практическая значимость работы заключается в том, что для получения сведений об устойчивости исследуемых течений не нужно применения сложного и трудоемкого теоретического анализа свойств как самой исследуемой математической модели, так и знать явный вид решений краевых и/либо смешанных задач, описывающих изучаемые течения и возмущения.
Таким образом, в работе прямым методом А. М. Ляпунова [6-8] будет показана абсолютная неустойчивость стационарных плоско-параллельных сдвиговых течений невязкой непрерывно стратифицированной жидкости, которая находится в поле силы тяжести между двумя неподвижными непроницаемыми твердыми параллельными бесконечными пластинами, относительно малых плоских возмущений. Кроме того, будет построена априорная оценка снизу, говорящая об экспоненциальном росте со временем исследуемых возмущений, и получены достаточные условия практической неустойчивости данных течений по отношению к малым плоским возмущениям.
Нами рассмотрены плоские течения идеальной несжимаемой неоднородной по плотности жидкости между двумя покоящимися непроницаемыми твердыми параллельными неограниченными стенками в поле силы тяжести.
Данные течения характеризуются эволюционными решениями начально-краевой задачи вида [1-3]
pDu = - рх, pDv = - р y - pg, Dp = 0, ux + vy = 0 bt ; (1) v = 0 на дт; и(x , у, 0) = u0 (x, у), v(x, у, 0) = v0 (x , у),
где p(x, у, t) - плотность жидкости; u(x, у, t\ v(x, у,t)
- составлжощие по ля скорости жидкости; р(х, у, t ) -поле давления; потесциальное поле силы тяжести
УФ = (0; g), g = const ; D = д / dt + ид /дх + vd / ду
- диффе(енциальный оператор; X, у - декартовы
координаты;Т = {(х У): — го < х < 0 < у < Н}
- область ' ' течения жидкости;
дт = {(х, у): - «з < X < +гс>, у = 0, И} - граница области течения; U q , Vq - начальные компоненты поля скорости жидкости; t - время; Н = const - ширина зазора между стенками. Нижними индексами из независимых переменных обозначаются частные производные искомых функций. Предполагается, что начальные составляющие и 0 и Vq поля скорости жидкости обращают в тождество четвертое соотношение смешанной задачи (1). Более того, считается, что начальная компонента VQ удовлетворяет пятому равенству этой задачи.
Начально-краевая задача (1) обладает интегралом энергии в форме
_ 7Н ( и2 + V2 )
К = III р------------+ рФ ауал = const
0 ^ ' б ’
когда ее решения либо периодичны, либо локализованы вдоль оси абсцисс.
Также у смешанной задачи (1) имеется еще один сохраняющийся на эволюционных решениях функционал. Действительно, если ввести функцию
Л(Х, у, t) , сохраняющую свои значения в любой
жидкой частице при ее движениях ( DA = 0 ), то интеграл +^и
I = Ц F (X)dxdy = const, (3)
—оэ 0
где F(X) - произво льная функция своего аргумента, будет сохраняться на решениях начально-краевой задачи (1).
Смешанная задача (1) обладает точными стационарными решениями вида
Р = Ро (У), и = U{у), V = 0 ,
Р
Ро
(4)
(где р0, и - некоторые функции ординаты, р0 - постоянная величина, у1 - переменная интегрирования). Данные решения соответствуют стационарным плоско-параллельным сдвиговым течениям невязкой непрерывно стратифицированной жидкости между двумя неподвижными непроницаемыми твердыми параллельными бесконечными пластинами в поле силы тяжести.
Цель дальнейшего рассмотрения состоит в том, чтобы выяснить, могут ли стационарные решения (4) быть неустойчивыми по отношению к малым плоским возмущениям.
Для достижения данной цели производим линеаризацию смешанной задачи (1) около точных стационарных решений (4), приводящую к начальнокраевой задаче вида
( < тт -Ро и, +ии х +У ~Г \ йу
Р
Ро
dv' TJ dv' — + U — dt дх
-Р' y -Р
. (5)
dt
■ 0 U +v' = 0 , в Т ;
iy
V = 0 на дт;
и' (x, у, 0) = и'0 (x, у), v'{x, у, 0) = v'0 {x, у), где и'(x, y, t), v'(x, у , t\p(x, y, t\ p(x, у , ) )
- малые плоские (озмущения полей скоростей, возмущения плотноути, возмущения давления;
U 0, V 0 - начальные составлязощие возмзщенного
поля скорости жидкости. Полагается, что функция
V10 , с одной стороны, превращает в тождество пятое равенство смешанной задачи (5), а с другой -вместе с функцией и 0 удовлетворяет четвертому соотношению этой же задачи.
К сожалению, линейный аналог интеграла энергии для начально- краевой задачи (5) для любых малых плоских возмущений с помощью энергетических соображений найти нельзя.
Далее прямым методом А. М. Ляпунова будет показано, что стационарные течения (4) абсолютно неустойчивы относительно малых плоских возмущений. Для того чтобы продемонстрировать неустойчивость какого-либо стационарного решения (4) смешанной задачи (1) по отношению к малым плоским возмущениям, нужно суметь выделить среди этих возмущений всего только одно, но за то экспоненциально быстро нарастающее во времени возмущение.
Исходя из данного обстоятельства, поиск требуемого возмущения осуществляется ниже в подклассе плоских течений, для него малые плоские возмущения (5) являются отклонениями траекторий движения жидких частиц от соответствующих линий тока установившихся течений (4). Эти возмущения могут быть описаны с помощью поля лагранжевых смещений £(х, у, ¿) = ) [9], которые определяются
уравнениями
dU
dy
& = v'-U^х . (6)
Посредством соотношений (6) перепишем начально-краевую задачу (5) в виде
р0 (#ltt + 2и+ и ) = -р\’ у = 0 > (7)
Р0 {^2„ + х, + и Ч2 XX )=- Ру +Ро у8%2’
%2 = о на дт;
£ {х, у, 0) = £ю (х у), £ {х, у, 0) =) (£ )0 {х, у), £ (Х у, о) = £о (х у), £, (х у, 0) = (&: )0 (х, у).
(где ^ 0, 0 - начальные составляющее поля лагран-
жевых смещений £ , а )0 и \<^ь )0 - начальные
компоненты его частной производной первого порядка по времени). Предполагается, что функции £ , д20 обращают в тождеетво второе уравнение из системы (7).
Этору подклассу линейным аналогом интеграла энергии для начально-краевой задачи (7) служит функционал вида
Ро fat + UL )2 +(%2, + U%2x )2 )
- 2PoU ( + £х^2, )
- 2Р U 2 (1x )2 + (£x )2 )'- Р
О у
8 (£ )2
йуйх. (8)
Непосредственной проверкой можно убедиться, что первая вариация 8] интеграла ] = Е1 + I (2), (3) зануляется на стационарных решениям (4) лишь тогда, ко гда выбирается связь функции X = X (р), и вариация возмущения плотности на границе равна нулю Зр\^г = 0. В то же время вторая вар нация §2 3 функционала 3 совпадает по форме с интегралом Е (8) [5].
Знакоопределенность подынтегрального выражения (8) имеется только на состояниях покоя и при следующем распределении плотности рассматриваемой жидкости: (1р01 < 0. Наложение же сколь угодного
слабого сдвига скорости «дестабилизирует» течение в смысле нарушения знакоопределенности подынтегрального выражения (8). В результате условий устойчивости точных стационарных решений (4) относительно малых плоских возмущений нет. Особо следует заметить, что в коэффициентах формы подынтегрального выражения (8) локальное число Ричардсона Ш = -^ /р0 )(йр0 /йу)(йи / йу)-2 не появляется. Этотфакт говорит о том, что локальное число Ричардсона как критерий устойчивости из энергетических соображений не возникает [5].
Таким образом, осталось исследовать неустойчивость стационарных решений (4) смешанной задачи (1) по отношению к малым плоским возмущениям.
Ниже интеграл Е (8) будет использоваться в форме
Е = Т + Т + П = const,
(9)
где
(10)
Для произвольной постоянной величины Я заметим, что функционал
2Т,= 2Т-Я — + Я2 М =
Я &
= Ял, (( + и£и-Я£)2+(2, + и^2х - )2 \lydx > 0 (12)
-« 0
неотрицателен. Используя соотношение (11) для функционала Т и для положительной величины
а = g max
OSySH
dp0 / dy
Po
из неравенства (12) можно
извлечь дифференциальное неравенство й 2Ш . йШ .2
—— - 2X--------+ 2(X + а)Ш >0 . (13)
В самом деле,
й~М - 2Л— + 2Л2М > -4П > - 2а(р0(£ )2йхйу > йт2 йт
— 2а| р0 (%2 )2 йхйу - 2а| р0 (^ )2 йхйу = -2аМ (14)
Откуда и вытекает искомое дифференциальное неравенство (13).
Оказывается, если постоянную величину X подчинить требованию Л > 0 и добавить дополнительные условия [10] к неравенству (13)
М (п /2 4 М2 +2а) >0. п = 0,1,2,... ;
dm
dt
т=1 +як ( + и£1х )2 +р0 ( + и%2 х )2 \lydx >0,
2 -« о
Т1 =~ Я \р0и (бЛ + #2х#2, ) + Рои 2 (( )2 + (#2х )2 )]№
-~0
П = -- \\р0йуйх.
2 0
В интересах последующего изучения удобно ввести в рассмотрение вспомогательный функционал вида
(ли/24 Л2 + 2а) > 2(Л + а / Л)М (ли/2-Л2 + 2а) ;
М (ли /2-2Я2 + 2 а )=М (0 )ехр(пг/1/2 л/22 + 2 а) , (15)
(ли / 2 •/Л2 + 2а) = (0)екр(лпЛ / 2 л/Л2 + 2а) ;
М (0) >0, 2м (0) > 2(1 + а/ Х)Ы (0) ,
йг
то из него в силу этих условий будет вытекать искомая априорная экспоненциальная по времени оценка снизу вида , ч ,
М (т ) > С ехр(Д/), (16)
где С - известная положительная постоянная величина.
Действительно, неравенство (13) может быть формально проинтегрир ов ано на полуинтервалах
г п < к п/( 2дД2 + 2а) +
Если дважды продифференцировать функционал М (10) по его аргументу ? и сделать несколько преобразований получившегося в результате интеграла, используя связи (6), (7) и (9), то несложно прийти к вириальному равенству в виде
= 4 (Т - П).
Л Т ) (11)
( tn = 2жп/^Л2 + 2а; п = 0,1, 2, ... ), дпя че]^о над10]вып0лнить несколько упрощающих замен искомого функционала М:
а) Мj (í ) = exp (- At)M(t): M"(t) + (Л2 + 2а)М( > 0;
б) M 2 (t) = M j (t У cos (t ^ + 2a):
0
М 2 (t )cos (t-j % + 2аf (t) -^Л2 + 2а М 2 (t )х
х sin (д/ Л2 + 2а )> 0;
в) М3(t) = М2(t)cos2 (tjА2 + 2a) M'(t)> 0.
Интегрирование последнего неравенства и
выполнение обратных замен приводят к неравенству
М (t )> (cin cos( t^j % + 2a)+
+ C2n sin ( t4Á2 + 2a ))exp (Át), (17)
где С1n и C2n - произвольные постоянные величины.
Учитывая нестрогость оценки (17), постоянные
C1n и C2n нетрудно связать со значениями функцион)1ла м: (10) и его первой производной М'(t) в моменты времени tn ( п = 0,1,2,...). В итоге соотношение (17) окончательно может быть записано в виде
м (Р)^/(0; (18)
f (f) ее (м (tn )cos( t4 Л2 + 2a) +
+ [M \tn )-ЛЫ (tn)] sin (t^jЛ2 + 2a)/4 Л + 2a )x x e^гpЛ(t -O-
Для того чтобы обосновать процедуру интегрирования нерав е нства (13) на промежутках tn < t <п/(2^Л2 + 2а) + tn ( п = 0,1,2,...), приведшую в результате к нюкней оценке (17), нужно вычислить производную первого порадка функции f по её аргумент) t:
f '(t) = (m \tn )cos (t^/A2 + 2a) +
+ Л [m '(. )-vm (tn + 2a - 4Л2 + 2aM (tn ))x
x sin(f ^Л2 + 2a ))exp(A(f-fn)). (i 9)
П]рин имая в о в ниман ие соотно шен ия (18) и (19), можно сдеаать заключение, что функция f (t) будет положительной а строго возрастаюше й на полуинтервалах tN< t<n(2.4Л2 + 2а + tn( п = 0,1,2,... )
в случае, всли истинны неравенства
М (t„ ) > Ia, М '(tn) > 2{Л + alЛ)м (tn ). (20)
Эти оеоавенства как раз и служат требуемыми гарантами ()равомернорти осуществлённой вышр про(1е(ауры интегрирования соотношения (13).
Поскольку промежутки tn <t < я/ 2т) Л2 + 2а + tn ( п = 0, 1,2,...) друг с другом не птресекаотсяб значения функционала М и его первой произ)одной M'(t) на левых концам настоящих промежутков могут задаваться любыми, беро всяких ограничений. В частно сти, их можно взять в форме
М {tn) = M (°)exp(Aí„), М\tn) = М /(°)exp(2T„).
Тогда нер авенства (2 0) будут выполнены, если М(0)>0, М'(о)> 2(Л +а) Л)М (&).
Функция же f (t) предстанев в виде
f{t)=(M (0) cos ( 14 л2 + 2а) + [M '(0)-
- ЛИ(0)]sin (t*JЛ2 + 2а)/д¡Л2 + 2а)ехр(Лг).
Аналогичные рассуждения могут быть проведены и в том случае, когда соотношение (13) надо будет интегрировать на остальных временных полуинтервалах. Учитывая этот факт, далее итоги интегрирования неравенства (13) на промежутках
Tkn ^ Т <П W Л + 2а +
(tkn =nk12л[Лг+2а + 2лп/л]Лг+2а;
к = 1,2,3; п = 0,1,2,...
представляются в форме кратких иллюстрирующих выкладок, без подробных комментариев:
1) Мt (t) =exp( — At)M ( t): MД?) + (Л2 + 2a)M1 > 0;
2) M 2 (t ) = Mj (t)/cos(^Л2 + 2a):
M '2{t )cos(^A2 + 2a)] (t ) --sJá2 + 2a M '2{t )x x sin (t^j Л2 + 2a) > 0;
3) M 3 (t ) = M ¡{t )cos2(^ A2 + 2a):
M'3)t)<o )k = 1,2), M'3)t)> 0 )k ^3)
4) M (t) > if3n cos(I2 + 2a) + C4n sin (t (/)2 + 2a) )sxp (At);
С, , С. — const;
3n ’ 4n
5) M (l) >ft (t )n
а) ,
fi() = \M(fin)sin{tjxr+2a) - [м'(tln)-.
1 - AM Ptu )]cos(^A2 + 2a)s¡A2 + 2a)x x exp(l(i - tj),
f1(t) = [m'(tln )sin (tVA2 + 2a) - {a[m \tln) -
— vM (tln)]/ A2 + 2a — ^ A2 + 2a M (tln )]x x co s (t(A2 + 2a))exp(A(t - tu));
б)
f2 [t) = -(M (t 2n )cos (tyl A2 + 2a) + [M '{t 2n)
-AM[í2и)]sin (t^¡I2 + 2a^y¡A2 + 2a)x
xexp(A(t -t2n))
f',(t) = - (м '(и, )cos (t4 X + la) + (x[m '(t, п) -- АИ (t2n )]/д/А2 + 2a -^A2 + 2a x x M (t 2л ))sin (t^j A2 + 2a) )exp(A(t -t2n));
f3[t )=[-M ХЧп)sin {t^[X~+2a) + [a\m '(t3n) -) - vM (t3n )]/л]Л2 + 2a -t]A2 + 2a x x M (t 3n ))cos [t^j 22 + 2a)]exp)2([ - ^J);
6)M(ife )>0, M '(^ )>2(А + а/А)м(ты);
7) M (tM)E) M (o)ex p Л , M \tl() s M '(0 )ехрЛ(ы;
M (0) >0, M '(0)> 2{А + а/ Л )m(0);
F (t) = (m (0)sin (tjI2 + 2a) - [m'(0) -
- vM (0)]c( os (t^jA2 + 2a)/iJA2 + 2a)exp(At),
f,_ (t )=-(m (0)c os (ty[¥+2a) + [M Xo)-
- vM(o)jsin(^A + 2a)/*JA2 + 2a)exp(Ai),
f3 (t) = (- M (0)sin (tj A2 + 2a) + [M '(0) -1 - vM (0)]cos (t*JA2 + 2a)/-yjA2 + 2a)exp(At).
Если проанализировать финальные выражения для функций /( ), f (t) (к = 1, 2, 3), то несложно уви-
деть, что графиками данных функций на отвечающих им полуинтервалах времени будут являться кривые, которые лежат поперёк полуполосы, экспоненциально быстро уходящей на бесконечность, причём их левые концы опираются сверху на нижнюю границу этой полуполосы
g(t) = М (0)exp(li),
а правые примыкают к её верхней границе
g2 (t) = [m '(0) - AM (o)]exp(At)/д/A2 + 2a.
Что говорит о том, что интеграл М (10) растет со временем, по крайнмй мере, не медленнее, чем экспоненциально. Тем самым продемонстрировано, что при наличии условий (15) из соотношения (14) действительно вытекает априорная экспоненциальная оценка снизу (16).
Таким обр(зо м, путем выбора специальных начальных утлотий на левых концах рассматриваемых временных промежутков (третье и четвертое выражения (15)) удалось доказать существование
единых начальных данных (два последние неравенства в системе соотношений (15)) для малых плоских возмущений (6), (7) точных стационарных решений (4) смешанной задачи (1), которые обеспечивают справедливость условий положительности и строгого возрастания функций/^ ), /) ) , где (к =1,2,3 ), (первые два неравенства из системы соотношений
(15)) на всех рассматриваемых полуинтервалах времени. Значит, согласно определению неустойчивости по А. М. Ляпунову на полубесконечном временном промежутке решения системы дифференциальных уравнений [10], этим показана принципиальная возможность возникновения и последующей эволюции во времени неограниченно нарастающих малых плоских возмущений (6), (7) точных стационарных решений (4) начально-краевой задачи (1).
Так как соотношение (16) построено без предъявления каких бы то ни было требований ограничительного характера к установившимся плоско-параллельным течениям (4), то теоретическая линейная неустойчивость последних течений относительно малых плоских возмущений (6) и (7) будет абсолютная. Кстати, для экспоненциально растущих во времени малых плоских возмущений (6), (7) счетные наборы условий (16) удовлетворяются тождественно. Более того, два первых неравенства в системе соотношений (15) будут желаемые достаточные условия практической линейной неустойчивости установившихся течений (4) стратифицированной идеально проводящей жидкости относительно малых плоских возмущений (6), (7). Следует отметить,
что достаточные условия практической линейной неустойчивости носят конструктивный характер, поскольку их истинность может быть проверена как в физических экспериментах, так и в численных экспериментах.
Отметим, что именно интеграл М (16) является искомым функционалом А. М. Ляпунова, который нарастает во времени в силу уравнений смешанной задачи (6), (7). Характерной чертой данного роста является большой произвол постоянной величины в показателе экспоненты из правой части неравенства
(16). Данное обстоятельство позволяет интерпретировать любое решение начально-краевой задачи (6), (7), (15), которое нарастает со временем согласно найденной априорной экспоненциальной оценке снизу (16), как аналог примера некорректности по Адамару [11].
Наконец, детально описанная выше процедура интегрирования соотношения (13) наглядно демонстрирует, что сведения о начальных условиях (15) для растущих во времени малых плоских возмущений (6), (7) могут быть извлечены и при изучении класса кусочно-гладких функций.
В данной работе рассмотрена задача линейной
устойчивости установившихся плоско-параллельных сдвиговых течений неоднородной по плотности невязкой несжимаемой жидкости в поле силы тяжести в прослойке между двумя покоящимися непроницаемыми твердыми параллельными неограниченными поверхностями.
Прямым методом А. М. Ляпунова продемонстрировано, что данные течения абсолютно неустойчивы к малым плоским возмущениям, а также получены достаточные условия практической неустойчивости этих течений относительно малых плоских возмущений, которые носят конструктивных характер, и их можно проверять непосредственно в ходе проведения численных расчетов и физических экспериментов. Построена априорная нижняя оценка, свидетельствующая об экспоненциальном во времени росте исследуемых малых возмущений, при этом инкремент представляет собой некий положительный параметр произвольного характера.
Следует подчеркнуть, что с математической точки зрения представленные в данной статье результаты служат, в большинстве своем, априорными, поскольку теоремы существования решений изучавшихся смешанных задач для систем дифференциальных уравнений с частными производными не доказаны.
Л и т е р а т у р а
1. Dikii L. A. Гидродинамическая устойчивость и динамика атмосферы. Гидрометеоиздат. - Leningrad. 1976. - 106 с.
2. Howard L. N. Note on a paper of John Miles // J. Fluid. Mech. - 1961, - V. 10, - № 4, - P. 509-512.
3. Miles J. W. On the stability of heterogeneous shear flows // J. Fluid. Mech. - 1961. - V. 10. - № 4. - P. 496-508.
4. Дразин Ф. Введение в теорию гидродинамической устойчивости. Пер. с англ. Г. Г. Цыпкина; Под ред. А. Т. Ильича. - М.: Физматлит., 2005. - 288 с. ISBN 5-9221-0629-5.
5. Владимиров В. А. Об интегралах плоских движений идеальной несжимаемой неоднородной по плотности жидкости // Механика жидкости и газа. - 1987. - № 3, - C. 16-20.
6. Gubarev Yu. G. The development of Lyapunov’s direct method in the application to new types of problems of hydrodynamic stability theory // In: Progress in nonlinear analysis research / Ed. Erik T. Hoffmann. Chapter 7. New York: Nova science publishers, inc., - 2009. - P. 137-181. (ISBN 978-1-60456359-7).
7. Lyapunov A. M. The general problem of the stability of motion. Taylor&Francis. - London. 1992. - P. 242.
8. Chetaev N. G. Stability of motion. Nauka. - Moscow. 1990. - P. 207.
9. S. Chandrasekhar: Ellipsoidal figures of equilibrium. - Yale University Press, New Haven. 1969. - P. 264.
10. Демидович Б. П. Лекции по математической теории устойчивости. - М.: Наука, 1967. - 472 с.
11. Годунов С. К. Уравнения математической физики. - М.: Наука, 1979. - 392 с.
