Reaksiyon-difüzyon manifoldları yöntemi (Reaction-Diffusion Manifolds, REDIM), ayrıntılı yanma kimyasındaki çok sayıda tür ve reaksiyonun oluşturduğu yüksek boyutlu termokimyasal durum uzayını, kimyasal reaksiyonlarla moleküler taşınım arasındaki etkileşimi koruyarak daha düşük boyutlu bir manifold üzerinde temsil etmeyi amaçlayan bir indirgenmiş kimya yöntemidir. Bu çalışma REDIM yöntemini OpenFOAM tabanlı yeni bir CFD çözücüsü olan REDIMFoam içine uygulamakta ve çözücüyü 1D serbest yayılan metan/hava alevleri ile 2D eksensimetrik karşı akışlı metan/hava difüzyon alevleri üzerinde ayrıntılı kimya kullanan sayısal çözümlerle karşılaştırmaktadır.

Çalışmanın amacı yalnız REDIM'in bir alev profilini yaklaşık olarak üretip üretemediğini görmek değildir. Araştırmacılar genelleştirilmiş ve fiziksel koordinat temsillerini, farklı ağ çözünürlüklerini, iki farklı progress-variable seçimini, iki ayrıntılı reaksiyon mekanizmasını, farklı giriş hızlarını ve kararlı durumdan sönmeye geçişi ayrı ayrı incelemiştir. Karşılaştırma referansı olarak ayrıntılı sonlu hızlı kimyayı çözen OpenFOAM tabanlı EBIdnsFoam kullanılmıştır.

1D serbest alev testinde ayrıntılı EBIdnsFoam hesabı 4,49 s CPU süresi kullanırken REDIMFoam genelleştirilmiş koordinatlarda 0,51 s, fiziksel koordinatlarda 0,55 s kullanmıştır. Böylece REDIM hesap süresi ayrıntılı mekanizma süresinin yaklaşık %11–12'sine düşmüştür. Bu, aynı test koşullarında yaklaşık %88–89 daha düşük CPU süresine karşılık gelir. Kaynak sonuç bölümünde genel avantajı yaklaşık bir büyüklük mertebesi hesaplama maliyeti azalması olarak özetlemektedir.

2D karşı akışlı alevlerde sıcaklık ve ana türlerin kütle kesirleri ayrıntılı çözümlerle çok yakın sonuçlar vermiştir. GRI3.0 mekanizmasında 1,0 m/s giriş hızında ayrıntılı çözümün tepe sıcaklığı 2016 K, REDIM sonucu 2015 K'dir. SanDiego-2014 için karşılık gelen değerler 2031 K ve 2032 K olarak raporlanmış ve REDIM kaynaklı tepe sıcaklığı hatasının %0,1'in altında olduğu belirtilmiştir. Küçük radikal türlerde, özellikle yüksek global gerinim hızlarında OH gibi türler için farklar daha görünür hale gelebilmektedir.

Çalışmanın önemli sonucu, REDIM'in yalnız manifold oluşturulurken açıkça verilen kararlı alev durumlarını değil, zamanla gerinim oranının artırılması sonucu gerçekleşen sönme sürecini de izleyebilmesidir. Bununla birlikte bu bulgular sayısal doğrulama sonuçlarıdır. Çalışma gerçek brülör üzerinde deney gerçekleştirmemekte, endüstriyel motor veya gaz türbini performansı doğrulamamakta ve bu makalede doğrudan türbülanslı alev simülasyonu sunmamaktadır.

REDIM Nedir ve Ayrıntılı Yanma Kimyasını Nasıl Azaltır?

REDIM, ayrıntılı kimyasal kinetiğin yüksek boyutlu durum uzayındaki çözümünün düşük boyutlu bir manifoldun yakınında bulunduğu varsayımından yararlanarak, her kimyasal tür için ayrı taşınım denklemi çözmek yerine sınırlı sayıdaki manifold koordinatını çözen ve sıcaklık ile tür bileşimlerini önceden oluşturulmuş REDIM tablosundan geri kazanan bir model indirgeme yöntemidir.

Ayrıntılı yanma kimyası neden pahalıdır?

Yanma kimyasında reaksiyonların karakteristik zaman ölçekleri çok farklı olabilir. Kaynak, kimyasal süreçlerin zaman ölçeklerinin yaklaşık 10⁻¹⁰ s ile 10⁰ s arasında değişebildiğini belirtmektedir. Çok hızlı ve çok yavaş kimyasal süreçlerin aynı denklem sisteminde bulunması sistemi stiff, yani sayısal olarak katı hale getirir.

Ayrıntılı mekanizma kullanıldığında akış denklemlerine ek olarak çok sayıda kimyasal türün taşınım denklemi ve yüzlerce kimyasal reaksiyonun kaynak terimi hesaplanabilir. Örneğin bu çalışmada kullanılan GRI3.0 mekanizması 53 tür ve 325 reaksiyon, SanDiego-2014 ise 48 tür ve 247 reaksiyon içermektedir.

Ayrıntılı çözümde hangi denklemler kullanılıyor?

Kaynak, ayrıntılı yanma modelinin Newtonyen akışkan varsayımı altında toplam kütle, türler, momentum ve enerji korunumuyla bir durum denkleminden oluştuğunu belirtir.

Denklem (1): süreklilik denklemi

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

Bu denklem toplam kütlenin korunmasını ifade eder. \(\rho\) yoğunluk, \(t\) zaman ve \(\mathbf{U}\) hız vektörüdür. İlk terim belirli bir noktadaki yoğunluğun zamanla değişimini, ikinci terim ise akışla taşınan kütlenin uzaysal dengesini temsil eder. Denklemde kimyasal reaksiyon nedeniyle toplam kütleyi yaratan veya yok eden bir kaynak terimi bulunmaz; kimyasal reaksiyonlar türleri birbirine dönüştürür fakat toplam kütle korunur.

Denklem (2): tür kütle kesri denklemi

\[ \frac{\partial(\rho Y_k)}{\partial t} +\nabla\cdot\left[\rho Y_k(\mathbf{U}+\mathbf{U}_c)\right] +\nabla\cdot \mathbf{j}_k = \dot{\omega}_k, \qquad k=1,\ldots,n_s-1 \]

\(Y_k\), \(k\) kimyasal türünün kütle kesridir. \(\mathbf{U}_c\) düzeltme hızını, \(\mathbf{j}_k\) difüzif kütle akısını, \(\dot{\omega}_k\) kimyasal reaksiyonların ilgili türü üretme veya tüketme hızını, \(n_s\) ise toplam tür sayısını ifade eder.

Toplam kütle kesirlerinin 1 olması gerektiğinden son tür ayrıca çözülmez ve:

\[ Y_s=1-\sum_{k=1}^{n_s-1}Y_k \]

ilişkisiyle hesaplanır.

Denklem (3): momentum denklemi

\[ \frac{\partial(\rho\mathbf{U})}{\partial t} = -\nabla\cdot(\rho\,\mathbf{U}\otimes\mathbf{U}) +\nabla\cdot \left[ \mu \left( \nabla\otimes\mathbf{U} + (\nabla\otimes\mathbf{U})^T -\frac{2}{3}(\nabla\cdot\mathbf{U})\mathbf{I} \right) \right] -\nabla p +\rho\mathbf{g} \]

Bu denklem akışkan momentumunun korunmasını ifade eder. \(\mathbf{U}\otimes\mathbf{U}\) hız vektörünün dyadik çarpımıdır; \(\mu\) dinamik viskozite, \(\mathbf{I}\) birim matris, \(p\) basınç ve \(\mathbf{g}\) yerçekimi ivmesidir.

İlk sağ taraf terimi momentumun taşınımını, viskoz terim akışkan içi momentum difüzyonunu, \(-\nabla p\) basınç kuvvetini ve \(\rho\mathbf{g}\) hacimsel yerçekimi kuvvetini temsil eder. Bununla birlikte çalışmanın 2D karşı akış doğrulama hesaplarında yerçekimi ihmal edilmiştir.

Denklem (4): enerji denklemi

\[ \frac{\partial(\rho h_s)}{\partial t} +\nabla\cdot(\rho\mathbf{U}h_s) = -\frac{\partial(\rho K)}{\partial t} -\nabla\cdot(\rho\mathbf{U}K) -\nabla\cdot\mathbf{q} +\dot{\omega}_T +\frac{\partial p}{\partial t} +\nabla\cdot\mathbf{q}_{rad} \]

\(h_s\) duyulur entalpi, \(K=\frac{1}{2}|\mathbf{U}|^2\) özgül kinetik enerji, \(\mathbf{q}\) ısı akısı, \(\dot{\omega}_T\) kimyasal reaksiyonlarla ilişkili ısı açığa çıkma hızı ve \(\mathbf{q}_{rad}\) radyatif ısı akısıdır.

Denklem, akışla taşınan duyulur entalpi ile kinetik enerji, iletim/difüzyon yoluyla ısı taşınımı, kimyasal ısı açığa çıkışı, basınç değişimi ve radyasyon arasındaki enerji dengesini temsil eder.

Denklem (5): ideal gaz durum denklemi

\[ \rho=\frac{p\overline{M}}{RT} \]

\(\overline{M}\) karışımın ortalama molar kütlesi, \(R\) evrensel gaz sabiti ve \(T\) sıcaklıktır. Aynı basınç ve bileşim altında sıcaklık arttıkça yoğunluk matematiksel olarak azalır; aynı sıcaklık ve bileşim altında basınç arttıkça yoğunluk artar.

Denklem (6): tür difüzyon akısı

Hirschfelder–Curtiss yaklaşımı altında kaynak tür difüzyon akısını:

\[ \mathbf{j}_k=-\rho D_k\nabla Y_k \]

şeklinde verir. \(D_k\), \(k\) türünün karışımın geri kalanına göre difüzyon katsayısıdır. Negatif işaret, difüzyonun kütle kesrinin yüksek olduğu bölgeden düşük olduğu bölgeye yöneldiğini gösterir.

Denklem (7): Lewis sayısı 1 varsayımı

Lewis sayısı bir kabul edildiğinde kaynak:

\[ D_k=\frac{\lambda}{\rho C_p} \]

ilişkisini kullanır. \(\lambda\) ısıl iletkenlik, \(C_p\) sabit basınçtaki özgül ısıdır. Lewis sayısının 1 kabulü, ısıl difüzyon ile kütle difüzyonu arasında belirli bir eşdeğerlik varsayımı getirir.

Bu, REDIM yönteminin her durumda yalnız Le=1 ile çalışabileceği anlamına gelmez. Yazarlar REDIM'in ayrıntılı moleküler taşınımı da ele alabildiğini önceki çalışmalara dayanarak belirtmektedir. Bu makalenin belirli karşı akış testlerinde ise karşılaştırmayı kimyasal indirgemeye odaklamak için Le=1 kullanılmıştır.

Denklem (8): optik olarak ince radyasyon modeli

\[ \nabla\cdot\mathbf{q}_{rad} = 4\sigma \left(T^4-T_b^4\right) \sum_i\left(p_i a_{p_i}\right) \]

\(\sigma\) Stefan–Boltzmann sabiti, \(T_b\) çevre sıcaklığı, \(p_i\) radyasyona katılan türün kısmi basıncı ve \(a_{p_i}\) Planck ortalama absorbsiyon katsayısıdır. Kaynakta H₂O, CO ve CO₂ radyatif türler olarak ele alınmaktadır.

Radyasyon denklemi genel EBIdnsFoam denklem takımında bulunmasına rağmen çalışmanın 2D karşı akış doğrulamasında radyasyon ve yerçekimi ihmal edilmiştir.

EBIdnsFoam ne yapıyor?

EBIdnsFoam, OpenFOAM ile Cantera'yı birleştiren ayrıntılı kimya çözücüsüdür. OpenFOAM toplam kütle, momentum, enerji ve tür korunum denklemlerini çözer; basınç, sıcaklık ve tür kütle kesirlerini Cantera'ya aktarır. Cantera ise moleküler viskozite, ısıl iletkenlik, tür difüzyon katsayıları ve kimyasal reaksiyon hızları gibi termokimyasal/taşınım özelliklerini hesaplar.

REDIM'in temel durum vektörü

Kaynak ayrıntılı reaksiyon-akış sistemini:

\[ \Psi= (h,p,\phi_1,\ldots,\phi_{n_s})^T \]

durum vektörüyle ifade eder. Burada:

\[ \phi_k=\frac{Y_k}{M_k} \]

özgül mol sayısıdır; \(M_k\) ilgili türün molar kütlesidir.

Denklem (9): ayrıntılı durum uzayının evrimi

\[ \frac{\partial\Psi}{\partial t} = F(\Psi) -\mathbf{U}\cdot\nabla\Psi -\frac{1}{\rho}\nabla\cdot(D\cdot\nabla\Psi) \equiv \Phi(\Psi) \]

Bu denklem REDIM'in matematiksel başlangıç noktasıdır. \(F(\Psi)\) kimyasal kaynak terimlerini, \(\mathbf{U}\cdot\nabla\Psi\) konvektif taşınımı ve \(D\) matrisi moleküler difüzyonu temsil eder. \(\Phi(\Psi)\), durum vektörünün bütün değişim hızlarını tek bir vektör fonksiyonunda toplar.

Yüksek boyutlu sistem düşük boyutlu manifolda nasıl indirgeniyor?

REDIM, sistem çözümünün yüksek boyutlu durum uzayında rastgele her noktayı kullanmadığını, fiziksel/kimyasal olarak erişilen durumların daha düşük boyutlu bir yapı çevresinde yoğunlaştığını kullanır.

Denklem (10) bu manifoldu:

\[ M= \left\{ \Psi:\Psi=\Psi(\theta), \quad \Psi:\mathbb{R}^{m_s}\rightarrow\mathbb{R}^{n} \right\} \]

biçiminde tanımlar. \(n\) ayrıntılı durum uzayının boyutu, \(m_s\) indirgenmiş manifoldun boyutudur ve amaç:

\[ m_s\ll n \]

olacak biçimde çok sayıda termokimyasal değişkeni birkaç koordinatla temsil etmektir.

Denklem (11): manifold değişmezlik koşulu

\[ \left( I-\Psi_\theta\Psi_\theta^+ \right)\cdot\Phi(\Psi)=0 \]

\(\Psi_\theta\), \(\Psi\)'nin indirgenmiş koordinat \(\theta\)'ya göre türevlerinden oluşan Jacobian benzeri matristir. \(\Psi_\theta^+\) onun psödo-tersidir.

\(\Psi_\theta\Psi_\theta^+\) manifoldun teğet uzayına projeksiyon yapar; dolayısıyla:

\[ I-\Psi_\theta\Psi_\theta^+ \]

manifoldun normal yönünü seçer. Denklem, sistemin değişim vektörünün manifold dışına doğru bileşeninin sıfır olması gerektiğini söyler. Başka bir ifadeyle manifold gerçekten değişmez hale geldiğinde sistem dinamiği manifold üzerinde kalır.

Psödo-ters neden gerekli?

İndirgenmiş koordinat sayısı ile ayrıntılı durum değişkeni sayısı aynı olmadığı için \(\Psi_\theta\) genellikle kare bir matris değildir. Bu nedenle sıradan matris tersi yerine psödo-ters kullanılır.

Denklem (12):

\[ \Psi_\theta^+\Psi_\theta = I_{m_s\times m_s} \]

ilişkisini verir.

Denklem (13): REDIM manifoldunun oluşturulması

Eşit difüzivite varsayımı \(D=dI\) altında kaynak manifold evrimini:

\[ \frac{\partial\Psi(\theta)}{\partial t} = \left( I-\Psi_\theta(\theta)\Psi_\theta^+(\theta) \right) \cdot \left( F(\Psi(\theta)) - \frac{1}{\rho} \left[ \left( d\Psi_\theta\cdot\chi(\theta) \right)_\theta \cdot\chi(\theta) \right] \right) \]

şeklinde verir. Burada:

\[ \chi(\theta)=\operatorname{grad}(\theta) \]

manifold için gradyan kestirimidir.

Bu denklem gerçek CFD alanındaki zaman evrimiyle karıştırılmamalıdır. Manifold üretiminde başlangıç koşulu, sınır koşulları ve gradyan kestirimleri tanımlanır; Denklem (13) matematiksel olarak durağan manifold elde edilinceye kadar ilerletilir.

REDIM tablosu nasıl oluşturuluyor?

Kaynağın verdiği süreç şöyledir:

1. Başlangıç durumu seçilir.
2. Uygun \(\theta\) parametrizasyon vektörü belirlenir.
3. Denklem (13) için sınır koşulları tanımlanır.
4. \(\chi(\theta)\) gradyan kestirimi sağlanır.
5. Denklem (13), \(t\rightarrow\infty\) yönünde durağan \(\Psi(\theta)\) çözümüne kadar entegre edilir.
6. \(\theta\), \(\Psi(\theta)\), \(\Psi_\theta\), \(\Psi_\theta^+\), \(D\Psi_\theta\) ve \(\Psi_\theta^+F\) tablodaki ağ noktalarında saklanır.

İlk tahmin neden hesap süresini etkiliyor?

2D REDIM tablosunun oluşturulmasında kaynak iki başlangıç yaklaşımını karşılaştırmıştır. Bir manifold 80, 100, 120, 140, 160, 180, 200 ve 300 s⁻¹ gerinim hızlarındaki flamelet'ler ve karışım doğrusu kullanılarak; diğer manifold yalnız 80 s⁻¹ flamelet'i ile karışım doğrusundan oluşturulmuştur.

Nihai manifoldlar arasında yalnız küçük farklar görülmesine rağmen yalnız iki başlangıç flamelet'i kullanılan kötü başlangıç tahmininin yakınsama süresi daha iyi başlangıç tahmininin yaklaşık üç katı olmuştur. Yani başlangıç tahmini teorik nihai çözümden çok çözümün ne kadar hızlı yakınsadığını etkilemiştir.

Genelleştirilmiş ve Fiziksel Koordinatlar Arasındaki Fark Nedir?

Genelleştirilmiş koordinatlar REDIM manifoldunun düzenli, yerel olarak ortogonal ve eş aralıklı tablo indekslerini doğrudan kullanarak interpolasyonu kolaylaştırırken fiziksel koordinatlar gerçek tür bileşimlerinin doğrusal kombinasyonlarını kullanarak sonuçları daha anlaşılır hale getirir; ancak tek-değerli ve monoton uygun fiziksel parametrizasyonu bulmak daha zordur.

Denklem (14): genelleştirilmiş koordinatın evrimi

Manifold oluşturulduktan sonra ayrıntılı durum vektörü \(\Psi\), düşük boyutlu \(\theta\) koordinatının fonksiyonudur. Kaynak zincir kuralı ve psödo-ters kullanarak:

\[ \frac{\partial\theta}{\partial t} = \Psi_\theta^+\cdot F(\Psi(\theta)) - \mathbf{U}\cdot\nabla\theta - \frac{1}{\rho} \Psi_\theta^+ \cdot \nabla\cdot \left( D\Psi_\theta\cdot\nabla\theta \right) \]

denklemini verir.

Buradaki önemli değişim, artık bütün türlerin doğrudan çözülmesi yerine \(m_s\) adet indirgenmiş koordinatın çözülmesidir.

REDIMFoam kaç denklem çözüyor?

Kaynak REDIM yaklaşımında toplam:

\[ 3+m_s \]

denklem çözüldüğünü belirtir: süreklilik, momentum, durum denklemi ve \(m_s\) adet REDIM koordinat denklemi. Örneğin 2D manifold için \(m_s=2\) olduğundan iki adet \(\theta\) denklemi vardır.

Bu sayı ayrıntılı kimyanın bütün tür denklemlerinden çok daha küçüktür; hesaplama maliyetindeki azalmanın temel nedenlerinden biri budur.

Denklem (15): 1D karşı akıştaki gerinim oranı

\[ \alpha= \sqrt{\frac{J}{\rho_{ub}}} \]

Burada \(J\) teğetsel basınç gradyanını, \(\rho_{ub}\) yanmamış yakıt karışımının yoğunluğunu temsil eder. \(\alpha\) birimi s⁻¹'dir ve karşı karşıya gelen akışların alevi ne ölçüde gerdiğini karakterize eder.

Denklem (16): genelleştirilmiş koordinatlarla REDIMFoam

\[ \frac{\partial(\rho\theta)}{\partial t} + \nabla\cdot(\rho\mathbf{U}\cdot\theta) + \Psi_\theta^+ \cdot \nabla\cdot \left( D\Psi_\theta\cdot\nabla\theta \right) = \rho\Psi_\theta^+\cdot F(\Psi(\theta)) \]

Bu, REDIMFoam içinde çözülen azaltılmış taşınım denklemidir. Denklemde konveksiyon ve difüzyon korunurken ayrıntılı kimyasal kaynak teriminin etkisi düşük boyutlu koordinat uzayına projekte edilir.

Fiziksel koordinatlar nasıl tanımlandı?

Kaynak fiziksel koordinatı sabit bir \(C\) parametrizasyon matrisiyle:

\[ \xi=C\Psi \]

şeklinde tanımlar. Örnek iki boyutlu parametrizasyon:

\[ C= \begin{bmatrix} 0 & 0 & 1 & 0 & 0 & \cdots & 0\\ 0 & 0 & 0 & 1 & 0.5 & \cdots & 0 \end{bmatrix} \]

ve:

\[ \xi= \begin{pmatrix} \phi_{N_2}\\ \phi_{CO_2}+0.5\phi_{CO} \end{pmatrix} \]

şeklindedir.

Böylece:

\[ \xi_1=\phi_{N_2} \]

ve:

\[ \xi_2=\phi_{CO_2}+0.5\phi_{CO} \]

olur. Bu koordinatların doğrudan kimyasal bileşimle ilişkili olması post-processing ve fiziksel yorumlamayı kolaylaştırır.

Denklem (18): fiziksel koordinat evrimi

\[ \frac{\partial\xi}{\partial t} = \Psi_\xi^+\cdot F(\Psi(\xi)) - \mathbf{U}\cdot\nabla\xi - \frac{1}{\rho} \Psi_\xi^+ \cdot \nabla\cdot \left( D\Psi_\xi\cdot\nabla\xi \right) \]

Denklem (19): fiziksel koordinatların REDIMFoam denklem biçimi

\[ \frac{\partial(\rho\xi)}{\partial t} + \nabla\cdot(\rho\mathbf{U}\cdot\xi) + \Psi_\xi^+ \cdot \nabla\cdot \left( D\Psi_\xi\cdot\nabla\xi \right) = \rho\Psi_\xi^+\cdot F(\Psi(\xi)) \]

Denklem (20): fiziksel koordinat projeksiyon matrisi

\[ \Psi_\xi^+ = C\Psi_\theta\Psi_\theta^+ \]

Denklem (21): fiziksel koordinat Jacobian'ı

\[ \Psi_\xi = \Psi_\theta\theta_\xi = \Psi_\theta \left(C\Psi_\theta\right)^{-1} \]

Bu ilişkiler, genelleştirilmiş koordinatlar üzerinde önceden oluşturulan REDIM tablosunun fiziksel koordinat biçimine dönüştürülebilmesini sağlar.

Hangisi daha iyi çıktı?

1D serbest yayılan metan/hava alevinde iki koordinat sistemi sıcaklık ve tür profilleri bakımından birbirine ve ayrıntılı çözümlere çok yakın sonuçlar vermiştir. CPU süreleri de 0,51 s ve 0,55 s ile birbirine yakındır.

Yazarlar genelleştirilmiş koordinatların fiziksel anlamının daha az doğrudan olmasına rağmen tablo ağının yerel olarak ortogonal ve eş aralıklı olması nedeniyle interpolasyonun daha verimli ve sağlam olduğunu belirtmekte ve çalışmanın sonraki simülasyonlarında genelleştirilmiş koordinatları kullanmaktadır.

REDIMFoam Hesaplama Maliyetini Ne Kadar Azalttı?

1D adyabatik metan/hava serbest alevi doğrulamasında REDIMFoam, aynı zaman adımı ve ağ boyutu için ayrıntılı EBIdnsFoam hesap süresinin yaklaşık %11–12'sini kullanmış; böylece bu testte CPU süresini yaklaşık sekiz ila dokuz kat azaltmıştır.

Denklem (22): göreli CPU süresi

Kaynak göreli CPU süresini:

\[ t_{CPU,rel} = \frac{ t_{CPU,\mathrm{REDIM}} }{ t_{CPU,\mathrm{detailed}} } \]

şeklinde tanımlar.

Karşılaştırmalar aynı zaman adımı, aynı mesh boyutu ve tek CPU çekirdeğinde yapılmıştır. Bu koşul önemlidir; farklı mesh, farklı işlemci veya farklı paralelleştirme düzenleri arasındaki ham süreler doğrudan bu oranla karşılaştırılamaz.

REDIM'in maliyet avantajının iki temel nedeni vardır. Birincisi, \(n_s\) tür denklemi yerine yalnız \(m_s\) indirgenmiş koordinat denklemi çözülmesidir. İkincisi, yüksek derecede doğrusal olmayan kimyasal kaynak terimlerinin her CFD hücresinde yeniden çözülmesi yerine termokimyasal büyüklüklerin REDIM lookup tablosundan lineer interpolasyonla geri alınabilmesidir.

CPU avantajı hangi sınır içinde yorumlanmalıdır?

4,49 s, 0,51 s ve 0,55 s değerleri belirli 1D serbest alev doğrulamasına aittir. Bunları her 2D, 3D, türbülanslı veya endüstriyel yanma probleminde aynı hızlanmanın garanti edildiği şeklinde genellemek doğru değildir.

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

1D serbest alev doğrulaması

Yeni REDIMFoam çözücüsünün ilk doğrulaması adyabatik, bir boyutlu, ön karışımlı metan/hava serbest aleviyle yapılmıştır. Yakıt karışımı stokiyometriktir:

\[ \Phi=1 \]

ve unity Lewis number varsayımı uygulanmıştır.

Bu test için \(m_s=1\) olan tek boyutlu REDIM hazırlanmış ve progress variable olarak \(\phi_{CO_2}\) seçilmiştir.

SanDiego-2014 mekanizmasıyla EBIdnsFoam ve INSFLA tarafından hesaplanan ana türler ile küçük türlerin profilleri birbirleriyle uyum göstermiş; hem genelleştirilmiş hem fiziksel koordinatlı REDIM sonuçları da bu ayrıntılı çözümleri yakından izlemiştir.

2D karşı akışlı alev geometrisi

İkinci ana doğrulama iki boyutlu, eksensimetrik, laminer karşı akışlı difüzyon alevinde yapılmıştır. Yakıt tarafında hacimce %25 CH₄ + %75 hava, diğer tarafta %21 O₂ + %79 N₂ içeren hava akışı bulunmaktadır.

Yakıt karışımı için kaynak:

\[ \Phi=3.17 \]

değerini vermektedir.

Her iki karşılıklı nozülün çapı 2 cm, nozüller arası mesafe 2 cm ve hesaplama bölgesinin merkez çizgiden yukarı doğru yüksekliği 4 cm'dir.

2D hesaplamalarda hangi varsayımlar kullanıldı?

• Akış iki boyutlu ve eksensimetriktir.
• Alev laminer karşı akışlı difüzyon alevidir.
• Her iki nozülün giriş hızlarının büyüklüğü eşit ve yönleri zıttır.
• Unity Lewis number varsayılmıştır.
• Yerçekimi ihmal edilmiştir.
• Radyasyon ihmal edilmiştir.
• 2D REDIM için \(m_s=2\) kullanılmıştır.
• Dış coflow hava hızı 0,01 m/s'dir.

Mesh bağımlılığı nasıl test edildi?

Mesh 2'den Mesh 3'e geçildiğinde merkez bölgesindeki x doğrultusu daha da inceltilirken y doğrultusundaki temel çözünürlük aynı kalmıştır.

Ayrıntılı kimyada mesh sonucu

REDIM'de mesh sonucu

Kaba Mesh 1'de küçük türlerin OH ve O profillerinde daha belirgin farklar görülmüş; Mesh 2 ve Mesh 3 sonuçları birbirine yakın çıkmıştır. REDIM'de Mesh 1'in alev kalınlığı da ince ağlardan biraz daha büyük olmasına rağmen kaynak bu farkı %2'nin altında raporlamıştır.

Bu sonuçlara dayanarak araştırmacılar sonraki hesaplamalarda Mesh 2'yi yeterli kabul etmiştir.

Progress Variable Seçimi REDIM Sonucunu Ne Kadar Değiştiriyor?

Çalışmanın test ettiği iki geçerli parametrizasyon sıcaklık, H₂O ve OH gibi büyüklüklerde neredeyse aynı sonuçları vermiş; bu nedenle bir-bir eşlemeyi bozmayan uygun koordinat seçimleri altında REDIM çözümünün progress-variable seçimine büyük ölçüde değişmez kaldığı gösterilmiştir.

Karşılaştırılan iki parametrizasyon

Birinci seçenek:

\[ \phi= (\phi_{N_2},\phi_{CO_2})^T \]

İkinci seçenek:

\[ \phi= (\phi_{N_2},\phi_{CO_2}+0.5\phi_{CO})^T \]

olarak seçilmiştir.

O kütle kesrinde ayrıntılı çözümden küçük bir fark gözlenirken iki REDIM parametrizasyonunun kendi aralarındaki profilleri neredeyse aynıdır.

Bu sonuç “hangi progress variable seçilirse seçilsin sonuç aynıdır” anlamına gelmez. Kaynağın şartı, parametrizasyonun manifold üzerinde bir-bir eşlemeyi sağlaması ve katlanmış bir bölge oluşturmamasıdır.

Reaksiyon mekanizması değiştirildiğinde ne oldu?

GRI3.0 ve SanDiego-2014 mekanizmaları için ayrı REDIM tabloları oluşturulmuştur.

Kaynak REDIM'in getirdiği tepe sıcaklığı hatasının %0,1'in altında olduğunu belirtmektedir. CO₂ gibi ana türlerde de fark ihmal edilebilir düzeydedir.

Buna karşılık O ve OH gibi küçük türlerde GRI3.0 ile SanDiego-2014 mekanizmalarının kendi aralarındaki fark, aynı mekanizmanın ayrıntılı ve REDIM çözümü arasındaki farktan daha büyüktür. Bu bulgu REDIM model indirgeme hatasının, bu testlerde ayrıntılı mekanizma seçimine bağlı model farkından daha küçük olduğunu göstermektedir.

Giriş hızı ve global gerinim oranı

2D karşı akışta global strain rate \(a_g\) kaynakta:

\[ a_g= \frac{2U_rU_F}{d} \left( 1+ \frac{\sqrt{\rho_F}} {U_r\sqrt{\rho_O}} \right) \]

şeklinde tanımlanır.

Burada:

• \(U_r=U_O/U_F\): oksitleyici ve yakıt giriş hızlarının oranı,
• \(U_O\): hava/oksitleyici giriş hızı,
• \(U_F\): yakıt tarafı giriş hızı,
• \(d\): iki nozül arasındaki mesafe,
• \(\rho_F\): yakıt karışımının yoğunluğu,
• \(\rho_O\): havanın yoğunluğudur.

Bu çalışmada iki giriş hızının büyüklüğü eşit olduğu için:

\[ U_r=1 \]

olur.

Gerinim oranının birimi s⁻¹'dir. Gerinim oranı büyüdükçe karşı akışın alev tabakasını daha güçlü zorladığı düşünülebilir; yeterince yüksek değerlerde kimyasal zaman ölçekleri ile taşınım arasındaki denge bozulabilir ve alev sönmeye yaklaşabilir.

Üç giriş hızındaki sonuçlar

Giriş hızı ve global gerinim oranı yükseldikçe her iki çözümde de tepe sıcaklığı azalmaktadır. REDIM aynı eğilimi ayrıntılı kimyayla birlikte takip etmektedir.

Kaynak global gerinim oranını 97–972 s⁻¹ aralığında genişlettiğinde sıcaklık ve H₂O sonuçlarının bütün aralıkta ayrıntılı çözümlerle iyi uyum gösterdiğini belirtmektedir. Yaklaşık 600 s⁻¹'in üzerinde OH gibi küçük türlerde farklar büyümeye başlamaktadır.

Yazarlar küçük türler için görülen bu farkların daha yüksek boyutlu manifoldlar kullanılarak azaltılabileceğini belirtmektedir. Bu, çalışmanın önemli bir model sınırıdır: düşük boyutlu REDIM ana termokimyasal yapıyı iyi yakalarken bütün küçük radikal türlerini her koşulda aynı kesinlikle temsil etmeyebilir.

REDIM Sönme Rejimini Manifold Oluşturulurken Sönme Flamelet'i Verilmeden Tahmin Edebiliyor mu?

Bu çalışmanın zamanla değişen 2D karşı akış testinde REDIM, manifold oluşturulurken açıkça sönme rejimine ait flamelet'ler verilmemesine rağmen sıcaklık, H₂O, O₂ ve genel sönme yörüngesini ayrıntılı çözüme yakın biçimde izlemiş; en belirgin fark OH gibi küçük bir radikal türde görülmüştür.

Sönme testi nasıl başlatıldı?

Kararlı durum deneylerinin ardından global gerinim oranı zamanla değişen basamak fonksiyonuyla 97 s⁻¹'den sönme sınırının üzerine yükseltilmiştir.

Kaynak aynı alt bölümde iki farklı quenching-limit değeri kullanmaktadır:

• ilk paragrafta 1457 s⁻¹,
• hemen ardından 1550 s⁻¹.

Makale bu iki değerin farklı tanımlar, farklı çözüm yöntemleri veya farklı aşamalar olup olmadığını açıklamamaktadır. Bu nedenle Verianla bunlardan birini sessizce düzeltmemektedir.

Kaynağa göre global gerinim oranının sönme sınırını iki milisaniyeden daha kısa süre aşması alevin sönmesine yol açmıştır.

Sönme yolu iki aşamaya ayrılıyor

Yazarlar sönme sürecini iki bölüm halinde açıklar:

1. artan gerinim oranının kontrol ettiği yavaş değişim dönemi,
2. quenching limitine ulaşıldıktan sonraki hızlı sönme dönemi.

REDIM çözümü tanımı gereği REDIM yavaş manifoldunun üzerinde hareket eder. Buna rağmen sıcaklık, H₂O ve O₂'nin durum uzayındaki gelişimi ayrıntılı EBIdnsFoam çözümüyle neredeyse üst üste gelmiştir.

OH kütle kesrinde daha görünür bir fark vardır. Kaynak buna rağmen REDIM'in 2D karşı akış alevinin sönme evrimini kabul edilebilir doğrulukla gösterdiği sonucuna varmaktadır.

Bu sonuç neden önemli?

Klasik sabit flamelet tabanlı yöntemlerde sönme rejimini doğru temsil edebilmek için sönmeye ait ek flamelet bilgileri gerekebilir. Kaynak, REDIM'in kimyasal reaksiyon ile moleküler taşınım arasındaki etkileşimi manifold üretiminde hesaba katması nedeniyle kararlı ve geçici rejimleri aynı düşük boyutlu yapı içinde daha esnek biçimde temsil edebildiğini savunmaktadır.

Çalışmanın doğrudan desteklediği sonuçlar

• REDIM yöntemi OpenFOAM tabanlı REDIMFoam çözücüsüne uygulanmıştır.
• Genelleştirilmiş ve fiziksel koordinat uygulamaları 1D serbest alevde birbirine ve ayrıntılı kimya sonuçlarına yakın sonuç vermiştir.
• Genelleştirilmiş koordinatlar interpolasyon verimliliği ve sağlamlığı nedeniyle çalışmanın sonraki hesaplarında tercih edilmiştir.
• 1D doğrulamada REDIMFoam CPU süresi ayrıntılı mekanizmanın yaklaşık %11–12'sine düşmüştür.
• Mesh 2, kaynak tarafından 2D sıcaklık ve küçük tür yapılarını çözmek için yeterli kabul edilmiştir.
• İki farklı geçerli progress-variable parametrizasyonu çok benzer REDIM sonuçları vermiştir.
• GRI3.0 ve SanDiego-2014 için REDIM sonuçlarının her biri kendi ayrıntılı mekanizmasıyla yakın uyum göstermiştir.
• Giriş hızı arttıkça tepe sıcaklığındaki düşüş REDIM tarafından da yakalanmıştır.
• REDIM, manifold üretiminde açık sönme flamelet'leri verilmemesine rağmen zamanla gelişen sönme sürecini takip edebilmiştir.

Çalışmanın doğrudan desteklemediği sonuçlar

• Çalışma fiziksel bir brülör deneyi değildir.
• REDIMFoam'ın gerçek motor veya gaz türbininde aynı hata ve hızlanmayla çalıştığı gösterilmemiştir.
• Bu makalede türbülanslı alev için yeni doğrudan simülasyon doğrulaması yapılmamıştır; sonuç bölümündeki türbülans bağlantısı önceki literatüre dayanmaktadır.
• Unity Lewis number kullanılan test sonuçları bütün ayrıntılı diferansiyel difüzyon durumlarına doğrudan genellenmemelidir.
• 2D karşı akış hesaplarında radyasyon ve yerçekimi ihmal edildiği için sonuçlar bu etkilerin baskın olduğu yanma problemlerini doğrudan doğrulamaz.
• Küçük radikal türlerde yüksek gerinim oranlarında görülen farkların tamamen ortadan kalktığı gösterilmemiştir.
• Bir boyutlu CPU karşılaştırmasındaki %11–12 göreli süre bütün CFD problemleri için garanti edilmiş bir hızlanma oranı değildir.

Çalışmanın temel yöntemsel sınırları

REDIM'in doğruluğu yalnız “REDIM yöntemi iyi veya kötü” şeklinde tek bir sabitle belirlenmez. Manifold boyutu, başlangıç ve sınır koşulları, gradyan kestirimi, seçilen ayrıntılı reaksiyon mekanizması, koordinat parametrizasyonu, mesh çözünürlüğü ve incelenen termokimyasal rejim sonuçları etkileyebilir.

Çalışma özellikle küçük türlerin yüksek gerinim oranlarındaki farklarının daha yüksek boyutlu manifoldlarla iyileştirilebileceğini belirtmektedir. Daha yüksek manifold boyutu ise azaltmanın hesaplama avantajıyla doğruluk arasında yeni bir denge oluşturabilir.

Genelleştirilmiş koordinatların interpolasyon açısından avantajlı olması, fiziksel koordinatlardan bilimsel olarak “daha doğru” oldukları anlamına gelmemektedir. 1D doğrulamada iki temsil neredeyse aynı doğruluk ve CPU maliyetini göstermiştir; öneri daha çok interpolasyon uygulamasının pratik sağlamlığına dayanmaktadır.

Kaynak ve Yöntem Notu

Özgün çalışma: Implementation and Validation of a REDIM-Based CFD Solver for Combustion Applications

Yazarlar: Ningyi Li; Thorsten Zirwes; Ulrich Maas.

Sorumlu yazar: Ningyi Li.

E-posta: ningyi.li@partner.kit.edu

Kurum 1: Institute of Technical Thermodynamics, Karlsruhe Institute of Technology, Engelbert-Arndold-Strasse 4, Karlsruhe, 76131, Germany.

Kurum 2: Institute for Reactive Flows, University of Stuttgart, Pfaffenwaldring 31, Stuttgart, 70569, Germany.

DOI: 10.18321/ectj1681.

Dergi: Eurasian Chemico-Technological Journal.

Sayfalar: 25–43.

Resmî çevrimiçi yayın tarihi: 10 Nisan 2026.

Makale işlem tarihleri: Alındı: 16 Ekim 2025; revize sürüm: 26 Aralık 2025; kabul: 18 Şubat 2026.

Kaynak türü: Hesaplamalı akışkanlar dinamiği, ayrıntılı kimyasal kinetik ve indirgenmiş kimya modeli doğrulamasını birleştiren özgün hesaplamalı araştırma makalesi.

Hakemlik: Derginin resmî politikası hakemli yayın olduğunu belirtmektedir. Araştırma makalelerinin kabulü için en az iki olumlu hakem değerlendirmesi gerektiği derginin yazar rehberinde açıklanmaktadır.

Lisans: Creative Commons Attribution 4.0 International — CC BY 4.0.

Yayınevi bilgisi: PDF'nin telif satırı “Published by al-Farabi Kazakh National University” ifadesini kullanmaktadır. Derginin resmî About sayfası yayıncıyı Institute of Combustion Problems, Almaty olarak tanımlamaktadır. Bu iki kurumsal ifade birbirinin yerine geçirilmemiştir.

Cilt bilgisi hakkında kaynak farkı: Yüklenen PDF üstbilgisinde “Eurasian Chem.-Technol. J. 28 (2026) 25–43” yazmaktadır. Derginin resmî çevrimiçi makale ve sayı kaydı ise “Vol. 26 No. 1 (2026), 25–43” biçimindedir. KIT ve University of Stuttgart kurumsal bibliyografik kayıtları da cilt 26'yı kullanmaktadır. Verianla bibliyografik notu bu farkı açıkça korumaktadır.

Temel yazılım: REDIMFoam, OpenFOAM tabanlıdır ve kaynak OpenFOAM v2006 temelini belirtmektedir. Ayrıntılı referans çözümler EBIdnsFoam ile oluşturulmuş; EBIdnsFoam OpenFOAM ile Cantera'yı birleştirmiştir. REDIM manifold üretiminde INSFLA ve HOMREA kodlarının değiştirilmiş/genişletilmiş uygulamalarından yararlanılmaktadır.

Ayrıntılı kimyasal mekanizmalar: GRI3.0 — 53 tür ve 325 reaksiyon; SanDiego-2014 — 48 tür ve 247 reaksiyon.

1D doğrulama: Adyabatik, ön karışımlı metan/hava serbest alev; Φ=1; Le=1; tek boyutlu REDIM, \(m_s=1\), progress variable \(\phi_{CO_2}\).

1D CPU sonuçları: EBIdnsFoam 4,49 s; REDIMFoam genelleştirilmiş koordinat 0,51 s; fiziksel koordinat 0,55 s.

2D doğrulama: Eksensimetrik laminer karşı akışlı difüzyon alevi; yakıt akışı hacimce %25 CH₄ + %75 hava, karşı akış hava %21 O₂ + %79 N₂; nozül çapı 2 cm; nozül mesafesi 2 cm; hesaplama yüksekliği 4 cm; coflow 0,01 m/s.

2D varsayımlar: Unity Lewis number; yerçekimi ihmal edilmiş; radyasyon ihmal edilmiş; \(m_s=2\) REDIM kullanılmıştır.

Seçilen mesh: Mesh 2; merkez bölgede Δx=0,05 mm ve Δy=0,125 mm.

Progress-variable testi: \((\phi_{N_2},\phi_{CO_2})\) ve \((\phi_{N_2},\phi_{CO_2}+0.5\phi_{CO})\).

Progress-variable tepe sıcaklıkları: Ayrıntılı GRI3.0 2016 K; ilk REDIM parametrizasyonu 2015 K; ikinci REDIM parametrizasyonu 2016 K.

Mekanizma karşılaştırması: GRI3.0 ayrıntılı 2016 K / REDIM 2015 K; SanDiego-2014 ayrıntılı 2031 K / REDIM 2032 K.

Giriş hızı sonuçları: 0,5 m/s → 97 s⁻¹ → 2043/2038 K; 1,0 m/s → 194 s⁻¹ → 2016/2015 K; 1,5 m/s → 291 s⁻¹ → 1997/1987 K; ilk sıcaklık ayrıntılı, ikinci REDIM sonucudur.

Kaynak içi sayı farkı: Sönme bölümünde quenching limit için hem 1457 s⁻¹ hem 1550 s⁻¹ ifadeleri bulunmaktadır. Kaynak bu iki değerin ilişkisini ayrıca açıklamamaktadır.

Finansman: İncelenen ana PDF içinde ayrı bir finansman beyanı bulunmamaktadır.

Veri erişilebilirliği: İncelenen ana PDF içinde ayrı bir Data Availability beyanı bulunmamaktadır.

Çıkar çatışması: İncelenen ana PDF içinde ayrı bir competing interests / conflict of interest beyanı bulunmamaktadır.

Yazar katkıları: İncelenen ana PDF içinde ayrı bir CRediT/yazar katkıları bölümü bulunmamaktadır. Eşit katkı/eş birinci yazarlık beyanı görülmemiştir.

Teşekkür: Yazarlar REDIM uygulamasına ilişkin yararlı tartışmalar için Dr. Viatcheslav Bykov ve Dr. Chunkan Yu'ya teşekkür etmektedir.

Temel yorum sınırı: Çalışmadaki doğrulama, ayrıntılı kimyasal kinetik kullanan başka sayısal çözümlerle karşılaştırmadır. Bu çalışma kendi kapsamında fiziksel alev deneyleriyle yeni bir deneysel doğrulama yapmamaktadır.

Türbülans sınırı: Sonuç bölümünde benzer çıkarımların önceki çalışmalardaki türbülanslı alevler için de geçerli olduğuna atıf yapılmaktadır; ancak bu makalenin yeni doğrulama senaryoları 1D ve 2D laminer alevlerdir.