
Ushbu tadqiqot og‘ir yuk toifasidagi, aralashtirgichli (mikserli) ta’minotga ega vodorodli uchqundan alangalanish dvigatelida yuzaga keladigan teskari chaqnashni (backfire) faqat kiritish (so‘rish) kollektori bosim signali yordamida ayni bir sikl ichida aniqlash mumkin yoki yo‘qligini o‘rganadi. Teskari chaqnash — silindr ichida yonishi kerak bo‘lgan vodorod-havo alangasining kiritish takti vaqtida orqaga qarab harakatlanib, kiritish klapani orqali kollektorga o‘tishi va u yerdagi yonuvchi aralashmani yondirib yuborishidir. Mikserli tizimlarda vodorod va havo kiritish klapanidan oldinroq aralashtirilgani sababli kollektorning katta qismi doimiy yonuvchi aralashma bilan to‘la bo‘ladi; shu bois teskari chaqnash taqsimlangan purkashli (PFI) tizimlarga nisbatan ancha yuqori bosim zarbalarini keltirib chiqarishi mumkin.
Tadqiqotchilar 11 litrli olti silindrli og‘ir yuk dizel dvigatelini bitta faol silindrli vodorodli uchqundan alangalanish tajriba dvigateliga aylantirdilar. Faol silindr hajmi 1 841 cm³, siqish darajasi 13:1 va sinov aylanishlar soni 1 800 ayl/daqiqadir. Kiritish kanalining taxminan 50 mm oldiga o‘rnatilgan yuqori chastotali piezoelektrik bosim datchigi 200 kHz chastotada signallarni qayd etdi; teskari chaqnash signali vaqt sohasidagi tezkor bosim ortishi hamda chastota sohasidagi keng polosali energiya o‘sishini birgalikda baholash orqali aniqlandi.
O‘rganilgan beshta teskari chaqnash hodisasida kiritish kollektori cho‘qqi bosimi 251–635 kPa oralig‘ida bo‘ldi, hisoblangan bosim ortishi kollektor ichidagi vodorod zaryadining taxminan 79,6–89,3% qismi yonganiga to‘g‘ri keldi. Teskari chaqnash yuz bergan sikl odatiy sikllar orasida yaqqol ekstremal chetlanish sifatida ajralib turdi. Moslashuvchan CUSUM/Page–Hinkley detektori ogohlantirish signalini ayni o‘sha sikl ichida va kollektor bosimi maksimumga yetishidan 4–16 millisoniya oldin hosil qildi.
Silindr bosimi egri chizig‘idan qayta tiklangan chiqarish klapani ochilish paytidagi (EVO) qoldiq gazlar harorati aralashma boyigan sari taxminan 600 K dan 784–791 K darajasiga ko‘tarildi. Bu harorat oralig‘i vodorod-havo aralashmasining o‘z-o‘zidan alangalanish chegarasi bo‘lgan 773–793 K bilan to‘liq mos tushadi. Tadqiqotchilar ushbu natijani issiq qoldiq gazlar yoki chiqarish klapani va sham atrofidagi qizigan mahalliy zonalar kiritish kollektoridagi aralashmani qayta yondirib yuboradigan termik mexanizmning dalili deb baholaydilar.
Shu bilan birga, natijalar beshta chaqnash hodisasi, bitta dvigatel geometriyasi va faqat 1 800 ayl/daq rejimi bilan cheklangan. Tavsiya etilgan aniqlash tizimi seriyali dvigatel boshqaruv blokida (ECU) sinab ko‘rilmagan, soxta ogohlantirishlar ehtimoli mustaqil ma’lumotlar bazasida tekshirilmagan va qoldiq gaz harorati to‘g‘ridan-to‘g‘ri harorat datchigi bilan o‘lchanmagan. Tadqiqot taqrizdan o‘tmagan preprint hisoblanadi.
Teskari chaqnash (backfire) nima?
Odatiy uchqundan alangalanish dvigatelida vodorod-havo aralashmasining yonishi silindr ichida, sham uchquni bilan boshlanishi va yonish kamerasida tugashi kerak. Teskari chaqnashda esa alanga yoki o‘t oldirish manbai kiritish klapani orqali teskari yo‘nalishda harakatlanib, kiritish kanali va kollektoriga o‘tadi. Agar kollektorda yonuvchi aralashma bo‘lsa, yonish silindrdan tashqarida davom etadi.
Bu hodisa quyidagi oqibatlarga olib kelishi mumkin:
- Kiritish kollektorida to‘satdan kuchli bosim zarbasi hosil bo‘lishi,
- Drossel qopqog‘i, mikser, datchiklar va quvurlarda mexanik zo‘riqish paydo bo‘lishi,
- Kollektordagi toza aralashmaning siqib chiqarilishi va hajmiy FIK ning pasayishi,
- Keyingi sikllarda detonatsiya yoki super-detonatsiya xavfining ortishi,
- Dvigatel erisha oladigan yuklama va quvvat zichligining cheklanishi,
- Takroriy hodisalarda uskunaning sinishi va xavfsizlik buzilishi.
Vodorodning juda past alangalanish energiyasi (0,02 mJ), yuqori laminar alanga tezligi (3,51 m/s) va keng yonish oralig‘i (hajm bo‘yicha 4,7–75%) bu xavfni keskin oshiradi.
Mikserli (aralashtirgichli) tizimlar nima sababdan sezgirroq?
Vodorod dvigatellarida uchta asosiy yoqilg‘i uzatish sxemasi mavjud:
- To‘g‘ridan-to‘g‘ri purkash (DI): Vodorod to‘g‘ridan-to‘g‘ri silindr ichiga purkaladi.
- Portga purkash (PFI): Vodorod kiritish kanalida, klapanga yaqin joyda purkaladi.
- Tashqi aralashtirish (mikserli): Vodorod kiritish klapanidan ancha yuqorida havo bilan oldindan aralashtiriladi.
PFI tizimlarida purkash vaqtini boshqarish orqali vodorodning kanalda bo‘lish vaqtini qisqartirish mumkin. Mikserli tizimda esa vodorod-havo aralashmasi kollektorning butun hajmida doimiy mavjud bo‘ladi. Orqaga qaytgan alanga yetib borishi mumkin bo‘lgan hajm ancha katta bo‘lgani uchun bosim to‘lqini ham kuchliroq bo‘ladi.
Bunday tizimlar sodda va arzon bo‘lgani sababli og‘ir yuk tabiiy gaz dvigatellarini vodorodga moslashtirishda dolzarbligini saqlab qolmoqda. Tadqiqotning amaliy ahamiyati aynan shu xavfli arxitekturani og‘ir yuk dvigateli miqyosida o‘rganganidadir.
Adabiyotda qanday bo‘shliq maqsad qilingan?
Oldingi ishlarda teskari chaqnash kollektor bosimi, silindr bosimi, tebranish, akustik signal, ion toki va optik OH/OH* tasvirlash usullari bilan o‘rganilgan. Optik usullar alangani ko‘rsatsa-da, shaffof oynalar va laboratoriya sharoitini talab qiladi. Silindr ichidagi piezoelektrik datchiklar juda aniq, ammo seriyali ishlab chiqarish uchun qimmat va resursi cheklangan bo‘lishi mumkin.
Standart kollektor bosimi (MAP) datchiklari arzon bo‘lsa-da, ularning chastota kengligi tezkor zarba va akustik to‘lqinlarni to‘liq qayd etishga yetmaydi. Tadqiqotchilar 200 kHz chastotada o‘lchangan tezkor kollektor bosimi orqali chaqnashning vaqt va chastota xususiyatlarini yagona sikl ichida birlashtirishni maqsad qilganlar.
Tadqiqotning uch asosiy hissasi:
- Bitta kollektor datchigi bilan ayni sikl ichida chaqnashni aniqlash,
- Bosim ortishidan foydalanib hajmga bog‘liq bo‘lmagan nisbiy shiddat ko‘rsatkichini hisoblash,
- Teskari chaqnash chegarasini qoldiq gazlar harorati bilan bog‘lovchi termik mexanizmni ko‘rsatib berish.
Tajriba dvigateli qanday tayyorlandi?
Tajriba stendi 6 silindrli, turbopurkovchili, 11 046 cm³ hajmli og‘ir yuk dizel dvigateliga asoslangan. 6 ta silindrdan faqat bittasi faol qoldirilgan. Qolgan 5 ta porshen teshilgan, ularning klapanlari chiqarib tashlanib, taze aralashma kirmasligi ta’minlangan.
| Dvigatel parametri | Qiymati |
|---|---|
| Bazaviy dvigatel | Qatordagi 6 silindrli, turbopurkovchili og‘ir yuk dizeli |
| Faol silindr | 6 tadan 1 tasi |
| Faol silindr hajmi | 1 841 cm³ |
| Dvigatelning umumiy hajmi | 11 046 cm³ |
| Silindr diametri | 123 mm |
| Porshen yo‘li (strok) | 155 mm |
| Shatun uzunligi | 275 mm |
| Siqish darajasi | 13:1 (asl dvigatelda 17:1) |
| Klapanlar soni | Silindr boshiga 1 kiritish va 1 chiqarish klapani |
| Yonish turi | Dizel siqilishidan uchqunli alangalanishga o‘zgartirilgan |
| Yonish kamerasi | Qayta yo‘nilgan vannasimon porshen chuqurchasi |
| Sham joylashuvi | Silindr markazidan 3,8 mm siljigan |
| Yoqilg‘i | Kamida 99,9% soflikdagi gazsimon vodorod |
| Yoqilg‘i uzatish | Kiritish klapanidan oldin havo bilan aralashtirish |
Dizel forsunkasi chiqarilib, o‘rniga o‘t oldirish shami o‘rnatildi. Siqish darajasi porshen chuqurchasini qayta ishlash orqali 17:1 dan 13:1 ga tushirildi. Sovitish suvi harorati 80°C da doimiy saqlandi.
Stend qanday datchiklar va tizimlarni o‘z ichiga olgan?
Tajriba sxemasida vodorod balloni, regulyator, massaviy sarf nazoratchisi, havo filtri, elektron drossel, kiritish kollektori, bir silindrli motor, gaz analizatori, boshqaruv bloki va ma’lumot yig‘ish kompyuteri aks ettirilgan.
| O‘lchov | Uskuna yoki usul | Xususiyati |
|---|---|---|
| Silindr bosimi | Kistler 6043ACP piezoelektrik datchigi | Silindr kallagiga o‘rnatilgan |
| O‘rtacha kollektor bosimi (MAP) | Kistler 4045A5 piezorezistiv datchigi | Past chastotali MAP o‘lchovi |
| Tezkor kollektor bosimi | Kistler 6052C piezoelektrik datchigi | Kiritish kanalidan 50 mm oldinda, 200 kHz |
| Yuqori tezlikdagi ma’lumot yig‘ish | AVL Indicom | Krank burchagi bilan sinxron bosim o‘lchovi |
| Chiqindi gazdagi H₂ va O₂ | ECOM MK6000 | H₂: 0–20 000 ppm; O₂: %0–21 |
| NOx | AVL i60 xemiluminessent analizatori | 0–10 000 ppm |
| HC | Olov-ionizatsion detektor (FID) | 0–10 000 ppm |
| CO va CO₂ | NDIR analizatorlari | CO: 0–10 000 ppm; CO₂: %0–25 |
| Chiqindi gazdagi O₂ | Paramagnit kislorod detektori | Yoqilg‘idan mustaqil λ tekshiruvi |
Tezkor datchikning 200 kHz namunalash chastotasi 100 kHz Naykvist chegarasini berib, chaqnash to‘lqinining nafaqat sekin ko‘tarilishini, balki yuqori chastotali zarba va akustik tebranishlarini ham yozib olish imkonini yaratdi.
Kiruvchi havo sarfi qanday o‘lchandi?
Kuchli bosim tebranishlarida ham havo sarfini barqaror ushlab turish uchun kritik oqimli drossellar (orifislar) to‘plami qo‘llanildi. O‘lchovlar etalon sarf o‘lchagich bilan kalibrlanib, quyidagi empirik formula olindi:
\[ \dot{m}_a = 4.444\times10^{-7}\left(\frac{P_0}{T_0^{0.5}}\right)+1.189\times10^{-5} \]
- \(\dot{m}_a\): Havo massaviy sarfi, kg/s,
- \(P_0\): Drosseldan oldingi mutlaq bosim, Pa,
- \(T_0\): Drosseldan oldingi mutlaq harorat, K.
Sinovdan o‘tkazilgan ish rejimlari
Barcha tajribalarda motor tezligi 1 800 ayl/daqiqada saqlandi. To‘rtta tajriba to‘plami o‘tkazildi:
| To‘plam | Tezlik | Yalpi IMEP | λ oralig‘i | O‘t oldirish burchagi | Maqsadi |
|---|---|---|---|---|---|
| Set 1 | 1 800 ayl/daq | 210–830 kPa | 2,05–4,40 | −32 dan −14 °aTDC gacha | Yuklama tahlili, FIK va λ tasdig‘i |
| Set 2 | 1 800 ayl/daq | 400 kPa | 1,80–2,05 | −22 dan +10 °aTDC gacha | Detonatsiya va chaqnash orasidagi o‘t oldirish oynasi |
| Set 3 | 1 800 ayl/daq | 500–820 kPa | 1,69–3,00 | −24 dan −1 °aTDC gacha | Boyitishda chaqnash chegarasi va harorat prekursori |
| Set 4 | 1 800 ayl/daq | 210–410 kPa | 2,05–4,40 | −32 dan −15 °aTDC gacha | CA50 taxminan 5 °aTDC da turganda λ ta’siri |
Manfiy burchaklar YuQNu dan oldingi (bTDC), musbat burchaklar esa YuQNu dan keyingi (aTDC) o‘t oldirishni bildiradi. λ > 1 bo‘lishi aralashmaning stexiometrik holatdan kambag‘alligini anglatadi. «Boyitish» iborasi λ ning 3,0 dan 2,2 ga tushirilishini anglatadi (aralashma baribir kambag‘al ligicha qoladi).
Buyurilgan λ ko‘rsatkichi qanday tekshirildi?
λ qiymati vodorod sarfi nazoratchisi va havo drosseli ko‘rsatkichlaridan hisoblangan. Standart keng polosali kislorod datchigi (UEGO) uglevodorodlarga sozlanganligi tufayli vodorodda xato ko‘rsatishi mumkin. Shu sababli tadqiqotchilar quruq chiqindi gaz tarkibiga asoslangan atomar balans formulasini tuzdilar:
\[ 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}\): Quruq gazdagi O₂ molyar ulushi,
- \(x_{H_2}\): Chiqindi gazdagi yonmagan H₂ molyar ulushi,
- \(\phi\): Havodagi kislorod ulushi (0,2095),
- \(\lambda\): Havo ortiqlik koeffitsiyenti.
2-to‘plamda kislorod balansi bilan hisoblangan λ bilan buyurilgan λ orasidagi rms farq atigi 0,03 bo‘ldi. UEGO datchigida esa bu farq 0,07 va 0,66 ga yetgan. Juda kambag‘al aralashmalarda kislorod sezgirligi pasaygani uchun 4-to‘plamda farq 1,09 gacha chiqqan. Tadqiqotchilar buyurilgan λ ni asosiy qiymat deb olganlar.
Dvigatel FIK ko‘rsatkichi yuklama bilan qanday o‘zgardi?
λ = 3,05–3,25 oralig‘ida va optimal o‘t oldirish burchaklarida sof indikatorli termik FIK past yuklamada 45,5% dan oshib, 264–372 kPa IMEP atrofida 48,1–48,2% lik maksimal darajaga yetdi va 815 kPa yuklamada 44,7% ga tushdi.
Bu ko‘rsatkich sof indikatorli termik FIK bo‘lib, tormoz FIK yoki avtomobil tejamkorligi emas. Past yuklamada issiqlik yo‘qotishlarining nisbiy ulushi katta bo‘ladi. Yuklama oshgani sari bu yo‘qotishlarning nisbiy ulushi kamayadi. Yuqori yuklamalarda esa silindr bosimi, harorat, devorga issiqlik o‘tishi va o‘t oldirishni kechiktirish majburiyati FIK ni yana pasaytiradi.
IMEP variatsiya koeffitsiyenti 2,4–3,0% oralig‘ida bo‘lgani sababli yuqori yuklamadagi pasayish yonish beqarorligidan emas, termodinamik yo‘qotishlardan ekani aytilgan.
Aralashma kambag‘allashganda yonish qanday kechadi?
λ oshib, aralashma suyulgan sari vodorod alangasining tarqalish tezligi pasayadi. CA50 nuqtasini taxminan 5 °aTDC da saqlash uchun o‘t oldirishni ancha oldinga surish kerak bo‘ldi. Kichik yuklamada λ 3,20 dan 4,40 ga chiqqanda o‘t oldirish −25 dan −32 °aTDC ga surildi va 10–90% yonish davomiyligi 37,4° dan 52,6° krank burchagiga yetdi.
| λ | Cho‘qqi silindr bosimi | Maksimal bosim ortish tezligi | 10–90% yonish davomiyligi | CA50 | O‘t oldirish burchagi | NOx |
|---|---|---|---|---|---|---|
| 3,20 | 2,20 MPa | 5,83 MPa/°CA | 37,4° | 5,43 °aTDC | −25 °aTDC | 1 ppm |
| 3,45 | 2,29 MPa | 6,08 MPa/°CA | 42,4° | 5,04 °aTDC | −27 °aTDC | 0 ppm |
| 4,04 | 2,48 MPa | 6,51 MPa/°CA | 48,6° | 5,19 °aTDC | −30 °aTDC | 0 ppm |
| 4,40 | 2,59 MPa | 6,79 MPa/°CA | 52,6° | 5,10 °aTDC | −32 °aTDC | 1 ppm |
Aralashma λ ≈ 2 gacha boyitilganda yonish tezlashib, bosim ortish sur’ati va NOx keskin oshdi. Yuqori yuklamada λ ≈ 2,55 bo‘lganda NOx 87–158 ppm ga yetgan bo‘lsa, λ ≥ 3,2 bo‘lgan kambag‘al rejimlarda 1 ppm yoki undan past bo‘ldi.
O‘t oldirish burchagining ikki tomonlama chegarasi
O‘t oldirish burchagini o‘zgartirish tajribasi xavfsiz ishlash oynasi ikki xil g‘ayritabiiy yonish bilan cheklanganini ko‘rsatdi:
- Haddan tashqari erta o‘t oldirish: Silindr bosimida yuqori chastotali tebranishlar va detonatsiya (knock),
- Haddan tashqari kech o‘t oldirish: Yonishning kiritish taktiga cho‘zilib ketishi va teskari chaqnash (backfire).
Masalan, −22 °aTDC da detonatsiya boshlangan bo‘lsa, −2 °aTDC da teskari chaqnash yuz bergan. Bu holat chaqnash faqat aralashma tarkibiga emas, balki oldingi sikldagi yonish qanchalik kech tugaganiga ham bog‘liqligini isbotlaydi.
Teskari chaqnashning bosim profili qanday bo‘ladi?
Odatiy sikllarda kiritish kollektori bosimi juda tor diapazonda tebranadi. Chaqnash siklida esa kiritish takti vaqtida bosim to‘satdan keskin ko‘tarilib, kuchli yakka zarba to‘lqinini hosil qiladi.
Ketma-ket 5 ta sikl taqqoslanganda chaqnash siklida kollektor bosimi 460 kPa ga yetgani holda silindr ichidagi bosim anomal yuqori bo‘lmagan. Bu farq hodisaning silindr ichidagi detonatsiya emas, balki kiritish tizimidagi tashqi yonish ekanini tasdiqlaydi.
300 ta sikl ustma-ust chizilganda yagona chaqnash sikli normal taqsimotdan butunlay ajralib chiqqan:
- Bitta holatda 46 kPa bazaviy bosimga nisbatan 251 kPa cho‘qqi,
- Boshqasida 100 kPa bazaviy bosimga nisbatan 458 kPa cho‘qqi qayd etilgan.
Barcha 5 ta hodisada maksimal cho‘qqi bosimi 635 kPa ga yetdi.
Vaqt-chastota imzosi nimani ko‘rsatadi?
Teskari chaqnash nafaqat bosim amplitudasi, balki spektral tarkibi bilan ham ajralib turadi. Vaqt sohasida tik o‘sish profili paydo bo‘lsa, spektral quvvat zichligi juda keng chastota polosasida keskin ko‘tariladi.
Spektrogrammalarda yuqori energiya kiritish taktiga to‘g‘ri keluvchi tor krank burchagi oralig‘ida to‘plangan. Hodisa bitta rezonans chastotasidan iborat bo‘lmay, keng polosali o‘tkinchi zarba xususiyatiga ega.
Matnning turli joylarida yuqori chastotali polosa chegarasi 50 Hz, 200 Hz va 500 Hz deb turlicha ko‘rsatilgan. Rasmlarda yuqori chegara taxminan 26–27 kHz ni tashkil etadi.
Bitta datchikli detektor qanday ishlaydi?
Tavsiya etilgan asosiy algoritm tezkor kollektor bosimini 500 Hz chastotali 2-tartibli Batervort (Butterworth) yuqori o‘tkazuvchi filtridan o‘tkazish bilan boshlanadi. Filtrlangan signal kvadratga ko‘tarilib, yuqori chastota energiyasi olinadi:
\[ e_n=\left(x_n^{HP}\right)^2 \]
Ushbu energiyaning joriy o‘rtachasi va tarqoqligi eksponensial harakatlanuvchi o‘rtacha (EWMA) bilan kuzatiladi. Standartlashtirilgan energiya quyidagicha yoziladi:
\[ z_n=\frac{e_n-\mu_n}{v_n} \]
- \(e_n\): Kvadratga ko‘tarilgan yuqori chastotali signal,
- \(\mu_n\): Energiyaning joriy o‘rtachasi,
- \(v_n\): Matnda dispersiya deb atalgan joriy ko‘rsatkich,
- \(z_n\): Standartlashtirilgan energiya.
So‘ngra Page–Hinkley shaklidagi CUSUM yig‘uvchisi ishlatiladi:
\[ g_n=\max\left(0,\;g_{n-1}+z_n-k\right) \]
\[ g_n>h\quad \Rightarrow\quad alarm \]
Tadqiqotda \(k=1\) va ogohlantirish chegarasi \(h=10\) deb olingan. Normal holatda yig‘uvchi nolga yaqin turadi, chaqnash boshlanishi bilan keskin o‘sib chegaradan oshadi.
O‘lchamsiz vaqt parametrlari
Bosim amplitudasi va ortish tezligi mahalliy standart chetlanishga normallashtirilgan:
\[ \Pi_{swing}=\frac{P_{max}-P_{min}}{\sigma_P} \]
\[ \Pi_{grad}=\frac{\max\left|dP/dt\right|}{\sigma_P} \]
Chaqnash sikllarida ushbu parametrlar normal o‘rtacha qiymatdan taxminan 23 standart chetlanishga yuqori bo‘lgan.
Spektral tekislikni hisoblash
Spektral tekislik o‘lchovi (SFM) spektral qiymatlarning geometrik o‘rtachasining arifmetik o‘rtachasiga nisbatidir:
\[ SFM=\frac{\exp\left(\frac{1}{N}\sum_k\ln S_k\right)}{\frac{1}{N}\sum_k S_k} \]
Chaqnashni yakuniy tasdiqlash uchun vaqt sohasidagi mezon bilan chastota sohasidagi spektral basiklik (kurtosis) mezonining ayni siklda mos kelishi shart qilingan. Bu soxta ogohlantirishlarning oldini oladi.
Teskari chaqnash shiddati qanday baholandi?
Tadqiqotchilar chaqnash shiddatini o‘lchangan bosim ortishining ayni kollektor zaryadining doimiy hajmda, adiabatik va to‘liq yonishi natijasida hosil bo‘lishi mumkin bo‘lgan nazariy bosim ortishiga nisbati sifatida aniqladilar:
\[ \eta=\frac{\Delta P_{meas}}{\Delta P_{pred}} \]
- \(\eta\): Yonish to‘liqligi ulushi,
- \(\Delta P_{meas}\): O‘lchangan bosim ortishi,
- \(\Delta P_{pred}\): To‘liq yonish uchun bashorat qilingan bosim ortishi.
Nazariy maksimal bosim ortishi:
\[ \Delta P_{pred}=P_0\left(\frac{N_bT_{ad}}{N_uT_0}-1\right) \]
- \(P_0\): Boshlang‘ich kollektor bosimi,
- \(N_u\): Yonmagan aralashmaning mol soni,
- \(N_b\): Yonish mahsulotlarining mol soni,
- \(T_0\): Boshlang‘ich harorat,
- \(T_{ad}\): Adiabatik alanga harorati.
Termodinamik hisoblar NASA-7 ko‘phadlari orqali bajarilgan. Bosimlar nisbati olingani sababli kollektor hajmi formuladan qisqarib ketadi va hajmga bog‘liq bo‘lmagan taqqoslama ko‘rsatkich beradi.
5 ta hodisaning shiddat ko‘rsatkichlari
| Tajriba | λ | Maksimal kollektor bosimi | O‘lchangan bosim ortishi | Nazariy bosim ortishi | Adiabatik alanga harorati | Yonish to‘liqligi (η) |
|---|---|---|---|---|---|---|
| Set 1, 017 | 2,30 | 458 kPa | 361 kPa | 454 kPa | 1 929 K | %79,6 |
| Set 2, 011 | 2,05 | 251 kPa | 205 kPa | 230 kPa | 2 066 K | %89,3 |
| Set 3, 004 | 2,275 | 455 kPa | 366 kPa | 426 kPa | 1 942 K | %85,9 |
| Set 3, 008 | 2,20 | 547 kPa | 442 kPa | 509 kPa | 1 981 K | %86,7 |
| Set 3, 015 | 2,03 | 635 kPa | 513 kPa | 616 kPa | 2 078 K | %83,3 |
Barcha hodisalarda kollektor ichidagi gazning 79,6–89,3% qismi yonib ketgan.
Qoldiq gazlar harorati qanday hisoblandi?
Silindr bosimi izidan ikkita harorat qayta tiklandi. Kiritish klapani yopilish paytidagi (IVC) massasi:
\[ mR=\frac{P_{IVC}V_{IVC}}{T_{IVC}} \]
Maksimal bosim paytidagi yongan gaz harorati:
\[ T_b=T_{IVC}\frac{P_{max}V\left(\theta_{P_{max}}\right)}{P_{IVC}V_{IVC}} \]
Chiqarish klapani ochilish paytidagi (EVO) qoldiq gaz harorati:
\[ T_{EVO}=T_{IVC}\frac{P_{EVO}V_{EVO}}{P_{IVC}V_{IVC}} \]
Termik mexanizm natijalari
| Yalpi IMEP | λ | Oddiy yongan gaz harorati | NASA-7 harorati | EVO qoldiq gaz harorati | Chaqnash yuz berdimi? |
|---|---|---|---|---|---|
| 510 kPa | 3,00 | 1 730 K | 1 563 K | 601 K | Yo‘q |
| 525 kPa | 2,60 | 1 905 K | 1 715 K | 704 K | Yo‘q |
| 520 kPa | 2,50 | 1 959 K | 1 764 K | 723 K | Yo‘q |
| 500 kPa | 2,275 | 2 070 K | 1 858 K | 791 K | Ha |
| 602 kPa | 2,80 | 1 753 K | 1 586 K | 647 K | Yo‘q |
| 608 kPa | 2,60 | 1 832 K | 1 656 K | 687 K | Yo‘q |
| 614 kPa | 2,40 | 1 940 K | 1 750 K | 767 K | Yo‘q |
| 610 kPa | 2,20 | 2 034 K | 1 828 K | 784 K | Ha |
Aralashma boyishi bilan qoldiq gaz harorati ko‘tarilib bordi. Chaqnash yuz bergan nuqtalarda \(T_{EVO}\) 784–791 K ga yetdi. Bu vodorodning o‘z-o‘zidan alangalanish harorati (773–793 K) bilan to‘liq mos tushadi.
Chiqarish klapani va sham atrofidagi mahalliy qizigan zonalar o‘rtacha haroratdan 100–200 K issiqroq bo‘lishi mumkinligi qayd etilgan. Mexanizm zanjiri quyidagicha:
- λ kamayishi bilan aralashma boyiydi va yonish harorati ko‘tariladi.
- Silindrda qolgan qoldiq gazlar yanada qiziydi.
- Gaz almashinuvi paytida issiq qoldiq gaz yoki qizigan sirt kiritish tizimiga tegadi.
- Kollektorda tayyor turgan vodorod-havo aralashmasi o‘z-o‘zidan alangalanadi.
- Alanga kollektor bo‘ylab orqaga qarab portlaydi.
Detektorning ishlash ko‘rsatkichlari
| Ko‘rsatkich | Normal sikllar | Chaqnash sikli |
|---|---|---|
| CUSUM cho‘qqi qiymati | taxminan 10 atrofida | taxminan 570–1 100 |
| Spektral basiklik | taxminan 0–5 | taxminan 16–20 |
| Normallashtirilgan bosim ortishi | Normal taqsimot | 23 standart chetlanishdan yuqori |
| Spektral quvvat o‘sishi | Baza daraja | 20–40 dB yuqoriroq |
CUSUM chegarasi kiritish taktida (−273 dan −214 °aTDC gacha) buzilgan. Ogohlantirish signali kollektor bosimi cho‘qqiga chiqishidan 4–16 ms (1 800 ayl/daqiqada 44–172° krank burchagi) oldin shakllangan.
Bu vaqt oralig‘i quyidagi himoya choralarini ko‘rish imkonini berishi mumkin:
- Vodorod uzatishni to‘xtatish,
- Uchqunni o‘chirish yoki kechiktirish,
- Drosselni yoki xavfsizlik klapanini boshqarish,
- Kollektor bosimini bo‘shatish,
- Dvigatel yuklamasini zudlik bilan tushirish.
Tadqiqotning kuchli tomonlari
- Teskari chaqnashni haqiqiy og‘ir yuk dvigateli o‘lchamida tadqiq etgan.
- Mikserli tizimning katta kollektor zaryadi xavfini bevosita o‘rgangan.
- 200 kHz datchik orqali vaqt va chastota axborotini birga tahlil qilgan.
- Aniqlash, shiddat o‘lchovi va mexanizmni yagona tajriba bazasida birlashtirgan.
- Mutlaq bosim o‘rniga normallashtirilgan moslashuvchan parametrlardan foydalangan.
- Ogohlantirishning bosim cho‘qqisidan oldin paydo bo‘lish vaqtini aniq ko‘rsatgan.
- Hajmga bog‘liq bo‘lmagan yonish to‘liqligi (\(\eta\)) mezonini ishlab chiqqan.
- Qoldiq gaz harorati bilan chaqnash chegarasi orasidagi moslikni tasdiqlagan.
Tadqiqotning cheklovlari
- Bor-yo‘g‘i 5 ta teskari chaqnash hodisasi o‘rganilgan.
- Barcha tajribalar faqat 1 800 ayl/daqiqada o‘tkazilgan; boshqa tezliklarda tekshirilmagan.
- Motor 6 silindrli dvigatelning bitta faol silindrli laboratoriya shaklidir.
- Seriyali ECU da real vaqtda sinov o‘tkazilmagan.
- Hisoblash yuki va protsessor vaqti hisobot qilinmagan.
- Qoldiq gaz harorati to‘g‘ridan-to‘g‘ri datchik bilan o‘lchanmagan, hisoblangan.
- Chaqnash siklida tutilgan massa yaqin rejimdan masshtablangan.
- Xom ma’lumotlar va kodlar ochiq taqdim etilmagan.
- Tadqiqot taqrizdan o‘tmagan.
Tadqiqot nimani tasdiqlaydi?
- Mikserli vodorod dvigatelida teskari chaqnash kollektorda katta bosim zarbasi hosil qilishini tasdiqlaydi.
- Chaqnash keng polosali va vaqt bo‘yicha lokalizatsiyalangan spektral imzo berishini ko‘rsatadi.
- Yuqori chastotali CUSUM usuli chaqnashni bosim cho‘qqisidan 4–16 ms oldin aniqlay olishini isbotlaydi.
- Aralashma boyishi bilan qoldiq gaz harorati vodorodning alangalanish haroratiga yetishini ko‘rsatadi.
- Erta o‘t oldirish detonatsiya, kech o‘t oldirish chaqnash chegarasini belgilashini tasdiqlaydi.
- λ ≈ 3 va o‘rta yuklama zonasida yuqori sof FIK (48%) bilan chaqnashsiz ishlash mumkinligini ko‘rsatadi.
Tadqiqot nimani isbotlamaydi?
- Detektor barcha vodorod dvigatellarida aynan shu parametrlar bilan ishlashini kafolatlamaydi.
- Oddiy arzon MAP datchigi laboratoriya datchigi kabi natija berishini isbotlamaydi.
- Chaqnash har doim faqat qoldiq gazdan kelib chiqishini isbotlamaydi (moy tomchisi, uglerod qoldig‘i kabi boshqa manbalar ham bo‘lishi mumkin).
- Signal chiqqandan keyin yoqilg‘ini to‘xtatish mexanik shikastlanishni amalda to‘liq bartaraf etishini eksperimental ko‘rsatmagan.
- 48% sof indikatorli FIK tormoz yoki avtomobil FIK ekanini anglatmaydi.
Tadqiqot usuli va natijalari
Texnik uslubiyot xulosasi
| Uslubiy element | Tadqiqotda qo‘llanilgan yondashuv |
|---|---|
| Tadqiqot turi | Eksperimental dvigatel va signallarni qayta ishlash tadqiqoti |
| Dvigatel | Bitta faol silindrli, og‘ir yuk vodorodli uchqunli dvigatel |
| Silindr hajmi | 1 841 cm³ |
| Siqish darajasi | 13:1 |
| Yoqilg‘i uzatish | Kiritish klapanidan oldin tashqi aralashtirish (mikser) |
| Tezlik | 1 800 ayl/daq |
| Yuklama | taxminan 200–800 kPa yalpi IMEP |
| λ oralig‘i | 1,69–4,40 |
| O‘t oldirish burchagi | −32 °aTDC dan +10 °aTDC gacha |
| Asosiy datchik | Kistler 6052C tezkor kollektor bosimi datchigi (200 kHz) |
| Filtr | 500 Hz, 2-tartibli Batervort yuqori o‘tkazuvchi filtri |
| Asosiy detektor | EWMA standartlashtirishi va Page–Hinkley/CUSUM |
| CUSUM parametrlari | k = 1, h = 10 |
| Shiddat mezoni | O‘lchangan bosim ortishining to‘liq yonish chegarasiga nisbati |
| Termik mexanizm | Silindr bosimidan hisoblangan Tb va TEVO |
| Tahlil qilingan chaqnashlar | Batafsil ko‘rilgan 5 ta hodisa |
Asosiy natijalar
| Ko‘rsatkich | Natija | Izoh |
|---|---|---|
| Maksimal kollektor bosimi | 251–635 kPa | Kollektorda og‘ir bosim yuki hosil bo‘lgan |
| Yonish to‘liqligi (η) | %79,6–89,3 | Kollektor zaryadining katta qismi yongan |
| Qoldiq gaz harorati | Chaqnashda 784–791 K | Vodorodning alangalanish harorati bilan mos |
| Spektral quvvat o‘sishi | taxminan 20–40 dB | Keng polosali zarba xususiyati |
| CUSUM chaqnash bali | taxminan 570–1 100 | Normal sikllardan (≈10) ancha yuqori |
| Aniqlash vaqti yutug‘i | 4–16 ms | Bosim cho‘qqisiga chiqmasdan oldin signal berilgan |
| Maksimal sof FIK | %48,1–48,2 | λ ≈ 3 va o‘rta yuklamada |
| Kambag‘al rejimdagi NOx | λ ≥ 3,2 da 0–1 ppm | Dvigateldan to‘g‘ridan-to‘g‘ri o‘lchangan |
Manba va metodologiya eslatmasi
Tadqiqotning asl nomi: Backfire in a heavy-duty mixer-type hydrogen spark-ignition engine: detection, severity quantification, and the residual-gas thermal mechanism
Mualliflar va ularning tartibi: Seongsu Kim; Hyunwook Park; Seungmook Oh; Changup Kim; Yonggyu Lee; Junsun Lee; Junghwan Kim.
Teng hissa: PDF da teng birinchi mualliflik bayonoti berilmagan.
Mas’ul muallif: Junghwan Kim (jkim77@cau.ac.kr).
Muassasalar:
- Energiya tizimlari muhandisligi kafedrasi, Chung-Ang universiteti, Seul, Koreya Respublikasi.
- Dvigatel tadqiqotlari bo‘limi, Koreya mashinasozlik va materiallar instituti (KIMM), Tedjon, Koreya Respublikasi.
DOI: 10.2139/ssrn.6945147
Manba turi: Og‘ir yuk vodorod dvigatelida o‘tkazilgan eksperimental tadqiqot va signallarga ishlov berish preprinti.
Taqriz holati: Tadqiqot ilmiy taqrizdan o‘tmagan.
Nashr platformasi: SSRN.
SSRN sanasi: 2026-yil 15-iyun.
Sahifalar soni: 37.
Rasmiy havolalar:SSRN maqola sahifasi va DOI havolasi.
Moliyalashtirish: Koreya Savdo, sanoat va energetika vazirligi huzuridagi Texnologik innovatsiyalar dasturi (20018148) tomonidan moliyalashtirilgan.
Manfaatlar to‘qnashuvi: Mualliflar manfaatlar to‘qnashuvi yo‘qligini bildirganlar.
Ushbu ilmiy sharh yuklangan 37 sahifalik PDF hujjatiga asoslanib tayyorlandi. PDF da bo‘lmagan yangi eksperimental natija yoki tijoriy kafolat qo‘shilmadi. Tashqi manbalardan faqat DOI, SSRN sanasi va muassasa ma’lumotlarini tekshirish uchun foydalanildi.
PDF dagi ichki tafovutlar va noaniqliklar
- Tezlikdan mustaqillik da’vosi: Matnda chegara motor tezligiga bog‘liq emas deyilgan, ammo barcha tajribalar faqat 1 800 ayl/daqiqada o‘tkazilgan.
- Z-skor formulasi: 6-formulada farq dispersiyaga bo‘lingan; standart statistikada u standart chetlanishga bo‘linishi kerak.
- Spektral tekislik talqini: Geometrik/arifmetik o‘rtacha nisbati tekis spektrda 1 ga yaqinlashadi, matnda esa teskarisi aytilgan.
- STFT parametrlari: Oyna uzunligi va qoplash darajasi ko‘rsatilmay, joyida qoralama matn qolib ketgan.
- Chastota polosasi tafovuti: Turli bo‘limlarda 50 Hz–26 kHz, 200 Hz–26 kHz va 0,5–26 kHz oraliqlari aralash berilgan.
- Cho‘qqi bosimi hisoboti: Xulosada 250–460 kPa deyilgan, 7-jadvalda esa 547 va 635 kPa qiymatlari ham bor.
- Tahrir xatosi: Algoritm tushuntirishining oxirida gap o‘rtasidan boshlangan takroriy matn parchasi qolgan.
Ushbu kamchiliklar kollektor bosimi orqali teskari chaqnashni erta aniqlash g‘oyasining ilmiy ahamiyatini yo‘qqa chiqarmaydi, biroq tizimni sanoatda qo‘llashdan oldin turli tezliklarda va motorlarda qo‘shimcha sinovlar o‘tkazilishi shart.

Izoh qoldiring
E-pochta manzilingiz chop etilmaydi. Majburiy maydonlar * bilan belgilangan