Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation and analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Pilgrim Hot Springs: GEOPHIRES Inputs and Outputs for Direct-Use Geothermal District Heating and Cooling

This dataset includes files for a techno-economic analysis conducted using the GEOPHIRES simulator to examine the feasibility of expanding a larger district heating site in a remote location: Pilgrim Hot Springs, Alaska. Files included here are GEOPHIRES inputs and outputs for five different scenarios with varying demand, cycle, and system design characteristics to analyze. Also included is the link to the GEOPHIRES GitHub, as well as a link to the dataset that contains the energy modelling used to determine the heating demand for the district. For a list of the differences between scenarios, see the included "Input Overview.txt" file. Fields included in the input files are: subsurface technical parameters, surface technical parameters, financial parameters, capital and O&M parameters, as well as simulation parameters. The output files are case reports that summarize all equipment, reservoir characteristics, costs, and heating profiles.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Lambda-PFLOTRAN 1.0: a workflow for incorporating organic matter chemistry informed by ultra high resolution mass spectrometry into biogeochemical modeling

Abstract. Organic matter (OM) composition plays a central role in microbial respiration of dissolved organic matter and subsequent biogeochemical reactions. Here, a direct connection of organic matter chemistry and thermodynamics to reactive transport simulators has been achieved through the newly developed Lambda-PFLOTRAN workflow tool that succinctly incorporates carbon chemistry data generated from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) into reaction networks to simulate organic matter degradation and the resulting biogeochemistry. Lambda-PFLOTRAN is a Python-based workflow, executed through a Jupyter notebook interface, that digests raw FTICR-MS data, develops a representative reaction network based on substrate-explicit thermodynamic modeling (also termed lambda modeling due to its key thermodynamic parameter λ used therein), and completes a biogeochemical simulation with the open source, reactive flow and transport code PFLOTRAN. The workflow consists of the following five steps: configuration, thermodynamic (lambda) analysis, sensitivity analysis, parameter estimation, and simulation output and visualization. Two test cases are provided to demonstrate the functionality of the Lambda-PFLOTRAN workflow. The first test case uses laboratory incubation data of temporal oxygen depletion to fit lambda parameters (i.e., maximum utilization rate and microbial carrying capacity). A slightly more complex second test case fits multiple lambda formulation and soil organic matter release parameters to temporal greenhouse gas generation measured during a soil incubation. Overall, the Lambda-PFLOTRAN workflow facilitates upscaling by using molecular-scale characterization to inform biogeochemical processes occurring at larger scales.

58 GEOSCIENCES↗
DOE OSTI2025

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the pebble-bed high-temperature gas-cooled reactor (PB-HTGR) reference plant models for the United States Nuclear Regulatory Commission. These models serve as the foundation for the future development of detailed design evaluation models based on license applications. The reference plant models have been developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, which is the code suite proposed for non-light-water reactor systems safety analysis. It incorporates various simulation tools developed by the Nuclear Energy Advanced Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of the work that was performed was to assess BlueCRAB’s level of readiness for modeling a PB-HTGR. To do so, we first developed numerical models in BlueCRAB that include the key physics for this technology to ensure an adequate level of fidelity for modeling PB-HTGR core performance and for performing multiphysics simulations for equilibrium core conditions and different accident scenarios. Then we simulated transient scenarios, including depressurized and pressurized loss of forced cooling accidents, over-cooling, and control rod withdrawal events with delayed and prompt supercritical reactivity insertions. The analysis in this report includes comparisons of the 2D thermal fluid porous media models in Pronghorn and SAM, and comparisons of coupled SAM/Griffin/SAM and coupled Pronghorn/Griffin for depressurized and pressurized loss of forced cooling, over-cooling, and control rod withdrawal events. In addition, we compare 3D, 2D, and 0D/PKE neutronic models for the two control rod withdrawal scenarios with coupled Pronghorn/Griffin. The comparisons show that the BlueCRAB models lead to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields, such as neutron flux, power, temperature, density, pressure, and velocity, are within the expected ranges, and their distributions can be explained by the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. This work further highlights the need for flexible tools with various levels of fidelity to cover the breadth and depth of needs that may arise in future technical evaluations of the PB-HTGR. We believe that the BlueCRAB capabilities will be a significant asset for confirmatory analyses in order to resolve important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Topological Data Analysis for Particulate Gels

Soft gels, formed via the self-assembly of particulate materials, exhibit intricate multiscale structures that provide them with flexibility and resilience when subjected to external stresses. Here, this work combines particle simulations and topological data analysis (TDA) to characterize the complex multiscale structure of soft gels. Our TDA analysis focuses on the use of the Euler characteristic, which is an interpretable and computationally scalable topological descriptor that is combined with filtration operations to obtain information on the geometric (local) and topological (global) structure of soft gels. We reduce the topological information obtained with TDA using principal component analysis (PCA) and show that this provides an informative low-dimensional representation of the gel structure. We use the proposed computational framework to investigate the influence of gel preparation (e.g., quench rate, volume fraction) on soft gel structure and to explore dynamic deformations that emerge under oscillatory shear in various response regimes (linear, nonlinear, and flow). Our analysis provides evidence of the existence of hierarchical structures in soft gels, which are not easily identifiable otherwise. Moreover, our analysis reveals direct correlations between topological changes of the gel structure under deformation and mechanical phenomena distinctive of gel materials, such as stiffening and yielding. In summary, we show that TDA facilitates the mathematical representation, quantification, and analysis of soft gel structures, extending traditional network analysis methods to capture both local and global organization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The impact of cloud microphysics and ice nucleation on Southern Ocean clouds assessed with single-column modeling and instrument simulators

Abstract. Supercooled liquid clouds are common at higher latitudes (especially over the Southern Ocean) and are critical for constraining climate projections. We take advantage of the Macquarie Island Cloud and Radiation Experiment (MICRE) to perform an analysis of observed and simulated cloud processes over the Southern Ocean in a region and season dominated by supercooled liquid clouds. Using a single-column version of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS), we compare two different cloud microphysical schemes to ground-based observations of cloud, precipitation, and radiation over a 2.5-month period (1 January–17 March 2017). Both schemes are able to reproduce aspects of the cloud and radiation observations during MICRE to within the uncertainty of the data when the thermodynamic profile is prescribed with relaxation. There are differences in water mass and representation of reflectivity between the schemes. A sensitivity study of the cloud microphysics schemes, one a bulk one-moment scheme and the other a two-moment scheme with prediction of mass and number, indicates that several key processes create differences between the schemes. Surface radiative fluxes and total water path are highly sensitive to the formation and fall speed of precipitation. The prediction of hydrometeor number with the two-moment scheme yields a better comparison with observed reflectivity and radiative fluxes, despite predicting higher liquid water contents than observed. With the two-moment scheme, we are also able to test the sensitivity of the results to the input of liquid cloud condensation nuclei (CCN) and ice nuclei (IN). The cloud properties and resulting radiative effects are found to be sensitive to the CCN and IN concentrations. More CCN and IN increase liquid and ice water paths, respectively. Thus, both the dynamic environment and aerosols, integrated through the cloud microphysics, are important for properly representing Southern Ocean cloud radiative effects.

58 GEOSCIENCES↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2025

Jet Drift in Heavy-Ion Collision: Acoplanarity and v 2

We introduce a sub-eikonal anisotropic contribution to jet-broadening, “jet drift,” that couples to the flow of the nuclear medium, showing that this effect results in a deflection of hard partons, and thus jets, in the direction of the medium flow. We study two-jet observables, i.e. v 2 and acoplanarity for √s = 5.02 TeV PbPb collisions at the LHC. We show that jet drift leads to an enhancement of the elliptic flow (v 2 ) of charged hadrons and modification of the acoplanarity of dihadrons and discuss the implications for jet substructure and medium response effects. We also demonstrate that we can obtain independent information about the event geometry encoded in the drift modification to acoplanarity vs v 2 enhancement, whereas the temperature dependence of these observables is qualitatively different. By disentangling this temperature vs geometry dependence of these observables, we can add discriminatory power to our analysis. Next, we study acoplanarity and v 2 across different system sizes and collisional energies i.e. 200 GeV AuAu collision at RHIC using our in-house simulation package (APE). Our analysis will be critical for understanding and generalizing dynamical jet-QGP interactions for different collisional systems and energies which will aid event engineering applications in the future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. This study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in Open-FOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS simulation includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

Ge, Wenjun↗
DOE OSTI2024

SCALE 6.3 Validation: Reactor Physics

This study was performed to validate the SCALE/Polaris v6.3.0–PARCS v3.4.2 code procedure with the ENDF/B-VII.1 AMPX 56-group library for light-water reactor analysis by comparing the simulated results with the measured data for critical experiments and operating light-water reactors. Uncertainties of the SCALE/Polaris–PARCS code procedure for light-water reactor physics analysis were evaluated in the validation for key nuclear parameters such as reactivity, control bank work, temperature coefficients, and pin and assembly power peaking factors. In addition, the SCALE/TRITON v6.3.1 procedure with the ENDF/B-VII.1 and VIII.0 252-group and continuous-energy cross sections was validated for non-lightwater reactors including the HTR-10 reactor, the High-Temperature Test Reactor, the Molten Salt Reactor Experiment, and the Experimental Breeder Reactor II.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Global 3D Simulations of the Tokamak Edge Region

The goal of this project is to better understand magnetized plasma dynamics at the edge region of the tokamak, including both closed-flux-surface region and the scrape-off-layer (SOL), through analytical analysis and advanced numerical simulations. Under the auspices of the U.S. DOE, the PI, along with students and collaborators have made a significant progress. These results were disseminated by 13 peer-reviewed article and many presentations and posters in domestic and international conferences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Failure Analysis–Informed Risk Assessment Framework for Geological Carbon Storage Using Numerical Simulation and Machine Learning

Geological carbon storage (GCS) is recognized as a critical technology for achieving large-scale reductions in anthropogenic carbon dioxide (CO 2 ) emissions. Ensuring long-term containment and safety requires robust risk assessment frameworks that account for geological uncertainty and identify potential failure scenarios. Among various indicators, the area of review (AoR) serves as a key metric for evaluating storage performance, regulatory compliance, and monitoring design, as it delineates the spatial extent impacted by pressure buildup and plume migration. However, conventional AoR-based risk assessments typically perturb parameters within narrow uncertainty bounds, potentially overlooking rare but high-impact events arising from extreme geological conditions. In this study, we present a failure analysis–informed risk assessment framework for large-scale GCS projects to improve site prescreening and monitoring design. A suite of 300 numerical simulations was generated using stochastic geological models that vary five key parameters: net-to-gross ratio, anisotropy azimuth, porosity multiplier, permeability multiplier, and vertical-to-horizontal permeability ratio. Among these, 200 realizations represent normal geological uncertainty, while 100 additional cases explore extreme yet plausible conditions for failure-case analysis. The AoR was simulated and computed from pressure and CO 2 saturation fields, where the baseline AoR boundary, representing the extent predicted under typical geological uncertainty, was defined as the union of 200 normal-range simulations, and failure was identified when extreme-range cases exceeded this baseline. Results show that incorporating broader parameter uncertainty produces significantly larger AoR extents, underscoring the potential underestimation of risk under conventional uncertainty ranges. Furthermore, spatial probability maps derived from failure-induced AoR exceedance identify regions requiring enhanced monitoring attention. Various machine learning (ML)–based classifiers were developed to predict failure occurrence from geological parameters, with the random forest model achieving the highest performance (F1-score of 0.986). Consistent findings from correlation coefficient, feature importance, and Sobol sensitivity analyses reveal that low net-to-gross ratios and permeability multipliers are the dominant risk drivers, reflecting reduced reservoir connectivity and limited pressure dissipation. Altogether, these results provide a novel framework for risk-informed site prescreening and monitoring design that explicitly considers rare but high-impact geological scenarios in GCS projects.

25 ENERGY STORAGE↗
DOE OSTI2024

The power reliability event simulator tool (PRESTO): A novel approach to distribution system reliability analysis and applications

The growing interest in onsite solar photovoltaic and energy storage systems is partially motivated by customer concerns regarding grid reliability. However, accurately assessing the effectiveness of PVESS in mitigating these interruptions requires a comprehensive understanding of location-specific outage patterns and the ability to simulate realistic scenarios. To address the gap, we introduce the Power Reliability Event Simulation TOol (PRESTO), the first publicly available tool that simulates location-specific power interruptions at the county level. PRESTO allows for a more realistic assessment of system reliability by considering the unpredictability and location-specific patterns of power interruptions. We applied PRESTO in a case study of a single-family home across three U.S. counties, examining the performance of a solar photovoltaic system with 10kWh of battery storage during short-duration power interruptions. Our findings show that this system reliably met 93% of energy demand for essential non-heating and cooling loads, fully serving these loads in 84% of events, despite the constraints of daily time-of-use bill management which limits the battery's state-of-charge reserve. However, when heating and cooling loads were included, system performance decreased significantly, with only 70% of demand met and full service in 43% of events. These results highlight the challenges of using solar photovoltaic and energy storage systems for short-duration outages, emphasizing the need to consider factors like battery size and grid charging strategies to improve reliability. Our study demonstrates the practical applications of PRESTO, providing valuable insights into potential mitigation strategies including grid charging and optimizing battery size.

14 SOLAR ENERGY↗
DOE OSTI2024

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗
DOE OSTI2024

Exponential Improvements in the Simulation of Lattice Gauge Theories Using Near-Optimal Techniques

We report a first-of-its-kind analysis on post-Trotter simulation of U(1), SU(2), and SU(3) lattice gauge theories including fermions in arbitrary spatial dimension. We provide explicit circuit constructions as well as T-gate counts and logical qubit counts for Hamiltonian simulation. We find a reduction of up to 25 orders of magnitude in space-time volume over Trotter methods for simulations of non-Abelian lattice gauge theories relevant to the standard model. This improvement results from our algorithm having polynomial scaling with the number of colors in the gauge theory, achieved by utilizing oracle constructions relying on the sparsity of physical operators, in contrast to the exponential scaling seen in state-of-the-art Trotter methods, which employ explicit mappings onto Pauli operators. Our work demonstrates that the use of advanced algorithmic techniques leads to dramatic reductions in the cost of simulating fundamental interactions, bringing it in step with resources required for first-principles quantum simulation of chemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗