SEARCH · Engineering Papers

Results for “Simulation framework”

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 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

DOE OSTI2022

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Photoelectrode Durability in Two- versus Three-Electrode Configurations: Understanding the Impact of Circuit Configuration on Water-Splitting Stability

Device durability remains a significant challenge in photoelectrochemical (PEC) water splitting under ambient conditions. Yet, a lack of understanding of the test configuration and applied bias effects continue to hinder progress. In this study, we differentiate two-electrode (2E) and three-electrode (3E) configurations for evaluating PEC material durability, focusing particularly on their impacts on photoabsorber solid-state operating conditions. Our results underscore the fallacy of inferring 2E device stability from durability measurements performed solely in 3E configurations. Unmeasured and often misunderstood total circuit bias in 3E tests moderates material degradation, leading to the overestimation of photoelectrode stability compared to short-circuit operation. We demonstrate how the photoabsorber's operating voltage critically governs charge separation, surface stability, and degradation mechanisms during PEC operation. With these findings, we propose a standardized framework for conducting more reliable 3E durability experiments that simulate unassisted performance to help accelerate the development of robust, stable materials for solar-driven water splitting.

08 HYDROGEN↗

DOE OSTI2023

Computational and Experimental Characterization of the Ligand Environment of a Ni-Oxo Catalyst Supported in the Metal–Organic Framework NU-1000

Heterogeneous catalysts exhibit significant changes in composition due to the influence of operating conditions, and these compositional changes can have dramatic effects on catalytic performance. For traditional bulk metal heterogeneous catalysts, relationships between composition and catalytic operating conditions are well documented. However, the influence of operating conditions on the compositions of single-site heterogeneous catalysts remains largely unresolved. To address this, we report a combined computational and experimental characterization of a Ni oxo catalyst under catalytic hydrogenation conditions. Specifically, pair distribution function (PDF) analysis is combined with ab initio thermodynamic modeling to investigate ligand environments present on a Ni oxo cluster supported in the metal–organic framework NU-1000. Comparisons of the experimentally observed and simulated Ni–O coordination numbers and Ni–O, Ni···Ni, and Ni···Zr distances provide insight into the Ni ligand environment under H 2 (g). These comparisons suggest significant OH and H 2 O content and, further, that different Ni ions within the cluster and/or NU-1000 structure may comprise subtly different numbers of these ligands. Further, the observation of significant H 2 O content under H 2 (g) suggests that the NU-1000 support supplies H 2 O to the cluster. Examples of ligand environments that could lead to the observed PDFs are provided. Furthermore, the combination of simulations and experiments provides new insights into the ligand environment for Ni-NU-1000 catalysts that will be useful for understanding the ligand environments of other single-site Ni catalysts as well.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Evaluation and Optimization of Snow Albedo Scheme in Noah-MP Land Surface Model Using In Situ Spectral Observations in the Colorado Rockies

The Biosphere-Atmosphere Transfer Scheme (BATS) ground snow albedo algorithm is commonly used in land-surface models (LSM), weather forecasting and research applications. This study addresses key uncertainties in BATS simulated ground snow albedo within the Noah-MP LSM framework through evaluation and optimization of the Noah-MP BATS ground snow albedo formulation using 2-band (visible and near-infrared (NIR)) in situ albedo observations at Rocky Mountain field stations. The Noah-MP BATS ground snow albedo scheme is extremely sensitive to its input parameters. Namely, an ensemble generated by varying BATS input parameters within potentially plausible ranges provides an average daily range (maximum ensemble member minus minimum ensemble member) of ground snow albedo exceeding 0.45 in visible and NIR bands. Parameter optimization improves agreement between simulated and in situ observed ground snow albedo in visible, NIR and broadband spectrums. Importantly, optimized parameters result in reduced biases relative to observed fresh-snow albedo and better agreement with observed albedo decay. Our analysis across different sites supports that the optimized BATS ground snow albedo parameters are appropriate to transfer in space and time, at least within the region studied (the central-southern Rocky Mountains). The primary error source remaining after parameter optimization is that observed fresh-snow albedo is highly variable, particularly in the NIR spectrum, whereas BATS fresh-snow albedo is constant, an issue which requires further investigation. This study shows significant correlations between observed fresh-snow albedo and surface meteorological conditions (e.g., downward shortwave radiation and temperature) which can support future model development that attempts to include a time-varying formulation for fresh-snow albedo.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Decoding defect statistics from diffractograms via machine learning

Abstract Diffraction techniques can powerfully and nondestructively probe materials while maintaining high resolution in both space and time. Unfortunately, these characterizations have been limited and sometimes even erroneous due to the difficulty of decoding the desired material information from features of the diffractograms. Currently, these features are identified non-comprehensively via human intuition, so the resulting models can only predict a subset of the available structural information. In the present work we show (i) how to compute machine-identified features that fully summarize a diffractogram and (ii) how to employ machine learning to reliably connect these features to an expanded set of structural statistics. To exemplify this framework, we assessed virtual electron diffractograms generated from atomistic simulations of irradiated copper. When based on machine-identified features rather than human-identified features, our machine-learning model not only predicted one-point statistics (i.e. density) but also a two-point statistic (i.e. spatial distribution) of the defect population. Hence, this work demonstrates that machine-learning models that input machine-identified features significantly advance the state of the art for accurately and robustly decoding diffractograms.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Quantitative x ray phase contrast imaging of oblique shock wave–interface interactions

Oblique shock wave–interface interactions of gases and liquids have been extensively studied in shock tubes using optical methods to measure equation-of-state (EOS) parameters. However, this is difficult with solids due to their opaqueness to visible light. X ray phase contrast imaging (XPCI) has the penetrative strength to probe solids while still being sensitive to mass density and enhancing the visibility of material boundaries. We investigate the accuracy and repeatability of measuring the mean value of the average mass density (areal density divided by thickness) over region S (B S ) and flow deflection angle (θ) from XPCI images of a sample. To that end, a Hough transform-based method for measuring θ is developed. To measure B S , the XPCI image intensity probability density function (PDF) is modeled accounting for the spatial distribution of x ray energy, scintillator response, and pulse-to-pulse variation in the x ray intensity. In addition, a Monte Carlo-based algorithm for computing the B S PDF is developed. Both methods are validated on an impact-generated oblique shock wave interacting at a solid polymer-to-polymer interface. Furthermore this is accomplished through a comparison to hydrodynamic simulations using well-established EOS. Under the modeling framework for the XPCI image intensity, B S is computed with an accuracy of < 0.1% and precision of 3%–5%, while θ has an uncertainty of 0.2°, respectively. This shows that the XPCI-based model that is developed here could be an invaluable tool for high-fidelity testing of hydrodynamic models in shock polar configurations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

A biology-informed similarity metric for simulated patches of human cell membrane

Complex scientific inquiries rely increasingly upon large and autonomous multiscale simulation campaigns, which fundamentally require similarity metrics to quantify ‘sufficient’ changes among data and/or configurations. However, subject matter experts are often unable to articulate similarity precisely or in terms of well-formulated definitions, especially when new hypotheses are to be explored, making it challenging to design a meaningful metric. Furthermore, the key to practical usefulness of such metrics to enable autonomous simulations lies in in situ inference, which requires generalization to possibly substantial distributional shifts in unseen, future data. Here, we address these challenges in a cancer biology application and develop a meaningful similarity metric for ‘patches’—regions of simulated human cell membrane that express interactions between certain proteins of interest and relevant lipids. In the absence of well-defined conditions for similarity, we leverage several biology-informed notions about data and the underlying simulations to impose inductive biases on our metric learning framework, resulting in a suitable similarity metric that also generalizes well to significant distributional shifts encountered during the deployment. We combine these intuitions to organize the learned embedding space in a multiscale manner, which makes the metric robust to incomplete and even contradictory intuitions. Our approach delivers a metric that not only performs well on the conditions used for its development and other relevant criteria, but also learns key spatiotemporal relationships without ever being exposed to any such information during training.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

D–MOPH–25: diverse MOF–molecule pairs for Henry’s constants prediction

Computational methods like grand-canonical Monte Carlo simulations and machine learning (ML) have accelerated metal–organic frameworks (MOF) exploration but are typically limited to a narrow range of adsorbates due to data availability and force field constraints. In this study, we introduce a dataset of diverse MOF–molecule pairs for Henry’s constant prediction, D–MOPH–25, which systematically explores a diverse chemical space by combining 113 molecular adsorbates with over 5000 MOF structures through an active learning process. D–MOPH–25 constitutes the most diverse adsorbate dataset used in any ML study of molecular adsorption in MOFs to date. Our workflow builds a benchmark for predicting Henry’s constants at 300 K, leveraging conformal prediction for uncertainty quantification. Assessment through Shannon entropy and uniform manifold approximation and projection confirms the comprehensiveness of D–MOPH–25 while highlighting the importance of robust classification to filter out unphysical data points in regression tasks. Although future enhancements in model architecture and sampling criteria could improve predictive performance, our dataset already spans the target space using only 2.31% of total possibilities. This comprehensive dataset facilitates assessment of model generalizability across adsorbate species and can establish a foundation for high-throughput MOF screening and ML-driven separation processes.

active learning↗

DOE OSTI2026

PAC in DESI. II. Galaxy-halo connection into the $10^{6}{\rm M}_{\odot}$ frontier

Understanding dwarf galaxy formation is crucial for testing dark matter models and reionization physics. However, constructing stellar-mass complete spectroscopic samples at low masses is increasingly difficult, and the potential existence of a local void complicates studies in an average environment. The Photometric object Around Cosmic webs (PAC) method, which combines deep photometric and spectroscopic data to measure the excess surface density $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ of photometric objects around spectroscopic tracers, offers a promising path forward. We model 349 $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ measurements from DESI Y1 BGS and DECaLS, reaching $M_*=10^{6.4}\,{\rm M}_{\odot}$, using a stellar mass-halo mass relation (SHMR)-based subhalo abundance matching framework applied to two high-resolution $N$-body simulations from the Jiutian suite. The resulting SHMR is constrained down to $M_{\rm h}\simeq10^{8.0}\,h^{-1}{\rm M}_{\odot}$, revealing a clear upturn at $\sim10^{10.0}\,h^{-1}{\rm M}_{\odot}$ toward lower masses, indicating rising star-formation efficiency (SFE) in small haloes. This feature persists under extensions of the model that allow mass-dependent scatter, reionization-induced suppression of the halo occupation fraction, galaxy assembly bias, and alternative cosmologies. Combining with the results from Paper I, we find that central red galaxies dominate the low-mass regime. Our results motivate a hypothesis in which SFE is significantly higher than previously thought prior to reionization, enabling relatively massive galaxies to form in small haloes. These systems are subsequently quenched by the UV background, producing the central red dwarf galaxies observed. Finally, we obtain $3σ$ and $5σ$ upper mass bounds of $10^{8.80}\,h^{-1}{\rm M}_{\odot}$ and $10^{10.24}\,h^{-1}{\rm M}_{\odot}$ on the smallest haloes required to exist.

Xu, Kun [Pennsylvania U.; Durham U., ICC; Tsung-Da↗

DOE OSTI2026

Impact of QCD energy evolution on observables in heavy-ion collisions

Here, we study how the inclusion of energy dependence as dictated by quantum chromodynamic (QCD) small-x evolution equations affects key observables in ultra-relativistic heavy-ion collisions. Specifically, we incorporate JIMWLK evolution into the IP-Glasma framework, which serves as the initial condition for a simulation pipeline that includes viscous relativistic hydrodynamics and a hadronic afterburner. This approach enables a consistent modeling of highly energetic nuclei across varying Bjorken-x values, which are relevant for different collision energies and rapidity regions. In comparison to the standard IP-Glasma setup without small-x evolution, we observe pronounced changes in particle multiplicities and spectral distributions, especially in smaller systems and at the highest available energies. We further explore effects on anisotropic flow observables and correlations between mean transverse momentum and elliptic flow. Our findings underscore the critical role of nonlinear QCD evolution in accurately modeling the early stages of heavy-ion collisions, as well as its implications for extracting transport properties of the quark-gluon plasma.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Asynchronous Grid Connections Providing Fast-Frequency Response: System Integration Study

This paper presents an integration study for the recent power electronic-based fast-frequency response technology, "asynchronous grid connection" which operates as an aggregator for behind-the-meter resources and distributed generators. Both technical feasibility and techno-economic viability studies are presented. The fast-frequency response characteristics, validated against Power Hardware-in-the-Loop experiments, are integrated into an IEEE 9- bus system in DigSilent PowerFactory for system-level dynamic analysis. It demonstrates that droop-based control enhancements to local distributed generators allow their aggregation to provide grid-supporting functionalities and participate in the ancillary service markets. To this end, a long-term simulation embedding the system within the ancillary service market framework of PJM has been performed. The fast-frequency response regulation is subsequently used to calculate the potential revenue and project the results on a 15-year investment horizon. Finally, the techno-economic analysis provides recommendations for enhancements to access the full potential of distributed generators on a technical and regulatory level.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

MAGIC: M arching Cubes Isosurface Uncertainty Visualization for G auss i an Uncertain Data With Spatial C orrelation

Here, in this paper, we study the propagation of data uncertainty through the marching cubes algorithm for isosurface visualization for correlated uncertain data. Consideration of correlation has been shown paramount for avoiding errors in uncertainty quantification and visualization in multiple prior studies. Although the problem of isosurface uncertainty with spatial data correlation has been previously addressed, there are two major limitations to prior treatments. First, there are no analytical formulations for uncertainty quantification of isosurfaces when the data uncertainty is characterized by a Gaussian distribution with spatial correlation. Second, as a consequence of the lack of analytical formulations,existing techniques resort to a Monte Carlo sampling approach, which is expensive and difficult to integrate into visualization tools. To address these limitations, we present a closed-form framework to efficiently derive uncertainty in marching cubes level-sets for Gaussian uncertain data with spatial correlation (MAGIC). To derive closed-form solutions, we leverage the Hinkley's derivation on the ratio of Gaussian distributions. With our analytical framework, we achieve a significant speed-up and enhanced accuracy of uncertainty quantification over classical Monte Carlo methods. We further accelerate our analytical solutions using many-core processors to achieve speed-ups up to 585× and integrability with production visualization tools for broader impact. We demonstrate the effectiveness of our correlation-aware uncertainty framework through experiments on meteorology, urban flow, and astrophysics simulation datasets.

Gaussian↗

DOE OSTI2022

Determining spatially varying profit‐maximizing management practices for miscanthus and switchgrass production in the rainfed United States

Abstract Determining optimal management practices for the profitable production of perennial energy crops is critical for scaling up production beyond experimental levels. Although many experimental field studies have examined the effects of management practices on the performance of miscanthus and switchgrass, there are no recommendations for economically optimal nitrogen (N) application rates and how they should vary spatially and with the age of the energy crop as well as on optimal rotation age of the energy crop to maximize profits. We develop a modeling framework to determine economically optimal crop management decisions and simulate the variability under various scenarios for miscanthus and switchgrass production across 2287 counties in the rainfed United States. We find that profit‐maximizing N recommendations for these crops vary across maturity stages and regions and can increase the landowner's profits compared with a uniform N rate across ages and regions. We also find that the optimal rotation for these crops is shorter than the productive physical lifespan (15–20 and 10 years for miscanthus and switchgrass, respectively). Specifically, the N rate that maximizes the economic returns is negligible for miscanthus and 111 kg ha −1 for switchgrass production at age 2. The mean profit‐maximizing N rate increases with age for miscanthus, peaking at 151 kg ha −1 at age 11 before declining to 114 kg ha −1 at the optimal rotation age of 13 years while that for switchgrass is 150 kg ha −1 for middle‐aged stands and declines to 114 kg ha −1 at the optimal rotation of 8–9 years. We find that miscanthus is the most profitable energy crop in the northern region of the rainfed United States while switchgrass is most profitable in the south of the rainfed United States. Our findings are useful for improving assessments of the profitability of energy crops and guiding future management decisions by landowners.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Targeted Adaptive Design

Modern advanced manufacturing and advanced materials design often require searches of relatively high-dimensional process control parameter spaces for settings that result in optimal structure, property, and performance parameters. The mapping from the former to the latter must be determined from noisy experiments or from expensive simulations. Here, we abstract this problem to a mathematical framework in which an unknown function from a control space to a design space must be ascertained by means of expensive noisy measurements, which locate control settings generating desired design features within specified tolerances, with quantified uncertainty. We describe targeted adaptive design (TAD), a new algorithm that performs this sampling task efficiently. TAD creates a Gaussian process surrogate model of the unknown mapping at each iterative stage, proposing a new batch of control settings to sample experimentally and optimizing the updated expected log-predictive probability density of the target design. TAD either stops upon locating a solution with uncertainties that fit inside the tolerance box or uses a measure of expected future information to determine that the search space has been exhausted with no solution. TAD thus embodies the exploration-exploitation tension in a manner that recalls, but is essentially different from, Bayesian optimization and optimal experimental design.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Baycal

Bayesian Model Calibration (BayCal) toolkit is a software plugin for Risk Analysis Virtual Environment (RAVEN) framework, arming at inversely quantifying the uncertainties associated with simulation model parameters based on available experiment data. BayCal seeks statistical inference of the uncertain input parameters that are consistent with the available measurement data or observed data. The unique feature of BayCal is the capability to be linked with RAVEN to build corresponding calibration workflows for complex multi-physics simulations. In addition to be able to use the machine learning capability of RAVEN to significantly reduce the computational cost of expensive simulation models, another distinctive feature of BayCal is the capability to deal with high-dimensional correlated model outputs, such as time series observations at multiple locations, via principal component analysis (PCA) technique.

Wang, Congjian↗

DOE OSTI2021

Development of Distributed Microgrid Controllers Using HELICS [SWR-21-94]

Development of Distributed Microgrid Controllers Using HELICS uses HELICS framework develop distributed microgrid controllers (DMCs) for coordinated grid service. Development of Distributed Microgrid Controllers Using HELICS includes thirteen agents, including OpenDSS agent, data manager agent, and eleven DMC agents. Development of Distributed Microgrid Controllers Using HELICS can be used as a standard framework to develop fully distributed control solutions that dispatches simulated DERs in simulation platform.

Wang, Jing↗

DOE OSTI2024

SAGESim

SAGESim: Scalable Agent-based GPU-Enabled Simulator is a distributed GPU agent-based modeling framework in Python. SAGESim enables researchers to study complex adaptive systems at large scales on high-performance computing systems

Gunaratne, Chathika↗

DOE OSTI2025

Simulation-Based Inference for Neutrino Interaction Model Tuning

This project demonstrates, for the first time, the application of simulation-based inference (SBI) techniques to tune neutrino–nucleus interaction models. Using a mock dataset based on the MicroBooNE tuning of the GENIE event generator, our approach employs a Neural Posterior Estimator (NPE) with Masked Autoregressive Flows (MAF) to infer the posterior distributions of key GENIE parameters directly from simulated histograms. The workflow provides a scalable and amortized framework for performing likelihood-free inference in high-dimensional parameter spaces, offering a pathway to more efficient and uncertainty-aware model tuning for next-generation neutrino experiments such as DUNE and SBND.

Tame-Narvaez, KarlaMaria [Fermi National Accelerat↗