There is an increasing aspiration to utilize machine learning (ML) for various tasks of relevance to national security. ML models have thus far been mostly applied to tasks and domains that, while impactful, have sufficient volume of data. For predictive tasks of national security relevance, ML models of great capacity (ability to approximate nonlinear trends in input-output maps) are often needed to capture the complex underlying physics. However, scientific problems of relevance to national security are often accompanied by various sources of sparse and/or incomplete data, including experiments and simulations, across different regimes of operation, of varying degrees of fidelity, and include noise with different characteristics and/or intensity. State-of-the-art ML models, despite exhibiting superior performance on the task and domain they were trained on, may suffer detrimental loss in performance in such sparse data environments. This report summarizes the results of the Laboratory Directed Research and Development project entitled Trust-Enhancing Probabilistic Transfer Learning for Sparse and Noisy Data Environments. The objective of the project was to develop a new transfer learning (TL) framework that aims to adaptively blend the data across different sources in tackling one task of interest, resulting in enhanced trustworthiness of ML models for mission- and safety-critical systems. The proposed framework determines when it is worth applying TL and how much knowledge is to be transferred, despite uncontrollable uncertainties. The framework accomplishes this by leveraging concepts and techniques from the fields of Bayesian inverse modeling and uncertainty quantification, relying on strong mathematical foundations of probability and measure theories to devise new uncertainty-aware TL workflows.
Direct Numerical Simulations (DNS) are performed to investigate the process of spontaneous ignition of hydrogen flames at laminar, turbulent, adiabatic and non-adiabatic conditions. Mixtures of hydrogen and vitiated air at temperatures representing gas-turbine reheat combustion are considered. Adiabatic spontaneous ignition processes are investigated first, providing a quantitative characterization of stable and unstable flames. Results indicate that, in hydrogen reheat combustion, compressibility effects play a key role in flame stability and that unstable ignition and combustion are consistently encountered for reactant temperatures close to the mixture's characteristic crossover temperature. Furthermore, it is also found that the characterization of the adiabatic processes is also valid in the presence of non-adiabaticity due to wall heat-loss. Finally, a quantitative characterization of the instantaneous fuel consumption rate within the reaction front is obtained and of its ability, at auto-ignitive conditions, to advance against the approaching turbulent flow of the reactants, for a range of different turbulence intensities, temperatures and pressure levels.
In the present work, premixed combustion in a turbulent boundary layer under auto-ignitive conditions is investigated using direct numerical simulation (DNS). The turbulent inflow of the reactive DNS is obtained by temporal sampling of a corresponding inert DNS of a turbulent boundary layer at a location with Reτ= 360, where Reτ is the friction Reynolds number. The reactants of the DNS are determined by mixing the products of lean natural gas combustion and a H2/N2 fuel jet, resulting in a lean mixture of high temperature with a short ignition delay time. In the free stream the reaction front is stabilized at a streamwise location which can be predicted using the free stream velocity U∞ and the ignition delay time τig. Inside the boundary layer, combustion modifies the near-wall coherent turbulent structures considerably and turbulence results in reaction front wrinkling. The combustion modes in various regions were examined based on the results of displacement velocity, species budget and chemical explosive mode analysis (CEMA). It was indicated that flame propagation prevails in the near-wall region and auto-ignition becomes increasingly important as the wall-normal distance increases. The interactions of turbulence and combustion were studied through statistics of reaction front normal vector and strain rate tensor. It was found that the reaction front normal preferentially aligns with the most compressive strain rate in regions where the effects of heat release on the strain rate are minor and with the most extensive strain rate where its effects are significant. Negative correlations between the wall heat flux and flame quenching distance were observed. A new quenching mode, back-on quenching, was identified. It was found that the heat release rate at the wall is the highest when head-on quenching occurs and lowest when back-on quenching occurs.
It has been previously demonstrated that thermal gas expansion might have a role in boundary layer flashback of premixed turbulent flames [Gruber et al., J Fluid Mech 2012], inducing local flow-reversal in the boundary layer's low-velocity streaks on the reactants’ side of the flame and facilitating its upstream propagation. We perform a two-dimensional numerical investigation of the interaction between a periodic shear flow and a laminar premixed flame. The periodic shear is a simplified model for the oncoming prolonged streamwise velocity streaks with alternating regions of high and low velocities found in turbulent boundary layers in the vicinity of the walls. The parametric study focuses on the amplitude and wavelength of the periodic shear flow and on the gas expansion ratio (unburnt-to-burnt density ratio). With the increase of the amplitudes of the periodic shear flow and of the gas expansion, the curved flame velocity increases monotonically. The flame velocity dependence on the periodic shear wavelength is non-monotonic, which is consistent with previous theoretical studies of curved premixed flame velocity. The flame shape that is initially formed by the oncoming periodic shear appears to be metastable. At a later stage of the flame propagation, the flame shape transforms into the stationary one dominated by the Darrieus-Landau instability.
A flamelet analysis of a highly resolved direct numerical simulation (DNS) of a multi-injection flame with both auto-ignition and ignition induced by flame-flame interaction was conducted. A novel method was proposed to identify the different combustion modes of ignition processes using generalized flamelet equations. A state-of-the-art DNS database for a multi-injection n-dodecane flame in a diesel engine environment was investigated. Three-dimensional flamelets were extracted from the DNS at different time instants with a focus on auto-ignition and interaction-ignition processes. The influences of mixture field interactions and the scalar dissipation rate on the ignition process were examined by varying the species composition boundary conditions of the transient flamelet equations. Results showed that auto-ignition is delayed if the burned products are added to the oxidizer side of the flamelet, and the ignition delay time is sensitive to the scalar dissipation rate. The significance of mass diffusion in the flame-normal direction is reduced due to the existence of burned products in the oxidizer stream. Budget analyses of the generalized flamelet equations revealed that the transport along the mixture fraction iso-surface is insignificant during the auto-ignition process, but becomes important when interaction-ignition occurs, which is further confirmed through a flamelet regime classification method.
In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.
Accurate modeling of mixing in large-eddy simulation (LES)/transported probability density function (PDF) modeling of turbulent combustion remains an outstanding issue. The issue is particularly salient in turbulent premixed combustion under extreme conditions such as high-Karlovitz number Ka. The present study addresses this issue by conducting an a priori analysis of a power-law scaling based mixing timescale model for the transported PDF model. A recently produced DNS dataset of a high-Ka turbulent jet flame is used for the analysis. A power-law scaling is observed for a scaling factor used to model the sub-filter scale mixing timescale in this high-Ka turbulent premixed DNS flame when the LES filter size is much greater than the characteristic thermal thickness of a laminar premixed flame. The sensitivity of the observed power-law scaling to the different viewpoints (local or global) and to the different scalars for the data analysis is examined and the dependence of the model parameters on the dimensionless numbers Ka and Re (the Reynolds number) is investigated. Different model formulations for the mixing timescale are then constructed and assessed in the DNS flame. The proposed model is found to be able to reproduce the mixing timescale informed by the high-Ka DNS flame significantly better than a previous model.
Krisman, Alex; Meagher, Patrick; Zhao, Xinyu; Park, Ji-Woong; Lu, Tianfeng; Chen, Jacqueline H.
The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.
In this paper, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.
In this paper, the first flamelet analysis is conducted of a highly resolved DNS of a multi-injection flame with both auto-ignition and ignition induced by flame-flame interaction. A novel method is proposed to identify the different combustion modes of ignition processes using generalized flamelet equations. The state-of-the-art DNS database generated by Rieth et al. (US National Combustion Meeting, 2019) for a multi-injection flame in a Diesel engine environment is investigated. Three-dimensional flamelets are extracted from the DNS at different time instants with a focus on auto-ignition and interaction-ignition processes. The influences of mixture field interactions and the scalar dissipation rate on the ignition process are investigated by varying the species composition boundary conditions of the transient flamelet equations. Budget analyses of the generalized flamelet equations show that the transport along the mixture fraction iso-surface is insignificant during the auto-ignition process, but becomes important when interaction-ignition occurs, which is further confirmed through a flamelet regime classification method.
It has been predicted that several changes will occur when premixed flames are subjected to the extreme levels of turbulence that can be found in practical combustors. This paper is a review of recent experimental and DNS results that have been obtained for the range of extreme turbulence, and it includes a discussion of cases that agree or disagree with predictions. “Extreme turbulence” is defined to correspond to a turbulent Reynolds number (ReΤ, based on integral scale) that exceeds 2800 or a turbulent Karlovitz number that exceeds 100, for reasons that are discussed in Section 2.1. Several data bases are described that include measurements made at Lund University, the University of Sydney, the University of Michigan and the U.S. Air Force Research Lab. The data bases also include DNS results from Sandia National Laboratory, the University of New South Wales, Newcastle University, the California Institute of Technology and the University of Cambridge. Several major observations are: (a) DNS now can be achieved for a realistic geometry (of the Lund University jet burner) even for extreme turbulence levels, (b) state relations (conditional mean profiles) from DNS and experiments do tend to agree with laminar profiles, at least for methane-air and hydrogen-air reactants that are not preheated, and (c) regime boundaries have been measured and they do not agree with predicted boundaries. These findings indicate that the range of conditions for which flamelet models should be valid is larger than what was previously believed. Additional parameters have been shown to be important; for example, broken reactions occur if the “back-support” is insufficient due to the entrainment of cold gas into the product gas. Turbulent burning velocity measurements have been extended from the previous normalized turbulence levels (u’/SL) of 24 up to a value of 163. Turbulent burning velocities no longer follow the trend predicted by Shchelkin but they tend to follow the trend predicted by Damköhler. The boundary where flamelet broadening begins was measured to occur at ReTaylor = 13.8, which corresponds to an integral scale Reynolds number (ReT) of 2800. This measured regime boundary can be explained by the idea that flame structure is altered when the turbulent diffusivity at the Taylor scale exceeds a critical value, rather than the idea that changes occur when Kolmogorov eddies just fit inside a flamelet. A roadmap to extend DNS to complex chemistry and to higher Reynolds numbers is discussed. Exascale computers, machine learning, adaptive mesh refinement and embedded DNS show promise. Some advances are reviewed that have extended the use of line and planar PLIF and CARS laser diagnostics to studies that consider complex hydrocarbon fuels and harsh environments.
A three-dimensional direct numerical simulation (DNS) is performed for a turbulent hydrogen-air flame, represented with detailed chemistry, stabilized in a model gas-turbine combustor. The combustor geometry consists of a mixing duct followed by a sudden expansion and a combustion chamber, which represents a geometrically simplified version of Ansaldo Energia's GT26/GT36 sequential combustor design. In this configuration, a very lean blend of hydrogen and vitiated air is prepared in the mixing duct and convected into the combustion chamber, where the residence time from the inlet of the mixing duct to the combustion chamber is designed to coincide with the ignition delay time of the mixture. The results show that when the flame is stabilized at its design position, combustion occurs due to both autoignition and flame propagation (deflagration) modes at different locations within the combustion chamber. A chemical explosive mode analysis (CEMA) reveals that most of the fuel is consumed due to autoignition in the bulk-flow along the centerline of the combustor, and lower amounts of fuel are consumed by flame propagation near the corners of the sudden expansion, where the unburnt temperature is reduced by the thermal wall boundary layers. An unstable operating condition is also identified, wherein periodic auto-ignition events occur within the mixing duct. These events appear upstream of the intended stabilization position, due to positive temperature fluctuations induced by pressure waves originating from within the combustion chamber. The present DNS investigation represents the initial step of a comprehensive research effort aimed at gaining detailed physical insight into the rate-limiting processes that govern the sequential combustor behavior and avoid the insurgence of the off-design auto-ignition events.
Shin, Dong H.; Richardson, Edward S.; Aparece-Scutariu, Vlad; Minamoto, Yuki; Chen, Jacqueline H.
The link between the distribution of fluid residence time and the distribution of reactive scalars is analysed using Direct Numerical Simulation data. Information about the reactive scalar distribution is needed in order to model the reaction terms that appear in Large Eddy and Reynolds-Averaged simulations of turbulent reacting flows. The lifted flame is simulated taking account of multi-step chemistry for dimethyl-ether fuel. Due to autoignition and flame propagation, the reaction progress increases with residence time. The variation of fluid residence time is evaluated by solving an Eulerian transport equation for the fluid age. The fluid age is a passive scalar with a spatially-uniform source term, meaning that its moments and dissipation rates in turbulent flows can be modelled using closures already established for conserved scalars such as mixture fraction. In combination with the mixture fraction, the fluid age serves as a useful mapping variable to distinguish younger less-reacted fluid near the inlet from older more-reacted fluid downstream. The local fluctuations of mixture fraction and fluid age have strong negative correlation and, building upon established presumed-pdf models for mixture fraction, this feature can be used to construct an accurate presumed-pdf model for the joint mixture fraction/fluid age pdf. It is demonstrated that the double-conditional first-order moment closure combined with the proposed presumed model for the joint pdf of mixture fraction and fluid age gives accurate predictions for unconditional reaction rates - both for pre-ignition radical species produced by low-temperature processes upstream of the flame base, and for major species that are produced at the flame front.
The ignition process in diesel engines is highly complex and incompletely understood. In the present study, two-dimensional direct numerical simulations are performed to investigate the ignition dynamics and their sensitivity to thermochemical and mixing parameters. The thermochemical and mixing conditions are matched to the benchmark Spray A experiment from the Engine Combustion Network. The results reveal a complex ignition process with overlapping stages of: low-temperature ignition (cool flames), rich premixed ignition, and nonpremixed ignition, which are qualitatively consistent with prior experimental and numerical investigations, however, this is the first time that fully-resolved simulations have been reported at the actual Spray A thermochemical condition. Parametric variations are then performed for the Damkohler number Da, oxidiser temperature, oxygen concentration, and peak mixture fraction (a measure of premixedness), to study their effect on the ignition dynamics. It is observed that with both increasing oxidiser temperature and decreasing oxygen concentration, that the cool flame moves to richer mixtures, the overlap in the ignition stages decreases, and the (nondimensional) time taken to reach a fully burning state increases. With increasing Da, the cool-flame speed is decreased due to lower mean mixing rates, which causes a delayed onset of high-temperature ignition. With increasing peak mixture fraction, the onset of each stage of ignition is not affected, but the overall duration of the ignition increases leading to a longer burn duration. Overall, the results suggest that turbulence-chemistry interactions play a significant role in determining the timing and location in composition space of the entire ignition process.
Laminar natural gas flames are investigated at engine-relevant thermochemical conditions where the ignition delay time τ is short due to very high ambient temperatures and pressures. At these conditions, it is not possible to measure or calculate well-defined values for the laminar flame speed sl, laminar flame thickness δl, and laminar flame time scale due to the explosive thermochemical state. The corresponding reference values, sR, δR, and that account for the effects of autoignition, are numerically estimated to investigate the enhancement of flame propagation, and the competition with autoignition that arises under nominally autoignitive conditions (characterised here by the number τ/τR). Large values of τ/τR indicate that autoignition is unimportant, values near or below unity indicate that flame propagation is not possible, and intermediate values indicate that a combination of both flame propagation and autoignition may be important, depending upon factors such as device geometry, turbulence, stratification, et cetera. The reference quantities are presented for a wide range of temperatures, equivalence ratios, pressures, and hydrogen concentrations, which includes conditions relevant to stationary gas turbine reheat burners and boosted spark ignition engines. It is demonstrated that the transition from flame propagation to autoignition is only dependent on residence time, when the results are non-dimensionalised by the reference values. The temporal evolution of the reference values are also reported for a modelled boosted SI engine. It is shown that the nominally autoignitive conditions enhance flame propagation, which may be an ameliorating factor for the onset of engine knock. The calculations are performed using a recently-developed, detailed 177 species mechanism for C0–C3 chemistry that is derived from theoretical chemistry and is suitable for a wide range of thermochemical conditions as it is not tuned or optimised for a particular operating condition.