Hipersonik araçların çevresindeki akışı hesaplayan sayısal yöntemlerin güçlü şok dalgalarını yakalaması, yoğunluk ve basıncı pozitif tutması, kararlı çalışması ve fiziksel korunum yasalarını mümkün olduğunca doğru temsil etmesi gerekir. Ancak yaygın Rusanov, HLL, HLLC ve Roe gibi Godunov tipi çözücüler, kararlı Euler akışlarında akım çizgileri boyunca sabit kalması gereken toplam entalpiyi sayısal olarak korumayabilir. Şoktan geçen akışta oluşan yapay entalpi değişimi, daha sonra akış alanının diğer bölgelerine taşınarak enerji ve yoğunluk dağılımını bozabilir.

Araştırmada, Gallice tarafından geliştirilen Simple Riemann çözücülerinin enerji bileşeni değiştirilerek MGallice ve MGallice-2D adlı entalpi koruyan tek ve çok boyutlu şemalar oluşturuldu. Değişiklik, enerji akısındaki sayısal yayınımın kütle akısındaki yayınım ile toplam entalpi çarpımına eşit olmasını sağlayan Häenel koşuluna dayanıyor.

Mach 2 hızındaki üçgen çıkıntılı kanal testinde dört şema da yaklaşık 4.000 iterasyonda kararlı çözüme ulaştı. Klasik Gallice çözücüleri durağan şokların üzerinden geçerken küçük fakat belirgin toplam entalpi sıçramaları üretti. Değiştirilmiş MGallice çözücüleri ise giriş toplam entalpisi olan 6,3 değerini bütün hesaplama alanında korudu.

Mach 20 hızında yarım silindir çevresindeki daha zor testte tek boyutlu çözücüler karbunkül adı verilen şok kararsızlığına maruz kaldı. Çok boyutlu çözücüler bu kararsızlığı önledi; ancak klasik Gallice-2D, şok arkasındaki entalpi kaybı nedeniyle durma noktasındaki yoğunluk dağılımını yanlış hesapladı. MGallice-2D toplam entalpiyi makine hassasiyetine kadar sabit tutarak beklenen yoğunluk alanını geri kazandırdı.

Çalışma, kararlı hipersonik Euler hesaplarında toplam entalpi korunumunun yalnızca teorik bir özellik olmadığını; güçlü şokların arkasındaki enerji ve yoğunluk alanını doğrudan değiştirebildiğini gösteriyor. Bununla birlikte önerilen şemaların entropi kararlılığı matematiksel olarak kanıtlanmış değildir ve yöntem henüz viskoz Navier–Stokes problemleri üzerinde sınanmamıştır.

Hipersonik akışların sayısal çözümü neden zordur?

Hipersonik hızlarda ilerleyen bir aracın önünde güçlü ve araçtan ayrılmış bir yay şoku oluşur. Şok boyunca basınç, yoğunluk ve sıcaklık kısa bir mesafede büyük ölçüde değişir. Aracın arkasında genişleme dalgaları meydana gelirken yüzeye yakın bölgelerde hız ve sıcaklık gradyanları oluşur.

Gerçek hipersonik akışlar viskozite, ısı iletimi, kimyasal reaksiyonlar, termokimyasal denge dışı etkiler, türbülans, yüzey ısınması ve aşınma gibi çok sayıda süreci içerebilir. Çalışma bu bütün sistemi değil, sıkıştırılabilir Navier–Stokes denklemlerinin şok oluşumundan sorumlu viskozitesiz bölümünü, yani Euler denklemlerini ele alıyor.

Sayısal bir şema güçlü şoku kararlı biçimde yakalasa bile fiziksel olarak korunması gereken bir büyüklükte yapay değişim üretebilir. Araştırmanın odaklandığı temel sorun, kararlı Euler akışlarında toplam entalpinin sayısal olarak korunmamasıdır.

Toplam entalpi nedir?

İki boyutlu Euler denklemlerinde korunan değişkenler yoğunluk, momentum ve toplam enerjidir:

\[ \mathbf{U}= \begin{pmatrix} \rho \\ \rho\mathbf{u} \\ \rho e \end{pmatrix} \]

Burada ρ yoğunluğu, u hız vektörünü ve e birim kütle başına toplam enerjiyi gösterir. Toplam enerji, iç enerji ile kinetik enerjinin toplamıdır:

\[ e=\varepsilon+\frac{1}{2}\mathbf{u}\cdot\mathbf{u} \]

Toplam entalpi ise:

\[ h=e+\frac{p}{\rho} \]

bağıntısıyla tanımlanır. p basınçtır. Çalışmada mükemmel gaz denklemi kullanılmıştır:

\[ p=\rho\varepsilon(\gamma-1) \]

γ politropik veya özgül ısı oranıdır. Sayısal testlerde γ = 7/5 olarak alınmıştır.

Kararlı Euler akışında entalpi neden sabit kalır?

İki boyutlu Euler denklemleri kütle, momentum ve toplam enerji korunumunu şu genel biçimde ifade eder:

\[ \frac{\partial\mathbf{U}}{\partial t} + \nabla\cdot\mathbf{F}(\mathbf{U}) =0 \]

Kararlı durumda zamana bağlı terim ortadan kalkar. Kütle ve enerji denklemleri birlikte değerlendirildiğinde:

\[ \mathbf{u}\cdot\nabla h=0 \]

sonucu elde edilir. Bu ifade, viskozite ve ısı iletimi bulunmayan kararlı bir Euler akışında toplam entalpinin akım çizgisi boyunca değişmediğini gösterir.

Akım çizgileri giriş sınırına bağlıysa toplam entalpi bütün hesaplama alanında giriş değeri h_∞ ile aynı kalmalıdır. Şok basıncı, sıcaklığı ve yoğunluğu değiştirir; ancak ideal kararlı Euler çözümünde toplam entalpi şok boyunca korunur.

Sayısal şemalar bu özelliği neden kaybedebilir?

Sonlu hacim yöntemlerinde hesaplama alanı hücrelere ayrılır. Her hücrenin kütle, momentum ve enerjisindeki değişim, hücre yüzeylerinden geçen sayısal akılarla belirlenir. Yüzeyin iki tarafındaki durumlar arasında bir Riemann problemi kurulur ve bu problemin yaklaşık çözümünden sayısal akı elde edilir.

Rusanov, HLL, HLLC ve Roe gibi klasik Godunov tipi şemalar korunumlu olmalarına rağmen kütle ve enerji akılarına aynı fiziksel ilişkiyi sağlayacak biçimde sayısal yayınım eklemez. Bu nedenle durağan bir şok boyunca enerji akısındaki yapay yayınım, kütle akısındaki yapay yayınımın toplam entalpiyle çarpımına eşit olmayabilir.

Sonuçta şokun bir tarafında hesaplanan toplam entalpi diğer tarafındakinden farklı hale gelir. Bu hata daha sonra akımla birlikte taşınır ve özellikle şok arkasındaki düşük hızlı bölgelerde enerji dağılımını etkileyebilir.

Entalpi hatasının mühendislikteki önemi

Araştırmacı, bazı Navier–Stokes hesaplarında sınır tabakasını tanımlayan ölçütlerin toplam entalpiye dayanabileceğini belirtiyor. Şoktan kaynaklanan yapay entalpi izi duvar boyunca taşındığında, fiziksel sınır tabakası ile sayısal enerji hatasının birbirinden ayrılması güçleşebilir.

Hipersonik araçlarda durma noktası çevresindeki yoğunluk, basınç, sıcaklık ve ısı akısı büyük önem taşır. Çalışmanın Mach 20 testi, toplam entalpi hatasının yalnızca yardımcı bir değişkeni bozmadığını; yoğunluk alanının biçimini de değiştirebildiğini gösteriyor.

Simple Riemann çözücüsü nedir?

Simple Riemann çözücüleri, başlangıçtaki sol ve sağ durumlar arasında sabit ara durumlar tanımlar. Bu ara durumlar farklı hızlarda hareket eden süreksizliklerle birbirinden ayrılır:

\[ W(\xi)= \begin{cases} U_l, & \xi\leq\lambda_1 \\ U_k, & \lambda_{k-1}<\xi\leq\lambda_k \\ U_r, & \xi>\lambda_m \end{cases} \]

Burada ξ = x/t kendine benzer değişkeni, U_l ve U_r sol ve sağ durumları, λ_k ise dalga hızlarını gösterir.

Gallice yaklaşımında çözücü önce Lagrange koordinatlarında kuruluyor. Daha sonra Lagrange–Euler dönüşümüyle sabit uzay hücreleri üzerinde çalışan Euler biçimine çevriliyor. Bu yöntem ara durumların pozitif yoğunluk ve iç enerji koşullarının incelenmesini kolaylaştırıyor.

Tek boyutlu Gallice şeması

Tek boyutlu olarak adlandırılan şemada bir hücre yüzeyindeki sayısal akı yalnızca yüzeyin iki tarafındaki hücrelerden hesaplanıyor. İki boyutlu ağ kullanılsa bile her yüzeydeki Riemann problemi yüzey normaline göre tek boyutlu ele alınıyor.

Ara durumun normal hızı:

\[ u_f= \frac{ \lambda_lu_{n,l} + \lambda_ru_{n,r} - \Delta p }{ \lambda_l+\lambda_r } \]

olarak tanımlanıyor. λ_l ve λ_r Lagrange dalga hızlarını, u_n yüzeye normal hız bileşenini ve Δp iki hücre arasındaki basınç farkını gösteriyor.

Makale, Euler biçimindeki tek boyutlu Gallice şemasının özel dalga hızı kestirimleri kullanan HLLC çözücüsüne büyük ölçüde eşdeğer olduğunu belirtiyor.

Çok boyutlu Gallice-2D şeması

Çok boyutlu şemada akı yalnızca yüzey komşularından değil, aynı düğümü paylaşan hücrelerin tamamından etkileniyor. Sayfa 4’teki Şekil 1, klasik yüzey akıları ile bir düğüm çevresinde hesaplanan çok boyutlu akıları karşılaştırıyor.

Ara durumun normal hızı, düğüme ait ortak bir hız vektörünün yüzey normaline izdüşümü olarak tanımlanıyor:

\[ u_n^*=\mathbf{u}_p\cdot\mathbf{n} \]

Düğüm hızı u_p, düğüm çevresindeki bütün yüzeylerin korunum koşulunu sağlayan doğrusal sistemden elde ediliyor. Bu çok boyutlu şablon, sayısal yayınımı akış ve ağ geometrisine bağlı olarak yönlendirebiliyor.

Çalışmaya göre çok boyutlu yapı, güçlü yay şoklarında görülen karbunkül kararsızlığına karşı tek boyutlu çözücülerden daha dirençlidir.

Pozitiflik ve kararlılık koşulları

Bir gaz dinamiği çözümünde yoğunluk ve iç enerjinin negatif olmaması gerekir. Çalışmada tek boyutlu çözücü için Lagrange dalga hızlarının belirli alt sınırları aşması halinde ara durumların pozitif kalacağı gösteriliyor.

Önerilen genel koşul:

\[ \lambda_l> \max \left( \rho_lc_l,\, \sqrt{\rho_l|\Delta p|},\, -\rho_l\Delta u_n \right) \]

\[ \lambda_r> \max \left( \rho_rc_r,\, \sqrt{\rho_r|\Delta p|},\, -\rho_r\Delta u_n \right) \]

biçimindedir. c ses hızını gösterir. Hesaplamalarda her yüzeyde sol ve sağ koşullardan elde edilen en büyük eğim kullanılmıştır.

Zaman adımı için Courant–Friedrichs–Lewy sayısı 0,5 olarak seçilmiştir. Kararlı problemlerde yerel zaman adımı uygulanmıştır.

Häenel koşulu

Toplam entalpi koruyan şemaların temel koşulu, enerji akısındaki sayısal yayınımın kütle akısındaki yayınımla orantılı olmasıdır:

\[ D_{\rho e}=hD_\rho \]

Burada D_ρe enerji akısındaki, D_ρ kütle akısındaki sayısal yayınımı ifade eder. Kararlı akışta h = h_∞ olduğunda enerji ve kütle akısı arasındaki fiziksel bağlantı korunur.

Klasik Gallice çözücüsünün ara enerji durumları bu eşitliği sağlamıyor. Ara durumlar üzerinden hesaplanan enerji farkı, yoğunluk farkı ile toplam entalpi çarpımına eşit olmuyor.

MGallice değişikliği

Araştırmacı, Gallice çözücüsünün kütle ve momentum yapısını korurken enerjiyle ilgili ara durumları toplam entalpi üzerinden yeniden tanımlıyor. Değiştirilmiş Lagrange durum vektörü:

\[ \widetilde{\mathbf{V}}_s= \begin{pmatrix} v_s \\ \mathbf{u}_s \\ h_s \end{pmatrix} \]

olarak kuruluyor. v özgül hacmi, u hızı ve h toplam entalpiyi gösterir.

Euler çerçevesindeki değiştirilmiş ara durumun enerji bileşeni:

\[ \widetilde{\mathbf{U}}_s^* = \rho_s^* \begin{pmatrix} 1 \\ \mathbf{u}_s+(u_n^*-u_{n,s})\mathbf{n} \\ h_s \end{pmatrix} \]

şeklinde oluşturuluyor. Böylece her dalga boyunca enerji bileşenindeki değişim yoğunluk değişimiyle aynı toplam entalpi üzerinden bağlantılı hale geliyor.

Kararlı çözümde sol ve sağ toplam entalpiler h_∞ değerine eşitse:

\[ \rho_l^*e_l^*-\rho_lh_l = (\rho_l^*-\rho_l)h_\infty \]

ve benzer eşitlikler diğer ara dalgalar için de sağlanıyor. Bu yapı Häenel koşulunu hem tek hem çok boyutlu değiştirilmiş Gallice çözücüleri için yerine getiriyor.

Değişikliğin önemli bedeli

Enerji akısının Lagrange biçimindeki karşılığı, değişiklikten sonra yayınımsız merkezî bir akıya dönüşüyor:

\[ \overline{pu} = \frac{1}{2} \left( p_lu_{n,l} + p_ru_{n,r} \right) \]

Kütle ve momentum denklemlerinde yeterli sayısal yayınım bulunması testlerde doğru çözümün yakalanmasına yetmiştir. Ancak yazar, değiştirilmiş MGallice çözücülerinin entropi artıran veya entropi kararlı olduğunu garanti eden bir koşul sunamamaktadır.

Bu nedenle yöntem toplam entalpiyi korurken klasik şemanın kanıtlanmış entropi özelliğinin tamamını otomatik olarak devralmıyor. Bu durum çalışmanın en önemli teorik sınırlılıklarından biridir.

İkinci mertebe MUSCL uzantısı

Birinci mertebe sonlu hacim şemaları güçlü şoklarda kararlı olsa da düzgün bölgelerde yüksek sayısal yayınım üretir. Araştırmada ikinci mertebeli çözüm için MUSCL türü yeniden yapılandırma kullanılmıştır.

• Yoğunluk, hız ve basınç gradyanları düğümlerde ağırlıklı doğrusal en küçük kareler yöntemiyle hesaplanmıştır.
• Hücre merkezi gradyanları, düğüm gradyanlarının normlarının tersine dayalı ağırlıklarla birleştirilmiştir.
• MGallice şemalarında toplam entalpi de yeniden yapılandırılmıştır.
• Gradyanlar R3 Nishikawa sınırlayıcısıyla sınırlandırılmıştır.
• Sınırlayıcıda p = 3 ve hücre ölçeğine bağlı ε_c kullanılmıştır.

Karşılaştırılan dört şema

Birinci test: Mach 2 hızındaki üçgen çıkıntılı kanal

İlk problem, karşılıklı duvarlarında küçük üçgen çıkıntılar bulunan iki boyutlu bir kanaldaki süpersonik akıştır. Alan 2 metre uzunluğunda ve 1 metre yüksekliğindedir. İki çıkıntı kanalın alt ve üst duvarlarına yerleştirilmiştir.

Akış çıkıntılara çarptığında çıkıntıların önünde iki durağan şok oluşuyor. Arkalarında gelişen genişleme dalgaları birbirleriyle etkileşiyor ve üst-alt duvarlardan yansıyarak elmas biçimli bir dalga düzeni meydana getiriyor.

Sayfa 9’daki Şekil 3, bütün şemaların yaklaşık 4.000 iterasyon sonunda kararlı çözüme ulaştığını gösteriyor. Kalıntı eğrileri birbirine yakın olduğundan entalpi değişikliği yakınsama yeteneğini belirgin biçimde bozmamıştır.

Yoğunluk alanları birbirine yakındı

Sayfa 10’daki Şekil 4, dört çözücünün 0,63 ile 1,77 arasındaki yoğunluk izoçizgilerini gösteriyor. Bütün çözücüler şokları, genişleme dalgalarını ve elmas biçimli etkileşim bölgesini benzer biçimde yeniden üretiyor.

Tek boyutlu şemaların biraz daha az yayınımlı olduğu belirtilmiştir. Ancak yalnızca yoğunluk alanına bakıldığında klasik ve değiştirilmiş çözücüler arasındaki temel fark açıkça görülmemektedir.

Entalpi alanları belirgin biçimde farklıydı

Sayfa 10’daki Şekil 5, toplam entalpi alanını 6,2 ile 6,44 arasında gösteriyor. Klasik Gallice şemasında durağan şokların üzerinde küçük entalpi sıçramaları oluşuyor ve bu hatalar duvar boyunca aşağı akıma taşınıyor.

Gallice-2D şemasındaki entalpi bozulması daha belirgin hale geliyor. Buna karşılık MGallice ve MGallice-2D alanları tekdüze görünüyor; giriş toplam entalpisi olan 6,3 değerinin bütün alanda korunduğu bildiriliyor.

Bu test, yoğunluk alanları benzer görünse bile çözücülerin enerji koruma özelliklerinin önemli ölçüde farklı olabileceğini gösteriyor.

İkinci test: Mach 20 hızında yarım silindir

İkinci problem, yarıçapı 1 metre olan yarım silindirin çevresindeki kararlı hipersonik akıştır. Bu problem hem güçlü ayrılmış yay şoku hem de karbunkül kararsızlığı nedeniyle zorlayıcıdır.

İterasyon sayısındaki raporlama tutarsızlığı

Metin, birinci mertebe çözümlerin 10⁵ iterasyon boyunca çalıştırıldığını ve bu noktada kalıntının üç-dört mertebe azaldığını belirtiyor. Ancak sayfa 10’daki Şekil 6’nın yatay ekseni yaklaşık 11.000 iterasyonda sona eriyor.

Bu nedenle belirtilen 10⁵ değeri ile gösterilen grafik arasında preprint içi bir raporlama tutarsızlığı bulunmaktadır. Grafik, çok boyutlu şemaların gösterilen aralıkta daha düşük ve daha düzenli azalan kalıntılar ürettiğini göstermektedir.

Karbunkül kararsızlığı

Karbunkül, güçlü ve ağ çizgileriyle hizalanmış şokların bazı şok yakalama şemalarında fiziksel olmayan biçimde bozulmasıdır. Yarım silindirin önündeki yay şokunun durma çizgisi, bu kararsızlığın ortaya çıkabildiği klasik bir durumdur.

Çalışmada Gallice ve MGallice adlı tek boyutlu şemaların güçlü yay şokunu doğru biçimde çözemediği belirtiliyor. Entalpi koruma değişikliği tek başına karbunkülü ortadan kaldırmamıştır.

Gallice-2D ve MGallice-2D şemaları ise çok boyutlu akı yapıları sayesinde karbunkül içermeyen yay şokları üretmiştir. Bu sonuç, şok kararlılığı ile toplam entalpi korunumunun iki ayrı sayısal özellik olduğunu gösteriyor.

Entalpi korumak yoğunluk alanını düzeltti

Sayfa 10’daki Şekil 7, dört şemanın 1 ile 6,3 arasındaki yoğunluk izoçizgilerini gösteriyor. Klasik Gallice-2D şemasında en yüksek yoğunluk beklenen durma noktasında oluşmuyor. Bu, yay şoku arkasındaki enerji dağılımının fiziksel olmayan biçimde değiştiğini gösteriyor.

MGallice-2D aynı çok boyutlu akı şablonunu kullanmasına rağmen entalpi koruyan enerji değişikliği sayesinde yoğunluk maksimumunu durma noktası çevresine geri getiriyor.

Sayfa 10’daki Şekil 8’de klasik çözücülerin şok üzerinde toplam entalpi düşüşü oluşturduğu görülüyor. Bu yapay eksilme şok arkasındaki sesaltı bölgede enerji dağılımını değiştiriyor ve sonuçta yoğunluk alanını bozuyor.

MGallice ve MGallice-2D çözümlerinde toplam entalpinin bütün alanda makine hassasiyetine kadar sabit kaldığı belirtiliyor.

İkinci mertebe çözümün etkisi

Yarım silindir problemi Gallice-2D ve MGallice-2D şemalarının ikinci mertebe MUSCL uzantılarıyla yeniden hesaplanmıştır. Sayfa 11’deki Şekil 9 yoğunluk ve toplam entalpi alanlarını yan yana gösteriyor.

İkinci mertebe yeniden yapılandırma klasik Gallice-2D şemasındaki toplam enerji yayınımını azaltmıştır. Bu nedenle Gallice-2D ve MGallice-2D yoğunluk alanları birinci mertebe sonuçlarına göre birbirine daha çok yaklaşmıştır.

Bununla birlikte klasik Gallice-2D toplam entalpiyi tam olarak korumamaya devam ederken MGallice-2D alanı yine giriş değerinde sabit tutmuştur. Daha yüksek mertebe hata miktarını azaltmış fakat korunum özelliğinin yerini almamıştır.

Çalışmanın temel mühendislik sonucu

Önerilen değişiklik yalnızca grafiklerde daha düzgün bir toplam entalpi alanı üretmemiştir. Mach 20 testinde şok arkasındaki enerji dağılımını düzelterek durma bölgesinin yoğunluk alanını da değiştirmiştir.

Bu sonuç üç ayrı sayısal gerekliliği ortaya koyuyor:

• Şema güçlü şokları kararlı biçimde yakalamalıdır.
• Çok boyutlu şok kararsızlıklarına karşı yeterli direnç göstermelidir.
• Kararlı çözümün toplam entalpi gibi fiziksel değişmezlerini ayrık düzeyde korumalıdır.

MGallice-2D, incelenen Mach 20 testinde çok boyutlu şok kararlılığı ile toplam entalpi korunumunu birlikte sağlayan tek yapı olmuştur.

Türkiye açısından olası kullanım alanları

Yöntem, yüksek hızlı hava araçları ve atmosferik giriş problemleri için geliştirilen yerli veya akademik CFD yazılımlarında sayısal akı seçiminin değerlendirilmesine aktarılabilir. Özellikle aşağıdaki araştırma alanlarıyla ilişkilidir:

• Hipersonik ve süpersonik araç aerodinamiği,
• Roket ve yüksek hızlı uçuş araçlarının dış akış analizleri,
• Yay şoku ve durma noktası hesaplamaları,
• Yapısal olmayan ağlarda sonlu hacim çözücüleri,
• Şok kararsızlığına dayanıklı çok boyutlu Riemann çözücüleri,
• Navier–Stokes çözümlerinde sınır tabakası ve duvar ısı akısı hesapları.

Ancak araştırmanın sonuçları doğrudan bir aracın yüzey sıcaklığı, ısı kalkanı kalınlığı veya aerodinamik kuvvetleri olarak kullanılamaz. Bu büyüklükler için viskozite, ısı iletimi, gerçek gaz etkileri, kimyasal reaksiyonlar, türbülans ve katı malzeme davranışı ayrıca modellenmelidir.

Çalışmanın gösterdiği sonuçlar

• Gallice tipi Simple Riemann çözücüleri enerji ara durumları değiştirilerek toplam entalpi koruyan hale getirilebilir.
• Değişiklik tek ve çok boyutlu çözücülere uygulanabilir.
• Mach 2 üçgen çıkıntılı kanal testinde değiştirilmiş şemalar toplam entalpiyi bütün alanda korumuştur.
• Mach 20 yarım silindir testinde çok boyutlu şemalar karbunkül kararsızlığını önlemiştir.
• MGallice-2D, toplam entalpi kaybını ortadan kaldırarak durma bölgesindeki yoğunluk dağılımını düzeltmiştir.
• İkinci mertebe yeniden yapılandırma klasik şemanın entalpi hatasını azaltmış fakat tamamen ortadan kaldırmamıştır.

Çalışmanın göstermediği sonuçlar

• Önerilen çözücülerin bütün hipersonik geometrilerde en iyi yöntem olduğu gösterilmemiştir.
• Entropi kararlılığı matematiksel olarak kanıtlanmamıştır.
• Viskoz Navier–Stokes denklemleri üzerinde doğrulama yapılmamıştır.
• Duvar ısı akısı veya sınır tabakası doğruluğu doğrudan ölçülmemiştir.
• Gerçek gaz, kimyasal reaksiyon veya termokimyasal denge dışı etkiler incelenmemiştir.
• Üç boyutlu geometri veya üç boyutlu Riemann çözücüsü sınanmamıştır.
• Deneysel rüzgâr tüneli verileriyle karşılaştırma yapılmamıştır.
• Hesaplama süresi ve işlemci maliyeti açısından ayrıntılı performans karşılaştırması verilmemiştir.

Çalışmanın Yöntemi ve Bulguları

Araştırma tasarımı

Çalışma, kararlı sıkıştırılabilir Euler denklemleri için geliştirilen tek ve çok boyutlu Godunov tipi Simple Riemann çözücülerinin matematiksel olarak değiştirilmesine ve iki sayısal test problemiyle karşılaştırılmasına dayanan hesaplamalı mühendislik araştırmasıdır.

Matematiksel yöntem

Test sonuçlarının karşılaştırması

Sayısal değerler

• Mach 2 kanal ağı: 200 × 100 dörtgen hücre.
• Mach 2 giriş toplam entalpisi: 6,3.
• Mach 2 çözümlerinde yoğunluk kalıntısı yaklaşık 4.000 iterasyonda 10⁻⁸ düzeyine indi.
• Mach 2 yoğunluk görsel aralığı: 0,63-1,77.
• Mach 2 entalpi görsel aralığı: 6,2-6,44.
• Mach 20 yarım silindir ağı: 32 × 128 dörtgen hücre.
• Mach 20 giriş toplam entalpisi: 283,5.
• Mach 20 yoğunluk görsel aralığı: 1-6,3.
• Mach 20 entalpi görsel aralığı: 240-330.
• Mach 20 metninde kalıntının üç-dört mertebe düştüğü belirtilmiştir.

Şekillerin mühendislik yorumu

• Sayfa 4, Şekil 1: Tek boyutlu yüzey akıları ile düğüm tabanlı çok boyutlu akı şablonunu karşılaştırıyor.
• Sayfa 6, Şekil 2: Lagrange çerçevesindeki dört durumlu Simple Riemann çözücüsünü ve dalga hızlarını gösteriyor.
• Sayfa 9, Şekil 3: Mach 2 kanal probleminde dört şemanın kalıntı eğrilerinin yaklaşık 4.000 iterasyonda yakınsadığını gösteriyor.
• Sayfa 10, Şekil 4: Mach 2 problemindeki yoğunluk dalga düzenlerinin dört şemada genel olarak benzer olduğunu gösteriyor.
• Sayfa 10, Şekil 5: Entalpi hatasının klasik şemalarda şoklardan aşağı akıma taşındığını, değiştirilmiş şemalarda ise kaybolduğunu gösteriyor.
• Sayfa 10, Şekil 6: Mach 20 testinin kalıntı geçmişini ve çok boyutlu şemaların daha düzenli yakınsama eğilimini gösteriyor.
• Sayfa 10, Şekil 7: Yay şoku ve yoğunluk alanını; tek boyutlu çözücülerde karbunkülü ve Gallice-2D ile MGallice-2D arasındaki durma noktası farkını gösteriyor.
• Sayfa 10, Şekil 8: Klasik şemalardaki şok kaynaklı entalpi eksilmesini ve değiştirilmiş şemaların tekdüze entalpi alanını gösteriyor.
• Sayfa 11, Şekil 9: İkinci mertebe yeniden yapılandırmanın klasik şemadaki hatayı azalttığını, MGallice-2D’nin ise toplam entalpiyi tam olarak korumaya devam ettiğini gösteriyor.

Güçlü yönler

• Önerilen değişiklik doğrudan fiziksel bir kararlı akış değişmezine dayanmaktadır.
• Tek ve çok boyutlu çözücüler aynı matematiksel çerçevede ele alınmıştır.
• Şema yalnızca analitik olarak sunulmamış, güçlü şok içeren iki problemde sınanmıştır.
• Klasik ve değiştirilmiş çözücülerin yoğunluk ile entalpi alanları ayrı ayrı karşılaştırılmıştır.
• Mach 20 testi şok kararlılığı ile entalpi korunumunun birbirinden farklı gereklilikler olduğunu göstermiştir.
• Birinci ve ikinci mertebe çözümler karşılaştırılmıştır.
• MGallice-2D’nin yoğunluk alanındaki iyileşmesi, korunum özelliğinin akış çözümünü değiştirebildiğini göstermiştir.

Kritik sınırlılıklar

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Yalnızca viskozitesiz Euler denklemleri çözülmüştür.
• Viskozite, ısı iletimi ve doğrudan duvar ısı akısı bulunmamaktadır.
• Gerçek gaz ve yüksek sıcaklık kimyası modellenmemiştir.
• Entropi kararlılığı için matematiksel garanti sunulmamıştır.
• Çok boyutlu çözücüde pozitifliği garanti eden açık eğim koşulu türetilememiştir; tek boyutlu çözücüdeki eğimler kullanılmıştır.
• Yalnızca iki sayısal test problemi sunulmuştur.
• Üç boyutlu çözücü veya üç boyutlu araç geometrisi incelenmemiştir.
• Fiziksel deney veya başka bir bağımsız sayısal referansla nicel hata analizi yapılmamıştır.
• Mach 2 probleminde karşılaştırma ağırlıklı olarak alan görsellerine dayanmaktadır.
• Mach 20 probleminde metindeki 10⁵ iterasyon ifadesi ile kalıntı grafiğinin yaklaşık 11.000 iterasyonluk ekseni arasında tutarsızlık vardır.
• İşlemci süresi, bellek kullanımı ve klasik çözücülere göre ek hesaplama maliyeti raporlanmamıştır.
• MGallice-2D’nin duvar ısı akısı veya sınır tabakası belirleme performansı henüz sınanmamıştır.

Kaynak ve Yöntem Notu

Özgün çalışma: “Enthalpy preserving Simple Riemann solvers for steady hypersonic flows.”

Yazar: Lucas Tallois.

Sorumlu yazar: Lucas Tallois. Çalışma tek yazarlıdır ve eş katkı beyanı bulunmamaktadır.

Kurum: CEA-CESTA, 15 avenue des Sablières, CS 60001, 33116 Le Barp Cedex, Fransa.

Kaynak türü: Özgün hesaplamalı akışkanlar dinamiği araştırması ve preprint.

Preprint tarihi: 12 Haziran 2026.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Sonuçlar nihai ve hakemli bir dergi makalesinin bulguları olarak değerlendirilmemelidir.

Dergi durumu: Metinde “Preprint submitted to Elsevier” ifadesi yer almaktadır. Belirli dergi adı, kabul tarihi veya yayımlanmış nihai sürüm bilgisi bulunmamaktadır.

Platform: SSRN.

DOI:10.2139/ssrn.6948779

Resmî bağlantı:SSRN çalışma kaydı

Finansman: Çalışma metninde finansman beyanı bulunmamaktadır.

Çıkar çatışması: Çalışma metninde çıkar çatışması beyanı bulunmamaktadır.

Veri ve kod erişimi: Açık veri deposu, sayısal ağ dosyaları veya kaynak kod bağlantısı verilmemiştir.

Bu içerikteki bilimsel yöntem, denklemler, sayısal şemalar, test koşulları, grafik yorumları, bulgular ve sınırlılıklar incelenen çalışmaya dayanmaktadır. Dış kaynaklar yalnızca başlık, yazar, kurum, platform ve DOI bilgilerinin bibliyografik doğrulanması amacıyla kullanılmış; çalışmada bulunmayan bilimsel sonuç eklenmemiştir.