Bu çalışma, çift girdaplı bir gaz türbini enjektöründe metan-hidrojen birlikte yakma alevlerinin enjektöre tutunmuş veya enjektörden kalkmış biçimde kararlı hâle gelmesini deneysel ölçümler ve bir boyutlu karşı akış simülasyonlarıyla incelemiştir. Araştırmacılar, 36 çalışma koşulunda yüksek hızlı OH düzlemsel lazerle indüklenmiş floresans görüntülemesi, parçacık görüntülemeli hız ölçümü ve CHEMKIN tabanlı OPPDIF hesaplamaları kullanmıştır. Hidrojen oranının artırılması genel olarak alevin enjektöre tutunmasını desteklerken ana akış hızının 10 m/s’den 15 m/s’ye çıkarılması kalkmış alev bölgesini genişletmiştir. Alev genişliği ile pilot alev kalkma yüksekliği hidrojen oranıyla birlikte artmış ve özellikle pilot ile ana geçit arasındaki momentum oranıyla güçlü ilişkiler göstermiştir. Bununla birlikte çalışma hakem değerlendirmesinden geçmemiş bir preprinttir; bir boyutlu model üç boyutlu türbülanslı alevin mutlak boyutlarını değil, yalnızca değişim eğilimlerini temsil etmektedir.

Deneysel görüntüler iki temel alev yapısı göstermiştir. Tutunmuş alevde dış reaksiyon tabakası enjektör çıkışına kadar uzanırken kalkmış alevde bu tabakanın kökü yerel olarak sönmüş ve reaksiyon bölgesi daha aşağı akışta yeniden kararlı hâle gelmiştir. Kalkmış alevde zaman ortalamalı OH bölgesinin daha geniş görünmesi, tek tek alev tabakalarının fiziksel olarak kalınlaşmasından çok alev yüzeyinin dalgalanması, kırışması, aralıklı davranması ve iç-dış reaksiyon tabakalarının kısmen birleşmesiyle açıklanmıştır.

Araştırmanın temel sorusu nedir?

Araştırmanın temel sorusu, dış ana geçitten verilen ön karışımlı CH₄/hava akımı ile iç pilot geçitten oksitleyici olmadan verilen saf H₂ akımının çift girdaplı bir enjektörde nasıl etkileştiği ve bu etkileşimin alevi enjektöre tutunmuş ya da enjektörden kalkmış bir yapıya nasıl dönüştürdüğüdür. Çalışma ayrıca karmaşık üç boyutlu türbülanslı alevin başlıca reaksiyon tabakalarının, yerel birer alevcik gibi ele alınarak bir boyutlu karşı akış modeliyle açıklanıp açıklanamayacağını araştırmıştır.

Literatürde çift girdaplı alevlerin kararlılığı üzerine çalışmalar bulunmasına karşın ortalama OH-PLIF görüntülerinden alev genişliği ve pilot alev kalkma yüksekliği gibi tekrar üretilebilir geometrik göstergelerin çıkarılması, bu göstergelerin pilot-ana akış momentumu ile ilişkilendirilmesi ve aynı eğilimlerin bir boyutlu OPPDIF modeliyle karşılaştırılması sınırlı kalmıştır. Araştırma bu boşluğu deney, görüntü işleme ve indirgenmiş sayısal modellemeyi aynı çerçevede birleştirerek hedeflemiştir.

Çift girdaplı enjektör nasıl çalışmaktadır?

Yanma odası optik ölçümlere izin veren dört kuvars pencereli, 70 mm × 70 mm iç kesitli ve 260 mm yüksekliğinde kare bir yapıdan oluşmaktadır. Dış halkasal ana geçitten %99,999 saflıkta metanla hazırlanan ön karışımlı CH₄/hava akımı, iç halkasal pilot geçitten ise %99,999 saflıkta hidrojen gönderilmiştir. Pilot hidrojen akımına girişte hava veya başka bir oksitleyici eklenmemiştir.

Ana girdap oluşturucu teğetsel, pilot girdap oluşturucu ise altı kanatlı eksenel tiptedir. Her iki akış saat yönünün tersine döndürüldüğü için düzenek eş yönlü çift girdaplı yapı oluşturmuştur. Ana ve pilot girdap sayıları sırasıyla 0,51 ve 0,83 olarak verilmiştir. Geometrik şemada ana çıkış çapı 20,0 mm, pilot bölümle ilişkili çaplar ise 10,7 mm ve 7,0 mm olarak gösterilmektedir.

Girdaplı akış, yanmış sıcak gazların merkez bölgesinden enjektöre doğru geri taşındığı bir iç devridaim bölgesi oluşturur. Bu geri dönen sıcak gaz, taze yakıt-hava karışımının yeniden tutuşmasına yardım eden sürekli bir ısı ve aktif tür kaynağı gibi davranmaktadır. Ancak yalnızca bir iç devridaim bölgesinin bulunması alevin mutlaka enjektöre tutunacağı anlamına gelmemiştir. Devridaim bölgesinin enjektöre ne kadar yaklaştığı, ana akışın eksenel hızı, kesme tabakasındaki hız gradyanı ve pilot hidrojen jetinin momentumu birlikte alev kökünün konumunu belirlemiştir.

Tutunmuş ve kalkmış alev arasındaki fark nedir?

Çalışmanın önemli yorumlarından biri, kalkmış alevin bütün alev yapısının tek parça hâlinde yukarı taşınması olarak değerlendirilmemesidir. İlk aşamada dış reaksiyon tabakasının enjektöre yakın kökünde yerel reaksiyon-taşınım dengesi bozulmakta, yerel sönme meydana gelmekte ve alev daha uygun hız ve karışım koşullarının bulunduğu aşağı akış konumunda yeniden tutunmaktadır.

Kalkmış alevde ortalama OH bölgesinin daha kalın görünmesi de tek bir anlık alev cephesinin iç yapısının kalınlaştığını kanıtlamamaktadır. Saniyede 10.000 kareyle alınan anlık görüntüler, alev yüzeyinin zaman içinde farklı konumlara hareket ettiğini, daha güçlü biçimde kırıştığını ve iç-dış tabakaların zaman zaman üst üste geldiğini göstermiştir. Dört bin görüntünün ortalanması bu hareketlerin tamamını geniş bir reaksiyon bölgesi biçiminde göstermektedir.

Hidrojen oranı alev kararlılığını nasıl değiştirmiştir?

Hidrojen oranının yükselmesi genel olarak alevin enjektöre tutunmasını desteklemiştir. Hidrojen oranı arttığında yalnızca karışımın kimyasal reaktivitesi değişmemiş, pilot nozuldan çıkan hidrojenin hızı ve momentumu da artmıştır. Böylece merkezdeki pilot jet, dış CH₄/hava alev kökünü enjektör yakınında destekleyen daha güçlü bir aerodinamik ve kimyasal etki oluşturmuştur.

Örneğin ana akış hızı 10 m/s ve küresel eşdeğerlik oranı 0,6 iken %60 hidrojen içeriğinde tutunmuş, %50 hidrojen içeriğinde ise kalkmış alev gösterilmiştir. Hidrojen oranı %60’tan %50’ye düşürüldüğünde ana CH₄/hava karışımının yerel eşdeğerlik oranı yükselmiş ve bu değişim tek başına değerlendirildiğinde yanma hızına olumlu etki yapabilecek bir yön oluşturmuştur. Buna rağmen alevin kalkması, geçişin yalnızca ana alev kimyasına bağlı olmadığını; pilot hidrojen jetinin zayıflayan hızı, momentumu ve tutuşturma desteğinin daha belirleyici olduğunu göstermiştir.

Kararlılık haritasında yalnızca metan kullanılan durumda yaklaşık Φ_g ≤ 0,5 bölgesinde fakir karışım sönmesi gözlenmiştir. Küçük miktarda hidrojen eklendiğinde ise alev daha düşük küresel eşdeğerlik oranlarında sürdürülebilmiştir. Bu bulgu, pilot hidrojen akımının merkez bölgede ısı ve aktif radikal desteği sağlayarak fakir yanma sınırını genişletebildiğini göstermektedir. Çalışma bununla birlikte emisyon, yakıt tüketimi veya sistem verimi ölçmediği için bu kararlılık artışı doğrudan çevresel ya da ekonomik kazanç olarak yorumlanmamalıdır.

Ana akış hızının etkisi nedir?

Ana geçitteki yığın hızının 10 m/s’den 15 m/s’ye çıkarılması kalkmış alev bölgesini belirgin biçimde genişletmiştir. OH-PLIF ölçümü yapılan ve ana hızın 15 m/s olduğu %20-%60 hidrojen, Φ_g = 0,5-0,7 koşullarının tamamında dış alev kökü enjektörden ayrılmış ve kalkmış alev morfolojisi korunmuştur.

Yüksek ana hız, ön karışımlı akımın eksenel taşınımını ve ana kesme tabakasındaki hız gradyanını artırmıştır. Böylece alev yüzeyinin enjektör yakınında kalabilmesi için gereken yerel yanma ve yeniden tutuşma koşulları daha zor hâle gelmiştir. Hidrojen oranının yükselmesi pilot desteğini artırsa da 15 m/s koşulunda bu destek, ana akımın baskın taşınım etkisini yenerek alevi tekrar enjektöre bağlamaya yetmemiştir.

Küresel eşdeğerlik oranı neden tek başına yeterli değildir?

Küresel eşdeğerlik oranı Φ_g, sistemdeki toplam yakıtın toplam havaya göre stokiyometrik durumunu ifade etmektedir. Ancak deney düzeneğinde Φ_g değiştirildiğinde ana geçidin eşdeğerlik oranı, toplam ısıl güç, pilot hızı, hız oranı ve pilot-ana etkileşimi de birlikte değişmiştir. Bu nedenle daha yüksek Φ_g değerinin her zaman alevi enjektöre daha güçlü bağlayacağı biçiminde basit bir sonuç çıkarılamamıştır.

Görüntülerde Φ_g yükseldikçe reaksiyon bölgesinin daha geniş bir alana yayıldığı, ısı açığa çıkışı ve termal genleşmeyle birlikte alevin aşağı akışta daha gelişmiş bir yapı oluşturduğu görülmüştür. Bazı koşullarda daha zengin ana karışıma rağmen kalkmış morfolojinin sürmesi, alev kararlılığının yalnızca kimyasal yanma hızına değil, pilot ve ana akımlar arasındaki aerodinamik dengeye bağlı olduğunu göstermiştir.

OH-PLIF görüntüleri nasıl sayısallaştırılmıştır?

Araştırmacılar her koşul için 4.000 anlık OH-PLIF görüntüsünü ortalamış, lazer tabakası yoğunluk düzeltmesini uygulamış ve alev ile arka planı ayırmak için OTSU eşikleme yöntemini kullanmıştır. OTSU yöntemi, görüntünün parlaklık histogramını iki sınıfa ayıran ve sınıflar arası varyansı en yüksek yapan eşiği otomatik seçmektedir. Böylece araştırmacı tarafından elle belirlenecek sabit bir eşikten doğabilecek öznellik azaltılmıştır.

Farklı koşullarda seçilen normalize OTSU eşikleri maksimum OH yoğunluğunun yaklaşık 0,20-0,30’u arasında dağılmış, ortalama değer yaklaşık 0,25 olmuştur. Bu sonuç daha sonra sayısal modelde de kullanılmış ve maksimum OH mol kesrinin %25’ine karşılık gelen konum karşılaştırma sınırı olarak seçilmiştir. Bu %25 düzeyi gerçek alev yüzeyinin mutlak fiziksel sınırı değildir; deneysel ve sayısal şekil göstergelerini aynı referansla karşılaştırmaya yarayan ampirik bir ölçüttür.

İki geometrik gösterge tanımlanmıştır:

• Alev genişliği: Enjektör merkez çizgisi ile dış reaksiyon tabakasının OTSU sınırı arasındaki radyal mesafe.
• Pilot alev kalkma yüksekliği: Nozul çıkış düzlemi ile iç reaksiyon tabakasının tabanı arasındaki eksenel mesafe.

Grafiklerden yaklaşık olarak okunduğunda ana hızın 10 m/s olduğu koşullarda alev genişliği %20 hidrojen düzeyinde yaklaşık 7,3-7,8 mm iken %80 hidrojen düzeyinde Φ_g değerine bağlı olarak yaklaşık 14-17,7 mm’ye ulaşmıştır. Pilot alev kalkma yüksekliği aynı hidrojen aralığında yaklaşık 3,4-4,6 mm’den yaklaşık 14,8-18,8 mm’ye yükselmiştir. Bunlar grafik üzerinden yaklaşık okunan değerlerdir ve kesin ham veri olarak değerlendirilmemelidir.

Ana hızın 15 m/s olduğu ve bütün ölçüm koşullarının kalkmış alev oluşturduğu durumda genişlik artışı daha sınırlı kalmıştır. Alev genişliği %20 hidrojen düzeyinde yaklaşık 6,8-7,7 mm, %60 hidrojen düzeyinde yaklaşık 8,3-10,5 mm aralığındadır. Pilot kalkma yüksekliği ise Φ_g değerine bağlı olarak daha farklı mutlak düzeyler göstermiştir. Bu durum kalkma yüksekliğinin yalnızca tutunmuş-kalkmış geçişini değil, kalkmış rejim içindeki iç reaksiyon tabakası konumunu da temsil ettiğini göstermektedir.

Momentum oranı neden anahtar parametre olarak değerlendirilmiştir?

Hidrojen oranı, ana hız ve eşdeğerlik oranı birbirini etkileyen işletme değişkenleridir. Araştırmacılar bu değişkenlerin aerodinamik sonucunu tek bir büyüklükte toplamak için pilot ve ana geçitlerin momentum akısı oranını kullanmıştır:

\[ J = \frac{\rho_P U_P^2}{\rho_M U_M^2} \]

Burada J boyutsuz momentum oranını, ρ_P ve ρ_M pilot ve ana geçitlerdeki karışım yoğunluklarını, U_P ve U_M ise nozul çıkış hızlarını göstermektedir. Momentum akısı yoğunluk ile hızın karesinin çarpımına bağlı olduğu için hızdaki bir değişim, pilot ve ana akımlar arasındaki dengeyi güçlü biçimde değiştirmektedir.

Grafiklerin yatay ekseninde yazarların kullandığı J^0,16 - 0,23 dönüşümü gösterilmiştir. Deneysel sonuçlarda elde edilen doğrusal uyumlar şöyledir:

Alev genişliğinin kalkmış rejimde daha düşük R² göstermesi, tek parametrenin bu rejimdeki bütün morfolojiyi açıklayamadığını göstermektedir. Dış alev kökü enjektörden ayrıldıktan sonra yerel karışım, ters akış, kesme tabakası gerinimi, kısmi ön karışım derecesi ve zamansal dalgalanma daha güçlü roller üstlenmektedir. Buna karşılık pilot kalkma yüksekliği, merkez jetin devridaim bölgesiyle nerede karşılaştığını doğrudan yansıttığı için momentum oranına daha düzenli tepki vermiştir.

Bir boyutlu OPPDIF modeli neyi temsil etmektedir?

Sayısal model bütün yanma odasını veya üç boyutlu girdaplı türbülansı çözmemiştir. Araştırmacılar OH-PLIF görüntülerinde görülen iki temel reaksiyon tabakasını ayrı yerel karşı akış alevleri olarak idealize etmiştir.

• Durum I: Dış reaksiyon tabakasını temsil etmektedir. Bir taraftan ön karışımlı CH₄/hava, karşı taraftan saf H₂ gönderilmiştir. Her iki giriş sıcaklığı 298 K, CH₄/hava giriş hızı 10 m/s, hidrojen hızı ise deney koşuluna bağlı değişken olarak tanımlanmıştır.
• Durum II: İç pilot reaksiyon tabakasını temsil etmektedir. Alt taraftan 298 K sıcaklığında H₂, karşı taraftan Durum I hesabından elde edilen yanmış gaz bileşimi gönderilmiştir. Yanmış gaz sıcaklığı ilgili koşulun adyabatik alev sıcaklığına eşitlenmiştir.

Her iki modelde nozul aralığı 20 mm, basınç 1 atm ve uyarlanabilir hesaplama ağı 600-900 nokta olarak belirlenmiştir. Kimyasal tepkime mekanizması olarak GRI-Mech 3.0 kullanılmıştır.

Durum I’de maksimum OH mol kesrinin %25’ine karşılık gelen sınırın konumu, deneydeki alev genişliğinin karşılığı olarak tanımlanmıştır. Durum II’de aynı %25 ölçütü hidrojen giriş sınırı ile reaksiyon bölgesi arasındaki pilot kalkma yüksekliğini belirlemek için kullanılmıştır.

Bir boyutlu modelde hidrojen oranı arttıkça hem alev genişliği hem pilot kalkma yüksekliği düzenli biçimde artmıştır. Momentum göstergesiyle yapılan doğrusal uyumda Durum I ve Durum II için R² = 0,95 elde edilmiştir. Bu güçlü uyum, yerel reaksiyon tabakalarının ana eğilimlerinin indirgenmiş bir karşı akış modeliyle açıklanabileceğini göstermiştir.

Ancak simülasyondaki mutlak alev genişliği yaklaşık 2,1-3,9 mm, pilot kalkma yüksekliği yaklaşık 2,0-6,6 mm düzeyindedir ve deneysel makroskopik değerlerle doğrudan örtüşmemektedir. Araştırmacılar da modelin değerini mutlak boyut tahmininde değil, hidrojen oranı ve momentum değiştikçe alev şeklinin hangi yönde değiştiğini açıklamasında görmektedir.

Çalışmada kullanılan temel denklemler

Ana girdap sayısı

\[ S_M = \frac{\pi(r_o^2-r_P^2)(r_i-t/2)}{4\,t\,h}\,r_o \]

S_M ana girdap oluşturucunun boyutsuz girdap yoğunluğudur. r_o, r_P ve r_i sırasıyla ana nozul çıkış yarıçapını, pilot nozul yarıçapını ve ana geçidin iç yarıçapını; t teğetsel kanal genişliğini, h kanat yüksekliğini göstermektedir. Denklem dizgisinde son r_o teriminin kesrin dışında görünmesi, bu terimin çarpan mı yoksa paydanın devamı mı olduğu konusunda tipografik belirsizlik oluşturmaktadır. Çalışma bu parantezlemeyi ayrıca açıklamamıştır; ifade burada metindeki sembol sırasına sadık kalınarak aktarılmıştır.

Pilot girdap sayısı

\[ S_P = \frac{2}{3}\left[\frac{1-(d_h/d)^3}{1-(d_h/d)^2}\right]\tan\theta \]

S_P pilot girdap sayısıdır. d_h merkez göbek çapını, d pilot girdap oluşturucunun çapını ve θ pilot kanat açısını göstermektedir. Çap oranları ve tanjant terimi boyutsuz olduğundan S_P de boyutsuzdur.

Küresel eşdeğerlik oranı

\[ \Phi_g = s\frac{\dot{m}_{CH_4}+\dot{m}_{H_2}}{\dot{m}_a} \]

Φ_g küresel eşdeğerlik oranıdır. ṁ_CH₄ ve ṁ_H₂ yakıtların, ṁ_a ise havanın kütlesel debisidir. Debiler aynı birimle, örneğin kg/s ile kullanıldığında oran boyutsuzdur. s, karışık yakıtın stokiyometrik yakıt-hava oranını temsil etmektedir.

Stokiyometrik oran

\[ s = \frac{(2-1.5X_{H_2})(W_{O_2}+3.76W_{N_2})}{X_{CH_4}W_{CH_4}+X_{H_2}W_{H_2}} \]

W_i ilgili türün moleküler ağırlığını, X_CH₄ ve X_H₂ yakıt bileşenlerinin mol veya hacim temelli oranlarını göstermektedir. 3,76 katsayısı havadaki azot-oksijen mol oranını temsil eden model katsayısıdır.

Hidrojen içeriği

\[ X_{H_2} = \frac{\dot{V}_{H_2}}{\dot{V}_{CH_4}+\dot{V}_{H_2}} \]

V̇_H₂ ve V̇_CH₄ hidrojen ile metanın hacimsel debileridir. Aynı hacimsel debi birimi kullanıldığında X_H₂ boyutsuzdur ve çalışmada yüzde olarak ifade edilmiştir.

Pilot-ana momentum oranı

\[ J = \frac{\rho_PU_P^2}{\rho_MU_M^2} \]

Bu denklem pilot akımının momentum akısını ana akımın momentum akısına göre karşılaştırmaktadır. J büyüdükçe merkez pilot akımının ana akıma göre aerodinamik etkisi artmaktadır. Çalışmanın temel sonucu, alev şeklinin yalnızca yakıt bileşiminden çok bu göreli momentum dengesiyle düzenli biçimde açıklanabildiğidir.

Grafikler ve şekiller ne göstermektedir?

• Deney düzeneği şeması: Dış CH₄/hava geçidi, iç H₂ pilot geçidi, iki farklı girdap oluşturucu ve kuvars pencereli yanma odasının geometrik ilişkisini göstermektedir.
• OH-PLIF ve PIV şemaları: 283 nm uyarma ışığıyla OH radikali görüntüleme sistemini ve 532 nm lazerle parçacık hareketinden hız alanı çıkaran PIV sistemini göstermektedir.
• Tutunmuş-kalkmış karşılaştırması: Dış reaksiyon tabakasının enjektör yakınında devam ettiği yapı ile kökte sönerek aşağı akışta yeniden oluştuğu yapı arasındaki farkı göstermektedir.
• PIV hız alanları: Her iki rejimde iç devridaim bölgesi bulunduğunu, ancak kalkmış alevde ters akışın enjektöre daha fazla yaklaştığını ve ana kesme tabakasındaki hız gradyanının arttığını göstermektedir.
• Anlık ve ortalama OH görüntüleri: Kalkmış alevin daha geniş ortalama OH bölgesinin anlık alev kırışması ve konum değişkenliğinden kaynaklandığını göstermektedir.
• Kararlılık haritası: Hidrojen artışının tutunmuş bölgeyi, ana hız artışının ise kalkmış bölgeyi genişlettiğini göstermektedir.
• OTSU grafiği: Farklı koşullardaki otomatik eşiklerin maksimum OH sinyalinin yaklaşık %20-%30’u arasında toplandığını göstermektedir.
• Şekil parametresi grafikleri: Alev genişliği ve pilot kalkma yüksekliğinin hidrojenle birlikte genel olarak arttığını göstermektedir.
• Momentum korelasyonları: İç reaksiyon tabakası konumunun momentum oranına, kalkmış alev genişliğine kıyasla daha düzenli tepki verdiğini göstermektedir.
• OPPDIF sonuçları: Bir boyutlu modelin deneysel mutlak boyutları değil, değişim yönlerini ve momentumla olan ilişkileri yeniden üretebildiğini göstermektedir.

Hidrojen ve eşdeğerlik oranına göre verilen OH-PLIF görüntülerinin her biri kendi maksimum OH yoğunluğuna göre ayrı ayrı normalize edilmiştir. Bu nedenle görüntüler farklı koşullar arasındaki mutlak OH miktarını veya toplam reaksiyon şiddetini karşılaştırmak için kullanılamaz. Görüntüler yalnızca reaksiyon bölgesinin konumu, genişliği ve şekli hakkında karşılaştırma sunmaktadır.

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

• Hidrojen oranının artırılması, incelenen çift girdaplı enjektörde pilot desteğini güçlendirerek tutunmuş alev rejimini genişletmiştir.
• Ana akış hızının 15 m/s’ye çıkarılması, kalkmış alev koşullarını belirgin biçimde artırmıştır.
• Kalkmış alev, dış alev kökünün yerel sönmesi ve aşağı akışta yeniden kararlı hâle gelmesiyle açıklanmıştır.
• Alev genişliği ve pilot kalkma yüksekliği, hidrojen içeriğiyle birlikte genel olarak artmıştır.
• Pilot kalkma yüksekliği hem tutunmuş hem kalkmış rejimde momentum oranıyla güçlü ilişki göstermiştir.
• Bir boyutlu karşı akış modeli, üç boyutlu türbülanslı alevin temel şekil değişimi eğilimlerini nitel olarak yeniden üretmiştir.

Çalışmanın kanıtlamadığı sonuçlar

• Çalışma ticari bir gaz türbininde saha performansı veya uzun süreli işletme başarısı göstermemektedir.
• NOx, CO, yanmamış hidrokarbon veya toplam CO₂ emisyonları ölçülmemiştir.
• Hidrojen oranının sistem verimini, elektrik üretim maliyetini veya türbin ömrünü ne ölçüde değiştirdiği belirlenmemiştir.
• Flashback güvenliği deneysel olarak haritalanmamış ve malzeme dayanımı değerlendirilmemiştir.
• OH-PLIF görüntüleri mutlak ısı açığa çıkışı ölçümü değildir.
• Kalkmış alevde daha güçlü kırışma görülmesi, türbülans şiddetinin veya türbülanslı yanma hızının kesin olarak arttığını kanıtlamamaktadır.
• Bir boyutlu simülasyon bütün girdaplı akış alanını, üç boyutlu türbülansı veya zamana bağlı akustik etkileşimleri çözmemektedir.
• Korelasyonlar nedenselliği tek başına kanıtlamaz ve momentum oranı özellikle kalkmış alev genişliğindeki bütün değişkenliği açıklamamaktadır.

Türkiye açısından ne ifade etmektedir?

Araştırma, doğal gazla çalışan gaz türbinlerinin hidrojen katkılı yakıtlara uyarlanmasında yalnızca yakıt yüzdesine odaklanmanın yeterli olmadığını göstermektedir. Pilot ve ana geçitlerin geometrisi, hızları, yoğunlukları ve göreli momentumu birlikte tasarlanmadığında daha reaktif bir yakıt kullanılması alevin istenen noktada tutulmasını tek başına garanti etmemektedir.

Türkiye’de gaz türbini, yanma odası, enerji dönüşümü ve hidrojen teknolojileri üzerinde çalışan üniversiteler ile sanayi kuruluşları açısından araştırmanın en aktarılabilir yönü, yüksek maliyetli üç boyutlu hesaplamalara geçmeden önce OH-PLIF/PIV ölçümlerinin OTSU tabanlı şekil göstergeleriyle sayısallaştırılması ve yerel reaksiyon tabakalarının indirgenmiş karşı akış modelleriyle incelenebilmesidir. Bununla birlikte mevcut sonuçlar tek bir laboratuvar enjektörü ve sınırlı çalışma aralığına aittir. Türkiye’deki gerçek türbinlere aktarım için basınçlı yanma deneyleri, emisyon ölçümleri, flashback sınırları, termoakustik kararlılık, malzeme dayanımı ve uzun süreli sistem testleri gereklidir.

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

Deney koşullarının teknik özeti

OH-PLIF ölçümü

• 30 W darbeli Nd:YAG lazer, 10 kHz ve 532 nm çıkış kullanılmıştır.
• Boyar lazer çıkışı 283 nm’ye dönüştürülerek OH A-X (1,0) geçişinin Q₁ çizgisi uyarılmıştır.
• Darbe enerjisi yaklaşık 75 µJ/darbe, darbe süresi 10 ns’dir.
• OH sinyali 315 ± 10 nm bant geçiren filtreyle kaydedilmiştir.
• Her koşul için 0,4 saniyelik, 10 kHz örnekleme hızında toplam 4.000 görüntü alınmıştır.
• Lazer tabakası yoğunluk profili kararlı Bunsen alevi kullanılarak belirlenmiş ve kalibrasyon bir saatlik aralıklarla tekrarlanmıştır.
• OH-PLIF uzamsal çözünürlüğü yaklaşık 0,11 mm/piksel olarak bildirilmiştir.
• Görüntüler MATLAB tabanlı kodla işlenmiş ve 4.000 kare ortalanmıştır.

PIV hız ölçümü

• 20 W gücünde sürekli 532 nm MGL-W-532 lazer kullanılmıştır.
• Akış izleyicisi olarak katı alümina parçacıkları kullanılmıştır.
• Görüntüler NAC MEMRECAM HX-7s kamerayla 5 kHz hızında alınmıştır.
• Ardışık görüntüler arasındaki süre Δt = 200 µs’dir.
• Hız vektörleri PIVlab üzerinden çapraz korelasyon yöntemiyle hesaplanmıştır.
• Sorgulama penceresi 16 × 16 piksel, bindirme oranı %50 ve vektör aralığı yaklaşık 1,17 mm’dir.
• Alt piksel yer değiştirme hatası yaklaşık 0,1 piksel kabul edildiğinde hız ölçümünün göreli belirsizliği ortalama hızın yaklaşık %1-%3’ü olarak tahmin edilmiştir.
• Toplam 4.000 görüntüden anlık ve ortalama hız alanları elde edilmiştir.

Sayısal model koşulları

Ana bulgular

1. Çift girdaplı enjektörde tutunmuş ve kalkmış olmak üzere iki temel alev rejimi belirlenmiştir.
2. Hidrojen içeriği arttıkça tutunmuş alev bölgesi genişlemiş, ana hız arttıkça kalkmış alev bölgesi genişlemiştir.
3. Metan-onlyaklaşımında Φ_g ≤ 0,5 civarında fakir karışım sönmesi görülürken hidrojen eklenmesi daha fakir koşullarda alevin sürmesini sağlamıştır.
4. Kalkmış alevde dış reaksiyon tabakasının kökü enjektör yakınında yerel olarak sönmüş ve aşağı akışta yeniden kararlı hâle gelmiştir.
5. Kalkmış alevde ana kesme tabakasındaki hız gradyanı ve iç devridaim bölgesinin enjektöre doğru nüfuzu artmıştır.
6. Alev genişliği ile pilot alev kalkma yüksekliği hidrojen içeriğiyle birlikte genel olarak yükselmiştir.
7. Deneysel pilot kalkma yüksekliğinin momentum göstergesiyle uyumu tutunmuş alevde R² = 0,90, kalkmış alevde R² = 0,92 olmuştur.
8. Tutunmuş alev genişliği için R² = 0,92, kalkmış alev genişliği için daha düşük olan R² = 0,67 elde edilmiştir.
9. Bir boyutlu modelde alev genişliği ve pilot kalkma yüksekliği için R² = 0,95 düzeyinde momentum ilişkisi elde edilmiştir.
10. Bir boyutlu model deneysel mutlak boyutları vermemiş, ancak hidrojen ve momentum değişimlerine bağlı eğilimleri doğru yönde yeniden üretmiştir.

İstatistiksel ve ölçümsel değerlendirme

Çalışmada doğrusal ilişkiler R² değerleriyle sunulmuştur. Buna karşılık regresyon katsayılarının güven aralıkları, p değerleri, bağımsız deney tekrar sayıları veya şekil parametreleri için standart sapma ve hata çubukları verilmemiştir. Dört bin kare, alevin zamansal davranışını örneklemektedir; bunlar bağımsız deneysel tekrarlarla aynı şey değildir. PIV hız belirsizliği %1-%3 olarak tahmin edilmiş, fakat alev genişliği ve kalkma yüksekliği için ayrı bir belirsizlik yayılımı sunulmamıştır.

Yöntemde görülen iç tutarsızlıklar

• PIV görüntü ölçeği metinde 0,146 nm/piksel olarak belirtilmiştir. Aynı paragrafta 16 pikselin 2,34 mm’ye karşılık geldiği yazılmıştır. 16 × 0,146 mm hesabı yaklaşık 2,34 mm verdiği için “nm” ifadesinin “mm” olması muhtemel görünmektedir; ancak çalışma bunu doğrulamadığından birim sessizce düzeltilmemelidir.
• OH-PLIF kamera modeli yöntem metninde Photron HSS6, Şekil 2’de ise Photron Fastcam SA-X2 olarak adlandırılmıştır.
• Ana girdap sayısı denkleminde son r_o teriminin kesir yapısındaki konumu tipografik olarak açık değildir.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı: Experimental investigation of CH4/H2 co-firing flame structure with comparison to 1D simulation in dual-swirl injector

Yazarlar: Jinseong Kim; Hyunchang Lee; Keeman Lee; Matthew Dunn; Assaad Masri.

Yazar sıralaması: Yukarıdaki liste çalışmadaki özgün sırayı korumaktadır.

Eş birinci yazar: Çalışmada eş birinci yazarlık veya eş katkı beyanı yer almamaktadır.

Sorumlu yazarlar: Keeman Lee ve Assaad Masri.

Kurumlar:

• School of Mechanical and Aerospace Engineering / Center for Aerospace Engineering Research, Sunchon National University, Suncheon 57922, Republic of Korea.
• School of Aerospace, Mechanical and Mechatronic Engineering, The University of Sydney, NSW 2006, Australia.

DOI:10.2139/ssrn.6944015

Dergi: Hakemli bir dergi veya nihai dergi yayını bu sürüm üzerinden doğrulanmamıştır.

Yayın platformu: SSRN.

Özgün yayınevi/platform sağlayıcısı: SSRN bir hakemli dergi yayınevi değildir. Elsevier Inc. tarafından işletilen preprint ve erken aşama araştırma platformudur.

Yayın yılı: 2026.

Kaynak türü: Deneysel araştırma ve bir boyutlu sayısal modelleme içeren preprint araştırma makalesi.

Hakemlik durumu: Bu çalışma bir preprinttir ve hakem değerlendirmesinden geçmemiştir. Bulgular bu sınırlılık dikkate alınarak okunmalıdır.

Resmî kaynak bağlantısı:https://ssrn.com/abstract=6944015

Finansman: Çalışmada Korea Institute of Energy Technology Evaluation and Planning tarafından desteklenen RS-2023-00254668 numaralı 150 MW F-sınıfı gaz türbini hidrojen birlikte yakma dönüşüm projesi ile KEIT ve Kore Cumhuriyeti Ticaret, Sanayi ve Enerji Bakanlığı tarafından desteklenen RS-2025-07982969 numaralı hidrojen karışımlı gaz ve alev tanılama projesi belirtilmiştir. Matthew Dunn ve Assaad Masri için Australian Research Council desteği bildirilmiştir.

Bu Türkçe açıklama, yüklenen çalışmanın metni, denklemleri, tabloları, grafik ve şekilleri incelenerek hazırlanmıştır. Çalışmada bulunmayan bilimsel bulgular ana içeriğe eklenmemiştir. Dış kaynaklar yalnızca başlık, yazarlar, DOI, SSRN kaydı, platform niteliği ve hakemlik durumunun bibliyografik doğrulaması amacıyla kullanılmıştır.

Çalışmanın temel sınırlılıkları; hakem değerlendirmesinden geçmemesi, tek bir laboratuvar enjektör geometrisiyle sınırlı olması, atmosferik basınçta yürütülmesi, bağımsız deney tekrarları ve şekil parametresi belirsizliklerinin ayrıntılı verilmemesi, emisyon ile flashback ölçümlerinin bulunmaması ve bir boyutlu modelin üç boyutlu türbülanslı akışı mutlak değerleriyle temsil etmemesidir. OTSU ve %25 OH ölçütleri karşılaştırma için seçilmiş ampirik sınırlardır; fiziksel alev yüzeyinin kesin konumu olarak değerlendirilmemelidir.