Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reacting flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

Bibliography of Doctor Chul Park

This document contains a comprehensive bibliography of the published works, and a short biography, of Dr. Chul Park. The contents of this bibliography were compiled primarily from the NASA RECON data base. The RECON citations have been modified to appear in a uniform format with all other listed citations . These other citations were located by computer searches in the INSPEC, NTIS, COMPENDEX, and Chemical Abstracts data bases, as well as through the cooperation of Dr. Chul Park, and his associates in the Reacting Flow Environments Branch at NASA Ames Research Center. All citations are presented in an approximate reverse chronological order from the present date. This work was created to honor the occasion of Dr. Chul Park's retirement on December 14, 1994, after 27 years of distinguished government service at the NASA Ames Research Center.

Gochberg, Lawrence A.↗

Reliability enhancement of Navier-Stokes codes through convergence acceleration

Methods for enhancing the reliability of Navier-Stokes computer codes through improving convergence characteristics are presented. The improving of these characteristics decreases the likelihood of code unreliability and user interventions in a design environment. The problem referred to as a 'stiffness' in the governing equations for propulsion-related flowfields is investigated, particularly in regard to common sources of equation stiffness that lead to convergence degradation of CFD algorithms. Von Neumann stability theory is employed as a tool to study the convergence difficulties involved. Based on the stability results, improved algorithms are devised to ensure efficient convergence in different situations. A number of test cases are considered to confirm a correlation between stability theory and numerical convergence. The examples of turbulent and reacting flow are presented, and a generalized form of the preconditioning matrix is derived to handle these problems, i.e., the problems involving additional differential equations for describing the transport of turbulent kinetic energy, dissipation rate and chemical species. Algorithms for unsteady computations are considered. The extension of the preconditioning techniques and algorithms derived for Navier-Stokes computations to three-dimensional flow problems is discussed. New methods to accelerate the convergence of iterative schemes for the numerical integration of systems of partial differential equtions are developed, with a special emphasis on the acceleration of convergence on highly clustered grids.

Merkle, Charles L.↗

Heat release effects on the instability of parallel shear layers

The influence of time-dependent heat addition on the linear instablity of shear layers is of considerable interest in understanding the dynamic behavior of reacting flows and combustion-turbulence interactions. The approach is based upon the Bernoulli enthalpy aeroacoustics theory, which utilizes the specific enthalpy and specific entropy as the primary thermodynamic variables. In addition, velocity oscillations are split into Helmoholtz decomposition theorem.

Hegde, U.↗

Spurious Numerical Solutions Of Differential Equations

Paper presents detailed study of spurious steady-state numerical solutions of differential equations that contain nonlinear source terms. Main objectives of this study are (1) to investigate how well numerical steady-state solutions of model nonlinear reaction/convection boundary-value problem mimic true steady-state solutions and (2) to relate findings of this investigation to implications for interpretation of numerical results from computational-fluid-dynamics algorithms and computer codes used to simulate reacting flows.

Lafon, A.↗

Dynamics and structure of turbulent premixed flames

In earlier work (Mantel & Bilger, 1994) the structure of the turbulent premixed flame was investigated using statistics based on conditional averaging with the reaction progress variable as the conditioning variable. The DNS data base of Trouve and Poinsot (1994) was used in this investigation. Attention was focused on the conditional dissipation and conditional axial velocity in the flame with a view to modeling these quantities for use in the conditional moment closure (CMC) approach to analysis of kinetics in premixed flames (Bilger, 1993). Two remarkable findings were made: there was almost no acceleration of the axial velocity in the flame front itself; and the conditional scalar dissipation remained as high, or higher, than that found in laminar premixed flames. The first finding was surprising since in laminar flames all the fluid acceleration occurs through the flame front, and this could be expected also for turbulent premixed flames at the flamelet limit. The finding gave hope of inventing a new approach to the dynamics of turbulent premixed flames through use of rapid distortion theory or an unsteady Bernoulli equation. This could lead to a new second order closure for turbulent premixed flames. The second finding was contrary to our measurements with laser diagnostics in lean hydrocarbon flames where it is found that conditional scalar dissipation drops dramatically below that for laminar flamelets when the turbulence intensity becomes high. Such behavior was not explainable with a one-step kinetic model, even at non-unity Lewis number. It could be due to depletion of H2 from the reaction zone by preferential diffusion. The capacity of the flame to generate radicals is critically dependent on the levels of H2 present (Bilger, et al., 1991). It seemed that a DNS computation with a multistep reduced mechanism would be worthwhile if a way could be found to make this feasible. Truly innovative approaches to complex problems often come only when there is the opportunity to work close at hand with the (in this case numerical) experimental data. Not only can one spot patterns and relationships in the data which could be important, but one can also get to know the limitations of the technique being used, so that when the next experiment is being designed it will address resolvable questions. A three-year grant from the Australian Research Council has enabled us to develop a small capability at the University of Sydney to work on DNS of turbulent reacting flow, and to analyze data bases generated at CTR. Collaboration between the University of Sydney and CTR is essential to this project and finding a workable modus operandum for this collaboration, given the constraints involved, has been a major objective of the past year's effort. The overall objectives of the project are: (1) to obtain a quantitative understanding of the dynamics of turbulent premixed flames at high turbulence levels with a view to developing improved second order closure models; and (2) to carry out new DNS experiments on turbulent premixed flames using a carefully chosen multistep reduced mechanism for the chemical kinetics, with a view to elucidating the laser diagnostic findings that are contrary to the findings for DNS using one-step kinetics. In this first year the objectives have been to make the existing CTR data base more accessible to coworkers at the University of Sydney, to make progress on understanding the dynamics of the flame in this existing CTR data base, and to carefully construct a suitable multistep reduced mechanism for use in a new set of DNS experiments on turbulent premixed flames.

Bilger, R. W.↗

Study of Turbulent Premixed Flame Propagation using a Laminar Flamelet Model

The laminar flamelet concept in turbulent reacting flows is considered applicable to many practical combustion systems (Linan & Williams 1993). For turbulent premixed combustion, the laminar flamelet regime is valid when turbulent Karlovitz number is less than unity, which is equivalent to stating that the characteristic thickness of the flame is less than that of a Kolmogorov eddy; this is known as the Klimov-Williams criterion (Williams 1985). In such a case, the flame maintains its laminar structure, and the effect of turbulent flow is merely to wrinkle and strain the flame front. The propagating wrinkled premixed flame can then be described as an infinitesimally thin surface dividing the unburnt fresh mixture and the burnt product.

Im, H. G.↗

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗

Study of Buoyancy Effects in Diffusion Flames Using Rainbow Schlieren Deflectometry

Diffusion flames are extensively encountered in many domestic and industrial processes. Even after many decades of research, a complete understanding of the diffusion flame structure is not available. The structure and properties of the flames are governed by the mixing (laminar or turbulent), chemical kinetics, radiation and soot processes. Another important phenomenon that affects flame structure in normal gravity is buoyancy. The presence of buoyancy has long hindered the rational understanding of many combustion processes. In gas jet diffusion flames, buoyancy affects the structure of the shear layer, the development of fluid instabilities, and formation of the coherent structures in the near nozzle region of the gas jets. The buoyancy driven instabilities generate vorticial structures outside the flame resulting in flame flicker. The vortices also strongly interact with the small-scale structures in the jet shear layer. This affects the transitional and turbulence characteristics of the flame. For a fundamental understanding of diffusion flames it is essential to isolate the effects of buoyancy. This is the primary goal of the experiments conducted in microgravity. Previous investigations, have shown dramatic differences between the jet flames in microgravity and normal gravity. It has been observed that flames in microgravity are taller and more sooty than in normal gravity. The fuels used in these experiments were primarily hydrocarbons. In the absence of buoyancy the soot resides near the flame region, which adversely affects the entrainment of reactants. It is very important to eliminate the interference of soot on flame characteristics in microgravity. The present work, therefore, focuses on the changes in the flame structure due to buoyancy without the added complexities of heterogeneous reactions. Clean burning hydrogen is used as the fuel to avoid soot formation and minimize radiative losses. Because of the low luminosity of hydrogen flames, we use rainbow schlieren deflectometry for visualization. The visualized images are digitized for quantification.The work reported here is divided into three sections; rainbow schlieren deflectometry (RSD), microgravity experiments and sub-atmospheric pressure experiments. The first section demonstrates the application of RSD for quantitative measurements in non-reacting and reacting flow systems. A computational effort to complement the experimental work is also included. In the second section, the experiments conducted at the 2.2s NASA Lewis Drop tower facility are described. The experiments were conducted to study the behavior of laminar, transitional and turbulent hydrogen flames in microgravity. The ability of RSD technique to provide quantitative data is highlighted. The final section deals with the sub-atmospheric pressure tests, which demonstrate that buoyancy in hydrogen diffusion flames can be scaled with pressure at normal gravity.

Agrawal, Ajay K.↗

The Effects of Buoyancy and Dilution on the Structure and Lift-off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of the factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. The goal of our characterizations of coflow laminar diffusion flames is to bring to microgravity the multidimensional diagnostic tools available in normal gravity, and in so doing provide a broader understanding of the successes and limitations of current combustion models. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in both buoyant and nonbuoyant environments. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2)delta, denoted CH*), and excited-state OH (A(exp 2)Sigma, denoted OH*) are measured in mu-g and normal gravity. Two-dimensional CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC-135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. Although CH* and OH* kinetics are not well understood, the CH*, OH*, and ground-state CH distributions are spatially coincident in the flame anchoring region. Therefore, the ground-state CH distribution, which is easily computed, and the readily measured CH*/OH* distributions can be used to provide a consistent and convenient way of measuring lift-off height and flame shape in the diffusion flame under investigation. Given that the fuel composition affects flame chemistry and that buoyancy influences the velocity profile of the flow, we have the opportunity to computationally and experimentally study the roles of fluids and chemistry. In performing this microgravity study, improvements to the computational model have been made and new calculations performed for a range of gravity and flow conditions. Furthermore, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented in the following sections.

Walsh, Kevin T.↗

A Pseubo-Temporal Multi-Grid Relaxation Scheme for Solving the Parabolized Navier-Stokes Equations

A multi-grid, flux-difference-split, finite-volume code, VULCAN, is presented for solving the elliptic and parabolized form of the equations governing three-dimensional, turbulent, calorically perfect and non-equilibrium chemically reacting flows. The space marching algorithms developed to improve convergence rate and or reduce computational cost are emphasized. The algorithms presented are extensions to the class of implicit pseudo-time iterative, upwind space-marching schemes. A full approximate storage, full multi-grid scheme is also described which is used to accelerate the convergence of a Gauss-Seidel relaxation method. The multi-grid algorithm is shown to significantly improve convergence on high aspect ratio grids.

Morrison, J. H.↗

Design of a Novel Gaseous Hydrogen-Oxygen Rocket Injector Element

An overview of activities supporting the design of a gaseous hydrogen-oxygen rocket injector element is presented in viewgraph form. The purpose of the research was to find a viable design for a rocket gas-gas injector that mixes fuel and oxidizer thoroughly and quickly. Computational fluid dynamics analyses were used with reacting flow to evaluate design options for mixing, temperature distribution, and combustion efficiency. A design was found that is an improvement over designs derived from liquid systems and is far better than traditional shear-coax.

Glenn, Dennis↗

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

Trace Chemistry

The goals of the trace chemistry group were to identify the processes relevant to aerosol and aerosol precursor formation occurring within aircraft gas turbine engines; that is, within the combustor, turbine, and nozzle. The topics of discussion focused on whether the chemistry of aerosol formation is homogeneous or heterogeneous; what species are important for aerosol and aerosol precursor formation; what modeling/theoretical activities to pursue; what experiments to carry out that both support modeling activities and elucidate fundamental processes; and the role of particulates in aerosol and aerosol precursor formation. The consensus of the group was that attention should be focused on SO2, SO3, and aerosols. Of immediate concern is the measurement of the concentration of the species SO3, SO2, H2SO4 OH, HO2, H2O2, O, NO, NO2, HONO, HNO3, CO, and CO2 and particulates in various engines, both those currently in use and those in development. The recommendation was that concentration measurements should be made at both the combustor exit and the engine exit. At each location the above species were classified into one of four categories of decreasing importance, Priority I through IV, as follows: Combustor exit: Priority I species - SO3:SO2 ratio, SO3, SO2, and particulates; Priority II species: OH and O; Priority III species - NO and NO2; and Priority IV species - CO and CO2. For the Engine exit: Priority I species - SO3:SO2 ratio, SO3, SO2,H2SO4, and particulates; Priority II species: OH,HO2, H2O2, and O; Priority III species - NO, NO2, HONO, and HNO3; and Priority IV species - CO and CO2. Table I summarizes the anticipated concentration range of each of these species. For particulate matter, the quantities of interest are the number density, size distribution, and composition. In order to provide data for validating multidimensional reacting flow models, it would be desirable to make 2-D, time-resolved measurements of the concentrations of the above species and, in addition, of the pressure, temperature, and velocity. A near term goal of the experimental program should be to confirm the nonlinear effects of sulfur speciation, and if present, to provide an explanation for them. It is also desirable to examine if the particulate matter retains any sulfur. The recommendation is to examine the effects on SOx production of variations in fuel-bound sulfur and aromatic content (which may affect the amount of particulates formed). These experiments should help us to understand if there is a coupling between particulate formation and SO, concentration. Similarly, any coupling with NOx can be examined either by introducing NOx into the combustion air or by using fuel-bound nitrogen. Also of immediate urgency is the need to establish and validate a detailed mechanism for sulfur oxidation/aerosol formation, whose chemistry is concluded to be homogeneous, because there is not enough surface area for heterogeneous effects. It is envisaged that this work will involve both experimental and theoretical programs. The experimental work will require, in addition to the measurements described above, fundamental studies in devices such as flow reactors and shock tubes. Complementing this effort should be modeling and theoretical activities. One impediment to the successful modeling of sulfur oxidation is the lack of reliable data for thermodynamic and transport properties for several species, such as aqueous nitric acid, sulfur oxides, and sulfuric acid. Quantum mechanical calculations are recommended as a convenient means of deriving values for these properties. Such calculations would also help establish rate constants for several important reactions for which experimental measurements are inherently fraught with uncertainty. Efforts to implement sufficiently detailed chemistry into computational fluid dynamic codes should be continued. Zero- and one-dimensional flow models are also useful vehicles for elucidating the minimal set of species and reactions that must be included in two- and three-dimensional modeling studies.

Radhakrishnan, Krishnan↗

The Effects of Buoyancy and Dilution on the Structure and Lift-Off of Coflow Laminar Diffusion Flames

The ability to predict the coupled effects of complex transport phenomena with detailed chemical kinetics in diffusion flames is critical in the modeling of turbulent reacting flows and in understanding the processes by which soot formation and radiative transfer take place. In addition, an understanding of those factors that affect flame extinction in diffusion flames is critical in the suppression of fires and in improving engine efficiency. A goal of this work is to bring to microgravity flame studies the detailed experimental and numerical tools that have been used to study ground-based systems. This will lead to a more detailed understanding of the interaction of convection, diffusion and chemistry in a nonbuoyant environment. To better understand these phenomena, experimental and computational studies of a coflow laminar diffusion flame have been carried out. To date, these studies have focused on a single set of flow conditions, in which a nitrogen-diluted methane fuel stream (65% methane by volume) was surrounded by an air coflow, with exit velocities matched at 35 cm/s. Of particular interest is the change in flame shape due to the absence of buoyant forces, as well as the amount of diluent in the fuel stream and the coflow velocity. As a sensitive marker of changes in the flame shape, the number densities of excited-state CH (A(exp 2 delta) denoted CH*), and excited-state OH (A(exp 2 sigma, denoted OH*) are measured. CH* and OH* number densities are deconvoluted from line-of-sight chemiluminescence measurements made on the NASA KC135 reduced-gravity aircraft. Measured signal levels are calibrated, post-flight, with Rayleigh scattering. In extending the study to microgravity conditions, improvements to the computational model have been made and new calculations performed for a range of gravity conditions. In addition, modifications to the experimental approach were required as a consequence of the constraints imposed by existing microgravity facilities. Results from the computations and experiments are presented.

Walsh, Kevin T.↗

Experimental and Modeling Study of the Burning of an Ethanol Droplet in Microgravity

The microgravity ethanol droplet combustion experiments were performed aboard the STS-94/MSL-1 Shuttle mission within the Fiber-Supported Droplet Combustion-2 (FSDC-2) program. The burning histories and flame standoffs for pure ethanol and ethanol/water droplets were obtained from the images recorded with two 8 mm videocameras. The obtained results show that average gasification rate is related to the initial droplet size in a manner similar to n-alkanes and methanol and consistent with the results of Hara and Kumagai and the data taken recently in the NASA-Lewis 2.2 s droptower. A transient, moving finite-element chemically reacting flow model applied previously to sphero-symmetric combustion of methanol, methanol/water, n-alkane, and n-alkane binary mixture droplets was adopted for the problem of ethanol droplet combustion. The model includes detailed description of gas-phase reaction chemistry and transport, a simplified description of liquid phase transport, and non-luminous radiative heat transfer. Gas-phase chemistry was described with the detailed reaction mechanism of Norton and Dryer, which consists of 142 reversible elementary reactions of 33 species. Another recently published reaction mechanism of high-temperature ethanol oxidation was also considered. The model predictions were found to compare favorably with the experimental data. The model analysis also indicates that water condensation in the case of ethanol has smaller effect on average droplet gasification rate as compared with previously studied methanol cases. This effect is explained by non-ideal (azeotropic) behavior of binary ethanol-water mixtures. Further analysis of computational results and ethanol droplet radiative extinction behavior will be discussed.

Kazakov, Andrei↗

Computational Modeling in Structural Materials Processing

High temperature materials such as silicon carbide, a variety of nitrides, and ceramic matrix composites find use in aerospace, automotive, machine tool industries and in high speed civil transport applications. Chemical vapor deposition (CVD) is widely used in processing such structural materials. Variations of CVD include deposition on substrates, coating of fibers, inside cavities and on complex objects, and infiltration within preforms called chemical vapor infiltration (CVI). Our current knowledge of the process mechanisms, ability to optimize processes, and scale-up for large scale manufacturing is limited. In this regard, computational modeling of the processes is valuable since a validated model can be used as a design tool. The effort is similar to traditional chemically reacting flow modeling with emphasis on multicomponent diffusion, thermal diffusion, large sets of homogeneous reactions, and surface chemistry. In the case of CVI, models for pore infiltration are needed. In the present talk, examples of SiC nitride, and Boron deposition from the author's past work will be used to illustrate the utility of computational process modeling.

Meyyappan, Meyya↗

Prebiotic Processing induced by Comet and Meteor Impact

In their study of organic synthesis from impact shocks using the laser-induced-plasma (LIP) technique, McKay and Borucki(l) found that organic synthesis preferentially occurred in a reducing gas mixture rich in methane, and not in a mixture rich in carbon dioxide. This result means chemical models based on the thermodynamical equilibrium approach do not apply to shock chemistry. In this study, we employ the technique of reacting flow, i.e., chemical kinetics in a fluid flow, to simulate the chemistry occurring in LIP and in the wake region from comet or meteor impact. Three different air compositions have been used: (1) 1/3 CO2 and 2/3 H2, (2) pure CH4, and (3) 1/4 CH4, 1/4 CO2, and 1/2 H2O. The stoichiometric ratio of gas mixtures (1) and (3) are kept the same. For (1) we obtain equal mole fractions of CO and H2O as the major products and for (2) C2H2 is the major product. In both cases our results are in agreement with Ref. (1). For (3) we find an interesting case where the nature of chemicals produced to be critically dependent on the flowfield temperature. At the higher temperature part of the wake region, CO and H2O are the dominant products, whereas in the cooler region C2H2 is the dominant product. Further studies of these reactions, as well as for the gas mixture including N2, are being pursued.

Dateo, Christopher E.↗

Simultaneous Laser Raman-rayleigh-lif Measurements and Numerical Modeling Results of a Lifted Turbulent H2/N2 Jet Flame in a Vitiated Coflow

An experiment and numerical investigation is presented of a lifted turbulent H2/N2 jet flame in a coflow of hot, vitiated gases. The vitiated coflow burner emulates the coupling of turbulent mixing and chemical kinetics exemplary of the reacting flow in the recirculation region of advanced combustors. It also simplifies numerical investigation of this coupled problem by removing the complexity of recirculating flow. Scalar measurements are reported for a lifted turbulent jet flame of H2/N2 (Re = 23,600, H/d = 10) in a coflow of hot combustion products from a lean H2/Air flame ((empty set) = 0.25, T = 1,045 K). The combination of Rayleigh scattering, Raman scattering, and laser-induced fluorescence is used to obtain simultaneous measurements of temperature and concentrations of the major species, OH, and NO. The data attest to the success of the experimental design in providing a uniform vitiated coflow throughout the entire test region. Two combustion models (PDF: joint scalar Probability Density Function and EDC: Eddy Dissipation Concept) are used in conjunction with various turbulence models to predict the lift-off height (H(sub PDF)/d = 7,H(sub EDC)/d = 8.5). Kalghatgi's classic phenomenological theory, which is based on scaling arguments, yields a reasonably accurate prediction (H(sub K)/d = 11.4) of the lift-off height for the present flame. The vitiated coflow admits the possibility of auto-ignition of mixed fluid, and the success of the present parabolic implementation of the PDF model in predicting a stable lifted flame is attributable to such ignition. The measurements indicate a thickened turbulent reaction zone at the flame base. Experimental results and numerical investigations support the plausibility of turbulent premixed flame propagation by small scale (on the order of the flame thickness) recirculation and mixing of hot products into reactants and subsequent rapid ignition of the mixture.

Cabra, R.↗