Search

There is considerable technological interest in developing new fuel-flexible combustion systems that can burn fuels such as hydrogenor syngas. Lean premixed systems have the potential to burn these types of fuels with high efficiency and low NOx emissions due to reduced burnt gas temperatures. Although traditional scientific approaches based on theory and laboratory experiment have played essential roles in developing our current understanding of premixed combustion, they are unable to meet the challenges of designing fuel-flexible lean premixed combustion devices. Computation, with itsability to deal with complexity and its unlimited access to data, hasthe potential for addressing these challenges. Realizing this potential requires the ability to perform high fidelity simulations of turbulent lean premixed flames under realistic conditions. In this paper, we examine the specialized mathematical structure of these combustion problems and discuss simulation approaches that exploit this structure. Using these ideas we can dramatically reduce computational cost, making it possible to perform high-fidelity simulations of realistic flames. We illustrate this methodology by considering ultra-lean hydrogen flames and discuss how this type of simulation is changing the way researchers study combustion.

This paper introduces an adaptive mesh and algorithm refinement method for fluctuating hydrodynamics. This particle-continuum hybrid simulates the dynamics of a compressible fluid with thermal fluctuations. The particle algorithm is direct simulation Monte Carlo (DSMC), a molecular-level scheme based on the Boltzmann equation. The continuum algorithm is based on the Landau-Lifshitz Navier-Stokes (LLNS) equations, which incorporate thermal fluctuations into macroscopic hydrodynamics by using stochastic fluxes. It uses a recently-developed solver for LLNS, based on third-order Runge-Kutta. We present numerical tests of systems in and out of equilibrium, including time-dependent systems, and demonstrate dynamic adaptive refinement by the computation of a moving shock wave. Mean system behavior and second moment statistics of our simulations match theoretical values and benchmarks well. We find that particular attention should be paid to the spectrum of the flux at the interface between the particle and continuum methods, specifically for the non-hydrodynamic (kinetic) time scales.

Using examples from active research areas in combustion and astrophysics, we demonstrate a computationally efficient numerical approach for simulating multiscale low Mach number reacting flows. The method enables simulations that incorporate an unprecedented range of temporal and spatial scales, while at the same time, allows an extremely high degree of reaction fidelity. Sample applications demonstrate the efficiency of the approach with respect to a traditional time-explicit integration method, and the utility of the methodology for studying the interaction of turbulence with terrestrial and astrophysical flame structures.

Type Ia supernovae (SNe Ia) are the largest thermonuclear explosions in the Universe. Their light output can be seen across great stances and has led to the discovery that the expansion rate of the Universe is accelerating. Despite the significance of SNe Ia, there are still a large number of uncertainties in current theoretical models. Computational modeling offers the promise to help answer the outstanding questions. However, even with today's supercomputers, such calculations are extremely challenging because of the wide range of length and time scales. In this paper, we discuss several new algorithms for simulations of SNe Ia and demonstrate some of their successes.

We present three-dimensional, time-dependent simulations of the flowfield of a laboratory-scale slot burner. The simulations are performed using an adaptive time-dependent low Mach number combustion algorithm based on a second-order projection formulation that conserves both species mass and total enthalpy. The methodology incorporates detailed chemical kinetics and a mixture model for differential species diffusion. Methane chemistry and transport are modeled using the DRM-19 mechanism along with its associated thermodynamics and transport databases. Adaptive mesh refinementdynamically resolves the flame and turbulent structures. Detailedcomparisons with experimental measurements show that the computational results provide a good prediction of the flame height, the shape of the time-averaged parabolic flame surface area, and the global consumption speed (the volume per second of reactants consumed divided by the area of the time-averaged flame). The thickness of the computed flamebrush increases in the streamwise direction, and the flamesurface density profiles display the same general shapes as the experiment. The structure of the simulated flame also matches the experiment; reaction layers are thin (typically thinner than 1 mm) and the wavelengths of large wrinkles are 5--10 mm. Wrinkles amplify to become long fingers of reactants which burn through at a neck region, forming isolated pockets of reactants. Thus both the simulated flame and the experiment are in the "corrugated flamelet regime."

The convective period leading up to a Type Ia supernova (SN Ia) explosion is characterized by very low Mach number flows, requiring hydrodynamical methods well-suited to long-time integration. We continue the development of the low Mach number equation set for stellar scale flows by incorporating the effects of heat release due to external sources. Low Mach number hydrodynamics equations with a time-dependent background state are derived, and a numerical method based on the approximate projection formalism is presented. We demonstrate through validation with a fully compressible hydrodynamics code that this low Mach number model accurately captures the expansion of the stellar atmosphere as well as the local dynamics due to external heat sources. This algorithm provides the basis for an efficient simulation tool for studying the ignition of SNe Ia.

Recent progress in simulation methodologies and high-performance parallel computers have made it is possible to perform detailed simulations of multidimensional reacting flow phenomena using comprehensive kinetics mechanisms. As simulations become larger and more complex, it becomes increasingly difficult to extract useful information from the numerical solution, particularly regarding the interactions of the chemical reaction and diffusion processes. In this paper we present a new diagnostic tool for analysis of numerical simulations of reacting flow. Our approach is based on recasting an Eulerian flow solution in a Lagrangian frame. Unlike a conventional Lagrangian view point that follows the evolution of a volume of the fluid, we instead follow specific chemical elements, e.g., carbon, nitrogen, etc., as they move through the system . From this perspective an "atom" is part of some molecule of a species that is transported through the domain by advection and diffusion. Reactions cause the atom to shift from one chemical host species to another and the subsequent transport of the atom is given by the movement of the new species. We represent these processes using a stochastic particle formulation that treats advection deterministically and models diffusion and chemistry as stochastic processes. In this paper, we discuss the numerical issues in detail and demonstrate that an ensemble of stochastic trajectories can accurately capture key features of the continuum solution. The capabilities of this diagnostic are then demonstrated by applications to study the modulation of carbon chemistry during a vortex-flame interaction, and the role of cyano chemistry in rm NO_x production for a steady diffusion flame.

Many turbulent premixed flames of practical interest are statistically stationary. They occur in combustors that have anchoring mechanisms to prevent blow-off and flashback. The stabilization devices often introduce a level of geometric complexity that is prohibitive for detailed computational studies of turbulent flame dynamics. As a result, typical detailed simulations are performed in simplified model configurations such as decaying isotropic turbulence or inflowing turbulence. In these configurations, the turbulence seen by the flame either decays or, in the latter case, increases as the flame accelerates toward the turbulent inflow. This limits the duration of the eddy evolutions experienced by the flame at a given level of turbulent intensity, so that statistically valid observations cannot be made. In this paper, we apply a feedback control to computationally stabilize an otherwise unstable turbulent premixed flame in two dimensions. For the simulations, we specify turbulent inflow conditions and dynamically adjust the integrated fueling rate to control the mean location of the flame in the domain. We outline the numerical procedure, and illustrate the behavior of the control algorithm. We use the simulations to study the propagation and the local chemical variability of turbulent flame chemistry.

We present three-dimensional, time-dependent simulations of a full-size laboratory-scale rod-stabilized premixed turbulent V-flame. The computations use an adaptive projection method based on a low Mach number formulation that incorporates detailed chemical kinetics and transport. The simulations are performed without introducing models for turbulence or turbulence chemistry interaction. We outline the numerical procedure and experimental setup, and compare computed results to mean flame location and surface wrinkling statistics gathered from experiment.

In this paper we discuss the application of a new diagnostic tool for analysis of flame simulations. This methodologogy is based on following specific chemical elements, e.g., carbon or nitrogen, as they move through the system. From this perspective an "atom" is a component of a molecule that is being transported through the simulation domain by advection and diffusion. Reactions cause the atom to shift from one species to another with the subsequent transport of the atom determined by the movement of the new species.

In this paper we use numerical simulation to investigate shock-induced ignition and combustion of a hydrocarbon gas. The focus of this paper is on quantifying the effect of fidelity in the chemical kinetics on the overall solution. We model the system using the compressible Navier Stokes equations for a reacting mixture. These equations express conservation of species mass, momentum, total energy.

With adaptive-grid computational methodologies and judicious use of compressible and low Mach number combustion models, we are carrying out three-dimensional, time-dependent direct numerical simulations of a laboratory-scale turbulent premixed methane burner. In the laboratory experiment, turbulence is generated by a grid located in the throat of a 50mm diameter circular nozzle; swirl is be introduced by four tangential air jets spaced uniformly around the circumference of the nozzle just above the turbulence grid. A premixed methane flame is stabilized above the nozzle in the central core region where a velocity deficit is induced by the swirling flow. The time-dependent flow field inside the nozzle, from the turbulence grid and the high-speed jets, to the nozzle exit plane is simulated using an adaptive-grid embedded-boundary compressible Navier-Stokes solver. The compressible calculation then provides time-dependent boundary conditions for an adaptive low Mach number model of the swirl-stabilized premixed flame. The low Mach model incorporates detailed chemical kinetics and species transport using 20 species and 84 reactions. Laboratory diagnostics available for comparisons include characterizations of the flow field just down stream of the nozzle exit plane, and flame surface statistics, such as mean location, wrinkling and crossing frequencies.

This conference paper considers how to use reaction path diagrams to better understand the output of reacting flow simulations. Briefly, these diagrams have long been used to depict the reactants and products in networks of chemical reactions. The diagrams can be generated in several ways from computer simulations of chemically reacting fluids to depict how the fluid moderates the chemistry by determining which species are brought into contact to react in quantity. The concept of a conditional diagram is introduced which depicts the reactions occurring in only a portion of the fluid domain, thus enabling comparisons between different regions of the fluid and the overall reaction network. Several examples are provided of the paths occurring in methane diffusion flames.

This paper reports on a combined experimental and modeling investigation of NOx formation in nitrogen-diluted laminar methane diffusion flames seeded with ammonia. The methane-ammonia mixture is a surrogate for biomass fuels which contain significant fuel-bound nitrogen. The experiments use flue-gas sampling to measure the concentration of stable species in the exhaust gas, including NO, O2, CO, and CO2. The computations evolve a two-dimensional low Mach number model using a solution-adaptive projection algorithm to capture fine-scale features of the flame. The model includes detailed thermodynamics and chemical kinetics, differential diffusion, buoyancy, and radiative losses. The model shows good agreement with the measurements over the full range of experimental NH3 seeding amounts. As more NH3 is added, a greater percentage is converted to N2 rather than to NO. The simulation results are further analyzed to trace the changes in NO formation mechanisms with increasing amounts of ammonia in the fuel.

In this paper we study the behavior of a premixed turbulent methane flame in three dimensions using numerical simulation. The simulations are performed using an adaptive time-dependent low Mach number combustion algorithm based on a second-order projection formulation that conserves both species mass and total enthalpy. The species and enthalpy equations are treated using an operator-split approach that incorporates stiff integration techniques for modeling detailed chemical kinetics. The methodology also incorporates a mixture model for differential diffusion. For the simulations presented here, methane chemistry and transport are modeled using the DRM-19 (19-species, 84-reaction) mechanism derived from the GRIMech-1.2 mechanism along with its associated thermodynamics and transport databases. We consider a lean flame with equivalence ratio 0.8 for two different levels of turbulent intensity. For each case we examine the basic structure of the flame including turbulent flame speed and flame surface area. The results indicate that flame wrinkling is the dominant factor leading to the increased turbulent flame speed. Joint probability distributions are computed to establish a correlation between heat release and curvature. We also investigate the effect of turbulent flame interaction on the flame chemistry. We identify specific flame intermediates that are sensitive to turbulence and explore various correlations between these species and local flame curvature. We identify different mechanisms by which turbulence modulates the chemisry of the flame.

In this paper we study the formation of NO in laminar, nitrogen diluted methane diffusion flames that are seeded with ammonia in the fuel stream. We have performed numerical simulations with detailed chemistry as well as laser-induced fluorescence imaging measurements for a range of ammonia injection rates. For comparison with the experimental data, synthetic LIF images are calculated based on the numerical data accounting for temperature and fluorescence quenching effects. We demonstrate good agreement between measurements and computations. The LIF corrections inferred from the simulation are then used to calculate absolute NO mole fractions from the measured signal.The NO formation in both doped and undoped flames occurs in the flame sheet. In the undoped flame, four different mechanisms including thermal and prompt NO appear to contribute to NO formation. As the NH3 seeding level increases, fuel-NO becomes the dominant mechanism and N2 shifts from being a net reactant to being a net product. Nitric oxide in the undoped flame as well as in the core region of the doped flames are underpredicted by the model; we attribute this mainly to inaccuracies in the NO recycling chemistry on the fuel-rich side of the flame sheet.

Recent progress in simulation methodologies and new, high-performance parallel architectures have made it is possible to perform detailed simulations of multidimensional combustion phenomena using comprehensive kinetics mechanisms. However, as simulation complexity increases, it becomes increasingly difficult to extract detailed quantitative information about the flame from the numerical solution, particularly regarding the details of chemical processes. In this paper we present a new diagnostic tool for analysis of numerical simulations of combustion phenomena. Our approach is based on recasting an Eulerian flow solution in a Lagrangian frame. Unlike a conventional Lagrangian viewpoint in which we follow the evolution of a volume of the fluid, we instead follow specific chemical elements, e.g., carbon, nitrogen, etc., as they move through the system. From this perspective an "atom" is part of some molecule that is transported through the domain by advection and diffusion. Reactions cause the atom to shift from one species to another with the subsequent transport given by the movement of the new species. We represent these processes using a stochastic particle formulation that treats advection deterministically and models diffusion as a suitable random-walk process. Within this probabilistic framework, reactions can be viewed as a Markov process transforming molecule to molecule with given probabilities. In this paper, we discuss the numerical issues in more detail and demonstrate that an ensemble of stochastic trajectories can accurately capture key features of the continuum solution. We also illustrate how the method can be applied to studying the role of cyanochemistry on NOx production in a diffusion flame.

The dominant computational cost in modeling turbulent combustion phenomena numerically with high fidelity chemical mechanisms is the time required to solve the ordinary differential equations associated with chemical kinetics. One approach to reducing that computational cost is to develop an inexpensive surrogate model that accurately represents evolution of chemical kinetics. One such approach, PRISM, develops a polynomial representation of the chemistry evolution in a local region of chemical composition space. This representation is then stored for later use. As the computation proceeds, the chemistry evolution for other points within the same region are computed by evaluating these polynomials instead of calling an ordinary differential equation solver. If initial data for advancing the chemistry is encountered that is not in any region for which a polynomial is defined, the methodology dynamically samples that region and constructs a new representation for that region. The utility of this approach is determined by the size of the regions over which the representation provides a good approximation to the kinetics and the number of these regions that are necessary to model the subset of composition space that is active during a simulation. In this paper, we assess the PRISM methodology in the context of a turbulent premixed flame in two dimensions. We consider a range of turbulent intensities ranging from weak turbulence that has little effect on the flame to strong turbulence that tears pockets of burning fluid from the main flame. For each case, we explore a range of sizes for the local regions and determine the scaling behavior as a function of region size and turbulent intensity.

The interaction of a vortex pair with a premixed flame serves as an important prototype for premixed turbulent combustion. In this study, we investigate the interaction of a counter-rotating vortex pair with an initially flat premixed methane flame. We focus on characterizing the mechanical nature of the flame-vortex interaction and on the features of the interaction strongly affected by fuel equivalence ratio. We compare computational solutions obtained using a time-dependent, two-dimensional adaptive low Mach number combustion algorithm that incorporates GRI-Mech 1.2 for the chemistry, thermodynamics and transport of the chemical species. We find that the circulation around the vortex scours gas from the preheat zone in front of the flame, making the interaction extremely sensitive to equivalence ratio. For nearly stoichiometric cases, the peak mole fraction of CH across the flame is relatively insensitive to the vortex whereas for richer flames we observe a substantial and rapid decline in the peak CH mole fraction, commencing early in the flame-vortex interaction. The peak concentration of HCO is found to correlate, in both space and time, with the peak heat release across a broad range of equivalence ratios. The model also predicts a measurable increase in C2H2 as a result of interaction with the vortex, and a marked increase in the low temperature chemistry activity.