Płomienie Przedprzapalne: Praca Sprzężenia Reakcji-Difuzji

Ogień to samonaprawiająca się fala reakcji-dyfuzji

Zapala się mieszanka wstępna – paliwo i utleniacz już zmieszane przed zapłonem, jak w piecu Bunsenowym lub skompresowanym ładunku silnika wewnętrznego spalania – rozprzestrzeniając się, ponieważ ciepło i reaktywne radykały dyfuzyjnie rozchodzą się z gorącej, spalonej strefy do zimniejszej, niespalonej mieszanki tu przed nią, podnosząc temperaturę tej ostatniej, aż do jej zapłonu. Cała frontalna część wtedy porusza się ze stałą prędkością, ustalona w wyniku równowagi między szybkością dyfuzji ciepła i szybkością zużycia paliwa przez reakcję chemiczną, gdy tylko dotrze – dokładnie taka sama równowaga reakcji-dyfuzji, która rządzi propagacją impulsu nerwowego lub rozprzestrzenianiem się epidemii w populacji, z tym samym mechanizmem spalania zastępującym człon reakcji.

Rządy równania

Minimalny model łączy pole temperatury (lub zmienne postępu) z polem koncentracji paliwa poprzez współdzielony, silnie zależny od temperatury termin reakcji:

∂T/∂t = α∇²T + (Q/cp) · ω(T, Y) ∂Y/∂t = D∇²Y - ω(T, Y) ω(T,Y) = A · Y · exp( -Ea / (R·T) ) (Szybkość reakcji Arrhenius) T = temperatura, Y = masa paliwa w procentach α = ciepło właściwe, D = dyfuzja gatunkowa Q = uwalniana energia na jednostkę reakcji, cp = ciepło właściwe A = czynnik bezwzględny (częstotliwość), Ea = energia aktywacji

Termin reakcji ω jest kluczowy dla całego problemu: ponieważ zależy w sposób wykładniczy od temperatury, jest zasadniczo wyłączony w zimnym, nie spalonym gazie i natychmiast włącza się dopiero po tym, jak dyfuzja wstępnie podgrzeje cienką strefę przed frontem — ta ekstremalna wrażliwość na temperaturę jest dokładnie powodem, dla którego płomień utrzymuje się jako cienka, dobrze zdefiniowana front, a nie rozmyta plama, oraz dlaczego samo zapłon zachowuje się prawie jak fenomen umowny.

∂T/∂t = α∇²T + (Q/cp) · ω(T, Y)
∂Y/∂t = D∇²Y - ω(T, Y)

ω(T,Y) = A · Y · exp( -Ea / (R·T) )      (Arrhenius reaction rate)

T  = temperature, Y = fuel mass fraction
α  = thermal diffusivity, D = species diffusivity
Q  = heat release per unit reaction, cp = specific heat
A  = pre-exponential (frequency) factor, Ea = activation energy

Laminarne prędkość płomienia: jak szybko front zjada mieszaninę

Dla prostego laminarnego, wstecznie zmiksowanego płomienia, klasyczny argument skalujący (równoważający strefę upłynnionego podgrzewania z strefą reakcji) daje prędkość płomienia jako w przybliżeniu:

S_L ≈ sqrt( α · A · exp(-Ea/RT_ad))

S_ad = temperatura płomienia adiabatycznego osiągnięta przez pełne spalanie

Zależność pierwiastkowa od Arrheniusa jest powodem, dla którego prędkość płomienia jest tak wrażliwa na skład powietrza i temperaturę początkową: mieszanki oddalone od stechiometrii spalają się znacznie chłodniej (niższa T_ad), co drastycznie tłumi wyrazu potęgowego i może spowolnić lub całkowicie ugaszyć płomień – podstawa granic palności, zakres składu powietrza poza którym płomień po prostu nie rozprzestrzenia się.

S_L ≈ sqrt( α · A · exp(-Ea/RT_ad) )

T_ad = adiabatic flame temperature reached by full combustion

Turbulencja: zaginanie krawędzi, nie zmieniająca się chemii

Rzeczywiste spalanie – w cylindrze silnika spalinowego, turbinie gazowej, palniku przemysłowym – jest turbulentne, a turbulencja zazwyczaj nie przyspiesza bezpośrednio podstawowej szybkości reakcji chemicznej. Zamiast tego, zaginia i rozciąga cienką, warstwową krawędź płomienia w złożoną powierzchnię, znacznie zwiększając całkowitą powierzchnię krawędzi w danym wolumenie. Skala spalania zależy od powierzchni płomienia, dlatego skuteczna turbulencyjna prędkość spalania S_T może przekroczyć prędkość warstwową S_L o duży czynnik, często modelowaną empirycznie jako S_T/S_L rosnącą wraz z intensywnością turbulencji u′ w stosunku do S_L samej – szybsza turbulencja bardziej agresywnie zaginia płomień, aż do momentu, w którym intensywność turbulencji lokalnie gaśnie lub pogrubia strefę reakcji, zamiast jej jedynie zaginąć.

Dlaczego symulacja tego jest naprawdę trudna

Rejon reakcji w prawdziwym płomieniu często ma grubość mniejszą niż jeden milimetr, podczas gdy turbulencyjne języki, które go zniekształcają, mogą sięgać metrów – to rozdzielenie skali o co najmniej trzech rzędach wielkości, które sprawia, że rozwiązywanie obu jednocześnie w pojedynczej symulacji jest dla większości praktycznych geometrii inżynieryjnych zbyt kosztowne. Dlatego modelowanie spalania opiera się na dużym stopniu na modele redukowane: podejście z płomieniami, które traktuje spękaną turbulencyjną żarę jako zbiór lokalnie laminarnych elementów płomieniowych, oraz zamknięcia szybkości reakcji, które szacują wpływ terminu Arrheniusa bez rozwiązywania pełnej wieloosobowej kinetyki chemicznej, która rzeczywiście występuje w prawdziwym paliwie, takim jak benzyna lub paliwo lotnicze (często setki gatunków i tysiące podstawowych reakcji w mechanizmie).

Frequently asked questions

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Premixed Flame Front i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.