Engineering Papers⌕ Search

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.

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

At least 1,729 records · Page 96

GDSA Repository Systems Analysis Investigations in FY 2024

The Disposal Research and Development (R&D) Program of the US Department of Energy (DOE) office of Nuclear Energy (NE-8) Spent Fuel and Waste Science and Technology (SFWST) Campaign is to conduct R&D on disposal of spent nuclear fuel (SNF) and high-level waste (HLW). The goal of the Geologic Disposal Safety Assessment (GDSA) within this project is to develop a disposal system modeling and analysis capability that supports the integrated modeling of coupled processes controlling disposal system performance of deep geologic repositories, including uncertainty. This report describes specific activities in the Fiscal Year (FY) 2024 associated with the GDSA Repository Systems Analysis (RSA) work package. The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and repository system performance models in crystalline, argillite, salt, and unsaturated alluvium potential host-rock environments, and to simulate and analyze these generic repository concepts and models using GDSA Framework toolkit, and other tools as needed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An FSI Reweighter for 𝑣 Interactions on Heavy Nuclei

Final State Interaction (FSI) processes in $\nu$--Pb differ significantly from those in neutrino--carbon ($\nu$--C), and MINER$\nu$A is particularly sensitive to these differences, especially in pion absorption. This work develops an FSI reweighting method for $\nu$--Pb (and C) interactions in MINER$\nu$A. We compare GENIE~2, the currently used generator version, with GENIE~3 predictions of the initial hadron energy as a function of FSI fate, and apply these weights to simulated observables. We identify changes in the FSI composition, providing a framework for estimating FSI-related systematic uncertainties and interpreting their effects in the data.

Moreno-Palacios, Oscar [Coll. William and Mary]↗

White Paper: Scalable Digital Twin Capabilities for Aging and Surveillance of Engineered Systems

This white paper presents a multi-year initiative to develop practical, secure, and scalable digital twin capabilities for engineered systems in aging and surveillance contexts—an approach pioneered at the National Nuclear Security Administration (NNSA) Lawrence Livermore National Laboratory (LLNL) that maps directly onto the needs and ambitions of the Navy for ship- and fleet-level digital twins. LLNL’s work in building part- and process-level digital twins for advanced manufacturing, with a vision to scale up to entire factory floors and, ultimately, enterprise-wide digital twins, offers an adaptable pathway for the Navy as it seeks to modernize lifecycle management, readiness, and predictive maintenance across ships and fleets. For our application, we integrate physics-based modeling with automated data ingestion, processing, and AI-driven calibration, creating hybrid models that are both interpretable and data responsive. We modernized legacy workflows, established centralized data infrastructure, automated experimental pipelines, and demonstrated end-to-end coupling of accelerated aging data with finite element simulations via optimization and surrogate modeling. The result is a generalizable framework that supports part-level digital twins today and lays the groundwork for future system-level twins suitable for Navy applications.

36 MATERIALS SCIENCE↗

Core-corona hadronization model and muon puzzle in extensive air showers

The discrepancy between models and measurements concerning the muon content in air showers produced by ultra-high energy cosmic rays, the so-called muon puzzle, remains unsolved. String fragmentation models currently used in hadronic interactions fail to reproduce the observations, while recent measurements at the LHC hint towards the existence of production mechanisms, such as collective statistical hadronization, which lead to an increase in the muon production after a hadronic interaction. The core-corona model of heavy ions implements both production mechanisms, where the large-density region of interaction (core) hadronizes statistically, while the low-density (corona) hadronizes following string fragmentation. In this contribution, we present a modification of the CONEX framework to implement the core-corona model in air shower simulations. We demonstrate a significant impact of the core effect, as observed at the LHC, on the muon content in air showers produced by ultra-high energy cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗

Development and Validation of a Latent Thermal Energy Storage Model Using Modelica

An abundance of research has been performed to understand the physics of latent thermal energy storage with phase change material. Some analytical and numerical findings have been validated by experiments, but there are few free and open-source models available to the general public for use in systems simulation and analysis. The Modelica programming language is a good avenue to make such models available, because it is object-oriented, equation-based, declarative, and acausal. These characteristics have enabling the creation of component model libraries that can be used to build larger system simulations for design analysis. The authors have previously developed a numerical framework to model phase change thermal storage and have validated model predictions with experiments. The objectives of this paper are to describe the transfer of the numerical framework to an implementation in a Modelica component model and to validate the Modelica model with data from the experiment and the original numerical framework.

25 ENERGY STORAGE↗

flat10MIP: an emissions-driven experiment to diagnose the climate response to positive, zero and negative CO2 emissions

Abstract. The proportionality between global mean temperature and cumulative emissions of CO2 predicted in Earth system models (ESMs) is the foundation of carbon budgeting frameworks. Deviations from this behavior could impact estimates of required net-zero timings and negative emissions requirements to meet the Paris Agreement climate targets. However, existing ESM diagnostic experiments do not allow for direct estimation of these deviations as a function of defined emissions pathways. Here, we perform a set of climate model diagnostic experiments for the assessment of transient climate response to cumulative CO2 emissions (TCRE), the Zero Emissions Commitment (ZEC), and climate reversibility metrics in an emissions-driven framework. The emissions-driven experiments provide consistent independent variables simplifying simulation, analysis and interpretation, with emissions rates more comparable to recent levels than existing protocols using model-specific compatible emissions from the CMIP DECK 1pctCO2 experiment, where emissions rates tend to increase during the experiment, such that at the time of CO2 doubling in year 70, emissions are much greater than present-day values. A base experiment, “esm-flat10”, has constant emissions of CO2 of 10 GtC per year (near-present-day values), and initial results show that the TCRE estimated in this experiment is about 0.1 K less than that obtained using 1pctCO2. A subset of ESMs exhibit land carbon sinks that saturate during this experiment. A branch experiment, esm-flat10-zec, illustrates that both positive and negative ZEC effects are less pronounced under esm-flat10 than under 1pctCO2 – the magnitude of ZEC50 in ESMs is, on average, reduced by 30 % compared with 1pctCO2 branch experiments. A final experiment, esm-flat10-cdr, assesses climate reversibility under negative emissions, where we find that peak warming may occur before or after net zero and that the asymmetry in temperature at a given level of cumulative emissions between the positive and negative emissions phases is well described by ZEC in most models. Further, we find that existing probabilistic simple climate model (SCM) ensembles tend to overestimate temperature reversibility compared with ESMs, highlighting the need for additional constraints. We propose a set of climate diagnostic indicators to quantify various aspects of climate reversibility. These experiments were suggested as potential candidates in CMIP7 and have since been adopted as “fast track” simulations.

Sanderson, Benjamin M↗

Exploring Building Retrofit Strategies Using AutoBEM Under Future Weather Scenarios

This study evaluates the long-term effectiveness of energy conservation measures (ECMs) on building energy consumption using AutoBEM, a scalable modeling framework driven by the high-resolution Model America dataset. We simulated 18,951 buildings in Flagstaff, Arizona under four climate scenarios using Future Typical Meteorological Year (fTMY) weather files for six time periods spanning from 1980 to 2099. Six ECMs were analyzed across electricity and gas usage, including HVAC fuel-switching, insulation upgrades, and infiltration control. While some measures, such as reducing space infiltration by percentage, showed minimal or even negative impact on total energy savings at the aggregate level, they proved highly effective for specific building types. Conversely, HVAC electrification offers high gas reduction but shifts demand to electricity, highlighting critical trade-offs under different climate trajectories. Building-type-specific analysis under SSP5-RCP8.5 (2080–2099) revealed significant variation in ECM performance, underscoring the need for targeted retrofit strategies. This study demonstrates the power of combining fTMY projections with large-scale simulations to inform data-driven retrofit planning.

Chowdhury, Shovan [ORNL]↗

Applications of Multibody for Everybody (M4E) in Marine Energy

Multibody for Everybody (M4E) is an open-source symbolic dynamics modeling framework designed to automate the derivation of equations of motion and simulation of constrained multibody systems using the joint coordinate formulation.

16 TIDAL AND WAVE POWER↗

Photoproduction of the b 1 (1235) Meson off the proton at E gamma = 6-12 GeV

The GlueX Experiment at Jefferson Lab directs a linearly polarized photon beam on a liquid hydrogen target surrounded by an almost hermetic detector. The experiment aims to study the meson spectrum in the light-quark sector and search for exotic spin-parity states predicted by lattice QCD calculations. The lightest exotic candidate is the π 1 meson, which is predicted to decay dominantly to b 1 π. In this thesis, we present efforts to characterize the ωπ 0 decay channel of the b 1 meson as precursor to an analysis of the sought-after π 1 exotic state. We present an introduction to the physics involved and an overview of the detector subsystems with a focus on the Barrel Calorimeter gain monitoring system, as part of the service work expected by the GlueX Collaboration. We discuss the data analysis method used and present two frameworks, for angular moments and partial waves analyses. Details of the simulation used to calculate the detector acceptance. A proof of concept partial wave analysis is presented along with experimental angular moments as first step to a full analysis of the angular distribution Our partial wave analysis indicates that both S- and D-waves are needed, in qualitative agreement with theoretical expectations. The cross-section of the omega π 0 channel is extracted to be of the scale of 1 μb in agreement with previous measurements. The differential cross-section indicates the presence of two production processes with p 1 = -5:24 ± 0:04 for 0:25 < -t < 0:95 GeV 2 ~c 2 and p 2 = -1:24 ± 0:05 for 0:95 < -t < 2:0 GeV 2 ~c 2 , which does not agree with results from previous experiments. The s-channel helicity conservation and helicity amplitudes of the b 1 meson are presented as additional experimental observables. Our results confirm that the b1 photoproduction process does not conserve s-channel helicity though they do not align with previous measurements. The helicity amplitude of the omega are extracted to be |F 1 | 2 = 0:3037 ± 0:0003 which does not agree with the expected value of a b 1 decay or results from a previous experiment. Future steps in all these analyses will be continued by a new graduate student in the group, culminating in a publication.

Foda, Ahmed↗

Optimizing Transportation Networks for E-Waste Reverse Logistics: A Multi-Modal Cost Allocation and Pricing Strategy

The exponential growth of electronic waste (e-waste) poses critical challenges for sustainable reverse logistics and transportation network optimization. This study develops a dual-channel transportation framework for e-waste logistics that integrates dynamic freight pricing, cost allocation mechanisms, and game-theoretic coordination. The model captures interactions between centralized hubs and distributed processing networks, accounting for freight rate elasticity, volume allocation, and capacity constraints. Using Stackelberg game theory and cost-sharing strategies, the framework optimizes transportation efficiency and profit distribution across logistics channels. Numerical simulations show that the dual-channel structure increases centralized hub profit by 226.8% compared to baseline single-channel operations, while boosting total transported volume by 1.2% and nearly doubling freight collector profit under cost-sharing. Scenario analyses across regional infrastructures reveal that network density, policy incentives, and logistics costs shape routing efficiency and profit allocation. These findings suggest that coordinated strategies combining dynamic pricing, targeted infrastructure investment, and strategic cost allocation are needed to design efficient, resilient, and regionally adaptable e-waste transportation systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Small-scale fluctuating magnetic fields of order nG are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel 3D magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use 384 × 512 × 384 tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay (E kin ∝ t -1.15 ) without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short (~0.3 of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian–Lagrangian bridge.

79 ASTRONOMY AND ASTROPHYSICS↗

SimPoint++: Less Simulation Points

We propose SimPoint++, which replaces the BIC method by combining WCSS and Silhouette to find the optimal cluster. The new Python framework of SimPoint++ provides a dimension reduction pipeline for effective clustering and supports multi-thread application analysis. We evaluate SimPoint++ with Spec CPU 2017 benchmarks. SimPoint++ achieves comparable or higher accuracy with significantly fewer simulation points, resulting in a 5x speed-up in simulation time compared to state-of-the-art solutions.

Chen, Si↗

Modeling pre-Exascale AMR Parallel I/O Workloads via Proxy Applications

The present work investigates the modeling of preexascale input/output (I/O) workloads of Adaptive Mesh Refinement (AMR) simulations through a simple proxy application. We collect data from the AMReX Castro framework running on the Summit supercomputer for a wide range of scales and mesh partitions for the hydrodynamic Sedov case as a baseline to provide sufficient coverage to the formulated proxy model. The non-linear analysis data production rates are quantified as a function of a set of input parameters such as output frequency, grid size, number of levels, and the Courant-Friedrichs-Lewy (CFL) condition number for each rank, mesh level and simulation time step. Linear regression is then applied to formulate a simple analytical model which allows to translate AMReX inputs into MACSio proxy I/O application parameters, resulting in a simple “kernel” approximation for data production at each time step. Results show that MACSio can simulate actual AMReX nonlinear “static” I/O workloads to a certain degree of confidence on the Summit supercomputer using the present methodology. The goal is to provide an initial level of understanding of AMR I/O workloads via lightweight proxy applications models to facilitate autotune data management strategies in anticipation of exascale systems.

Godoy, William↗

Recovery of MOF-5 from Extreme High-Pressure Conditions Facilitated by a Modern Pressure Transmitting Medium

Here mechanisms underlying the mechanically induced amorphization of metal-organic frameworks (MOFs) are of current interest, and both high-pressure experimentation and molecular dynamics simulations have been used to reveal the fundamentals of load bearing, deformation, and pressure-induced amorphization (PIA) in these highly porous materials. Unfortunately, MOFs are typically highly susceptibility to amorphization, which limits the conditions under which they can be processed and used. However, their flexible structures can be stabilized to high pressures by incorporating guest species into the framework matrix. In this study, a large molecule pressure transmitting medium (DAPHNE 7575 ) is used as a structure-fortifying guest species to stabilize the prototypical MOF -5 to ultra-high pressures ( >9 GPa) and enable the recovery of crystalline material upon decompression. Structural changes associated with the penetration of the pressure transmitting medium on compression are examined using a combination of high pressure synchrotron powder diffraction and molecular dynamics simulations. This work enhances the understanding of PIA in MOFs while showcasing a potential route for the stabilization of MOFs to surprisingly high-pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating Resiliency of Transportation Network under Targeted and Potential Climate Change Disruptions

Ensuring robustness and resilience in intermodal transportation systems is essential for the continuity and reliability of global logistics. These systems are vulnerable to various disruptions, including natural disasters and technical failures. Despite significant research on freight transportation resilience, investigating the robustness of the system after targeted and climate-change-driven disruption remains a crucial challenge. Drawing on network science methodologies, this study models the interdependencies within the rail and water transport networks and simulates different disruption scenarios to evaluate system responses. Here, we use the data from the U.S. Department of Energy Volpe Center for network topology and tonnage projections. The proposed framework quantifies deliberate, stochastic, and climate-driven infrastructure failure, using higher resolution downscaled multiple Earth System Models’ simulations from Coupled Model Intercomparison Project Phase version 6. We show that the disruptions of a few nodes could have a larger impact on the total tonnage of freight transport than on network topology. For example, the removal of targeted 20 nodes can bring the total tonnage carrying capacity to 30% with about 75% of the rail freight network intact. This research advances the theoretical understanding of transportation resilience and provides practical applications for infrastructure managers and policymakers. By implementing these strategies, stakeholders and policymakers can better prepare for and respond to unexpected disruptions, ensuring sustained operational efficiency in transportation networks.

Climate Change Disruptions↗

Exploring scenarios for enhanced fuel compression and performance on the National Ignition Facility with machine-learning-aided design techniques

Recent fusion experiments on the National Ignition Facility (NIF) have achieved ignition, producing multi-MJ fusion yields for input laser energies of roughly 2 MJ [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)]. Building on the success of the target designs that have achieved ignition, we explore new implosion scenarios predicted to generate significantly more compression of the dense DT ice layer and correspondingly higher yields while preserving many of the key physics characteristics of present-day ignition designs. Our main result is a novel 3-shock implosion scheme that effectively minimizes the shock-induced entropy in the dense, accelerating DT shell and maximizes the resulting fuel compression subject to a fixed leading shock strength consistent with present-day ignition experiments, which is necessary to melt the crystalline high-density carbon ablator. Compared to the first NIF experiment to fulfill Lawson's ignition criterion, shot N210808 [Abu-Shawareb et al., Phys. Rev. Lett. 129, 075001 (2022)], our design exhibits a 40% increase in simulated peak areal density (ρR) and a 5× increase in 1D fusion yield using a 4% lighter ablator and identical DT payloads. We also present a complete integrated 2D hohlraum design and laser pulse specifications capable of generating the desired 3-shock drive and maintaining control of the low-mode capsule implosion symmetry, where the increase in simulated 2D yield relative to N210808 is > 10×. This new implosion regime was discovered with help from a machine-learning-enabled capsule design optimization framework. We outline the workflow this automated tool uses to identify improved design candidates by running several rounds of capsule simulations, constructing a surrogate model mapping input variations to key physics output quantities, and querying the resulting statistical model to propose adjustments to the x-ray drive and capsule to reach a set of physics objectives prescribed by the designer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Emulation of quantum correlations by classical dynamics in a spin-$\frac{1}{2}$ Heisenberg chain

We simulate the dynamical spin structure factor (DSSF) 𝒮⁡(𝑞,𝜔) of the spin-1/2 Heisenberg antiferromagnetic chain using classical simulations. By employing Landau-Lifshitz Dynamics, we emulate quantum correlations through temperature-dependent corrections, including rescaling of magnetic dipoles and renormalization of exchange interactions. Here, our results closely match Quantum Monte-Carlo calculations for 𝑘 B⁢ 𝑇/𝐽≳1, extending the applicability of classical dynamics to the challenging case of gapless excitations. At higher temperatures, our simulations comply with general predictions for uncorrelated paramagnetic fluctuations in the infinite temperature limit. Entanglement witnesses derived from the quantum-equivalent DSSF act as sensitive diagnostics for the quantum-to-classical crossover. Their reliability stems from their dependence on spectral features alone, enabling classical dynamics to emulate quantum thresholds without genuine entanglement. This framework also reproduces transverse spin correlations in finite magnetic fields, in agreement with quantum simulations. Together, our results establish quantum-corrected classical dynamics as a scalable and predictive tool for interpreting scattering experiments and exploring quantum correlations in strongly correlated spin systems.

Inelastic neutron scattering↗

Sources of stochasticity in the growth of cloud droplets: supersaturation fluctuations versus turbulent transport

The role played by fluctuations of supersaturation in the growth of cloud droplets is examined in this study. The stochastic condensation framework and the three regimes of activation of cloud droplets – namely, mean-dominant, fluctuation-influenced, and fluctuation-dominant, are used for analyzing the data from high-resolution large-eddy simulations of the Pi convection-cloud chamber. Based on a detailed budget analysis the significance of all the terms in the evolution of the droplet size distribution equation is evaluated in all three regimes. The analysis indicates that the mean-growth rate is a dominant process in shaping the droplet size distribution in all three regimes. Turbulence introduces two sources of stochasticity, turbulent transport and particle lifetime, and supersaturation fluctuations. The transport of cloud droplets plays an important role in all three regimes, whereas the direct effect of supersaturation fluctuations is primarily related to the activation and growth of the small droplets in the fluctuation-influenced and fluctuation-dominant regimes. Here, we compare our results against the previous studies (experimental and theory) of the Pi chamber, and discuss the limitations of the existing models based on the stochastic condensation framework. Furthermore, we extend the discussion of our results to atmospheric clouds, and in particular focus on recent adiabatic turbulent cloud parcel simulations based on the stochastic condensation framework, and emphasize the importance of entrainment/mixing and turbulent transport in shaping the droplet size distribution.

54 ENVIRONMENTAL SCIENCES↗