SEARCH · Search NASA

Results for “engine simulation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Zero Order Reactioin Kinetics: Enabling the use of Detailed Chemical Kinetics in Combustion Simulations (Final CRADA Report)

This was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS), as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Gamma Technologies, LLC (GT or Participante), to incorporate the ability to access LLNL chemical kinetics technologies while using GT-SUITE, GT’s market leading engine simulation software. At the end of the project, LLNL has released Zero-RK version 3.5 with zero- and one-dimensional (0-D and 1-D) solver functionality that interfaces with GT’s GT-SUITE v2023 and later releases. GT has tested its product to assure their customers that the interface can provide reduction in chemistry solution time for detailed chemistry simulations. The process has also positioned GT to easily benefit from future improvements of the Zero-RK suite of tools developed under the DOE Vehicle Technologies Office.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Commercialization of a Comprehensive Spark-Ignition Model for Automotive Engine Applications in CONVERGE CFD

The work completed under this CRADA commercializes an advanced hybrid Lagrangian-Eulerian spark-ignition (LESI) model for multi-dimensional engine simulations. In this report, LESI is coupled with combustion models (WSR, TFM, G-equation) to model spark channel elongation and kernel and flame development in a constant volume combustion vessel, and a four-valve direct-injection spark-ignition (4VDISI) engine in two operating conditions: stoichiometric E30 with no dilution, and stoichiometric PACE-20 with dilution using CONVERGE CFD solver.

42 ENGINEERING↗

DOE OSTI2021

Integration of Open-Source URBANopt and Dragonfly Energy Modeling Capabilities into Practitioner Workflows for District-Scale Planning and Design

High-performance districts and communities offer opportunities for reducing energy use, emissions, and costs, and can be instrumental in helping cities achieve their climate goals. The design of such communities requires identification of opportunities early on and their re-evaluation throughout the planning process. There is a need for energy modeling tools that connect 3D Computer-Aided Design (CAD) platforms to simulation engines, enabling detailed energy analysis of districts within the workflows and tools used by practitioners. This paper introduces the Dragonfly and URBANoptTM combined toolset that supports the creation of urban models from a range of geometry formats typically used by designers and planners, and provides an integrated pathway to simulate district-scale energy systems. The toolset is piloted by a global architecture and master planning firm to evaluate several key urban-scale technical questions for the design of a district in Chicago. The findings indicate that, while energy savings can be achieved through traditional architectural studies and enhancements to individual building efficiency, the modeling toolset helps identify additional savings and insights that can be achieved when considering district-scale energy systems. Finally, this study demonstrates how the Dragonfly/URBANopt toolset can integrate with master planning workflows, thereby enabling an iterative performance-based design process.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2020

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

DOE OSTI2024

Effect of Cyclo-Pentane Impurities on the Autoignition Reactivity and Properties of a Gasoline Surrogate Fuel

Surrogate fuels that reproduce the characteristics of full-boiling range fuels are key tools to enable numerical simulations of fuel-related processes and ensure reproducibility of experiments by eliminating batch-to-batch variability. Within the PACE initiative, a surrogate fuel for regular-grade E10 (10%vol ethanol) gasoline representative of a U.S. market gasoline, termed PACE-20, was developed and adopted as baseline fuel for the consortium. Although extensive testing demonstrated that PACE-20 replicates the properties and combustion behavior of the full-boiling range gasoline, several concerns arose regarding the purity level required for the species that compose PACE-20. This is particularly important for cyclo-pentane, since commercial-grade cyclo-pentane typically shows 60%–85% purity. Here, in the present work, the effects of the purity level of cyclo-pentane on the properties and combustion characteristics of PACE-20 were studied. Chemical kinetic simulations were performed to predict the effects of cyclo-pentane impurities on the properties, octane rating, and autoignition reactivity under homogeneous charge compression-ignition conditions of PACE-20. From the numerical results, cyclo-pentane with 85% purity or higher is required to reasonably match both the research octane number and motor octane number of the target gasoline. Finally, homogeneous charge compression-ignition engine simulations show that impurities have only a modest effect on reactivity at naturally aspirated conditions, but cyclo-pentane purity is critical to properly replicate the pressure dependency of the reactivity.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

TRUST Nonlinear Dynamics (ND) Report: FY22

The objective of the Delivery Environments (DE) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) work package, is to quantify and help increase confidence in specific areas of computational and experimental capabilities that are applicable to current and future delivery environments. More complete quantification of confidence in experimental and computational capabilities and the sufficient increase of confidence in those capabilities is critical to improving weapons engineering design, qualification, and assessment efforts that are relevant to the current and future stockpile. This work will use and provide feedback on analysis tools and experimental results databases for efficient and responsive engineering which are currently under development: engineering common model framework (ECMF), engineering quantification of margins and uncertainties (EQMU), and the test information management system (TIMS). Additionally, this work is being used to test a newly available framework in W-13, Weapons Analysis for Verified Engineering Simulation (WAVES).

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2022

aphBO-2GP-3B: a budgeted asynchronous parallel multi-acquisition functions for constrained Bayesian optimization on high-performing computing architecture

High-fidelity complex engineering simulations are often predictive, but also computationally expensive and often require substantial computational efforts. The mitigation of computational burden is usually enabled through parallelism in high-performance cluster (HPC) architecture. Optimization problems associated with these applications is a challenging problem due to the high computational cost of the high-fidelity simulations. In this paper, an asynchronous parallel constrained Bayesian optimization method is proposed to efficiently solve the computationally expensive simulation-based optimization problems on the HPC platform, with a budgeted computational resource, where the maximum number of simulations is a constant. The advantage of this method are three-fold. Firstly, the efficiency of the Bayesian optimization is improved, where multiple input locations are evaluated parallel in an asynchronous manner to accelerate the optimization convergence with respect to physical runtime. This efficiency feature is further improved so that when each of the inputs is finished, another input is queried without waiting for the whole batch to complete. Second, the proposed method can handle both known and unknown constraints. Third, the proposed method samples several acquisition functions based on their rewards using a modified GP-Hedge scheme. The proposed framework is termed aphBO-2GP-3B, which means asynchronous parallel hedge Bayesian optimization with two Gaussian processes and three batches. The numerical performance of the proposed framework aphBO-2GP-3B is comprehensively benchmarked using 16 numerical examples, compared against other 6 parallel Bayesian optimization variants and 1 parallel Monte Carlo as a baseline, and demonstrated using two real-world high-fidelity expensive industrial applications. The first engineering application is based on finite element analysis (FEA) and the second one is based on computational fluid dynamics (CFD) simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Introduction to Finite Element Analysis: Abaqus Tutorial [Slides]

Abaqus is a suite of powerful engineering simulation programs, based on the finite element method, that can solve a variety of problems, ranging from simple to complex, including the fields of: Structural mechanics, Structural dynamics, Thermodynamics, Heat transfer, and more

42 ENGINEERING↗

DOE OSTI2024

Materials genome innovation for computational software (magics) center

Functional layered material (LM) architectures will dominate nanomaterials science in this century. We have developed theory, modeling, simulation, and software and data tools that enhance understanding and AI guide synthesis, enable characterization of complex structures, and improve capabilities in the predictive design and growth of LMs. Research at the Center has focused on: Computational synthesis and characterization: AI guided synthesis and experimental synthesis of stacked LMs with tailored properties via optimized chemical vapor deposition (CVD) growth and liquid-phase exfoliation; study defects, edges, grain boundaries, wrinkling of atomic layers and their effects on chemical, mechanical, electrical, and optical properties. Far-from-equilibrium processes: Joint experimental and simulation based probe of electronic processes with NAQMD and ultrafast X-ray free-electron laser (XFEL) and ultrafast electron diffraction (UED) facilities at Stanford. Experimentally validate NAQMD by ultrafast electron diffraction and X-ray spectroscopy studies of structural and excited state dynamics, shape fluctuations, and phonon dynamics. Scalable software: Simulation engines for desktop-to-exascale platforms using low-overhead, linear-scaling QMD algorithms; divide-conquer-recombine NAQMD with electronic excitations; extended-Lagrangian reactive molecular dynamics (RMD), machine learning (ML) based neural-network quantum molecular dynamics (NNQMD), and super-state accelerated molecular dynamics (AMD) and kinetic Monte Carlo codes; thermal and electrical transport software; and design 3D architectures of LMs with desired functionality using scalable software. Distribution of software and data, and training: Software and simulation-experimental data generated within the Center are distributed to the materials science community via Berkeley Materials Project (MP) framework. We have also organized three workshops for software distribution and training at USC (Nov. 2017, Mar. 2018) and Gaithersburg, MD (Nov. 2018) to train researchers, with the last one in focused on underrepresented groups, in collaboration with Howard University which is one of the largest HBCUs. The Center supported a total of 46 personnel and 6 undergraduate students. These include 14 faculty, 11 postdoctoral research associates, 20 graduate research assistants, and mentored 6 undergraduate students. This resulted in the publications of 63 research papers that include 46 publications on Reactive and Quantum Dynamics Simulations, 13 publications on Machine Learning for Quantum Materials, and 4 publications on Quantum Computing.

2D Materials↗

DOE OSTI2021

Pyrolysis of bio-derived dioxolane fuels: A ReaxFF molecular dynamics study

Alkyl-substituted 1,3-dioxolanes, including 4,5-dimethyl-2-pentan-3-yl-1,3-dioxolane (Fuel 1), 4,5-dimethyl-2-pentyl-1,3-dioxolane (Fuel 2), and 2-(heptan-3-yl)-4,5-dimethyl-1,3-dioxolane (Fuel 3), have been recently suggested as potential biodiesels. In this paper, we investigate the initial pyrolysis of the alkyl-substituted 1,3-dioxolanes at high temperatures using ReaxFF molecular dynamics (MD) simulations. We analyze the decomposition rate, reaction mechanism, and product distribution in the pyrolysis of the three alkyl 1,3-dioxlanes. The three fuels primarily decompose to 4,5-dimethyl-1,3-dioxolane radical and hydrocarbons derived from the alkyl side-chains. The further decomposition of 4,5-dimethyl-1,3-dioxolane radical primarily leads to 2-C 4 H 8 and CO 2 within a few decomposition steps. The hydrocarbon product distribution is significantly affected by the molecular structure of the alkyl side-chain, which would have a strong influence on the sooting tendency of these fuels. The ReaxFF simulations predict that the order of sooting tendency would be Fuel 3 > Fuel 1 > Fuel 2, which agrees with the measured sooting tendency trend. Based on the pyrolysis mechanism identified by ReaxFF, we propose a new alkyl dioxolane, 4-hexyl-5-methyl-1,3-dioxolane (Fuel 4), which might produce even less soot, by modifying the molecular structure of Fuel 2. Our ReaxFF simulation shows that Fuel 4 produce much less C 4 H 8 , an effective non-aromatic soot precursor, than Fuel 2. Moreover, more carbon atoms are bonded to each oxygen atom in Fuel 4 than Fuel 2, which would help reduce soot yield by removing more carbon atoms from the soot-producing pool of species. The major decomposition pathways identified in this work can be used to develop chemical kinetic models for 1,3-dioxolane based compounds, as biodiesel components, applicable to combustion engine simulations. We also demonstrate that the chemical kinetic insight offered by ReaxFF simulations can be used to design new fuel molecules with more desired properties.

09 BIOMASS FUELS↗

DOE OSTI2023

Genetic algorithm optimization of a chemical kinetic mechanism for propane at engine relevant conditions

Propane has demonstrated significant potential for reductions in greenhouse gas and pollutant emissions in medium- and heavy-duty engine applications, but further improvements require accurate, compact, and scalable chemical kinetic mechanisms to design the next generation of propane fueled engines, particularly at the boosted operating conditions necessary to meet the power density demand of medium- and heavy-duty applications. In this work, six key chemical reactions were identified in a reduced mechanism with 70 species and 352 reactions through a sensitivity analysis performed at conditions typical of thermodynamic trajectories observed in a high compression ratio, long stroke engine operated on propane from throttled to boosted operating conditions. While the original mechanism was validated against rapid compression machine (RCM) data, it was found to overpredict experimental autoignition tendencies in 2-zone, 0-D SI engine simulations performed in Chemkin Pro. Subsequently, a genetic algorithm approach was used to optimize the six reaction rate parameters within established uncertainty bounds by performing RCM simulations and comparing to two independent sets of literature ignition delay times for propane, thus generating two new kinetic mechanisms. The first optimization achieved a mean absolute percent error (MPE) reduction in 2nd stage ignition delay of 61.4% in seven generations, while the second optimization utilized a newer experimental RCM dataset, and achieved MPE reduction of 56.7% in seven generations, and further marginal improvement to 57.8% reduction in 34 generations. Finally, the two mechanisms were then evaluated again in the 2-zone 0-D SI engine model in Chemkin Pro comparing typical mean and knocking cycle trajectories, and it was found that the second optimized mechanism provided better prediction of knock onset at the representative conditions evaluated in this work, particularly for higher load operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Programming approaches for scalability, performance, and portability of combustion physics codes

Here, this paper presents the process, strategy, and results associated with porting a typical combustion physics flow solver to current state-of-the-art and future massively-parallel computer architectures. Major focus is placed on the distinct algorithmic structure of these types of codes and how it can be integrated with modern programming paradigms for heterogeneous platforms (i.e., distributed many-core systems with accelerators). An end-to-end case study is presented that exemplifies the process in a generic manner, which then serves as a clear guide with respect to the strategy and best practices leading to a robust and adaptable framework that performs well, is durable over time, is portable, and requires minimal human-effort. This end is accomplished beginning with the use of a mature, validated, structured, multiblock code framework optimized for application of both Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS). This code has been ported to a variety of platforms over the past decade, including most recently the Oak Ridge Leadership Computing Facility’s “Summit” Platform. The experience gained on these multiple platforms provides general insights and thus the results presented are not specific to any one code or platform other than the overarching trend toward distributed many-core systems with accelerators in order to move toward exascale performance. The resultant performance and scalability of the ported code is demonstrated on a real-world application; a state-of-the-art rotating detonation rocket engine simulation that matches the complex geometry and boundary conditions imposed as part of a companion experimental campaign.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Model Development for Multi-Component Fuel Vaporization and Flash Boiling

The objectives of this project are to improve the multi-component fuel droplet and film vaporization models used in internal combustion engine simulation, and to develop a comprehensive model to predict the characteristics of multi-component flash boiling spray. This work explores two approaches to fuel composition treatment for modeling multi-component fuel vaporization: one based on discretization and surrogates, and the other based on continuous thermodynamic distribution of fuel properties. The experimental data collected for model validation are done under three fuel form factors: droplet, spray, and thin film. The study on sprays also include experimentation under non-flash and flash boiling conditions, a phenomenon that enhances fuel vaporization. The main goals of this work are: Design and develop a multi-component fuel droplet and wall film vaporization model using both discrete and continuous thermodynamics methods. Design and develop an analytical model for multi-component flash boiling. Integrate the multi-component droplet and film model into multi-dimensional engine calculations to predict the fuel vaporization process under engine operation condition. Conduct multi-component droplet and fuel film vaporization experiments in a non-combusting chamber to verify the proposed vaporization models. Characterize flash boiling phenomena of multi-component fuel sprays by optical and laser diagnostic techniques. This report will detail the experimental setup and the numerical basis for developing a model to achieve the main goals listed above. Key features of observed multi-component fuel vaporization will be summarized at the end of each experimental sections, and corresponding model performance evaluation will be presented at the end of each model development section.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

URBANopt: An Open-Source Software Development Kit for Community and Urban District Energy Modeling: Preprint

Urban building modeling tools are developing rapidly; these tools use emerging simulation workflows for specific urban environmental design tasks, such as assessing the impacts of energy efficiency technologies at a district scale. However, with the emergence of new environmental design tasks, addressing all possible use cases and tasks is challenging and cannot be covered by a single tool. Urban-scale analysis at this level of complexity often requires linking multiple emerging tools, rather than using a single tool, to adequately evaluate a variety of possible fields in urban environmental design. To achieve this, flexible platforms are needed to support multiple input formats (e.g., geometric and non-geometric building properties), enabling the mapping of such inputs to underlying simulation engines. This paper provides an overview of the open-source URBANopt Software Development Kit (SDK) for modeling high-performance buildings and energy systems at a district scale. URBANopt's flexible SDK is composed of several modules that can be customized to integrate with other tools and generate new workflows to perform urban environmental design tasks, such as capturing interactions between individual buildings, district energy systems, distributed energy resources, and the electric distribution grid. We describe the functionality of the core SDK modules in URBANopt (called Core Gem, GeoJSON Gem, and Scenario Gem) and discuss the flexibility of these modules as a means of integration with a variety of tools. We also document and demonstrate technical details of writing and combining new modules to create customized workflows. Finally, we present a case study that uses the URBANopt SDK to model a hypothetical mixed-use urban project and simulate various scenarios to meet district energy performance goals.

buildings↗

DOE OSTI2023

ResStock™ v3.2.0 [SWR-19-15 and SWR-20-07]

The ResStock™ analysis tool was built on NREL's OpenStudio® platform, and is a project geared at modeling existing residential building stocks at national, regional, or local scales with a high-degree of granularity (e.g., one physics-based simulation model for every 200 dwelling units), using the EnergyPlus® simulation engine. Information about ComStock™, a sister tool for modeling the commercial building stock, can be found here: https://www.nrel.gov/buildings/comstock.html This repository contains: Housing characteristics of the U.S. residential building stock, in the form of conditional probability distributions stored as tab-separated value (.tsv) files. Comments at the bottom of each file document data sources and assumptions for each. A library of housing characteristic "options" that translate high-level characteristic parameters into arguments for OpenStudio measures, and which are referenced by the housing characteristic .tsv files and building energy upgrades defined in project definition files Project definition files: v2.3.0 and later: buildstockbatch YML files openable in any text editor v2.2.5 and prior: Project folder openable in PAT Unit-level OpenStudio Measures for automatically constructing OpenStudio Models of each representative dwelling unit model: v3.0.0 and later: OpenStudio-HPXML Measures v2.5.0 and prior: OpenStudio Measures Higher-level OpenStudio Measures for controlling simulation inputs and outputs This repository does not contain software for running ResStock simulations, which can be found as follows: Versions 2.3.0 and later only support the use of buildstockbatch for deploying simulations on high-performance or cloud computing. Version 2.3.0 also removed separate projects for single-family detached and multifamily buildings, in lieu of a combined project_national representing the U.S. residential building stock. See the changelog for more details. Versions 2.2.5 and prior support the use of the publicly available OpenStudio-PAT software as an interface for deploying simulations on cloud computing. Read the documentation for v2.2.5.

Horowitz, Scott↗

DOE OSTI2021

Experimental validation of an organic rankine-vapor compression cooling cycle using low GWP refrigerant R1234ze(E)

There is a significant global opportunity to capture and utilize low grade waste heat to reduce fossil fuel consumption, greenhouse gas emissions, and improve energy efficiency across a wide range of industries. In this work, an advanced type of thermally activated cooling system, an organic Rankine-vapor compression cycle (ORVC) with novel heat integration strategies, was designed and tested at a relevant scale for industrial waste heat recovery (300 kW th cooling capacity). The ORVC linked an organic Rankine power cycle and a vapor compression cooling cycle using a turbine and compressor that shared a single shaft. The ORVC test facility absorbed waste heat from a liquid stream at 91 °C to simulate engine coolant in diesel generator sets, rejected heat to a glycol stream at 30 °C, and generated chilled water at 7 °C. A cooling capacity of 264 kW ± 3.5 kW was experimentally validated with a COP of 0.56 ± 0.01 during steady-state operation at the design temperatures. The thermal efficiency, accounting for pump work, of the Rankine cycle was 7.7% ± 0.22% and the COP of the vapor compression cycle was 5.25 ± 0.09. The centrifugal turbo-compressor operated at 31.5 kRPM ± 0.3 kRPM, with a turbine and compressor isentropic efficiencies of 76.7% ± 0.90% and 84.8% ± 0.54%, respectively, with near-perfect power transmission between these components. The pressure drop in the piping and heat exchangers were significantly larger than expected which had a detrimental impact on the performance of the ORVC. In addition, the condenser on the cooling cycle could not deliver the subcooling as specified from the design point modeling. Furthermore, the results from the sensitivity analysis showed that the higher condenser glycol outlet temperature had the largest impact on performance, which is consistent with other analytical models in the literature. When the ORVC simulations were updated with experimental values for isentropic efficiencies of the turbomachinery, the thermal COP was 0.66 which represents an estimate of the predicted performance if test facility limitations are overcome.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2021

Influence of NOx Chemistry on the Prediction of Natural Gas End-Gas Autoignition

Research on the influence of NOx (NO/NO2) chemistry on HC oxidation has shown that NOx plays an important role in the promotion/inhibition of HC autoignition [1-6]. Current data for n-heptane, iso-octane, ethanol, and toluene [3, 4] have shown that NO chemistry inhibits autoignition at low temperatures and concentrations. However, at medium and high temperatures and low concentrations, autoignition is promoted, decreasing the Negative Temperature Coefficient (NTC) effect. This autoignition promotion, however, becomes weaker as NO concentrations are increased. Additionally, fundamental kinetics studies have suggested that NOx chemistry also promotes the autoignition of smaller HC, such as CH4 [7-10], C2H6 [11-14], and C3H8 [15]. Additionally to the chemical kinetics studies, a growing body of evidence supports the importance of NOx chemistry on SI engine simulations [16, 17]. Foong et al. [18] suggested that the addition of NO in their models significantly improved the agreement between simulations and experiments where the onset of knock was advanced due to the presence of NO, and Morganti et al. [19] showed that NO addition in the residual-gas was necessary for their models to obtain good agreement with experimental data. Lastly, Mohr et al. [20, 21] investigated the effect of Natural Gas (NG) reactivity and EGR substitution rate and composition on homogeneous ignition delay, flame speeds, and EGAI for stoichiometric NG/oxidizer/EGR blends. They observed that the addition of EGR composed of only inert species (CO2 and Ar) suppressed homogeneous autoignition in all NG fuels tested under all conditions. However, for all NG fuels tested, increasing EGR rates with reactive species (Ar, CO2, CO, and NO) promoted homogeneous autoignition under all conditions, thus reproducing the same trends for various NG fuel compositions that were observed in other studies and fuels. For this reason, a study is being conducted to analyze the effects of NOx chemistry on the prediction of NG/air/EGR homogeneous autoignition. This study has been divided into two tasks. First, the addition of NOx chemistry to the ARIES82 mechanism [20] is being validated using homogeneous ignition delay data collected by Mohr [20, 21]. Second, once the new modified mechanism with NOx reactions has been validated, it will be employed on engine multi-dimensional modeling presented in [22, 23] to analyze the influence of NOx chemistry on NG EGAI. Results thus far have shown that homogeneous ignition delay calculations are sensitive to NOx chemistry, where the modified mechanisms captured well all trends and closely matched the homogeneous ignition delays observed by Mohr et al. Also, the results in this work suggest that NOx chemistry is necessary to correctly capture the trends in NG homogeneous autoignition when EGR with reactive species is used since the opposite behavior is observed when only inert species are considered for the EGR composition.

03 NATURAL GAS↗

DOE OSTI2021

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗