Ushbu tadqiqot gipertovushli parvoz vaqtida kremniy karbidning (SiC) qizishi, oksidlanishi va yuzadan material yo‘qotishi xatti-harakatini reaktiv molekulyar dinamika simulyatsiyalari bilan ATAC-1 aerotermik yechuvchisini birlashtirib o‘rganadi. Tadqiqotchilar LAMMPS va ReaxFF yordamida hisoblangan issiqlik o‘tkazuvchanligi, solishtirma issiqlik, sirt reaksiyalari va uchuvchan turlarni B-prime ablatsiya jadvaliga o‘tkazgan; bu jadvaldan Mach 5, 10, 15 va 20 sharoitlaridagi konussimon burun geometriyasining harorati va sirt chekinishini hisoblashda foydalangan. Simulyatsiyalar harorat oshgani sayin SiC ning issiqlik o‘tkazuvchanligi kamayishini, solishtirma issiqligi oshishini va taxminan 1500 °C dan yuqorida sirt kimyosi murakkabroq uchuvchan turlarni hosil qilishi mumkinligini ko‘rsatadi. Biroq tadqiqot peer-reviewdan o‘tmagan preprintdir; ishlatilgan ATAC-1 yechuvchisi to‘liq Navier–Stokes hisobini bajarmaydi va ayrim yuqori Mach sonlaridagi natijalar fizik geometriyadan oshib ketadigan, mualliflarning o‘zi ham yaroqsiz deb hisoblagan sirt yo‘qotish qiymatlarini hosil qiladi.

Tadqiqotning asosiy yangiligi atom darajasidagi kimyoviy va termomekanik ma’lumotni bevosita tezkor muhandislik loyihalash vositasiga uzatishga urinishdir. Shu tariqa faqat xona haroratidagi tajribalardan kelib chiqqan empirik jadvallar o‘rniga bog‘larning uzilishi va yangi molekulalar hosil bo‘lishi kabi jarayonlarni ifodalay oladigan molekulyar dinamika natijalaridan foydalanish maqsad qilingan. Biroq usul hali validatsiya yoki yakuniy loyihalash vositasi emas, balki tadqiqotchilar ta’biri bilan aytganda dastlabki bosqich va amalga oshirish mumkinligini ko‘rsatuvchi namoyishdir.

Tadqiqotning asosiy savoli nima?

Tadqiqotning asosiy savoli gipertovushli parvoz vaqtida kremniy karbid yuzasidan qaysi atom va molekulyar turlar ajralishini hamda bu kimyoviy jarayonlarni transport vositasining sirt harorati va material chekinishiga qanday o‘tkazish mumkinligini aniqlashdir. Tadqiqotchilar ayniqsa atom miqyosidagi oksidlanish va bog‘ uzilishi jarayonlari katta masshtabli aerotermik hisoblarda ishlatiladigan B-prime jadvallarini oziqlantira oladimi-yo‘qmi, shuni sinaydi.

Gipertovushli parvozda burun, qanot old qirrasi va boshqa tashqi yuzalar zarba to‘lqini, aerodinamik qarshilik, sirt ishqalanishi va yuqori haroratli chegara qatlami tufayli juda katta issiqlik yukiga duch keladi. Termik himoya tizimi bu energiyaning bir qismini o‘zida ushlab qolsa, bir qismi materialning parchalanishi, gazlanishi yoki yuzadan uzoqlashishi orqali sarflanadi. Yuzadan chiqadigan gazlarning chegara qatlamiga aralashib, devorga yetib keladigan issiqlik oqimini kamaytirishi tadqiqotda “to‘sish ta’siri” yoki “blocking effect” deb ataladi.

Mavjud modellarda qanday bo‘shliq nishonga olingan?

Keng tarqalgan B-prime modellari sirt kimyosi muvozanat holatiga yetadi deb faraz qiladigan jadvallarga tayanadi. Bu yondashuv hisoblash jihatidan tez; biroq zarba chegara qatlamidagi qisqa muddatli, metastabil va kinetik boshqariladigan reaksiyalarni to‘liq ifodalay olmasligi mumkin. Finite-rate kimyo modellari batafsilroq bo‘lsa-da, qimmat tajribalar va keng reaksiyalar mexanizmini talab qiladi. Tadqiqotda mavjud modellar orasida turli sirt turlarining hosil bo‘lishi bo‘yicha 300 foizdan ortiq farqlar qayd etilgani ta’kidlanadi.

Tadqiqotchilarning yechim taklifi kimyoviy reaksiyalarni klassik molekulyar dinamika darajasida hisoblash va shu natijalardan tezkor B-prime jadvali yaratishdir. Bu jadval yuqori aniqlikdagi oqim yechuvchisiga emas, balki dastlabki loyihalash va parametrik tekshiruvlar uchun ishlatiladigan ATAC-1 kodiga uzatilgan. Demak, tadqiqot atomistik kimyo bilan muhandislik masshtabidagi aerotermik yechim o‘rtasida zaif bog‘langan ma’lumot uzatish zanjirini yaratadi.

Zarba chegara qatlamida qaysi jarayonlar birgalikda sodir bo‘ladi?

Tadqiqotning 1-shaklida gipertovushli sirt yaqinidagi asosiy o‘zaro ta’sirlar birgalikda ko‘rsatilgan. Erkin oqim hosil qilgan zarbadan keyin gaz molekulalari dissotsiatsiyalanadi, ichki energiya darajalari qo‘zg‘aladi va sirt bilan kimyoviy o‘zaro ta’sirga kirishadi. Ayni hududda radiatsiya, materialdan gaz chiqishi, uglerodli turlarning diffuziyasi va zarba buzilishlari sodir bo‘ladi. Bu jarayonlarning barchasi bir xil vaqt va uzunlik masshtabida emas: bog‘ uzilishi femtosekund yoki nanosekund darajasida sodir bo‘lsa, sirt yeyilishi va geometrik o‘zgarish soniyalar davomida yig‘ilishi mumkin.

Atomistik model qanday qurilgan?

Kremniy karbidning termomekanik va oksidlanish xatti-harakati Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) dasturi bilan bajarilgan klassik molekulyar dinamika simulyatsiyalari orqali modellashtirilgan. Kimyoviy bog‘larning uzilishi, yangi bog‘larning hosil bo‘lishi va atomlar orasida zaryad uzatilishiga imkon berish uchun reaktiv kuch maydoni ReaxFF va zaryad muvozanatlash modeli QEq ishlatilgan.

ReaxFF tizimining umumiy potensial energiyasi tadqiqotda quyidagi komponentlar bilan ifodalanadi:

\[ 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: Bog‘langan atomlarning bog‘ energiyasi hissasi.
• U_over: Atomlarning fizik bo‘lmagan ortiqcha koordinatsiyasini cheklaydi.
• U_under: Yetarlicha koordinatsiyalanmagan atomlarning barqarorligiga hissa qo‘shadi.
• U_valence: Bog‘ burchaklariga oid energiya hadi.
• U_lp: Yolg‘iz elektron jufti energiyasi.
• U_torsion: To‘rt atomli buralish burchaklaridan kelib chiqadigan energiya.
• U_vdW: Qisqa masofali van der Waals o‘zaro ta’sirlarini ifodalaydi.
• U_Coul: Uzoq masofali Coulomb o‘zaro ta’siri.

Bu umumiy energiya molekulyar dinamikada har bir atomga ta’sir qiluvchi kuchlar va vaqt bo‘yicha harakatni hisoblashda ishlatiladi. Modelning ahamiyati shundaki, u bog‘ tuzilishi sobit bo‘lgan an’anaviy potensiallardan farqli ravishda oksidlanish vaqtida yangi Si–O va C–O bog‘larining hosil bo‘lishiga imkon beradi.

Issiqlik o‘tkazuvchanligi qanday hisoblangan?

Polikristall va monokristall α(6H)-SiC tuzilmalarining issiqlik o‘tkazuvchanligi Green–Kubo munosabati bilan hisoblangan:

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

Bu yerda k_ij issiqlik o‘tkazuvchanlik tensorini, V simulyatsiya hujayrasining hajmini, k_B Boltzmann doimiysini, T mutlaq haroratni va J issiqlik oqimini ifodalaydi. Burchak qavs ichidagi had boshlang‘ich issiqlik oqimi bilan keyingi vaqtdagi issiqlik oqimi orasidagi avtokorrelyatsiyadir. Korrelyatsiyaning uzoq vaqt saqlanishi energiya tuzilma ichida samaraliroq tashilishini anglatadi.

Termomekanik hisoblarda taxminan 20 × 20 × 50 nm o‘lchamdagi simulyatsiya hujayrasi ishlatilgan. X va Y yo‘nalishlarida davriy chegara shartlari qo‘llanib, issiqlik energiyasi 2-shaklda ko‘rsatilganidek issiq hududdan sovuq hududga uzatilgan. Tasvirdagi turli ranglar polikristall tuzilmadagi alohida donalarni ko‘rsatadi.

Solishtirma issiqlik va bosim hisoblari

Solishtirma issiqlik sig‘imi uchun LAMMPS ichidagi garmonik moslashtirilgan o‘rtachalash usuli (harmonically mapped averaging, HMA) ishlatilgan. Maqoladagi energiya, bosim va doimiy hajmdagi solishtirma issiqlik ifodalari quyidagicha:

\[ \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 \]

Bu tenglamalarda d tizim o‘lchamlari sonini, N atomlar sonini, U umumiy energiyani, F atom kuchlarini, Δr atomlarning panjara holatlaridan siljishini, Φ Hessian matritsasini, ρ son zichligini va P_vir virial bosimni ko‘rsatadi. Energiya tebranishlarining kattaligi doimiy hajmdagi solishtirma issiqlik hisobiga bevosita hissa qo‘shadi.

Oksidlanish va ablatsiya simulyatsiyasi qanday bajarilgan?

Oksidlanish modeli uchun gorizontal o‘lchami 20 × 20 nm va vertikal o‘lchami 15 nm bo‘lgan alohida hujayra qurilgan. Hujayrada 10 nm qalinlikdagi SiC qatlami va uning ustida chegara qatlamini ifodalovchi 5 nm bo‘shliq mavjud. Molekulyar va atomar kislorod turli tezlik va haroratlarda SiC yuzasiga yo‘naltirilgan.

3-shaklda pastki qismdagi qizil atomlar qotirilgan, yashil atomlar harorati termostat bilan boshqarilgan va kislorod atomlari tizimga yuqori qismdan qo‘shilgan. Hujayra ustida yuzadan yetarlicha uzoqlashgan turlar simulyatsiyadan o‘chirib tashlanadigan hudud ta’riflangan. Shu tariqa kislorodning yuzaga urilishi, kimyoviy bog‘ hosil qilishi va ayrim mahsulotlarning yuzadan ajralishi atom darajasida kuzatilgan.

Langevin termostatidagi umumiy kuch quyidagicha ajratilgan:

\[ \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 potensialidan keladigan konservativ kuchni, F_f ishqalanish kuchini va F_r tasodifiy kuch hadini ifodalaydi. m atom massasi, V tezlik, τ_damp harorat relaksatsiyasi vaqti va dt vaqt qadamidir. Tasodifiy kuch tenglamasi bu yerda maqolada yozilgan shakli bilan aynan uzatilgan.

Molekulyar dinamika qaysi uchuvchan turlarni bashorat qilgan?

Simulyatsiyalarda 1500 °C dan pastda SiC yuzasidan ajraladigan asosiy turlar karbon monoksid (CO), atomar kislorod (O) va karbonat angidrid (CO₂) ekani bildirilgan. Yuqoriroq haroratlarda C₂O₂, C₂O₃ va SiO₃ kabi murakkabroq va metastabil turlar hosil bo‘lishi mumkinligi bashorat qilingan.

Bunga qarshi, ATAC-1 ga uzatilgan dastlabki B-prime jadvali hisoblash yukini kamaytirish maqsadida faqat CO, SiO₂ va CO₂ turlari bilan cheklangan. Demak, molekulyar dinamikada kuzatilgan barcha reaksiya mahsulotlari makromasshtabli sirt chekinishi hisobiga ko‘chirilmagan. Azot bilan va ultra yuqori haroratli keramika tarkibida ishlatiladigan o‘tish metallari bilan yuz berishi mumkin bo‘lgan reaksiyalar ham joriy versiyada mavjud emas.

ATAC-1 qanday ishlaydi?

ATAC-1 tadqiqotda ba’zan CFD kodi deb atalgan bo‘lsa-da, to‘liq Navier–Stokes tenglamalarini yechmaydi. Kod zarbadan keyingi oqimdan boshlanadigan inviscid oqim chiziqlarini kuzatib, bosim va tezlik taqsimotini chiqaradi; keyin chegara qatlami yaqinlashuvlari, yarim empirik korrelyatsiyalar va chekli ayirmalar bilan ishlaydigan termik yechuvchidan foydalanadi. Sirt harorati va material chekinish tezligi har bir geometrik nuqtadagi nochiziqli energiya balansi Newton–Raphson usuli bilan yechilib aniqlanadi.

Bu yondashuv yuqori aniqlikdagi CFD yechuvchilaridan tezroq; biroq viskoz oqim tafsilotlarini, turbulentlikni, murakkab zarba o‘zaro ta’sirlarini va ajralgan oqim hududlarini bevosita yecha olmaydi. Tadqiqotning maqsadi ham yakuniy aniqlikdagi parvoz bashoratini berish emas, balki yangi B-prime jadvali tezkor muhandislik vositasiga integratsiya qilinishi mumkinligini ko‘rsatishdir.

Gipertovushli burun konusi uchun qaysi sharoitlar o‘rganilgan?

ATAC-1 hisoblarida 1976 AQSh Standart Atmosferasi va havo uchun Clark Mollier jadvali ishlatilgan. Yaw va roll burchaklari nol deb qabul qilingan. Ikki asosiy balandlik o‘rganilgan:

• 15.240 m: 11.664,53 Pa atmosfera bosimi va −56,5 °C harorat.
• 21.336 m: 4.487,68 Pa atmosfera bosimi va −55,23 °C harorat.

Geometriya 7,5° samarali yarim konus burchagiga ega 10 cm masshtabli simmetrik burun konusidir. Vaqt qadami 0,01 soniya deb belgilangan, sirtda P1–P11 deb nomlangan 11 batafsil tahlil nuqtasi tanlangan. 4-shakl bu nuqtalar burun uchida va konusning shamolga qaragan hamda shamol orti yuzalarida qanday taqsimlanganini ko‘rsatadi.

Asosiy parametrik qator Mach 5, 10, 15 va 20 tezliklarida bajarilgan. Hujum burchagi 5° bilan birga 10°–20° oralig‘ida bir darajalik qadamlar bilan o‘zgartirilgan va asosiy hisoblarning parvoz davomiyligi 32 soniya deb olingan. Bundan tashqari, Mach 5 va Mach 10 uchun 19° hujum burchagida 900 soniyalik uzoq davomli ssenariylar hisoblangan.

Turg‘unlik harorati qanday baholangan?

Tadqiqotda turli balandlik va Mach sonlari uchun soddalashtirilgan turg‘unlik harorati munosabati ham ishlatilgan:

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

T₀ turg‘unlik haroratini, T_∞ tegishli balandlikdagi erkin oqim haroratini, M Mach sonini, Pr Prandtl sonini va γ solishtirma issiqlik nisbatini ifodalaydi. Hisobda havo ideal gaz deb qabul qilingan va γ qiymati 1,4 da sobit tutilgan. Bu sodda munosabat tezlikning kinetik energiyasi oqim to‘xtatilganda qay darajada issiqlik energiyasiga aylanishi mumkinligini ko‘rsatadi; biroq haqiqiy sirt reaksiyalari, yuqori bosimdagi faza xatti-harakati va muvozanatdan tashqari kimyoni o‘zi yolg‘iz tushuntira olmaydi.

Termik xususiyat natijalari nimani ko‘rsatgan?

5-shaklda MD bilan hisoblangan α-SiC issiqlik o‘tkazuvchanligi taxminan 750 °C da 48 W/m·°C darajadan 2250 °C atrofida taxminan 13 W/m·°C gacha pasayadi. Pasayish dastlabki harorat oralig‘ida tezroq, taxminan 1600 °C dan keyin esa tekisroq. Tadqiqotchilar 1250 °C va 1500 °C yaqinidagi kichik tebranishlarni oksidlanish yoki faza o‘tishi chegaralari bilan bog‘laydi, biroq natijalarga simulyatsiya o‘lchami, chegara shartlari, termostat relaksatsiyasi va e’tibordan chetda qolgan kvant ta’sirlari ta’sir qilishi mumkinligini tan oladi.

6-shakldagi solishtirma issiqlik qiymatlari taxminan 750 °C dan 2250 °C gacha muntazam oshib, grafikda taxminan 5,4 dan 6,15 gacha ko‘tariladi. Biroq grafikning vertikal o‘qi “J/kg·°C” deb belgilangan bo‘lsa-da, berilgan sonli masshtab bu birlik bilan ochiq ravishda tushuntirilmagan. Shu sabab grafikdagi qiymatlarni mutlaq muhandislik ma’lumoti sifatida bevosita ishlatishdan oldin birlik va masshtab validatsiyasi zarur.

Tezlik, balandlik va hujum burchagi haroratga qanday ta’sir qilgan?

7-shaklda umumiy tendensiya Mach soni oshgani sari maksimal sirt harorati ko‘tarilishidir. 15.240 m balandlikdagi Mach 20 qatori ko‘pchilik hujum burchaklarida SiC uchun modelda ishlatilgan 2850 °C sublimatsiya chegarasiga yaqinlashadi yoki undan oshadi. Eng aniq oshish 17° hujum burchagida kuzatiladi va grafikda maksimal harorat taxminan 3200 °C gacha chiqadi.

Yuqoriroq 21.336 m balandlikda atmosfera zichligi past bo‘lgani sabab Mach 5 haroratlari 15.240 m natijalaridan pastroq. Biroq Mach 15 va Mach 20 kabi yuqori tezliklarda haroratlarning yana 2600–2850 °C oralig‘iga yetgani ko‘rinadi. Shu sabab balandlik ta’sirini alohida baholab bo‘lmaydi; u tezlik, sirt bosimi, hujum burchagi va ishlatilgan oqim yaqinlashuvi bilan birga ko‘rib chiqilishi kerak.

Material chekinishi natijalari nega muammoli?

8-shakl fizik tendensiyalar bilan birga modelning amal qilish chegarasidan chiqadigan natijalarni ham hosil qilgan. Mach 20, 15.240 m va 17° hujum burchagi qatorida sirt chekinishi taxminan 3,10 × 10¹⁴ cm deb ko‘rsatilgan. Bu qiymat hech qanday haqiqiy transport vositasi geometriyasiga mos kelmaydi. Matnning bir qismida 17°–18° anomaliyasi Mach 15 ga tegishli deb berilgan bo‘lsa, grafik legendasi favqulodda qiymatni Mach 20 va 15 km qatoriga joylashtiradi. Demak, natija bayoni bilan grafik o‘rtasida Mach raqami bo‘yicha nomuvofiqlik mavjud.

Tadqiqotchilar bu anomaliya harorat va bosim B-prime jadvali kalibrlangan soha tashqarisiga chiqqanda ATAC-1 material xatti-harakatini eksponensial tarzda ekstrapolyatsiya qilishidan kelib chiqishi mumkinligini tushuntiradi. Sirt harorati 2850 °C sublimatsiya chegarasidan oshganda, model qattiq material butunlay yo‘qolgan bo‘lsa ham chekinish hisobini davom ettiradi. Boshlang‘ich material qalinligi sarflanib bo‘lgandan keyingi vaqt qadamlarining fizik ma’nosi yo‘q.

Sirtda eng katta yeyilish qayerda hosil bo‘lgan?

9-shakldagi Mach 15 sirt xaritalari material yo‘qotishining katta qismi burun uchida va parvoz yo‘nalishiga qaragan shamol oldi tomonida to‘planganini ko‘rsatadi. Hujum burchagi o‘zgargani sari yuqori chekinish hududi simmetrik burun markazidan yon tomonga siljiydi. Bu natija turli hujum burchaklarida bir xil umumiy harorat hosil bo‘lsa ham, sirt bo‘ylab mahalliy yuk taqsimoti o‘zgarishi mumkinligini ko‘rsatadi.

Biroq xaritaning rang shkalasi 32 soniya ichida 80 cm gacha umumiy chekinishni ko‘rsatadi. O‘rganilgan burun geometriyasi 10 cm masshtabda ekanini hisobga olsak, bu qiymatlar ham haqiqiy material hisoblash muddati tugashidan oldin tugab ketishini bildiradi. Bu xaritalarni miqdoriy xizmat muddati bashorati emas, balki model sirtning qaysi hududlarida yuqori yukni bashorat qilishini ko‘rsatadigan tendensiya xaritalari sifatida o‘qish kerak.

Mach 5 va Mach 10 uzoq davomli ssenariylari

15.240 m balandlik, 19° hujum burchagi va 900 soniya parvoz davomiyligida Mach 5 bilan Mach 10 taqqoslangan. Mach 5 ssenariysida sirt harorati taxminan 900 °C atrofida bo‘lgan va tadqiqotda SiC uchun ko‘rsatilgan taxminan 1400 °C chekinish boshlanish darajasidan past qolgan. 10-shakldagi sirt yo‘qotishi juda kichik sonli masshtab ichida qoladi.

Mach 10 ssenariysida esa harorat 2500 °C dan yuqorida tez tebranishlar ko‘rsatgan va hisoblangan jami chekinish taxminan 101,6 cm ga yetgan. Bu qiymat boshlang‘ich geometriyadan oshadi. Mualliflar tez harorat tebranishlari vaqt qadami aniqligi va ablatsiya hisobi bilan bog‘liq cheklovlarni aks ettirishi mumkinligini; geometriya butunlay yo‘qolgandan keyingi hisoblar yaroqsizligini ta’kidlaydi.

Mach 7 misoli nimani ko‘rsatgan?

Soddalashtirilgan turg‘unlik harorati hisobi 15.240 m balandlikda Mach 7 tezlikda SiC yuzasi ablatsiya boshlanish haroratidan oshishi mumkinligini ko‘rsatgan. Bunga qarshi, ATAC-1 bilan 0° hujum burchagi va 20 soniya uchun bajarilgan hisobda sirt bosimi 2,43 MPa ga yetgan va faqat yaxlitlash xatosi darajasida material yo‘qotilishi yuz bergan.

Bu nomuvofiqlik ablatsiyani faqat harorat chegarasi bilan aniqlab bo‘lmasligini ko‘rsatadi. Yuqori bosimdagi erish, bug‘lanish va sublimatsiya chegaralarini ishonchli faza diagrammalari bilan ta’riflash kerak. Tadqiqotda ishlatilgan sodda ideal gaz va sobit xususiyat farazlari bu o‘tishlarni barcha tafsilotlari bilan ifodalamaydi.

Tadqiqotning kuchli tomonlari nimalar?

• Atomistik reaksiyalarni tezkor aerotermik loyihalash vositasiga uzatadigan aniq ko‘p masshtabli ish oqimini taqdim etadi.
• Bog‘ uzilishi, bog‘ hosil bo‘lishi va zaryad uzatilishini ifodalay oladigan ReaxFF ishlatilgan.
• Issiqlik o‘tkazuvchanligi, solishtirma issiqlik, sirt kimyosi, uchuvchan turlar, harorat va geometrik chekinish bir xil model zanjirida ko‘rib chiqilgan.
• Turli Mach sonlari, ikki balandlik, ko‘plab hujum burchaklari va ikki xil parvoz davomiyligi taqqoslangan.
• Mualliflar modelning fizik bo‘lmagan natijalarini yashirmagan; ekstrapolyatsiya va geometriyaning tugab qolishi muammolarini ochiq ko‘rsatgan.
• Usul yuqori xarajatli tajribalardan oldin keng parametr sohasini ko‘zdan kechirish uchun potensial dastlabki loyihalash yondashuvini taklif qiladi.

Tadqiqotning cheklovlari nimalar?

• Ushbu tadqiqot peer-reviewdan o‘tmagan preprintdir.
• ATAC-1 to‘liq Navier–Stokes tenglamalarini yechmaydi va ajralgan oqim hududlarini ifodalay olmaydi.
• B-prime jadvali faqat Si–C–O o‘zaro ta’sirlarini o‘z ichiga oladi; azot va boshqa keramika komponentlarini tashqarida qoldiradi.
• Molekulyar dinamikada bashorat qilingan turlarning barchasi makromasshtabli modelga uzatilmagan.
• Emissivlik va spektral yutilish barcha haroratlarda sobit 0,87 deb qabul qilingan.
• Issiqlik o‘tkazuvchanligi natijalarida tajribaviy tendensiyalardan og‘ishlar borligini mualliflar tan oladi.
• 6-shakldagi solishtirma issiqlik birligi yoki masshtablanishi yetarlicha aniq emas.
• Matnda havola qilingan 1-jadval yuklangan preprint versiyasida mavjud emas.
• Ayrim sirt chekinishi natijalari boshlang‘ich geometriya va fizik jihatdan mumkin bo‘lgan qiymatlardan oshib ketadi.
• 8-shakl bilan izoh matni o‘rtasida anomal qiymat Mach 15 yoki Mach 20 qatoriga tegishlimi, degan masalada nomuvofiqlik mavjud.
• Model natijalari yuqori aniqlikdagi CFD, gipertovushli zarba tunneli yoki haqiqiy parvoz ma’lumoti bilan validatsiya qilinmagan.
• Noaniqlik tarqalishi va miqdoriy xato oraliqlari taqdim etilmagan.

Tadqiqot nimani qo‘llab-quvvatlaydi?

Tadqiqot reaktiv molekulyar dinamikadan olingan termik va kimyoviy ma’lumot B-prime shakliga aylantirilib, aerotermik yechuvchiga uzatilishi mumkinligini ko‘rsatadi. Shuningdek, harorat oshgani sari SiC ning hisoblangan issiqlik o‘tkazuvchanligi kamayishini, solishtirma issiqligi oshishini va uchuvchan sirt turlari murakkablashishi mumkinligini ko‘rsatadi. Model yuqori harorat va yuqori Mach sharoitlarida burun uchi va shamol oldi yuzasi yuqoriroq ablatsiya xavfiga ega ekanini ko‘rsatuvchi tendensiya hosil qiladi.

Tadqiqot nimani isbotlamaydi?

Tadqiqot berilgan chekinish qiymatlari haqiqiy gipertovushli transport vositasining xizmat muddatini to‘g‘ri bashorat qilishini isbotlamaydi. Simulyatsiyalar biror parvoz vositasi muvaffaqiyatli ishlab chiqarilganini, SiC himoya tizimi ma’lum missiya davomiyligiga bardosh berishini yoki Mach 20 sharoitida xavfsiz ishlashini ko‘rsatmaydi. Fizik bo‘lmagan katta natijalar joriy versiyadan sertifikatsiya, xavfsizlik chegarasi yoki yakuniy o‘lchamlash maqsadida foydalanib bo‘lmasligini ko‘rsatadi.

Turkiya nuqtai nazaridan mumkin bo‘lgan ahamiyati nima?

Tadqiqotning metodologik yondashuvi Turkiyada gipertovushli havo vositalari, raketa va raketa burun geometriyalari, qayta kirish tizimlari, yuqori haroratli keramika va termik himoya materiallari ustida ishlaydigan tadqiqot guruhlari uchun ko‘chirilishi mumkin. Ayniqsa qimmat zarba tunneli yoki ark-jet tajribalaridan oldin atomistik hisoblar bilan reaksiya turlarini aniqlash va bu ma’lumotni tezkor loyihalash kodlariga uzatish tajriba sohasini toraytirishga yordam berishi mumkin.

Biroq tadqiqotning o‘z natijalari tezkor va past aniqlikdagi modellar albatta yuqori aniqlikdagi CFD, material tajribalari va mos faza diagrammalari bilan sinovdan o‘tkazilishi kerakligini ham ko‘rsatadi. Turkiyaga ko‘chirilishi mumkin bo‘lgan eng kuchli saboq ma’lum bir chekinish qiymati emas; atomistik model, muhandislik yechuvchisi va tajribaviy validatsiya o‘rtasida bosqichma-bosqich model validatsiyasi zanjiri qurilishi zarurligidir.

Tadqiqot usuli va natijalari

Asosiy sonli natijalar

• Hisoblangan issiqlik o‘tkazuvchanligi taxminan 750 °C da 48 W/m·°C dan 2250 °C da taxminan 13 W/m·°C gacha pasaygan.
• 6-shakldagi solishtirma issiqlik egri chizig‘i taxminan 5,4 dan 6,15 gacha ko‘tarilgan; biroq grafikdagi “J/kg·°C” birlik yorlig‘i bilan sonli masshtab o‘rtasidagi bog‘lanish tushuntirilmagan.
• 1500 °C dan pastda asosiy uchuvchan turlar CO, O va CO₂ deb bashorat qilingan.
• 1500 °C dan yuqorida C₂O₂, C₂O₃ va SiO₃ turlari hosil bo‘lishi mumkinligi bildirilgan.
• SiC uchun ATAC-1 da ishlatilgan sublimatsiya harorati 2850 °C.
• Mach 20, 15.240 m va 17° hujum burchagida maksimal sirt harorati grafikda taxminan 3200 °C ga yetgan.
• Ayni parametr hududida 8-shakl taxminan 3,10 × 10¹⁴ cm kabi fizik bo‘lmagan sirt chekinishini hosil qilgan.
• Mach 15, 15.240 m va 32 soniyalik sirt xaritalarida eng katta chekinish burun uchi va shamol oldi hududida to‘plangan.
• Mach 5, 19° hujum burchagi va 900 soniyada sirt harorati taxminan 900 °C darajada qolgan va chekinish e’tiborsiz darajada hisoblangan.
• Mach 10, ayni balandlik va davomiylikda taxminan 101,6 cm chekinish hosil qilgan; mualliflar bu qiymat boshlang‘ich geometriyadan oshgani uchun fizik jihatdan yaroqsiz ekanini ta’kidlagan.
• Mach 7, 0° hujum burchagi va 20 soniyada sirt bosimi 2,43 MPa deb hisoblangan va harorat chegarasi bashoratiga qaramay sezilarli ablatsiya yuz bermagan.

Modelni muhandislikda ishlatish uchun zarur keyingi qadamlar

1. B-prime jadvalini azot va ultra yuqori haroratli keramikadagi boshqa elementlarni ham qamrab oladigan tarzda kengaytirish.
2. Emissivlik va spektral yutilishni haroratga bog‘liq o‘zgaruvchilar sifatida ta’riflash.
3. MD hisoblarini kattaroq hujayralar, uzoqroq relaksatsiya vaqtlar va noaniqlik tahlillari bilan takrorlash.
4. Jadvalni ATAC-2, US3D, LAURA yoki DPLR kabi yuqori aniqlikdagi yechuvchilarga uzatish.
5. Zarba tunneli, ark-jet yoki nazoratli ablatsiya tajribalari bilan sirt harorati va chekinish qiymatlarini validatsiya qilish.
6. Material butunlay tugaganda hisobni to‘xtatadigan geometrik va fizik chegaralarni qo‘llash.
7. Yuqori harorat va bosimlarda ishonchli SiC faza diagrammalarini modelga qo‘shish.

Manba va usul bo‘yicha izoh

Tadqiqotning asl nomi: Hypersonic Flight Aerothermodynamics and Ablation: A Theoretical Model

Mualliflar va tartib: Robert Slapikas, Luis Bravo, Anindya Ghoshal.

Teng birinchi muallif: Tadqiqotda teng hissa yoki teng birinchi mualliflik ma’lumoti mavjud emas.

Aloqa muallifi: SSRN yozuvida Anindya Ghoshal aloqa muallifi sifatida ko‘rsatilgan.

Institutsional aloqalar: Robert Slapikas DEVCOM Army Research Laboratory va SURVICE Engineering Company bilan bog‘langan. Luis Bravo va Anindya Ghoshal DEVCOM Army Research Laboratory bilan bog‘langan.

Muallif hissalari: Robert Slapikas tadqiqot, usul, dasturiy ta’minot, vizualizatsiya, validatsiya, ma’lumotlarni tartibga solish, formal tahlil va dastlabki qoralamani yozishni bajargan. Luis Bravo usul, dastlabki qoralamani yozish, maslahat va resurslar ta’minotiga hissa qo‘shgan. Anindya Ghoshal usul, yozish, maslahat, resurslar, loyiha boshqaruvi va moliyalashtirishni ta’minlashga hissa qo‘shgan.

DOI:10.2139/ssrn.6952638

Nashr platformasi: SSRN.

Asl nashriyot yoki platforma operatori: Tadqiqot peer-reviewli jurnalda nashr etilmagan. SSRN Elsevier tomonidan boshqariladigan preprint va dastlabki tadqiqot platformasidir.

Nashr yili: 2026.

SSRN ga yuklangan sana: 16 iyun 2026.

Peer-review holati: Ushbu tadqiqot preprint bo‘lib, peer-reviewdan o‘tmagan.

Manba turi: Molekulyar dinamika va aerotermik hisoblarni o‘z ichiga olgan modellashtirish/simulyatsiya preprinti.

Rasmiy manba:SSRN rasmiy tadqiqot sahifasi

Moliyalashtirish: Tadqiqot DEVCOM Army Research Laboratory tomonidan qo‘llab-quvvatlangan va W15P7T-19-D-0126 raqamli hamkorlik kelishuvi doirasida bajarilgan.

Manfaatlar to‘qnashuvi: Mualliflar ma’lum moliyaviy manfaatlar to‘qnashuvi yoki tadqiqotga ta’sir qilishi mumkin bo‘lgan shaxsiy munosabat mavjud emasligini bildirgan.

Ushbu turkcha izoh yuklangan 36 sahifalik tadqiqot; uning tenglamalari, simulyatsiya sharoitlari, grafiklari, sirt xaritalari, natijalari va mualliflar ko‘rsatgan cheklovlar o‘rganilib tayyorlangan. Tadqiqotda bo‘lmagan ilmiy natijalar qo‘shilmagan. Ilmiy mazmun uchun tashqi manba ishlatilmagan; tashqi tekshiruv faqat DOI, platforma, yuklangan sana va aloqa muallifi kabi bibliografik identifikatsiya ma’lumotlari bilan cheklangan.

Yuklangan versiyada matn havola qilgan 1-jadval mavjud emas. 6-shakldagi solishtirma issiqlik birligi yoki masshtabi yetarlicha tushuntirilmagan. Bundan tashqari, 8-shakldagi nihoyatda katta chekinish qiymati bilan izoh matnidagi Mach raqami o‘rtasida nomuvofiqlik mavjud. Shu sabablarga ko‘ra tadqiqot yakuniy transport vositasi xizmat muddati yoki material xavfsizligi bashorati sifatida emas, atomistik ma’lumotni tezkor muhandislik yechuvchisiga uzatishni ko‘rsatuvchi dastlabki bosqichdagi hisoblash tadqiqoti sifatida baholanishi kerak.

Ushbu tadqiqot haqiqiy parvoz muvaffaqiyati, dala validatsiyasi yoki sertifikatsiya dalilini taqdim etmaydi. Modelning miqdoriy natijalari yuqori aniqlikdagi CFD, nazoratli gipertovushli tajribalar va noaniqlik tahlili bilan validatsiya qilinmasdan haqiqiy transport vositasi loyihasida bevosita ishlatilmasligi kerak.