Главная » Просмотр файлов » Диссертация

Диссертация (1150751), страница 5

Файл №1150751 Диссертация (Скорость физико-химической релаксации в вязких неравновесных течениях газов) 5 страницаДиссертация (1150751) страница 52019-06-29СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 5)

. . , L + Lm ,(1.61)26так как для гармонического осциллятора в силу равноудаленности энергий колебательных уровней справедливо соотношениеεci + εck = εci±m + εck∓m .(1.62)В этом случае VV 1 обмен является строго резонансным, и в быстрыхпроцессах колебательные степени свободы не обмениваются энергией сдругими модами.1.3.3. Система уравнений для макропараметровОператор быстрых процессов в многотемпературном приближениипредставляет собой сумму операторов упругих столкновений, столкновений с обменами вращательной энергией и VV 1 -обменов:rapV V1elrotJcij= Jcij+ Jcij+ Jcij,(1.63)а оператор медленных процессов содержит операторы столкновений сVV 2 -обменами колебательной энергией, VT-переходами, бимолекулярными химическими реакциями и реакциями диссоциации–рекомбинации:2⇄22⇄3V V2slVTJcij= Jcij+ Jcij+ Jcij+ Jcij.(1.64)На основании аддитивных инвариантов быстрых процессов (1.57), (1.58),(1.59) имеем следующий набор макропараметров, являющихся функциями координат r и времени t : числовая плотность частиц сорта c — nc ,макроскопическая скорость газа — v , температура газа — T , среднеечисло колебательных квантов в расчете на единицу массы в молекулахсорта c — Wc , определяемое как∑ ∫ρc Wc =i fcij duc .(1.65)ijФункция распределения нормируется относительно данных макропараметров, условия нормировки (1.36)–(1.39) принимают вид∑∫∑ ∫ (0)fcij duc =fcij duc = nc , c = 1, .

. . , L,(1.66)ijij∑cij∫mcuc fcij duc =∑cij∫mc(0)uc fcij duc = ρv,(1.67)27)∑ ∫ ( mc c2ccc+ εi + εj + εc fcij duc =2cij)∑ ∫ ( mc c23(0)c=+ εci + εcj + εc fcij duc = nkT + Erot + Evibr + Ef , (1.68)22cij∑ ∫∑ ∫ (0)i fcij duc =i fcij duc = ρc Wc .ij(1.69)ijЗамкнутую систему уравнений переноса и релаксации получаем,умножая уравнение (1.30) на аддитивные инварианты ψcij , интегрируяпо скорости, и суммируя по колебательным и вращательным уровням (вслучае инвариантов количества движения и полной энергии также производим суммирование по химическим сортам).

В результате имеем:dnc+ nc ∇ · v + ∇ · (nc Vc ) = Rcreact , c = 1, . . . , L,dtdvρ= ∇ · P,dtdU= −∇ · q + P : ∇v,ρdt(1.70)(1.71)(1.72)dWc+ ∇ · qw,c = Rcw − Wc mc Rcreact + Wc ∇ · (ρc Vc ), c = 1, . . . , Lm . (1.73)dtУравнение неразрывностиρcdρ+ ρ∇ · v = 0(1.74)dtявляется следствием уравнений (1.70) и получается, если умножить (1.70)на mc и просуммировать по c .В уравнении (1.73) qw,c — поток числа колебательных квантов молекул сорта c , определяемый как∑ ∫qw,c =i cc fcij duc ,(1.75)ijрелаксационные члены Rcreact (описывающие изменение числа частиц сорта c за счет химических реакций), Rcw (описывающие изменение числаколебательных квантов молекул сорта c вследствие всех медленных процессов) определяются как∑ ∫∑∫slreactRcw =i Jcijduc ,Rcreact =Jcijduc .(1.76)ijij28Произведем разбиение оператор медленных процессов согласно (1.24):∑N V V2slJcij=V V2 ,rνr,ci Jcijr=1+NV T∑∑N2⇄2V T,rνr,ci Jcij+r=1∑N2⇄32⇄2,rνr,ci Jcij+r=12⇄3,rνr,ci Jcij, (1.77)r=1а также введем для удобства следующие множества различных типов медленных процессов:V = {V V2 , V T },R = {2 ⇄ 2, 2 ⇄ 3},VR = {V V2 , V T, 2 ⇄ 2, 2 ⇄ 3}.(1.78)С учетом определения скорости медленных процессов (1.26), релаксационные члены (1.76) можно переписать в следующем виде [128]:Rcreact = NaNγ∑∑∑iRcwνr,ci ξ˙γ,r ,γ∈R r=1NγNγ∑ ∑∑ ∑ ∑νr,ci ξ˙γ,r = Na∆iξ˙γ,r ,= Naiiγ∈VR r=1(1.79)(1.80)γ∈VR r=1где ∆i = i′ − i , i′ — колебательный уровень после столкновения.Система уравнений (1.70)-(1.72) представляет собой систему уравнений гидродинамики, которая включает уравнения изменения числовыхплотностей в силу химических реакций, уравнение движения, уравнение переноса энергии.

Эта система дополняется уравнениями релаксации(1.73), которые включают в себя члены, связанные с диффузией колебательных квантов и диффузией частиц разных химических сортов. Такжеуравнения (1.73) включают в себя релаксационные члены Rcw , Rcreact , которые описывают изменение числа колебательных квантов в системе засчет переходов колебательной энергии и химических реакций.1.3.4. Нулевое приближениеВ нулевом приближении метода Энскога–Чепмена функция распределения представляет собой распределение Тринора [169], обобщенное наслучай произвольной смеси [85]:)(( m )3/2 nεcijmc c2ccc c(0)fcij =−− λc i ,(1.81)s exp −2πkTZcint ij2kTkT29Здесь scij — статистический вес состояния молекулы, находящейся на i -мколебательном и j -м вращательном уровнях, k — постоянная Больцмана, Zcint — неравновесная статистическая сумма по внутренним степенямсвободы, и()εc1 11λc =−,(1.82)k T1c TT1c — температура первого колебательного уровня молекул сорта c , вводимая по аналогии с [169] на основании отношения заселенностей первогои нулевого колебательных уровней:()sc1εc1nc1= c exp − c ,(1.83)nc0s0kT1где sci — статистический вес колебательного уровня i .

Отметим, что длядвухатомных молекул sci = 1 , поэтому в дальнейшем sci будем опускать.Статистическая сумма по внутренним степеням свободы определяется как(()c )c∑∑εεjZcint = Zcrot Zcvibr , Zcrot =scj exp −, Zcvibr =exp − i − λc i .kTkTji(1.84)На основании вида функции распределения (1.81) заселенности колебательных уровней вычисляются следующим образом:( c)iεc1ncεi − iεc1c−nci (T, T1 ) = vibrexp −.(1.85)Zc (T, T1c )kTkT1cПодставим функцию распределения нулевого приближения в интегральные операторы медленных процессов.

С учетов микроскопическихзаконов сохранения получим следующее представление оператора процесса r :∑ ∫ (0) (0)r(0)rJcij = Γfcij fdkl geσr dud .(1.86)lj ′ l′В случае VT-переходов необходимо опустить суммирование по l′ , а в случае реакций диссоциации–рекомбинации — по j ′ l′ . Величины Γr определяются в виде( )ArrΓ = 1 − exp,(1.87)kT30где Ar — аналоги сродства химических реакций, обобщенные на случайколебательной неравновесности [132]. Величины Ar определяются следующим образом:(())TTcdAV V2 ,r = ∆iε1 1 − c + ∆kε1 1 − d ,(1.88)T1T1()TcAV T,r = ∆iε1 1 − c ,(1.89)T1A2⇄2,r3mc mdZcint Zdintnc nd= kT ln+ kT ln int int − kT ln+2mc′ md′nc′ nd′Zc′ Zd′′′+ εc′ + εd′ − εc − εd + εc0 + εd0 − εc0 − εd0 +()()()()TTTT′′+ iε1c− 1 − i′ εc1− k ′ εd1− 1 + kεd1, (1.90)′ − 1′ − 1cdT1cT1T1T1d3mc3− ln(2πkT ) + 3kT ln h + kT ln Zcint −A2⇄3,r = kT ln2mc′ ,a mf ′ ,a 2)(ncTcc− kT ln+ εc′ ,a + εf ′ ,a − εc − ε0 + iε1− 1 , (1.91)nc′ ,a nf ′ ,aT1cгде c′ , a и f ′ , a обозначают атомарные продукты реакций диссоциации,а εc0 = 21 hνc − 14 hνc αc — энергия нулевого колебательного уровня, рассчитанная по формуле (1.41).Подставляя функцию распределения (1.81) в выражения для теплового потока, скоростей диффузии, потока числа колебательных квантови тензора напряжений, имеем:q(0) = 0,V(0)c = 0,q(0)w,c = 0,P(0) = −pU,(1.92)где p = nkT , U — единичный тензор.Таким образом, уравнения переноса в нулевом (невязком) приближении метода Энскога–Чепмена принимают видdnc+ nc ∇ · v = Rcreact(0) , c = 1, .

. . , L,dtdv= −∇p,dtdUρ= −p∇ · v,dtρ(1.93)(1.94)(1.95)31dWc= Rcw(0) − Wc mc Rcreact(0) , c = 1, . . . , Lm .(1.96)dtСистема содержит уравнение движения Эйлера (1.94) и уравнениеэнергии (1.95), дополненные уравнениями многотемпературной химической кинетики (1.93) и релаксационными уравнениями (1.96) для числаколебательных квантов молекулярных сортов.Скорость процессов в нулевом приближении линейно зависит от величины Γr :(r))νr,ciLc (L ∏∏nciξ˙r(0) = Γr kr(0)(1.97)Nac=1 i=1ρc(0)и от коэффициента скорости соответствующей реакции kr , получаемойподстановкой функции распределения (1.81) в определение (1.27). Очевидно, что для функции распределения (1.81) выражение (1.97) отличается от закона действующих масс, так как скорость реакции зависит нетолько от температуры газа T , но и от температуры первого колебательного уровня T1c (через заселенности колебательных уровней).

Когда течение достигает колебательного равновесия (т.е. T = T1c ), распределениезаселенностей колебательных уровней становится равновесным Больцмановским, и уравнение (1.97) сводится к закону действующих масс.Подставляя функцию распределения нулевого приближения в выражения для удельной полной энергии U и среднего числа колебательныхквантов Wc , имеемmm∑∑∑3ρc Erot,c (T ) +ρc Evibr,c (T, T1c ) +ρc Ef,c = ρU (nc , T, T1c ),ρU = nkT +2c=1c=1c=1(1.98))( ccc∑∑iε11ε−iε1Evibr,c (T, T1c ) =εci nci =− 1c ,εci exp − icvibrρc imc Zc (T, T1 ) ikTkT1)(1.99)( ccc∑∑11iε1εi − iε1Wc (T, T1c ) =inci =−.iexp−cρc imc Zcvibr (T, T1 ) ikTkT1c(1.100)Отсюда видно, что макроскопические уравнения переноса (1.70)-(1.73) могут быть переписаны в переменных nc , v , T и T1c вместо nc , v , U , Wc ,LLL32на основании следующих выражений:LmL∑∑∂U∂U∂UdU =dnc +dT +dT1c ,c∂nc∂T∂T1c=1c=1(1.101)∂Wc∂WcdT +dT1c .(1.102)c∂T∂T1Вычисляя производные, получаем, что⟨ c⟩εj rot ρc ⟨ c ⟩∂U3kT3nkT mc=−+− 2 εj rot +∂nc2ρ2ρ2ρρ⟨εc ⟩ρcεc ρc+ i vibr − 2 ⟨εci ⟩vibr + − 2 εc , (1.103)ρρρρdWc =mm∑ρc ∂Erot,c ∑ρc ∂Evibr,c13∂U=nk ++,∂Tρ2ρ∂Tρ∂Tc=1c=1LL(1.104)∂Uρc ∂Evibr,c=,(1.105)c∂T1ρ ∂T1c⟨ c ⟩2 )∂Erot,c1 (⟨( c )2 ⟩(1.106)=εj− εj rot ,rot∂Tmc kT 2)∂Evibr,c1 (⟨ c 2 ⟩c 2cccc=(εi )− ⟨εi ⟩vibr + ε1 ⟨i⟩vibr ⟨εi ⟩vibr − ε1 ⟨iεi ⟩vibr ,vibr∂Tmc kT 2(1.107)cε1∂Evibr,ccc=(1.108)2 (⟨iεi ⟩vibr − ⟨i⟩vibr ⟨εi ⟩vibr ) ,cc∂T1mc k (T1 ))⟨ 2⟩∂Wc1 ( c2ccc=⟨iεi ⟩vibr − ε1 i vibr − ⟨i⟩vibr ⟨εi ⟩vibr + ε1 ⟨i⟩vibr , (1.109)∂Tmc kT 2(⟨ ⟩)εc1∂Wc22=i vibr − ⟨i⟩vibr ,(1.110)∂T1cmc k (T1c )2где угловые скобки означают операции осреднения по вращательной иколебательной энергиям, которые определяются следующим образом:)(εcj1 ∑ c,(1.111)⟨ζj ⟩rot = rotsj ζj exp −ZckTj⟨ζi ⟩vibr =1∑Zcvibri)( ccciεε−iε1− 1c .sci ζi exp − ikTkT1(1.112)Таким образом, видно, что наборы макроскопических переменных( nc , v , U , Wc ) и ( nc , v , T , T1c ) эквивалентны.33Релаксационные уравнения (1.96) отличаются от обычно применяемого в вычислительной аэродинамике уравнения релаксации колебательной энергии, т.к.

в случае ангармонического колебательного спектра инвариантом является именно число колебательных квантов. Релаксационныеreact(0)w(0)также учитывают изменение числа колебательныхчлены Rc , Rcквантов в системе за счет химических реакций. При этом стоит отметить, что уравнения (1.96) выведены из строгих определений кинетической теории и не делают предположений относительно вида релаксационw(0)react(0)ных членов Rc , Rc. Таким образом, связь колебательной релаксации и химических реакций вводится естественным образом и не требуетдополнительных предположений, привлекаемых в феноменологическихмоделях [118,143]. С учетом представления (1.97) для скорости реакций внулевом приближении метода Энскога–Чепмена, релаксационные членыпредставляют собой линейные функции коэффициентов скорости реакций.

Характеристики

Список файлов диссертации

Скорость физико-химической релаксации в вязких неравновесных течениях газов
Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6549
Авторов
на СтудИзбе
300
Средний доход
с одного платного файла
Обучение Подробнее