Bu çalışma, ağır hizmet tipi ve karıştırıcı beslemeli bir hidrojen kıvılcım ateşlemeli motorda meydana gelen geri tepmenin yalnızca emme manifoldu basınç sinyali kullanılarak aynı çevrim içinde algılanıp algılanamayacağını incelemektedir. Geri tepme, silindir içinde yanması gereken hidrojen-hava alevinin emme zamanı sırasında ters yönde ilerleyerek emme manifolduna ulaşması ve manifold içinde bulunan yanıcı karışımı tutuşturmasıdır. Karıştırıcı tipi sistemlerde hidrojen ve hava emme supabından önce karıştırıldığı için manifoldun büyük bölümü sürekli yanıcı karışım içermekte; bu nedenle geri tepme, port enjeksiyonlu sistemlere göre daha yüksek basınç yükleri oluşturabilmektedir.

Araştırmacılar, 11 litrelik altı silindirli ağır hizmet tipi bir dizel motoru tek aktif silindirli hidrojen kıvılcım ateşlemeli deney motoruna dönüştürmüştür. Aktif silindir hacmi 1.841 cm³, sıkıştırma oranı 13:1 ve deney devri 1.800 dev/dakikadır. Emme portunun yaklaşık 50 mm önüne yerleştirilen yüksek bant genişlikli piezoelektrik basınç sensörü 200 kHz ile örneklenmiş; geri tepme sinyali zaman alanındaki hızlı basınç yükselişiyle frekans alanındaki geniş bant enerji artışının birlikte değerlendirilmesiyle belirlenmiştir.

İncelenen beş geri tepme olayında emme manifoldu tepe basıncı 251–635 kPa arasında değişmiş, hesaplanan basınç artışı manifold içindeki hidrojen yükünün yaklaşık yüzde 79,6–89,3’ünün yanmasına karşılık gelmiştir. Geri tepme çevrimi, normal çevrimler arasında tek bir aşırı uç değer olarak belirgin biçimde ayrılmıştır. Uyarlanabilir CUSUM/Page–Hinkley algılayıcısı alarmı, aynı çevrim içinde ve manifold basıncı maksimuma ulaşmadan 4–16 milisaniye önce üretmiştir.

Basınç izinden yeniden oluşturulan egzoz supabı açılma anındaki artık gaz sıcaklığı, karışım zenginleştikçe yaklaşık 600 K’den 784–791 K düzeyine yükselmiştir. Bu sıcaklık aralığı, çalışmada kullanılan hidrojen-hava kendiliğinden tutuşma sınırı olan yaklaşık 773–793 K ile örtüşmektedir. Araştırmacılar bu sonucu, sıcak artık gazların veya egzoz supabı ve buji çevresindeki daha sıcak yerel bölgelerin emme manifoldundaki karışımı yeniden tutuşturduğu termal mekanizmanın kanıtı olarak yorumlamaktadır.

Bununla birlikte sonuçlar beş geri tepme olayıyla, tek motor geometrisiyle ve yalnızca 1.800 dev/dakika koşuluyla sınırlıdır. Önerilen algılama sistemi üretim tipi bir motor kontrol ünitesinde uygulanmamış, yanlış alarm oranı bağımsız bir doğrulama veri kümesinde ölçülmemiş ve artık gaz sıcaklığı doğrudan sıcaklık sensörüyle doğrulanmamıştır. Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.

Geri tepme nedir?

Normal bir kıvılcım ateşlemeli motorda hidrojen-hava karışımının yanması silindir içinde, bujinin ateşlemesiyle başlamalı ve yanma odasında tamamlanmalıdır. Geri tepmede ise alev veya tutuşma kaynağı emme supabı üzerinden ters yönde ilerleyerek emme portuna ve manifolda ulaşır. Manifold içinde yanıcı karışım bulunuyorsa yanma silindir dışında devam eder.

Bu olay aşağıdaki sonuçlara yol açabilir:

• Emme manifoldunda ani ve yüksek basınç oluşması,
• Gaz kelebeği, karıştırıcı, sensör ve kanal bileşenlerinde mekanik yük meydana gelmesi,
• Manifold içindeki taze karışımın yer değiştirmesi ve hacimsel verimin düşmesi,
• Sonraki çevrimlerde vuruntu veya süper vuruntu riskinin artması,
• Motorun erişebileceği yük ve güç yoğunluğunun sınırlandırılması,
• Tekrarlanan olaylarda donanım hasarı ve güvenlik riski oluşması.

Hidrojenin düşük tutuşma enerjisi, yüksek alev hızı ve geniş yanabilirlik aralığı bu riski artırmaktadır. Çalışmada hidrojen için minimum tutuşma enerjisi 0,02 mJ, stokiyometrik laminer alev hızı 3,51 m/s ve yanabilirlik aralığı hacimsel olarak yüzde 4,7–75 olarak verilmektedir.

Karıştırıcı tipi sistemler neden daha hassastır?

Hidrojen motorlarında üç temel yakıt besleme düzeni bulunmaktadır:

1. Doğrudan enjeksiyon: Hidrojen doğrudan silindir içine gönderilir.
2. Port yakıt enjeksiyonu: Hidrojen emme portunda ve emme supabına yakın bir noktada püskürtülür.
3. Harici karıştırma veya ön karışımlı emme: Hidrojen, emme supabından daha yukarıda havayla karıştırılır.

Port enjeksiyonunda püskürtme zamanlaması, hidrojenin emme portunda bulunduğu süreyi kısaltabilir. Harici karıştırmalı sistemde ise hidrojen-hava karışımı manifoldun önemli bölümünde sürekli bulunmaktadır. Geriye ilerleyen alevin erişebileceği yanıcı karışım hacmi daha büyük olduğundan basınç darbesi de daha ağır olabilir.

Bu sistemler basit ve görece düşük maliyetli olmaları nedeniyle ağır hizmet tipi doğal gaz motorlarının hidrojene dönüştürülmesinde önemini korumaktadır. Çalışmanın endüstriyel önemi, geri tepme riskinin yüksek olduğu bu mimariyi gerçek ağır hizmet tipi motor boyutlarında incelemesinden kaynaklanmaktadır.

Literatürde hangi boşluk hedeflenmiştir?

Önceki araştırmalarda geri tepme; manifold basıncı, silindir basıncı, motor titreşimi, akustik sinyal, iyon akımı ve optik OH/OH* görüntüleme yöntemleriyle incelenmiştir. Optik yöntemler tutuşmanın fiziksel konumunu gösterebilmekte, fakat şeffaf erişim pencereleri ve laboratuvar donanımı gerektirmektedir. Silindir içi piezoelektrik basınç sensörleri yüksek doğruluk sağlasa da seri üretim ve uzun süreli saha kullanımı için pahalı ve dayanım açısından sınırlı olabilir.

Standart manifold basınç sensörleri üretim motorlarında uygulanabilir olmakla birlikte bant genişlikleri hızlı şok ve akustik bileşenleri ayrıntılı biçimde kaydetmeye yetmeyebilir. Araştırmacılar, 200 kHz örneklenen hızlı manifold basıncıyla geri tepmenin hem zaman hem frekans özelliklerini tek çevrim içinde birleştirmeyi amaçlamıştır.

Çalışmanın üç temel katkı iddiası şöyledir:

• Tek manifold sensörüyle aynı çevrim içinde geri tepme algılamak,
• Basınç artışından olayın şiddetini hacimden bağımsız bir oranla hesaplamak,
• Geri tepme sınırını artık gaz sıcaklığıyla ilişkilendiren termal mekanizmayı göstermek.

Deney motoru nasıl hazırlanmıştır?

Deney platformunun temeli, sıralı altı silindirli, turboşarjlı ve doğrudan enjeksiyonlu 11.046 cm³ ağır hizmet tipi dizel motordur. Altı silindirden yalnızca biri aktif bırakılmıştır. Diğer beş piston üzerinde delikler açılmış ve bu silindirlerin emme-egzoz supapları çıkarılarak taze karışımın devre dışı silindirlere girmesi engellenmiştir.

Özgün dizel enjektör çıkarılarak yerine buji yerleştirilmiştir. Sıkıştırma oranı, piston çanağı yeniden işlenerek 17:1’den 13:1’e düşürülmüştür. Motor soğutma suyu çıkış sıcaklığı harici eşanjörle 80°C’de tutulmuştur.

Deney düzeneği hangi akımları ve sensörleri içeriyordu?

Deney düzeneğini gösteren şemada hidrojen tüpü, basınç regülatörü, kütlesel debi kontrolörü, hava filtresi, elektrikli gaz kelebeği, emme manifoldu, tek silindirli motor, egzoz analiz sistemi, motor kontrolörü ve veri toplama bilgisayarı aynı akış üzerinde gösterilmektedir. Hidrojen ile hava, emme manifoldundan önce birleşmektedir.

Hızlı manifold sensörünün 200 kHz örnekleme hızı 100 kHz Nyquist frekansı sağlamıştır. Bu bant genişliği, geri tepme basınç darbesinin yalnızca yavaş basınç yükselişini değil, akustik ve şok cephesi bileşenlerini de incelemek amacıyla seçilmiştir.

Emme havası nasıl ölçülmüştür?

Yanma sırasında oluşan büyük manifold basınç dalgalanmalarına rağmen hedeflenen hava debisini kararlı tutmak için kritik akış orifislerinden oluşan bir demet kullanılmıştır. Orifisler referans debimetreyle kalibre edilmiş ve hava kütlesel debisi için şu regresyon bağıntısı elde edilmiştir:

\[ \dot{m}_a = 4.444\times10^{-7}\left(\frac{P_0}{T_0^{0.5}}\right)+1.189\times10^{-5} \]

• \(\dot{m}_a\): Hava kütlesel debisi, kg/s,
• \(P_0\): Orifis öncesindeki mutlak basınç, Pa,
• \(T_0\): Orifis öncesindeki mutlak sıcaklık, K’dir.

Hangi çalışma koşulları denenmiştir?

Motor devri bütün deneylerde 1.800 dev/dakikada tutulmuştur. Yük, hava fazlalık katsayısı ve ateşleme zamanlamasını kapsayan dört deney seti kullanılmıştır.

Tablodaki negatif kıvılcım açıları üst ölü noktadan önceki ateşlemeyi, pozitif değerler ise üst ölü noktadan sonraki ateşlemeyi ifade etmektedir. λ değerinin 1’den büyük olması karışımın stokiyometrik karışımdan daha fakir olduğunu gösterir. Çalışmada “zenginleştirme”, λ değerinin örneğin 3,0’dan 2,2’ye düşürülmesidir; karışım yine de stokiyometriye göre fakirdir.

Komut edilen λ değeri nasıl doğrulanmıştır?

Deney düzeneğinde λ, hidrojen kütlesel debi kontrolörü ile kritik orifis hava debisinden hesaplanmıştır. Geniş bant oksijen sensörünün doğrusal hesaplaması hidrokarbon yakıta göre ayarlandığından hidrojen için yanlı olabilmektedir. Araştırmacılar bu nedenle kuru egzoz oksijenini ve egzoz hidrojen kaçağını kullanan bağımsız bir atom dengesi oluşturmuştur:

\[ A=x_{O_2}-0.5x_{H_2} \]

\[ B=1-1.5x_{H_2} \]

\[ \lambda=\frac{0.5(A-B)}{cA-0.5B} \]

\[ c=\frac{K+1}{2}=2.3866,\qquad K=\frac{1-\phi}{\phi},\qquad \phi=0.2095 \]

• \(x_{O_2}\): Kuru egzozdaki oksijen mol kesri,
• \(x_{H_2}\): Kuru egzozdaki yanmamış hidrojen mol kesri,
• \(\phi\): Havadaki oksijen mol kesri,
• \(\lambda\): Gerçek havanın stokiyometrik hava gereksinimine oranıdır.

Temiz veriye sahip Set 2’de egzoz oksijeniyle hesaplanan λ ile komut edilen λ arasındaki rms fark 0,03; Set 3’te 0,10’dur. Hidrokarbon yakıta göre kalibre edilmiş UEGO sensörünün karşılık gelen farkları 0,07 ve 0,66’dır. Çok fakir koşullarda egzoz oksijeninin λ’ya duyarlılığı azaldığından Set 4’te fark 1,09’a kadar yükselmiştir. Araştırmacılar komut edilen λ’yı ana referans olarak kullanmış, oksijen dengesini λ yaklaşık 3,5’in altında bağımsız kontrol olarak değerlendirmiştir.

Motor verimi yükle nasıl değişmiştir?

Yaklaşık λ = 3,05–3,25 aralığında ve kıvılcım zamanlaması her yük için ayarlanırken net indike termik verim düşük yükte yüzde 45,5’ten yükselmiş, 264–372 kPa net IMEP civarında yüzde 48,1–48,2 ile geniş bir maksimum oluşturmuş ve 815 kPa yükte yüzde 44,7’ye gerilemiştir.

Buradaki değer net indike termik verimdir; motordan alınan fren gücünü, bütün yardımcı sistem kayıplarını veya araç verimini temsil etmez. En düşük yükte sabit ısı ve yanma kayıplarının üretilen işe oranı yüksektir. Yük arttıkça bu kayıpların göreli etkisi azalır. Daha yüksek yüklerde ise silindir basıncı, sıcaklık, duvar ısı kaybı ve ateşleme zamanını geciktirme gereksinimi verimi yeniden düşürmektedir.

IMEP değişim katsayısı yüzde 2,4–3,0 aralığında kaldığı için yüksek yükteki verim azalmasının yanma kararsızlığından ziyade termodinamik kayıplarla ilişkili olduğu yorumlanmıştır.

Karışım fakirleştikçe yanma nasıl değişmektedir?

λ yükselip karışım fakirleştikçe hidrojenin alev yayılma hızı düşmektedir. CA50’nin yaklaşık 5 °aTDC’de tutulabilmesi için buji daha erken ateşlenmiştir. Hafif yükte λ 3,20’den 4,40’a yükselirken kıvılcım zamanlaması −25’ten −32 °aTDC’ye ilerletilmiş ve yüzde 10–90 yanma süresi 37,4 dereceden 52,6 krank açısı derecesine çıkmıştır.

Karışım zenginleştirilerek λ yaklaşık 2’ye düşürüldüğünde yanma hızlanmış, basınç yükselme hızı ve NOx artmıştır. Yüksek yükte λ yaklaşık 2,55 civarında NOx 87–158 ppm düzeyine çıkarken λ 3,2 ve üzerindeki fakir koşullarda 1 ppm veya altında kalmıştır.

Ateşleme zamanlamasının iki tarafındaki sınırlar

Sabit yük ve λ altında yapılan kıvılcım taraması, kullanılabilir ateşleme penceresinin iki farklı anormal yanma biçimiyle sınırlandığını göstermiştir:

• Çok erken ateşleme: Silindir basıncı tepesine binen yüksek frekanslı salınımlar ve vuruntu,
• Çok geç ateşleme: Yanmanın emme zamanına kadar uzaması ve geri tepme.

Örnek taramada −22 °aTDC kıvılcım zamanında vuruntu görülürken −2 °aTDC civarında geri tepme meydana gelmiştir. Bu sonuç, geri tepmenin yalnızca karışımın zenginliğiyle değil, önceki çevrim yanmasının ne kadar geç tamamlandığıyla da ilişkili olduğunu göstermektedir.

Geri tepmenin basınç izi nasıldır?

Normal çevrimlerde emme manifoldu basıncı dar bir bantta kalmıştır. Geri tepme çevriminde ise basınç emme zamanı sırasında aniden yükselmiş ve tek çevrimlik büyük bir darbe oluşturmuştur.

Beş ardışık çevrimin karşılaştırıldığı grafikte geri tepme çevriminin manifold basıncı yaklaşık 460 kPa’ya çıkarken aynı çevrimin silindir basıncı olağanüstü yüksek değildir. Bu ayrım, olayın silindir içindeki klasik vuruntu değil, emme tarafında meydana gelen yanma olduğunu göstermektedir.

İki farklı çalışma noktasında 300 çevrim üst üste çizildiğinde tek geri tepme çevrimi normal dağılımın çok dışında kalmıştır:

• Bir olayda yaklaşık 46 kPa taban basıncına karşı 251 kPa tepe,
• Diğerinde yaklaşık 100 kPa taban basıncına karşı 458 kPa tepe görülmüştür.

Beş olayın tamamını içeren şiddet tablosunda en yüksek tepe 635 kPa’ya ulaşmıştır.

Zaman-frekans imzası ne göstermektedir?

Geri tepme yalnızca basınç genliğiyle değil, spektral içeriğiyle de normal çevrimlerden ayrılmıştır. Zaman alanında keskin başlangıçlı ve büyük genlikli bir basınç yükselişi oluşurken güç spektral yoğunluğu geniş bir frekans bandında yükselmiştir.

Spektrogramlarda yüksek enerji, emme zamanına karşılık gelen dar bir krank açısı aralığında yoğunlaşmaktadır. Olay belirli tek bir rezonans frekansından oluşmamakta; geniş bantlı ve zamansal olarak yerelleşmiş bir geçici darbe niteliği göstermektedir.

Çalışmanın farklı bölümlerinde yüksek frekans bandının alt sınırı 50 Hz, 200 Hz ve yaklaşık 500 Hz olarak farklı biçimlerde ifade edilmektedir. Şekillerde üst sınır yaklaşık 26–27 kHz’tir. Bu bant tanımlarındaki farklılık, uygulama öncesinde netleştirilmelidir.

Tek sensörlü algılayıcı nasıl çalışmaktadır?

Önerilen ana algılayıcı, hızlı manifold basıncının 500 Hz kesim frekanslı ikinci dereceden nedensel Butterworth yüksek geçiren filtreden geçirilmesiyle başlamaktadır. Filtrelenmiş örnek karesi alınarak yüksek frekans enerjisi oluşturulur:

\[ e_n=\left(x_n^{HP}\right)^2 \]

Bu enerjinin çalışan ortalaması ve yayılımı üstel ağırlıklı hareketli ortalamayla izlenir. Makalede standartlaştırılmış enerji şu biçimde yazılmıştır:

\[ z_n=\frac{e_n-\mu_n}{v_n} \]

• \(e_n\): Karesi alınmış yüksek geçiren manifold sinyali,
• \(\mu_n\): Enerjinin çalışan ortalaması,
• \(v_n\): Makalede çalışan varyans olarak tanımlanan terim,
• \(z_n\): Standartlaştırılmış enerji olarak adlandırılan değer.

Ardından Page–Hinkley biçimindeki CUSUM biriktiricisi kullanılır:

\[ g_n=\max\left(0,\;g_{n-1}+z_n-k\right) \]

\[ g_n>h\quad \Rightarrow\quad alarm \]

Çalışmada boyutsuz kayma parametresi \(k=1\), alarm eşiği ise \(h=10\) olarak seçilmiştir. Normal sinyalde biriktirici sıfıra yakın kalırken geri tepmenin yükselen kenarında hızla artarak eşiği geçmektedir.

Boyutsuz zaman alanı özellikleri

Basınç genliği ve yükselme hızı, çalışma noktasına bağlı mutlak eşikler yerine yerel standart sapmaya göre normalize edilmiştir:

\[ \Pi_{swing}=\frac{P_{max}-P_{min}}{\sigma_P} \]

\[ \Pi_{grad}=\frac{\max\left|dP/dt\right|}{\sigma_P} \]

Geri tepme çevrimlerinde normalize basınç salınımı ve basınç yükselme gradyanı yerel normal ortalamanın yaklaşık 23 standart sapma üzerinde bulunmuştur.

Spektral düzlüğün hesaplanması

Tek çevrimin spektral düzlüğü, spektral değerlerin geometrik ortalamasının aritmetik ortalamasına oranıdır:

\[ SFM=\frac{\exp\left(\frac{1}{N}\sum_k\ln S_k\right)}{\frac{1}{N}\sum_k S_k} \]

Çalışmada ayrıca dürtüsel ve durağan olmayan sinyallere duyarlı spektral basıklık kullanılmıştır. Bir çevrimin geri tepme olarak işaretlenmesi için zaman alanındaki basınç/rise-rate testinin ve frekans alanındaki testin aynı çevrimde anlaşması amaçlanmıştır.

Bu ortak koşulun gerekçesi şöyledir:

• Yavaş yük değişimleri büyük basınç değişimi oluştursa bile yüksek yükselme hızı ve geniş bant enerji testini geçmeyebilir.
• Tek örnekli elektriksel veya sensör kaynaklı darbeler zaman testini geçse bile beklenen fiziksel spektral yapıyı taşımayabilir.
• İki bağımsız özelliğin aynı anda kullanılması yanlış alarm olasılığını azaltabilir.

Geri tepme şiddeti nasıl hesaplanmıştır?

Araştırmacılar olayın şiddetini, ölçülen manifold basınç artışının aynı manifold yükünün sabit hacimde, adyabatik ve tamamen yanması hâlinde oluşacak kuramsal basınç artışına oranı olarak tanımlamıştır:

\[ \eta=\frac{\Delta P_{meas}}{\Delta P_{pred}} \]

• \(\eta\): Yanma tamamlanma oranı,
• \(\Delta P_{meas}\): Ölçülen manifold basınç yükselişi,
• \(\Delta P_{pred}\): Tam yanma için tahmin edilen basınç yükselişi.

Tam yanma basınç tavanı şu denklemle hesaplanmıştır:

\[ \Delta P_{pred}=P_0\left(\frac{N_bT_{ad}}{N_uT_0}-1\right) \]

• \(P_0\): Başlangıç manifold basıncı,
• \(N_u\): Yanmamış karışımın mol sayısı,
• \(N_b\): Yanmış ürünlerin mol sayısı,
• \(T_0\): Başlangıç karışım sıcaklığı,
• \(T_{ad}\): Adyabatik alev sıcaklığıdır.

Mol sayıları ve adyabatik alev sıcaklığı NASA-7 termodinamik polinomlarıyla hesaplanmıştır. Basınçların oranı kullanıldığı için manifold hacmi denklemden düşmektedir. Bu özellik, gerçek manifold hacminin kesin bilinmediği durumlarda karşılaştırmalı bir şiddet göstergesi sağlamaktadır.

Ancak \(\eta\), doğrudan kimyasal tür ölçümünden elde edilen gerçek yanmış hidrojen fraksiyonu değildir. Sabit hacim, adyabatik yanma, tekdüze başlangıç durumu ve basınç dalgasının ölçüm noktasındaki temsili gibi varsayımlara dayanan model tabanlı bir basınç vekilidir.

Beş geri tepme olayının şiddeti

Beş olayda manifold yükünün büyük bölümü yanmıştır. Bununla birlikte tamamlanma oranı λ ile tekdüze biçimde değişmemektedir. Bunlar kontrollü bir şiddet taramasından ziyade farklı koşullarda oluşan beş bağımsız olaydır.

Artık gaz sıcaklığı nasıl hesaplanmıştır?

Geri tepmenin termal öncülünü araştırmak için silindir basınç izinden iki sıcaklık yeniden oluşturulmuştur. Silindirde hapsolmuş kütle, emme supabının kapanma durumundan hesaplanmıştır:

\[ mR=\frac{P_{IVC}V_{IVC}}{T_{IVC}} \]

Tepe basınç anındaki yanmış gaz sıcaklığı:

\[ T_b=T_{IVC}\frac{P_{max}V\left(\theta_{P_{max}}\right)}{P_{IVC}V_{IVC}} \]

Egzoz supabının açılma anındaki ortalama artık gaz sıcaklığı:

\[ T_{EVO}=T_{IVC}\frac{P_{EVO}V_{EVO}}{P_{IVC}V_{IVC}} \]

• \(T_b\): Tepe basınç anındaki iki bölgeli yanmış gaz sıcaklığı,
• \(T_{EVO}\): Egzoz supabı açılırken silindirdeki kütle ortalamalı artık gaz sıcaklığı,
• \(T_{IVC}\): Emme supabı kapanırken karışım sıcaklığı,
• \(P\) ve \(V\): İlgili krank açısındaki silindir basıncı ve hacmidir.

Basit ideal gaz hesabı, iki bölgeli ayırma ve bileşim düzeltilmiş NASA-7 hesabıyla çapraz kontrol edilmiştir. NASA-7 sonuçları yaklaşık yüzde 10 daha düşük olmakla birlikte λ ile değişim eğilimi aynıdır.

Isıl mekanizma için bulunan sonuçlar

Karışım zenginleştikçe hem yanmış gaz hem artık gaz sıcaklığı düzenli biçimde yükselmiştir. Geri tepme olaylarında kütle ortalamalı \(T_{EVO}\) değeri 784–791 K’ye ulaşmıştır. Çalışmada verilen hidrojen-hava minimum kendiliğinden tutuşma sıcaklığı yaklaşık 773–793 K’dir.

Egzoz supabı veya buji yakınındaki yerel sıcak noktaların kütle ortalamasından 100–200 K daha sıcak olabileceği; bunların yaklaşık 950 K olarak verilen yerel sıcak nokta sınırına ulaşabileceği belirtilmiştir. Buna göre önerilen mekanizma şöyledir:

1. λ azaldıkça karışım zenginleşir ve yanma sıcaklığı yükselir.
2. Yanma sonunda silindirde kalan artık gaz daha sıcak hâle gelir.
3. Egzoz ve emme zamanları arasındaki gaz değişimi sırasında sıcak artık gaz veya yerel sıcak nokta emme tarafına ulaşır.
4. Manifoldda hazır bulunan hidrojen-hava karışımı kendiliğinden tutuşur.
5. Alev manifold boyunca ters yönde ilerler ve büyük basınç darbesi oluşturur.

Bu sonuç mekanizmayı kuvvetle desteklemektedir; ancak sıcak nokta doğrudan görüntülenmemiş veya hızlı termometreyle ölçülmemiştir. Geri tepme çevrimlerinde yavaş veri toplama kaydı bulunmadığı için hapsolmuş kütle en yakın kararlı çalışma noktasından λ ile ölçeklenmiştir. Bu nedenle sıcaklıklar doğrudan ölçüm değil, basınç temelli yeniden oluşturmadır.

Algılayıcının performansı

Tek parametre seti beş geri tepme çalışmasına uygulanmıştır. Her çalışmada 300 çevrim bulunmaktadır. Geri tepme çevrimi hem CUSUM hem de spektral basıklıkta en yüksek skoru üretmiştir.

CUSUM eşiği geri tepme çevriminin emme zamanında, −273 ile −214 °aTDC arasında aşılmıştır. Alarm manifold basıncı tepeye ulaşmadan 4–16 ms, yani 1.800 dev/dakikada yaklaşık 44–172 krank açısı derecesi önce üretilmiştir.

Bu süre kuramsal olarak aşağıdaki müdahalelere olanak sağlayabilir:

• Hidrojen akışının kesilmesi,
• Kıvılcımın kesilmesi veya yeniden zamanlanması,
• Gaz kelebeği ya da emniyet valfinin yönetilmesi,
• Manifold basınç tahliye yolunun açılması,
• Motor yükünün düşürülmesi.

Ancak bu müdahaleler deneysel olarak uygulanmamış ve geri tepme hasarını gerçekten azalttıkları gösterilmemiştir. Çalışma algılama için kullanılabilecek zaman penceresini belirlemektedir.

Grafiklerin birlikte anlattığı fiziksel süreç

• Deney düzeneği şeması: Hidrojen ve havanın emme manifoldundan önce birleştiğini, hızlı manifold basıncı ile silindir basıncının aynı krank açısı temelinde kaydedildiğini gösterir.
• λ doğrulama grafiği: Egzoz oksijeniyle hesaplanan λ değerlerinin temiz deney setlerinde komut edilen değerlere yakın olduğunu gösterir.
• Verim grafiği: Yaklaşık λ = 3’te net indike verimin düşük ve çok yüksek yükler arasında part-load bölgesinde maksimum oluşturduğunu gösterir.
• Yanma tanımlayıcıları: Fakirleşmeyle yanma süresinin uzadığını ve kıvılcımın daha erkene alınması gerektiğini gösterir.
• Beş çevrimlik basınç grafiği: Geri tepmenin yalnızca bir çevrimde manifold tarafında büyük basınç artışı oluşturduğunu gösterir.
• 300 çevrim topluluğu: Geri tepme çevriminin yükten bağımsız biçimde olağan dağılımın çok dışında kaldığını gösterir.
• Ateşleme taraması: Erken ateşlemede vuruntu, geç ateşlemede geri tepme oluşarak güvenli zamanlama penceresini iki taraftan sınırlar.
• PSD ve spektrogramlar: Geri tepmenin dar bantlı rezonans yerine geniş bantlı ve emme zamanında yerelleşmiş bir darbe olduğunu gösterir.
• Sıcaklık grafiği: λ azaldıkça \(T_b\) ve \(T_{EVO}\) değerlerinin yükseldiğini, geri tepme noktalarının hidrojen kendiliğinden tutuşma bandına ulaştığını gösterir.
• CUSUM grafikleri: Beş olayın tamamında geri tepme çevriminin alarm eşiğini çok büyük farkla geçtiğini ve alarmın basınç tepesinden önce oluştuğunu gösterir.

Çalışmanın güçlü yönleri

• Geri tepmeyi yapay olarak oluşturulan küçük bir düzenek yerine gerçek ağır hizmet tipi motor geometrisinde incelemektedir.
• Karıştırıcı tipi hidrojen beslemesindeki büyük manifold yükünü doğrudan değerlendirmektedir.
• 200 kHz manifold basıncı sayesinde zaman ve frekans bilgisini birlikte kullanmaktadır.
• Algılama, şiddet ölçümü ve mekanizma incelemesini aynı deney veri kümesinde birleştirmektedir.
• Mutlak basınç eşiği yerine normalize edilmiş ve uyarlanabilir özellikler kullanmayı amaçlamaktadır.
• Algılayıcının yalnızca çevrimi değil, çevrim içindeki alarm zamanını da raporlamaktadır.
• Yanma tamamlanması oranı manifold hacmine cebirsel olarak bağlı olmayan bir şiddet göstergesi sunmaktadır.
• Artık gaz sıcaklığı ile geri tepme sınırı arasında iki ayrı yük taramasında aynı eğilimi göstermektedir.
• Algılayıcının beş ayrı geri tepme olayında aynı parametrelerle güçlü ayrım sağladığını göstermektedir.
• Yakıt türüne bağımlı UEGO ölçümünü egzoz oksijen atom dengesiyle çapraz kontrol etmektedir.

Çalışmanın başlıca sınırlılıkları

• Yalnızca beş geri tepme olayı ayrıntılı olarak değerlendirilmiştir.
• Bütün deneyler 1.800 dev/dakikada yürütülmüştür; farklı devirlerde doğrulama bulunmamaktadır.
• Motor, altı silindirli üretim motorunun yalnızca bir aktif silindirli laboratuvar sürümüdür.
• Karıştırıcı tipi tek manifold geometrisi, tek sıkıştırma oranı ve tek buji konumu kullanılmıştır.
• Algılayıcı üretim tipi motor kontrol ünitesinde gerçek zamanlı olarak uygulanmamıştır.
• Hesaplama süresi ve işlemci yükü sayısal olarak raporlanmamıştır.
• STFT pencere uzunluğu, pencere tipi ve örtüşme oranı belirtilmemiştir.
• Algılayıcı için duyarlılık, özgüllük, kesinlik, yanlış alarm oranı veya ROC eğrisi verilmemiştir.
• Eşikler bağımsız test veri kümesinde doğrulanmamıştır.
• “Kalibrasyonsuz” ifadesine rağmen eşik, çalışma zarfındaki en kötü normal çevrim skoruna göre seçilmektedir.
• Artık gaz sıcaklığı doğrudan ölçülmemiştir.
• Geri tepme çevrimlerindeki hapsolmuş kütle, en yakın kararlı koşuldan ölçeklenmiştir.
• Yanma tamamlanması sabit hacim ve adyabatik tam yanma varsayımlarına dayanmaktadır.
• Basınç tabanlı şiddet göstergesi gerçek parça hasarı, gerilme veya yorulma ölçümleriyle ilişkilendirilmemiştir.
• Ham yüksek hızlı veriler ve analiz kodu paylaşılmamış; yalnızca talep üzerine sağlanacağı belirtilmiştir.
• Çalışma hakem değerlendirmesinden geçmemiştir.

Çalışma neyi desteklemektedir?

• Karıştırıcı tipi ağır hizmet hidrojen motorunda geri tepmenin emme manifoldu basıncında tek çevrimlik büyük bir darbe oluşturduğunu desteklemektedir.
• Geri tepmenin geniş bantlı ve zaman içinde yerelleşmiş spektral imzaya sahip olduğunu göstermektedir.
• Uyarlanabilir yüksek frekans enerjisi temelli CUSUM yönteminin incelenen beş olayda geri tepmeyi basınç tepesinden önce algılayabildiğini göstermektedir.
• Manifold basınç yükselişinin olay şiddeti için karşılaştırmalı bir vekil olarak kullanılabileceğini desteklemektedir.
• Karışım zenginleştikçe artık gaz sıcaklığının yükselerek hidrojen kendiliğinden tutuşma bölgesine ulaştığını göstermektedir.
• Çok erken kıvılcımın vuruntu, çok geç kıvılcımın ise geri tepme sınırı oluşturduğunu göstermektedir.
• Yaklaşık λ = 3 ve orta yük bölgesinin incelenen motorda yüksek net indike verimle geri tepmesiz çalışma sağlayabildiğini göstermektedir.

Çalışma neyi kanıtlamamaktadır?

• Algılayıcının bütün hidrojen motorlarında aynı eşikle çalışacağını kanıtlamamaktadır.
• Eşiğin farklı motor devirlerinden gerçekten bağımsız olduğunu göstermemektedir.
• Standart seri üretim MAP sensörünün 200 kHz laboratuvar sensörüyle aynı algılama başarısını sağlayacağını göstermemektedir.
• Geri tepmenin bütün durumlarda yalnızca artık gaz tarafından başlatıldığını kanıtlamamaktadır.
• Yağ damlacığı, karbon birikintisi, buji elektrodu veya sıcak metal kenar gibi diğer yerel tutuşma kaynaklarını dışlamamaktadır.
• Alarm sonrasında yakıt veya kıvılcım kesmenin olayın mekanik etkisini önleyeceğini deneysel olarak göstermemektedir.
• Yanma tamamlanması oranının gerçek mekanik hasarla doğrusal ilişkili olduğunu kanıtlamamaktadır.
• Yüzde 48 civarındaki net indike verimin fren termik verimi veya araç verimi olduğunu göstermemektedir.
• Hidrojen motorlarının genel olarak sıfır yaşam döngüsü karbon emisyonuna sahip olduğunu göstermemektedir.
• Hakemli bilimsel uzlaşı veya ticari sistem doğrulaması sunmamaktadır.

Endüstriyel ve teknolojik anlamı

Çalışmanın uygulama açısından en önemli sonucu, manifold basıncının yalnızca geri tepme olduktan sonra kayıt yapan bir sensör değil, olay büyürken müdahale sinyali sağlayabilecek bir kontrol girdisi olabileceğidir. Basınç tepesinden önceki 4–16 ms kısa bir süre olsa da yüksek hızlı motor kontrolü için anlamlı olabilir.

Yöntemin üretim sistemine taşınabilmesi için daha düşük maliyetli sensörlerde bant genişliği, sensör yerleşimi, titreşim bağışıklığı, elektromanyetik gürültü, sürekli çalışma dayanımı ve ECU hesaplama yükü değerlendirilmelidir. Farklı motor devirleri, çok silindirli ateşleme sıraları, turboşarj basıncı, EGR, farklı manifold hacimleri ve port enjeksiyonlu sistemler üzerinde bağımsız doğrulama yapılmalıdır.

Gelecekte algılama sinyali kapalı çevrim kontrolle birleştirilebilir. Sistem alarm sonrası hidrojen akışını azaltabilir, ateşlemeyi geçici olarak durdurabilir, λ hedefini yükseltebilir veya yük talebini sınırlandırabilir. Ancak bu müdahalelerin zamanlaması ve güvenliği ayrı deneylerle doğrulanmalıdır.

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

Teknik yöntem özeti

Beş olayın algılama sonuçları

Temel nicel bulgular

Algılama ve mekanizma zinciri

1. Hidrojen ve hava manifold içinde önceden karışmış durumdadır.
2. Karışım zenginleştikçe silindir yanma ve artık gaz sıcaklıkları yükselir.
3. Geç yanma veya sıcak yerel gaz, emme olayı sırasında manifolda ulaşır.
4. Manifold karışımı tutuşarak geniş bantlı ve hızlı bir basınç darbesi üretir.
5. Yüksek geçiren manifold enerjisi CUSUM eşiğini aşar.
6. Zaman alanındaki hızlı yükseliş ile frekans alanındaki geniş bant enerji aynı çevrimde doğrulanır.
7. Alarm basınç tepesinden önce üretilir.
8. Basınç artışının tam yanma tavanına oranı olay şiddeti için hesaplanır.

Uygulamadan önce gerekli doğrulamalar

• Farklı motor devirlerinde eşik aktarılabilirliği,
• Çok silindirli motorda silindir kaynağının belirlenmesi,
• Üretim tipi düşük maliyetli basınç sensörleriyle bant genişliği yeterliliği,
• Bağımsız veri kümesinde yanlış pozitif ve yanlış negatif oranları,
• Motor titreşimi, supap darbeleri ve elektriksel parazit altında dayanıklılık,
• Ardışık veya tekrarlanan geri tepme olaylarında uyarlanabilir tabanın davranışı,
• Yakıt ve kıvılcım kesme komutunun gerçek koruma etkisi,
• Farklı manifold hacimleri ve turboşarj basınçlarında şiddet ölçeği,
• Doğrudan sıcaklık veya optik ölçümle termal mekanizmanın doğrulanması.

Kaynak ve Yöntem Notu

Çalışmanın özgün adı: Backfire in a heavy-duty mixer-type hydrogen spark-ignition engine: detection, severity quantification, and the residual-gas thermal mechanism

Yazarlar ve sıraları: Seongsu Kim; Hyunwook Park; Seungmook Oh; Changup Kim; Yonggyu Lee; Junsun Lee; Junghwan Kim.

Eş katkı bilgisi: PDF’de eş katkı veya eş birinci yazarlık beyanı bulunmamaktadır.

Sorumlu yazar: Junghwan Kim.

Sorumlu yazar e-postası: jkim77@cau.ac.kr

Kurumlar:

1. Department of Energy Systems Engineering, Chung-Ang University, 84 Heukseok-ro, Seoul, Republic of Korea 06974.
2. Department of Engine Research, Korea Institute of Machinery and Materials, 156 Gajungbuk-ro, Daejeon, Republic of Korea 34103.

Yazar-kurum eşleştirmesi: Seongsu Kim ve Junghwan Kim, Chung-Ang University; Hyunwook Park, Seungmook Oh, Changup Kim, Yonggyu Lee ve Junsun Lee, Korea Institute of Machinery and Materials.

DOI: 10.2139/ssrn.6945147

Kaynak türü: Ağır hizmet tipi hidrojen motorunda yürütülen deneysel araştırma ve sinyal işleme preprinti.

Hakemlik durumu: Bu çalışma bir preprinttir ve hakem değerlendirmesinden geçmemiştir.

Yayın platformu: SSRN.

SSRN yüklenme tarihi: 15 Haziran 2026.

Sayfa sayısı: 37.

Dergi: Hakemli bir dergi adı veya kabul edilmiş nihai yayın bilgisi doğrulanamamıştır.

Yayınevi durumu: Çalışma bir dergi makalesi olarak yayımlanmamıştır. SSRN, Elsevier bünyesinde işletilen bir erken araştırma ve preprint paylaşım platformudur; bu durum çalışmanın hakemli bir Elsevier dergi makalesi olduğu anlamına gelmez.

Resmî bağlantılar:SSRN kayıt sayfası ve DOI bağlantısı.

Finansman: Çalışma, Kore Ticaret, Sanayi ve Enerji Bakanlığı tarafından finanse edilen Technology Innovation Program 20018148 kapsamında desteklenmiştir.

Çıkar çatışması: Yazarlar, çalışmayı etkileyebilecek bilinen finansal çıkar veya kişisel ilişki bulunmadığını beyan etmiştir.

Veri erişimi: Verilerin talep üzerine sağlanacağı belirtilmiştir. Ana PDF ile birlikte ham 200 kHz basınç verisi, analiz kodu veya gerçek zamanlı kontrol yazılımı paylaşılmamıştır.

Yazar katkıları: Seongsu Kim deney, biçimsel analiz ve ilk taslak yazımını; Hyunwook Park deney, veri düzenleme ve görselleştirmeyi; Seungmook Oh finansman ve proje yönetimini; Changup Kim kavramsallaştırmayı; Yonggyu Lee yöntemi; Junsun Lee deney ve kaynakları; Junghwan Kim kavramsallaştırma, danışmanlık ve metin incelemesini üstlenmiştir.

Bu Türkçe makale, yüklenen 37 sayfalık PDF’nin metni, denklemleri, deney şemaları, basınç grafikleri, spektral yoğunluk grafikleri, spektrogramları, CUSUM sonuçları ve bütün tabloları incelenerek hazırlanmıştır. PDF’de bulunmayan yeni deneysel sonuç, güvenlik garantisi, ticari uygulanabilirlik veya motor dayanıklılığı iddiası eklenmemiştir. Dış kaynaklar yalnızca DOI, SSRN yüklenme tarihi, iletişim yazarı ve platform kimliğini doğrulamak için kullanılmıştır.

Metinde belirlenen iç tutarsızlıklar ve yeniden üretilebilirlik sorunları

• Hızdan bağımsız eşik iddiası: Özet ve girişte algılama eşiğinin yük ve motor devrinden bağımsız olduğu belirtilmektedir. Ancak bütün deneyler 1.800 dev/dakikada yürütülmüştür. Farklı devirlerde deneysel doğrulama yapılmadan hızdan bağımsızlık gösterilmiş sayılamaz.
• “Kalibrasyonsuz” ifadesi: Algılayıcının çalışma noktasına göre ayrı kalibrasyon tablosu gerektirmediği doğrudur. Bununla birlikte boyutsuz eşiklerin çalışma zarfındaki en kötü normal çevrim skoruna göre seçildiği belirtilmektedir. Bu nedenle yöntem tamamen kalibrasyonsuz değil, çalışma noktasına özgü olmayan ortak eşikli bir yaklaşımdır.
• Tek sensörle neden teşhisi iddiası: Geri tepme algılama ve şiddet hesabında manifold basıncı kullanılmaktadır. Ancak artık gaz sıcaklığı ve neden analizi silindir içi basınç izine dayanmaktadır. Bu nedenle manifold sensörünün tek başına hem algılama, hem şiddet, hem de neden teşhisi sağladığı ifadesi yöntemle tam uyumlu değildir.
• Z-skoru denklemindeki varyans: Nomenklatürde \(v_n\) çalışan varyans olarak tanımlanmış ve Denklem 6’da enerji farkı \(v_n\)’ye bölünmüştür. Standart bir z-skorunda farkın varyansa değil standart sapmaya, yani \(\sqrt{v_n}\)’ye bölünmesi beklenir. \(v_n\) gerçekte standart sapmayı temsil ediyorsa adlandırma; varyansı temsil ediyorsa denklem boyutsal açıdan düzeltilmelidir.
• Spektral düzlük yorumu: Denklemde verilen geometrik ortalama/aritmetik ortalama oranı, genel olarak düz ve geniş bantlı spektrumda 1’e; tonal veya keskin tepeli spektrumda 0’a yaklaşır. Metinde bunun tersi ifade edilmektedir. Geri tepme verilerinde kullanılan dB değeri ve bant seçimi bu genel tanımdan farklı bir davranış oluşturuyorsa açıklanmalıdır.
• Eksik STFT parametreleri: Pencere uzunluğu ve örtüşme oranı için “parametre optimizasyonu tamamlandığında belirtilecek” ifadesi bulunmaktadır. Bu açık bir taslak yer tutucusudur ve spektrogramların bağımsız yeniden üretimini engeller.
• Frekans bandı tutarsızlığı: Metnin farklı bölümlerinde geniş bant analizi için yaklaşık 50 Hz–26 kHz, 200 Hz–26/27 kHz ve 0,5–26 kHz aralıkları kullanılmaktadır. Yüksek geçiren filtre ise 500 Hz’dir. Nihai algılama bandı tek ve açık biçimde tanımlanmamıştır.
• Tepe basıncı özeti: Sonuç bölümünde manifold basıncının yaklaşık 250–460 kPa’ya ulaştığı yazmaktadır. Beş olayın yer aldığı Tablo 7’de ise 547 ve 635 kPa tepe değerleri de bulunmaktadır. Sonuç cümlesi yalnızca iki örnek grafiği yansıtmakta, bütün olay aralığını yansıtmamaktadır.
• Tablo 7’de bulunmayan hidrojen kütlesi: Metin, Tablo 7’de 1,5 litrelik manifold varsayımıyla hidrojen kütlelerinin de raporlandığını belirtmektedir. Görünen tabloda böyle bir sütun bulunmamaktadır.
• Gerçek zamanlı uygulama eksikliği: Karar akış şeması ve çevrim başına hesaplama maliyeti sonraki revizyona bırakılmıştır. Bu nedenle ECU üzerinde gerçek zamanlı uygulanabilirlik doğrudan gösterilmemiştir.
• Sınıflandırma performansı: Beş olayda güçlü ayrım gösterilse de bağımsız test kümesi, yanlış alarm sayısı, duyarlılık, özgüllük ve güven aralığı raporlanmamıştır.
• Termal nedenin dolaylı oluşu: \(T_{EVO}\) sıcaklığı doğrudan ölçüm değil, ideal gaz ve basınç izi temelli kütle ortalamasıdır. Geri tepme çalışmaları için kütle en yakın kararlı koşuldan ölçeklenmiştir. Sonuç artık gaz mekanizmasını destekler, ancak yerel tutuşma kaynağını doğrudan gözlemlemez.
• Şiddet ölçütünün fiziksel varsayımları: Manifold yanması gerçekte sabit hacimli ve tam adyabatik değildir; supap, kanal, ısı kaybı ve basınç dalgası etkileri vardır. Bu nedenle \(\eta\), doğrudan yanmış hidrojen fraksiyonu yerine model tabanlı şiddet vekili olarak yorumlanmalıdır.
• Düzenleme hatası: Algılayıcı yönteminin sonunda cümlenin ortasında başlayan ve önceki açıklamayı tekrar eden bozuk bir metin parçası bulunmaktadır. Bu durum preprintin nihai editoryal kontrolden geçmediğini göstermektedir.

Bu sorunlar, manifold basıncı sinyalinin geri tepme algılamasında güçlü bir bilgi kaynağı olduğunu gösteren deneysel bulguları ortadan kaldırmamaktadır. Ancak farklı hızlarda, farklı motorlarda, daha fazla olayla ve bağımsız doğrulama verileriyle çalışılmadan yöntemin evrensel veya üretime hazır olduğu sonucuna varılmamalıdır.