Mitigating chemical stiffness in turbulent reacting flow simulations via an adaptive reduced-order model with time-dependent bases
Abstract not provided.
Abstract not provided.
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.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Combustion and Flame
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.
Combustion and Flame
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.
Combustion and Flame
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.
Proceedings of the Combustion Institute
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.
Proceedings of the Combustion Institute
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.
Proceedings of the Combustion Institute
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.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Combustion and Flame
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.
Abstract not provided.
Proceedings of the Combustion Institute
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.