CC BY
55
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2012 Математика и механика № 1(17)
УДК 532.5:532.517.4
И.В. Ершов, К.И. Зырянов
ДИССИПАЦИЯ ВОЛН КЕЛЬВИНА - ГЕЛЬМГОЛЬЦА В КОЛЕБАТЕЛЬНО НЕРАВНОВЕСНОМ ДВУХАТОМНОМ ГАЗЕ1
На основе уравнений двухтемпературной аэродинамики численно исследовано влияние колебательной релаксации на подавление неустойчивости Кельвина - Гельмгольца в развивающемся во времени сдвиговом слое колебательно неравновесного двухатомного газа.
Ключевые слова: неустойчивость Кельвина - Гельмгольца, колебательная релаксация, кинетическая энергия возмущений, диссипация.
В работе [1] рассматривалось влияние умеренного возбуждения внутренних степеней свободы молекул газа на развитие инерционной неустойчивости Кельвина - Гельмгольца в свободном слое сдвига. Расчеты были выполнены на основе полных уравнений Навье - Стокса сжимаемого теплопроводного газа, в которых с помощью коэффициента объемной вязкости учитывается только возбуждение вращательных степеней свободы молекул газа, а колебательные моды молекул предполагались невозбужденными.
Несмотря на инерционный характер возбуждения волн Кельвина - Гельмгольца результаты отчетливо показали, что с возрастанием объемной вязкости усиливается диссипация кинетической энергии пульсаций, а скорость ее производства уменьшается. Вместе с тем было высказано предположение, что при дополнительном возбуждении колебательных уровней энергии эти эффекты должны усилиться. В этом случае релаксационный процесс может быть использован для управления течением, например для затягивания ламинарно-турбулентного перехода, так как колебательная неравновесность легко создается искусственным путем.
В данной работе нелинейное развитие неустойчивости Кельвина - Гельмгольца рассматривается в рамках уравнений двухтемпературной аэрогидродинамики. В них предполагается, что поступательные и вращательные степени свободы молекул образуют квазиравновесный термостат, характеризуемый статической температурой потока, а релаксация колебательных степеней свободы молекул к равновесию описывается уравнением Ландау - Теллера для колебательной температуры газа [2-4]. В качестве начальных возмущений использовались невязкие моды с максимальным инкрементом нарастания, рассчитанные в рамках линеаризованной системы уравнений двухтемпературной газовой динамики.
Постановка задачи
В координатной плоскости (х1, х2) рассматривается плоскопараллельное сдвиговое течение колебательно неравновесного двухатомного газа, в котором стационарный (несущий) поток равномерной плотности р0 направлен вдоль оси
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (коды проектов 08-01-00116; 11-01-00064).
абсцисс х и имеет профиль скорости с точкой перегиба при х2 = 0: П5 (х2) = и0 Ш(х2 / 50), где параметр 50, имеющий размерность длины, определя-
ется соотношением
Sc - и0 ^dUs/dx2 |x2=0
]-1
и характеризует максимальный наклон профиля скорости.
Считается, что в невозмущенном течении все внутренние степени свободы молекул находятся в равновесии и описываются единой статической температурой потока T0 = const. Отметим, что в отсутствии внешней накачки энергии в колебательную моду молекул это единственно возможное стационарное решение энергетических уравнений используемой системы уравнений двухтемпературной аэрогидродинамики.
В качестве характерных величин для обезразмеривания были выбраны длина 50, асимптотическое значение скорости U0, постоянные плотность р0, температура Т0 и образованные из них характерные время т0 = S0 / U0 и давление p0 =р0 Uq . В обезразмеренных переменных система уравнений двухтемпературной аэрогидродинамики примет вид
др дрщ
С д u■ д u ^
—- + u,—:
dt 1 д хп
v 1У
(дТ дТ Л -------+ и,------
д t 1 д Xi ,
1 У
д t
д p 1 д +
д X,- Re д X,-
д X,
- = 0,
П(Т)
д ut
д X:
+ (Y- 1)р Т
д u.
Y
Re
а, +-
д д X,
П(Т)
д uj
д X
i ]
д X, RePr д X,
П(Т )
дТ
д X,
Y - Р (V)
Y( Y- 1)П(Т )M2
2Re
д u, + д ui д X, д X,
v i , У
+ 21 а, --
дu,
д X,
(1)
Y vP
ГдТ дТ Л
__L + u ____-
. дt J дX, ,
v 1 У
Y а 2 д
Re Pr д x,
П(Т)
дТ,
д X,
Y - Р(Т- - Т)
уМ2р = рТ, п(Т) = Т4/5, уу = ууй/(1 -упЬ), /,] = 1,2,
где по повторяющимся индексам подразумевается суммирование.
В уравнения системы (1) включена зависимость коэффициентов переноса от температуры, что позволило более полно учесть диссипативный эффект «сброса» энергии колебательной моды на поступательные степени свободы молекул, сопровождающийся возрастанием вязкости газа. Выбранная зависимость соответствует условиям относительно холодного несущего потока (мягким потенциалам межмолекулярного взаимодействия). Также принято, что удельные теплоемкости не зависят от статической и колебательной температур потока и постоянны.
Для обезразмеривания коэффициентов переноса использованы их значения при температуре Т0 , отмеченные соответствующим индексом. В коэффициенте
теплоемкости при постоянном объеме выделены парциальные теплоемкости с1г, сго1 и сУ1Ъ, связанные соответственно с поступательными, вращательными и колебательными степенями свободы молекул. Коэффициент теплопроводности также разбит на сумму составляющих Х0 = \ 0 + Xг0 и Xv 0, характеризующих молекулярный перенос тепла поступательными, вращательными и колебательными степенями свободы молекул.
Параметры, входящие в уравнения системы (1), определяются следующим образом. Коэффициенты а1 = Пь 0 / П0 - отношение объемной и сдвиговой вязкостей
и а2 = Xу 0 / Х0; у = (с1г + сго + Я)/(с1г + сго1) = ср /(с1г + сго1) - показатель адиабаты, Я - газовая постоянная; параметр уУ1Ъ = сУ1Ъ /(с1г + сго + сУ1Ъ) равен отношению парциальных теплоемкостей отдельных степеней свободы молекул и косвенно характеризует степень неравновесности колебательной моды молекул, т - время релаксации колебательной моды молекул. Параметры Яе = и080р0 /п0,
М = и0/л1 УЯТ0 и Рг = П0 ср / Х0 - соответственно числа Рейнольдса, Маха и
Прандтля несущего потока.
Система (1) описывает распространенную в аэродинамике ситуацию, когда характерные времена микроскопических процессов энергообмена между различными степенями свободы молекул газа оцениваются системой неравенств [2-4]
т« ~ Тг <<Т™ <<т* ~ т0.
Причем в этом случае поступательные и вращательные степени свободы молекул с малыми соизмеримыми временами релаксации тй ~ тЛ на временах порядка
характерного времени течения т0 образуют квазиравновесный термостат с температурой потока Т . В подсистеме же колебательных уровней энергии по истечению времени тт устанавливается квазиравновесное распределение с колебательной температурой ТV. Обмен энергией между колебательной модой и квази-равновесными степенями свободы молекул газа описывается релаксационным уравнением Ландау - Теллера с характерным временем т^ = т .
Поскольку поступательные и вращательные степени свободы молекул находятся в состоянии равновесия, то в качестве значений теплоемкостей с1г и сго1, связанных с поступательным и вращательным движением молекул, могут быть выбраны их равновесные значения: с1г = 3Я / 2 , сго1 = Я. В результате, используя
соотношения Эйкена [2, 3]: Х( 0 = 5пс / 2, Xг 0 = 6п0сго1 / 5, Xv 0 = 6п0сУ1Ъ / 5, получаем, что параметр а2 может быть записан в удобном для дальнейших вычислений виде
а = = 12 Уу1Ъ (ctг + Cгot) = 20 УУЪ
2 Х0 (1 -у У1Ъ) (25 с* +12 сго1) 33(1 -у ^
Нижний предел уУ1Ъ = 0 соответствует случаю невозбуждения колебательной моды молекул. С другой стороны, равнораспределение энергии по степеням свободы молекул не является здесь верхним пределом для параметра уУ1Ъ , поскольку закон равнораспределения энергии неприменим в неравновесной ситуации, опи-
сываемой системой уравнений (1), когда разрыв между статической температурой потока Т и колебательной температурой Т, может быть достаточно велик. В монографии [3] показано, что при Т = 300 K неравновесная теплоемкость cvib и 1,8 R . Используя равнораспределение энергии в состоянии термодинамического квазиравновесия по поступательным и вращательным модам молекул, получаем, что параметр yvlb и 0,42 . С ростом разрыва между температурами Т, и Т значение Yvlb увеличивается, приближаясь в пределе к единице, когда энергия
колебательной моды молекул существенно превышает температуру квазиравно-весного термостата, определяемого поступательными и вращательными степенями свободы молекул. В расчетах максимальное значение параметра y vlb было выбрано равным: y vlb = 0,4 с тем, чтобы остаться в рамках используемой модели, избежав возбуждения высоких колебательных уровней энергии.
В принятых безразмерных переменных стационарный несущий поток задается соотношениями
US (x2) — th X2 , Т = Тv, S = pS = 1, pS -1j (Y M ). (2)
В момент времени t = 0 на основной поток накладывалось двумерное возмущение с длиной волны X и волновым вектором k = (Р, 0), где р = 2п / X. Исходная краевая задача ставится в бесконечной полосе, центр которой совпадает с началом координат: х1 е [-х1 0; х1 0 ], х2 е (-да; да). Ширина полосы по координате
х1 0 выбиралась равной длине волны возмущения р = 2п / X, при этом х10 = п / р .
В расчетах асимптотические условия при х2 ^ ±да переносились на х2 = ±х2 0,
где ордината х2 0 определялась из условия достижимой компьютерной точности
| х2 0 -1| < 10-12. В итоге было принято значение х2 0 = 20 . На границах расчетной
области во все моменты времени при х1 = ± х10 ставились периодические условия,
а при х2 = ±х2 0 - условия невозмущенного потока (2). Начальные условия для
поля скорости и термодинамических величин, включающие возмущения, определялись в виде
u1(0, х1, х2) - thх2 + u1(0, х1, х2), u2(0, х1, х2) - u'2(0, х1, х2), р(0, х1, х2) = 1 + р'(0, х1, х2), Т(0, х1, х2) - 1, Т, (0, х1, х2) = £ Т(0, х1, х2), где параметр £ задает амплитуду возмущения колебательной температуры. В работах [5] показано, что при течениях в соплах, недорасширенных струях или для умеренной лазерной накачки колебательных мод, когда диссоциацией можно пренебречь, значения параметра £ для двухатомных газов лежат в интервале
£ = 1 - 5 . В расчетах максимальное значение параметра £ принималось равным
£ = 3.
В качестве вводимых в основной поток начальных возмущений компонент вектора скорости u{, u2 и плотности р' использовались собственные линейные невязкие колебания с наибольшими инкрементами нарастания, которые рассчи-
тывались в [6, 7] на основе линеаризованной системы уравнений двухтемпературной газовой динамики.
В табл. 1 приведены использованные в настоящей работе собственные значения фазовых скоростей с, соответствующих волновых чисел р и инкрементов
роста Рс наиболее неустойчивых невязких мод для чисел Маха несущего потока M = 0 -1, времени колебательной релаксации т = 1 и значений параметра Y vlb = 0,0,4. На рис. 1 представлены примеры изолиний поля завихренности в начальный момент времени t = 0 :
( du i ( du i
Ip v, +—r I cos p x +1 P vr- —- I sin p x ,
\ dx2 ) \ dx2 ) _
где u1 = ur + i щ , u2 = vr + i v, - комплексные компоненты возмущений поля скорости.
Т аблица 1
Спектральные характеристики и инкременты роста наиболее неустойчивых невязких мод при т = 1 и ут1Ь = 0,0,4
M Р Ci Рс,-
0 0,4 0 0,4 0 0,4
0 0,4446 0,4446 0,4266 0,4266 0,1897 0,1897
0,2 0,4260 0,4377 0,4255 0,4115 0,1813 0,1801
0,5 0,3970 0,3890 0,3556 0,3449 0,1413 0,1341
0,8 0,2790 0,2895 0,2790 0,2142 0,0778 0,0620
1 0 0 0 0 0 0
Рис. 1. Картина изолиний поля завихренности <в в момент времени Г = 0 для М = °,5 и т = 1 (а - уУ1Ь = 0 ; б - Уу1ъ = 0,4 )
Для расчета эволюции возмущений уравнения (1) аппроксимировались весовой конечно-разностной схемой Ковени - Яненко [8] с расщеплением по физическим процессам и пространственным переменным. Ее исследование применительно к системе (1) было проведено в [5]. Вычисления велись на равномерной сетке с шагом к по обеим пространственным переменным. По периодической координате х1 выполнялась циклическая прогонка. Шаги по времени и пространству выбирались равными к = 0,025 и Д Г = 0,01. Режимы течения определялись следующими значениями параметров: Яе = 100; М = 0,2, 0,5; Рг = 0,75 ; у = 1,4 ; а1 = 0 - 2; £ = 1 - 3, уУ1Ь = 0,0,4; т = 1 -3.
ю (0) = -sech2 х2 -
Эволюция возмущений в колебательно неравновесном двухатомном газе
Эволюция вихревого возмущения прослеживалась по поведению поля изолиний полной завихренности потока ю . Сеточная функция ю рассчитывалась с помощью аппроксимации со вторым порядком трансверсальной компоненты вектора вихря центральными разностями от компонент поля скорости. Расчеты позволили детально воспроизвести известную картину нелинейной динамики крупной вихревой структуры «са^-еуе» в процессе возникновения и развития неустойчивости Кельвина - Гельмгольца [9]. Структура достаточно быстро при Г = 3 достигает своего максимального размера, после чего возмущение начинает затухать. К моменту времени Г = 6, который можно условно принять за время «жизни», структура минимизируется. При этом ее размер несколько превышает начальный. Примеры изолиний завихренности в вихре в момент достижения максимального размера, для различных значений степени колебательной неравновесности £ приведены на рис. 2.
а б
Рис. 2. Изолинии поля завихренности ю при М = 0,5, а1 = 0 , уу1ъ = 0,4 и т = 1 в момент времени Г = 3 (а - £ = 1; б - £ = 3)
Сравнение картины изолиний на рис. 2 показывает, что при одинаковой начальной амплитуде возмущения поля скорости (завихренности) большее начальное возмущение колебательной моды (рис. 2, б) уменьшает градиент ю в ядре сформировавшейся вихревой структуры. В данном случае это объясняется более высокой температурой газа в ядре и соответственно большим значением сдвиговой вязкости п(Т), что приводит к большей диффузии завихренности.
В процессе приближения к термическому равновесию перетекание энергии из возбужденной колебательной моды в поступательные степени свободы молекул газа приводит к повышению статической температуры потока. Этот результат релаксационного процесса отражен на графиках рис. 3, где приведены профили колебательной и статической температур в поперечном сечении вихря. Сравниваются распределения температур для трех амплитуд возмущений £ колебательной моды молекул в начальный момент времени Г = 0 (кривые 1) и при Г = 3 (кривые 2), когда структура достигает максимального размера. Видно, что при большей амплитуде возмущения £ колебательной моды процесс релаксации идет более интенсивно. Статическая температура в ядре вихря, условно выделяемом замкну-
тыми изолиниями завихренности (см. рис. 2), оказывается выше, чем на периферии области течения. Это связано с тем, что газ, попавший в ядро вихря, остается в нем в течение всей «жизни» вихря, в то время как вне ядра газ выносится конвекцией из расчетной области и релаксация продолжается за ее пределами.
Т
1,35
---5=1,5 ^ - - 5=2 \ —5=3 А
а
1 | \ \ А \ \ \ N \ \
О
Х2
Рис. 3. Профили статической Т(х2) (а) и колебательной Ту(х2) (б) температуры в сечении х1 =0 для Яе = 100 , М = 0,5 , а1 = 0 , ууЛ = 0,4 и т = 1 при Г = 0 (кр. 1) и Г = 3 (кр. 2)
В свою очередь, нагрев приводит к возрастанию вязкости газа и уменьшению плотности в ядре вихря. Следствием этого должно быть дополнительное затухание пульсаций скорости потока. Этот эффект был зафиксирован в работе [10], где моделировалась дорожка Кармана за цилиндром в колебательно возбужденном газе. В наших расчетах выделить в чистом виде этот канал подавления возмущений не удалось. Возможное объяснение этому заключается в том, что в [10] характерное время релаксации, амплитуда возмущений и время похождения вихрем расчетной области существенно превышали соответствующие величины, использованные в данной работе.
Процесс релаксации, свободный от влияния конвекции, можно наблюдать в центре вихря. Соответствующая эволюция во времени статической и колебательной температур к термодинамическому равновесию показана на рис. 4.
На рис. 4, а приведены релаксационные кривые для трех значений начального возбуждения. Рис. 4, б представляет процесс релаксации для трех значений характерного времени релаксации т. Как следует из уравнения Ландау - Теллера, понижение колебательной температуры в центре вихря определяется не только релаксацией, но и молекулярной теплопроводностью. Вместе с тем кривые изменения Ту, особенно на рис. 4, б, демонстрируют преимущественно экспоненциальное убывание. Это обстоятельство вполне согласуется с оценками динамики обоих процессов, которые можно получить из уравнения для Ту. Собственно релаксация приводит к оценке Ту (Г) ~ Ту (0)ехр(-Г / т). В то же время процесс теплопроводности оценивается как Ту (Г) ~ (ЯеРг/ Г)ехр(-сош1 • ЯеРг/ Г). Отсюда видно, что при рассматриваемых параметрах режима диффузионный процесс существенно медленнее релаксационного.
ТТ
ТТу
Рис. 4. Временные зависимости статической Т(Г) и колебательной Ту(Г) температур в центре вихревой структуры для Яе = 100 , М = 0,5, а1 = 0 , уу1Ь =0,4: а - т = 1 (1 - 5 =1,5;
2 -5 = 2; 3 -5 = 3); б - 5 = 3 (1 - т = 1; 2 - т = 2; 3 - т = 3)
С точки зрения возможных приложений интерес представляет влияние колебательной релаксации на диссипацию кинетической энергии возмущений. Рассматривалась эволюция во времени кинетической энергии возмущения и абсолютной величины рейнольдсовых напряжений:
1 х1,0 х2,0 х1,0 х2,0 Е(Г)= — | а х1 | а х2 р(12 + м'22), ст12(Г)= | ёх1 | ёх2 |рм'1м'2|.
Уравнение производства энергии возмущений В (Г) имеет вид [11]
ёЕ 1 , ч
В(Г) = ~Т7 = -Л + —2 — Т;-( -3 +а1 -4), аГ Яе
(3)
,/1 = - | ёх1 | ёх2 р м'1м'2
а х2
■4,0 ■*■2,0 / з ,
Г Л Г Л I д и 1 д м 2
./2= I ёх! I ёх2 р|—— + -
-1 -1 I д х, д х2
—х, Л — х^ ^1 2
-3 = | ёх1 | ёх2 п(Т)
*1,0 Л2,0
д х,
+ 1
I д х2
3 ^ д х1 д х2
1 / 2 х1,0 х2,0
д м '2 д х2
С,;, ( д м '1 д м '2
-4 = I ёх1 I ёх2П(Т) | — + -—
-1 -1 ^ д х1 д х2
■4,0 2,0
Здесь величины м1, м2, - пульсации компонент скорости и давления,
х10 = п / р , х2 0 = 20 . Волновые числа р для различных значений режимных параметров уу, М и т брались из табл. 1.
х
Интегралы -1 и —2 описывают соответственно обмен энергией между возмущением и основным потоком и работу при пульсационном расширении (сжатии) газа. Интегралы -3 и —4, определяющие диссипацию кинетической энергии возмущений, в данном случае включают зависящий от температуры коэффициент сдвиговой вязкости. Последнее естественно сохраняет их положительную определенность. Таким образом, в уравнении (3) вклад релаксации колебательной моды в диссипацию кинетической энергии возмущений осуществляется не только через плотность и давление газа, связанные уравнением состояния с температурой, но и непосредственно через температурную зависимость п(Т) = Т 4/5.
Пульсационные характеристики течения определялись следующим образом:
Т' (Г, х1, х2 ) = Т (Г, х1, х2 ) — Т 5 (Г, х1, х2 ),
где компоненты вектор-функция Т представляют собой мгновенные значения характеристик возмущенного течения, а Т5 - соответствующие характеристики
несущего потока. Поскольку в данном случае невозмущенное течение (2) не является точным стационарным решением системы (1), то его мгновенные характеристики рассчитывались параллельно с расчетом возмущенного потока. Интегралы в (3) вычислялись по квадратурным формулам трапеций с шагом к = 0,025 на сетке, использованной в расчетах.
Примеры временных зависимостей кинетической энергии возмущений Е (Г, 5) и производства пульсационной энергии В(Г, 5) для некоторых режимов представлены соответственно на рис. 5, 6. Видно, что расслоение кривых Е (Г) по параметру неравновесности колебательной энергии 5 достаточно заметное, а эволюция структуры для всех значений 5 носит универсальный характер — рост и достижение максимального значения энергии структуры при Г и 3 , затем ее спад и стабилизация при Г > 5 к некоторому значению, несколько превышающему начальное.
Из графиков рис. 5 следует, что чем больше значение параметра 5 , тем меньше кинетическая энергия структуры на всем временном интервале. Возрастание глубины возбуждения колебательной моды 5 , как и увеличение объемной вязкости Пь
(параметра р
Iдм. дм^ др 1 д
—-+м .—- = —-+-
дГ 1 д х.
К- - —1 у д х Яе д х1
п(т ) дм. д х1.
Яе К 3 у д х.
д м7-
п(Т) Т
д х
),
].
приводит к возрастанию диссипации кинетической энергии возмущений. Временные зависимости модуля рейнольдсовых напряжений | ст12 | (Г, 5, а1) практически повторяют приведенные графики Е(Г, 5, а1). Вместе с тем расчеты показали, что варьирование времени колебательной релаксации т в диапазоне, принятом в вычислениях, слабо влияет на поведение временных зависимостей Е(Г) и | ст12 | (Г).
Сопоставление графиков на рис. 5, 6 показывает, что на временном интервале, где происходит рост кинетической энергии возмущений, ее производство сначала положительно и возрастает, достигая некоторого максимума, а затем начинает убывать и далее становится отрицательным, что соответствует убыванию энергии структуры. Точки перехода кривых В(Г) через нуль соответствуют максимумам на кривых Е (Г), а точки максимума и минимума графиков В (Г) отвечают точкам перегиба на графиках Е(Г).
Е(Г) Е(Г)
Рис. 5. Временные зависимости энергии возмущений Е(Г) для Яе = 100 , М = 0,5 ,
у У1Ь =0,4 и т = 2 ( а - а1 = 0; б - а1 = 2 ; 1 - 5 = 1; 2 - 5 = 2; 3 - 5 = 3)
ДО до
Рис. 6. Временные зависимости производства пульсационной энергии В(Г) для Яе = 100 ,
М = 0,5 , у У1Ь =0,4 и т = 2 ( а - а1 = 0; б - а1 = 2 ; 1 - 5 = 1; 2 - 5 = 2; 3 - 5 = 3)
Для количественного сравнения вклада колебательной релаксации в диссипацию кинетической энергии возмущений находились относительные отклонения
< Е(Г, 5) > — < Е(Г, 1) >
г(5) =
•100%, 5 = 1,5 ^ 3,
< Е (Г, 1) >
где угловые скобки <... > обозначают осреднение по условному времени «жизни» структуры © = 6. Результаты расчетов относительных отклонений е Е (5) для некоторых наборов параметров представлены в табл. 2.
Из табл. 2 следует, что рост параметра £ при фиксированных значениях числа Маха М и параметрах а1, т приводит к большей диссипации кинетической энергии возмущения. В частности, в отсутствие объемной вязкости, когда а1 = 0, относительное уменьшение средней по времени энергии < Е(/, £) > для £ = 3, М = 0,5 достигает еЕ = 9,19% при т = 1 и еЕ = 12,56% при т = 3. Полученное
значение хорошо согласуется с результатами работы [5] для диссипации энергии возмущений < Е(/, £) > за счет только релаксационного процесса колебательной
моды молекул при отсутствии в течение других диссипативных процессов, хотя в [5] рассматривалась модельная задача, в которой воспроизводилось только затухание структуры, что в данном случае соответствует убывающим ветвям кривых на рис. 5. В расчетах для максимальных значений объемной вязкости при а1 = 2, числе Маха М = 0,5, параметре £ = 3 среднее относительное подавление кинетической энергии возмущений составило е Е = 10,07 % для времени колебательной релаксации т = 1 и еЕ = 14,62 % для т = 3 .
Т аблица 2
Относительные отклонения %, для Ке = 100, = 0,4
М т =1
а!=0 а!=1 а!=2
£ = 1,5 £ =2 £ = 3 £ =1,5 £ = 2 £ =3 £ = 1,5 £ = 2 £ = 3
0,2 2,17 3,62 5,15 2,32 4,03 6,06 2,54 4,46 6,68
0,5 3,02 6,04 9,19 3,29 6,31 9,69 4,10 7,18 10,07
М т =2
а!=0 а!=1 а!=2
£ = 1,5 £ =2 £ = 3 £ =1,5 £ = 2 £ =3 £ = 1,5 £ = 2 £ = 3
0,2 3,07 5,12 7,28 3,27 5,69 7,96 3,59 6,31 10,67
0,5 4,24 8,52 11,99 4,45 8,89 12,58 3,10 7,18 12,07
М т =3
а!=0 а!=1 а!=2
£ = 1,5 £ =2 £ = 3 £ =1,5 £ = 2 £ =3 £ = 1,5 £ = 2 £ = 3
0,2 5,32 7,87 9,61 5,66 8,86 11,26 6,23 9,92 12,62
0,5 7,35 9,75 12.56 7,71 10,73 12,31 8,12 11,18 14,62
Можно отметить отраженную в табл. 2 однонаправленность воздействия сжимаемости, объемной вязкости и возбуждения колебательной степени свободы молекул, возрастание которых вызывает увеличение диссипации энергии возмущений. Вместе с тем видно, что допустимые возбуждения колебательной мод молекул приводит к существенно большему диссипативному эффекту по сравнению с возможным вкладом объемной вязкости. Относительная доля последнего не превышает 15 %. Это позволяет рассматривать принудительную накачку колебательных степеней свободы молекул в качестве реального способа управления течениями молекулярных газов.
Заключение
Выполнено численное моделирование нелинейного развития дозвуковых вихревых возмущений в эволюционирующем во времени сдвиговом слое колебательно возбужденного газа. Диапазон параметров течения, в частности отношения коэффициентов объемной и сдвиговой вязкостей ai и параметра неравновесности колебательной моды £, соответствовал реально достижимым значениям для двухатомных газов. Полученные результаты позволяют сделать следующие выводы.
Релаксация неравновесной колебательной моды при уровнях возбуждения, которые можно получить в двухатомных газах при течениях в соплах, недорасши-ренных струях или умеренной лазерной накачке, сопровождается заметным подавлением вихревых возмущений. Вызванное ей относительное увеличение диссипации кинетической энергии крупной вихревой структуры, осредненное по времени ее «жизни», в отсутствие объемной вязкости достигает примерно 13 %. Этот фактор может существенно повлиять на тягу сопел реактивных двигателей, сопротивление тракта газодинамических лазеров и т.п.
Допустимые возбуждения колебательной моды в пределах 1 <£< 3 приводят к существенно большему диссипативному эффекту по сравнению с возможным вкладом объемной вязкости при а1 < 2.
В заключение авторы выражают благодарность профессору Ю.Н. Григорьеву за постоянное внимание к работе и обсуждение ее результатов.
ЛИТЕРАТУРА
1. Григорьев Ю.Н., Ершов И.В., Зырянов К.В. Численное моделирование волн Кельвина-Гельмгольца в слабо неравновесном молекулярном газе // Вычисл. технологии. 2008. Т. 13. № 5. С. 25-40.
2. Жданов В.М. Процессы переноса и релаксации в молекулярных газах / В.М. Жданов, М.Е. Алиевский. М.: Наука, 1989. 335с.
3. Нагнибеда Е.А. Кинетическая теория процессов переноса и релаксации в потоках неравновесных реагирующих газов / Е.А. Нагнибеда, Е.В. Кустова. СПб.: Изд-во С.-Пе-тербурского ун-та, 2003. 272 с.
4. Осипов А.И., Уваров А.В. Кинетические и газокинетические процессы в неравновесной молекулярной физике // УФН. 1992. Т. 162. № 11. С. 1-42.
5. Григорьев Ю.Н., Ершов И.В., Ершова Е.Е. Влияние колебательной релаксации на пуль-сационную активность в течениях возбужденного двухатомного газа // ПМТФ. 2004. Т. 45. № 3. С. 15-23.
6. Григорьев Ю.Н., Ершов И.В. Устойчивость течений колебательно возбужденных газов. Энергетический подход // Вестник Нижегородского университета им. Н.И. Лобачевского. 2011. № 4 (3). С. 735-737.
7. Григорьев Ю.Н., Ершов И.В. Линейная устойчивость невязкого сдвигового течения колебательно возбужденного двухатомного газа // ПММ. 2011. Т. 45. Вып. 4. С. 581-593.
8. Ковеня В.М. Метод расщепления в задачах газовой динамики / В.М. Ковеня, Н.Н. Яненко. Новосибирск: Наука, 1981. 320 с.
9. Patnaik P.C., Sherman F.S., Coreos G.M. A numerical simulation of Kelvin-Helmholtz waves of finite amplitude // J. Fluid Mech. 1976. V. 73. Part 2. P. 215-239.
10. Винниченко Н.А., Никитин Н.В., Уваров А.В. Вихревая дорожка Кармана в колебательно-неравновесном газе // Изв. РАН. МЖГ. 2005. № 5. С. 107-114.
11. Григорьев Ю.Н., Ершов И.В. Подавление вихревых возмущений релаксационным процессом в течениях возбужденного молекулярного газа // ПМТФ. 2003. Т. 44. № 4. С. 22-34.
Статья поступила 06.10.2011 г.
Ershov I.V., Zyryanov K.I. DISSIPATION OF KELVIN-HELMHOLTZ WAVES IN VIBRATIONAL NON-EQUILIBRIUM DIATOMIC GAS. On the basis of the equations of two-temperature aerodynamics, the influence of vibrational relaxation on suppression of the Kelvin-Helmholtz instability in a developing in time shift layer of vibrational non-equilibrium diatomic gas is numerically investigated.
Keywords: Kelvin-Helmholtz instability, vibrational relaxation, kinetic energy of the disturbances, dissipation.
ERSHOV Igor Valer’evich (Novosibirsk State University of Civil Engineering)
E-mail: i_ershov@ngs.ru
ZYRYANOV Kirill Igorevich (Novosibirsk State University of Civil Engineering)
E-mail: k-zyryanov@ngs.ru
