Бул изилдөө Жердин тартылуу күчү шартында көлөм боюнча %6 H₂ камтыган өтө-жарды суутек–аба аралашмасында пайда болгон жалын өзөгүнүн өсүшүн, көтөрүлүшүн, формасынын өзгөрүшүн жана бөлүнүшүн сандык түрдө изилдейт. Өтө-жарды суутек жалын өзөгү — тутандыргандан кийин суутектин артыкчылыктуу жана термиялык диффузиясы, химиялык жылуулук бөлүнүшү, радиациялык муздоо жана көтөрүүчү күч биргелешип аныктаган чакан жана убакытка көз каранды реакциялык түзүлүш. Изилдөөдө суутек химиясынын деталдуу механизми жана төмөн-Mach-сандуу газ динамикасы колдонулган; Соре диффузиясы менен радиациялык жылуулук жоготуусу өз алдынча күйгүзүлүп-өчүрүлүп, төрт физикалык моделдин айкалышы салыштырылган.

Тутандыргандан кийинки алгачкы этап бир өлчөмдүү сфералык симуляциялар менен изилденип, болжол менен алгачкы 200 ms ичинде жалын фронтунун радиусунун квадраты убакыт менен болжолдуу сызыктуу өзгөрөрү, башкача айтканда \(R_f^2 \simeq Kt\) түрүндөгү диффузия башкарган өсүү режими байкалганы аныкталган. Соре диффузиясы натыйжалуу өсүү коэффициенти \(K\)-ны жогорулатса, радиациялык жылуулук жоготуусу аны төмөндөтөт. Бул жүрүм-турум туруктуу ылдамдыктагы сфералык жалындын таралышынан айырмаланат; өзөк чоңойгон сайын фронттун учурдагы радиалдык ылдамдыгы азаят.

Эки өлчөмдүү симуляцияларда өзөк чоңойгон сайын көтөрүүчү күч үстөм боло баштап, ысык продукт аймагы жогору көтөрүлүп, өзүнүн тороидалдык куюнун түзгөн. Башында болжолдуу сфералык болгон түзүлүш кийин жалпак эллипстик калпакчага айланган. Соре диффузиясы чоңураак жана ылдамыраак көтөрүлгөн өзөктөрдү түзүп, үстүңкү калпакчаны тегерегирээк сактаган. Радиация көтөрүлүү траекториясына азыраак таасир эткен, бирок салыштырмалуу каптал жалпактоону күчөткөн.

Бардык физикалык-модель айкалыштарында жалындын бөлүнүү жолу негизинен бирдей болгон. Тороидалдык куюн үстүңкү калпакча менен ылдый ийилген жалын бутактын ортосундагы тар реакция аймагын созуп, жергиликтүү өчүүгө алып келген. Соре диффузиясы бөлүнүүнү кечиктирген; радиация болсо эртелеткен. Эки таасир тең активдүү болгондо бөлүнүү 305 msде, экөө тең өчүк болгондо 306 msде болгон жана 1 ms айырма изилдөөнүн убакыт аныктоо тактыгында калгандыктан маанилүү физикалык жылыш катары чечмеленген эмес.

Изилдөө жаңы физикалык күйүү эксперименти эмес. Натыйжалар бир өлчөмдүү сфералык, эки өлчөмдүү тегиз жана эки өлчөмдүү ок-симметриялуу сандык моделдерге негизделет. Ок-симметриялуу натыйжалар мурда жарыяланган эксперименттик алдыңкы-чекит траекториялары менен салыштырылган; толук үч өлчөмдүү асимметриялуу бөлүнүү жана турбуленттүү жалындын жүрүм-туруму бул изилдөөдө чечилген эмес.

Изилдөө маселеси

Суутек–аба аралашмалары күйүү чегине жакындаганда классикалык, жылмакай жана болжолдуу тегиз premix жалын жүрүм-турумунан олуттуу четтей алат. Химиялык жылуулук бөлүнүшү жайлаган сайын молекулалык ташуу, жылуулук жоготуусу, ийрилик, жалындын созулушу жана табигый конвекциянын салыштырмалуу мааниси өсөт. Молекулалык суутектин жогорку диффузиялуулугунан улам суутек ысык реакция аймагына башка көптөгөн түрлөргө караганда тезирээк ташылышы мүмкүн.

Жердин тартылуу күчү шартында буга көтөрүүчү күч да кошулат. Тутандыргандан кийин пайда болгон ысык жана тыгыздыгы төмөн продукт өзөгү жогору карай жылат. Бул кыймыл жалындын ордун гана өзгөртпөйт; ошол эле учурда жалындын айланасында ири масштабдуу куюн жаратат. Демек, байкалган өйдө багытталган жалын кыймылы бир гана ламинардык күйүү ылдамдыгынын натыйжасы эмес. Химиялык фронттун кыймылы, ысык газ массасынын көтөрүлүшү, радиалдык жана октук конвекция жана үзгүлтүксүз форма өзгөрүшү чогуу роль ойнойт.

Булак изилдеген негизги боштук — Соре термиялык диффузиясы менен радиациялык жылуулук жоготуусунун эркин көтөрүлгөн өтө-жарды суутек жалын өзөгүнө биргелешкен таасири Жердин тартылуу күчү шартында системалуу түрдө сандык бааланбаганы. Мурунку изилдөөлөр бул механизмдерди көбүнчө микрогравитациядагы жалын шарлары, дубал менен чектелген жалындар, күйгүзгүч менен турукташтырылган жалындар же табигый конвекциясыз клеткалык жалындар контекстинде караган.

Соре Диффузиясы Өтө-Жарды Суутек Жалында Эмне Кылат?

Соре диффузиясы температура градиентинин жеңил H жана H₂ түрлөрүндө кошумча масса агымын жаратуусуна алып келет; изилдөөнүн өтө-жарды суутек аралашмасында бул механизм суутектин ысык реакция аймагына ташылышын күчөтөт, өзөктүн өсүшүн ылдамдатат, көбүрөөк ысык продукт пайда болушуна көмөктөшөт жана куюндан келип чыккан созулууга каршы байланышкан жалын түзүлүшүнүн узагыраак сакталуусуна жардам берет.

Соре таасири жана кадимки молекулалык диффузия

Булакта аралашма-орточо түр агымынын оңдолбогон формасы төмөнкүдөй берилген:

\[ \mathbf{j}_k^{(0)} = -\rho D_k \nabla Y_k - D_k^T \nabla \ln T \tag{10} \]

Бул жерде \(\mathbf{j}_k^{(0)}\), \(k\) түрүнүн оңдолбогон диффузиялык масса агымы; \(\rho\) аралашманын тыгыздыгы, \(D_k\) аралашма-орточо молекулалык диффузия коэффициенти, \(Y_k\) масса үлүшү, \(D_k^T\) термиялык-диффузия коэффициенти жана \(T\) температура.

Биринчи мүчө түрдүн концентрациясынын мейкиндиктеги өзгөрүшүнөн келип чыккан кадимки диффузияны, экинчи мүчө температура градиенти жараткан Соре ташылышын көрсөтөт. Соре таасири өчүрүлгөндө бардык \(D_k^T\) маанилери нөлгө коюлган.

Жалпы диффузиялык масса агымы нөл болушу үчүн оңдоо ылдамдыгы колдонулган:

\[ \mathbf{j}_k= \mathbf{j}_k^{(0)} - Y_k\sum_{l=1}^{N_s}\mathbf{j}_l^{(0)}, \qquad \sum_{k=1}^{N_s}\mathbf{j}_k=0. \tag{11} \]

Термиялык диффузия H жана H₂ үчүн гана Schlup–Blanquart кыскартылган модели менен колдонулган. Булактагы коэффициенттер \(\alpha_H=0{,}895\) жана \(\alpha_{H_2}=0{,}910\). Оор түрлөрдүн түздөн-түз термиялык-диффузия таасири эске алынган эмес. Өз ара Dufour жылуулук агымы premix суутек жалындарында Соре салымына салыштырмалуу кичине деп кабыл алынгандыктан моделге кошулган эмес.

Радиациялык жылуулук жоготуусу кантип моделденген?

Радиация оптикалык жактан жука эмиссия ыкмасы менен көрсөтүлгөн жана радиация чыгаруучу жалгыз түр суу буусу болгон. Көлөмдүк радиациялык булак:

\[ \dot q_{\mathrm{rad}} = -4\sigma_{\mathrm{SB}} \left(T^4-T_b^4\right)\kappa_P, \qquad T_b=300\ \mathrm{K} \tag{12} \]

түрүндө. \(\sigma_{\mathrm{SB}}\) Stefan–Boltzmann туруктуусу, \(T_b\) фондук радиация температурасы жана \(\kappa_P\) Planck-орточо сиңирүү коэффициенти.

Planck-орточо коэффициенти булакта Barlow жана кесиптештеринин полиномдук корреляциясы боюнча эсептелген:

\[ \kappa_P = X_{H_2O} \frac{p_0}{p_{\mathrm{ref}}} \sum_{n=0}^{5} a_n \left(\frac{1000}{T}\right)^n, \qquad p_{\mathrm{ref}}=101325\ \mathrm{Pa}. \tag{13} \]

Оптикалык жука ыкма продукттардын өз радиациясын кайра сиңиришин эске албайт. Оптикалык калыңдык көбөйгөндө бул жоромол таза жылуулук жоготуусун жогору баалоого ыкташат. Ошондуктан изилдөө деталдуу спектралдык радиациялык берүүнү эмес, сантиметр масштабындагы бул жалын өзөгүндө радиациялык муздоонун негизги таасирин чечет.

Баштапкы аралашма жана тутандыруу

Сфералык эталон маселеде борбордогу радиусу 0,4 mm болгон аймак заматта 2000 Kге көтөрүлгөн, бирок баштапкы химиялык курамы өзгөртүлгөн эмес. Ошентип, алгачкы этаптагы өзөктүн өсүшү тартылуу күчүнөн келип чыккан форма өзгөрүшүнөн бөлүп изилденген.

Эки өлчөмдүү маселелерде тутандыруу борбору 300 K изотермиялык төмөнкү дубалдан 6,4 mm жогору жайгашкан. Эсептөө аймагынын жарым туурасы же ок-симметриялуу радиусу 51,2 mm, бийиктиги 102,4 mm. Үстүңкү жана тышкы чек аралар эркин чыгуу, ал эми ок симметрия чеги.

Төрт физикалык модель айкалышы

Бул матрица эки физикалык механизмдин өзүнчө таасирлерин жана алар бирге иштегендеги натыйжаны бир эле баштапкы аралашмада салыштырууга мүмкүндүк берет.

Төмөн-Mach математикалык модели

Реактивдүү аралашма термиялык жактан кемчиликсиз идеал газ катары каралган. Төмөн-Mach ыкмасында басым мейкиндикте болжол менен бир калыптагы термодинамикалык бөлүк жана гидродинамикалык пертурбацияга бөлүнөт:

\[ p(\mathbf{x},t)=p_0(t)+\pi(\mathbf{x},t). \tag{1} \]

Бул изилдөөдө \(p_0=101325\) Pa туруктуу кармалган.

Массанын сакталуусу:

\[ \frac{\partial \rho}{\partial t} +\nabla\cdot(\rho\mathbf{u})=0. \tag{2} \]

Түрлөрдүн сакталуусу:

\[ \frac{\partial(\rho Y_k)}{\partial t} +\nabla\cdot(\rho Y_k\mathbf{u}) = -\nabla\cdot\mathbf{j}_k+\dot\omega_k, \qquad k=1,\ldots,N_s. \tag{3} \]

Импульс теңдемеси:

\[ \frac{\partial\mathbf{u}}{\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u} = -\frac{1}{\rho}\nabla\pi + \frac{1}{\rho}\nabla\cdot\boldsymbol{\tau} + \frac{\rho-\rho_0}{\rho}\mathbf{g}. \tag{4} \]

Акыркы мүчө тыгыздык айырмасынан пайда болгон көтөрүү таасирин көрсөтөт. Гидростатикалык басым компоненти басым талаасынан алынып салынгандыктан, эталон тыгыздык \(\rho_0\) ушул мүчөдө көрүнөт.

Сезилүүчү энтальпия теңдемеси жылуулук өткөрүүнү, түрлөрдүн энтальпия ташылышын, химиялык жылуулук бөлүнүшүн, илешкектүүлүк мүчөсүн жана радиацияны бир энергия теңдемесинде бириктирет:

\[ \frac{\partial(\rho h_s)}{\partial t} +\nabla\cdot(\rho h_s\mathbf{u}) = \frac{dp_0}{dt} +\nabla\cdot(\lambda\nabla T) -\nabla\cdot \left( \sum_{k=1}^{N_s}h_{s,k}\mathbf{j}_k \right) -\sum_{k=1}^{N_s}\Delta h_{f,k}^{0}\dot\omega_k +\boldsymbol{\tau}:\nabla\mathbf{u} +\dot q_{\mathrm{rad}}. \tag{5} \]

Суутектин кычкылданышы үчүн Kéromnès жана кесиптештеринин деталдуу химиялык механизми колдонулуп, азот инерттүү түр катары каралган. Температурага көз каранды жылуулук сыйымдуулуктары жана энтальпиялар JANAF термохимиялык маалыматтарына негизделген полиномдор менен көрсөтүлгөн.

Сандык чечкич жана grid

Теңдемелер төмөн-Mach FDS алгоритминин ички NRG ишке ашырылышы менен чечилген. Тыгыздык жана түрлөр химия жана ташуу булактары менен бирге predictor–corrector projection схемасы аркылуу алдыга жылдырылган; абал теңдемеси талап кылынган ылдамдык дивергенциясын берген; гидродинамикалык басым үчүн өзгөрмө коэффициенттүү Poisson теңдемеси multigrid ыкмасы менен чечилген.

Агым ламинардуу деп кабыл алынып, турбуленттүүлүк же sub-grid күйүү модели колдонулган эмес. Убакыт кадамы адаптивдүү болуп, төмөн-Mach чечкичинин сандык туруктуулук критерийлери менен чектелген.

Эки өлчөмдүү эсептөөлөрдө негизги уяча аралыгы 0,4 mm жана чечилген жарым аймак 128 × 256 уячаны камтыйт. Соре күйүк, радиация өчүк абал эң жука жалын фронтун бергендиктен grid жакындашуу тести ушул абалда жүргүзүлгөн. 0,1, 0,2 жана 0,4 mm gridдерде алдыңкы-чекит траекториясы менен 50, 100 жана 150 msдеги \(T=1000\) K изотермалары дээрлик дал келип, кийинки эсептөөлөрдө 0,4 mm grid колдонулган.

\(R_f^2 \simeq Kt\) Режими Эмнени Билдирет?

\(R_f^2 \simeq Kt\) режими изилденген тутандыргандан кийинки алгачкы аралыкта сфералык жалын өзөгүнүн радиус квадраты болжол менен туруктуу ылдыйыш менен өсөрүн жана ошондуктан фронттун ылдамдыгы туруктуу эмес экенин билдирет; \(dR_f/dt=K/(2R_f)\) болгондуктан өзөк чоңойгон сайын радиалдык алдыга жылуу ылдамдыгы болжол менен \(R_f^{-1}\) менен азаят. Булак бул байланышты алгачкы болжол менен 200 ms үчүн натыйжалуу диффузия-башкарылган режим катары чечмелейт жана узак убакыт үчүн универсалдуу асимптотикалык мыйзам катары бербейт.

Натыйжалуу өсүү коэффициенти K

Жалын фронтунун fit кылынган жүрүм-туруму:

\[ R_f^2(t) = R_f^2(t^*) + K(t-t^*) \tag{14} \]

түрүндө. Туундасы:

\[ U_f = \frac{dR_f}{dt} = \frac{K}{2R_f}. \tag{15} \]

Ошондуктан \(K\) молекулалык диффузия коэффициенти же ламинардык күйүү ылдамдыгы менен бирдей эмес. Булакта \(K\) артыкчылыктуу диффузия, жылуулук өткөрүү, чектүү ылдамдыктагы химия, термиялык кеңейүү, Соре ташылышы жана радиациялык жоготуунун жалпы таасирин камтыган натыйжалуу өзөк-өсүү коэффициенти катары аныкталат.

Таблица Соре диффузиясы радиациялуу да, радиациясыз да учурда \(K\)-ны жогорулатарын, ал эми радиация Соре абалынын экөөндө тең \(K\)-ны төмөндөтөрүн көрсөтөт.

Эмне үчүн \(R_f^2\) болжол менен сызыктуу өсөт?

Булактын масштабдоо аргументинде жетишпеген реактантынын сфералык реакция аймагына диффузиялык ташылышы болжол менен:

\[ \dot m_D \sim 4\pi R_f^2 \rho D_{\mathrm{eff}} \frac{\Delta Y}{R_f} \propto R_f \tag{16} \]

түрүндө масштабданат. Сфералык фронттун өсүшү үчүн зарыл керектөө болсо:

\[ \dot m_F \sim 4\pi R_f^2\rho_u\dot R_f. \tag{17} \]

Эки масштабды тең салмактоо \(\dot R_f\propto R_f^{-1}\) натыйжасын жана ошентип \(d(R_f^2)/dt\) маанисинин болжол менен туруктуу болушун түшүндүрөт. Бул химия маанисиз дегенди билдирбейт; химиялык реакция диффузияны камсыз кылган температура жана концентрация градиенттерин түзүүнү улантат.

Ысык продукт көлөмү жана көтөрүүчү күч

Сфералык ысык аймактын көлөмү:

\[ V_h=\frac{4}{3}\pi R_f^3. \tag{18} \]

\(R_f^2\sim Kt\) болгондо көлөм болжол менен:

\[ V_h\sim\frac{4\pi}{3}(Kt)^{3/2} \tag{19} \]

түрүндө масштабданат. Ошондуктан \(K\)-дагы салыштырмалуу орточо өзгөрүү кийин көтөрүлүү кыймылын камсыз кылган ысык газ көлөмүндө чоңураак салыштырмалуу өзгөрүү жаратышы мүмкүн.

Булактын интегралдык көтөрүүчү күч туюнтмасы:

\[ F_b = g\int_{\Omega_h}(\rho_u-\rho)\,dV \approx g(\rho_u-\rho_h)V_h \tag{20} \]

түрүндө. Ошентип, алгачкы сфералык этапта Соре же радиация өзгөрткөн ысык продукт көлөмү кийинки көтөрүлүү жүрүм-турумуна өтөт.

Радиация 200 msде эмнени өзгөрттү?

Булакта 200 msдеги радиалдык профилдер радиация продукт температурасын болжол менен 61 Kге төмөндөтөрүн көрсөтөт. Бул төмөндөө айлана-чөйрөдөн жогору температура өсүшүнүн болжол менен %7 көрсөткүчүнө барабар. Мунун ордуна продукттун тыгыздыгы болжол менен 0,015 kg/m³ көбөйгөн жана бул өсүш айланадагы аралашманын тыгыздыгынын болжол менен %1,89 көрсөткүчүнө барабар болгон.

Төмөн-Mach ыкмасында болжол менен туруктуу термодинамикалык басым шартында тыгыздык:

\[ \rho=\frac{p_0M}{R_uT} \tag{21} \]

менен берилет. Кичине өзгөрүүлөр үчүн:

\[ \frac{\delta\rho}{\rho} \approx \frac{\delta M}{M} - \frac{\delta T}{T}. \tag{22} \]

Бул туюнтма температура төмөндөшүнүн тыгыздыкка таасири курам жана орточо молярдык массадагы өзгөрүү менен жарым-жартылай теңделиши мүмкүн экенин көрсөтөт.

Изилдөөнүн Ыкмасы жана Жыйынтыктары

Талдоо чынжыры

1. %6 көлөмдүк H₂–аба аралашмасы үчүн төрт Соре/радиация айкалышы аныкталды.
2. Бир өлчөмдүү сфералык эсептөөлөр менен тартылуу күчүнөн көз каранды эмес алгачкы өзөк өсүшү изилденди.
3. \(R_f^2(t)\) ийри сызыктары аркылуу натыйжалуу өсүү коэффициенти \(K\) аныкталды.
4. Температура, тыгыздык, ылдамдык жана H₂/O₂/H₂O профилдери аркылуу Соре менен радиациянын физикалык таасирлери ажыратылды.
5. Эки өлчөмдүү тегиз жана ок-симметриялуу эсептөөлөр менен табигый конвекция жана көтөрүлүү чечилди.
6. Алдыңкы-чекит траекториялары мурда жарыяланган эки эксперименттик маалымат топтому менен салыштырылды.
7. Температура жана OH талаалары менен ысык продукт аймагы жана активдүү реакция катмары ажыратылды.
8. Үстүңкү жалын калпакчасы эллипске fit кылынып, форма катышы \(a_r/a_z\) көзөмөлдөндү.
9. Тороидалдык куюн менен реакция аймагынын өз ара аракеттенүүсү жана биринчи байланыштын үзүлүшү боюнча бөлүнүү убактысы аныкталды.

Тартылуу күчү качан маанилүү боло баштайт?

Алгачкы этапта өзөк болжол менен сфералык жана диффузия башкарган. Булак жалындын өсүү убактысы менен көтөрүүчү күч аркылуу көтөрүлүү убактысын салыштырган критерийди колдонгон. Davies–Taylor негизиндеги көтөрүлүү ылдамдыгы:

\[ V_{\mathrm{rise}} = \frac{2}{3}\sqrt{g'R}, \qquad g'=g\left(1-\frac{\rho_b}{\rho_u}\right) \tag{23} \]

түрүндө жазылат.

Бул изилдөөнүн радиуска көз каранды өсүү ылдамдыгы колдонулганда убакыт масштабдарынын катышы:

\[ \Pi_{bK} = \frac{t_{\mathrm{growth}}}{t_{\mathrm{rise}}} = \frac{V_{\mathrm{rise}}}{\dot R_f} = \frac{4}{3} \sqrt{\frac{g'R_f^3}{K^2}} \tag{24} \]

болот. \(\Pi_{bK}\ll1\) диффузия башкарган болжолдуу изотроптук өсүүнү, \(\Pi_{bK}=O(1)\) болсо байкаларлык көтөрүүчү күчтөн келип чыккан кыймыл жана деформациянын башталышын көрсөтөт.

Буга туура келген мүнөздүү радиус:

\[ R_{bK}^{*} = \left( \frac{9K^2}{16g'} \right)^{1/3}. \tag{25} \]

Соре да, радиация да күйүк болгондо \(K=1{,}113\times10^{-3}\) m²/s жана булакта колдонулган өкүлчүлүк тыгыздыктар менен \(R_{bK}^{*}\) болжол менен 4,6 mm табылган. Бул маани так эксперименттик чек эмес, көтөрүүчү күч маанилүү боло баштаган масштаб үчүн тартиптик баа.

Тегиз жана ок-симметриялуу натыйжалардын айырмасы

Эки модель тең өзөктүн убакыт өткөн сайын ылдамдап көтөрүлөрүн жана төрт Соре/радиация айкалышынын бирдей салыштырмалуу иретин көрсөткөн. Бирок тегиз модель алдыңкы-чекиттин жылышын системалуу түрдө төмөн баалаган. Ок-симметриялуу эсептөөлөр, айрыкча Volodin жана кесиптештеринин чоң көлөмдүү экспериментиндеги траекторияга кыйла жакын жыйынтык берген.

Бул ок-симметриялуу модель “толук үч өлчөмдүү туура чечим” дегенди билдирбейт. Чыныгы эксперименттик жалын үч өлчөмдүү объект. Булак ок-симметриялуу геометрия учурдагы эксперимент менен салыштыруу үчүн жакшыраак жакындашуу экенин, ал эми чыныгы 3Б таасирлерди көтөрүүчү күч, кошумча масса, сүйрөө жана форма салымдарына ажыратуу үчүн толук үч өлчөмдүү симуляция керек экенин белгилейт.

250 msдеги температура талаасы

Булактын 6-сүрөтүндө төрт абалдын баарында баштапкы компакттуу өзөк 250 msде үстү кең ысык калпакча жана ылдый созулган wake түзүлүшүнө айланган. Тороидалдык куюн үстүңкү чек араны өйдө жана капталга алып барганда, төмөнкү чек араны ичкери жыгып калпакча формасын түзгөн.

Соре күйүк абалдар чоңураак жана жогору көтөрүлгөн ысык аймактарды жараткан. Радиация продукт температурасын жана өзөктүн өлчөмүн бир аз азайткан; бирок алгачкы жана орто көтөрүлүү этабында радиациянын траекторияга таасири Соре таасиринен кичине.

OH талаасы эмне үчүн маанилүү?

Температура талаасы ысык продукттардын кең аймагын көрсөтсө, OH масса үлүшү активдүү реакция аймагын тандалма көрсөтөт. Булактын 7-сүрөтүндө олуттуу OH негизинен үстүңкү бетте жана эки каптал ийинде жайгашкан, ал эми төмөнкү ысык чек аранын чоң бөлүгү химиялык жактан алсыз же активдүү эмес.

Бул натыйжа көтөрүлүп жаткан түзүлүш бардык багытта бирдей жылган жабык сфералык жалын фронту эмес экенин көрсөтөт. Төмөнкү ок боюнча куюн тартып кирген муздак жаңы аралашма төмөнкү реакция аймагын муздатып жана созот; ошентип активдүү үстүңкү калпакча менен негизинен активдүү эмес төмөнкү wake ортосунда ачык түзүлүш өнүгөт.

Үстүңкү калпакча кантип сандык бааланды?

Үстүңкү жалын чек арасы окко туураланган эллипс менен fit кылынган:

\[ \frac{(x-x_c)^2}{a_r^2} + \frac{(z-z_c)^2}{a_z^2} = 1. \tag{26} \]

Ок-симметриялуу учурда \(x\) ордуна радиалдык координата \(r\) колдонулат. \(a_r\) каптал, \(a_z\) тик жарым огу.

Форма катышы:

\[ \chi=\frac{a_r}{a_z}. \tag{27} \]

\(\chi=1\) меридиан тегиздигинде тегерек үстүңкү калпакчаны, \(\chi>1\) болсо горизонталдык багытта кеңейген жалпак калпакчаны көрсөтөт.

Булакта форма катышы убакыттын ордуна алдыңкы-чекит бийиктиги \(z_{lp}\) менен салыштырылган. Ошентип, ар башка көтөрүлүү ылдамдыктары бар абалдар окшош тик өнүгүү деңгээлдеринде салыштырылган. Жалпысынан өзөк көтөрүлгөн сайын \(\chi\) көбөйүп, сфералык жакындашуу бара-бара азыраак ылайыктуу болгон.

Форма катышынын убакыт боюнча өзгөрүшү:

\[ \frac{d}{dt} \ln\left(\frac{a_r}{a_z}\right) = \frac{\dot a_r}{a_r} - \frac{\dot a_z}{a_z}. \tag{28} \]

Демек, форма катышынын өсүшү салыштырмалуу каптал өсүү салыштырмалуу тик өсүүдөн ылдамыраак экенин билдирет.

Соре жана радиация формага бир багытта таасир этеби?

Жок. Бирдей алдыңкы-чекит бийиктигинде радиация адатта \(a_r/a_z\) катышын жогорулатып, калпакчанын жалпагыраак болушуна алып келген. Соре диффузиясы болсо өзөктүн жалпы өлчөмүн жана көтөрүлүшүн жогорулатканына карабастан форма катышын көбүнчө төмөндөтүп, үстүңкү калпакчаны тегерегирээк сактаган.

Бул натыйжа маанилүү айырманы көрсөтөт: алдыңкы-чекит бийиктигин же ылдамдыгын гана кароо жалын морфологиясын толук түшүндүрбөйт. Бирдей бийиктикке жеткен эки өзөк ар башка туурасы-бийиктиги катышына, ийрилигине жана химиялык активдүүлүк бөлүштүрүлүшүнө ээ болушу мүмкүн.

Көтөрүлүп Жаткан Суутек Жалын Өзөгү Эмне Үчүн Бөлүнөт?

Өзөк көтөрүлгөн ысык продукттар түзгөн тороидалдык куюн үстүңкү активдүү калпакча менен куюндун айланасына ийилген төмөнкү жалын бутагынын ортосундагы тар реакция аймагын катуу созгондуктан бөлүнөт; өсүп жаткан ийрилик жана созулуу реактантардын ысык аймакта калуу убактысын азайтып, конвективдүү муздоону күчөткөндө химиялык жылуулук бөлүнүшү бул жоготууларды толуктай албай, байланыш аймагында жергиликтүү өчүү пайда болот.

Бөлүнүү убактысынын аныктамасы

Бөлүнүү убактысы \(t_b\), мурда үзгүлтүксүз болгон жана үстүңкү калпакча менен куюн орогон төмөнкү бутакты бириктирген реакция аймагы биринчи жолу байланышын жоготкон учур катары аныкталган. Убакыт болжол менен 1 ms тактык менен көзөмөлдөнгөн.

Радиация өчүк болгондо Соре диффузиясы бөлүнүүнү 306 msден 327 msге, башкача айтканда 21 ms кечиктирген. Радиация күйүк болгондо ошол эле таасир 282 msден 305 msге, башкача айтканда 23 ms кечигүү. Булак муну байланышкан жалын түзүлүшүнүн өмүрүндө болжол менен %7–8 узаруу катары берет.

Радиация тескери багытта иштеген. Соре өчүк болгондо бөлүнүүнү 306 msден 282 msге, 24 ms эртелеткен; Соре күйүк болгондо 327 msден 305 msге, 22 ms эртелеткен.

Эки механизм бир учурда күйүк болгондо 305 ms, экөө тең өчүк болгондо 306 ms табылган. Бул 1 ms айырма изилдөөдө чечилген физикалык өзгөрүү катары кабыл алынбайт. Эки абал бөлүнүү убактысы боюнча жакын болсо да, температура, курам, өзөк өлчөмү жана көтөрүлүү тарыхы бирдей эмес.

Бөлүнүү толук өчүүнү билдиреби?

Жок. Биринчи байланыш үзүлгөндөн кийин үстүңкү калпакчада да, ийилген төмөнкү бутакта да химиялык энергия бөлүнүшү уланат. Баштапкы жалгыз өзөк жоголуунун ордуна бир нече реакция түзүлүшүнө айланат. Булак бул бөлүктөр кийин термиялык жана химиялык wake ичинде жылып, өз ара аракеттениши же экинчи өзөктөрдү түзүшү мүмкүн экенин билдирет.

Бөлүнүү жолун агымбы, Соре менен радиациябы аныктайт?

Сапаттык бөлүнүү ордун жана жолун негизинен гидродинамикалык агым топологиясы аныктайт. Төрт моделдин баарында бирдей негизги тороидалдык куюн–жалын байланышы байкалган. Соре диффузиясы жана радиация негизинен жалын куюндан келип чыккан созулууга жана муздоого канчалык чыдай аларын жана критикалык шартка качан жетерин өзгөрткөн.

Турбуленттүүлүк жагынан эмне айтууга болот?

Бул изилдөө турбуленттүү жалын симуляциясы эмес. Агым ламинардуу жана бир үстөм тороидалдык куюн тарабынан башкарылат. Ошентсе да жергиликтүү жалын–куюн өз ара аракеттенүүсүндө беттик деформация, ийрилик, созулуу, реактантардын калуу убактысы, конвективдүү ташуу жана жергиликтүү өчүү сыяктуу турбуленттүү жалын–eddy өз ара аракеттеринде да кездешкен физикалык процесстер байкалат.

Ошондуктан изилдөө агым башкарган жергиликтүү өчүү чеги Соре ташылышы жана жылуулук жоготуусу менен өзгөрүшү мүмкүн экенин көрсөтөт. Бирок турбуленттүү суутектин күйүү чектери сандык түрдө эсептелди деген жыйынтык чыгарууга болбойт; бул үчүн ар кандай өлчөмдөргө жана созулуу катыштарына ээ көптөгөн куюндарды камтыган үч өлчөмдүү симуляциялар талап кылынат.

Изилдөө колдогон жыйынтыктар

• %6 көлөмдүк H₂–аба аралашмасында алгачкы өзөк өсүшү изилденген убакыт аралыгында болжол менен \(R_f^2\simeq Kt\) режимин көрсөткөн.
• Соре диффузиясы \(K\)-ны жогорулатып, өзөктүн өсүшүн ылдамдаткан.
• Радиациялык жылуулук жоготуусу \(K\)-ны төмөндөтүп, продукттарды муздаткан.
• Соре таасири чоңураак жана ылдамыраак көтөрүлгөн ысык продукт өзөгүн түзгөн.
• Ок-симметриялуу эсептөөлөр бар эксперименттик алдыңкы-чекит траекториясын тегиз эсептөөлөргө караганда жакшыраак кайра чыгарган.
• Радиация бирдей бийиктикке карата үстүңкү калпакчанын жалпагыраак болушун колдогон.
• Соре диффузиясы үстүңкү калпакчанын тегерегирээк калышын колдогон.
• Бөлүнүү бардык абалдарда тороидалдык куюн жараткан жергиликтүү созулуу жана өчүү аркылуу болгон.
• Соре бөлүнүүнү кечиктирген, радиация эртелеткен.
• Соре менен радиациянын бөлүнүү убактысына карама-каршы таасирлери экөө тең күйүк абалда дээрлик бири-бирин теңдеген.

Изилдөө колдобогон жыйынтыктар

• Изилдөө %6 H₂ аралашмасынын бардык реалдуу объект геометрияларындагы жүрүм-турумун аныктабайт.
• Эсептелген алдыңкы-чекит ылдамдыгы ламинардык күйүү ылдамдыгына түздөн-түз тең эмес.
• \(R_f^2\simeq Kt\) режими чексиз убакыт үчүн универсалдуу өсүү мыйзамы катары далилденген эмес.
• Ок-симметриялуу модель толук үч өлчөмдүү жалын бөлүнүшүн чечпейт.
• Изилдөө турбуленттүү жалын же турбуленттүү күйүү чеги эсептөөсү эмес.
• Оптикалык жука радиация модели деталдуу спектралдык радиациялык берүүнү чечпейт.
• Изилдөөнүн өзүндө жаңы эксперимент жүргүзүлгөн эмес; эксперименттик маалыматтар мурдагы басылмалардан салыштыруу максатында алынган.
• Бөлүнүү толук күйүү өчүшүн билдирбейт.

Негизги методологиялык чектөөлөр

Тегиз жана ок-симметриялуу эки өлчөмдүү геометриялар чыныгы үч өлчөмдүү жалындын ок-симметриялуу эмес деформацияларын көрсөтө албайт. Айрыкча бөлүнүүдөн кийин ар башка жалын бөлүктөрүнүн үч өлчөмдүү өз ара аракеттери жана бир нече куюндун бир убактагы таасири камтылбайт.

Агым ламинардуу бойдон калат деп болжолдонуп, турбуленттүүлүк же subgrid күйүү модели колдонулган эмес. Ошондуктан натыйжаларды түздөн-түз турбуленттүү суутек булуттарына колдонууга болбойт.

Радиация H₂O эмиссиясы жана оптикалык жука ыкма шартында гана моделденген. Өзүн-өзү сиңирүү эсепке алынбагандыктан, таза радиациялык жоготуу оптикалык калыңдык өскөн шарттарда ашыкча бааланышы мүмкүн.

Бар эксперименттер менен салыштыруу алдыңкы-чекит кыймылы жагынан гана жүргүзүлгөн жана эксперименттердин геометриялары да бири-биринен айырмаланат. Чоң көлөмдүү эксперимент болжол менен эркин өзөктү көрсөтсө, экинчи маалымат топтому жабык түтүктөн алынган. Геометриялык чектөө жана дубал өз ара аракеттери айрыкча кийинки этаптарда маанилүү боло алат.

Булак жана Ыкма Жөнүндө Эскертүү

Түпнуска изилдөө: Effects of Soret Diffusion and Radiative Heat Loss on the Evolution of Buoyant Flame Kernels in Ultra-Lean Hydrogen-Air Mixture.

Авторлор: Ivan S. Yakovenko; Alexey D. Kiverin.

Жооптуу автор: Ivan S. Yakovenko.

Мекеме: Joint Institute for High Temperatures of the Russian Academy of Sciences, Moscow, Russia.

Журнал жана басылма идентификациясы: Fire, 2026, Том 9, Макала 383. Басмакана: MDPI. DOI: 10.3390/fire9090383.

Басылма процесси: Макала 9 август 2026-жылы алынган, 29 август 2026-жылы кайра каралган, 4 сентябрь 2026-жылы кабыл алынган жана 5 сентябрь 2026-жылы жарыяланган.

Булак түрү: Рецензияланган, ачык жеткиликтүү, деталдуу сандык күйүү жана суюктуктар динамикасы изилдөө макаласы. Изилдөөнүн өзү жаңы физикалык күйүү экспериментин камтыбайт.

Лицензия: Creative Commons Attribution (CC BY).

Негизги ыкма: %6 көлөмдүк H₂–аба аралашмасы; деталдуу Kéromnès суутек кычкылдануу механизми; төмөн-Mach газ динамикасы; аралашма-орточо молекулалык ташуу; H жана H₂ үчүн кыскартылган Соре модели; оптикалык жука H₂O радиация модели; 1B сфералык, 2B тегиз жана 2B ок-симметриялуу симуляциялар; ички NRG төмөн-Mach FDS чечкичи.

Grid: Негизги 2B эсептөөлөрдө 0,4 mm бирдей уяча аралыгы жана чечилген жарым аймакта 128 × 256 уяча. Жакындашуу текшерүүсү 0,1, 0,2 жана 0,4 mm gridдер менен жүргүзүлгөн.

Эксперименттик текшерүү контексти: Сандык алдыңкы-чекит траекториялары Volodin жана кесиптештери менен Anikin жана кесиптештеринин мурда жарыяланган эксперименттик маалыматтары менен салыштырылган. Эксперименттер бул изилдөөнүн алкагында жүргүзүлгөн эмес.

Каржылоо: Изилдөө Russian Federation Ministry of Science and Higher Education тарабынан State Assignment No. 075-00270-26-00 алкагында колдоого алынган.

Маалыматтардын жеткиликтүүлүгү: Авторлор макаланын жыйынтыктарын колдогон чийки маалыматтар суроо-талап боюнча берилерин билдиришкен.

Кызыкчылыктардын кагылышы: Авторлор кызыкчылыктардын кагылышы жок экенин билдиришкен.

Авторлордун салымдары: Концептуалдаштыруу, ыкма, программалык камсыздоо, текшерүү, изилдөө, маалыматтарды башкаруу, алгачкы долбоор жана визуалдаштыруу Ivan S. Yakovenko тарабынан аткарылган. Формалдык талдоо Ivan S. Yakovenko жана Alexey D. Kiverin тарабынан; ресурстар, карап чыгуу жана редакциялоо, көзөмөл, долбоорду башкаруу жана каржылоону камсыздоо эки автор тарабынан биргеликте жүргүзүлгөн.

Илимий чечмелөө чеги

Изилдөөнүн негизги жыйынтыгы — өтө-жарды суутек жалындарын Жердин тартылуу күчү шартында бир гана “күйүү ылдамдыгы” менен жетиштүү түшүндүрүү мүмкүн эмес. Молекулалык ташуу алгачкы өсүүгө, ысык продукт көлөмү көтөрүүчү күчкө, көтөрүүчү күч өз куюн талаасына, куюн болсо жалын морфологиясына жана жергиликтүү өчүүгө таасир этет. Соре диффузиясы жана радиация бул чынжырдын убакыт масштабын жана жалындын созулууга чыдамдуулугун өзгөртөт; бирок бөлүнүүнүн негизги топологиялык жолу изилденген төрт абалда тең гидродинамикалык жактан окшош.

Натыйжаларды реалдуу үч өлчөмдүү, ок-симметриялуу эмес же турбуленттүү суутек өрттөрүнө түздөн-түз колдонуу туура эмес. Авторлордун жыйынтык бөлүмүндө да бул камтуу үчүн толук үч өлчөмдүү симуляциялар зарыл экени белгиленген.