About the Premixed Flame Front Simulation
This simulation solves the one-dimensional reaction-diffusion equation ∂T/∂t = D·∂²T/∂x² + ω(T) for the normalised temperature T (where T = 0 represents cold unburnt reactants and T = 1 represents the adiabatic flame temperature). The source term ω(T) = A·exp(−Ea/RT) is the Arrhenius reaction rate, capturing the sharp exponential sensitivity of combustion chemistry to temperature. At low temperatures the rate is negligible; once the gas is heated past the ignition threshold the rate increases explosively, converting reactants to hot products and sustaining the propagating flame front.
The numerical scheme uses a 1D finite-difference grid of N = 700 cells with Neumann (zero-flux) boundary conditions. An explicit Euler time-stepper advances the solution with step dt chosen to satisfy the diffusion stability limit dt ≤ dx²/(2D). Several sub-steps are taken per animation frame so the simulation runs smoothly at 60 fps while maintaining accuracy. The main canvas shows the temperature profile T(x) as a glowing filled area coloured from deep red (cold) through orange-white (peak flame). A kymograph inset records successive snapshots, whose diagonal slope gives the flame speed directly. The HUD displays the measured flame speed, theoretical laminar flame speed, flame thickness, and peak temperature in real time.
Frequently Asked Questions
What is a premixed flame front?
A premixed flame front is a thin reaction zone in which fuel and oxidizer, mixed before ignition, react and release heat. The hot combustion products heat adjacent unburnt gas by thermal conduction, which in turn ignites it — creating a self-sustaining propagating wave. In a 1D model this appears as a sharp travelling temperature profile.
What is the reaction-diffusion equation for flame propagation?
The governing equation is ∂T/∂t = D∇²T + ω(T), where T is the normalised temperature (0 = cold reactants, 1 = adiabatic products), D is the thermal diffusivity, and ω(T) is the Arrhenius reaction rate. The diffusion term D∇²T transports heat ahead of the front; the source term ω(T) releases chemical energy.
What is the Arrhenius reaction rate?
The Arrhenius rate is ω = A·exp(−Ea/RT), where A is the pre-exponential frequency factor, Ea is the activation energy, R is the gas constant, and T is temperature. At low temperatures the exponential is tiny and the reaction is frozen; once T rises past the ignition threshold the rate increases explosively, producing the sharp heat-release zone.
How is the laminar flame speed derived?
The laminar flame speed s_L = √(D·A) · exp(−Ea / 2R·T₀) comes from a travelling-wave analysis of the reaction-diffusion equation. Balancing diffusive heat transport with Arrhenius heat release gives this square-root dependence on diffusivity D and pre-exponential factor A, and an exponential sensitivity to activation energy Ea and inlet temperature T₀.
What does activation energy Ea control in this simulation?
Activation energy Ea sets how strongly temperature-sensitive the chemistry is. High Ea produces a very steep flame front that only reacts in a narrow high-temperature zone (thin flame). Low Ea allows the reaction to proceed at lower temperatures, broadening the front and increasing the flame speed. In real flames, Ea/R is typically 10 000–20 000 K.
Why does increasing thermal diffusivity D make the flame faster?
Thermal diffusivity D measures how quickly heat spreads from the hot reaction zone into the cold unburnt gas. A higher D pre-heats a wider layer of fresh mixture per unit time, shortening the ignition delay and accelerating the flame. The formula s_L ∝ √D captures this square-root sensitivity.
What is the preheat zone?
The preheat zone is the region immediately upstream of the reaction zone where cold reactants are warmed by thermal conduction from the hot flame. In the 1D model it appears as the gently rising temperature ramp in front of the sharp heat-release peak. Its thickness scales as δ ∼ D / s_L — thinner for faster-burning mixtures.
How does initial temperature affect flame propagation?
Raising the initial reactant temperature T₀ reduces the energy needed to ignite the mixture, effectively lowering the apparent activation energy. The laminar flame speed increases exponentially with T₀ through the Arrhenius factor exp(−Ea/2RT₀). This is why pre-heated intake air or exhaust gas recirculation dramatically accelerates combustion.
What is the difference between laminar and turbulent flames?
A laminar flame propagates as a smooth planar front at speed s_L. A turbulent flame has a wrinkled or corrugated front: turbulent eddies increase the flame surface area and mix reactants with hot products, raising the global burning rate to a turbulent flame speed s_T much greater than s_L. This simulation models the underlying laminar 1D dynamics that are the building block of turbulent combustion models.
Where does premixed flame propagation appear in real applications?
Premixed flames occur in spark-ignition (petrol) engines, gas-turbine combustors running in premixed mode, domestic burners, and hydrogen fuel cells. Understanding s_L and flame thickness is critical for designing efficient, low-emission combustors, predicting engine knock, and preventing flashback or blowout in gas turbines.