Combustion Science and Technology, Vol.128, No.1-6, 99-130, 1997
One-dimensional simulations of freely propagating turbulent premixed flames
The propagation rate and the structure of freely propagating premixed turbulent flames are investigated using a one-dimensional simulation model based on a new version of the linear-eddy model (LEM) of Kerstein (1991, 1992). This model explicitly includes thermo-diffusive, finite-rate kinetic, and heat release effects. Reasonably good quantitative agreement in predictions of turbulent flame speed with fan-stirred bomb experiments of Abdel-Gayed et al. (1984a) is obtained over most of the reported u'/S-L range. LEM predicts a rapid increase in u(t)/S-L with u'/S-L for low u' followed by a bending slope of u(t)/S-L with increasing u' that was also observed in the experiments. Here. u(t) and S-L are, respectively, the turbulent and stretch free planar laminar flame speeds and u' is the r.m.s, turbulence intensity. Comparisons with an earlier model based on the G-equation (Menon and Kerstein, 1992) for flamelet combustion are also made. The resulting propagation speeds are also in good agreement. Comparisons with weak-swirl burner experiments of stationary flames by Bedat and Cheng (1995) show that the model underpredicts the reported u(t)/S-L with u'/S-L. However, progress variable probability density functions at different locations within the flame reveal the onset of distributed combustion which is predicted by the location of the flame on the Borghi combustion phase diagram (Bedat and Cheng, 1995). Finally, constant Reynolds number simulations for a range in S-L/u' compare well with experiments by Abdel-Gayed er al. (1979) for low u', but predict a plateau in u(t)/S-L as u increases, and decreasing u(t)/S-L with further decrease in S-L/u'. This behavior is interpreted as correct based on physical arguments.