Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulent flames”

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 109 records · Page 6

100 kHz OH PLIF of High-Pressure H2/NH3 and H2/CH4 Flames

The reaction zones in a multi-element, multi-stage, micromix injector are investigated using planar laser-induced fluorescence (PLIF) imaging of the hydroxyl (OH) radical at a 100 kHz repetition rate. The combustor is operated at a pressure of ∼10 bar. Multiple fuel blends of hydrogen, methane, and ammonia are considered to study the effect of changing fuel composition on the flame structure. Compact regions of high OH signal are observed in each flame for fuel mixtures with high hydrogen content, indicating rapid reactant consumption in the shear layers. Instantaneous images exhibit wrinkled flame fronts and the formation of unreacted pockets due to highly turbulent burning. The addition of methane is observed to reduce the reactivity and flame speed, resulting in more diffuse regions of OH between each flame. The shape of the reactant core is altered, and the length scale of the flame surface wrinkling is increased. Ammonia addition has a more severe effect on the distribution of OH and the local flame structure than methane addition. A time series of images shows the presence of OH vortical structures and significant turbulent fluctuations in the flames near the wall, suggesting the presence of a recirculation zone that could drive the observed disparity in flame lengths.

Hodge, Alexander J. [Purdue University]↗

Simulations of Multi-Mode Combustion Regimes Realizable In a Gasoline Direct Injection Engine

Lean and dilute gasoline compression ignition (GCI) operation in spark ignition (SI) engines are an attractive strategy to attain high fuel efficiency and low NOx levels. However, this combustion mode is often limited to low-load engine conditions due to the challenges associated with autoignition controllability. In order to overcome this constrain, multi-mode (MM) operating strategies, consisting of advanced compression ignition (ACI) at low load and conventional SI at high load, have been proposed. In this three-dimensional computational fluid dynamics study, the concept of multi-mode combustion using two RON98 gasoline fuel blends (Co-Optima Alkylate and E30) in a gasoline direct injection (GDI) engine were explored. To this end, a new reduced mechanism for simulating the kinetics of E30 fuel blend is introduced in this study. To cover the varying engine load demands for multi-mode engines, primary combustion dynamics observed in ACI and SI combustion modes was characterized and validated against experimental measurements. In order to implement part-load conditions, a strategy of mode transition between SI and ACI combustion (i.e., mixed-mode combustion) was then explored numerically by creating a virtual test condition. The results obtained from the mixed-mode simulations highlight an important feature that deflagrative flame propagation regime coexists with ignition-assisted end-gas autoignition. This study also identifies a role of turbulent flow property adjacent to premixed flame front in characterizing the mixed-mode combustion. The employed hybrid combustion model was verified to perform simulations aiming at suitable range of multi-mode engine operations.

Kim, Sayop↗

MACHINE LEARNING BASED CHEMICAL EXPLOSIVE MODE ANALYSIS (ML-CEMA)

The software consists of a Machine Learning based Chemical Explosive Mode Analysis (ML-CEMA) tool for advanced computational flame diagnostics. CEMA, originally based on eigen-analysis of the local thermochemical system, is capable of identifying reaction fronts and limit phenomena such as auto-ignition and extinction in practical combustion systems, such as internal combustion engines and gas turbine combustors. However, the original CEMA is computationally expensive for large reaction mechanisms that are typically needed to describe fuel chemistry of practical large-hydrocarbon fuels. This novel ML-CEMA tool employs a ML technique to accelerate the eigen-analysis of the basic CEMA approach by orders of magnitude, thus making it suitable for practical fuels. In ML-CEMA,zero-dimensional (0D) reactors and one-dimensional (1D) premixed flames are first used to generate a large number of data points for neural network based ML training. The trained ML model is then used to perform CEMA prediction. This ML-CEMA tool has been demonstrated in canonical 0D and 1D configurations as well as highly-transient three-dimensional spray flames exhibiting multi-mode turbulent combustion, showing promising results. ML-CEMA, as a standalone tool, can be used for computationally-efficient diagnostics of massive datasets generated from both experiments and simulations. For example, based on spatially resolved measurements of a small set of reactive scalars(such as temperature, hydroxyl radical and formaldehyde), ML-CEMA can effectively identify flame fronts and rare events. ML-CEMA also provides a robust online or offline flame feature detection tool. When used for on-the-fly simulations, ML-CEMA further enables zone-adaptive combustion modeling, in which the predicted eigenvalue is used as a robust mode indicator for judicious assignment of locally-valid combustion models. This ML-CEMA based zone-adaptive model can lead to substantial computational cost savings when used for large-scale simulations of multi-mode combustion systems. Third Party Code Web Page to Download Code Web Page Location of Third Party License

Xu, Chao↗

Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine

The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.

42 ENGINEERING↗

A direct numerical simulation study of the dilution tolerance of propane combustion under spark-ignition engine conditions

Modern spark ignition internal combustion (IC) engines rely on highly diluted fuel-air mixtures to achieve high brake thermal efficiencies. To support this, new engine designs have introduced high stroke-to-bore ratios and cylinder head designs that promote high tumble flow and turbulence intensities. However, mixture dilution through exhaust gas recirculation (EGR) is limited by combustion instabilities manifested in the form of cycle-to-cycle variability. Propane has been observed to have superior EGR dilution tolerance than gasoline, which makes it a very competitive low-carbon fuel for the new IC engines without sacrificing efficiency. Two-dimensional direct numerical simulations (DNS) are performed with detailed chemistry to study and contrast the effect of turbulence intensity and dilution on propane and iso-octane premixed flames at high pressure conditions similar to those in-cylinder. A new reduced mechanism for propane consisting of 53 transported species and 17 quasi-steady state species is developed based on a previously published mechanism and used in these simulations. Three levels of turbulence intensity and two levels of exhaust gas dilution are chosen based on conditions relevant to IC engine operation. The DNS results are analyzed based on the evolution of the flame surface area and the statistics of its driving terms, which are found to be similar for both fuels when there is no dilution but considerably different under high dilution. The analysis of the DNS data provides fundamental insights into the underlying mechanisms for improved stability under dilution.

42 ENGINEERING↗

Simultaneous 10 kHz three-dimensional CH 2 O and tomographic PIV measurements in a lifted partially-premixed jet flame

High-speed, three-dimensional (3D) scalar-velocity field measurements are demonstrated in a lifted partially-premixed dimethyl-ether/air jet flame using simultaneous laser-induced fluorescence (LIF) of formaldehyde (CH 2 O) and tomographic particle image velocimetry (TPIV). The 3D LIF measurements are performed by raster scanning the laser beam from a 100 kHz pulse-burst laser across the probe volume using an acousto-optic deflector. The volumetric reconstruction of the LIF signal from ten parallel planes provides quasi-instantaneous 3D LIF measurements that are synchronized with 10 kHz TPIV measurements. Here, the temporally resolved CH 2 O-LIF and velocity field data are used to analyze Lagrangian particle trajectories and displacement speeds at the base of the lifted flame. The particle trajectories reveal flow structures that are difficult to observe in an Eulerian reference frame. Positive and negative displacement speeds are observed at the CH 2 O-LIF surfaces at the inner and outer regions of the jet flame with a maximum displacement speed of approximately eight times the laminar flame speed of a stoichiometric dimethyl-ether/air mixture. The observations in present study support formation of lean premixed flames near the flame base whose role in flame stabilization should be investigated in more detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flashback Studies of High-Hydrogen Flames using High-Speed OH Planar Laser-Induced Fluorescence

Flame flashback is a crucial challenge during the application of high-hydrogen fuels for the development of next-generation lean premixed combustors. The fundamentals of stable flame configuration, flashback phenomena, and spatiotemporal information of flame dynamics in a premixed, swirl-stabilized burner are investigated experimentally at atmospheric pressure conditions using high-repetition-rate nanosecond (ns)-duration hydroxyl radical planar laser-induced fluorescence (OH-PLIF) measurements. The inlet conditions are varied systematically with respect to equivalence ratio (ϕ), hydrogen percentage in a hydrogen-methane mixture (50%–90%) and pre-mixer velocities and the experiments are performed on the burner modified to provide optical access to the premixing section. Signal interferences such as saturation effects, background chemiluminescence, and laser non-uniformity issues were minimized or corrected during data acquisition and subsequent data processing steps. A stable flame characterization was performed at 20 kHz-repetition-rate for varied pre-mixer velocities, hydrogen content and ϕ in the reactant mixture and lift-off length (L) investigation showed dependence on each of these parameters. Flame dynamics such as flame flashback events, local ignition, local extinction, flame curvature as well as the main reaction zone of the combustion mechanism are well characterized. Flashback ϕ showed an expected linearly increasing trend with increasing pre-mixer velocity and decreasing hydrogen content and the conclusions drawn aligned well with detailed L investigation. Such spatially and temporally resolved kHz-rate OH-PLIF is a promising technique to observe rapidly occurring flashback dynamics and can further validate turbulence-chemistry interaction models of swirling flames.

Parajuli, Pradeep↗

DNS of ignition and flame stabilization in a simplified gas turbine premixer

With the increasing need for fuel flexibility, mitigation of auto-ignition (AI) inside gas turbine (GT) premixers becomes crucial. They must be designed to yield a sufficiently homogeneous fuel-air mixture to achieve low emissions while at the same time avoiding the occurrence of AI and subsequent flame stabilization. This challenge requires a detailed understanding of turbulent mixing and chemistry interactions. In the present work, a direct numerical simulation (DNS) of an array of jets in crossflow (JICF), representative of an industrial GT premixer, is reported to shed light on these complex phenomena. It is found that AI kernels form in the aft part of the premixer and coalesce into a flame front that then propagates upstream, mainly through the boundary layer, and successively engulfs the jets. This, therefore, suggests a significant role of the jet array pattern on the flame stabilization. It is noted that AI kernels continue to form independently during the whole time of the simulation. To clarify the contribution of AI and diffusion in the ignition kernels and the main flame, chemical explosive mode analysis (CEMA) is employed jointly with a kernel tracking algorithm. It is found that during the initial formation of the flame, many ignition kernels form in mixtures with low scalar dissipation rate and large contribution from AI mode. As they quickly grow, they merge into a single flame front that becomes increasingly more diffusion-assisted over time, balancing the AI mode. Turbulence is shown to have a significant enhancing effect in lean premixed flames, but further analysis is required to fully characterize it. These findings are relevant for the industrial premixer studied, and also for novel micromix concepts that may be used in the next generation of GT combustion systems.

ADVANCED PROPULSION SYSTEMS↗

A two-equation soot-in-flamelet modeling approach applied under Spray A conditions

Soot production (including formation and oxidation) is studied in the transient, high-pressure and turbulent n-dodecane Spray A flames from the Engine Combustion Network (ECN) using computational fluid dynamics (CFD) simulations. A two-equation soot-in-flamelet modeling approach is applied within the framework of the Unsteady Flamelet Progress Variable (UFPV) model and results are validated against experimental data. Equations for soot mass fraction and soot number density derived in the mixture fraction space are solved in the context of detailed flamelet calculations. Source terms for the different steps in the soot chemistry are tabulated and incorporated in the flamelet manifold. For the reference condition, the modeling approach based on the tabulated flamelet manifold reduces the computational cost of a CFD calculation by approximately 40 times compared to a non-tabulated well-mixed (WM) modeling approach. The soot-in-flamelet approach is then extended to study the effect of ambient oxygen concentration, ambient mixture composition and ambient temperature on soot production. Results show that the modeling approach is able to capture the experimental trends for the soot volume fraction (SVF) with good quantitative agreement, especially in the soot ramp-up region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗

Optical Investigation of a Partial Fuel Stratification Strategy to Stabilize Overall Lean Operation of a DISI Engine Fueled with Gasoline and E30

This paper offers new insights into a partial fuel stratification (PFS) combustion strategy that has proven to be effective at stabilizing overall lean combustion in direct injection spark ignition engines. To this aim, high spatial and temporal resolution optical diagnostics were applied in an optically accessible engine working in PFS mode for two fuels and two different durations of pilot injection at the time of spark: 210 µs and 330 µs for E30 (gasoline blended with ethanol by 30% volume fraction) and gasoline, respectively. In both conditions, early injections during the intake stroke were used to generate a well-mixed lean background. The results were compared to rich, stoichiometric and lean well-mixed combustion with different spark timings. In the PFS combustion process, it was possible to detect a non-spherical and highly wrinkled blue flame, coupled with yellow diffusive flames due to the combustion of rich zones near the spark plug. The initial flame spread for both PFS cases was faster compared to any of the well-mixed cases (lean, stoichiometric and rich), suggesting that the flame propagation for PFS is enhanced by both enrichment and enhanced local turbulence caused by the pilot injection. Different spray evolutions for the two pilot injection durations were found to strongly influence the flame kernel inception and propagation. PFS with pilot durations of 210 µs and 330 µs showed some differences in terms of shapes of the flame front and in terms of extension of diffusive flames. Yet, both cases were highly repeatable.

33 ADVANCED PROPULSION SYSTEMS↗

On the evolution of scalar iso-surface area density in a turbulent mixing layer

In this paper, the transport equation for scalar iso-surface area density (Σ) in a turbulent, temporally developing mixing layer is examined. Exploring the spatial and temporal evolution of the terms in the Σ transport equation is vital to improving our understanding of turbulent flows characterized by distinct interfaces, e.g. the flame surface or the turbulent/non-turbulent interface. Previous work reported by the authors identified that Σ exhibits self-similar behaviour consistent with the development of the temporal mixing layer. Accordingly, each of the terms in the Σ transport equation is found to behave in a self-similar manner, though there are notable differences in the self-similar behaviours for each term. Based on the results presented herein, it is suggested that the rate of change of Σ (∂Σ/∂t) and the advection term scale with hλ Φ /ΔU, where h is the width of the mixing layer, λ Φ is the scalar Taylor length scale and ΔU is the velocity difference. The production and destruction terms are found to scale with an additional factor (ReSc) 1/2 . In contrast, the molecular diffusion term is found to scale with a factor (ReSc) –1/2 compared to ∂Σ/∂t. Importantly, it is found that the difference between the production and destruction terms, or net surface ‘stretch’, scales with the same factor as ∂Σ/∂t and the advection term, which may have a significant impact on how the evolution of Σ is understood and modelled in turbulent flows.

42 ENGINEERING↗

A co-kurtosis PCA based dimensionality reduction with nonlinear reconstruction using neural networks

For turbulent reacting flow systems, identification of low-dimensional representations of the thermo-chemical state space is vitally important, primarily to significantly reduce the computational cost of device-scale simulations. Principal component analysis (PCA), and its variants, are a widely employed class of methods. Recently, an alternative technique that focuses on higher-order statistical interactions, co-kurtosis PCA (CoK-PCA), has been shown to effectively provide a low-dimensional representation by capturing the stiff chemical dynamics associated with spatiotemporally localized reaction zones. While its effectiveness has only been demonstrated based on a priori analyses with linear reconstruction, in this work, we employ nonlinear techniques to reconstruct the full thermo-chemical state and evaluate the efficacy of CoK-PCA compared to PCA. Specifically, we combine a CoK-PCA-/PCA-based dimensionality reduction (encoding) with an artificial neural network (ANN) based reconstruction (decoding) and examine, a priori, the reconstruction errors of the thermo-chemical state. In addition, we evaluate the errors in species production rates and heat release rates, which are nonlinear functions of the reconstructed state, as a measure of the overall accuracy of the dimensionality reduction technique. We employ four datasets to assess CoK-PCA/PCA coupled with ANN-based reconstruction: zero-dimensional (homogeneous) reactor for autoignition of an ethylene/air mixture that has conventional single-stage ignition kinetics, a dimethyl ether (DME)/air mixture which has two-stage (low and high temperature) ignition kinetics, a one-dimensional freely propagating premixed ethylene/air laminar flame, and a two-dimensional dataset representing turbulent autoignition of ethanol in a homogeneous charge compression ignition (HCCI) engine. Finally, results from the analyses demonstrate the robustness of the CoK-PCA based low-dimensional manifold with ANN reconstruction in accurately capturing the data, specifically from the reaction zones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗

Towards stabilization of pressurized oxy-coal diffusion flames

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. Turbulent jet diffusion burners are adopted for the pressurized oxy-combustor in this research and a bluff-body is employed to create a recirculation zone to stabilize the flame in pressurized oxy-combustor. The objective of this numerical work is to perform a systematic analysis of the characteristics of this pressurized non-premixed flame. Specifically, a 15-bar POFC combustor of power 100 kWth is modeled employing the Ansys FLUENT commercial platform, using Reynolds-averaged Navier-Stokes (RANS) modeling. The present computational work focuses on identifying the aerothermodynamic features of the isothermal and reacting flow with a disk-shape bluff-body. Results show that fuel-oxidizer stream momentum ratio has great impacts on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. Parametric study of the width of the burner tip identifies an optimal range for this pressurized burner.

Li, Lei↗