The effects of combustion on the strain rate field are investigated in turbulent premixed CH4/air Bunsen flames using simultaneous tomographic PIV and OH LIF measurements. Tomographic PIV provides three-dimensional velocity measurements, from which the complete strain rate tensor is determined. The OH LIF measurements are used to determine the position of the flame surface and the flame-normal orientation within the imaging plane. This combination of diagnostic techniques enables quantification of divergence as well as flame-normal and tangential strain rates, which are otherwise biased using only planar measurements. Measurements are compared in three lean-to-stoichiometric flames that have different amounts of heat release and Damköhler numbers greater than unity. The effects of heat release on the principal strain rates and their alignment relative to the local flame normal are analyzed. The extensive strain rate preferentially aligns with the flame normal in the reaction zone, which has been indicated by previous studies. The strength of this alignment increases with increasing heat release and, as a result, the flame-normal strain rate becomes highly extensive. These effects are associated with the gas expansion normal to the flame surface, which is largest for the stoichiometric flame. In the preheat zone, the compressive strain rate has a tendency to align with the flame normal. Away from the flame front, the flame – strain rate alignment is arbitrary in both the reactants and products. The flame-tangential strain rate is on average positive across the flame front, and therefore the turbulent strain rate field contributes to the enhancement of scalar gradients as in passive scalar turbulence. As a result, increases in heat release result in larger positive values of the divergence as well as flame-normal and tangential strain rates, the tangential strain rate has a weaker dependence on heat release than the flame-normal strain rate and the divergence.
Stabilization methods of turbulent flames often involve mixing of reactants with hot products of combustion. The stabilizing effect of combustion product enthalpy has been long recognized, but the role played by the chemical composition of the product gases is typically overlooked. We employ a counterflow system to pinpoint the effects of the combustion product stoichiometry on the structure of turbulent premixed flames under conditions of both stable burning and local extinction. To that end, a turbulent jet of lean-to-rich, CH4/O2/N2-premixed reactants at a turbulent Reynolds number of 1050 was opposed to a stream of hot products of combustion that were generated in a preburner. While the combustion product stream temperature was kept constant, its stoichiometry was varied independently from that of the reactant stream, leading to reactant-to-product stratification of relevance to practical combustion systems. The detailed structure of the turbulent flame front was analyzed in two series of experiments using laser-induced fluorescence (LIF): joint CH2O LIF and OH LIF measurements and joint CO LIF and OH LIF measurements. Results revealed that a decrease in local CH2O+OH and CO+OH reaction rates coincide with the depletion of OH radicals in the vicinity of the combustion product stream. These critical combustion reaction rates were more readily quenched in the presence of products of combustion from a stoichiometric flame, whereas they were favored by lean combustion products. As a result, stoichiometric combustion products contributed to a greater occurrence of local extinction. Furthermore, they limited the capacity of premixed reactants to ignite and of the turbulent premixed flames to stabilize. In contrast, lean and rich combustion products facilitated flame ignition and stability and reduced the rate of local extinction. The influence of the combustion product stream on the turbulent flame front was limited to a zone of approximately two millimeters from the gas mixing layer interface (GMLI) of the product stream. Flame fronts that were separated from the GMLI by larger distances were unaffected by the product stream stoichiometry.
In turbulent flows, the interaction between vorticity, ω, and strain rate, s, is considered a primary mechanism for the transfer of energy from large to small scales through vortex stretching. The ω-s coupling in turbulent jet flames is investigated using tomographic particle image velocimetry (TPIV). TPIV provides a direct measurement of the three-dimensional velocity field from which ω and s are determined. The effects of combustion and mean shear on the ω-s interaction are investigated in turbulent partially premixed methane/air jet flames with high and low probabilities of localized extinction as well as in a non-reacting isothermal air jet with Reynolds number of approximately 13 000. Results show that combustion causes structures of high vorticity and strain rate to agglomerate in highly correlated, elongated layers that span the height of the probe volume. In the non-reacting jet, these structures have a more varied morphology, greater fragmentation, and are not as well correlated. The enhanced spatiotemporal correlation of vorticity and strain rate in the stable flame results in stronger ω-s interaction characterized by increased enstrophy and strain-rate production rates via vortex stretching and straining, respectively. The probability of preferential local alignment between ω and the eigenvector of the intermediate principal strain rate, s2, which is intrinsic to the ω-s coupling in turbulent flows, is larger in the flames and increases with the flame stability. The larger mean shear in the flame imposes a preferential orientation of ω and s2 tangential to the shear layer. The extensive and compressive principal strain rates, s1 and s3, respectively, are preferentially oriented at approximately 45° with respect to the jet axis. The production rates of strain and vorticity tend to be dominated by instances in which ω is parallel to the s1 - s2 plane and orthogonal to s3.
In this LDRD project, we developed a capability for quantitative high - speed imaging measurements of high - pressure fuel injection dynamics to advance understanding of turbulent mixing in transcritical flows, ignition, and flame stabilization mechanisms, and to provide e ssential validation data for developing predictive tools for engine combustion simulations. Advanced, fuel - efficient engine technologies rely on fuel injection into a high - pressure, high - temperature environment for mixture preparation and com bustion. Howe ver, the dynamics of fuel injection are not well understood and pose significant experimental and modeling challenges. To address the need for quantitative high - speed measurements, we developed a Nd:YAG laser that provides a 5ms burst of pulses at 100 kHz o n a robust mobile platform . Using this laser, we demonstrated s patially and temporally resolved Rayleigh scattering imaging and particle image velocimetry measurements of turbulent mixing in high - pressure gas - phase flows and vaporizing sprays . Quantitativ e interpretation of high - pressure measurements was advanced by reducing and correcting interferences and imaging artifacts.
Turbulent dimethyl ether (DME) jet flames provide a canonical flame geometry for studying turbulence-flame interactions in oxygenated fuels and for developing predictive models of these interactions. The development of accurate models for DME/air flames would establish a foundation for studies of more complex oxygenated fuels. We present a joint experimental and computational investigation of the velocity field and OH and CH2O distributions in a piloted, partially-premixed turbulent DME/air jet flame with a jet exit Reynolds number, ReD, of 29,300. The turbulent DME/air flame is analogous to the well-studied, partially-premixed methane/air jet flame, Sandia Flame D, with identical stoichiometric mixture fraction, ξst = 0.35, and bulk jet exit velocity, Vbulk = 45.9 m/s. Measurements include particle image velocimetry (PIV) and simultaneous CH2O and OH laser-induced fluorescence (LIF) imaging. Simulations are performed using a large eddy simulation combined with conditional moment closure (LES-CMC) on an intermediate size grid of 1.3 million cells. Overall, the downstream evolution of the mean and RMS profiles of velocity, OH, and CH2O are well predicted, with the largest discrepancies occurring for CH2O at x/D = 20-25. LES-CMC simulations employing two different chemical reaction mechanisms (Kaiser et al., 2000 [20] and Zhao et al., 2008 [21]) show approximately a factor of two difference in the peak CH2O mole fractions, whereas OH mole fractions are in good agreement between the two mechanisms. The single-shot LIF measurements of OH and CH2O show a wide range of separation distances between the spatial distributions of these intermediate species with gaps on the order of millimeters. The intermittency in the overlap between these species indicates that the consumption rates of formaldehyde by OH in the turbulent DME/air jet flame may be highly intermittent with significant departures from flamelet models.
The effects of combustion on the strain rate field in turbulent jets were studied using 10 kHz tomographic particle image velocimetry (TPIV). Measurements were performed in three turbulent jets: a well-studied, piloted partially-premixed methane/air jet flame, Sandia flame C, with low probability of localized extinction; a second piloted jet flame, analogous to flame C but with a reduced pilot flow rate and a high probability of localized extinction; and a non-reacting air jet. Since the jet exit Reynolds number of approximately 13000 was nearly identical in the three jets, differences in the strain rate fields were attributed to the effects of combustion. Spatiotemporal characteristics of the strain rate field were analyzed. Overall, the strain rate norm was larger in the flames than in the non-reacting jet with the most stable flame having the largest values. In all three jets, the compressive strain rate was on average the largest of the three principal strain rates. At high strain rates, the ratios of the compressive and extensive strain rate to the intermediate strain rate were similar to those found in isotropic incompressible turbulent flows. The three-dimensional velocity measurements were used to analyze the spatial distribution of strain rate clusters, defined as singly-connected groups of voxels where the strain rate magnitude exceeded a threshold value. The presence of a stable flame significantly attenuated the number of clusters of intermediate strain rate. Strain rate bursts, corresponding to sudden increases in the number of clusters, were identified in the three jets. Bursts in the non-reacting jet and the unstable flame contained up to twice as many clusters as in the stable flame. The temporal intermittency of intense strain rate clusters was analyzed using the time-series measurements. Clusters with strain rates greater than five times the standard deviation of the strain rate norm were highly intermittent.
The University of Michigan and Sandia National Laboratories collaborated on the initial development of a compact single-camera approach for simultaneously measuring 3-D gasphase velocity and temperature fields at high frame rates. A compact diagnostic tool is desired to enable investigations of flows with limited optical access, such as near-wall flows in an internal combustion engine. These in-cylinder flows play a crucial role in improving engine performance. Thermographic phosphors were proposed as flow and temperature tracers to extend the capabilities of a novel, compact 3D velocimetry diagnostic to include high-speed thermometry. Ratiometric measurements were performed using two spectral bands of laser-induced phosphorescence emission from BaMg2Al10O17:Eu (BAM) phosphors in a heated air flow to determine the optimal optical configuration for accurate temperature measurements. The originally planned multi-year research project ended prematurely after the first year due to the Sandia-sponsored student leaving the research group at the University of Michigan.
Recent advances in high frame rate complementary metal-oxide-semiconductor (CMOS) cameras coupled with high repetition rate lasers have enabled laser-based imaging measurements of the temporal evolution of turbulent reacting flows. This measurement capability provides new opportunities for understanding the dynamics of turbulence-chemistry interactions, which is necessary for developing predictive simulations of turbulent combustion. However, quantitative imaging measurements using high frame rate CMOS cameras require careful characterization of the their noise, non-linear response, and variations in this response from pixel to pixel. We develop a noise model and calibration tools to mitigate these problems and to enable quantitative use of CMOS cameras. We have demonstrated proof of principle for image de-noising using both wavelet methods and Bayesian inference. The results offer new approaches for quantitative interpretation of imaging measurements from noisy data acquired with non-linear detectors. These approaches are potentially useful in many areas of scientific research that rely on quantitative imaging measurements.
This LDRD was a Sandia Fellowship that supported Andrea Hsu's PhD research at Texas A&M University and her work as a visitor at Sandia's Combustion Research Facility. The research project at Texas A&M University is concerned with the experimental characterization of hypersonic (Mach>5) flowfields using experimental diagnostics. This effort is part of a Multidisciplinary University Research Initiative (MURI) and is a collaboration between the Chemistry and Aerospace Engineering departments. Hypersonic flight conditions often lead to a non-thermochemical equilibrium (NTE) state of air, where the timescale of reaching a single (equilibrium) Boltzmann temperature is much longer than the timescale of the flow. Certain molecular modes, such as vibrational modes, may be much more excited than the translational or rotational modes of the molecule, leading to thermal-nonequilibrium. A nontrivial amount of energy is therefore contained within the vibrational mode, and this energy cascades into the flow as thermal energy, affecting flow properties through vibrational-vibrational (V-V) and vibrational-translational (V-T) energy exchanges between the flow species. The research is a fundamental experimental study of these NTE systems and involves the application of advanced laser and optical diagnostics towards hypersonic flowfields. The research is broken down into two main categories: the application and adaptation of existing laser and optical techniques towards characterization of NTE, and the development of new molecular tagging velocimetry techniques which have been demonstrated in an underexpanded jet flowfield, but may be extended towards a variety of flowfields. In addition, Andrea's work at Sandia National Labs involved the application of advanced laser diagnostics to flames and turbulent non-reacting jets. These studies included quench-free planar laser-induced fluorescence measurements of nitric oxide (NO) and mixture fraction measurements via Rayleigh scattering.
The catalytic effect of nitric oxide (NO) on the dynamics of extinction and re-ignition of a vortex-perturbed non-premixed hydrogen-air flame is studied in a counterflow burner. A diffusion flame is established with counterflowing streams of nitrogen-diluted hydrogen at ambient temperature and air heated to a range of temperatures that brackets the auto-ignition temperature. Localized extinction is induced by impulsively driving a fuel-side toroidal vortex into the steady flame, and the recovery of the extinguished region is monitored by planar laser-induced fluorescence (PLIF) of the hydroxyl radical (OH). The dynamics of flame recovery depend on the air temperature and fuel concentration, and four different recovery modes are identified. These modes involve combinations of edge-flame propagation and the expansion of an auto-ignition kernel that forms within the extinguished region. The addition of a small amount of NO significantly alters the re-ignition process by shifting the balance between chain-termination and chain-propagation reactions to enhance auto-ignition. The ignition enhancement by this catalytic effect causes a shift in the conditions that govern the recovery modes. In addition, the effects of NO concentration and vortex strength on the flame recovery are examined. Direct numerical simulations of the flame-vortex interaction with and without NO doping show how the small amount of OH produced by NO-catalyzed reactions has a significant impact on the development of an auto-ignition kernel. This joint experimental and numerical study provides detailed insight into the interaction between transient flows and ignition processes.
Combined experimental and numerical studies of the transient response of ignition to strained flows require a well-characterized ignition trigger. Laser deposition of a small radical pool provides a reliable method for initiating ignition of mixtures that are near the ignition limit. Two-dimensional direct numerical simulations are used to quantify the sensitivity of ignition kernel formation and subsequent edge-flame propagation to the oxidizer temperature and the initial width and amplitude of O-atom deposition used to trigger ignition in an axisymmetric counterflow of heated air versus ambient hydrogen/nitrogen. The ignition delay and super-equilibrium OH concentration in the nascent ignition kernel are highly sensitive to variations in these initial conditions. The ignition delay decreases as the amplitude of the initial O-atom deposition increases. The spatial distribution and the magnitude of the OH overshoot are governed by multi-dimensional effects. The degree of OH overshoot near the burner centerline increases as the diameter of the initial O-atom deposition region decreases. This result is attributed to preferential diffusion of hydrogen in the highly curved leading portion of the edge flame that is established following thermal runaway. The edge-flame speed and OH overshoot at the leading edge of the edge flame are relatively insensitive to variations in the initial conditions of the ignition. The steady edge-flame speed is approximately twice the corresponding laminar flame speed. The rate at which the edge flame approaches its steady state is insensitive to the initial conditions and depends solely on the diffusion time scale at the edge flame. The edge flame is curved toward the heated oxidizer stream as a result of differences in the chemical kinetics between the leading edge and the trailing diffusion flame. The structure of the highly diluted diffusion flame considered in this study corresponds to Liñán's 'premixed flame regime' in which only the oxidizer leaks through the reaction zone such that the flame is located at fuel lean rather than stoichiometric mixture fraction conditions.
