A three-dimensional direct numerical simulation (DNS) of a turbulent lifted dimethyl ether (DME) slot jet flame was performed at elevated pressure to study interactions between chemical reactions with low-temperature heat release (LTHR), negative temperature coefficient (NTC) reactions and shear generated turbulence in a jet in a heated coflow. By conditioning on mixture fraction, local reaction zones and local heat release rate, the turbulent flame is revealed to exhibit a "pentabrachial" structure that was observed for a laminar DME lifted flame [Krisman et al., (2015)]. The propagation characteristics of the stabilization and triple points are also investigated. Potential stabilization points, spatial locations characterized by preferred temperature and mixture fraction conditions, exhibit autoignition characteristics with large reaction rate and negligible molecular diffusion. The actual stabilization point which coincides with the most upstream samples from the pool of potential stabilization points fovr each spanwise location shows passive flame structure with large diffusion. The propagation speed along the stoichiometric surface near the triple point is compared with the asymptotic value obtained from theory [Ruetsch et al., (1995)]. At stoichiometric conditions, the asymptotic and averaged DNS values of flame displacement speed deviate by a factor of 1.7. However, accounting for the effect of low-temperature species on the local flame speed increase, these two values become comparable. This suggests that the two-stage ignition influences the triple point propagation speed through enhancement of the laminar flame speed in a configuration where abundant low-temperature products from the first stage, low-temperature ignition are transported to the lifted flame by the high-velocity jet.
Zhao, Xin Y.; Bhagatwala, Ankit; Chen, Jacqueline H.; Haworth, Daniel C.; Pope, Stephen B.
In this study, the modeling of mixing by molecular diffusion is a central aspect for transported probability density function (tPDF) methods. In this paper, the newly-proposed shadow position mixing model (SPMM) is examined, using a DNS database for a temporally evolving di-methyl ether slot jet flame. Two methods that invoke different levels of approximation are proposed to extract the shadow displacement (equivalent to shadow position) from the DNS database. An approach for a priori analysis of the mixing-model performance is developed. The shadow displacement is highly correlated with both mixture fraction and velocity, and the peak correlation coefficient of the shadow displacement and mixture fraction is higher than that of the shadow displacement and velocity. This suggests that the composition-space localness is reasonably well enforced by the model, with appropriate choices of model constants. The conditional diffusion of mixture fraction and major species from DNS and from SPMM are then compared, using mixing rates that are derived by matching the mixture fraction scalar dissipation rates. Good qualitative agreement is found, for the prediction of the locations of zero and maximum/minimum conditional diffusion locations for mixture fraction and individual species. Similar comparisons are performed for DNS and the IECM (interaction by exchange with the conditional mean) model. The agreement between SPMM and DNS is better than that between IECM and DNS, in terms of conditional diffusion iso-contour similarities and global normalized residual levels. It is found that a suitable value for the model constant c that controls the mixing frequency can be derived using the local normalized scalar variance, and that the model constant a controls the localness of the model. A higher-Reynolds-number test case is anticipated to be more appropriate to evaluate the mixing models, and stand-alone transported PDF simulations are required to more fully enforce localness and to assess model performance.
The diffusive-thermal (D-T) instability of opposed nonpremixed tubular flames near extinction is investigated using two-dimensional (2-D) direct numerical simulations together with the linear stability analysis. Two different initial conditions (IC), i.e. the perturbed IC and the C-shaped IC are adopted to elucidate the effects of small and large amplitude disturbances on the formation of flame cells, similar to conditions found in linear stability analysis and experiments, respectively. The characteristics of the D-T instability of tubular flames are identified by a critical Damköhler number, DaC, at which the D-T instability first occurs and the corresponding number of flame cells for three different tubular flames with different flame radii. It is found that DaC predicted through linear stability analysis shows good agreement with that obtained from the 2-D simulations performed with two different ICs. The flame cell number, Ncell, from the 2-D simulations with the perturbed IC is also found to be equal to an integer close to the maximum wavenumber, kmax, obtained from the linear stability analysis. However, Ncell from the 2-D simulations with the C-shaped IC is smaller than kmax and Ncell found from the simulations with the perturbed IC. This is primarily because the strong reaction at the edges of the horseshoe-shaped cellular flame developed from the C-shaped IC is more likely to produce larger flame cells and reduce Ncell. It is also found that for cases with the C-shaped IC, once the cellular instability occurs, the number of flame cells remains constant until global extinction occurs by incomplete reaction manifested by small Da. It is also verified through the displacement speed, Sd, analysis that the two edges of the horseshoe-shaped cellular flame are stationary and therefore do not merge due to the diffusion-reaction balance at the edges. Moreover, large negative Sd is observed at the local extinction points while small positive or negative Sd features in the movement of flame cells as they adjust their location and size towards steady state.
This paper presents results from one and two-dimensional direct numerical simulations under Reactivity Controlled Compression Ignition (RCCI) conditions of a primary reference fuel (PRF) mixture consisting of n-heptane and iso-octane. RCCI uses in-cylinder blending of two fuels with different autoignition characteristics to control combustion phasing and the rate of heat release. These simulations employ an improved model of compression heating through mass source/sink terms developed in a previous work by Bhagatwala et al. (2014), which incorporates feedback from the flow to follow a predetermined experimental pressure trace. Two-dimensional simulations explored parametric variations with respect to temperature stratification, pressure profiles and n-heptane concentration. Statistics derived from analysis of diffusion/reaction balances locally normal to the flame surface were used to elucidate combustion characteristics for the different cases. Both deflagration and spontaneous ignition fronts were observed to co-exist, however it was found that higher n-heptane concentration provided a greater degree of flame propagation, whereas lower n-heptane concentration (higher fraction of iso-octane) resulted in more spontaneous ignition fronts. A significant finding was that simulations initialized with a uniform initial temperature and a stratified n-heptane concentration field, resulted in a large fraction of combustion occurring through flame propagation. It was also found that the proportion of spontaneous ignition fronts increased at higher pressures due to shorter ignition delay when other factors were held constant. For the same pressure and fuel concentration, the contribution of flame propagation to the overall combustion was found to depend on the level of thermal stratification, with higher initial temperature gradients resulting in more deflagration and lower gradients generating more ignition fronts. Statistics of ignition delay are computed to assess the Zel'dovich (1980) theory for the mode of combustion propagation based on ignition delay gradients.
Application resilience is a key challenge that has to be addressed to realize the exascale vision. Online recovery, even when it involves all processes, can dramatically reduce the overhead of failures as compared to the more traditional approach where the job is terminated and restarted from the last checkpoint. In this paper we explore how local recovery can be used for certain classes of applications to further reduce overheads due to resilience. Specifically we develop programming support and scalable runtime mechanisms to enable online and transparent local recovery for stencil-based parallel applications on current leadership class systems. We also show how multiple independent failures can be masked to effectively reduce the impact on the total time to solution. We integrate these mechanisms with the S3D combustion simulation, and experimentally demonstrate (using the Titan Cray-XK7 system at ORNL) the ability to tolerate high failure rates (i.e., node failures every 5 seconds) with low overhead while sustaining performance, at scales up to 262144 cores.
Proceedings of ISAV 2015: 1st International Workshop on In Situ Infrastructures for Enabling Extreme-Scale Analysis and Visualization, Held in conjunction with SC 2015: The International Conference for High Performance Computing, Networking, Storage and Analysis
Next generation architectures necessitate a shift away from traditional workflows in which the simulation state is saved at prescribed frequencies for post-processing analysis. While the need to shift to in situ workflows has been acknowledged for some time, much of the current research is focused on static workflows, where the analysis that would have been done as a post-process is performed concurrently with the simulation at user-prescribed frequencies. Recently, research efforts are striving to enable adaptive workflows, in which the frequency, composition, and execution of computational and data manipulation steps dynamically depend on the state of the simulation. Adapting the workflow to the state of simulation in such a data-driven fashion puts extremely strict efficiency requirements on the analysis capabilities that are used to identify the transitions in the workflow. In this paper we build upon earlier work on trigger detection using sublinear techniques to drive adaptive workflows. Here we propose a methodology to detect the time when sudden heat release occurs in simulations of turbulent combustion. Our proposed method provides an alternative metric that can be used along with our former metric to increase the robustness of trigger detection. We show the effectiveness of our metric empirically for predicting heat release for two use cases.
Computing distance fields is fundamental to many scientific and engineering applications. Distance fields can be used to direct analysis and reduce data. In this paper, we present a highly scalable method for computing 3D distance fields on massively parallel distributed-memory machines. A new distributed spatial data structure, named parallel distance tree, is introduced to manage the level sets of data and facilitate surface tracking over time, resulting in significantly reduced computation and communication costs for calculating the distance to the surface of interest from any spatial locations. Our method supports several data types and distance metrics from real-world applications. We demonstrate its efficiency and scalability on state-of-the-art supercomputers using both large-scale volume datasets and surface models. We also demonstrate in-situ distance field computation on dynamic turbulent flame surfaces for a petascale combustion simulation. Our work greatly extends the usability of distance fields for demanding applications.
Post-Moore's law scaling is creating a disruptive shift in simulation workflows, as saving the entirety of raw data to persistent storage becomes expensive. We are moving away from a post-process centric data analysis paradigm towards a concurrent analysis framework, in which raw simulation data is processed as it is computed. Algorithms must adapt to machines with extreme concurrency, low communication bandwidth, and high memory latency, while operating within the time constraints prescribed by the simulation. Furthermore, in- put parameters are often data dependent and cannot always be prescribed. The study of sublinear algorithms is a recent development in theoretical computer science and discrete mathematics that has significant potential to provide solutions for these challenges. The approaches of sublinear algorithms address the fundamental mathematical problem of understanding global features of a data set using limited resources. These theoretical ideas align with practical challenges of in-situ and in-transit computation where vast amounts of data must be processed under severe communication and memory constraints. This report details key advancements made in applying sublinear algorithms in-situ to identify features of interest and to enable adaptive workflows over the course of a three year LDRD. Prior to this LDRD, there was no precedent in applying sublinear techniques to large-scale, physics based simulations. This project has definitively demonstrated their efficacy at mitigating high performance computing challenges and highlighted the rich potential for follow-on re- search opportunities in this space.
Direct numerical simulations of three-dimensional spatially-developing turbulent Bunsen flames were performed at three different turbulence intensities. The simulations were performed using a reduced methane-air chemical mechanism which was specifically tailored for the lean premixed conditions simulated here. A planar-jet turbulent Bunsen flame configuration was used in which turbulent preheated methane-air mixture at 0.7 equivalence ratio issued through a central jet and was surrounded by a hot laminar coflow of burned products. The turbulence characteristics at the jet inflow were selected such that combustion occured in the thin reaction zones (TRZ) regime. At the lowest turbulence intensity, the conditions fall on the boundary between the TRZ regime and the corrugated flamelet regime, and progressively moved further into the TRZ regime by increasing the turbulent intensity. The data from the three simulations was analyzed to understand the effect of turbulent stirring on the flame structure and thickness. Statistical analysis of the data showed that the thermal preheat layer of the flame was thickened due to the action of turbulence, but the reaction zone was not significantly affected. A global and local analysis of the burning velocity of the flame was performed to compare the different flames. Detailed statistical averages of the flame speed were also obtained to study the spatial dependence of displacement speed and its correlation to strain rate and curvature.
The one-dimensional turbulence (ODT) model is applied to a reactant-to-product counterflow configuration and results are compared with DNS data. The model employed herein solves conservation equations for momentum, energy, and species on a one dimensional (1D) domain corresponding to the line spanning the domain between nozzle orifice centers. The effects of turbulent mixing are modeled via a stochastic process, while the Kolmogorov and reactive length and time scales are explicitly resolved and a detailed chemical kinetic mechanism is used. Comparisons between model and DNS results for spatial mean and root-mean-square (RMS) velocity, temperature, and major and minor species profiles are shown. The ODT approach shows qualitatively and quantitatively reasonable agreement with the DNS data. Scatter plots and statistics conditioned on temperature are also compared for heat release rate and all species. ODT is able to capture the range of results depicted by DNS. However, conditional statistics show signs of underignition.
A turbulent lifted slot-jet flame is studied using direct numerical simulation (DNS). A one-step chemistry model is employed with a mixture-fraction-dependent activation energy which can reproduce qualitatively the dependence of the laminar burning rate on the equivalence ratio that is typical of hydrocarbon fuels. The basic structure of the flame base is first examined and discussed in the context of earlier experimental studies of lifted flames. Several features previously observed in experiments are noted and clarified. Some other unobserved features are also noted. Comparison with previous DNS modelling of hydrogen flames reveals significant structural differences. The statistics of flow and relative edge-flame propagation velocity components conditioned on the leading edge locations are then examined. The results show that, on average, the streamwise flame propagation and streamwise flow balance, thus demonstrating that edge-flame propagation is the basic stabilisation mechanism. Fluctuations of the edge locations and net edge velocities are, however, significant. It is demonstrated that the edges tend to move in an essentially two-dimensional (2D) elliptical pattern (laterally outwards towards the oxidiser, then upstream, then inwards towards the fuel, then downstream again). It is proposed that this is due to the passage of large eddies, as outlined in Su et al. (Combust. Flame, vol. 144 (3), 2006, pp. 494-512). However, the mechanism is not entirely 2D, and out-of-plane motion is needed to explain how flames escape the high-velocity inner region of the jet. Finally, the time-averaged structure is examined. A budget of terms in the transport equation for the product mass fraction is used to understand the stabilisation from a time-averaged perspective. The result of this analysis is found to be consistent with the instantaneous perspective. The budget reveals a fundamentally 2D structure, involving transport in both the streamwise and transverse directions, as opposed to possible mechanisms involving a dominance of either one direction of transport. It features upstream transport balanced by entrainment into richer conditions, while on the rich side, upstream turbulent transport and entrainment from leaner conditions balance the streamwise convection.
Application resilience is a key challenge that must be ad-dressed in order to realize the exascale vision. Previous work has shown that online recovery, even when done in a global manner (i.e., involving all processes), can dramatically re-duce the overhead of failures when compared to the more traditional approach of terminating the job and restarting it from the last stored checkpoint. In this paper we suggest going one step further, and explore how local recovery can be used for certain classes of applications to reduce the over-heads due to failures. Specifically we study the feasibility of local recovery for stencil-based parallel applications and we show how multiple independent failures can be masked to effectively reduce the impact on the total time to solution.
This paper reports the results of a joint experimental and numerical study of the flow characteristics and flame structure of a hydrogen rich jet injected normal to a turbulent, vitiated crossflow of lean methane combustion products. Simultaneous high-speed stereoscopic PIV and OH PLIF measurements were obtained and analyzed alongside three-dimensional direct numerical simulations of inert and reacting JICF with detailed H2/CO chemistry. Both the experiment and the simulation reveal that, contrary to most previous studies of reacting JICF stabilized in low-to-moderate temperature air crossflow, the present conditions lead to a burner-attached flame that initiates uniformly around the burner edge. Significant asymmetry is observed, however, between the reaction zones located on the windward and leeward sides of the jet, due to the substantially different scalar dissipation rates. The windward reaction zone is much thinner in the near field, while also exhibiting significantly higher local and global heat release than the much broader reaction zone found on the leeward side of the jet. The unsteady dynamics of the windward shear layer, which largely control the important jet/crossflow mixing processes in that region, are explored in order to elucidate the important flow stability implications arising in the inert and reacting JICF. The paper concludes with an analysis of the ignition, flame characteristics, and global structure of the burner-attached flame. Chemical explosive mode analysis (CEMA) shows that the entire windward shear layer, and a large region on the leeward side of the jet, are highly explosive prior to ignition and are dominated by non-premixed flame structures after ignition. The predominantly mixing limited nature of the flow after ignition is examined by computing the Takeno flame index, which shows that ~70% of the heat release occurs in non-premixed regions.