Bu çalışma, karasal yerçekimi altında %6 hacimsel H₂ içeren ultra-zayıf hidrojen–hava karışımında oluşan alev çekirdeğinin büyümesi, yükselmesi, şekil değiştirmesi ve parçalanmasını sayısal olarak incelemektedir. Ultra-zayıf hidrojen alev çekirdeği, ateşleme sonrasında hidrojenin tercihli ve termal difüzyonu, kimyasal ısı açığa çıkışı, radyatif soğuma ve kaldırma kuvvetinin birlikte belirlediği küçük ve zamana bağlı bir reaksiyon yapısıdır. Araştırmada ayrıntılı hidrojen kimyası ve düşük-Mach-sayısı gaz dinamiği kullanılmış; Soret difüzyonu ile radyatif ısı kaybı bağımsız olarak açılıp kapatılarak dört fiziksel model kombinasyonu karşılaştırılmıştır.

Erken ateşleme-sonrası aşama bir boyutlu küresel simülasyonlarla incelenmiş ve yaklaşık ilk 200 ms boyunca alev cephesi yarıçapının karesinin zamanla yaklaşık doğrusal değiştiği, yani \(R_f^2 \simeq Kt\) biçiminde difüzyon kontrollü bir büyüme rejimi izlediği bulunmuştur. Soret difüzyonu etkin büyüme katsayısı \(K\)'yı artırırken radyatif ısı kaybı azaltmıştır. Bu davranış sabit hızlı bir küresel alev yayılımından farklıdır; çekirdek büyüdükçe ön cephenin anlık radyal hızı azalmaktadır.

İki boyutlu simülasyonlarda çekirdek büyüdükçe kaldırma kuvveti baskın hale gelmiş, sıcak ürün bölgesi yukarı yükselmiş ve kendi toroidal girdabını oluşturmuştur. Başlangıçta yaklaşık küresel olan yapı daha sonra basık eliptik bir kapağa dönüşmüştür. Soret difüzyonu daha büyük ve daha hızlı yükselen çekirdekler oluşturmuş ve üst kapağı daha dairesel tutmuştur. Radyasyon yükselme yörüngesini daha zayıf ölçüde etkilemiş, ancak göreli yanal basıklaşmayı artırmıştır.

Bütün fiziksel-model kombinasyonlarında alevin parçalanma yolu esas olarak aynıdır. Toroidal girdap üst kapak ile aşağı kıvrılmış alev kolu arasındaki dar reaksiyon bölgesini gererek yerel sönmeye neden olmuştur. Soret difüzyonu parçalanmayı geciktirmiş; radyasyon ise öne çekmiştir. Her iki etki birlikte aktif olduğunda parçalanma 305 ms'de, ikisi de kapalı olduğunda 306 ms'de gerçekleşmiş ve 1 ms'lik fark çalışmanın zaman belirleme hassasiyeti içinde kaldığından anlamlı fiziksel kayma olarak yorumlanmamıştır.

Çalışma yeni bir fiziksel yanma deneyi değildir. Sonuçlar bir boyutlu küresel, iki boyutlu düzlemsel ve iki boyutlu eksensimetrik sayısal modellere dayanmaktadır. Eksensimetrik sonuçlar daha önce yayımlanmış deneysel önder-nokta yörüngeleriyle karşılaştırılmıştır; tam üç boyutlu asimetrik parçalanma ve türbülanslı alev davranışı bu çalışmada çözülmemiştir.

Araştırma problemi

Hidrojen–hava karışımları yanabilirlik sınırına yaklaştığında klasik, düzgün ve yaklaşık düzlemsel premiks alev davranışından önemli ölçüde uzaklaşabilir. Kimyasal ısı açığa çıkışı yavaşladıkça moleküler taşınım, ısı kaybı, eğrilik, alev gerinimi ve doğal konveksiyonun göreli önemi artar. Moleküler hidrojenin yüksek difüzivitesi nedeniyle hidrojen, sıcak reaksiyon bölgesine diğer birçok türden daha hızlı taşınabilir.

Karasal yerçekimi altında buna kaldırma kuvveti de eklenir. Ateşleme sonrası oluşan sıcak ve düşük yoğunluklu ürün çekirdeği yukarı doğru hareket eder. Bu hareket yalnız alevin yerini değiştirmez; aynı zamanda alevin çevresinde büyük ölçekli bir girdap üretir. Dolayısıyla gözlenen yukarı yönlü alev hareketi, yalnız bir laminer yanma hızının sonucu değildir. Kimyasal ön cep hareketi, sıcak gaz kütlesinin yükselmesi, radyal ve eksenel konveksiyon ve sürekli şekil değişimi birlikte rol oynar.

Kaynağın araştırdığı temel boşluk, Soret termal difüzyonu ile radyatif ısı kaybının serbestçe yükselen ultra-zayıf bir hidrojen alev çekirdeği üzerindeki ortak etkisinin karasal yerçekimi altında sistematik olarak nicelendirilmemiş olmasıdır. Önceki çalışmalar bu mekanizmaları çoğunlukla mikrogravite alev topları, duvarla sınırlı alevler, brülörle stabilize edilmiş alevler veya doğal konveksiyonsuz hücresel alevler bağlamında ele almıştır.

Soret Difüzyonu Ultra-Zayıf Hidrojen Alevinde Ne Yapar?

Soret difüzyonu, sıcaklık gradyanının hafif H ve H₂ türlerinde ek bir kütle akısı oluşturmasına neden olur; çalışmanın ultra-zayıf hidrojen karışımında bu mekanizma hidrojenin sıcak reaksiyon bölgesine taşınımını artırır, çekirdeğin büyümesini hızlandırır, daha fazla sıcak ürün oluşmasını destekler ve girdap kaynaklı gerinime karşı bağlantılı alev yapısının daha uzun süre korunmasına katkı sağlar.

Soret etkisi ve sıradan moleküler difüzyon

Kaynakta karışım-ortalama tür akısının düzeltilmemiş biçimi şu şekilde verilmiştir:

\[ \mathbf{j}_k^{(0)} = -\rho D_k \nabla Y_k - D_k^T \nabla \ln T \tag{10} \]

Burada \(\mathbf{j}_k^{(0)}\), \(k\) türünün düzeltilmemiş difüzif kütle akısıdır; \(\rho\) karışım yoğunluğu, \(D_k\) karışım-ortalama moleküler difüzyon katsayısı, \(Y_k\) kütle kesri, \(D_k^T\) termal-difüzyon katsayısı ve \(T\) sıcaklıktır.

Birinci terim tür derişimindeki uzaysal değişimden kaynaklanan normal difüzyonu, ikinci terim ise sıcaklık gradyanının oluşturduğu Soret taşınımını temsil eder. Soret etkisi kapatıldığında bütün \(D_k^T\) değerleri sıfıra ayarlanmıştır.

Toplam difüzif kütle akısının sıfır olması için düzeltme hızı uygulanmıştır:

\[ \mathbf{j}_k= \mathbf{j}_k^{(0)} - Y_k\sum_{l=1}^{N_s}\mathbf{j}_l^{(0)}, \qquad \sum_{k=1}^{N_s}\mathbf{j}_k=0. \tag{11} \]

Termal difüzyon yalnız H ve H₂ için Schlup–Blanquart indirgenmiş modeliyle uygulanmıştır. Kaynaktaki katsayılar \(\alpha_H=0{,}895\) ve \(\alpha_{H_2}=0{,}910\)'dur. Daha ağır türlerin doğrudan termal-difüzyon etkisi ihmal edilmiştir. Karşılıklı Dufour ısı akısı ise premiks hidrojen alevlerinde Soret katkısına göre daha küçük kabul edildiğinden modele eklenmemiştir.

Radyatif ısı kaybı nasıl modellenmiştir?

Radyasyon optik olarak ince emisyon yaklaşımıyla temsil edilmiştir ve radyasyon yapan tek tür su buharıdır. Hacimsel radyatif kaynak:

\[ \dot q_{\mathrm{rad}} = -4\sigma_{\mathrm{SB}} \left(T^4-T_b^4\right)\kappa_P, \qquad T_b=300\ \mathrm{K} \tag{12} \]

biçimindedir. \(\sigma_{\mathrm{SB}}\) Stefan–Boltzmann sabiti, \(T_b\) arka plan radyasyon sıcaklığı ve \(\kappa_P\) Planck-ortalama soğurma katsayısıdır.

Planck-ortalama katsayısı kaynakta Barlow ve arkadaşlarının polinom korelasyonuna göre hesaplanmıştır:

\[ \kappa_P = X_{H_2O} \frac{p_0}{p_{\mathrm{ref}}} \sum_{n=0}^{5} a_n \left(\frac{1000}{T}\right)^n, \qquad p_{\mathrm{ref}}=101325\ \mathrm{Pa}. \tag{13} \]

Optik olarak ince yaklaşım ürünlerin kendi radyasyonlarını yeniden soğurmasını ihmal eder. Bu ihmal optik kalınlık arttığında net ısı kaybını yüksek tahmin etme eğilimindedir. Dolayısıyla çalışma ayrıntılı spektral radyasyon aktarımını değil, santimetre ölçekli bu alev çekirdeğinde radyatif soğumanın birincil etkisini çözmektedir.

Başlangıç karışımı ve ateşleme

Küresel referans probleminde merkezdeki 0,4 mm yarıçaplı bölge anlık olarak 2000 K'ye çıkarılmış, ancak başlangıç kimyasal bileşimi değiştirilmemiştir. Böylece ilk aşamadaki çekirdek büyümesi, yerçekimi kaynaklı şekil değişiminden ayrılarak incelenmiştir.

İki boyutlu problemlerde ateşleme merkezi 300 K'lik izotermal alt duvarın 6,4 mm üzerinde yer almaktadır. Hesap alanının yarı genişliği veya eksensimetrik yarıçapı 51,2 mm, yüksekliği 102,4 mm'dir. Üst ve dış sınırlar serbest çıkış, eksen ise simetri sınırıdır.

Dört fiziksel model kombinasyonu

Bu matris, iki fiziksel mekanizmanın ayrı etkilerinin ve birlikte çalıştıklarındaki sonucun aynı başlangıç karışımında karşılaştırılmasına izin vermektedir.

Düşük-Mach matematiksel model

Reaktif karışım termal olarak mükemmel ideal gaz olarak ele alınmıştır. Düşük-Mach yaklaşımında basınç, uzaysal olarak yaklaşık düzgün termodinamik bölüm ile hidrodinamik pertürbasyona ayrılır:

\[ p(\mathbf{x},t)=p_0(t)+\pi(\mathbf{x},t). \tag{1} \]

Bu çalışmada \(p_0=101325\) Pa sabit tutulmuştur.

Kütlenin korunumu:

\[ \frac{\partial \rho}{\partial t} +\nabla\cdot(\rho\mathbf{u})=0. \tag{2} \]

Türlerin korunumu:

\[ \frac{\partial(\rho Y_k)}{\partial t} +\nabla\cdot(\rho Y_k\mathbf{u}) = -\nabla\cdot\mathbf{j}_k+\dot\omega_k, \qquad k=1,\ldots,N_s. \tag{3} \]

Momentum denklemi:

\[ \frac{\partial\mathbf{u}}{\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u} = -\frac{1}{\rho}\nabla\pi + \frac{1}{\rho}\nabla\cdot\boldsymbol{\tau} + \frac{\rho-\rho_0}{\rho}\mathbf{g}. \tag{4} \]

Son terim yoğunluk farkından doğan kaldırma etkisini temsil eder. Hidrostatik basınç bileşeni basınç alanından çıkarıldığı için referans yoğunluğu \(\rho_0\) bu terimde görünür.

Duyulur entalpi denklemi ısı iletimini, türlerin entalpi taşınımını, kimyasal ısı açığa çıkışını, viskoz terimi ve radyasyonu aynı enerji denkleminde birleştirir:

\[ \frac{\partial(\rho h_s)}{\partial t} +\nabla\cdot(\rho h_s\mathbf{u}) = \frac{dp_0}{dt} +\nabla\cdot(\lambda\nabla T) -\nabla\cdot \left( \sum_{k=1}^{N_s}h_{s,k}\mathbf{j}_k \right) -\sum_{k=1}^{N_s}\Delta h_{f,k}^{0}\dot\omega_k +\boldsymbol{\tau}:\nabla\mathbf{u} +\dot q_{\mathrm{rad}}. \tag{5} \]

Hidrojen oksidasyonu için Kéromnès ve arkadaşlarının ayrıntılı kimyasal mekanizması kullanılmış, azot inert tür olarak ele alınmıştır. Sıcaklığa bağlı ısı kapasiteleri ve entalpiler JANAF termokimyasal verilerine dayanan polinomlarla temsil edilmiştir.

Sayısal çözücü ve grid

Denklemler, düşük-Mach FDS algoritmasının kurum içi NRG uygulamasıyla çözülmüştür. Yoğunluk ve türler kimya ve taşınım kaynaklarıyla birlikte predictor–corrector projection şemasıyla ilerletilmiş; durum denklemi gerekli hız diverjansını sağlamış; hidrodinamik basınç için değişken katsayılı Poisson denklemi multigrid yöntemiyle çözülmüştür.

Akış laminer kabul edilmiş ve türbülans veya alt-grid yanma modeli kullanılmamıştır. Zaman adımı adaptiftir ve düşük-Mach çözücünün sayısal kararlılık ölçütleriyle sınırlandırılmıştır.

İki boyutlu hesaplarda temel hücre aralığı 0,4 mm'dir ve çözülen yarı alan 128 × 256 hücre içerir. Soret açık, radyasyon kapalı durum en ince alev cephesini ürettiği için grid yakınsama testi bu durumda yapılmıştır. 0,1, 0,2 ve 0,4 mm gridlerde önder-nokta yörüngesi ile 50, 100 ve 150 ms'deki \(T=1000\) K izotermleri neredeyse çakışmış ve sonraki hesaplarda 0,4 mm grid kullanılmıştır.

\(R_f^2 \simeq Kt\) Rejimi Ne Anlama Geliyor?

\(R_f^2 \simeq Kt\) rejimi, incelenen erken ateşleme-sonrası aralıkta küresel alev çekirdeğinin yarıçap karesinin yaklaşık sabit eğimle büyüdüğünü ve dolayısıyla ön cep hızının sabit olmadığını ifade eder; \(dR_f/dt=K/(2R_f)\) olduğundan çekirdek büyüdükçe radyal ilerleme hızı yaklaşık \(R_f^{-1}\) ile azalır. Kaynak bu ilişkiyi ilk yaklaşık 200 ms için etkili bir difüzyon-kontrollü rejim olarak yorumlar ve uzun zaman için evrensel asimptotik yasa olarak sunmaz.

Etkin büyüme katsayısı K

Alev cephesinin fit edilen davranışı:

\[ R_f^2(t) = R_f^2(t^*) + K(t-t^*) \tag{14} \]

şeklindedir. Türevi:

\[ U_f = \frac{dR_f}{dt} = \frac{K}{2R_f}. \tag{15} \]

Bu nedenle \(K\), moleküler bir difüzyon katsayısıyla veya laminer yanma hızıyla özdeş değildir. Kaynakta \(K\); tercihli difüzyon, ısı iletimi, sonlu hızlı kimya, termal genleşme, Soret taşınımı ve radyatif kaybın toplam etkisini içeren etkin çekirdek-büyüme katsayısı olarak tanımlanır.

Tablo, Soret difüzyonunun hem radyasyonlu hem radyasyonsuz durumda \(K\)'yı yükselttiğini, radyasyonun ise her iki Soret durumunda \(K\)'yı düşürdüğünü göstermektedir.

Neden \(R_f^2\) yaklaşık doğrusal büyüyor?

Kaynağın ölçekleme argümanında eksik reaktanın küresel reaksiyon bölgesine difüzif taşınımı yaklaşık:

\[ \dot m_D \sim 4\pi R_f^2 \rho D_{\mathrm{eff}} \frac{\Delta Y}{R_f} \propto R_f \tag{16} \]

şeklinde ölçeklenir. Küresel ön cephenin büyümesi için gerekli tüketim ise:

\[ \dot m_F \sim 4\pi R_f^2\rho_u\dot R_f. \tag{17} \]

İki ölçeğin dengelenmesi \(\dot R_f\propto R_f^{-1}\) sonucunu ve dolayısıyla \(d(R_f^2)/dt\) değerinin yaklaşık sabit olmasını açıklar. Bu, kimyanın önemsiz olduğu anlamına gelmez; kimyasal reaksiyon, difüzyonu besleyen sıcaklık ve derişim gradyanlarını üretmeye devam eder.

Sıcak ürün hacmi ve kaldırma kuvveti

Küresel sıcak bölgenin hacmi:

\[ V_h=\frac{4}{3}\pi R_f^3. \tag{18} \]

\(R_f^2\sim Kt\) olduğunda hacim yaklaşık:

\[ V_h\sim\frac{4\pi}{3}(Kt)^{3/2} \tag{19} \]

olarak ölçeklenir. Bu nedenle \(K\)'daki görece orta bir değişiklik, sonraki kaldırma hareketini besleyen sıcak gaz hacminde daha büyük bir göreli değişim oluşturabilir.

Kaynağın bütünleşik kaldırma kuvveti ifadesi:

\[ F_b = g\int_{\Omega_h}(\rho_u-\rho)\,dV \approx g(\rho_u-\rho_h)V_h \tag{20} \]

şeklindedir. Böylece erken küresel aşamada Soret veya radyasyon tarafından değiştirilen sıcak ürün hacmi, daha sonraki yükselme davranışına taşınır.

Radyasyon 200 ms'de neyi değiştirdi?

Kaynakta 200 ms'deki radyal profiller, radyasyonun ürün sıcaklığını yaklaşık 61 K azalttığını göstermektedir. Bu düşüş, ortam üzerindeki sıcaklık artışının yaklaşık %7'sidir. Buna karşılık ürün yoğunluğu yaklaşık 0,015 kg/m³ yükselmiş ve bu artış ortam karışımı yoğunluğunun yaklaşık %1,89'u kadar olmuştur.

Düşük-Mach yaklaşımında yaklaşık sabit termodinamik basınç altında yoğunluk:

\[ \rho=\frac{p_0M}{R_uT} \tag{21} \]

ile verilir. Küçük değişimler için:

\[ \frac{\delta\rho}{\rho} \approx \frac{\delta M}{M} - \frac{\delta T}{T}. \tag{22} \]

Bu ifade, sıcaklık düşüşünün yoğunluk üzerindeki etkisinin bileşim ve ortalama molar kütledeki değişim tarafından kısmen dengelenebileceğini göstermektedir.

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

Analiz zinciri

1. %6 hacimsel H₂–hava karışımı için dört Soret/radyasyon kombinasyonu tanımlandı.
2. Bir boyutlu küresel hesaplarla yerçekiminden bağımsız erken çekirdek büyümesi incelendi.
3. \(R_f^2(t)\) eğrileri üzerinden etkin büyüme katsayısı \(K\) belirlendi.
4. Sıcaklık, yoğunluk, hız ve H₂/O₂/H₂O profilleriyle Soret ve radyasyonun fiziksel etkileri ayrıştırıldı.
5. İki boyutlu düzlemsel ve eksensimetrik hesaplarla doğal konveksiyon ve yükselme çözüldü.
6. Önder-nokta yörüngeleri daha önce yayımlanmış iki deneysel veri setiyle karşılaştırıldı.
7. Sıcaklık ve OH alanları kullanılarak sıcak ürün bölgesi ile aktif reaksiyon tabakası ayrıldı.
8. Üst alev kapağı elipsle fit edilerek biçim oranı \(a_r/a_z\) izlendi.
9. Toroidal girdap ile reaksiyon bölgesi arasındaki etkileşim ve ilk bağlantı kopması üzerinden parçalanma zamanı belirlendi.

Yerçekimi ne zaman önemli hale geliyor?

Erken evrede çekirdek yaklaşık küresel ve difüzyon kontrollüdür. Kaynak, alevin büyüme süresi ile kaldırmayla yükselme süresini karşılaştıran bir ölçüt kullanmıştır. Davies–Taylor tabanlı yükselme hızı:

\[ V_{\mathrm{rise}} = \frac{2}{3}\sqrt{g'R}, \qquad g'=g\left(1-\frac{\rho_b}{\rho_u}\right) \tag{23} \]

şeklinde yazılır.

Bu çalışmanın yarıçapa bağlı büyüme hızı kullanıldığında zaman ölçeği oranı:

\[ \Pi_{bK} = \frac{t_{\mathrm{growth}}}{t_{\mathrm{rise}}} = \frac{V_{\mathrm{rise}}}{\dot R_f} = \frac{4}{3} \sqrt{\frac{g'R_f^3}{K^2}} \tag{24} \]

olur. \(\Pi_{bK}\ll1\) difüzyon kontrollü yaklaşık izotropik büyümeyi, \(\Pi_{bK}=O(1)\) ise belirgin kaldırma kaynaklı hareket ve deformasyon başlangıcını temsil eder.

Buna karşılık gelen karakteristik yarıçap:

\[ R_{bK}^{*} = \left( \frac{9K^2}{16g'} \right)^{1/3}. \tag{25} \]

Hem Soret hem radyasyonun açık olduğu durumda \(K=1{,}113\times10^{-3}\) m²/s ve kaynakta kullanılan temsili yoğunluklarla \(R_{bK}^{*}\) yaklaşık 4,6 mm bulunmuştur. Bu değer hassas bir deneysel eşik değil, kaldırmanın önem kazanmaya başladığı ölçek için mertebe tahminidir.

Düzlemsel ve eksensimetrik sonuçların farkı

Her iki model de çekirdeğin zamanla hızlanarak yükseldiğini ve dört Soret/radyasyon kombinasyonunun aynı göreli sıralamasını göstermiştir. Ancak düzlemsel model önder-nokta yer değiştirmesini sistematik biçimde daha düşük tahmin etmiştir. Eksensimetrik hesaplar, özellikle büyük hacimli Volodin ve arkadaşları deneyinin yörüngesine çok daha yakın sonuç vermiştir.

Bu durum eksensimetrik modelin “tam üç boyutlu doğru çözüm” olduğu anlamına gelmez. Gerçek deneysel alev üç boyutlu bir nesnedir. Kaynak, eksensimetrik geometrinin mevcut deney karşılaştırması için daha iyi yaklaşım olduğunu, gerçek 3B etkileri kaldırma, ek kütle, sürükleme ve şekil katkılarına ayırmanın ise tam üç boyutlu simülasyon gerektirdiğini belirtmektedir.

250 ms'deki sıcaklık alanı

Kaynağın Şekil 6'sında bütün dört durumda başlangıçtaki kompakt çekirdek, 250 ms'de üst kısmı geniş bir sıcak kapak ve aşağı doğru uzanan wake yapısına dönüşmüştür. Toroidal girdap üst sınırı yukarı ve yana taşırken alt sınırı içeri doğru çökerterek kapak biçimini oluşturmuştur.

Soret açık durumlar daha büyük ve daha yukarı taşınmış sıcak bölgeler üretmiştir. Radyasyon ürün sıcaklığını ve çekirdek boyutunu bir miktar düşürmüştür; ancak erken ve orta yükselme aşamasında radyasyonun yörünge üzerindeki etkisi Soret etkisinden daha küçüktür.

OH alanı neden önemli?

Sıcaklık alanı sıcak ürünlerin geniş bölgesini gösterirken OH kütle kesri aktif reaksiyon bölgesini daha seçici biçimde ortaya koymaktadır. Kaynağın Şekil 7'sinde önemli OH esas olarak üst yüzey ve iki yan omuzda bulunmakta, alt sıcak sınırın önemli bölümü ise kimyasal olarak zayıf veya etkin değildir.

Bu sonuç, yükselen yapının her yönde düzgün biçimde ilerleyen kapalı küresel bir alev cephesi olmadığını göstermektedir. Alt eksen boyunca girdabın çektiği soğuk taze karışım, alt reaksiyon bölgesini soğutur ve gerer; böylece aktif üst kapak ile büyük ölçüde inaktif alt wake arasında açık bir yapı gelişir.

Üst kapak nasıl nicelleştirildi?

Üst alev sınırı eksene hizalanmış bir elipsle fit edilmiştir:

\[ \frac{(x-x_c)^2}{a_r^2} + \frac{(z-z_c)^2}{a_z^2} = 1. \tag{26} \]

Eksensimetrik durumda \(x\) yerine radyal koordinat \(r\) kullanılır. \(a_r\) yanal, \(a_z\) düşey yarı eksendir.

Biçim oranı:

\[ \chi=\frac{a_r}{a_z}. \tag{27} \]

\(\chi=1\) meridyen düzleminde dairesel üst kapağı, \(\chi>1\) ise yatay yönde genişlemiş basık kapağı gösterir.

Kaynakta biçim oranı zaman yerine önder-nokta yüksekliği \(z_{lp}\) ile karşılaştırılmıştır. Böylece farklı yükselme hızlarına sahip durumlar benzer düşey gelişim seviyelerinde kıyaslanmıştır. Genel olarak çekirdek yükseldikçe \(\chi\) artmış ve küresel yaklaşım giderek daha az uygun hale gelmiştir.

Biçim oranının zamansal değişimi:

\[ \frac{d}{dt} \ln\left(\frac{a_r}{a_z}\right) = \frac{\dot a_r}{a_r} - \frac{\dot a_z}{a_z}. \tag{28} \]

Dolayısıyla biçim oranının artması, göreli yanal büyümenin göreli düşey büyümeden daha hızlı olduğunu ifade eder.

Soret ve radyasyon şekli aynı yönde mi etkiliyor?

Hayır. Aynı önder-nokta yüksekliğinde radyasyon genellikle \(a_r/a_z\) oranını artırmış ve kapağın daha basık olmasına yol açmıştır. Soret difüzyonu ise çekirdeğin toplam boyutunu ve yükselmesini artırmasına rağmen biçim oranını genellikle azaltmış, üst kapağı daha dairesel tutmuştur.

Bu sonuç önemli bir ayrımı gösterir: yalnız önder-nokta yüksekliğine veya hızına bakmak alev morfolojisini tam açıklamaz. Aynı yüksekliğe ulaşmış iki çekirdek farklı en-boy oranına, eğriliğe ve kimyasal aktivite dağılımına sahip olabilir.

Yükselen Hidrojen Alev Çekirdeği Neden Parçalanıyor?

Çekirdek, yükselen sıcak ürünlerin oluşturduğu toroidal girdabın üst aktif kapak ile girdap etrafına kıvrılmış alt alev kolu arasındaki dar reaksiyon bölgesini güçlü biçimde germesi nedeniyle parçalanır; artan eğrilik ve gerinim reaktanların sıcak bölgede kalma süresini azaltıp konvektif soğumayı artırdığında kimyasal ısı açığa çıkışı bu kayıpları karşılayamaz ve bağlantı bölgesinde yerel sönme meydana gelir.

Parçalanma zamanının tanımı

Parçalanma zamanı \(t_b\), daha önce sürekli olan ve üst kapak ile girdap tarafından sarılmış alt kolu birbirine bağlayan reaksiyon bölgesinin ilk kez bağlantısını kaybettiği an olarak belirlenmiştir. Zaman yaklaşık 1 ms hassasiyetle izlenmiştir.

Radyasyon kapalıyken Soret difüzyonu parçalanmayı 306 ms'den 327 ms'ye, yani 21 ms geciktirmiştir. Radyasyon açıkken aynı etki 282 ms'den 305 ms'ye, yani 23 ms gecikmedir. Kaynak bunu bağlantılı alev yapısının ömründe yaklaşık %7–8 uzama olarak ifade etmektedir.

Radyasyon ters yönde davranmıştır. Soret kapalıyken parçalanmayı 306 ms'den 282 ms'ye, 24 ms öne çekmiş; Soret açıkken 327 ms'den 305 ms'ye, 22 ms öne çekmiştir.

İki mekanizma aynı anda açık olduğunda 305 ms, ikisi de kapalı olduğunda 306 ms bulunmuştur. Bu 1 ms fark çalışma tarafından çözümlenmiş fiziksel değişiklik olarak kabul edilmemektedir. İki durum aynı parçalanma zamanına yakın olsa bile sıcaklık, bileşim, çekirdek boyutu ve yükselme geçmişleri aynı değildir.

Parçalanma tam sönme anlamına geliyor mu?

Hayır. İlk bağlantı koptuktan sonra hem üst kapakta hem de kıvrılmış alt kolda kimyasal enerji açığa çıkışı devam etmektedir. Başlangıçtaki tek çekirdek ortadan kaybolmak yerine birden fazla reaksiyon yapısına dönüşür. Kaynak bu parçaların daha sonra termal ve kimyasal wake içinde hareket edebileceğini, etkileşebileceğini veya ikincil çekirdekler oluşturabileceğini belirtmektedir.

Akış mı, Soret ve radyasyon mu parçalanma yolunu belirliyor?

Nitel parçalanma konumunu ve yolunu esas olarak hidrodinamik akış topolojisi belirlemektedir. Bütün dört modelde aynı temel toroidal girdap–alev bağlantısı görülmüştür. Soret difüzyonu ve radyasyon esas olarak alevin girdap kaynaklı gerinim ve soğumaya ne kadar dayanabileceğini ve kritik koşula ne zaman ulaşacağını değiştirmiştir.

Türbülans açısından ne söylenebilir?

Bu çalışma türbülanslı alev simülasyonu değildir. Akış laminer ve tek baskın toroidal girdap tarafından kontrol edilmektedir. Bununla birlikte yerel alev–girdap etkileşiminde yüzey deformasyonu, eğrilik, gerinim, reaktan kalış süresi, konvektif taşınım ve yerel sönme gibi türbülanslı alev–eddy etkileşimlerinde de bulunan fiziksel süreçler görülmektedir.

Bu nedenle çalışma, akış tarafından kontrol edilen bir yerel sönme eşiğinin Soret taşınımı ve ısı kaybıyla değişebileceğini göstermektedir. Ancak türbülanslı hidrojen yanabilirlik sınırlarının nicel olarak hesaplandığı sonucu çıkarılamaz; bunun için farklı boyut ve gerinim oranlarına sahip çok sayıda girdabı içeren üç boyutlu simülasyonlar gerekir.

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

• %6 hacimsel H₂–hava karışımındaki erken çekirdek büyümesi incelenen zaman aralığında yaklaşık \(R_f^2\simeq Kt\) rejimi göstermiştir.
• Soret difüzyonu \(K\)'yı yükseltmiş ve çekirdeğin büyümesini hızlandırmıştır.
• Radyatif ısı kaybı \(K\)'yı azaltmış ve ürünleri soğutmuştur.
• Soret etkisi daha büyük ve daha hızlı yükselen sıcak ürün çekirdeği oluşturmuştur.
• Eksensimetrik hesaplar mevcut deneysel önder-nokta yörüngesini düzlemsel hesaplardan daha iyi yeniden üretmiştir.
• Radyasyon aynı yüksekliğe göre üst kapağın daha basık olmasını desteklemiştir.
• Soret difüzyonu üst kapağın daha dairesel kalmasını desteklemiştir.
• Parçalanma bütün durumlarda toroidal girdabın oluşturduğu yerel gerinim ve sönmeyle gerçekleşmiştir.
• Soret parçalanmayı geciktirmiş, radyasyon öne çekmiştir.
• Soret ve radyasyonun parçalanma zamanı üzerindeki karşıt etkileri birlikte açık durumda neredeyse birbirini dengelemiştir.

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

• Çalışma %6 H₂ karışımının bütün gerçek tesis geometrilerindeki davranışını belirlemez.
• Hesaplanan önder-nokta hızı laminer yanma hızıyla doğrudan özdeş değildir.
• \(R_f^2\simeq Kt\) rejimi sonsuz zaman için evrensel büyüme yasası olarak kanıtlanmamıştır.
• Eksensimetrik model tam üç boyutlu alev parçalanmasını çözmez.
• Çalışma türbülanslı alev veya türbülanslı yanabilirlik sınırı hesaplaması değildir.
• Optik olarak ince radyasyon modeli ayrıntılı spektral radyatif transferi çözmez.
• Çalışmanın kendisinde yeni deney yapılmamıştır; deneysel veriler önceki yayınlardan karşılaştırma amacıyla alınmıştır.
• Parçalanma tam yanma sönmesi anlamına gelmez.

Temel yöntemsel sınırlılıklar

Düzlemsel ve eksensimetrik iki boyutlu geometriler, gerçek üç boyutlu alevin non-axisymmetric deformasyonlarını temsil edemez. Özellikle parçalanmadan sonra farklı alev parçalarının üç boyutlu etkileşimleri ve birden fazla girdabın eşzamanlı etkisi kapsam dışıdır.

Akışın laminer kalacağı varsayılmış ve türbülans ya da subgrid yanma modeli kullanılmamıştır. Dolayısıyla sonuçlar doğrudan türbülanslı hidrojen bulutlarına aktarılamaz.

Radyasyon yalnız H₂O emisyonu ve optik olarak ince yaklaşım altında modellenmiştir. Kendiliğinden soğurma hesaba katılmadığından net radyatif kayıp optik kalınlığın arttığı durumlarda fazla tahmin edilebilir.

Mevcut deneylerle karşılaştırma yalnız önder-nokta hareketi açısından yapılmaktadır ve deneylerin geometrileri de birbirinden farklıdır. Büyük hacimli deney yaklaşık serbest bir çekirdeği temsil ederken ikinci veri seti kapalı tüpten gelmektedir. Geometrik sınırlama ve duvar etkileşimleri özellikle daha geç evrelerde önem kazanabilir.

Kaynak ve Yöntem Notu

Özgün çalışma: Effects of Soret Diffusion and Radiative Heat Loss on the Evolution of Buoyant Flame Kernels in Ultra-Lean Hydrogen-Air Mixture.

Yazarlar: Ivan S. Yakovenko; Alexey D. Kiverin.

Sorumlu yazar: Ivan S. Yakovenko.

Kurum: Joint Institute for High Temperatures of the Russian Academy of Sciences, Moscow, Russia.

Dergi ve yayın kimliği: Fire, 2026, Cilt 9, Makale 383. Yayınevi: MDPI. DOI: 10.3390/fire9090383.

Yayın süreci: Makale 9 Ağustos 2026'da alınmış, 29 Ağustos 2026'da revize edilmiş, 4 Eylül 2026'da kabul edilmiş ve 5 Eylül 2026'da yayımlanmıştır.

Kaynak türü: Hakemli, açık erişimli, ayrıntılı sayısal yanma ve akışkanlar dinamiği araştırma makalesidir. Çalışmanın kendisi yeni fiziksel yanma deneyi içermemektedir.

Lisans: Creative Commons Attribution (CC BY).

Temel yöntem: %6 hacimsel H₂–hava karışımı; ayrıntılı Kéromnès hidrojen oksidasyon mekanizması; düşük-Mach gaz dinamiği; karışım-ortalama moleküler taşınım; H ve H₂ için indirgenmiş Soret modeli; optik olarak ince H₂O radyasyon modeli; 1B küresel, 2B düzlemsel ve 2B eksensimetrik simülasyonlar; kurum içi NRG düşük-Mach FDS çözücüsü.

Grid: Ana 2B hesaplarda 0,4 mm uniform hücre aralığı ve çözülen yarı alanda 128 × 256 hücre. Yakınsama kontrolü 0,1, 0,2 ve 0,4 mm gridlerle gerçekleştirilmiştir.

Deneysel doğrulama bağlamı: Sayısal önder-nokta yörüngeleri Volodin ve arkadaşları ile Anikin ve arkadaşlarının daha önce yayımlanmış deneysel verileriyle karşılaştırılmıştır. Deneyler bu araştırma kapsamında gerçekleştirilmemiştir.

Finansman: Çalışma Russian Federation Ministry of Science and Higher Education tarafından State Assignment No. 075-00270-26-00 kapsamında desteklenmiştir.

Veri erişilebilirliği: Yazarlar makalenin sonuçlarını destekleyen ham verilerin talep üzerine sağlanacağını belirtmiştir.

Çıkar çatışması: Yazarlar çıkar çatışması olmadığını beyan etmiştir.

Yazar katkıları: Kavramsallaştırma, yöntem, yazılım, doğrulama, araştırma, veri düzenleme, ilk taslak ve görselleştirme Ivan S. Yakovenko tarafından yürütülmüştür. Biçimsel analiz Ivan S. Yakovenko ve Alexey D. Kiverin tarafından; kaynaklar, gözden geçirme ve düzenleme, gözetim, proje yönetimi ve finansman temini iki yazar tarafından ortak yürütülmüştür.

Bilimsel yorum sınırı

Çalışmanın ana sonucu, ultra-zayıf hidrojen alevlerinin karasal yerçekimi altında tek bir “yanma hızı” ile yeterince açıklanamadığıdır. Moleküler taşınım erken büyümeyi, sıcak ürün hacmi kaldırma kuvvetini, kaldırma kendi girdap alanını ve girdap da alev morfolojisi ile yerel sönmeyi etkilemektedir. Soret difüzyonu ve radyasyon bu zincirin zaman ölçeğini ve alevin gerinime dayanımını değiştirir; ancak parçalanmanın temel topolojik yolu incelenen dört durumda da hidrodinamik olarak benzerdir.

Sonuçların gerçek üç boyutlu, non-axisymmetric veya türbülanslı hidrojen yangınlarına doğrudan aktarılması uygun değildir. Yazarların sonuç bölümünde de bu kapsam için tam üç boyutlu simülasyonların gerekli olduğu belirtilmektedir.