Bu çalışma, şok dalgaları, kayma tabakaları ve şok-dalga/sınır-tabaka etkileşimleri gibi çok güçlü değişimlerin bulunduğu hipersonik akışları hesaplamak için Haar dalgacık eşyerleştirme yöntemini geliştirmiştir. Yöntem, üç boyutlu sıkıştırılabilir Navier–Stokes denklemlerine uygulanmış ve üçüncü mertebeden Runge–Kutta zaman integrasyonu kullanılmıştır. Araştırmacılar yöntemi uzun süre taşınan izentropik bir girdap, Mach 7,5 karışım tabakası, Mach 6 üç boyutlu sınır tabakası geçişi ve Mach 6 sıkıştırma köşesi üzerinden sınamıştır. Sonuçlar; yöntemin süreksizlik çevresinde Gibbs tipi yapay salınımlar oluşturmadan keskin değişimleri temsil edebildiğini, küçük genlikli kararsızlıkları örtük sayısal sönümle bastırmadığını ve seçilen hipersonik örneklerde girdap oluşumu, akış geçişi ve duvara yakın termal gradyanları yeniden üretebildiğini göstermektedir. Bununla birlikte değerlendirmeler ağırlıklı olarak nitel akış yapısı karşılaştırmalarına dayanmaktadır; deneysel ölçüm, ayrıntılı hata normu ve güncel şok yakalama yöntemleriyle eşit maliyetli nicel kıyaslama sunulmamıştır.
Haar dalgacıklarının temel avantajı, basamak biçimli veya keskin değişen fonksiyonlara doğal olarak uyum sağlayan parçalı sabit bir tabana sahip olmasıdır. Klasik küresel spektral tabanlar bir süreksizlik çevresinde gerçek çözümde bulunmayan yüksek frekanslı dalgalanmalar oluşturabilir. Haar tabanı ise çalışmadaki Heaviside basamak fonksiyonu örneğinde birinci ve ikinci türevleri Gibbs salınımları olmadan temsil etmiştir. Yöntem güçlü gradyan bulunmayan bölgelerde örtük sayısal yayınım eklememekte; şok etkileşimleri gibi kararlılık gerektiren bölgelerde ise ağ inceldikçe azalan açık yapay viskozite kullanmaktadır.
Türkiye açısından değerlendirme: Yöntem; hipersonik hava araçları, yeniden giriş sistemleri, yüksek hızlı aerotermodinamik, sınır tabakası geçişi, şok-dalga/yüzey etkileşimleri ve termal koruma sistemi araştırmalarında Türkiye’de geliştirilen hesaplamalı akışkanlar dinamiği altyapılarına uyarlanabilir. Üniversiteler, araştırma merkezleri ve havacılık-uzay kuruluşlarında kullanılabilmesi için algoritmanın yerel çözücülere uygulanması, açık referans testleriyle yeniden üretilmesi, rüzgâr tüneli veya yayımlanmış deney verilerine karşı doğrulanması, paralel hesaplama performansının ölçülmesi ve WENO, TENO, sonlu hacim ya da süreksiz Galerkin yöntemleriyle aynı ağ ve hata düzeyinde karşılaştırılması gerekir. Çalışma Türkiye’de geliştirilmiş bir araç geometrisi, yerli uçuş deneyi, malzeme sistemi veya operasyonel tasarım incelemediğinden Türkiye’ye özgü performans, güvenlik ya da maliyet sonucu doğrudan çıkarılamaz.
Hipersonik akışların sayısal çözümü neden zordur?
Hipersonik akış alanlarında birbirinden çok farklı uzaysal ölçekler aynı anda bulunabilir. Bir aracın toplam uzunluğu metrelerle ifade edilirken şok tabakasındaki değişim mikrometre ölçeğinde gerçekleşebilir. Çalışmada tipik bir şok genişliğinin yaklaşık 7–9 mikrometre düzeyinde olabileceği ve termodinamik değişkenlerin bu çok dar bölgede ani bir sıçrama gösterdiği belirtilmektedir.
Aynı akış alanında aşağıdaki yapılar birlikte bulunabilir:
- Şok dalgaları,
- Sıkıştırma ve genişleme bölgeleri,
- Kayma tabakaları,
- Şok-şok etkileşimleri,
- Şok-dalga/sınır-tabaka etkileşimleri,
- Küçük genlikli kararsızlık dalgaları,
- Laminer-türbülanslı geçiş yapıları,
- Duvara yakın çok güçlü sıcaklık ve yoğunluk gradyanları.
Bir sayısal yöntemin şoku kararlı biçimde çözebilmesi için keskin değişimler çevresindeki yapay salınımları bastırması gerekir. Ancak fazla sayısal yayınım, geçişi başlatan çok küçük genlikli dalgaları da sönümleyebilir. Bu nedenle hipersonik geçiş hesaplarında iki gereksinim birbiriyle yarışır: şok çevresinde kararlılık ve düzgün bölgelerde düşük sönüm.
Mevcut yöntemlerin temel güçlükleri
Sonlu fark ailesindeki ENO, WENO ve TENO yöntemleri şoklu akışlarda yaygın biçimde kullanılmaktadır. Bu yöntemler çözümün düzgünlüğüne göre hesaplama kalıbını veya ağırlıkları değiştirerek süreksizlik çevresindeki salınımları sınırlar. Kompakt sonlu fark yöntemleri ise düzgün bölgelerde spektrale yakın çözünürlük sağlayabilir.
Sonlu eleman ve süreksiz Galerkin yöntemleri karmaşık geometrilere ve düzensiz ağlara daha kolay uygulanabilir. Ancak yüksek mertebeli elemanlarda şok yakalama için sensöre bağlı yapay viskozite veya başka kararlılık mekanizmaları gerekir. Ayrıca yüksek Mach sayılı geçiş dalgalarını çözmek için her elemanda çok sayıda serbestlik derecesi gerekebilir.
Klasik spektral yöntemler, düzgün fonksiyonlarda az sayıda serbestlik derecesiyle yüksek doğruluk sağlayabilir. Fakat küresel trigonometrik veya polinom tabanı, süreksizlik çevresinde Gibbs olgusu nedeniyle yapay dalgalanmalar oluşturur. Güçlü filtreleme veya yapay yayınım bu dalgalanmaları azaltabilse de küçük fiziksel kararsızlıkların kaybolmasına yol açabilir.
Dalgacık yaklaşımı neyi değiştirmektedir?
Dalgacık yöntemleri tek bir küresel düzgünlük varsaymak yerine, bir ana dalgacığın farklı ölçek ve konumlara taşınmış biçimlerini kullanır. Böylece hem uzaysal konum hem de ölçek bilgisi korunabilir. Bir çözümdeki yerel yapı, kullanılan dalgacık tabanına benziyorsa az sayıda katsayıyla temsil edilebilir.
Haar dalgacığı, en basit dalgacık ailesidir. Temel biçimi belirli bir aralıkta +1, bitişik aralıkta −1 ve diğer bölgelerde sıfırdır. Bu parçalı sabit ve süreksiz yapı, şok benzeri keskin geçişlere doğal olarak uygundur.
Çalışmanın özgün katkıları üç başlık altında sunulmuştur:
- Haar dalgacık eşyerleştirme yönteminin üç boyutlu sıkıştırılabilir Navier–Stokes denklemleri için tam matematiksel formülasyonunun geliştirilmesi,
- Yöntemin örtük sayısal sönüm oluşturmadığının uzun süreli girdap testiyle gösterilmesi,
- Mach 7,5 kayma tabakası, Mach 6 sınır tabakası geçişi ve Mach 6 sıkıştırma köşesi gibi hipersonik akışlara uygulanması.
Gibbs olgusu nedir?
Gibbs olgusu, süreksiz bir fonksiyonun düzgün ve küresel tabanlarla yaklaşık olarak temsil edilmesi sırasında sıçrama noktasının çevresinde oluşan yapay salınımlardır. Ağ veya terim sayısı artırıldığında salınımların kapladığı bölge daralabilir; ancak tepe genliği tamamen ortadan kalkmayabilir.
Çalışmada şok boyunca ani değişen bir termodinamik değişkeni temsil etmek için Heaviside basamak fonksiyonu kullanılmıştır:
\[ H_s(x)= \begin{cases} 1, & x\geq 0{,}5,\\ 0, & x<0{,}5. \end{cases} \]
Fonksiyonun ikinci türevi Haar dalgacık tabanında:
\[ \frac{d^2H_s}{dx^2}=\sum_{j=1}^{2M}a_jh_j(x) \]
biçiminde ifade edilmiştir. Burada \(h_j(x)\) Haar taban fonksiyonlarını, \(a_j\) ise dalgacık katsayılarını göstermektedir.
İki kez integrasyon yapıldığında:
\[ H_s=Q\mathbf{a}+Ax+B \]
ifadesi elde edilir. Sınır koşulları kullanıldığında katsayı vektörü:
\[ \mathbf{a}=\left[Q-xQ\big|_{x=1}\right]^{-1}(H_s-x) \]
olarak hesaplanır.
Şekil 1’de Haar yöntemiyle Chebyshev eşyerleştirme yöntemi karşılaştırılmıştır. Chebyshev çözümü süreksizliğin iki tarafında ve özellikle ikinci türevde çok büyük yapay dalgalanmalar üretirken Haar yöntemi bu dalgalanmaları göstermemektedir. Çalışmanın hipersonik akışlar açısından dayandığı temel matematiksel avantaj budur.
Haar dalgacık tabanı nasıl tanımlanmıştır?
\(x\in[0,1]\) aralığında bir Haar dalgacığı:
\[ h_i(x)= \begin{cases} 1, & x\in[\xi_1,\xi_2),\\ -1, & x\in[\xi_2,\xi_3),\\ 0, & \text{diğer durumlarda} \end{cases} \]
biçimindedir. Alt aralık sınırları:
\[ \xi_1=\frac{k}{m},\qquad \xi_2=\frac{k+0{,}5}{m},\qquad \xi_3=\frac{k+1}{m} \]
olarak verilir. \(m=2^j\) genişleme seviyesini, \(k\) ise dalgacığın uzaydaki yerini belirler. Çözünürlük seviyesi \(J\) yükseldikçe daha küçük ölçekli Haar fonksiyonları eklenir ve temsil edilen toplam dalgacık sayısı artar.
Şekil B.8, ilk üç Haar dalgacığını göstermektedir. İlk fonksiyon bütün aralıkta sabitken sonraki fonksiyonlar daha küçük aralıklarda +1 ve −1 değerleri alır. Bu hiyerarşik yapı, büyük ölçekli ve küçük ölçekli özelliklerin aynı tabanda temsil edilmesini sağlar.
En yüksek türev neden doğrudan dalgacıklarla temsil edilmiştir?
Haar fonksiyonları süreksiz olduğundan klasik anlamda türevlerini almak uygun değildir. Çalışmada bu sorun ters yönde çözülmüştür: Diferansiyel denklemdeki en yüksek mertebeli türev Haar tabanında yazılmış, daha sonra bu ifade integrasyon yoluyla asıl fonksiyona dönüştürülmüştür.
Haar fonksiyonunun \(r\) kez integrali:
\[ h_{r,i}(x)=\int_0^x\int_0^x\cdots\int_0^x h_i(t)\,dt^r \]
olarak tanımlanmıştır. Böylece süreksiz tabanın doğrudan türevini almak yerine, analitik olarak tanımlı integral matrisleri kullanılmıştır.
Üç boyutlu Navier–Stokes denklemleri
Yöntem, kalorik olarak mükemmel bir gaz için üç boyutlu sıkıştırılabilir Navier–Stokes denklemlerine uygulanmıştır:
\[ \frac{\partial\mathbf{w}}{\partial t} +\frac{\partial\mathbf{F}_i}{\partial x} +\frac{\partial\mathbf{G}_i}{\partial y} +\frac{\partial\mathbf{H}_i}{\partial z} = \frac{\partial\mathbf{F}_v}{\partial x} +\frac{\partial\mathbf{G}_v}{\partial y} +\frac{\partial\mathbf{H}_v}{\partial z}. \]
Korunan değişkenler vektörü:
\[ \mathbf{w}=(\rho,\rho u,\rho v,\rho w,E)^T \]
olarak verilmiştir. Burada \(\rho\) yoğunluk, \(u\), \(v\) ve \(w\) Kartezyen hız bileşenleri, \(E\) ise toplam enerji yoğunluğudur.
Basınç:
\[ p=(\gamma-1)\left[E-\frac{1}{2}\rho(u^2+v^2+w^2)\right] \]
ilişkisiyle elde edilir.
Viskoz gerilme tensörü:
\[ \sigma_{ij}=\mu\left( \frac{\partial u_i}{\partial x_j} +\frac{\partial u_j}{\partial x_i} \right) -\frac{2}{3}\mu\delta_{ij} \frac{\partial u_k}{\partial x_k} \]
ve ısı akısı:
\[ q_i=-k\frac{\partial T}{\partial x_i} \]
olarak tanımlanmıştır.
Isıl iletkenlik ve sıcaklık:
\[ k=\frac{C_p\mu}{Pr},\qquad T=\frac{p}{R\rho} \]
ilişkileriyle hesaplanmıştır. Viskozite için sıcaklığa bağlı kuvvet yasası kullanılmıştır:
\[ \frac{\mu}{\mu_{ref}}= \left(\frac{T}{T_{ref}}\right)^\omega \]
ve kaynakta \(\mu_{ref}=1{,}716\times10^{-5}\), \(T_{ref}=273\) ve \(\omega=0{,}7\) değerleri verilmiştir. Viskozite birimi kaynak denkleminde ayrıca yazılmamıştır.
Uzaysal ayrıklaştırma nasıl yapılmıştır?
Hesaplama \((\xi,\eta,\zeta)\) koordinatlarında yürütülmüş ve üç yön için bağımsız çözünürlük seviyeleri \(I\), \(J\) ve \(K\) kullanılmıştır. Her zaman adımında korunan değişkenlerin Haar katsayıları üç yönde ayrı ayrı hesaplanmaktadır.
\(\xi\) yönündeki katsayı vektörü örnek olarak:
\[ \mathbf{a}_{\xi} =MQ_{\xi}^{-1} \left[ \mathbf{w}-\mathbf{w}_L -\xi(\mathbf{w}_R-\mathbf{w}_L) \right] \]
biçiminde verilmiştir. Burada \(\mathbf{w}_L\) ve \(\mathbf{w}_R\), ilgili yöndeki sol ve sağ sınır değerleridir.
Birinci türev:
\[ \frac{\partial\mathbf{w}}{\partial\xi} =MP_\xi\mathbf{a}_\xi +(\mathbf{w}_R-\mathbf{w}_L) \]
ikinci türev ise:
\[ \frac{\partial^2\mathbf{w}}{\partial\xi^2} =MH_\xi\mathbf{a}_\xi \]
olarak hesaplanmıştır. Hesaplama koordinatlarındaki türevler daha sonra kullanılan ağ dönüşümü aracılığıyla fiziksel \(x\), \(y\) ve \(z\) türevlerine dönüştürülmüştür.
Zaman integrasyonu
Çalışmada üçüncü mertebeden Runge–Kutta yöntemi kullanılmıştır:
\[ \mathbf{w}_1= \Delta t\frac{\partial\mathbf{w}_0}{\partial t} +\mathbf{w}_0 \]
\[ \mathbf{w}_2= \frac{1}{4} \left( \Delta t\frac{\partial\mathbf{w}_1}{\partial t} +\mathbf{w}_1 \right) +\frac{3}{4}\mathbf{w}_0 \]
\[ \mathbf{w}_3= \frac{2}{3} \left( \Delta t\frac{\partial\mathbf{w}_2}{\partial t} +\mathbf{w}_2 \right) +\frac{1}{3}\mathbf{w}_0. \]
\(\mathbf{w}_0\), önceki zaman adımındaki çözümü temsil etmektedir. Her doğrulama örneğinde zaman adımı, belirtilen CFL sayısına göre belirlenmiştir.
Yöntem gerçekten hiç yapay yayınım kullanmıyor mu?
Çalışmada iki ayrı kavramın ayrılması gerekir:
- Haar ayrıklaştırmasının kendi yapısından kaynaklanan örtük sayısal yayınım bulunmadığı ileri sürülmektedir.
- Şok veya güçlü gradyan içeren hesaplarda kararlılık için açık yapay viskozite eklenmektedir.
Açık yapay viskozitenin katsayısı ağ çözünürlüğüyle ters orantılı seçilmiştir. Bu nedenle ağ inceltildikçe yapay viskozitenin azalması amaçlanmaktadır. Ancak kaynakta bütün testler için kullanılan katsayıların ayrıntılı listesi, sensör biçimi veya çözüm kalitesine duyarlılık analizi verilmemiştir. Dolayısıyla “hiç sayısal yayınım yoktur” ifadesi yalnızca yöntemin örtük davranışı için geçerlidir; güçlü gradyan testlerinde açık kararlılık terimi kullanılmaktadır.
H, P ve Q matrislerinin önemi
Haar yaklaşımında \(H\) taban değerlerini, \(P\) birinci integralleri, \(Q\) ise ikinci integralleri temsil etmektedir. Navier–Stokes denklemlerindeki en yüksek uzaysal türev ikinci mertebe olduğundan \(Q\) matrisi gerekli en düşük integral seviyesidir.
Şekil 2 ve Şekil 3, bu matrislerin ve terslerinin sıfır olmayan elemanlarını göstermektedir. \(H\) ve \(P\) ileri matrisleri oldukça seyrektir. \(Q\) ve özellikle katsayı hesaplamasında kullanılan ters matrisler daha yoğundur. \(H\) matrisi ortogonaldir; integrasyon sonrasında bu özellik kaybolmaktadır.
Mevcut uygulama, doğruluğu göstermek amacıyla yoğun matris-vektör çarpımı kullanmıştır. Bu uygulamanın her uzaysal yöndeki hesaplama maliyeti yaklaşık:
\[ O(N_i^2) \]
ölçeğindedir ve paralel verimliliği sınırlıdır. Araştırmacılar, hızlı dalgacık dönüşümlerinin bu maliyeti teorik olarak:
\[ O(N_i) \]
mertebesine düşürebileceğini belirtmektedir. Ancak hızlı algoritma bu çalışmada uygulanıp ölçülmemiştir; gelecekteki geliştirme yönü olarak sunulmuştur.
Birinci test: Uzun süreli izentropik girdap
İlk test, yöntemin örtük sayısal yayınımını ölçmek amacıyla durağan bir izentropik girdabın uzun süre korunup korunamadığını incelemiştir. Düzgün ve tam çözümü bilinen bu problemde girdabın zamanla yayılması veya tepe değerinin azalması doğrudan sayısal yayınım göstergesi olarak yorumlanabilir.
Ortalama akış:
\[ \rho=1,\qquad p=1,\qquad u=U_\infty,\qquad v=V_\infty \]
olarak tanımlanmış ve durağan durum için \(U_\infty=V_\infty=0\) alınmıştır.
Girdap hız pertürbasyonları:
\[ (\delta u,\delta v)= \frac{\beta}{2\pi} e^{0{,}5(1-r^2)} (-\bar{y},\bar{x}) \]
ve sıcaklık pertürbasyonu:
\[ \delta T= -\frac{(\gamma-1)\beta^2} {8\gamma\pi^2} e^{1-r^2} \]
olarak verilmiştir. Entropi pertürbasyonu sıfırdır.
| Girdap testi parametresi | Değer |
|---|---|
| Hesaplama alanı | \([0,50]\times[-5,5]\) |
| Girdap merkezi | \((25,0)\) |
| Girdap şiddeti | \(\beta=5\) |
| Çözünürlük seviyeleri | \(I=7,\ J=5\) |
| Akış yönü ağ aralığı | \(\Delta x=0{,}1953\) |
| Düşey ağ aralığı | \(\Delta y=0{,}1562\) |
| CFL sayısı | \(10^{-2}\) |
| Karşılaştırma zamanları | \(t=0,\ 50,\ 100\) |
Şekil 4’te üç zamandaki merkez çizgisi yoğunluk profilleri neredeyse üst üste gelmektedir. Araştırmacılar bunu yöntemin 100 zaman birimi boyunca girdabı örtük yayınımla bozmadığının kanıtı olarak yorumlamıştır.
Bu sonuç önemli olmakla birlikte çalışmada \(L_1\), \(L_2\) veya \(L_\infty\) hata normları verilmemiştir. Eğrilerin görsel olarak çakışması fiziksel yapının korunduğunu göstermekte, fakat sayısal hatanın kesin büyüklüğünü vermemektedir.
İkinci test: Mach 7,5 hipersonik kayma tabakası
İkinci örnek, farklı hızlara ve yoğunluklara sahip iki akışın arasında gelişen iki boyutlu bir karışım tabakasıdır. Giriş hız profili hiperbolik tanjant biçimindedir:
\[ u(y)= \frac{U_1+U_2}{2} +\frac{U_1-U_2}{2} \tanh\left(\frac{2y}{\delta_\omega}\right). \]
| Kayma tabakası parametresi | Değer |
|---|---|
| Üst akım hızı | \(U_1=4\) |
| Alt akım hızı | \(U_2=2\) |
| Girdap kalınlığı | \(\delta_\omega=1\) |
| Serbest akım basıncı | \(p_\infty=0{,}3327\) |
| Akım yoğunlukları | \(\rho_1=1{,}6374,\ \rho_2=0{,}3626\) |
| Mach sayıları | \(M_1=7{,}5,\ M_2=1{,}76\) |
| Konvektif Mach sayısı | 1,2 |
| Hesaplama alanı | \(0\leq x/\delta_\omega\leq200\), \(-20\leq y/\delta_\omega\leq20\) |
| Çözünürlük seviyeleri | \(I=7,\ J=5\) |
| CFL sayısı | 0,45 |
Geçiş, duvara dik hız pertürbasyonlarıyla başlatılmıştır:
\[ v'= \sum_{k=1}^{2} a_k \cos\left(\frac{2\pi kt}{T}+\phi_k\right) \exp\left(-\frac{y^2}{b}\right). \]
Dalga boyu \(\lambda=30\), iki modun genlikleri \(a_1=a_2=0{,}05\), fazları \(\phi_1=0\) ve \(\phi_2=\pi/2\), genişlik parametresi ise \(b=10\) olarak seçilmiştir.
Konvektif Mach sayısının 1,2 olması, kayma tabakası içinde küçük şok yapılarının oluşmasına elverişli bir rejim meydana getirmiştir. Bu testte açık yapay viskozite kullanılmış ve katsayı ağ çözünürlüğü yükseldikçe azaltılmıştır.
Şekil 5’te yoğunluk alanı boyunca Kelvin–Helmholtz kararsızlığının gelişimi görülmektedir:
- Başlangıç pertürbasyonları akış yönünde büyümektedir.
- Kayma tabakası dalgalanarak yuvarlanan girdaplara dönüşmektedir.
- Komşu girdaplar eşleşmekte ve daha büyük yapılar meydana getirmektedir.
- Karışım tabakasının kalınlığı akış yönünde artmaktadır.
Bu sonuç yöntemin şok benzeri keskin değişimlerle düşük genlikli fiziksel kararsızlıkları aynı hesapta temsil edebildiğini göstermektedir. Ancak çalışmada girdap büyüme oranı, karışım tabakası kalınlığı veya baskın frekans için referans çözümle nicel hata karşılaştırması verilmemiştir.
Üçüncü test: Mach 6 sınır tabakası geçişi
Üçüncü test, düz bir yüzey üzerindeki üç boyutlu hipersonik sınır tabakasının kararsızlık dalgalarından geçiş yapılarına dönüşmesini incelemektedir. Problem, daha önce yayımlanmış Mach 6 referans yapılandırmasına dayandırılmıştır.
| Sınır tabakası parametresi | Değer |
|---|---|
| Serbest akım Mach sayısı | 6 |
| Giriş Reynolds sayısı | \(Re_{\delta_0^*}=3000\) |
| Alan uzunluğu | \(1000\delta_0^*\) |
| Alan yüksekliği | \(37{,}5\delta_0^*\) |
| Açıklık genişliği | \(20\pi\delta_0^*\) |
| Duvar/serbest akım sıcaklık oranı | \(T_w=6{,}5T_\infty\) |
| Çözünürlük seviyeleri | \(I=9,\ J=5,\ K=6\) |
Kaynak, \(T_w/T_\infty=6{,}5\) durumunu hipersonik uçuş koşullarına göre güçlü biçimde soğutulmuş duvar olarak tanımlamaktadır. Burada “soğutulmuş” ifadesi serbest akım statik sıcaklığına değil, çok daha yüksek olabilen adyabatik veya geri kazanım sıcaklığına göre yorumlanmalıdır.
Geçişi uyarmak için yüzeyde emme-üfleme biçiminde hız zorlaması kullanılmıştır:
\[ v_w=f(x)g(z) \sum_{i=1}^{2} A_i\sin(\omega_it-\beta_iz). \]
Zorlama bölgesi:
\[ 15\leq\frac{x-x_0}{\delta_0^*}\leq20 \]
aralığındadır. Akış yönündeki Gauss zarfının merkezi \(17{,}5\delta_0^*+x_0\), genişliği \(0{,}75\delta_0^*\) olarak seçilmiştir. İki eğik modun:
- Genliği \(A=0{,}05U_\infty\),
- Frekansı \(\omega=0{,}9c_\infty/\delta_0^*\),
- Açıklık yönündeki dalga sayıları \(\beta=\pm0{,}3/\delta_0^*\)
olarak verilmiştir.
Şekil 6, duvara paralel düzlemde normalize edilmiş yoğunluk gradyanlarını ve Q ölçütü izoyüzeylerini göstermektedir. Görselde ilk bölümde düzenli ve akış yönüne yakın yapılar bulunurken ilerleyen bölgede girdap yapıları çoğalmakta, bükülmekte ve üç boyutlu karmaşık bir örgüye dönüşmektedir.
Araştırmacılar bu gelişimi şu mekanizmalarla ilişkilendirmiştir:
- Eğik kararsızlık dalgalarının büyümesi,
- Doğrusal olmayan mod etkileşimleri,
- Akış yönlü girdapların oluşumu,
- Üç boyutlu geçiş yapılarının ortaya çıkması,
- Bozulmanın erken evresinin başlaması.
Çalışma tam gelişmiş türbülansın istatistiksel özelliklerini veya yüzey ısı akısı aşımını nicel olarak incelememektedir. Gösterilen sonuç, yöntemin geçişin temel akış yapılarını üretebildiğine ilişkin nitel bir doğrulamadır.
Dördüncü test: Mach 6 sıkıştırma köşesi
Son test, düz bir levhanın ardından 8° eğimli bir sıkıştırma rampası bulunan üç boyutlu akış alanıdır. Rampanın başlangıcı küçük bir yarıçapla yuvarlatılmıştır. Bu yapılandırma, şok-dalga/sınır-tabaka etkileşiminin ve duvara yakın güçlü termal gradyanların birlikte bulunduğu klasik bir yüksek hızlı aerodinamik problemidir.
| Sıkıştırma köşesi parametresi | Değer |
|---|---|
| Serbest akım Mach sayısı | 6 |
| Rampa açısı | 8° |
| Toplam sıcaklık | 648 K |
| Serbest akım statik basıncı | 1145 Pa |
| Duvar/serbest akım sıcaklık oranı | 2,05 |
| Referans uzunluk | \(\delta_0=1\) mm |
| Akış yönü alanı | \(x/\delta_0\in[-60,140]\) |
| Duvara dik alan | \(y/\delta_0\in[0,30]\) |
| Açıklık yönü alanı | \(z/\delta_0\in[0,8\pi]\) |
| Çözünürlük seviyeleri | \(I=9,\ J=5,\ K=6\) |
Giriş koşulu ayrı bir düz levha öncül çözümünden üretilmiştir. Bu öncül alan, üç boyutlu bozuntuları girişte sürekli olarak sağlamaktadır. Başlangıçtaki aşırı geçici davranışı azaltmak için iç alan, basınçla uyumlu ve düzgün bir sınır tabakası benzeri profille başlatılmıştır.
Şekil 7’de simetri düzlemindeki sıcaklık alanı gösterilmektedir. Düz levha üzerinde duvara yakın bir termal tabaka gelişmektedir. Akış köşeye ve rampaya ulaştığında sıkıştırma nedeniyle duvar çevresindeki sıcaklık yükselmektedir. Rampanın devamında sıcak tabaka kalınlaşmakta ve kuvvetli sıcaklık gradyanları yüzeye yakın bölgede yoğunlaşmaktadır.
Bu test, yöntemin karmaşık geometride aşağıdaki nitel özellikleri üretebildiğini göstermektedir:
- Rampanın oluşturduğu sıkıştırma etkisi,
- Duvar yakınında sıcaklık yükselmesi,
- Termal sınır tabakasının bozulması,
- Köşe sonrasında kuvvetli gradyanların devam etmesi,
- Üç boyutlu giriş bozuntularının aşağı akışa taşınması.
Bununla birlikte çalışma yüzey basınç katsayısı, Stanton sayısı, ısı akısı, ayrılma uzunluğu veya şok konumu için deneysel ya da yüksek çözünürlüklü referansla nicel karşılaştırma sunmamaktadır.
Çalışmanın güçlü yönleri
- Haar dalgacık yöntemi üç boyutlu sıkıştırılabilir Navier–Stokes denklemlerine sistematik olarak uygulanmıştır.
- Süreksizliklerin Gibbs salınımları olmadan temsil edilebilmesi doğrudan Heaviside fonksiyonu üzerinden gösterilmiştir.
- Doğrulama problemleri düzgün ve durağan bir girdaptan üç boyutlu hipersonik şok/sınır-tabaka etkileşimine kadar kademeli olarak zorlaştırılmıştır.
- Uzun süreli girdap testi, düşük genlikli yapıların örtük yayınımla bastırılmadığını göstermektedir.
- Mach 7,5 kayma tabakasında Kelvin–Helmholtz girdap oluşumu ve eşleşmesi üretilmiştir.
- Mach 6 sınır tabakasında eğik modların büyümesi ve akış yönlü girdap yapıları elde edilmiştir.
- Mach 6 sıkıştırma köşesinde duvara yakın kuvvetli termal gradyanlar çözülmüştür.
- Haar matrislerinin hiyerarşik ve seyrek yapısı gelecekteki hızlandırma için açık bir algoritmik yol sunmaktadır.
- Yazar katkıları, çıkar çatışması, veri erişimi ve finansman bilgileri açıkça verilmiştir.
Çalışmanın sınırlılıkları
- Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
- Fiziksel rüzgâr tüneli veya uçuş deneyiyle doğrulama yapılmamıştır.
- Dört testin büyük bölümü nitel görsel akış yapısı değerlendirmesine dayanmaktadır.
- Hata normları, gözlenen yakınsama mertebesi ve kapsamlı ağ yakınsama tabloları verilmemiştir.
- WENO, TENO, kompakt fark veya süreksiz Galerkin yöntemleriyle aynı ağ ve zaman adımında doğrudan doğruluk-maliyet karşılaştırması yapılmamıştır.
- Güçlü gradyan testlerinde açık yapay viskozite kullanılmış, fakat bütün katsayılar ve duyarlılık sonuçları ayrıntılı biçimde raporlanmamıştır.
- Mevcut yoğun matris uygulaması \(O(N_i^2)\) maliyetli ve paralel verimliliği sınırlıdır.
- Önerilen hızlı \(O(N_i)\) algoritma uygulanmamış ve performansı ölçülmemiştir.
- Uyarlamalı ağ veya gerçek zamanlı çok çözünürlüklü etkinleştirme bu çalışmada uygulanmamıştır.
- Karmaşık endüstriyel geometriler veya düzensiz ağlar üzerinde doğrulama sunulmamıştır.
- Kimyasal tepkime, gerçek gaz etkileri, titreşimsel enerji ve yüksek sıcaklık gaz fiziği modele dâhil edilmemiştir.
- Veriler yalnızca talep üzerine sağlanacaktır; açık kod deposu ve yeniden üretim paketi verilmemiştir.
Çalışma neyi desteklemektedir?
- Haar tabanı basamak tipi süreksizliklerde Gibbs salınımlarını önleyebilir.
- Haar eşyerleştirme yöntemi düzgün akış yapılarını uzun süre düşük örtük sönümle koruyabilir.
- Yöntem şok veya güçlü gradyanlarla geçiş dalgalarının birlikte bulunduğu örnekleri kararlı biçimde hesaplayabilir.
- Seçilen çözünürlüklerde Kelvin–Helmholtz girdapları ve hipersonik sınır tabakası geçiş yapıları üretilebilir.
- Haar matrislerinin düzenli seyreklik yapısı hızlı algoritmalar için kullanılabilir.
- Yöntem hipersonik geçiş hesapları açısından araştırılmaya değer bir alternatif oluşturmaktadır.
Çalışma neyi kanıtlamamaktadır?
- Haar yönteminin WENO, TENO veya diğer güncel yöntemlerden genel olarak daha doğru olduğunu kanıtlamamaktadır.
- Yöntemin aynı doğrulukta daha düşük hesaplama maliyetine sahip olduğunu göstermemektedir.
- Her türlü şok şiddeti, geometri veya sınır koşulunda kararlı çalışacağını kanıtlamamaktadır.
- Hızlı dalgacık dönüşümü kullanıldığında gerçek uygulamanın \(O(N)\) hızına ulaşacağını deneysel olarak göstermemektedir.
- Uyarlamalı çözünürlüğün otomatik ve güvenilir biçimde çalıştığını göstermemektedir.
- Gerçek gaz, kimyasal reaksiyon veya plazma içeren hipersonik rejimleri doğrulamamaktadır.
- Belirli bir hava aracının aerodinamik, ısıl veya uçuş güvenliği performansını hesaplamamaktadır.
- Endüstriyel tasarım veya sertifikasyon için hazır bir çözücü sunmamaktadır.
Gelecekte hangi doğrulamalar gereklidir?
Yöntemin araştırma düzeyinden genel amaçlı hipersonik çözücü düzeyine ilerlemesi için aşağıdaki çalışmalar önem taşımaktadır:
- Bilinen analitik çözümlerde sistematik \(L_1\), \(L_2\) ve \(L_\infty\) hata ölçümü,
- Farklı çözünürlüklerde gözlenen yakınsama mertebesinin hesaplanması,
- Yapay viskozite katsayısının şok kalınlığı ve geçiş dalgalarına etkisinin incelenmesi,
- Aynı ağ ve CFL koşullarında WENO/TENO yöntemleriyle karşılaştırma,
- CPU/GPU süreleri, bellek tüketimi ve paralel ölçeklenmenin ölçülmesi,
- Hızlı Haar dönüşümlerinin uygulanması,
- Uyarlamalı çok çözünürlüklü ağ stratejisinin geliştirilmesi,
- Rüzgâr tüneli basıncı, ısı akısı ve geçiş konumu verilerine karşı doğrulama,
- Gerçek gaz ve yüksek sıcaklık termokimyasıyla genişletme,
- Karmaşık üç boyutlu araç geometrilerine uygulama.
Çalışmanın Yöntemi ve Bulguları
Yöntemsel çerçeve
| Bileşen | Çalışmada kullanılan yaklaşım |
|---|---|
| Denklem sistemi | Kalorik olarak mükemmel gaz için üç boyutlu sıkıştırılabilir Navier–Stokes denklemleri |
| Uzaysal yöntem | Haar dalgacık eşyerleştirme yöntemi |
| Temel temsil | En yüksek türevin Haar tabanında yazılması ve integrasyonla çözümün elde edilmesi |
| Zaman integrasyonu | Üçüncü mertebeden Runge–Kutta |
| Koordinatlar | Hesaplama alanında \((\xi,\eta,\zeta)\), dönüşümle fiziksel \((x,y,z)\) |
| Matrisler | Haar tabanı için H, birinci integral için P, ikinci integral için Q |
| Şok kararlılığı | Ağ çözünürlüğü arttıkça azalan açık yapay viskozite |
| Mevcut hesaplama maliyeti | Yoğun matris işlemleri nedeniyle yön başına yaklaşık \(O(N_i^2)\) |
| Önerilen gelecek maliyet | Hızlı dalgacık dönüşümüyle teorik \(O(N_i)\) |
| Doğrulama sayısı | Dört test |
| Deneysel doğrulama | Bulunmamaktadır |
| Açık kod | Bulunmamaktadır |
| Veri erişimi | Talep üzerine |
Dört testin karşılaştırmalı özeti
| Test | Akış rejimi | Amaç | Temel gözlem |
|---|---|---|---|
| İzentropik girdap | Durağan, düzgün ve viskozitesiz | Örtük sayısal sönümü sınamak | Yoğunluk profili \(t=0\), 50 ve 100’de görsel olarak korunmuştur |
| Karışım tabakası | Mach 7,5 ve Mach 1,76 akımlar | Hipersonik viskoz geçiş ve şokçukları incelemek | Kelvin–Helmholtz yuvarlanması ve girdap eşleşmesi yakalanmıştır |
| Sınır tabakası geçişi | Üç boyutlu Mach 6 | Eğik kararsızlıkların doğrusal olmayan bozulmasını incelemek | Akış yönlü girdaplar ve erken geçiş yapıları üretilmiştir |
| Sıkıştırma köşesi | Üç boyutlu Mach 6, 8° rampa | Şok/sınır-tabaka etkileşimi ve termal gradyanları incelemek | Rampada sıcaklık artışı ve duvar yakınında güçlü gradyanlar oluşmuştur |
Şekillerin bilimsel anlamı
| Şekil | Gösterilen içerik | Makaledeki işlevi |
|---|---|---|
| Şekil 1 | Heaviside fonksiyonunun birinci ve ikinci türevlerinde Haar ve Chebyshev karşılaştırması | Haar yönteminin Gibbs salınımı üretmediğini göstermek |
| Şekil 2 | H, P ve Q matrislerinin seyreklik desenleri | Haar tabanının hiyerarşik ve seyrek hesaplama yapısını göstermek |
| Şekil 3 | \(H^{-1}\), \(P^{-1}\) ve \(Q^{-1}\) matrislerinin desenleri | Katsayı hesaplamasında ters matrislerin daha yoğun olduğunu göstermek |
| Şekil 4 | Durağan girdabın üç zamandaki yoğunluk profili | Uzun süreli örtük yayınımın düşük olduğunu göstermek |
| Şekil 5 | Hipersonik karışım tabakasının yoğunluk alanı | Kelvin–Helmholtz yuvarlanması ve girdap eşleşmesini göstermek |
| Şekil 6 | Mach 6 sınır tabakasında yoğunluk gradyanı ve Q ölçütü yapıları | Akış yönlü girdapların ve üç boyutlu geçişin gelişimini göstermek |
| Şekil 7 | Mach 6 sıkıştırma köşesindeki sıcaklık alanı | Rampaya bağlı ısınma ve güçlü duvar yakınındaki gradyanları göstermek |
| Şekil B.8 | İlk üç Haar dalgacığı | Parçalı sabit ve çok ölçekli tabanın biçimini açıklamak |
Ana bulgular
- Haar tabanı, çalışmadaki süreksiz test fonksiyonunda Gibbs tipi sayısal dalgalanmalar oluşturmamıştır.
- Durağan izentropik girdap 100 zaman birimi boyunca görünür profil kaybı olmadan korunmuştur.
- Mach 7,5 akım içeren karışım tabakasında fiziksel kararsızlık büyümesi, girdap yuvarlanması ve girdap eşleşmesi hesaplanmıştır.
- Mach 6 sınır tabakasında eğik modların doğrusal olmayan etkileşimi ve akış yönlü girdapların oluşumu yeniden üretilmiştir.
- Mach 6 sıkıştırma rampasında duvara yakın sıcaklık artışı ve güçlü termal gradyanlar çözülmüştür.
- Yöntem güçlü gradyanlarda açık yapay viskoziteyle kararlı hâle getirilmiş, viskozite katsayısının ağ inceldikçe azalması sağlanmıştır.
- Mevcut yoğun matris uygulaması fiziksel doğruluğa odaklanmış, yüksek hesaplama performansı henüz gösterilmemiştir.
- Haar matrislerinin yapısı gelecekte hızlı dönüşüm ve uyarlamalı çözünürlük için uygun görülmüştür.
Bulguların dikkatli yorumu
Çalışma, Haar eşyerleştirme yaklaşımının hipersonik geçiş ve güçlü gradyan içeren akışlar için uygulanabilir olduğunu göstermektedir. “Uygulanabilir” ifadesi, seçilen testlerin kararlı biçimde çalıştırıldığı ve beklenen temel fiziksel yapıların üretildiği anlamındadır.
Sonuçlar yöntemin genel olarak üstün, daha hızlı veya daha doğru olduğunu henüz göstermemektedir. Bunun için eşit çözünürlükte güncel yöntemlerle nicel hata, hesaplama süresi ve bellek karşılaştırması gerekir. Ayrıca geçiş ve sıkıştırma köşesi sonuçlarının deneysel basınç, ısı akısı veya geçiş konumu ölçümleriyle sınanması gereklidir.
Kaynak ve Yöntem Notu
Çalışmanın tam özgün adı: A Haar wavelet collocation method for hypersonic flows with strong gradients
Yazarlar: Yu Zhang, Diego Exposito-Brioso ve Xuerui Mao.
Çalışmada kullanılan kısaltılmış adlar: Y. Zhang, D. Exposito ve X. Mao.
Yazar sıralaması: Kaynakta verildiği sırayla korunmuştur.
Eş birinci yazar: Eşit katkı veya eş birinci yazarlık beyanı yer almamaktadır.
Sorumlu yazarlar: Diego Exposito-Brioso ve Xuerui Mao.
İletişim adresleri: dieexbr17@gmail.com ve xmao@bit.edu.cn.
Kurum 1: School of Interdisciplinary Science, Beijing Institute of Technology, Beijing 100081, Çin.
Kurum 2: School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, Çin.
Kurum 3: Beijing Institute of Technology (Zhuhai), Zhuhai 519088, Çin.
Kurum 4: State Key Laboratory of Explosion Science and Safety Protection, Beijing Institute of Technology, Beijing 100081, Çin.
Resmî kaynak bağlantısı:SSRN resmî kayıt sayfası
Yayın platformu: SSRN.
Yayın yılı: 2026.
Dergi: Belirli bir hakemli dergi adı veya dergi kabulü bu sürümde yer almamaktadır.
Hakemli yayınevi: Hakemli bir dergi yayınevi bilgisi bu sürümde yer almamaktadır.
Kaynak türü: Üç boyutlu sıkıştırılabilir akışlar için sayısal yöntem geliştirme ve doğrulama çalışması niteliğinde preprint araştırma makalesi.
Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. Sonuçlar, preprint niteliği ve belirtilen doğrulama sınırlılıkları dikkate alınarak değerlendirilmelidir.
Yu Zhang’ın katkıları: Yöntem, yazılım, doğrulama, biçimsel analiz, araştırma, veri düzenleme, özgün taslak yazımı ve görselleştirme.
Diego Exposito-Brioso’nun katkıları: Kavramsallaştırma, yöntem, biçimsel analiz, inceleme ve düzenleme, danışmanlık ve finansman temini.
Xuerui Mao’nun katkıları: Kaynak sağlama, inceleme ve düzenleme, danışmanlık, proje yönetimi ve finansman temini.
Finansman: National Natural Science Foundation of China, Hibe No. W2433020.
Çıkar çatışması: Yazarlar çalışmayı etkileyebilecek bilinen bir mali çıkar veya kişisel ilişki bulunmadığını beyan etmiştir.
Veri erişimi: Verilerin talep üzerine sağlanacağı belirtilmiştir. Açık veri deposu bağlantısı verilmemiştir.
Kod erişimi: Açık kaynak kod, sürüm etiketi veya yeniden üretim paketi sunulmamıştır.
Yöntemsel sınır: Çalışma kalorik olarak mükemmel gaz varsayımına dayanmaktadır. Gerçek gaz, termokimyasal dengesizlik, kimyasal reaksiyon, iyonlaşma ve radyasyon modellenmemiştir. Hipersonik örneklerin doğrulaması ağırlıklı olarak beklenen akış yapılarının nitel olarak yeniden üretilmesine dayanmaktadır.
Hesaplama sınırı: Mevcut uygulama yoğun matris işlemleri nedeniyle yaklaşık \(O(N_i^2)\) maliyet taşımaktadır. Hızlı \(O(N_i)\) dalgacık algoritmaları ve uyarlamalı çok çözünürlüklü uygulama, gerçekleştirilmiş sonuç değil gelecek çalışma önerisidir.
Yapay viskozite uyarısı: Yöntem örtük sayısal sönüm oluşturmadığını ileri sürmekle birlikte, güçlü şok etkileşimlerinde açık yapay viskozite kullanmıştır. Bu nedenle sonuçlar “tamamen viskozitesiz sayısal çözüm” olarak yorumlanmamalıdır.
Bu Türkçe içerik yüklenen çalışmanın metni, denklemleri, şekilleri, ekleri ve sonuçları temel alınarak hazırlanmıştır. Çalışmada bulunmayan bir deney başarısı, uçuş performansı, araç tasarımı, endüstriyel üstünlük veya başka yöntemlere karşı kanıtlanmış genel üstünlük iddiası eklenmemiştir.
Bir yorum bırakın
E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir