[ Flash Home ]

[ Validation and
Basic Physics
Group Home
]

[ Overview ]

[ Gallery ]

[ Publications ]

DNS Studies of Burning Fronts in
Passive-Reactive Diffusion

N.Vladimirova, F.Cattaneo, A.Malagoli, A.Oberman, R.Rosner, O.Ruchayskiy

Problem 3: Quenching by Cellular Flow, Le=1

Initial Temperature The initial band of hot material is distorted by cellular flow. Reaction term is step function of temperature, so that no burning occur if temperature is lower than critical. For temperatures higher than critical burning rate is constant and chosen to make laminar speed equal to unity.

colormap
colormap


Some Visual Observations

  • Quenching by cellular flow requires much higher velocities compared to shear flow, but it is possible.

  • The process of flame propagation is as follows. In any unburned row, hot material from a neighboring row is quickly dispersed, uniformly preheating the cell. Before ignition in that cell, preheating of the next cell begins.

  • Quenching occurs if a row of unburned cells loses too much heat to its cooler neighbors before igniting.

  • It takes significantly longer to ignite the first row, because it has not been preheated (note the logarithmic time scale for images below).

  • Quenching is not observed if the first unburned row ignites.

  • We did not investigate the large cell size limit, corresponding to the geometric optics regime. In that regime for shear flow, the band may break into isolated flamelets, which may or may not reconnect later, depending on U. In cellular flow we did not observe isolated flamelets.

Evolution of flame with the amplitude of shear above critical.
U=13
Temperature movie dt=0.0100, tmax=3.00:   quicktime (1.3M),   avi (1.2M)
Reaction rate movie dt=0.0100, tmax=3.00:   quicktime (1.4M),   avi (1.0M)
Temperature movie dt=0.0001, tmax=0.03:   quicktime (1.6M),   avi (1.4M)
Reaction rate movie dt=0.0001, tmax=0.03:   quicktime (1.5M),   avi (1.2M)

Evolution of flame with the amplitude of shear below critical.
U=14
Temperature movie dt=0.0100, tmax=3.00:   quicktime (0.7M),   avi (0.6M)
Reaction rate movie dt=0.0100, tmax=3.00:   quicktime (0.7M),   avi (0.5M)
Temperature movie dt=0.0001, tmax=0.03:   quicktime (1.6M),   avi (1.5M)
Reaction rate movie dt=0.0001, tmax=0.03:   quicktime (1.4M),   avi (1.0M)


Critical velocity as function of initial band width
Ucrit vs W plot

Critical amplitude of the shear flow velocity is proportional to the forth power of initial band width. The computer simulations were suggested by P. Constantin, A. Kiselev and L. Ryzhik, and show good agreement with their theory.

The dependence of critical velocity on initial band width can be explained from the point of view of the effective flame. For the flame to be quenched, band width should be of the order of effective flame thickness. Effective flame thickness is proportional to the U1/4, which leads to the critical velocity proportional to the forth power of initial band width.


Advection-reaction-diffusion problems:

Problem 1: Burning in Shear Flow, Le=1, KPP
Problem 2: Burning in Cellular Flow, Le=1, KPP
Problem 3: Quenching by Shear Flow, Le=1, step function reaction
Problem 4: Quenching by Cellular Flow, Le=1, step function reaction
Problem 5: Quenching by Shear Flow, Le>1, step function reaction

[Validation and Basic Physics Group Home] [Overview] [Gallery] [Publications]
[Flash Home]


nata@flash.uchicago.edu
This file was last modified on 8 March 2001.
The ASCI Flash Center is based at the University of Chicago under U. S. Department of Energy contract B341495. All material on these pages is Copyright © 1999 ASCI Flash Center.