Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation fidelity”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Toward better hydraulic fracturing fluids and their application in energy production: A review of sustainable technologies and reduction of potential environmental impacts

Recent advances in hydraulic fracturing, in conjunction with horizontal drilling, have enabled large-scale extraction of natural gas and oil from shale formations. Despite its advances and enormous economic benefits, opportunities remain to increase hydraulic fracturing efficiency and minimize potential environmental impacts. This review specifically examines three key themes associated with development and utilization of hydraulic fracturing fluids: 1) characteristics and behavior of fracturing fluids, 2) understanding and predicting migration and fate of fracturing fluids, 3) technologies to reduce environmental impact of fracturing fluids. Furthermore, the paper discusses key and new techniques and findings on rheology of hydrogel-based fluids, high fidelity simulation of propagation transport, potential environmental impacts, geosynthetics in mitigating contamination, and greener fracturing fluids. It is indicated that future development relies on advances in understanding of physical processes, modeling capabilities, and monitoring techniques.

02 PETROLEUM↗
DOE OSTI2018

Discrete element modeling of deformable pinewood chips in cyclic loading test

The design of efficient lignocellulosic biomass feedstock systems is challenging, as current laboratory characterization and design methods were developed primarily for fine powders with relatively low compressibility. Here, the discrete element method (DEM) is gaining prominence as an alternative method for modeling the bulk flow and transport of particulate materials in hoppers and feeders. However, prior DEM simulations investigating the flow of wood chips modeled the particles as simple rigid geometries such as spheres, rods or blocks, and neglected the effects of particle deformability and irregular shapes. As a consequence, those simplified DEM models may not provide enough key diagnosis to help improve the design of biomass feeding and handling equipment. This work presents a bonded-sphere DEM approach for characterizing the mechanical behavior of bulk flexible, deformable pinewood chips in a cyclic stress loading test. Clustered spheres that can bend and twist via elastic bonds are used to model irregular-shape particles in real pinewood chip samples. An axial compression tester, which contains 0.06 million bonded-sphere particles (1.35 million spheres) in a quarter cylinder, is simulated with a domain size similar to the experiments. With careful calibration, the simulations have delivered the bulk densities and the bulk moduli of elasticity that are in good agreement with those measured in the corresponding experiments. However, it is also been found challenging for the present DEM model to accurately predict the overall stress-strain behavior of bulk pinewood chips, especially the large sustained plastic deformation during unloading. Additional numerical tests have shown that the adjustment in certain contact parameters (e.g. bond stiffness) can lead to profound solution improvement, but meanwhile will induce extra challenges such as increased computing time. Future work will include an elasto-plastic particle bond model to enhance the simulation fidelity.

42 ENGINEERING↗
DOE OSTI2021

Scaling of oxygen-methane reacting coaxial jets using x-ray fluorescence to measure mixture fraction

Stoichiometric mixing length L s of reacting coaxial jet flames is a critical scaling parameter for liquid rocket engine combustors. Previous studies have shown that L s for shear coaxial flames can be scaled like their non-reacting counterparts using a non-dimensional momentum flux ratio J. In addition, stoichiometric mixing lengths of reacting and nonreacting coaxial jets collapse upon a single line by altering J using an effective outer flow gas density. This effective density is calculated from a modified version of the equivalence principle, originally developed by Tacina and Dahm [1, 2] and accounts for the effects of heat release on mixing. However, previous studies also required a second nonphysical scaling constant S c for the reacting jets, which is not predicted by the equivalence principle [3]. It was originally hypothesized that S c is attributed to the limitation of hydroxyl (OH) planar laser-induced fluorescence, which only infers L s . Direct quantitative measurement of conserved scalar fields using conventional optical diagnostics is difficult due to the lack of a tracer that easily fluoresces, survives high temperature oxygen flames, and is not dominated by quenching effects. To measure a conserved scalar field, this work implements x-ray fluorescence of Kr and Ar tracers to obtain quantitative mixture fraction fields. From these mixture fraction fields, stoichiometric mixing lengths for two CH 4 /O 2 flames are calculated and scaled against nonreacting coaxial mixing lengths using the equivalence principle. By directly measuring the stoichiometric mixing length, it is established that the additional constant is a byproduct of the OH measurement technique and the equivalence principle fully captures the scaling. Finally, comparison with high-fidelity simulation of the flame further supports this conclusion. In addition to further strengthening this scaling method, this work represents the first use of x-ray fluorescence to make quantitative conserved scalar measurements in turbulent flames.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Mixing and detonation structure in a rotating detonation engine with an axial air inlet

High-fidelity simulations of an experimental rotating detonation engine with an axial air inlet were conducted. The system operated with hydrogen as fuel at globally stoichiometric conditions. Instantaneous data showed that the detonation front is highly corrugated, and is considerably weaker than an ideal Chapman–Jouguet wave. Regions of deflagration are present ahead of the wave, caused by mixing with product gases from the previous cycle, as well as the injector recovery process. It is found that as the post-detonation high pressure flow expands, the injectors recover unsteadily, leading to a transient mixing process ahead of the next cycle. The resulting flow structure not only promotes mixing between product and reactant gases, but also increases likelihood of autoignition. These results show that the detonation process is very sensitive to injector design and the transient behavior during the detonation cycle. Phase-averaged statistics and conditionally averaged data are used to understand the overall reaction structure. As a result, comparisons with available experimental data on this configuration show remarkable good agreement of the predicted reacting flow structure.

42 ENGINEERING↗
DOE OSTI2016

Reynolds averaged turbulence modelling using deep neural networks with embedded invariance

There exists significant demand for improved Reynolds-averaged Navier–Stokes (RANS) turbulence models that are informed by and can represent a richer set of turbulence physics. This paper presents a method of using deep neural networks to learn a model for the Reynolds stress anisotropy tensor from high-fidelity simulation data. A novel neural network architecture is proposed which uses a multiplicative layer with an invariant tensor basis to embed Galilean invariance into the predicted anisotropy tensor. It is demonstrated that this neural network architecture provides improved prediction accuracy compared with a generic neural network architecture that does not embed this invariance property. Furthermore, the Reynolds stress anisotropy predictions of this invariant neural network are propagated through to the velocity field for two test cases. For both test cases, significant improvement versus baseline RANS linear eddy viscosity and nonlinear eddy viscosity models is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Modeling Single-Crystal Battery Materials: From Fundamental Understanding to Performance Evaluation

The performance of rechargeable batteries is fundamentally influenced by the physicochemical properties and microstructural features of their key material components. Recent experimental advancements have highlighted the potential of single-crystal (SC) morphologies to address inherent limitations of polycrystalline (PC) electrodes and solid-state electrolytes, offering tunable charge transport kinetics and improved cell cycling performance. Here, this review examines how state-of-the-art computational modeling, from atomistic and mesoscale to continuum-level approaches, including machine learning methodologies, has been utilized to investigate the critical factors governing the electrochemical behavior of SC battery materials. We explore how predictive modeling can elucidate the processing–structure–property–performance relationships of SC cathodes, anodes, and solid-state electrolytes, with a focus on unique SC characteristics such as crystallographic anisotropy, size effects, and facet-dependent properties. Additionally, we identify limitations in commonly used modeling techniques and discuss strategies to address these challenges. By integrating high-fidelity simulations with experimental insights, this review aims to outline a clear path for the rational design and optimization of SC battery components, paving the way for accelerated advancements in energy storage technologies.

Materials science↗
DOE OSTI2025

A Computational Framework for Simulations of Dissipative Nonadiabatic Dynamics on Hybrid Oscillator-Qubit Quantum Devices

Here, we introduce a computational framework for simulating nonadiabatic vibronic dynamics on circuit quantum electrodynamics (cQED) platforms. Our approach leverages hybrid oscillator-qubit quantum hardware with midcircuit measurements and resets, enabling the incorporation of environmental effects such as dissipation and dephasing. To demonstrate its capabilities, we simulate energy transfer dynamics in a triad model of photosynthetic chromophores inspired by natural antenna systems. We specifically investigate the role of dissipation during the relaxation dynamics following photoexcitation, where electronic transitions are coupled to the evolution of quantum vibrational modes. Our results indicate that hybrid oscillator-qubit devices, operating with noise levels below the intrinsic dissipation rates of typical molecular antenna systems, can achieve the simulation fidelity required for practical computations on near-term and early fault-tolerant quantum computing platforms.

Hamiltonians↗
DOE OSTI2018

Capturing Phase Behavior of Ternary Lipid Mixtures with a Refined Martini Coarse-Grained Force Field

Whether lipid rafts are present in the membranes of living cells remains hotly disputed despite their incontrovertible existence in liposomes at 298 K. In attempts to resolve this debate, molecular dynamics (MD) simulations have been extensively used to study lipid phase separation at high resolution. However, computation has been of limited utility in this respect because the experimental distributions of phases in lamellar lipid mixtures are poorly reproduced by simulations. In particular, all-atom (AA) approaches suffer from restrictions on accessible time scales and system sizes whereas the more efficient coarse-grained (CG) force fields remain insufficiently accurate to achieve correspondence with experiment. We refine the CG Martini parameters for the high- and low-melting temperature (T m ) lipids 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) and 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC). Our approach involves the modification of bonded Martini parameters based on fitting to atomistic simulations conducted with the CHARMM36 lipid force field. The resulting CG parameters reproduce experimental structural and thermodynamic properties of homogeneous lipid membranes while concurrently improving simulation fidelity to experimental phase diagrams of DPPC, DOPC, and cholesterol lipid mixtures. Importantly, the refined parameters provide much better phase accuracy for regions near the critical point that mimic the lipid concentrations under physiological conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Screening of Therapeutic Agents for COVID-19 Using Machine Learning and Ensemble Docking Studies

The current pandemic demands a search for therapeutic agents against the novel coronavirus SARS-CoV-2. In this work, we present an efficient computational strategy that combines machine learning (ML)-based models and high-fidelity ensemble docking studies to enable rapid screening of possible therapeutic ligands. Targeting the binding affinity of molecules for either the isolated SARS-CoV-2 S-protein at its host receptor region or the S-protein:human ACE2 interface complex, we screen ligands from drug and biomolecule data sets that can potentially limit and/or disrupt the host-virus interactions. Top scoring one hundred eighty-seven ligands (with 75 approved by the Food and Drug Administration) are further validated by all atom docking studies. Important molecular descriptors ( 2 χ n , topological surface area, and ring count) and promising chemical fragments (oxolane, hydroxy, and imidazole) are identified to guide future experiments. Overall, this work expands our knowledge of small-molecule treatment against COVID-19 and provides a general screening pathway (combining quick ML models with expensive high-fidelity simulations) for targeting several chemical/biochemical problems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

The KPP Boundary Layer Scheme for the Ocean: Revisiting Its Formulation and Benchmarking One‐Dimensional Simulations Relative to LES

Abstract We evaluate the Community ocean Vertical Mixing project version of the K‐profile parameterization (KPP) for modeling upper ocean turbulent mixing. For this purpose, one‐dimensional KPP simulations are compared across a suite of oceanographically relevant regimes against horizontally averaged large eddy simulations (LESs). We find the standard configuration of KPP consistent with LES across many forcing regimes, supporting its physical basis. Our evaluation also motivates recommendations for KPP best practices within ocean circulation models and identifies areas where further research is warranted. The original treatment of KPP recommends the matching of interior diffusivities and their gradients to the KPP‐predicted values computed in the ocean surface boundary layer (OSBL). However, we find that difficulties in representing derivatives of rapidly changing diffusivities near the base of the OSBL can lead to loss of simulation fidelity. To mitigate this difficulty, we propose and evaluate two computationally simpler approaches: (1) match to the internal predicted diffusivity alone and (2) set the KPP diffusivity to 0 at the OSBL base. We find the KPP entrainment buoyancy flux to be sensitive to vertical grid resolution and details of how to diagnose the KPP boundary layer depth. We modify the KPP turbulent shear velocity parameterization to reduce resolution dependence. Additionally, an examination of LES vertical turbulent scalar flux budgets shows that the KPP‐parameterized nonlocal tracer flux is incomplete due to the assumption that it solely redistributes the surface tracer flux. This result motivates further studies of the nonlocal flux parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The GFDL‐CM4X Climate Model Hierarchy, Part II: Case Studies

This paper is Part II of a two‐part paper that documents the Climate Model version 4X (CM4X) hierarchy of coupled climate models developed at the Geophysical Fluid Dynamics Laboratory. Part I of this paper is presented in Griffies et al. (2025a, https://doi.org/10.1029/2024MS004861 ). Here we present a suite of case studies that examine ocean and sea ice features that are targeted for further research, which include sea level, eastern boundary upwelling, Arctic and Southern Ocean sea ice, Southern Ocean circulation, and North Atlantic circulation. The case studies are based on experiments that follow the protocol of version 6 from the Coupled Model Intercomparison Project. The analysis reveals a systematic improvement in the simulation fidelity of CM4X relative to its CM4.0 predecessor, as well as an improvement when refining the ocean/sea ice horizontal grid spacing from the 0.25° of CM4X‐p25 to the 0.125° of CM4X‐p125. Even so, there remain many outstanding biases, thus pointing to the need for further grid refinements, enhancements to numerical methods, and/or advances in parameterizations, each of which target long‐standing model biases and limitations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE

On the basis of several recent breakthroughs in fusion research, many activities have been launched around the world to develop fusion power plants on the fastest possible time scale. In this context, high-fidelity simulations of the plasma behavior on large supercomputers provide one of the main pathways to accelerating progress by guiding crucial design decisions. When it comes to determining the energy confinement time of a magnetic confinement fusion device, which is a key quantity of interest, gyrokinetic turbulence simulations are considered the approach of choice – but the question, whether they are really able to reliably predict the plasma behavior is still open. The present study addresses this important issue by means of careful comparisons between state-of-the-art gyrokinetic turbulence simulations with the GENE code and experimental observations in the ASDEX Upgrade tokamak for an unprecedented number of simultaneous plasma observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A robust synthetic data generation framework for machine learning in high-resolution transmission electron microscopy (HRTEM)

Machine learning techniques are attractive options for developing highly-accurate analysis tools for nanomaterials characterization, including high-resolution transmission electron microscopy (HRTEM). However, successfully implementing such machine learning tools can be difficult due to the challenges in procuring sufficiently large, high-quality training datasets from experiments. In this work, we introduce Construction Zone, a Python package for rapid generation of complex nanoscale atomic structures which enables fast, systematic sampling of realistic nanomaterial structures and can be used as a random structure generator for large, diverse synthetic datasets. Using Construction Zone, we develop an end-to-end machine learning workflow for training neural network models to analyze experimental atomic resolution HRTEM images on the task of nanoparticle image segmentation purely with simulated databases. Further, we study the data curation process to understand how various aspects of the curated simulated data—including simulation fidelity, the distribution of atomic structures, and the distribution of imaging conditions—affect model performance across three benchmark experimental HRTEM image datasets. Using our workflow, we are able to achieve state-of-the-art segmentation performance on these experimental benchmarks and, further, we discuss robust strategies for consistently achieving high performance with machine learning in experimental settings using purely synthetic data. Construction Zone and its documentation are available at https://github.com/lerandc/construction_zone.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-driven discovery of dynamics from time-resolved coherent scattering

Coherent X-ray scattering (CXS) techniques are capable of interrogating dynamics of nano- to mesoscale materials systems at time scales spanning several orders of magnitude. However, obtaining accurate theoretical descriptions of complex dynamics is often limited by one or more factors—the ability to visualize dynamics in real space, computational cost of high-fidelity simulations, and effectiveness of approximate or phenomenological models. In this work, we develop a data-driven framework to uncover mechanistic models of dynamics directly from time-resolved CXS measurements without solving the phase reconstruction problem for the entire time series of diffraction patterns. Our approach uses neural differential equations to parameterize unknown real-space dynamics and implements a computational scattering forward model to relate real-space predictions to reciprocal-space observations. This method is shown to recover the dynamics of several computational model systems under various simulated conditions of measurement resolution and noise. Moreover, the trained model enables estimation of long-term dynamics well beyond the maximum observation time, which can be used to inform and refine experimental parameters in practice. Finally, we demonstrate an experimental proof-of-concept by applying our framework to recover the probe trajectory from a ptychographic scan. Our proposed framework bridges the wide existing gap between approximate models and complex data.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Evaluation of machine learning algorithms for prediction of regions of high Reynolds averaged Navier Stokes uncertainty

Reynolds Averaged Navier Stokes (RANS) models are widely used in industry to predict fluid flows, despite their acknowledged deficiencies. Not only do RANS models often produce inaccurate flow predictions, but there are very limited diagnostics available to assess RANS accuracy for a given flow configuration. If experimental or higher fidelity simulation results are not available for RANS validation, there is no reliable method to evaluate RANS accuracy. This paper explores the potential of utilizing machine learning algorithms to identify regions of high RANS uncertainty. Three different machine learning algorithms were evaluated: support vector machines, Adaboost decision trees, and random forests. The algorithms were trained on a database of canonical flow configurations for which validated direct numerical simulation or large eddy simulation results were available, and were used to classify RANS results on a point-by-point basis as having either high or low uncertainty, based on the breakdown of specific RANS modeling assumptions. Classifiers were developed for three different basic RANS eddy viscosity model assumptions: the isotropy of the eddy viscosity, the linearity of the Boussinesq hypothesis, and the non-negativity of the eddy viscosity. It is shown that these classifiers are able to generalize to flows substantially different from those on which they were trained. As a result, feature selection techniques, model evaluation, and extrapolation detection are discussed in the context of turbulence modeling applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Data-driven prediction of scaling and ignition of inertial confinement fusion experiments

Recent advances in inertial confinement fusion (ICF) at the National Ignition Facility (NIF), including ignition and energy gain, are enabled by a close coupling between experiments and high-fidelity simulations. Neither simulations nor experiments can fully constrain the behavior of ICF implosions on their own, meaning pre- and postshot simulation studies must incorporate experimental data to be reliable. Linking past data with simulations to make predictions for upcoming designs and quantifying the uncertainty in those predictions has been an ongoing challenge in ICF research. We have developed a data-driven approach to prediction and uncertainty quantification that combines large ensembles of simulations with Bayesian inference and deep learning. The approach builds a predictive model for the statistical distribution of key performance parameters, which is jointly informed by past experiments and physics simulations. The prediction distribution captures the impact of experimental uncertainty, expert priors, design changes, and shot-to-shot variations. We have used this new capability to predict a 10× increase in ignition probability between Hybrid-E shots driven with 2.05 MJ compared to 1.9 MJ, and validated our predictions against subsequent experiments. We describe our new Bayesian postshot and prediction capabilities, discuss their application to NIF ignition and validate the results, and finally investigate the impact of data sparsity on our prediction results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Latent space mapping: Revolutionizing predictive models for divertor plasma detachment control

The inherent complexity of boundary plasma, characterized by multi-scale and multi-physics challenges, has historically restricted high-fidelity simulations to scientific research due to their intensive computational demands. Consequently, routine applications such as discharge control and scenario development have relied on faster but less accurate empirical methods. This work introduces DivControlNN, a novel machine-learning-based surrogate model designed to address these limitations by enabling quasi-real-time predictions (i.e., ~ 0.2 ms) of boundary and divertor plasma behavior. Trained on over 70,000 2D UEDGE simulations from KSTAR tokamak equilibria, DivControlNN employs latent space mapping to efficiently represent complex divertor plasma states, achieving a computational speed-up of over 10 8 compared to traditional simulations while maintaining a relative error below 20% for key plasma property predictions. During the 2024 KSTAR experimental campaign, a prototype detachment control system powered by DivControlNN successfully demonstrated detachment control on its first attempt, even for a new tungsten divertor configuration and without any fine-tuning. These results highlight the transformative potential of DivControlNN in overcoming diagnostic challenges in future fusion reactors by providing fast, robust, and reliable predictions for advanced integrated control systems.

Artificial neural networks↗