Radiative instability of a barotropic jet flow in a stratified rotating atmosphere

Cover Page

Abstract


The problem of the stability of a jet flow with a piecewise linear velocity profile in a stratified rotating atmosphere is considered. The linearized system of equations for perturbations is reduced to a single equation with respect to the amplitude of the longitudinal velocity component containing the turning points. In terms of Airy functions, an asymptotic solution of the equation is constructed that is valid for small values of the Rossby number. It is shown that the flow becomes unstable due to radiation of inertial-gravitational waves. An analytical expression is obtained for the growth rate of perturbations.


Full Text

ВВЕДЕНИЕ

Под излучательной неустойчивостью понимают неустойчивость гидродинамических течений, обусловленную излучением волн, уносящих энергию на бесконечность [1]. Этот тип неустойчивости достаточно широко исследовался как в теоретической акустике [2–6], так и в гидродинамике стратифицированной жидкости [7–10]. Основное внимание при этом уделялось рассмотрению локализованных осесимметричных и эллиптических вихрей. Струйные течения, формирующиеся на скачках потенциальной завихренности (PV) во вращающемся слое мелкой воды, рассматривались в работах [11, 12]. Волновые возмущения этих течений представляют собой поверхностные гравитационные волны и так называемые краевые квазигеострофические волны Россби, связанные со скачком PV. К излучательной неустойчивости течений приводит взаимодействие волн двух классов [12].

Описание излучательной неустойчивости требует сращивания асимптотических решений уравнений динамики и, в общем случае, достаточно громоздко. В настоящей работе, в рамках модели стратифицированной вращающейся атмосферы, рассмотрен пример излучательной неустойчивости, допускающий простое аналитическое исследование. В этом примере рассматривается двумерное струйное течение с кусочно-линейным профилем скорости. Линеаризованная система уравнений для возмущений течения сводится к одному дифференциальному уравнению, содержащему точки поворота. С использованием равномерных асимптотик решений этого уравнения показано, что течение, которое в модели несжимаемой однородной среды часто называют stable jet [13], становится неустойчивым за счет излучения внутренних инерционно-гравитационных волн (ИГВ). Применительно к излучению акустических волн аналогичный результат был получен в недавней работе [14].

Следует отметить, что практический интерес к излучательной неустойчивости стимулирован поиском динамических механизмов, обеспечивающих непрерывную (спонтанную) генерацию волн. С этими механизмами связывают наблюдаемые волновые проявления от атмосферных фронтов и струйных течений [15]. В последние годы активно исследовались два таких механизма — механизм генерации волн вихревыми возмущениями в сдвиговых течениях [16–18] и механизм нелинейного геострофического приспособления [19, 20]. Излучательная неустойчивость дополняет эти механизмы.

Отметим также, что, в теоретическом плане, излучательная неустойчивость тесно связана с концепцией о разделении движений стратифицированной вращающейся жидкости на быстрые волновые движения (ИГВ) и медленные вихревые (квазигеострофические) движения. Асимптотический анализ уравнений динамики показывает, что такое разделение является неполным даже при малых значениях числа Россби [21–24]. Этот факт приводит к спонтанному (вторичному) излучению волн в ходе медленной эволюции квазигеострофических движений.

ОСНОВНЫЕ УРАВНЕНИЯ И УСТОЙЧИВОСТЬ В КВАЗИГЕОСТРОФИЧЕСКОМ ПРИБЛИЖЕНИИ

Рассматривается горизонтальный слой стратифицированной вращающейся атмосферы высоты H с постоянными значениями частоты плавучести N и инерционной частоты f. Предполагается, что в этом слое имеется баротропное течение U=U0U(y/L) с характерной скоростью U0 и горизонтальным масштабом L. В приближении гидростатики поведение малых возмущений течения описывается системой уравнений

 

 Du/DtfU'(y)v=px,Dv/Dt+fu=py,   pz=σ,  Dσ/Dt+N2w=0,ux+vy+wz=0.(1)

 

Здесь u, v, w — компоненты скорости вдоль горизонтальных осей x, y и направленной вертикально вверх оси z соответственно, p — возмущение нормированного на среднюю плотность давления, σ — плавучесть, D/Dt=/t+U(y)/x — оператор полной производной. В отсутствие фонового течения система (1) описывает динамику длинных инерционно-гравитационных волн [13].

При задании условия w = 0 на твердых границах слоя z = 0, H, система (1) допускает разделение переменных — решение представляется в виде рядов по собственным функциям (модам) оператора d2/dz2:

u,v,p=n=1un,vn,pncosnπHz,        w,σ=n=1wn,σnsinnπHz.

Подставляя эти ряды wn, σn в (1) и исключая wn,σn, для каждой вертикальной моды получим систему уравнений (нижний индекс для краткости опускаем)

Du/DtfU'(y)v=px,Dv/Dt+fu=py,Dp/Dt+a2(ux+vy)=0.(2)

В этой системе, аналогичной линеаризованной системе уравнений мелкой воды, a=NH/πn — фазовая скорость n-й моды гравитационных волн в отсутствие вращения [13].

В анализе удобно использовать безразмерную форму (2). Принимая в качестве масштабов координат, времени, скорости и давления соответственно L, L/U0, U0, fU0L, из (2) получим систему

εDu/Dt1εU'(y)v=pxεDv/Dt+u=py,εF2Dp/Dt+ux+vy=0,, (3)

содержащую два безразмерных параметра: ε=U0/fL — число Россби, F=L/LR — отношение горизонтального масштаба течения к радиусу деформации Россби LR=a/f. Из системы (3) следует уравнение переноса потенциальной завихренности (PV) для возмущений

Dq/Dt=U''(y)v, q=vxuyF2(1εU'(y))p, (4)

играющее важную роль в последующем рассмотрении.

Далее будем рассматривать кусочно-линейный профиль скорости

U(y)=1y, (5)

описывающий неограниченное струйное течение с разрывом завихренности при y = 0 (размерная форма профиля U(y)=U0(1y/L)). Для профиля (5) U'(y)=sgn(y), U''(y)=2δ(y), где sgn(y), δ(y) функция выделения знака и дельта функция соответственно. С учетом второго равенства, уравнению (4) удовлетворяют возмущения с нулевой PV

vxuyF2(1εU'(y))p=0,y0. (6)

Для таких возмущений система уравнений (3) сводится к одному уравнению относительно продольного компонента скорости

(εF)2D2u/Dt2+F21εU'(y)u=uxx+uyy. (7)

Это уравнение получается применением оператора D/Dt к первому уравнению системы с последующим использованием (6). Отметим, что входящий в (7) безразмерный параметр εF=U0/a является некоторым аналогом числа Маха (отношение характерной скорости течения к скорости ИГВ).

Прежде чем переходить к дальнейшему анализу, исследуем устойчивость течения (5) в квазигеострофическом приближении ε<<1, когда асимптотически u=pyv=px и PV q=qg=pxx+pyyF2p. Динамика возмущений при этом описывается уравнением (4). Полагая в (4) q = qg и отыскивая решения в форме нормальных мод p=p~(y)eik(xct), получим уравнение

Ucp~yyκ2p~=2δ(y)p~, κ=k2+F2, (8)

где k — волновое число, c — комплексная фазовая скорость. Из (8) следует, что всюду вне линии разрыва p~yyκ2p~=0, откуда p~=rexpκy, p~y=rsgn(y)κexpκy, где r — произвольная постоянная. Для определения c используется условие на разрыве, которое получается интегрированием (8) по малой окрестности точки y = 0:

(1c)p~y(+0)p~y(0)=2p~(0). (9)

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

c=1κ1=1k2+F21/2, (10)

свидетельствующее об устойчивости струйного течения в рамках квазигеострофической модели. Возмущение давления, p=rexpκyexpik(xct), играющее роль функции тока возмущения, списывает прижатую к скачку завихренности волну, распространяющуюся с фазовой скоростью c (10). В квазигеострофической теории подобные волны часто называют краевыми волнами Россби [13]. Ниже будет показано, что за счет излучения ИГВ, амплитуда этих волн медленно нарастает, т. е. возникает неустойчивость. В математическом плане эта неустойчивость обусловлена появлением комплексной поправки к скорости (10) с положительной мнимой частью.

АМПЛИТУДНОЕ УРАВНЕНИЕ И УСЛОВИЕ ИЗЛУЧЕНИЯ

Рассмотрим теперь динамику возмущений в рамках полной (нефильтрованной) системы (3), (4). Будем представлять в возмущения форме нормальных мод

(u,v,p)=u~(y),v~(y),p~(y)eik(xct), (11)

где в круглых скобках правой части стоят амплитуды, зависящие только от y. Система уравнений для амплитуд сводится к одному уравнению относительно амплитуды продольной компоненты скорости

u~yy+κ2(εF~)2(yc~)21+εk2F~2sgn(y)u~=0, (12)

где обозначено c~=1c, κ=k2+F2F~=Fk/κ. Данное уравнение прямо следует из уравнения (7). Отметим, что параметр c~ есть разность осевой скорости течения U(0) = 1 и фазовой скорости c. При ε<<1 и вблизи линии разрыва (12) сводится к уравнению u~yyκ2u~=0, которое удовлетворяется в квазигеострофической теории.

Амплитуды поперечной компоненты скорости и возмущения давления связаны с амплитудой u~ диагностическими (поляризационными) соотношениями

v~=ikΔ1u~yεF2(1εU')(Uc)u~,p~=Δ1(1εU')u~y+εk2(Uc)u~,Δ=k2+F2(1εU')2, (13)

где, напомним, U'(y)=sgn(y). Соотношения (13) следуют из амплитудной формы первого уравнения системы (3) и уравнения (6)

ikε(Uc)u~(1εU'(y))+ikp~=0,ik,v~u~y+F21εU'(y)p~=0. (14)

Рассматривая (14) как систему линейных алгебраических уравнений относительно v~,p~, получим соотношения (13).

Для построения решения однородного уравнения (12) к нему нужно присоединить условия на разрыве y = 0 и условия на бесконечности. Условия на разрыве завихренности сводятся к непрерывности нормальной компоненты скорости и давления при переходе через разрыв: p~(+0)=p~(0), v~(+0)=v~(0). С учетом (13) эти условия дают соотношения, связывающие u~ и u~y по обе стороны от разрыва.

На бесконечности будем требовать выполнения принципиально важного условия излучения. Оно состоит в том, что поток волновой энергии Π=vp (в системе координат, связанной с течением) должен быть направлен от линии разрыва, т. е. быть положительным при y+, и отрицательным при y. Здесь угловыми скобками отмечена операция горизонтального осреднения ϕ=limL12LLLϕdx. Для математической формулировки этого условия представим c=cr+ici, и учтем, что вещественная часть решения для нормальной моды записывается в виде p=ekcitp~eiθ+p~eiθ, θ=k(xcrt), где звездочкой обозначено комплексное сопряжение (аналогично для остальных переменных). При этом, vp=e2kcitp~+p~=2e2kcitRep~, где символами Re и Im далее обозначаются действительная и мнимая части выражений. Используя соотношения (13) можно показать, что Rep~=εkΔ1Reic~yu~u~y. С учетом, Re(iz)=Im(z), для потока энергии на бесконечности асимптотически получим

Π=vp=2εkΔ1ye2kcitW, W=Imu~u~y, y. (15)

Непосредственно из (15) следует, что условию излучения удовлетворяют комплексные решения уравнения (12), для которых W < 0 при y+. и W > 0 при y-, Обратим внимание, что при y в уравнении (12) можно пренебречь c~ и величина W играет роль вронскиана уравнения с вещественными коэффициентами. Согласно общей теории [25], для таких уравнений W = const.

Уравнение (12) вместе с условием излучения и условиями на разрыве определяет спектральную задачу для нахождения комплексной фазовой скорости c=cr+ici. Течение неустойчиво, если мнимая часть ci>0.

ИЗЛУЧАТЕЛЬНАЯ НЕУСТОЙЧИВОСТЬ ПРИ МАЛЫХ ЧИСЛАХ РОССБИ

Замена переменной ξ=εF~y, F~=Fk/κ, приводит уравнение (12) к виду

 u~ξξ+λ2(ξεF~c~)21εk2F~2U'u~=0,λ2=(κ/εF~)2. (16)

Асимптотическая форма (16) при ε<<1 сводится к уравнению

u~ξξ+λ2ξ21u~=0, (17)

которое не содержит спектральный параметр c~=1c. Этот факт существенно упрощает решение спектральной задачи — после нахождения решения уравнения (17), удовлетворяющего условию излучения, параметр c~ находится из условий на разрыве. Асимптотическую форму этих условий проще всего получить, привлекая амплитудную форму уравнения для PV (4)

iku¯c)ikv~u~y+F21εU'(y)p~=2δ(y)v~. (18)

Интегрирование (18) по малой окрестности разрыва получим

ik(1c)u~(+0)u~(0)=2v~(0). (19)

Заметим далее, что уравнение (17) инвариантно относительно замены yy и при выполнении условия излучения при y+, оно автоматически удовлетворяется при y, если использовать правило сопряжения решений u~(y)=u~+(y). Здесь через u~+(y)u(y) обозначены решения, соответственно, в областях y > 0, y < 0. С учетом указанного правила, условия на разрыве сводятся к условию u~(+0)=u~(0) и условию     uξ(+0)=c~κλu(+0),  c~=1c.(20)

Последнее условие следует из (19) с использованием вытекающего из (13) асимптотического соотношения v~=ikκ1u~y. Таким образом, приходим к следующей редуцированной спектральной задаче — на полуоси y > 0 найти значения параметра c~=1c, при котором существуют нетривиальные решения уравнения (17), удовлетворяющие условию (20) и условию излучения W < 0 при y+.

Отметим, что правило сопряжения и условие (20) точно выполняются для решения u~=p~y в квазигеострофическом приближении. Это приближение дает асимптотику точного решения непосредственно вблизи линии разрыва y = 0.

Поскольку спектральная задача (17), (20) аналогична рассмотренной в [14] спектральной задаче для акустических волн, далее ограничимся кратким анализом. На полуоси ξ>0 уравнение (17) содержит точку поворота ξ=1, при переходе через которую монотонный (экспоненциальный) характер решения сменяется осциллирующим волновым поведением. Методы построения асимптотических решений уравнений с точкой поворота в случае λ>>1 (большой параметр) описаны в [25, 26]. Соответствующие решения представляются в терминах функций Эйри Ai(z),   Bi(z). Для уравнения (17) в качестве линейно независимых решений удобно использовать функции Φ1(z)=Bi(z)+iAi(z), Φ2(z)=Φ1(z)=Bi(z)iAi(z).

Функция Φ1(z) имеет асимптотические представления при z+

Φ1(z)=1πz1/4expi23z3/2+π4,Φ1(z)=1πz1/4exp23z3/2+i2exp23z3/2, (21)

вытекающие из известных асимптотических представлений для функций Эйри.

Асимптотическое разложение решения уравнения (17), равномерно пригодное для всех ξ>0, включая точку поворота ξ=1, дается выражением [26]

u~ξ=1s'ξC1Φ1λ2/3sξ+C2Φ2λ2/3sξξ1:       23s(ξ)3/2=1ξτ21dτ,ξ1:       23s(ξ)3/2=ξ11τ2dτ, (22)

где C1, C2 — произвольные постоянные. Чтобы удовлетворить условию излучения в этом выражении следует положить C2 = 0 (см. ниже). С использованием асимптотических представлений (21), из (22) найдем асимптотики решения в областях, разделенных точкой поворота

u~(ξ)=Ceiπ/4ξ214expiλα(ξ),C=C1π1/2λ1/6,   ξ>1, (23)

u~(ξ)=C1ξ24exp(λβ(ξ))+i2exp(λβ(ξ)),0<ξ<1. (24)

Здесь функции α(ξ),  β(ξ) определены выражениями

α(ξ)=12ξξ21lnξ+ξ21,β(ξ)=π412arcsinξ+ξ1ξ2, (25)

полученными вычислением входящих в (22) интегралов.

Из (23) следует асимптотика u~(ξ)=Ceiπ/4ξ1/2expiλξ2/2,  справедливая при ξ. Используя эту асимптотику легко показать, что при ξ вронскиан W=Imu~u~ξ=λC2<0, т. е. выполнено условие излучения на бесконечности. Для решения, отвечающего C20, вронскиан W > 0. По этой причине выше было положено C2 = 0.

Разлагая β(ξ) по степеням ξ, из (24) получим асимптотику при ξ0+0

u~(ξ)=Ceλπ/4eλξ+(i/2)eλπ/2eλξ. (26)

Определяя отсюда u~(+0),  u~ξ(+0) и подставляя в условие (20), найдем значение спектрального параметра

c~=κ1  1(i/2)eλπ/21+(i/2)eλπ/2κ1  1mieλπ/21+m, (27)

где m=(1/4)exp(λπ). Соответственно, для фазовой скорости c=1c~ асимптотически получим

c=cr+ici=1κ1+iκ1eλπ/2. (28)

Поскольку мнимая часть ci=κ1eλπ/2>0, имеет место неустойчивость с экспоненциально малым инкрементом нарастания

γ=kci=(k/κ)eλπ/2. (29)

Таким образом, излучательная неустойчивость при малых числах Россби развивается крайне медленно. Отметим, что инкремент (29) зависит от волнового числа и γ(k)0 при k0 и k.Максимум инкремента асимптотически достигается при k=km=F, соответствующее максимальное значение γm=21/2exp(πε1).

Остановимся на некоторых особенностях излучаемого волнового поля в дальней волновой зоне ξ+ (на больших расстояниях от линии разрыва). С учетом (23), для продольного компонента u=u~(y)eik(xct) скорости после выделения вещественной части получим

u=Ceγtξ1/2cosθ,  θ=k(xcrt)+λ(ξ2/2)+(π/4).(30)

Согласно (30), осцилляции поля u по координате ξ носят негармонический характер — локальная длина волны по этой координате стремится к нулю. Частота временных осцилляций

ωR=kcr=k1k2+F21/2, (31)

зависит от волнового числа и определяется квазигеострофическим решением для краевой волны Россби. В фиксированный момент времени с ростом ξ продольный компонент затухает как 1/ξ. Согласно соотношениям (13), поперечный компонент v (вместе с возмущением давления) пропорционален ξ, т. е. нарастает по корневому закону. Этот рост обеспечивает постоянство плотности потока импульса I=uv.

Теневая картина изолиний полного поля продольной скорости (в фиксированный момент времени) аналогична построенной в [14] картине для акустических волн (см. рис. 2 в [14]). На небольших расстояниях от линии разрыва имеем вихревую картину с замкнутыми изолиниями, определяемыми квазигеострофическим решением («кошачий глаз»). С ростом ξ эта картина трансформируется в волновую картину поля (30). Изолинии этого поля асимптотически совпадают с линиями равных фаз θ=const (линиями волновых гребней и ложбин). Обозначая x1=xcrt, для линий равных фаз из (30) получим уравнения x1=k1(θπ/4)0.5εFy2, описывающие однопараметрическое семейство парабол с верхней ветвью, наклоненной влево (по потоку).

Завершим данный раздел важным замечанием. Из представленного анализа следует, что необходимое условие излучательной неустойчивости состоит в наличии точки поворота у амплитудного уравнения. В рассмотренном случае малых значений числа Россби Ro=ε<<1, эта точка (ξ=1) представляет собой уровень, на котором частота ωR (31) волны Россби становится равной частоте ωg инерционно-гравитационных волн. Действительно, с учетом доплеровской поправки в размерных переменных (отмеченных звездой) частота длинных ИГВ определяется выражением [13]

ωg=kU(y)±a2(k2+l2)+f21/2. (32)

Нормируя (32) на U0/L, полагая l=0 и переходя к безразмерному волновому числу k=kL, для безразмерной частоты (32) получим

ωg=k1y±(κ/εFk). (33)

Приравнивая (33) (со знаком плюс) выражению (31), асимптотически найдем y=κ/εFk или ξ=1. Подчеркнем, что совпадение частот ωR=ωg происходит только за счет доплеровской поправки, приводящей к зависимости частоты ИГВ от пространственной координаты.

Сделанное замечание, в сочетании с дисперсионными соотношениями (31), (33), позволяет предположить, что для рассматриваемого неограниченного течения излучательная неустойчивость (совпадение или синхронизация частот ωR=ωg) будет иметь место и при произвольных значениях числа Россби. В этом состоит принципиальное отличие от рассмотренного в [12] ограниченного струйного течения на мелкой воде (с экспоненциальным профилем), для которого синхронизация частот ωR=ωg происходит только при Ro > 1.

ЧИСЛЕННЫЕ ОЦЕНКИ

Приведем некоторые численные оценки для слоя атмосферы толщиной H = 10 км с характерными частотами N=102c1, f=104c1.  Для первой вертикальной моды (n = 1) фазовая скорость a=NH/π=31.8  м/с, бароклинный радиус деформации Россби LR=a/f=318 км. Если рассматривать струйное течение с характерной скоростью U0=10ì/c и горизонтальным масштабом L=LR, то число Россби ε=U0/fL=U0/a=0.31 и параметр F = 1. Соответствующий адвективный временной масштаб T0=L/U0=8.33  ÷, инерционный масштаб f1=2.8  ÷ч. Число Россби, очевидно, представляет собой отношение этих масштабов. Для приведенных масштабов характеристики неустойчивых волновых возмущений определяются значением безразмерного волнового числа. При k = 5 и k = 10 (длины волн соответственно 496 и 198 км) размерный период колебаний T=2πT0/ω составляет соответственно 13 ч и 5.7 ч. С ростом k этот период неограниченно уменьшается. Инкременту нарастания возмущений (35) отвечает временной масштаб T=T0/γm, который многократно превышает адвективный временной масштаб T0. Для наименьшего значения этого масштаба T=T0/γm, достигаемого при k = 1, превышение составляет пять порядков. Отметим, что аналогичный результат был получен и в работе [12] — для ограниченного струйного течения на мелкой воде время нарастания возмущений многократно превышает период вращения системы. Приведем также оценку координаты точки поворота y=κ/εFk или, в размерном виде y=(LR/ε)(κ/k)., При k = 10 и F = 1 получим y1000​ км, т. е. волновая картина начинает проявляться на расстоянии тысячи километров от оси струи.

Подчеркнем, что приведенные оценки носят иллюстративный характер и, в зависимости от значений внешних параметров, могут существенно меняться. В частности, для океана с H = 5 км, N=103c1, f=104c1 фазовая скорость первой бароклинной моды a=NH/π=1.6м/c, бароклинный радиус деформации LR=a/f=50км. При L=LR и U0=0.3  м/c число Россби ε=0.18. Значения этих параметров отличаются от атмосферных значений, что приводит к отличию характеристик неустойчивых возмущений.

ЗАКЛЮЧЕНИЕ

В исследованиях устойчивости геофизических течений, как правило, используется квазигеострофическое приближение, фильтрующее быстрые волновые движения [13]. При таком рассмотрении упускается из виду неустойчивость, обусловленная излучением волн — излучательная неустойчивость. Как показано в настоящей работе, а также в работах [11, 12, 23, 24], этот тип неустойчивости развивается слишком медленно, чтобы приводить к видимым перестройкам структуры течений. Значимость неустойчивости, однако, состоит в том, что ее развитие сопровождается непрерывной (спонтанной) генерацией волн. Излучательная неустойчивость, таким образом, представляет физический механизм спонтанной генерации волновых движений, который следует учитывать при интерпретации данных дистанционного зондирования атмосферы.

Источник финансирования. Работа выполнена при поддержке РФФИ (проекты 18-05-00831, 18-05-00414).

About the authors

M. V. Kalashnik

Obukhov Institute of Atmospheric Physics, RAS; Schmidt Institute of Physics of the Earth, Russian Academy of Sciences

Author for correspondence.
Email: kalashnik-obn@mail.ru

Russian Federation, Pyzhevsky per., 3, Moscow 109017; B. Gruzinskaya Str., 10, building 1, Moscow, 123242

References

  1. Степанянц Ю.А., Фабрикант А.Л. Распространение волн в сдвиговых потоках. М.: Наука, 1996. 240 с.
  2. Broadbend E., Moore D. Acoustic destabilization of vortices // Phil. Trans. Roy. Soc. 1979. V. A290. P. 353–371.
  3. Копъев В.Ф., Леонтьев Е.А. Об акустической неустойчивости аксиального вихря // Акуст. журн. 1983. Т. 29. № 2. С. 192–198.
  4. Цейтлин В.Ю. Об акустическом излучении распределенных вихревых структур // Акуст. журн. 1988. Т. 34. № 2. С. 321–326.
  5. Гряник В.М. Излучение звука вихрями, акустическая неустойчивость и акустический коллапс // Изв. АН СССР. ФАО. 1988. Т. 24. № 5. С. 513–519.
  6. Zeitlin V. On the backreaction of acoustic radiation for distributed two-dimensional vortex structures // Phys. Fluids. 1991. V. 3. № . 6. P. 1677–1680.
  7. Le Dizès S. Inviscid waves on a Lamb-Oseen vortex in a rotating stratified fluid: consequanced for the elliptic instability // J. Fluid Mech. 2008. V. 597. P. 283–303.
  8. Le Dizès S., Billant P. Radiative instability in stratified vortices // Phys. Fluids. 2009. V. 21. № . 9. 096602.
  9. Schecter D., Montgomery M. Damping and pumping of vortex Rossby waves in a monotonic cyclone: critical layer stirring versus inertia-buoyancy emission // Phys. Fluids. 2004. V. 16. P. 1334–1348.
  10. Schecter D. The spontaneous imbalance of an atmospheric vortices at high Rossby number // J. Atm. Sci. 2008. V. 65. P. 2498–2521.
  11. Boss E., Paldor N., Thompson L. Stability of a potential vorticity front: from quasi-geostrophy to shallow water // J. Fluid Mech. 1996. V. 315. P. 65–84.
  12. Dritshel D.G., Vanneste J. Instability of a shallow-water potential-vorticity front // J. Fluid Mech. 2006. V. 561. P. 237–254.
  13. Vallis G.K. Atmospheric and oceanic fluid dynamics: fundamentals and large-scale circulation. Cambridge: Cambridge University Press, 2006. 738 p.
  14. Калашник М.В. Излучательная неустойчивость неограниченного струйного течения // Акуст. журн. 2018. Т. 64. № 2. С. 139–145.
  15. Plougonven R., Zhang F. Internal gravity waves from atmospheric jets and fronts // Rev. Geophys. 2014. V. 52. P. 33–76.
  16. Калашник М.В. Генерация акустических волн волной завихренности в сдвиговом потоке // Известия РАН. Механика жидкости и газа. 2014. № 4. С. 119–129.
  17. Lott F., Plougonven R., Vanneste J. Gravity waves generated by sheared potential vorticity anomalies // J. Atm. Sci. 2010. V. 67. P. 157–170.
  18. Kalashnik M.V., Chkhetiani O.G. Generation of gravity waves by singular potential-vorticity disturbances in shear flows // J. Atmos. Sci. 2017. V. 74. P. 293–307.
  19. Reznik G.M., Zeitlin V., Ben Jelloul M. Nonlinear theory of geostrophic adjustment. Part I: Rotating shallow water model // J. Fluid Mech. 2001. V. 445. P. 93–120.
  20. Zeitlin V., Reznik G.M., Ben Jelloul M. Nonlinear theory of geostrophic adjustment. Part II: Two-layer and continuosly stratified primitive equations // J. Fluid Mech. 2003. V. 491. P. 207–228.
  21. Ford R. The instability of an axisymmetric vortex with monotonic potential vorticity in rotating shallow water // J. Fluid Mech. 1994. V. 280. P. 303–334.
  22. Ford R., McIntyre M., Norton W. Balance and slow quasimanifold: some explicit results // J. Atm. Sci. 2000. V. 57. P. 1236–1254.
  23. Vanneste J., Yavneh I. Exponential small inertia-gravity waves and the breakdown of quasigeostrophic balance // J. Atm. Sci. 2004. V. 61. P. 211–223.
  24. Zeitlin V. Decoupling of balanced and unbalanced motions and inertia-gravity wave emission: small versus large Rossby number // J. Atm. Sci. 2008. V. 65. P. 3528–3542.
  25. Федорюк М.В. Асимптотические методы для линейных обыкновенных дифференциальных уравнений. М.: Наука, 1984. 344 с.
  26. Найфэ А. Методы возмущений. М.: Мир, 1976. 456 с.

Statistics

Views

Abstract - 122

PDF (Russian) - 101

Cited-By


PlumX

Refbacks

  • There are currently no refbacks.

Copyright (c) 2019 Russian academy of sciences

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies