Search NASA⌕ Search

SEARCH · Search NASA

Results for “Non-Equilibrium Flow”

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 145 records · Page 8

Electrified vapour deposition at ultrahigh temperature and atmospheric pressure for nanomaterials synthesis

Vapour-phase synthesis methods have shown promise for the scalable synthesis of nanomaterials and coatings. However, the vaporization of different precursors for the synthesis of a broad nanomaterial space, particularly at atmospheric pressure, while maintaining compositional and structural control of the final product is challenging. Here we report the generation of an ultrahigh-temperature atomic vapour at atmospheric pressure based on electrified heating, for the growth of multi-elemental nanomaterials and thin films. This process relies on a reactor design whereby solid-state precursors are vaporized within a semi-confined space beneath an electrified heater that can reach ~3,000 K. The proximity of the heater rapidly breaks down the bonds of metal salt precursors and decomposes them into an atomic vapour that expands into a high-temperature (>2,000 K), highly reactive and high-flux vapour (10 21 –10 22 atoms per cm 2 per second) that travels upwards in a directional flow. When mixed with entrained ambient gases, the highly reactive atomic species rapidly nucleate and grow into the desired final products, including alloys, oxides, sulfides and thin films, which can be deposited on a low-temperature substrate. This EVD approach can synthesize a broad range of functional nanomaterials at atmospheric pressure, including single-phase multi-elemental nanomaterials formed under thermodynamically non-equilibrium conditions.

36 MATERIALS SCIENCE↗

Transition in Hypersonic Flows Including High-temperature Gas Effects

Hypersonic transition poses a special challenge for direct numerical simulations. Comparable data from Wind-tunnel tests or free-flight testing are not available or not accurate enough for comparison. The wind-tunnel testing does not allow for the exact match to the free-flight conditions at such high Mach-numbers. Flat-plate boundary-layer transition at high Mach-numbers is investigated in this work. A simulation case was chosen where chemical non-equilibrium plays an important role but ionization can be neglected. The chosen case at an altitude of H=50Km lies close to one point on the descent path of the Space Shuttle. The failure of the Space Shuttle has shown that an improved vehicle for space transportation is imperative in the close future. Transition research for an improved space-transportation vehicle is crucial in order to estimate the heat load during re-entry.

Stemmer, Christian↗

A Technique for Determining the Nozzle-Flow Properties of Air in an Equilibrium, Nonequilibrium, or Frozen State

One of the problems associated with the development of high-enthalpy hypersonic test facilities hos been the determination of nozzle flow properties. This note presents a technique for determining the flow properties of air in an equilibrium, non-equilibrium, or frozen state using two test section measurements. A knowledge of two stagnation properties is required before the technique can be applied. The stagnation pressure can be measured, and reference 1 offers a method for determining stagnation enthalpy. Results of the calculations using the method of reference 1 are presented in reference 2 in a chart (chart 21) which can be readily used to determine the stagnation enthalpy from measurements of stagnation pressure, mass flow of air through the tunnel, and throat diameter.

Hayman, Lovick O.↗

RTE: A computer code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Naraghi, Mohammad H. N.↗

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Emerging supersolidity in photonic-crystal polariton condensates

A supersolid is a counter-intuitive phase of matter in which its constituent particles are arranged into a crystalline structure, yet they are free to flow without friction. This requires the particles to share a global macroscopic phase while being able to reduce their total energy by spontaneous, spatial self-organization. The existence of the supersolid phase of matter was speculated more than 50 years ago1, 2, 3–4. However, only recently has there been convincing experimental evidence, mainly using ultracold atomic Bose–Einstein condensates (BECs) coupled to electromagnetic fields. There, various guises of the supersolid were created using atoms coupled to high-finesse cavities5,6, with large magnetic dipole moments7, 8, 9, 10, 11, 12–13, and spin–orbit-coupled, two-component systems showing stripe phases14, 15–16. Here we provide experimental evidence of a new implementation of the supersolid phase in a driven-dissipative, non-equilibrium context based on exciton–polaritons condensed in a topologically non-trivial, bound state in the continuum (BiC) with exceptionally low losses, realized in a photonic-crystal waveguide. We measure the density modulation of the polaritonic state indicating the breaking of translational symmetry with a precision of several parts in a thousand. Direct access to the phase of the wavefunction allows us to also measure the local coherence of the supersolid. We demonstrate the potential of our synthetic photonic material to host phonon dynamics and a multimode excitation spectrum.

Trypogeorgos, Dimitrios↗

Experimental and analytical study of cryogenic propellant boiloff to develop and verify alternate pressurization concepts for Space Shuttle external tank using a scaled down tank

Self pressurization by propellant boiloff is experimentally studied as an alternate pressurization concept for the Space Shuttle external tank (ET). The experimental setup used in the study is an open flow system which is composed of a variable area test tank and a recovery tank. The vacuum jacketed test tank is geometrically similar to the external LOx tank for the Space Shuttle. It is equipped with instrumentation to measure the temperature and pressure histories within the liquid and vapor, and viewports to accommodate visual observations and Laser-Doppler Anemometry measurements of fluid velocities. A set of experiments were conducted using liquid Nitrogen to determine the temperature stratification in the liquid and vapor, and pressure histories of the vapor during sudden and continuous depressurization for various different boundary and initial conditions. The study also includes the development and calibration of a computer model to simulate the experiments. This model is a one-dimensional, multi-node type which assumes the liquid and the vapor to be under non-equilibrium conditions during the depressurization. It has been tested for a limited number of cases. The preliminary results indicate that the accuracy of the simulations is determined by the accuracy of the heat transfer coefficients for the vapor and the liquid at the interface which are taken to be the calibration parameters in the present model.

Akyuzlu, K. M.↗

H2OTSTUF: Appropriate Operating Regimes for Magnetohydrodynamic Augmentation

A trade study of magnetohydrodynamic (MHD) augmented propulsion reveals a unique operating regime at lower thrust levels. Substantial mass savings are realized over conventional chemical, solar, and electrical propulsion concepts when MHD augmentation is used to obtain optimal I(sub sp). However, trip times for the most conservative estimates of power plant specific impulse and accelerator efficiency may be prohibitively long. Quasi-one-dimensional calculations show that a solar or nuclear thermal system augmented by MHD can provide competitive performance while utilizing a diverse range of propellants including water, which is available from the Space Shuttle, the Moon, asteroids, and various moons and planets within our solar system. The use of in-situ propellants will reduce costs of space operations as well as enable human exploration of our Solar System. The following conclusions can be drawn from the results of the mission trade study: (1) There exists a maximum thrust or mass flow rate above which MHD augmentation increases the initial mass in low earth orbit (LEO); (2) Mass saving of over 50% can be realized for unique combination of solar/MHD systems; (3) Trip times for systems utilizing current power supply technology may be prohibitively long. Theoretical predictions of MHD performance for in space propulsion systems show that improved efficiencies can reduce trip times to acceptable levels; (4) Long trip times indicative of low thrust systems can be shortened by an increase in the MHD accelerator efficiency or a decrease in the specific mass of the power supply and power processing unit; and (5) As for all propulsion concepts, missions with larger (Delta)v's benefit more from the increased specific impulse resulting from MHD augmentation. Using a quasi-one-dimensional analysis, the required operating conditions for a MHD accelerator to reach acceptable efficiencies are outlined. This analysis shows that substantial non-equilibrium ionization is desirable.

Jones, Jonathan E.↗

High-temperature LDV seed particle development

The feasibility of developing a method for making monodisperse, unagglomerated spherical particles greater than 50 nm in diameter was demonstrated. Carbonaceous particles were made by pyrolyzing ethylene with a pulsed CO2 laser, thereby creating a non-equilibrium mixture of carbon, hydrogen, hydrocarbon vapors, and unpyrolyzed ethylene. Via a complex series of reactions, the carbon and hydrocarbon vapors quickly condensed into the spherical particles. By cooling and dispersing them in a supersonic expansion immediately after their creation, the hot newly-formed spheres were prevented from colliding and coalescing, thus preventing the problem of agglomeration which as plagued other investigators studying laser-simulated particle formation. The cold particles could be left suspended in the residual gases indefinitely without agglomerating. Their uniform sizes and unagglomerated nature were visualized by collecting the particles on filters that were subsequently examined using electron microscopy. It was found the mean particle size can be coarsely controlled by varying the initial ethylene pressure, and can be finely controlled by varying the fluence (energy/unit area) with which the laser irradiates the gas. The motivating application for this research was to manufacture particles that could be used as laser Doppler velocimetry (LDV) seeds in high-temperature high-speed flows. Though the particles made in this program will not evaporate until heated to about 3000 K, and thus could serve as LDV seeds in some applications, they are not ideal when the hot atmosphere is also oxidizing. In that situation, ceramic materials would be preferable. Research performed elsewhere has demonstrated that selected ceramic materials can be manufactured by laser pyrolysis of appropriate supply gases. It is anticipated that, when the same gases are used in conjunction with the rapid cooling technique, unagglomerated spherical ceramic particles can be made with little difficulty. Such particles would also be valuable to manufacturers of ceramic or abrasive products, and this technique may find its greatest commercial potential in those areas.

Frish, Michael B.↗

The nonlinear evolution of modes on unstable stratified shear layers

The nonlinear development of disturbances in stratified shear flows (having a local Richardson number of value less than one quarter) is considered. Such modes are initially fast growing but, like related studies, we assume that the viscous, non-parallel spreading of the shear layer results in them evolving in a linear fashion until they reach a position where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. Four different basic integro-differential amplitude equations are possible, including one due to a novel mechanism; the relevant choice of amplitude equation, at a particular instance, being dependent on the relative sizes of the disturbance amplitude, the growth rate of the disturbance, its wavenumber, and the viscosity of the fluid. This richness of choice of possible nonlinearities arises mathematically from the indicial Frobenius roots of the governing linear inviscid equation (the Taylor-Goldstein equation) not, in general, differing by an integer. The initial nonlinear evolution of a mode will be governed by an integro-differential amplitude equations with a cubic nonlinearity but the resulting significant increase in the size of the disturbance's amplitude leads on to the next stage of the evolution process where the evolution of the mode is governed by an integro-differential amplitude equations with a quintic nonlinearity. Continued growth of the disturbance amplitude is expected during this stage, resulting in the effects of nonlinearity spreading to outside the critical level, by which time the flow has become fully nonlinear.

Blackaby, Nicholas↗

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗

Supersonic flow of chemically reacting gas-particle mixtures. Volume 1: A theoretical analysis and development of the numerical solution

A numerical solution for chemically reacting supersonic gas-particle flows in rocket nozzles and exhaust plumes was described. The gas-particle flow solution is fully coupled in that the effects of particle drag and heat transfer between the gas and particle phases are treated. Gas and particles exchange momentum via the drag exerted on the gas by the particles. Energy is exchanged between the phases via heat transfer (convection and/or radiation). Thermochemistry calculations (chemical equilibrium, frozen or chemical kinetics) were shown to be uncoupled from the flow solution and, as such, can be solved separately. The solution to the set of governing equations is obtained by utilizing the method of characteristics. The equations cast in characteristic form are shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The particle distribution is represented in the numerical solution by a finite distribution of particle sizes.

Penny, M. M.↗

A survey of upwind methods for flows with equilibrium and non-equilibrium chemistry and thermodynamics

Several versions of flux-vector split and flux-difference split algorithms were compared with regard to general applicability and complexity. Test computations were performed using curve-fit equilibrium air chemistry for an M = 5 high-temperature inviscid flow over a wedge, and an M = 24.5 inviscid flow over a blunt cylinder for test computations; for these cases, little difference in accuracy was found among the versions of the same flux-split algorithm. For flows with nonequilibrium chemistry, the effects of the thermodynamic model on the development of flux-vector split and flux-difference split algorithms were investigated using an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Several numerical examples are presented, including nonequilibrium air chemistry in a high-temperature shock tube and nonequilibrium hydrogen-air chemistry in a supersonic diffuser.

Grossman, B.↗

Upwind schemes and bifurcating solutions in real gas computations

The area of high speed flow is seeing a renewed interest due to advanced propulsion concepts such as the National Aerospace Plane (NASP), Space Shuttle, and future civil transport concepts. Upwind schemes to solve such flows have become increasingly popular in the last decade due to their excellent shock capturing properties. In the first part of this paper the authors present the extension of the Osher scheme to equilibrium and non-equilibrium gases. For simplicity, the source terms are treated explicitly. Computations based on the above scheme are presented to demonstrate the feasibility, accuracy and efficiency of the proposed scheme. One of the test problems is a Chapman-Jouguet detonation problem for which numerical solutions have been known to bifurcate into spurious weak detonation solutions on coarse grids. Results indicate that the numerical solution obtained depends both on the upwinding scheme used and the limiter employed to obtain second order accuracy. For example, the Osher scheme gives the correct CJ solution when the super-bee limiter is used, but gives the spurious solution when the Van Leer limiter is used. With the Roe scheme the spurious solution is obtained for all limiters.

Suresh, Ambady↗

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.↗

MHD Energy Bypass Scramjet Engine

Revolutionary rather than evolutionary changes in propulsion systems are most likely to decrease cost of space transportation and to provide a global range capability. Hypersonic air-breathing propulsion is a revolutionary propulsion system. The performance of scramjet engines can be improved by the AJAX energy management concept. A magneto-hydro-dynamics (MHD) generator controls the flow and extracts flow energy in the engine inlet and a MHD accelerator downstream of the combustor accelerates the nozzle flow. A progress report toward developing the MHD technology is presented herein. Recent theoretical efforts are reviewed and ongoing experimental efforts are discussed. The latter efforts also include an ongoing collaboration between NASA, the US Air Force Research Laboratory, US industry, and Russian scientific organizations. Two of the critical technologies, the ionization of the air and the MHD accelerator, are briefly discussed. Examples of limiting the combustor entrance Mach number to a low supersonic value with a MHD energy bypass scheme are presented, demonstrating an improvement in scramjet performance. The results for a simplified design of an aerospace plane show that the specific impulse of the MHD-bypass system is better than the non-MHD system and typical rocket over a narrow region of flight speeds and design parameters. Equilibrium ionization and non-equilibrium ionization are discussed. The thermodynamic condition of air at the entrance of the engine inlet determines the method of ionization. The required external power for non-equilibrium ionization is computed. There have been many experiments in which electrical power generation has successfully been achieved by magneto-hydrodynamic (MHD) means. However, relatively few experiments have been made to date for the reverse case of achieving gas acceleration by the MHD means. An experiment in a shock tunnel is described in which MHD acceleration is investigated experimentally. MHD has several potential aerospace applications. The first is to improve the performance of hypersonic air-breathing engines for space launch and cruise vehicles. The second is to improve the performance of a high enthalpy wind tunnel. The third is to control a hypersonic vehicle. With such applications in mind, theoretical and experiments are being conducted at the NASA Ames Research Center to develop the MHD technology.

Mehta, Unmeel B.↗

Quantitative Species Measurements in Microgravity Combustion Flames

Flame-vortex interactions are canonical configurations that can be used to study the underlying processes occurring in complicated turbulent reacting flows. The elegant simplicity of the flame-vortex interaction permits the study of these complex interactions under relatively controllable experimental configurations, in contrast to direct measurements in turbulent flames. The ability to measure and model the fundamental phenomena that occur in a turbulent flame, but with time and spatial scales which are amenable to our diagnostics, permits significant improvements in the understanding of turbulent combustion under both normal and reduced gravity conditions. In this paper, we report absolute mole fraction measurements of methane in a reacting vortex ring. These microgravity experiments are performed in the 2.2-sec drop tower at NASA Glenn Research Center. In collaboration with Drs. Chen and Dahm at the University of Michigan, measured methane absorbances are incorporated into a new model from which the temperature and concentrations of all major gases in the flame can be determined at all positions and times in the development of the vortex ring. This is the first demonstration of the ITAC (Iterative Temperature with Assumed Chemistry) approach, and the results of these computations and analyses are presented in a companion paper by Dahm and Chen at this Workshop. We believe that the ITAC approach will become a powerful tool in understanding a wide variety of combustion flames under both equilibrium and non-equilibrium conditions.

Silver, Joel A.↗