Ин таҳқиқот рафтори гармшавӣ, оксидшавӣ ва аз даст рафтани мавод аз сатҳи карбиди кремнийро (SiC) ҳангоми парвози гиперсадоӣ тавассути якҷо кардани симулятсияҳои реактивии динамикаи молекулавӣ бо ҳалкунандаи аэротермалии ATAC-1 меомӯзад. Муҳаққиқон гузаронандагии гармӣ, гармиғунҷоиши хос, реаксияҳои сатҳӣ ва намудҳои бухоршавандаро, ки бо LAMMPS ва ReaxFF ҳисоб шудаанд, ба ҷадвали аблятсияи B-prime интиқол додаанд; сипас ин ҷадвалро барои ҳисоб кардани ҳарорат ва ақибнишинии сатҳи геометрияи бинии конусӣ дар шароити Mach 5, 10, 15 ва 20 истифода бурдаанд. Симулятсияҳо нишон медиҳанд, ки бо зиёд шудани ҳарорат гузаронандагии гармии SiC кам, гармиғунҷоиши хоси он зиёд мешавад ва химияи сатҳӣ дар болотар аз тақрибан 1500 °C метавонад намудҳои бухоршавандаи мураккабтар тавлид кунад. Бо вуҷуди ин, таҳқиқот preprint-и аз peer review нагузашта аст; ҳалкунандаи ATAC-1-и истифодашуда ҳисоби пурраи Navier–Stokes-ро анҷом намедиҳад ва баъзе натиҷаҳои Mach-и баланд қиматҳои талафоти сатҳиро тавлид мекунанд, ки аз геометрияи физикӣ зиёд буда, худи муаллифон низ онҳоро беэътибор мешуморанд.

Навоварии асосии таҳқиқот кӯшиши интиқоли мустақими иттилооти химиявӣ ва термомеханикии сатҳи атомӣ ба воситаи зудкори тарҳрезии муҳандисӣ мебошад. Ҳамин тавр, ба ҷойи ҷадвалҳои эмпирикие, ки танҳо аз таҷрибаҳои ҳарорати хонагӣ ҳосил шудаанд, ҳадаф истифодаи натиҷаҳои динамикаи молекулавӣ мебошад, ки метавонанд равандҳои канда шудани пайванд ва ташаккули молекулаҳои навро намояндагӣ кунанд. Аммо усул ҳанӯз воситаи тасдиқшуда ё тарҳрезии ниҳоӣ нест, балки ба таъбири муҳаққиқон марҳилаи аввал ва намоиши имконпазирӣ мебошад.

Саволи асосии таҳқиқот чист?

Саволи асосии таҳқиқот муайян кардани он аст, ки ҳангоми парвози гиперсадоӣ аз сатҳи карбиди кремний кадом намудҳои атомӣ ва молекулавӣ ҷудо мешаванд ва ин равандҳои химиявиро чӣ гуна ба ҳарорати сатҳ ва ақибнишинии маводи воситаи парвозӣ интиқол додан мумкин аст. Муҳаққиқон махсусан месанҷанд, ки оё равандҳои оксидшавӣ ва кандашавии пайванд дар миқёси атомӣ метавонанд ҷадвалҳои B-prime-ро, ки дар ҳисобҳои аэротермалии миқёси калон истифода мешаванд, таъмин кунанд ё не.

Дар парвози гиперсадоӣ бинӣ, канори пеши бол ва дигар сатҳҳои беруна аз сабаби мавҷи зарба, кашиши аэродинамикӣ, соиши сатҳ ва қабати сарҳадии ҳарорати баланд ба бори шадиди гармӣ дучор мешаванд. Системаи муҳофизати гармӣ як қисми ин энергияро дар худ нигоҳ медорад, қисми дигар бошад тавассути вайроншавӣ, газшавӣ ё дур шудани мавод аз сатҳ масраф мешавад. Дар таҳқиқот кам шудани ҷараёни гармӣ ба девор аз ҳисоби омехта шудани газҳои аз сатҳ хориҷшуда ба қабати сарҳадӣ “таъсири боздоранда” ё “blocking effect” номида мешавад.

Дар моделҳои мавҷуда кадом холигӣ ҳадаф гирифта шудааст?

Моделҳои маъмули B-prime ба ҷадвалҳое такя мекунанд, ки фарз мекунанд химияи сатҳӣ ба ҳолати мувозинат мерасад. Ин равиш аз ҷиҳати ҳисоббарорӣ зуд аст; аммо метавонад реаксияҳои кӯтоҳмуддат, metastable ва аз ҷиҳати кинетикӣ идорашавандаро дар қабати сарҳадии зарба пурра намояндагӣ накунад. Моделҳои химияи finite-rate муфассалтаранд, вале таҷрибаҳои гарон ва механизмҳои густурдаи реаксияро талаб мекунанд. Дар таҳқиқот таъкид мешавад, ки байни моделҳои мавҷуда дар ташаккули намудҳои гуногуни сатҳӣ фарқҳое зиёда аз 300 фоиз гузориш шудаанд.

Пешниҳоди муҳаққиқон ҳисоб кардани реаксияҳои химиявӣ дар сатҳи динамикаи классикии молекулавӣ ва сохтани ҷадвали зудкори B-prime аз ин натиҷаҳо мебошад. Ин ҷадвал на ба ҳалкунандаи ҷараёни дақиқии баланд, балки ба коди ATAC-1, ки барои тарҳрезии пешакӣ ва таҳқиқоти параметрӣ истифода мешавад, интиқол дода шудааст. Ҳамин тавр, таҳқиқот байни химияи атомистӣ ва ҳалли аэротермалии миқёси муҳандисӣ занҷири интиқоли додаи сустпайвастро ташкил мекунад.

Дар қабати сарҳадии зарба кадом равандҳо ҳамзамон ба амал меоянд?

Дар Шакли 1-и таҳқиқот ҳамкориҳои асосӣ дар наздикии сатҳи гиперсадоӣ якҷо нишон дода шудаанд. Пас аз зарбаи ба вуҷудовардаи ҷараёни озод молекулаҳои газ dissociate мешаванд, сатҳҳои энергияи дохилии онҳо ангезиш меёбанд ва бо сатҳ ба ҳамкории химиявӣ медароянд. Дар ҳамин минтақа радиатсия, хориҷшавии газ аз мавод, диффузияи намудҳои дорои карбон ва вайроншавии зарба ба амал меоянд. Ҳамаи ин равандҳо дар як миқёси вақт ва дарозӣ нестанд: кандашавии пайванд дар сатҳи femtosecond ё nanosecond рух медиҳад, дар ҳоле ки фарсудашавии сатҳ ва тағйири геометрӣ метавонанд дар тӯли сонияҳо ҷамъ шаванд.

Модели атомистӣ чӣ гуна сохта шудааст?

Рафтори термомеханикӣ ва оксидшавии карбиди кремний тавассути симулятсияҳои классикии динамикаи молекулавӣ, ки бо нармафзори Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) иҷро шудаанд, моделсозӣ шудааст. Барои иҷозат додан ба кандашавии пайвандҳои химиявӣ, ташаккули пайвандҳои нав ва интиқоли заряд байни атомҳо майдони қувваи реактивии ReaxFF ва модели мувозинати заряд QEq истифода шудаанд.

Энергияи умумии потенсиалии системаи ReaxFF дар таҳқиқот бо ҷузъҳои зерин ифода мешавад:

\[ U_{\mathrm{ReaxFF}} = U_{\mathrm{bond}} + U_{\mathrm{over}} + U_{\mathrm{under}} + U_{\mathrm{valence}} + U_{\mathrm{lp}} + U_{\mathrm{torsion}} + U_{\mathrm{vdW}} + U_{\mathrm{Coul}} \]

• U_bond: Ҳиссаи энергияи пайванди атомҳои пайвастшуда.
• U_over: Координатсияи аз ҳад зиёди ғайрифизикии атомҳоро маҳдуд мекунад.
• U_under: Ба устувории атомҳои нокифоя ҳамоҳангшуда саҳм мегузорад.
• U_valence: Ҳади энергия вобаста ба кунҷҳои пайванд.
• U_lp: Энергияи ҷуфти танҳоии электрон.
• U_torsion: Энергияе, ки аз кунҷҳои печиши чоратомӣ ба вуҷуд меояд.
• U_vdW: Ҳамкориҳои кӯтоҳмасофаи van der Waals-ро намояндагӣ мекунад.
• U_Coul: Ҳамкории дарозмасофаи Coulomb.

Ин энергияи умумӣ дар динамикаи молекулавӣ барои ҳисоб кардани қувваҳои ба ҳар атом таъсиркунанда ва ҳаракати онҳо дар тӯли вақт истифода мешавад. Аҳамияти модел дар он аст, ки бар хилофи потенсиалҳои анъанавии дорои сохтори собити пайванд, ба ташаккули пайвандҳои нави Si–O ва C–O ҳангоми оксидшавӣ иҷозат медиҳад.

Гузаронандагии гармӣ чӣ гуна ҳисоб шудааст?

Гузаронандагии гармии сохторҳои поликристаллӣ ва монокристаллии α(6H)-SiC бо муносибати Green–Kubo ҳисоб шудааст:

\[ k_{ij} = \frac{V}{k_B T^2}\int_0^\infty \langle J_i(0)J_i(t)\rangle\,dt \]

Дар ин ҷо k_ij тензори гузаронандагии гармӣ, V ҳаҷми ҳуҷайраи симулятсия, k_B собити Boltzmann, T ҳарорати мутлақ ва J ҷараёни гармиро ифода мекунанд. Ҳади дохили қавсҳои кунҷӣ автокоррелятсияи байни ҷараёни гармии ибтидоӣ ва ҷараёни гармӣ дар вақти баъдиро нишон медиҳад. Нигоҳ доштани коррелятсия барои муддати дароз маънои интиқоли самараноктари энергия дар дохили сохторро дорад.

Дар ҳисобҳои термомеханикӣ ҳуҷайраи симулятсионии тақрибан 20 × 20 × 50 nm истифода шудааст. Дар самтҳои X ва Y шартҳои сарҳадии даврӣ татбиқ шуда, энергияи гармӣ тавре ки дар Шакли 2 нишон дода шудааст, аз минтақаи гарм ба минтақаи сард интиқол ёфтааст. Рангҳои гуногун дар тасвир донаҳои алоҳидаи сохтори поликристаллиро нишон медиҳанд.

Ҳисобҳои гармиғунҷоиши хос ва фишор

Барои гармиғунҷоиши хос усули миёнагирии гармоникӣ хариташуда (harmonically mapped averaging, HMA) дар дохили LAMMPS истифода шудааст. Ифодаҳои энергия, фишор ва гармиғунҷоиши хос дар ҳаҷми собит дар мақола чунинанд:

\[ \langle U\rangle_{\mathrm{HMA}} = \frac{d}{2}(N-1)k_B T + \left\langle U+\frac{1}{2}\mathbf{F}\cdot\Delta\mathbf{r}\right\rangle \]

\[ \langle P\rangle_{\mathrm{HMA}} = \Delta\hat{P}+\left\langle P_{\mathrm{vir}}+\frac{\beta\Delta\hat{P}-\rho}{d(N-1)}\mathbf{F}\cdot\Delta\mathbf{r}\right\rangle \]

\[ \langle C_V\rangle_{\mathrm{HMA}} = \frac{d}{2}(N-1)k_B+\frac{\langle U_{\mathrm{HMA}}^2\rangle-\langle U_{\mathrm{HMA}}\rangle^2}{k_BT^2}+\frac{1}{4T}\left\langle \mathbf{F}\cdot\Delta\mathbf{r}+\Delta\mathbf{r}\cdot\Phi\cdot\Delta\mathbf{r}\right\rangle \]

Дар ин муодилаҳо d шумораи андозаҳои система, N шумораи атомҳо, U энергияи умумӣ, F қувваҳои атомӣ, Δr ҷойивазкунии атомҳо аз мавқеъҳои шабака, Φ матритсаи Hessian, ρ зичии ададӣ ва P_vir фишори virial-ро нишон медиҳанд. Андозаи лаппиши энергия ба ҳисоби гармиғунҷоиши хос дар ҳаҷми собит мустақиман саҳм мегузорад.

Симулятсияи оксидшавӣ ва аблятсия чӣ гуна иҷро шудааст?

Барои модели оксидшавӣ ҳуҷайраи ҷудогона бо андозаи уфуқии 20 × 20 nm ва андозаи амудии 15 nm сохта шудааст. Дар ҳуҷайра қабати SiC бо ғафсии 10 nm ва болои он фазои холии 5 nm, ки қабати сарҳадиро намояндагӣ мекунад, мавҷуд аст. Оксигени молекулавӣ ва атомӣ бо суръат ва ҳароратҳои гуногун ба сатҳи SiC равона шудааст.

Дар Шакли 3 атомҳои сурхи қисми поён собит карда шудаанд, ҳарорати атомҳои сабз бо термостат идора шудааст ва атомҳои оксиген аз қисми боло ба система илова шудаанд. Дар болои ҳуҷайра минтақае таъриф шудааст, ки намудҳои аз сатҳ ба қадри кофӣ дуршуда аз симулятсия нест карда мешаванд. Ҳамин тавр, бархӯрди оксиген ба сатҳ, ташаккули пайванди химиявӣ ва ҷудо шудани баъзе маҳсулот аз сатҳ дар миқёси атомӣ пайгирӣ шудааст.

Қувваи умумӣ дар термостати Langevin чунин ҷудо карда шудааст:

\[ \mathbf{F} = \mathbf{F}_c+\mathbf{F}_f+\mathbf{F}_r \]

\[ \mathbf{F}_f = -\left(\frac{m}{\tau_{\mathrm{damp}}}\right)\mathbf{V} \]

\[ \mathbf{F}_r = \frac{k_BT m}{dt\,\tau_{\mathrm{damp}}} \]

F_c қувваи консервативӣ аз потенсиали ReaxFF, F_f қувваи соиш ва F_r ҳади қувваи тасодуфиро ифода мекунад. m массаи атом, V суръат, τ_damp вақти оромиши ҳарорат ва dt қадами вақт мебошад. Муодилаи қувваи тасодуфӣ дар ин ҷо айнан ба шакле оварда шудааст, ки дар мақола навишта шудааст.

Динамикаи молекулавӣ кадом намудҳои бухоршавандаро пешгӯӣ кардааст?

Дар симулятсияҳо гузориш шудааст, ки дар поён аз 1500 °C намудҳои асосии аз сатҳи SiC ҷудошаванда оксиди карбон (CO), оксигени атомӣ (O) ва гази карбон (CO₂) мебошанд. Дар ҳароратҳои баландтар пешгӯӣ шудааст, ки намудҳои мураккабтар ва metastable, аз қабили C₂O₂, C₂O₃ ва SiO₃, метавонанд ба вуҷуд оянд.

Баръакс, ҷадвали ибтидоии B-prime, ки ба ATAC-1 интиқол дода шудааст, бо мақсади кам кардани бори ҳисоббарорӣ танҳо бо намудҳои CO, SiO₂ ва CO₂ маҳдуд шудааст. Аз ин рӯ, ҳамаи маҳсулоти реаксияе, ки дар динамикаи молекулавӣ мушоҳида шудаанд, ба ҳисоби ақибнишинии сатҳи миқёси макро интиқол дода нашудаанд. Реаксияҳои эҳтимолӣ бо нитроген ва металлҳои гузаришие, ки дар керамикаи ҳарорати ултрабаланд истифода мешаванд, низ дар версияи ҷорӣ вуҷуд надоранд.

ATAC-1 чӣ гуна кор мекунад?

Гарчанде ATAC-1 дар таҳқиқот баъзан ҳамчун коди CFD номида мешавад, он муодилаҳои пурраи Navier–Stokes-ро ҳал намекунад. Код хатҳои ҷараёни inviscid-ро, ки аз ҷараёни пас аз зарба оғоз мешаванд, пайгирӣ карда, тақсимоти фишор ва суръатро ба даст меорад; сипас аз наздиккуниҳои қабати сарҳадӣ, коррелятсияҳои нимэмпирикӣ ва ҳалкунандаи гармии ба фарқҳои ниҳоӣ асосёфта истифода мекунад. Ҳарорати сатҳ ва суръати ақибнишинии мавод тавассути ҳал кардани тавозуни ғайрихаттии энергия дар ҳар нуқтаи геометрӣ бо усули Newton–Raphson муайян мешаванд.

Ин равиш аз ҳалкунандаҳои CFD-и дақиқии баланд тезтар аст; аммо ҷузъиёти ҷараёни часпак, турбулентӣ, ҳамкориҳои мураккаби зарба ва минтақаҳои ҷараёни ҷудошударо мустақиман ҳал карда наметавонад. Ҳадафи таҳқиқот низ на пешгӯии парвоз бо дақиқии ниҳоӣ, балки нишон додани он аст, ки ҷадвали нави B-prime метавонад ба воситаи зудкори муҳандисӣ ҳамгиро карда шавад.

Барои конуси бинии гиперсадоӣ кадом шароитҳо омӯхта шудаанд?

Дар ҳисобҳои ATAC-1 Атмосфераи Стандартии ИМА-и 1976 ва ҷадвали Clark Mollier барои ҳаво истифода шудаанд. Кунҷҳои yaw ва roll сифр қабул шудаанд. Ду баландии асосӣ омӯхта шудааст:

• 15.240 m: 11.664,53 Pa фишори атмосфера ва −56,5 °C ҳарорат.
• 21.336 m: 4.487,68 Pa фишори атмосфера ва −55,23 °C ҳарорат.

Геометрия конуси симметрии бинии дорои нимкунҷи муассири конусии 7,5° ва миқёси 10 cm мебошад. Қадами вақт 0,01 сония муайян шуда, дар сатҳ 11 нуқтаи таҳлили муфассал бо номи P1–P11 интихоб шудаанд. Шакли 4 нишон медиҳад, ки ин нуқтаҳо дар нӯги бинӣ ва сатҳҳои бодхӯрдаву пуштибоди конус чӣ гуна паҳн шудаанд.

Силсилаи асосии параметрӣ дар суръатҳои Mach 5, 10, 15 ва 20 иҷро шудааст. Кунҷи ҳамла дар баробари 5° аз 10°–20° бо фосилаи якдараҷаӣ тағйир дода шуда, давомнокии парвози ҳисобҳои асосӣ 32 сония гирифта шудааст. Илова бар ин, барои Mach 5 ва Mach 10 дар кунҷи ҳамлаи 19° сенарияҳои дарозмуддати 900 сония иҷро шудаанд.

Ҳарорати рукуд чӣ гуна тахмин шудааст?

Дар таҳқиқот барои баландӣ ва Mach-ҳои гуногун муносибати содакардашудаи ҳарорати рукуд низ истифода шудааст:

\[ T_0(h,M)=T_\infty(h)\left[1+Pr^{1/3}\left(\frac{\gamma-1}{2}\right)M^2\right] \]

T₀ ҳарорати рукуд, T_∞ ҳарорати ҷараёни озод дар баландии дахлдор, M рақами Mach, Pr рақами Prandtl ва γ таносуби гармиҳои хосро ифода мекунанд. Дар ҳисоб ҳаво гази идеалӣ қабул шуда, қимати γ дар 1,4 собит нигоҳ дошта шудааст. Ин муносибати сода нишон медиҳад, ки энергияи кинетикии суръат ҳангоми боздоштани ҷараён то чӣ андоза метавонад ба энергияи гармӣ табдил ёбад; аммо реаксияҳои воқеии сатҳӣ, рафтори фаза дар фишори баланд ва химияи берун аз мувозинатро танҳо худаш шарҳ дода наметавонад.

Натиҷаҳои хосиятҳои гармӣ чиро нишон доданд?

Дар Шакли 5 гузаронандагии гармии α-SiC, ки бо MD ҳисоб шудааст, дар тақрибан 750 °C аз 48 W/m·°C то дар атрофи 2250 °C тақрибан 13 W/m·°C коҳиш меёбад. Коҳиш дар фосилаи аввали ҳарорат тезтар, пас аз тақрибан 1600 °C бошад ҳамвортар аст. Муҳаққиқон лаппишҳои хурди наздики 1250 °C ва 1500 °C-ро ба остонаҳои оксидшавӣ ё гузариши фаза вобаста мекунанд, вале эътироф мекунанд, ки ба натиҷаҳо андозаи симулятсия, шартҳои сарҳадӣ, вақти оромиши термостат ва таъсирҳои квантии сарфи назаршуда таъсир расонда метавонанд.

Қиматҳои гармиғунҷоиши хос дар Шакли 6 аз тақрибан 750 °C то 2250 °C мунтазам зиёд шуда, дар график аз тақрибан 5,4 ба 6,15 мерасанд. Аммо гарчанде меҳвари амудии график ҳамчун “J/kg·°C” нишон дода шудааст, миқёси ададии пешниҳодшуда бо ин воҳид ба таври равшан шарҳ дода нашудааст. Аз ин рӯ, пеш аз истифодаи мустақими қиматҳои график ҳамчун додаи мутлақи муҳандисӣ санҷиши воҳид ва миқёс лозим аст.

Суръат, баландӣ ва кунҷи ҳамла ба ҳарорат чӣ гуна таъсир карданд?

Дар Шакли 7 тамоюли умумӣ он аст, ки бо баланд шудани рақами Mach ҳарорати максималии сатҳ зиёд мешавад. Силсилаи баландии 15.240 m ва Mach 20 дар аксари кунҷҳои ҳамла ба ҳадди сублиматсияи 2850 °C, ки дар модел барои SiC истифода шудааст, наздик мешавад ё аз он мегузарад. Аз ҳама равшантар гузаштан дар кунҷи ҳамлаи 17° дида мешавад ва ҳарорати максималӣ дар график ба тақрибан 3200 °C мерасад.

Дар баландии баландтари 21.336 m аз сабаби зичии пасти атмосфера ҳароратҳои Mach 5 аз натиҷаҳои 15.240 m пасттаранд. Аммо дар суръатҳои баланд, аз қабили Mach 15 ва Mach 20, ҳарорат боз ба диапазони 2600–2850 °C мерасад. Аз ин рӯ, таъсири баландиро танҳо арзёбӣ кардан мумкин нест; онро бояд якҷо бо суръат, фишори сатҳ, кунҷи ҳамла ва наздиккунии ҷараёни истифодашуда баррасӣ кард.

Чаро натиҷаҳои ақибнишинии мавод мушкил доранд?

Шакли 8 дар баробари тамоюлҳои физикӣ натиҷаҳоеро низ тавлид кардааст, ки аз ҳадди эътибори модел берунанд. Дар силсилаи Mach 20, 15.240 m ва кунҷи ҳамлаи 17° ақибнишинии сатҳ тақрибан 3,10 × 10¹⁴ cm нишон дода шудааст. Ин қимат бо ҳеҷ геометрияи воқеии воситаи нақлиёт мувофиқат намекунад. Дар як қисми матн аномалияи 17°–18° ба Mach 15 нисбат дода мешавад, дар ҳоле ки легендаи график қимати ғайриоддиро дар силсилаи Mach 20 ва 15 km ҷойгир мекунад. Аз ин рӯ, байни баёни натиҷа ва график аз ҷиҳати рақами Mach номувофиқӣ вуҷуд дорад.

Муҳаққиқон шарҳ медиҳанд, ки ин аномалия метавонад аз он барояд, ки вақте ҳарорат ва фишор аз соҳаи калибровкаи ҷадвали B-prime берун мешаванд, ATAC-1 рафтори маводро ба таври экспоненсиалӣ extrapolate мекунад. Вақте ҳарорати сатҳ аз ҳадди сублиматсияи 2850 °C мегузарад, модел ҳатто агар маводи сахт пурра аз байн рафта бошад ҳам, ҳисоби ақибнишиниро идома медиҳад. Қадамҳои вақти пас аз тамом шудани ғафсии ибтидоии мавод маънои физикӣ надоранд.

Фарсудашавии бештар дар куҷои сатҳ ба вуҷуд омадааст?

Дар Шакли 9 харитаҳои сатҳи Mach 15 нишон медиҳанд, ки қисми зиёди талафоти мавод дар нӯги бинӣ ва тарафи бодхӯрдае, ки ба самти парвоз рӯ ба рӯст, ҷамъ мешавад. Бо тағйир ёфтани кунҷи ҳамла минтақаи ақибнишинии баланд аз маркази симметрии бинӣ ба канор ҳаракат мекунад. Ин натиҷа нишон медиҳад, ки ҳатто агар дар кунҷҳои гуногуни ҳамла ҳарорати умумии якхела ба вуҷуд ояд ҳам, тақсимоти маҳаллии бор дар рӯи сатҳ метавонад тағйир ёбад.

Аммо миқёси ранги харита дар давоми 32 сония то 80 cm ақибнишинии умумиро нишон медиҳад. Бо назардошти он ки геометрияи бинии омӯхташуда миқёси 10 cm дорад, ин қиматҳо низ нишон медиҳанд, ки маводи воқеӣ пеш аз анҷоми вақти ҳисоб пурра тамом мешуд. Ин харитаҳо бояд на ҳамчун пешгӯии миқдории умр, балки ҳамчун харитаҳои тамоюл хонда шаванд, ки нишон медиҳанд модел дар кадом минтақаҳои сатҳ бори баландро пешгӯӣ мекунад.

Сенарияҳои дарозмуддати Mach 5 ва Mach 10

Дар баландии 15.240 m, кунҷи ҳамлаи 19° ва давомнокии парвози 900 сония Mach 5 бо Mach 10 муқоиса шудааст. Дар сенарияи Mach 5 ҳарорати сатҳ ба тақрибан 900 °C расида, аз сатҳи оғози ақибнишинии тақрибан 1400 °C, ки дар таҳқиқот барои SiC нишон дода шудааст, паст мондааст. Талафоти сатҳӣ дар Шакли 10 дар миқёси хеле хурди ададӣ боқӣ мемонад.

Дар сенарияи Mach 10 бошад, ҳарорат дар болотар аз 2500 °C лаппишҳои зуд нишон дода, ақибнишинии умумии ҳисобшуда ба тақрибан 101,6 cm расидааст. Ин қимат аз геометрияи ибтидоӣ зиёд аст. Муаллифон мегӯянд, ки лаппишҳои тези ҳарорат метавонанд маҳдудиятҳои ҳалномаи қадами вақт ва ҳисоби аблятсияро инъикос кунанд; ҳисобҳои баъд аз пурра аз байн рафтани геометрия беэътиборанд.

Мисоли Mach 7 чиро нишон дод?

Ҳисоби содакардашудаи ҳарорати рукуд нишон дод, ки дар баландии 15.240 m ва суръати Mach 7 сатҳи SiC метавонад аз ҳарорати оғози аблятсия гузарад. Баръакс, дар ҳисоби ATAC-1 барои кунҷи ҳамлаи 0° ва 20 сония фишори сатҳ ба 2,43 MPa расида ва танҳо талафоти мавод дар сатҳи хатои гирдкунӣ ба вуҷуд омадааст.

Ин номувофиқӣ нишон медиҳад, ки аблятсияро танҳо бо остонаи ҳарорат муайян кардан мумкин нест. Ҳудудҳои гудозиш, бухоршавӣ ва сублиматсия дар фишори баланд бояд бо диаграммаҳои боэътимоди фаза муайян карда шаванд. Фарзияҳои соддаи гази идеалӣ ва хосиятҳои собит, ки дар таҳқиқот истифода шудаанд, ин гузаришҳоро бо тамоми ҷузъиёт намояндагӣ намекунанд.

Ҷиҳатҳои қавии таҳқиқот кадомҳоянд?

• Ҷараёни кори равшани бисёрмиқёсаеро пешниҳод мекунад, ки реаксияҳои атомистиро ба воситаи зудкори тарҳрезии аэротермалӣ интиқол медиҳад.
• ReaxFF истифода шудааст, ки метавонад кандашавии пайванд, ташаккули пайванд ва интиқоли зарядро ифода кунад.
• Гузаронандагии гармӣ, гармиғунҷоиши хос, химияи сатҳӣ, намудҳои бухоршаванда, ҳарорат ва ақибнишинии геометрӣ дар як занҷири модел баррасӣ шудаанд.
• Рақамҳои гуногуни Mach, ду баландӣ, шумораи зиёди кунҷҳои ҳамла ва ду давомнокии гуногуни парвоз муқоиса шудаанд.
• Муаллифон натиҷаҳои ғайрифизикии моделро пинҳон накарда, мушкилоти extrapolation ва тамом шудани геометрияро ошкоро зикр кардаанд.
• Усул ҳамчун равиши эҳтимолии тарҳрезии пешакӣ барои скан кардани соҳаи васеи параметрҳо пеш аз таҷрибаҳои гаронарзиш пешниҳод шудааст.

Маҳдудиятҳои таҳқиқот кадомҳоянд?

• Ин таҳқиқот preprint-и аз peer review нагузашта мебошад.
• ATAC-1 муодилаҳои пурраи Navier–Stokes-ро ҳал намекунад ва минтақаҳои ҷараёни ҷудошударо намояндагӣ карда наметавонад.
• Ҷадвали B-prime танҳо ҳамкориҳои Si–C–O-ро дар бар мегирад; нитроген ва дигар ҷузъҳои керамикиро берун мегузорад.
• Ҳамаи намудҳои пешгӯишуда дар динамикаи молекулавӣ ба модели миқёси макро интиқол дода нашудаанд.
• Эмиссивият ва ҷаббиши спектралӣ барои ҳамаи ҳароратҳо 0,87 собит қабул шудаанд.
• Муаллифон эътироф мекунанд, ки натиҷаҳои гузаронандагии гармӣ аз тамоюлҳои таҷрибавӣ фарқ мекунанд.
• Воҳид ё миқёспазии гармиғунҷоиши хос дар Шакли 6 ба қадри кофӣ равшан нест.
• Ҷадвали 1, ки дар матн ба он истинод шудааст, дар версияи preprint-и боршуда вуҷуд надорад.
• Баъзе натиҷаҳои ақибнишинии сатҳ аз геометрияи ибтидоӣ ва қиматҳои аз ҷиҳати физикӣ имконпазир зиёдтаранд.
• Байни Шакли 8 ва матни шарҳдиҳанда дар бораи он ки қимати аномалӣ ба силсилаи Mach 15 ё Mach 20 тааллуқ дорад, номувофиқӣ вуҷуд дорад.
• Натиҷаҳои модел бо CFD-и дақиқии баланд, нақби зарбаи гиперсадоӣ ё додаи воқеии парвоз тасдиқ нашудаанд.
• Паҳншавии номуайянӣ ва фосилаҳои миқдории хато пешниҳод нашудаанд.

Таҳқиқот чиро дастгирӣ мекунад?

Таҳқиқот нишон медиҳад, ки иттилооти гармӣ ва химиявии аз динамикаи реактивии молекулавӣ ҳосилшуда метавонад ба шакли B-prime табдил дода шуда, ба ҳалкунандаи аэротермалӣ интиқол дода шавад. Ҳамчунин нишон медиҳад, ки бо зиёд шудани ҳарорат гузаронандагии гармии ҳисобшудаи SiC кам, гармиғунҷоиши хос зиёд мешавад ва намудҳои бухоршавандаи сатҳӣ метавонанд мураккабтар гарданд. Модел тамоюле тавлид мекунад, ки дар шароити ҳарорати баланд ва Mach-и баланд нӯги бинӣ ва сатҳи бодхӯрда хавфи бештари аблятсия доранд.

Таҳқиқот чиро исбот намекунад?

Таҳқиқот исбот намекунад, ки қиматҳои ақибнишинии додашуда умри воқеии воситаи гиперсадоиро дуруст пешгӯӣ мекунанд. Симулятсияҳо нишон намедиҳанд, ки ягон воситаи парвозӣ бомуваффақият истеҳсол шудааст, системаи муҳофизати SiC ба давомнокии муайяни миссия тоб меорад ё дар шароити Mach 20 бехатар кор мекунад. Натиҷаҳои бузурги ғайрифизикӣ нишон медиҳанд, ки версияи ҷорӣ барои сертификатсия, ҳудуди бехатарӣ ё андозабандии ниҳоӣ истифода шуда наметавонад.

Аҳамияти эҳтимолӣ аз нигоҳи Туркия чист?

Равиши методологии таҳқиқот барои гурӯҳҳои таҳқиқотӣ дар Туркия, ки дар воситаҳои ҳавоии гиперсадоӣ, геометрияи бинии мушак ва ракета, системаҳои бозворидшавӣ, керамикаи ҳарорати баланд ва маводи муҳофизати гармӣ кор мекунанд, қобили интиқол аст. Хусусан пеш аз таҷрибаҳои гарони нақби зарба ё ark-jet муайян кардани намудҳои реаксия бо ҳисобҳои атомистӣ ва интиқоли ин иттилоот ба кодҳои зудкори тарҳрезӣ метавонад барои танг кардани фазои таҷриба кӯмак кунад.

Бо вуҷуди ин, натиҷаҳои худи таҳқиқот ҳам нишон медиҳанд, ки моделҳои зуд ва дақиқии паст бояд ҳатман бо CFD-и дақиқии баланд, таҷрибаҳои мавод ва диаграммаҳои мувофиқи фаза санҷида шаванд. Муҳимтарин дарсе, ки ба Туркия интиқол додан мумкин аст, қимати мушаххаси ақибнишинӣ нест; балки зарурати сохтани занҷири марҳилавии тасдиқи модел байни модели атомистӣ, ҳалкунандаи муҳандисӣ ва тасдиқи таҷрибавӣ мебошад.

Усул ва натиҷаҳои таҳқиқот

Натиҷаҳои асосии ададӣ

• Гузаронандагии гармии ҳисобшуда дар тақрибан 750 °C аз 48 W/m·°C то дар 2250 °C тақрибан 13 W/m·°C коҳиш ёфтааст.
• Каҷи гармиғунҷоиши хос дар Шакли 6 аз тақрибан 5,4 ба 6,15 баланд шудааст; аммо робитаи байни тамғаи воҳиди “J/kg·°C” дар график ва миқёси ададӣ шарҳ дода нашудааст.
• Дар поён аз 1500 °C намудҳои асосии бухоршаванда CO, O ва CO₂ пешгӯӣ шудаанд.
• Дар болотар аз 1500 °C гузориш шудааст, ки намудҳои C₂O₂, C₂O₃ ва SiO₃ метавонанд ташаккул ёбанд.
• Ҳарорати сублиматсияи барои SiC дар ATAC-1 истифодашуда 2850 °C мебошад.
• Дар Mach 20, 15.240 m ва кунҷи ҳамлаи 17° ҳарорати максималии сатҳ дар график ба тақрибан 3200 °C расидааст.
• Дар ҳамин минтақаи параметрӣ Шакли 8 ақибнишинии сатҳии ғайрифизикии тақрибан 3,10 × 10¹⁴ cm тавлид кардааст.
• Дар харитаҳои сатҳи Mach 15, 15.240 m ва 32 сония баландтарин ақибнишинӣ дар нӯги бинӣ ва минтақаи бодхӯрда ҷамъ шудааст.
• Дар Mach 5, кунҷи ҳамлаи 19° ва 900 сония ҳарорати сатҳ дар сатҳи тақрибан 900 °C монда, ақибнишинӣ дар миқёси ночиз ҳисоб шудааст.
• Mach 10 дар ҳамон баландӣ ва давомнокӣ тақрибан 101,6 cm ақибнишинӣ тавлид кардааст; муаллифон гуфтаанд, ки ин қимат аз геометрияи ибтидоӣ зиёд аст ва аз ин рӯ аз ҷиҳати физикӣ беэътибор мебошад.
• Дар Mach 7, кунҷи ҳамлаи 0° ва 20 сония фишори сатҳ 2,43 MPa ҳисоб шуда, сарфи назар аз пешгӯии остонаи ҳарорат, аблятсияи назаррас ба вуҷуд наомадааст.

Қадамҳои баъдӣ барои истифодаи модел дар муҳандисӣ

1. Васеъ кардани ҷадвали B-prime то нитроген ва дигар унсурҳои керамикаи ҳарорати ултрабаландро фаро гирад.
2. Таъриф кардани эмиссивият ва ҷаббиши спектралӣ ҳамчун тағйирёбандаҳои вобаста ба ҳарорат.
3. Такрор кардани ҳисобҳои MD бо ҳуҷайраҳои калонтар, вақтҳои дарозтари оромиш ва таҳлили номуайянӣ.
4. Интиқол додани ҷадвал ба ҳалкунандаҳои дақиқии баландтар, аз қабили ATAC-2, US3D, LAURA ё DPLR.
5. Тасдиқ кардани ҳарорати сатҳ ва қиматҳои ақибнишинӣ бо нақби зарба, ark-jet ё таҷрибаҳои назоратшавандаи аблятсия.
6. Татбиқи ҳудудҳои геометрӣ ва физикӣ, ки ҳангоми пурра тамом шудани мавод ҳисобро бозмедоранд.
7. Илова кардани диаграммаҳои боэътимоди фазаи SiC дар ҳарорат ва фишорҳои баланд ба модел.

Ёддошт оид ба манбаъ ва усул

Номи аслии таҳқиқот: Hypersonic Flight Aerothermodynamics and Ablation: A Theoretical Model

Муаллифон ва тартиб: Robert Slapikas, Luis Bravo, Anindya Ghoshal.

Муаллифи якуми баробар: Дар таҳқиқот маълумоти саҳми баробар ё муаллифии якуми баробар вуҷуд надорад.

Муаллифи тамос: Дар сабти SSRN Anindya Ghoshal ҳамчун муаллифи тамос нишон дода шудааст.

Пайвандҳои муассисавӣ: Robert Slapikas бо DEVCOM Army Research Laboratory ва SURVICE Engineering Company вобаста аст. Luis Bravo ва Anindya Ghoshal бо DEVCOM Army Research Laboratory вобастаанд.

Саҳми муаллифон: Robert Slapikas таҳқиқот, усул, нармафзор, визуализатсия, тасдиқ, тартиб додани дода, таҳлили формалӣ ва навиштани тарҳи аввалро анҷом додааст. Luis Bravo ба усул, навиштани тарҳи аввал, машварат ва таъмини захираҳо саҳм гузоштааст. Anindya Ghoshal ба усул, навиштан, машварат, захираҳо, идоракунии лоиҳа ва таъмини маблағгузорӣ саҳм гузоштааст.

DOI:10.2139/ssrn.6952638

Платформаи нашр: SSRN.

Ношири аслӣ ё оператори платформа: Таҳқиқот дар маҷаллаи peer-reviewed нашр нашудааст. SSRN платформаи preprint ва таҳқиқоти барвақт мебошад, ки аз ҷониби Elsevier идора мешавад.

Соли нашр: 2026.

Санаи боргузорӣ ба SSRN: 16 июни 2026.

Ҳолати peer review: Ин таҳқиқот preprint аст ва аз peer review нагузаштааст.

Навъи манбаъ: Preprint-и моделсозӣ/симулятсия, ки динамикаи молекулавӣ ва ҳисобҳои аэротермалиро дар бар мегирад.

Манбаи расмӣ:Саҳифаи расмии таҳқиқот дар SSRN

Маблағгузорӣ: Таҳқиқот аз ҷониби DEVCOM Army Research Laboratory дастгирӣ шуда, дар доираи созишномаи ҳамкорӣ бо рақами W15P7T-19-D-0126 иҷро шудааст.

Бархӯрди манфиат: Муаллифон гузориш додаанд, ки бархӯрди молиявии маълум ё муносибати шахсие, ки метавонад ба таҳқиқот таъсир расонад, вуҷуд надорад.

Ин шарҳи туркӣ бо баррасии таҳқиқоти боршудаи 36-саҳифагӣ, муодилаҳо, шароити симулятсия, графикҳо, харитаҳои сатҳ, натиҷаҳо ва маҳдудиятҳои зикркардаи муаллифон омода шудааст. Натиҷаҳои илмие, ки дар таҳқиқот вуҷуд надоранд, илова нашудаанд. Барои мундариҷаи илмӣ манбаи беруна истифода нашуда, санҷиши беруна танҳо бо маълумоти ҳувияти библиографӣ, аз қабили DOI, платформа, санаи боргузорӣ ва муаллифи тамос маҳдуд шудааст.

Дар версияи боршуда Ҷадвали 1, ки матн ба он истинод мекунад, вуҷуд надорад. Воҳид ё миқёси гармиғунҷоиши хос дар Шакли 6 ба қадри кофӣ шарҳ дода нашудааст. Илова бар ин, байни қимати фавқулода бузурги ақибнишинӣ дар Шакли 8 ва рақами Mach дар матни шарҳ номувофиқӣ вуҷуд дорад. Аз ин рӯ, таҳқиқот бояд на ҳамчун пешгӯии ниҳоии умри восита ё бехатарии мавод, балки ҳамчун таҳқиқоти ҳисоббарории марҳилаи барвақт арзёбӣ шавад, ки интиқоли додаи атомистиро ба ҳалкунандаи зудкори муҳандисӣ нишон медиҳад.

Ин таҳқиқот далели муваффақияти парвози воқеӣ, тасдиқи саҳроӣ ё сертификатсия пешниҳод намекунад. Натиҷаҳои миқдории модел набояд пеш аз тасдиқ бо CFD-и дақиқии баланд, таҷрибаҳои назоратшавандаи гиперсадоӣ ва таҳлили номуайянӣ мустақиман дар тарҳрезии воситаи воқеӣ истифода шаванд.