In various applications with recirculation, liquid droplets can be immersed in gases that may have a wide range of possible compositions, from pure air to hot combustion products. In order to gain fundamental understanding of the behaviour of individual droplets in vitiated air, numerical simulations of kerosene single droplet evaporation, autoignition, and combustion in conditions relevant to gas turbines have been performed. The droplet autoignition behaviour has been analysed in both physical and mixture fraction space for a wide range of vitiated air compositions and initial droplet diameters. Results show that the autoignition time delay decreases with increasing level of dilution with hot combustion products and decreasing initial droplet diameter. Chemistry is confined up to a radius of almost 10 initial droplet diameters and the location of autoignition is influenced by both the initial droplet diameter and the level of dilution. The time evolution of species in the gaseous phase after autoignition shows similar trends for all the diameters and dilution levels investigated here with the peak of temperature and OH mass fraction moving towards the droplet surface as a consequence of the balance between fuel production and consumption. In mixture fraction space, the location of the peaks of temperature and OH mass fraction after autoignition do not change in time whereas other intermediate species such as CH 2 O and pyrolysis products still exhibit a quite variable behaviour. The long-time flame structure has been compared with gaseous laminar counterflow simulations and, although qualitatively similar, the flame structure in the two configurations has differences with implications for flamelet combustion models used in spray combustion. The droplet evaporation, autoignition, and combustion behaviour has been summarized through a regime diagram showing the evaporation and autoignition time delays as a function of both initial droplet diameter and vitiated air dilution. This allows the identification of different states in the droplet combustion scenario and the introduction of critical values of dilution and initial droplet diameter beyond which single droplet rather than cloud combustion can occur, which can be exploited in the design of lean burn gas turbine combustion systems. In conclusion, the approach presented here can be easily extended to other conditions and fuels allowing the generation of regime diagrams for various operating conditions.

42 ENGINEERING↗

DOE OSTI2017

Scalability of Several Asynchronous Many-Task Models for In Situ Statistical Analysis.

This report is a sequel to [PB16], in which we provided a first progress report on research and development towards a scalable, asynchronous many-task, in situ statistical analysis engine using the Legion runtime system. This earlier work included a prototype implementation of a proposed solution, using a proxy mini-application as a surrogate for a full-scale scientific simulation code. The first scalability studies were conducted with the above on modestly-sized experimental clusters. In contrast, in the current work we have integrated our in situ analysis engines with a full-size scientific application (S3D, using the Legion-SPMD model), and have conducted nu- merical tests on the largest computational platform currently available for DOE science ap- plications. We also provide details regarding the design and development of a light-weight asynchronous collectives library. We describe how this library is utilized within our SPMD- Legion S3D workflow, and compare the data aggregation technique deployed herein to the approach taken within our previous work.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

A Unified Data-Driven Approach for Programming In Situ Analysis and Visualization: An Interim Report of Sandia Sub-Team Contributions

As we look ahead to next generation high performance computing platforms, the placement and movement of data is becoming the key-limiting factor on both performance and energy efficiency. Furthermore, the increased quantities of data the systems are capable of generating, in conjunction with the insufficient rate of improvements in the supporting I/0 infrastructure, is forcing applications away from the off-line post-processing of data towards techniques based on in ,situ analysis and visualization. Together, these challenges are shaping how we will both design and develop effective, performant and energy-efficient software. In particular, the challenges highlight the need for data and data-centric operations to be fundamental in the reasoning about, and optimization of, scientific workflows on extreme-scale architectures.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

A priori and a posteriori investigations for developing large eddy simulations of multi-species turbulent mixing under high-pressure conditions

A Direct Numerical Simulation (DNS) database was created representing mixing of species under high-pressure conditions. The configuration considered is that of a temporally evolving mixing layer. The database was examined and analyzed for the purpose of modeling some of the unclosed terms that appear in the Large Eddy Simulation (LES) equations. Several metrics are used to understand the LES modeling requirements. First, a statistical analysis of the DNS-database large-scale flow structures was performed to provide a metric for probing the accuracy of the proposed LES models as the flow fields obtained from accurate LESs should contain structures of morphology statistically similar to those observed in the filtered-and-coarsened DNS (FC-DNS) fields. To characterize the morphology of the large-scales structures, the Minkowski functionals of the iso-surfaces were evaluated for two different fields: the second-invariant of the rate of deformation tensor and the irreversible entropy production rate. To remove the presence of the small flow scales, both of these fields were computed using the FC-DNS solutions. It was found that the large-scale structures of the irreversible entropy production rate exhibit higher morphological complexity than those of the second invariant of the rate of deformation tensor, indicating that the burden of modeling will be on recovering the thermodynamic fields. Second, to evaluate the physical effects which must be modeled at the subfilter scale, an a priori analysis was conducted. This a priori analysis, conducted in the coarse-grid LES regime, revealed that standard closures for the filtered pressure, the filtered heat flux, and the filtered species mass fluxes, in which a filtered function of a variable is equal to the function of the filtered variable, may no longer be valid for the high-pressure flows considered in this study. The terms requiring modeling are the filtered pressure, the filtered heat flux, the filtered pressure work, and the filtered species mass fluxes. Improved models were developed based on a scale-similarity approach and were found to perform considerably better than the classical ones. These improved models were also assessed in an a posteriori study. Different combinations of the standard models and the improved ones were tested. At the relatively small Reynolds numbers achievable in DNS and at the relatively small filter widths used here, the standard models for the filtered pressure, the filtered heat flux, and the filtered species fluxes were found to yield accurate results for the morphology of the large-scale structures present in the flow. Analysis of the temporal evolution of several volume-averaged quantities representative of the mixing layer growth, and of the cross-stream variation of homogeneous-plane averages and second-order correlations, as well as of visualizations, indicated that the models performed equivalently for the conditions of the simulations. The expectation is that at the much larger Reynolds numbers and much larger filter widths used in practical applications, the improved models will have much more accurate performance than the standard one.

42 ENGINEERING↗

DOE OSTI2015

Irreversible entropy production rate in high-pressure turbulent reactive flows

Not Available

Borghesi, Giulio↗