Nitrogen Oxide Emissions from Premixed Reacting Jets in a Vitiated Crossflow
Not provided.
SEARCH · Search NASA
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.
Not provided.
Many complex turbulent flows in nature and engineering can be qualitatively regarded as being constituted of multiple simpler unit flows. The objective of this work is to characterize the coherent structures in such complex flows as a combination of constituent unitary flow structures for the purpose of reduced-order representation. While turbulence is clearly a non-linear phenomenon, we aim to establish the degree to which the optimally weighted superposition of unitary flow structures can represent the complex flow structures. The rationale for investigating such superposition stems from the fact that the large-scale coherent structures are generated by underlying flow instabilities that may be reasonably described using linear analysis. Clearly, the degree of validity of superposition will depend on the flow under consideration. In this work, we take the first step toward establishing a procedure for investigating superposition. Experimental data of single and triple tandem jets in crossflow are used to demonstrate the procedure. A composite triple tandem jet flow field is generated from optimal superposition of single jet data and compared against ‘true’ triple jet data. Direct comparisons between the true and composite fields are made for spatial, temporal, and kinetic energy content. The large-scale features (obtained from proper orthogonal decomposition or POD) of true and composite tandem jet wakes exhibit nearly 70% agreement in terms of modal eigenvector correlation. Corresponding eigenvalues reveal that the kinetic energy of the flow is also emulated with only a slight overprediction. Temporal frequency features are also examined in an effort to completely characterize POD modes. The proposed method serves as a foundation for more rigorous and robust dimensional reduction in complex flows based on unit flow modes.
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.
Micromix combustion presents a promising solution for integrating hydrogen into gas turbine combustors, offering rapid and compact air-fuel mixing that effectively shortens the combustion zone and significantly reduces nitrogen oxides (NOX) emissions. A major challenge in hydrogen combustion systems is flashback, which is mitigated through short mixing lengths and high flow velocities. This study investigates three fuel injector designs – zero swirl, low swirl, and high swirl – for a 6-mm inner diameter nozzle, where hydrogen is injected through six 0.3 mm diameter holes in a jet-in-crossflow configuration relative to the airstream. Numerical simulations indicate stable flame formation, moderate peak flame temperatures, and subsequent low NOx potential for all designs. These injectors have been fabricated using additive manufacturing, and upcoming experiments will explore flame structure and NOx behavior across varying fuel injection angles, equivalence ratios and pre-mixer velocities. These efforts aim to advance the development of low-emission hydrogen micromix combustion.
A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.
High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.
NO x emissions remain a primary concern for modern and future gas turbines, particularly as we progress towards a net-zero carbon future. It is critical to identify novel strategies to mitigate environmental pollutants for both modern turbines operated on natural gas and future turbines projected to use carbon-free fuels (such as hydrogen). In the current study, NO x emissions of a reacting methane-air jet in a vitiated crossflow are experimentally investigated in a model axially staged combustor at 5 atm. The combustion products from the main stage combustor flow into the axial stage at temperatures ranging from 1580 to 1650 °C. The axial stage contains a reacting jet in crossflow, which provides an overall temperature rise ranging from 50 to 220 °C to generate exit temperatures similar to gas turbine engines. The focus of this paper is to explore the effects of flame liftoff and ignition timescales on the NO x contribution of this secondary jet. This is done by varying the premixed level and preheat temperature of the axial jet. The effects are explored at different momentum flux ratios (J), jet equivalence ratios (φ jet ), exit temperatures (T exit ), and main stage (headend) temperatures (T HE ). High-speed CH* chemiluminescence imaging is employed for each test case to determine the flame stabilization point. For a given temperature rise and combustor exit temperature, an increase in ignition delay was observed with a decrease in jet temperature. This provided a NO x benefit, particularly for jets that are injected at or close to an ignitable mixture fraction. Similarly, the coaxial injector caused a significant ignition delay compared to the fully premixed injector; this also led to a NO x benefit. Overall, the results show that a thermally driven delay in the axial stage provides a greater NO x reduction than a mixing delay.
With the increasing need for fuel-flexibility, an understanding of auto-ignition inside gas turbine premixers is crucial. Direct numerical simulation (DNS) is a great tool to reveal the fine grain details of physics inside the premixer and guide the choice of the models to use for less expensive simulations. We present numerical simulations of a simplified and downscaled premixer, with multiple air jets in rich fuel crossflow. The simulations show intense fuel-air mixing and flame-turbulence interaction. We examine flow acceleration by the incoming jets, and the combustion within the device. We observe low-temperature and high-temperature combustion signatures, and compute the flame index to quantify local stabilization modes.
The current study explores the possibility of cooling the vanes and blades of a direct-fired sCO 2 turbine using film cooling geometries. The operating conditions of a direct-fired sCO 2 cycle and thermophysical properties of the fluid at those conditions can alter the flow field characteristics of the coolant jet and its mixing with the mainstream. Very little information is present in the literature regarding the performance of film cooling geometries employing a supercritical fluid. The objective of this study is to estimate the resulting film cooling effectiveness while also capturing the effects of the crossflow-to-mainstream velocity ratio on the coolant jet. A computational fluid dynamic model is used to study the coolant jet exiting a cylindrical hole located on a flat plate, with the coolant fed by an internal channel. Steady-state RANS equations were solved along with the SST kω model to provide the turbulence closure. The operating conditions for the direct-fired sCO 2 turbine are obtained using an in-house Cooled Turbine Model. Numerical predictions revealed that the crossflow effects and jet lift-off were more pronounced in the case of sCO 2 when compared to air. Spatial distribution of flow field and cooling effectiveness will be presented at different operating conditions. Recommendations for future studies will be given wherever applicable.
The current study explores the possibility of cooling the vanes and blades of a direct-fired sCO 2 turbine using film cooling. The operating conditions of a direct-fired sCO 2 cycle and thermophysical properties of the fluid at those conditions can alter the flow field characteristics of the coolant jet and its mixing with the mainstream. Very little information is present in the literature regarding the performance of film cooling geometries employing supercritical CO 2 . The objective of this study is to estimate the resulting film cooling effectiveness while also capturing the effects of the crossflow-to-mainstream velocity ratio on the coolant jet. Herein, a computational fluid dynamic model is used to study the coolant jet exiting a cylindrical hole located on a flat plate, with the coolant fed by an internal channel. Steady-state Reynolds-averaged Navier–Stokes equations were solved along with the (shear-stress transport) SST k–ω model to provide the turbulence closure. The operating conditions for the direct-fired sCO 2 turbine are obtained using an in-house Cooled Turbine Model. Numerical predictions revealed that the crossflow effects and jet lift-off were more pronounced in the case of sCO 2 when compared to air. Spatial distribution of flow field and cooling effectiveness are presented at different operating conditions.
The primary parameter of a standard k-ϵ model, C μ , was calculated from stereoscopic particle image velocimetry (PIV) data for a supersonic jet exhausting into a transonic crossflow. This required the determination of turbulent kinetic energy, turbulent eddy viscosity, and turbulent energy dissipation rate. Image interrogation was optimized, with different procedures used for mean strain rates and Reynolds stresses, to produce useful turbulent eddy viscosity fields. The eddy viscosity was calculated by a least-squares fit to all components of the three-dimensional strain-rate tensor that were available from the PIV data. This eliminated artifacts and noise observed when using a single strain component. Local dissipation rates were determined via Kolmogorov’s similarity hypotheses and the second-order structure function. The eddy viscosity and dissipation rates were then combined to determine C μ . Considerable spatial variation was observed in C μ , with the highest values found in regions where turbulent kinetic energy was relatively low but where turbulent mixing was important, e.g., along the high-strain jet edges and in the wake region. This suggests that use of a constant C μ in modeling may lead to poor Reynolds stress predictions at mixing interfaces. A data-driven modeling approach that can predict this spatial variation of C μ based on known state variables may lead to improved simulation results without the need for calibration.
Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.
Counter-rotating vortices, formed by the interaction of film-cooling jets and the hot gas flow, adversely affect the performance of conventional film-cooling designs. Downstream vortex generators have been shown to improve cooling effectiveness by mitigating the effects of the counter-rotating vortices and by deflecting the cooling jet laterally. In this study, computational and experimental methods were used to examine how cylindrical film-cooling holes (D = 3.2 mm, L/D = 6, p/D = 3, α = 30 deg) with and without downstream vortex generators perform when the coolant supply channel is perpendicular to the direction of the hot gas. For this study, the hot gas had a temperature of 650 K and an average Mach number of 0.23. Here, the hot-gas-to-coolant temperature ratio was 1.9, and two blowing ratios (0.75 and 1.0) were studied. Results from the computational fluid dynamics study show how crossflow affects the interaction between the film-cooling jet and hot gas flow with and without downstream vortex generators. The experimental measurements were based on infrared thermography in a conjugate heat transfer environment. Results were obtained for film-cooling performance in terms of overall effectiveness, film effectiveness, and local heat transfer coefficients. The downstream vortex generators can increase the laterally averaged effectiveness by a factor of 1.5 relative to cylindrical holes, but this higher performance is restricted to low crossflow velocities and higher blowing ratios.
This presentation describes experimental and CFD results of a novel idea to improve film cooling effectiveness. This concept involves downstream vortex generators to control counter-rotating vortices that can degrade discrete film cooling jets. The performance can be influenced significantly by how the coolant is supplied to the film cooling holes. In this paper, the effects of a coolant supply channel that is normal to the hot gas path are investigated and discussed.
Abstract Film cooling is used to protect turbine components from the extreme temperatures by ejecting coolant through arrays of holes to create an air buffer from the hot combustion gases. Limitations in traditional machining meant film cooling holes universally have sharp inlets, which create separation regions at the hole entrance. The present study uses experimental and computational data to show that these inlet separation are a major cause of performance variation in crossflow fed film cooling holes. Three-hole designs were experimentally tested by independently varying the coolant velocity ratio (VR) and the coolant channel velocity ratio (VRc) to isolate the effects of crossflow on hole performance. Leveraging additive manufacturing (AM) technologies, the addition of a 0.25D radius fillet to the inlet of a 7-7-7 shaped hole is shown to significantly improve diffuser usage and significantly reduce variation in performance with VRc. A second AM design used a very large radius of curvature inlet to reduce biasing caused by the inlet crossflow. Experiments showed that this “swept” hole design did minimize biasing of the coolant flow to one side of the shaped hole, and it significantly reduced variations due to varying VRc. RANS simulations at six VR and three VRc conditions were made for each geometry to better understand how the new geometries changed the velocity field within the hole. The sharp and rounded inlets were seen to have very similar tangential velocity fields and jet biasing. Both AM inlets created more uniform, slower velocity fields entering the diffuser. The results of this article indicate that large improvements in film cooling performance can be found by leveraging AM technology.