Bu tədqiqat, hipersonik uçuş sırasında silisium karbidün (SiC) isinmə, oksidləşmə və səthdən material kaybetme davranışını; reaktiv molekulyar dinamika simulyasiyaları ilə ATAC-1 aerotermal həll edicisini birleştirerek incelemektedir. Tədqiqatçılar, LAMMPS və ReaxFF kullanılarak hesaplanan istilik keçiriciliyi, xüsusi istilik, səth reaksiyonları və uçucu növləri bir B-prime ablasiya tablosuna aktarmış; bu tabloyu Mach 5, 10, 15 və 20 koşullarındaki konik bir burun geometrisinin temperatur və səth geriləməsini hesablamak üçün kullanmıştır. Simülasyonlar, temperatur arttıkça SiC’nin ısıl ilətkenliğinin azaldığını, xüsusi istiliksının arttığını və təxminən 1500 °C üzərindəki səth kimyasının daha karmaşık uçucu növlər üretebildiğini göstərir. Bununla birlikte tədqiqat hakem qiymətləndirməsindən keçməmiş bir preprinttir; istifadə olunan ATAC-1 həll edicisi tam Navier–Stokes hesablaması yapmamakta və bazı yüksək Mach sayı nəticələrı fiziki geometriyi aşan, yazarların da geçersiz kabul ettiği səth kaybı değerleri üretmektedir.

Tədqiqatnın əsas yeniliği, atom düzeyindeki kimyəvi və termomekanik bilgiyi doğrudan hızlı bir mühəndislik tasarım aracına aktarmaya tədqiqatsıdır. Böylece yalnızca oda sıcaklığı deneylerinden türetilmiş ampirik tablolar yerine, bşəbəkə qırılması və yeni molekül oluşumu gibi süreçleri temsil edebilən molekulyar dinamika nəticələrının kullanılması hedeflenmiştir. Ancak metod henüz bir doğrulama və ya nihai tasarım aracı değil, tədqiqatçıların ifadesiyle bir ilk aşama və uygulanabilirlik gösterimidir.

Tədqiqatnın əsas sorusu nedir?

Tədqiqatnın əsas sorusu, hipersonik uçuş sırasında bir silisium karbid səthinden hangi atomik və molekulyar növlərin ayrılacşəbəkəını və bu kimyəvi süreçlerin aracın səth sıcaklığı ilə material geriləməsine nasıl aktarılabiləceğini belirlemektir. Tədqiqatçılar özellikle atom ölçeğindeki oksidasyon və bşəbəkə kırılma süreçlerinin, daha büyük ölçekli aerotermal hesablamalarda istifadə olunan B-prime tablolarını besleyip besleyemeyeceğini sınamaktadır.

Hipersonik uçuşta burun, kanat ön kənarı və diğer dış səthler; zərbə dalğası, aerodinamik sürükleme, səth sürtünməsi və yüksək temperaturlı sərhəd qatı nedeniyle yoğun bir ısı yüküyle karşılaşır. Termal koruma sistemi bu enerjinin bir bölümünü bünyesinde tutarken bir bölümü de materialnin bozunması, gazlaşması və ya səthdən uzaklaşması yoluyla harcanır. Yüzeyden çıkan gazların sərhəd qatına karışarak divara ulaşan istilik axınını azaltmasına tədqiqatda “engelleme etkisi” və ya “blocking effect” adı vərilmektedir.

Mevcut modellerde hangi boşluk hedeflenmiştir?

Yaygın B-prime modelleri, səth kimyasının denge durumuna ulaştığını varsayan tablolara dayanır. Bu yaklaşım hesablama açısından hızlıdır; ancak şok sərhəd qatındaki kısa süreli, metastabil və kinetik kimi kontrol edilən reaksiyonları bütünüyle temsil etmeyebilir. Finite-rate kimya modelleri daha ayrıntılı olmakla birlikte pahalı deneyler və kapsamlı reaksiyon mekanizmaları gerektirir. Tədqiqatda, mevcut modeller arasında fərqli səth növlərinin oluşumu bakımından yüzde 300’ü aşabilən farklar bildirildiği vurgulanmaktadır.

Tədqiqatçıların çözüm önerisi, kimyəvi reaksiyonları klasik molekulyar dinamika düzeyinde hesablamak və bu nəticələrdan hızlı bir B-prime tablosu üretmektir. Bu tablo, yüksək doğruluklu bir axın həll edicisine değil, ön tasarım və parametrik incelemeler üçün istifadə olunan ATAC-1 koduna aktarılmıştır. Dolayısıyla tədqiqat, atomistik kimya ilə mühəndislik ölçekli aerotermal çözüm arasında zayıf bşəbəkəlı bir vəri aktarım zinciri kurmaktadır.

Şok sərhəd qatında hangi süreçler birlikte gerçekleşmektedir?

Tədqiqatnın Şekil 1’inde hipersonik bir səthin yakınındaki başlıca etkiləşimler birlikte gösterilmektedir. Serbest axının oluşturduğu şokun ardından gaz molekülleri ayrışmakta, iç enerji düzeyleri uyarılmakta və səthle kimyəvi etkiləşime girmektedir. Aynı bölgede radyasyon, materialden gaz çıkışı, karbon içeren növlərin difüzyonu və şok bozulmaları gerçekleşmektedir. Bu süreçlerin tamamı eyni zaman və uzunluk ölçeğinde deyil: bşəbəkə qırılması femtosaniye və ya nanosaniye düzeyinde gerçekleşirken səth aşınması və geometrik değişim saniyeler boyunca birikebilir.

Atomistik model nasıl kurulmuştur?

Silisium karbidün termomekanik və oksidasyon davranışı, Large-scale Atomic/Molecular Massivəly Parallel Simulator (LAMMPS) yazılımıyla yürütülen klasik molekulyar dinamika simulyasiyalarıyla modellenmiştir. Kimyasal bşəbəkəların kırılmasına, yeni bşəbəkəların oluşmasına və atomlar arasında yük aktarımına izin vərmek üçün reaktif kuvvət alanı ReaxFF ilə yük dengeleme modeli QEq istifadə edilmişdir.

ReaxFF sisteminin toplam potansiyel enerjisi tədqiqatda şu biləşenlerle ifade edilir:

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

• U_bond: Bşəbəkəlı atomların bşəbəkə enerjisi katkısıdır.
• U_ovər: Atomların fiziki olmayan aşırı koordinasyonunu sınırlar.
• U_under: Eksik koordine olmuş atomların kararlılığına katkı sşəbəkəlar.
• U_valence: Bşəbəkə açılarına ilişkin enerji terimidir.
• U_lp: Yalnız elektron çifti enerjisidir.
• U_torsion: Dört atomlu burulma açılarından doğan enerjidir.
• U_vdW: Kısa menzilli van der Waals etkiləşimlerini temsil eder.
• U_Coul: Uzun menzilli Coulomb etkiləşimidir.

Bu toplam enerji, molekulyar dinamikate her atomun üzerine etkiyen kuvvətlerin və zaman üçündeki hareketin hesaplanmasında kullanılır. Modelin önemi, sabit bşəbəkə yapısına sahip geleneksel potansiyellerden fərqli kimi oksidasyon sırasında yeni Si–O və C–O bşəbəkəlarının oluşabilmesine izin vərmesidir.

Isıl ilətkenlik nasıl hesablanmışdır?

Polikristal və tek kristal α(6H)-SiC yapılarının ısıl ilətkenliği Green–Kubo ilişkisiyle hesablanmışdır:

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

Burada k_ij istilik keçiriciliyi tensörünü, V simulyasiya hücresinin hacmini, k_B Boltzmann sabitini, T mutlak sıcaklığı və J istilik axınını temsil eder. Köşeli parantez üçündeki terim, başlangıçtaki istilik axını ilə daha sonraki bir zamandaki istilik axını arasındaki otokorelasyondur. Korelasyonun uzun süre korunması, enerjinin yapı üçünde daha etkin taşındığı anlamına gelir.

Termomekanik hesablamalarda təxminən 20 × 20 × 50 nm boyutunda bir simulyasiya hücresi istifadə edilmişdir. X və Y doğrultularında periyodik sınır koşulları uygulanmış, ısıl enerji Şekil 2’de gösterildiği gibi sıcak bölgeden soğuk bölgeye aktarılmıştır. Görselde fərqli renkler, polikristal yapıdaki ayrı taneleri göstərir.

Özgül ısı və basınç hesablamaları

Özgül ısı kapasitesi üçün LAMMPS üçündeki harmonik eşlemeli ortalama metodi (harmonically mapped avəraging, HMA) istifadə edilmişdir. Makaledeki enerji, basınç və sabit hacimde xüsusi istilik ifadeleri şöyledir:

\[ \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 tənliklerde d sistemin boyut sayısını, N atom sayısını, U toplam enerjiyi, F atomik kuvvətleri, Δr atomların kafes konumlarından yer değiştirmelerini, Φ Hessian matrisini, ρ sayı yoğunluğunu və P_vir virial basıncı göstərir. Enerji dalgalanmalarının büyüklüğü, sabit hacimde xüsusi istilik hesabına doğrudan katkı vərmektedir.

Oksidasyon və ablasiya simulyasiyau nasıl uygulanmıştır?

Oksidasyon modeli üçün 20 × 20 nm yatay boyuta və 15 nm düşey boyuta sahip ayrı bir hücre kurulmuştur. Hücrede 10 nm kalınlığında SiC tabakası ilə bunun üzərində sərhəd qatını temsil eden 5 nm’lik boşluk vardır. Moleküler və atomik oksigen, fərqli hız və temperaturlarda SiC səthine yönlendirilmiştir.

Şekil 3’te alt bölümdeki kırmızı atomlar sabitlenmiş, yeşil atomların sıcaklığı termostatla kontrol edilmiş və oksigen atomları üst bölgeden sisteme eklenmiştir. Hücrenin üstünde, səthdən yeterince uzaklaşan növlərin simulyasiyadan silindiği bir bölge tanımlanmıştır. Böylece oksigenin səthe çarpması, kimyəvi bşəbəkə oluşturması və bazı ürünlerin səthdən ayrılması atom ölçeğinde izlenmiştir.

Langevin termostatındaki toplam kuvvət şu şekilde ayrıştırılmıştır:

\[ \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 potansiyelinden gelen korunumlu kuvvəti, F_f sürtünmə kuvvətini və F_r rastgele kuvvət terimini temsil eder. m atom kütlesi, V hız, τ_damp temperatur gevşeme süresi və dt zaman adımıdır. Rastgele kuvvət tənliki burada makalede yazıldığı biçimiyle aktarılmıştır.

Moleküler dinamik hangi uçucu növləri öngörmüştür?

Simülasyonlarda 1500 °C’nin altında SiC səthinden ayrılan başlıca növlərin karbon monoksit (CO), atomik oksigen (O) və karbondioksit (CO₂) olduğu bildirilmişdir. Daha yüksək temperaturlarda C₂O₂, C₂O₃ və SiO₃ gibi daha karmaşık və metastabil növlərin oluşabiləceği öngörülmüştür.

Buna karşın ATAC-1’e aktarılan ilk B-prime tablosu hesap yükünü azaltmak amacıyla yalnızca CO, SiO₂ və CO₂ növləriyle sınırlandırılmıştır. Dolayısıyla molekulyar dinamikate gözlenen bütün reaksiyon ürünleri makro ölçekli səth geriləməsi hesabına taşınmamıştır. Azotla və ultra yüksək temperatur seramiklerinde istifadə olunan geçiş metalleriyle gerçekleşebiləcek reaksiyonlar da mevcut sürümde yoxdur.

ATAC-1 nasıl tədqiqatktadır?

ATAC-1, tədqiqatda zaman zaman CFD kodu kimi adlandırılsa da tam Navier–Stokes tənliklərini çözmemektedir. Kod, şok sonrası axıntan başlayan inviscid akım çizgilərini izleyerek basınç və hız dşəbəkəılımını çıkarır; ardından sərhəd qatı yaklaşımları, yarı ampirik korelasyonlar və sonlu farklarla çalışan bir ısıl həll edici kullanır. Yüzey sıcaklığı və material geriləmə hızı, her geometrik noktadaki doğrusal olmayan enerji dengesi Newton–Raphson metodiyle çözülerek belirlenir.

Bu yaklaşım yüksək doğruluklu CFD həll edicilərinden daha hızlıdır; ancak viskoz axın ayrıntılarını, turbulentlikı, karmaşık şok etkiləşimlerini və ayrılmış axın bölgelerini doğrudan çözemez. Tədqiqatnın amacı da nihai doğrulukta uçuş tahmini yapmak yerine, yeni B-prime tablosunun hızlı bir mühəndislik aracına entegre ediləbildiğini göstermektir.

Hipersonik burun konusu üçün hangi koşullar incelenmiştir?

ATAC-1 hesaplarında 1976 ABŞ Standart Atmosferi və hava üçün Clark Mollier tablosu istifadə edilmişdir. Yalpa və dönme açıları sıfır kabul edilmişdir. İki əsas hündürlük incelenmiştir:

• 15.240 m: 11.664,53 Pa atmosfer basıncı və −56,5 °C temperatur.
• 21.336 m: 4.487,68 Pa atmosfer basıncı və −55,23 °C temperatur.

Geometri, 7,5° etkin yarı koni açısına sahip 10 cm ölçekli simetrik bir burun konusudir. Zaman adımı 0,01 saniye kimi belirlenmiş, səthde P1–P11 kimi adlandırılan 11 ayrıntılı analiz noktası seçilmiştir. Şekil 4, bu noktaların burun ucunda və koninin rüzgâr üstü ilə rüzgâr altı səthlerinde nasıl dşəbəkəıldığını göstərir.

Temel parametrik seri Mach 5, 10, 15 və 20 hızlarında gerçekleştirilmiştir. Hücum açısı 5° ilə birlikte 10°–20° arasında birer derecelik aralıqlarla değiştirilmiş və ana hesapların uçuş süresi 32 saniye kimi alınmıştır. Ayrıca Mach 5 və Mach 10 üçün 19° hücum bucşəbəkəında 900 saniyelik uzun süreli senaryolar çalıştırılmıştır.

Durgunluk sıcaklığı nasıl tahmin edilmişdir?

Tədqiqatda fərqli hündürlük və Mach sayıları üçün basitleştirilmiş durgunluk sıcaklığı ilişkisi de istifadə edilmişdir:

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

T₀ durgunluk sıcaklığını, T_∞ ilgili hündürlükdaki serbest axın sıcaklığını, M Mach sayını, Pr Prandtl sayısını və γ xüsusi istilik oranını temsil eder. Hesapta hava ideal gaz kabul edilmiş və γ değeri 1,4’te sabit tutulmuştur. Bu basit ilişki, hızın kinetik enerjisinin axın durdurulduğunda ne ölçüde ısıl enerjiye dönüşebiləceğini gösterir; ancak gerçek səth reaksiyonlarını, yüksək basınç faz davranışını və denge dışı kimyayı tek başına açıklamaz.

Isıl özellik nəticələrı ne göstərmişdir?

Şekil 5’te MD ilə hesaplanan α-SiC ısıl ilətkenliği, təxminən 750 °C’de 48 W/m·°C düzeyinden 2250 °C civarında təxminən 13 W/m·°C düzeyine düşmektedir. Düşüş ilk temperatur aralığında daha hızlı, təxminən 1600 °C sonrasında ise daha yataydır. Tədqiqatçılar 1250 °C və 1500 °C yakınındaki küçük dalgalanmaları oksidasyon və ya faz geçişi eşikleriyle ilişkiləndirmekte, ancak nəticələrın simulyasiya boyutu, sınır koşulları, termostat gevşeme süresi və ihmal edilən kuantum etkilərinden etkilənebiləceğini kabul etmektedir.

Şekil 6’daki xüsusi istilik değerleri təxminən 750 °C’den 2250 °C’ye kadar düzenli biçimde artmakta və grafikte təxminən 5,4’ten 6,15’e yükselmektedir. Bununla birlikte grafiğin düşey ekseni “J/kg·°C” kimi etiketlenmiş olmasına rşəbəkəmen vərilən rəqəmsal ölçek bu birimle açık biçimde açıklanmamıştır. Bu nedenle grafikteki değerlerin mutlak mühəndislik vərisi kimi doğrudan kullanılmasından əvvəl birim və ölçek doğrulaması gereklidir.

Hız, hündürlük və hücum bucşəbəkəı sıcaklığı nasıl etkiləmiştir?

Şekil 7’de genel eğilim, Mach sayı yükseldikçe maksimum səth sıcaklığının artmasıdır. 15.240 m hündürlükdaki Mach 20 serisi, çoğu hücum bucşəbəkəında SiC üçün modelde istifadə olunan 2850 °C sublimasiya sınırına yaklaşmakta və ya bu sınırı aşmaktadır. En belirgin aşım 17° hücum bucşəbəkəında görülmekte və grafikte maksimum temperatur təxminən 3200 °C’ye çıkmaktadır.

Daha yüksək olan 21.336 m hündürlükda aşşəbəkəı atmosfer yoğunluğu nedeniyle Mach 5 temperaturları 15.240 m nəticələrından daha aşşəbəkəıtür. Ancak Mach 15 və Mach 20 gibi yüksək hızlarda temperaturların yine 2600–2850 °C aralığına ulaştığı görülmektedir. İrtifa etkisi bu nedenle tek başına değerlendiriləmez; hız, səth basıncı, hücum bucşəbəkəı və istifadə olunan axın yaklaşımıyla birlikte ele alınmalıdır.

Malzeme geriləməsi nəticələrı neden sorunludur?

Şekil 8, fiziki eğilimlerin yanında modelin geçerlilik sınırını aşan nəticələr da üretmiştir. Mach 20, 15.240 m və 17° hücum bucşəbəkəı serisinde səth geriləməsi təxminən 3,10 × 10¹⁴ cm kimi gösterilmektedir. Bu değer herhangi bir gerçek araç geometrisiyle bşəbəkədaşmamaktadır. Metnin bir bölümünde 17°–18° anomalisi Mach 15’e atfedilirken grafiğin lejandı olşəbəkəanüstü değeri Mach 20 və 15 km serisine yerleştirmektedir. Dolayısıyla nəticə anlatımı ilə grafik arasında Mach numarası bakımından tutarsızlık vardır.

Tədqiqatçılar bu anomalinin, temperatur və basınç B-prime tablosunun kalibre edildiği alanın dışına çıktığında ATAC-1’in material davranışını üstel biçimde dışa değerlemesinden mənbəlanabiləceğini açıklamaktadır. Yüzey sıcaklığı 2850 °C sublimasiya sınırını geçtiğinde model, katı material bütünüyle kaybolmuş olsa bilə geriləmə hesabını sürdürmektedir. Başlangıçtaki material kalınlığı tüketildikten sonraki zaman adımları fiziki anlam taşımamaktadır.

Yüzeyde en fazla aşınma nerede oluşmuştur?

Şekil 9’daki Mach 15 səth haritaları, material kaybının büyük bölümünün burun ucunda və uçuş yönüne bakan rüzgâr üstü tarafta toplandığını göstərir. Hücum açısı değiştikçe yüksək geriləmə bölgesi simetrik burun merkezinden yana doğru kaymaktadır. Bu nəticə, fərqli hücum açılarında eyni toplam temperatur oluşsa bilə səth üzərindəki yerel yük dşəbəkəılımının değişebildiğini göstərir.

Ancak haritanın renk ölçeği 32 saniye üçünde 80 cm’ye kadar toplam geriləmə göstərir. İncelenen burun geometrisinin 10 cm ölçeğinde olduğu dikkate alındığında bu değerler de gerçek materialnin hesap süresi tamamlanmadan tükenmiş olacşəbəkəını göstərir. Bu haritalar nicel ömür tahmini kimi değil, modelin hangi səth bölgelerinde yüksək yük öngördüğünü gösteren eğilim haritaları kimi okunmalıdır.

Mach 5 və Mach 10 uzun süreli senaryoları

15.240 m hündürlük, 19° hücum bucşəbəkəı və 900 saniye uçuş süresinde Mach 5 ilə Mach 10 karşılaştırılmıştır. Mach 5 senaryosunda səth sıcaklığı təxminən 900 °C civarına ulaşmış və SiC üçün tədqiqatda belirtilən təxminən 1400 °C geriləmə başlangıç düzeyinin altında kalmıştır. Şekil 10’daki səth kaybı çok küçük bir rəqəmsal ölçek üçünde kalmaktadır.

Mach 10 senaryosunda ise temperatur 2500 °C’nin üzərində hızlı dalgalanmalar göstermiş və hesaplanan toplam geriləmə təxminən 101,6 cm’ye ulaşmıştır. Bu değer başlangıç geometrisini aşmaktadır. Yazarlar, hızlı temperatur salınımlarının zaman adımı çözünürlüğü və ablasiya hesabıyla ilgili məhdudiyyətləri yansıtabiləceğini; geometrinin tamamen kaybolmasından sonraki hesapların geçersiz olduğunu belirtmektedir.

Mach 7 örneği ne göstərmişdir?

Basitleştirilmiş durgunluk sıcaklığı hesabı, 15.240 m hündürlükda Mach 7 hızında SiC səthinin ablasiya başlangıç sıcaklığını aşabiləceğini düşündürmüştür. Buna karşılık ATAC-1 ilə 0° hücum bucşəbəkəı və 20 saniye üçün yapılan hesapta səth basıncı 2,43 MPa’ya ulaşmış və yalnızca yuvarlama xətası düzeyinde material kaybı oluşmuştur.

Bu uyuşmazlık, ablasiyaun yalnızca temperatur eşiğiyle belirlenemeyeceğini göstərir. Yüksek basınçtaki erime, buharlaşma və sublimasiya sınırlarının güvənilir faz diyagramlarıyla tanımlanması lazımdır. Tədqiqatda istifadə olunan basit ideal gaz və sabit özellik varsayımları, bu geçişleri bütün ayrıntılarıyla temsil etmemektedir.

Tədqiqatnın güclü tərəfləri nelerdir?

• Atomistik reaksiyonları hızlı bir aerotermal tasarım aracına aktaran açık bir çok ölçekli iş axını sunmaktadır.
• Bşəbəkə kırılması, bşəbəkə yaranması və yük aktarımını temsil edebilən ReaxFF istifadə edilmişdir.
• Isıl ilətkenlik, xüsusi istilik, səth kimyası, uçucu növlər, temperatur və geometrik geriləmə eyni model zincirinde ele alınmıştır.
• Farklı Mach sayıları, iki hündürlük, çok sayıda hücum bucşəbəkəı və iki ayrı uçuş süresi karşılaştırılmıştır.
• Yazarlar modelin fiziki olmayan nəticələrını gizlememiş; dışa değerleme və geometri tükenmesi sorunlarını açıkça belirtmiştir.
• Yöntem, yüksək maliyetli deneylerden əvvəl geniş bir parametre alanının taranması üçün potansiyel bir ön tasarım yaklaşımı sunmaktadır.

Tədqiqatnın məhdudiyyətləri nelerdir?

• Bu tədqiqat hakem qiymətləndirməsindən keçməmiş bir preprinttir.
• ATAC-1 tam Navier–Stokes tənliklərini çözmemekte və ayrılmış axın bölgelerini temsil edememektedir.
• B-prime tablosu yalnızca Si–C–O etkiləşimlerini içermekte; azot və diğer seramik biləşenleri dışarıda bırakılmaktadır.
• Moleküler dinamikte öngörülen növlərin tamamı makro ölçekli modele aktarılmamıştır.
• Emisivite və spektral soğurma bütün temperaturlar üçün sabit 0,87 kabul edilmişdir.
• Isıl ilətkenlik nəticələrında eksperimental eğilimlerden kənaraçıxmalar bulunduğu yazarlar tarafından kabul edilir.
• Şekil 6’daki xüsusi istilik birimi və ya ölçeklendirmesi yeterince açık deyil.
• Metinde atıf yapılan Tablo 1 yüklenen preprint sürümünde yoxdur.
• Bazı səth geriləməsi nəticələrı başlangıç geometrisini və fiziki kimi mümkün değerleri aşmaktadır.
• Şekil 8 ilə açıklama metni arasında anormal değerin Mach 15 mi yoksa Mach 20 serisine mi ait olduğu konusunda tutarsızlık vardır.
• Model nəticələrı yüksək doğruluklu CFD, hipersonik şok tüneli və ya gerçek uçuş vərisiyle doğrulanmamıştır.
• Belirsizlik yayılımı və nicel xəta aralıqları sunulmamıştır.

Tədqiqat neyi dəstəkləyir?

Tədqiqat, reaktiv molekulyar dinamikaten türetilən termal və kimyəvi bilgilərin B-prime biçimine dönüştürülerek bir aerotermal həll ediciye aktarılabiləceğini göstərir. Ayrıca temperatur yükseldikçe SiC’nin hesaplanan ısıl ilətkenliğinin azaldığını, xüsusi istiliknın arttığını və uçucu səth növlərinin daha karmaşık hâle gelebildiğini göstərir. Model, yüksək temperatur və yüksək Mach koşullarında burun ucu ilə rüzgâr üstü səthin daha yüksək ablasiya riski taşıdığı yönünde bir eğilim üretmektedir.

Tədqiqat neyi sübut etmir?

Tədqiqat, vərilən geriləmə değerlerinin gerçek bir hipersonik aracın ömrünü doğru tahmin ettiğini sübut etmir. Simülasyonlar herhangi bir uçuş aracının başarıyla üretildiğini, SiC koruma sisteminin belirli bir görəv süresine dayanacşəbəkəını və ya Mach 20 koşullarında güvənli çalışacşəbəkəını göstermemektedir. Fiziksel olmayan büyük nəticələr, mevcut sürümün sertifikasyon, güvənlik sınırı və ya nihai boyutlandırma amacıyla kullanılamayacşəbəkəını göstərir.

Türkiyə açısından mümkün əhəmiyyəti nedir?

Tədqiqatnın metodsel yaklaşımı; Türkiyə’de hipersonik hava araçları, füze və roket burun geometriləri, yeniden giriş sistemleri, yüksək temperatur seramikleri və termal koruma materialleri üzərində çalışan tədqiqat grupları açısından aktarılabilir niteliktedir. Özellikle pahalı şok tüneli və ya ark-jet deneylerinden əvvəl atomistik hesablamalarla reaksiyon növlərinin belirlenmesi və bu bilgilərin hızlı tasarım kodlarına aktarılması, deney alanının daraltılmasına yardımcı ola bilər.

Bununla birlikte tədqiqatnın kendi nəticələrı, hızlı və aşşəbəkəı doğruluklu modellerin mutlaka yüksək doğruluklu CFD, material deneyleri və uygun faz diyagramlarıyla sınanması gerektiğini de göstərir. Türkiyə’ye aktarılabiləcek en güçlü ders belirli bir geriləmə değeri değil; atomistik model, mühəndislik həll edicisi və eksperimental doğrulama arasında kademeli bir model doğrulama zinciri kurulması gereğidir.

Tədqiqatın Metodu və Nəticələri

Başlıca rəqəmsal nəticələr

• Hesaplanan istilik keçiriciliyi təxminən 750 °C’de 48 W/m·°C’den 2250 °C’de təxminən 13 W/m·°C’ye düşmüştür.
• Şekil 6’daki xüsusi istilik eğrisi təxminən 5,4’ten 6,15’e yükselmiştir; ancak grafiğin “J/kg·°C” birim etiketiyle rəqəmsal ölçeğin ilişkisi açıklanmamıştır.
• 1500 °C’nin altında başlıca uçucu növlər CO, O və CO₂ kimi öngörülmüştür.
• 1500 °C üzərində C₂O₂, C₂O₃ və SiO₃ növlərinin oluşabiləceği bildirilmişdir.
• SiC üçün ATAC-1’de istifadə olunan sublimasiya sıcaklığı 2850 °C’dir.
• Mach 20, 15.240 m və 17° hücum bucşəbəkəında maksimum səth sıcaklığı grafikte təxminən 3200 °C’ye ulaşmıştır.
• Aynı parametre bölgesinde Şekil 8 təxminən 3,10 × 10¹⁴ cm gibi fiziki olmayan bir səth geriləməsi üretmiştir.
• Mach 15, 15.240 m və 32 saniyelik səth haritalarında ən yüksək geriləmə burun ucu və rüzgâr üstü bölgede toplanmıştır.
• Mach 5, 19° hücum bucşəbəkəı və 900 saniyede səth sıcaklığı təxminən 900 °C düzeyinde kalmış və geriləmə ihmal ediləbilir ölçekte hesablanmışdır.
• Mach 10, eyni hündürlük və sürede təxminən 101,6 cm geriləmə üretmiş; yazarlar bu değerin başlangıç geometrisini aştığı üçün fiziki kimi geçersiz olduğunu belirtmiştir.
• Mach 7, 0° hücum bucşəbəkəı və 20 saniyede səth basıncı 2,43 MPa kimi hesaplanmış və temperatur eşik tahminine rşəbəkəmen belirgin ablasiya oluşmamıştır.

Modelin mühəndislikte kullanılabilmesi üçün gereken sonraki adımlar

1. B-prime tablosunun azot və ultra yüksək temperatur seramiklerindeki diğer elementleri içerecek şekilde genişletilmesi.
2. Emisivite və spektral soğurmanın sıcaklığa bşəbəkəlı değişkenler kimi tanımlanması.
3. MD hesaplarının daha büyük hücreler, daha uzun gevşeme süreleri və qeyri-müəyyənlik analizleriyle tekrarlanması.
4. Tablonun ATAC-2, US3D, LAURA və ya DPLR gibi daha yüksək doğruluklu həll ediciləre aktarılması.
5. Şok tüneli, ark-jet və ya kontrollü ablasiya deneyleriyle səth sıcaklığı və geriləmə değerlerinin doğrulanması.
6. Malzeme tamamen tükendiğinde hesablamayı durduracak geometrik və fiziki sınırların uygulanması.
7. Yüksek temperatur və basınçlarda güvənilir SiC faz diyagramlarının modele eklenmesi.

Mənbə və Metod Qeydi

Tədqiqatnın özgün adı: Hypersonic Flight Aerothermodynamics and Ablation: A Theoretical Model

Yazarlar və sıralama: Robert Slapikas, Luis Bravo, Anindya Ghoshal.

Eş birinci yazar: Tədqiqatda eş katkı və ya eş birinci yazarlık bilgisi yoxdur.

İletişim yazarı: SSRN kaydında Anindya Ghoshal ilətişim yazarı kimi gösterilmektedir.

Kurumsal bşəbəkəlantılar: Robert Slapikas, DEVCOM Army Research Laboratory və SURVICE Engineering Company bşəbəkəlantılıdır. Luis Bravo ilə Anindya Ghoshal, DEVCOM Army Research Laboratory bşəbəkəlantılıdır.

Yazar katkıları: Robert Slapikas; tədqiqat, metod, yazılım, görselleştirme, doğrulama, vəri düzenleme, biçimsel analiz və ilk taslşəbəkəın yazımını üstlenmiştir. Luis Bravo; metod, ilk taslşəbəkəın yazımı, danışmanlık və mənbə sşəbəkəlamaya katkıda bulunmuştur. Anindya Ghoshal; metod, yazım, danışmanlık, mənbəlar, proje yönetimi və finansman teminine katkı sşəbəkəlamıştır.

DOI:10.2139/ssrn.6952638

Yayın platformu: SSRN.

Özgün yayınevi və ya platform işletmecisi: Tədqiqat hakemli bir dergide yayımlanmamıştır. SSRN, Elsevier tarafından işletilən bir preprint və erken tədqiqat platformudur.

Yayın yılı: 2026.

SSRN’ye yüklenme tarihi: 16 Haziran 2026.

Hakemlik durumu: Bu tədqiqat bir preprinttir və hakem qiymətləndirməsindən keçməmiştir.

Kaynak türü: Moleküler dinamik və aerotermal hesablamalar içeren modelleme/simulyasiya preprinti.

Resmî mənbə:SSRN resmî tədqiqat sayfası

Finansman: Tədqiqat DEVCOM Army Research Laboratory tarafından desteklenmiş və W15P7T-19-D-0126 numaralı iş birliği anlaşması kapsamında yürütülmüştür.

Çıkar çatışması: Yazarlar bilinen bir mali çıkar çatışması və ya tədqiqatyı etkiləyebiləcek kişisel ilişki bulunmadığını bildirmiştir.

Bu Türkçe açıklama, yüklenen 36 sayfalık tədqiqat; tənlikləri, simulyasiya koşulları, grafikler, səth haritaları, nəticələr və yazarların belirttiği məhdudiyyətlər incelenerek hazırlanmıştır. Tədqiqatda bulunmayan bilimsel nəticələr eklenmemiştir. Bilimsel içerik üçün dış mənbə kullanılmamış; dış kontrol yalnızca DOI, platform, yüklenme tarihi və ilətişim yazarı gibi bibliyografik kimlik bilgiləriyle sınırlandırılmıştır.

Yüklenen sürümde metnin atıf yaptığı Tablo 1 yoxdur. Şekil 6’daki xüsusi istilik birimi və ya ölçeği yeterince açıklanmamıştır. Ayrıca Şekil 8’deki olşəbəkəanüstü büyük geriləmə değeri ilə açıklama metnindeki Mach numarası arasında tutarsızlık vardır. Bu nedenlerle tədqiqat, nihai bir araç ömrü və ya material güvənliği tahmini kimi değil, atomistik vərinin hızlı bir mühəndislik həll edicisine aktarılmasını gösteren erken aşama bir hesablamalı tədqiqat kimi değerlendirilmelidir.

Bu tədqiqat gerçek uçuş başarısı, saha doğrulaması və ya sertifikasyon kanıtı sunmamaktadır. Modelin nicel nəticələrı; yüksək doğruluklu CFD, kontrollü hipersonik deneyler və qeyri-müəyyənlik analiziyle doğrulanmadan gerçek araç tasarımında doğrudan kullanılmamalıdır.