SEARCH · Search NASA

Results for “Large eddy 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.

522 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Atomistic Simulation of Glasses and Amorphous Materials: Challenges and Opportunities for the Next Decade

Atomistic simulations have become indispensable tools for understanding glass structure, dynamics, and properties, yet persistent challenges limit their predictive power. This perspective examines three interconnected issues, namely glass formation procedures, interatomic potential development, and machine learning applications, which emerged from the 5th International Workshop on Challenges of Atomistic Simulations of Glasses and Amorphous Materials. We identify convergent community priorities for (i) standardized validation protocols, (ii) curated benchmark datasets with complete metadata, and (iii) open repositories for glasses. A systematic was forward is provided by a hierarchical validation framework for assessing the structural fidelity, property prediction, and behavioral realism of simulation techniques. Looking ahead, transformative advances are promised by the fusion of classical techniques with machine learning based approaches, for instance, by integrating swap Monte Carlo with machine-learning (ML) potentials, leveraging foundation models through transfer learning, and finetuning ML potentials with experimental data. Progress depends on the community committing to validated models, reproducible protocols, and sustained data sharing.

Krishnan, N. M. Anoop↗

DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗

DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Wettability of Two-Dimensional Carbon Allotropes from Molecular Simulations

Force-field Monte Carlo and Molecular Dynamics simulations are used to compare wetting behaviors of model carbon sheets mimicking neat graphene, its saturated derivative, graphane, and related planar allotropes penta-graphene, γ-graphyne, and ψ-graphene in contact with aqueous droplets or an aqueous film confined between parallel carbon sheets. Atomistic and area-integrated surface/water potentials are found to be essentially equivalent in capturing moderate differences between the wetting free energies of tested substrates. Despite notable differences in mechanical and electric properties of distinct allotropes, the predicted allotrope/water contact angles span a narrow window of weakly hydrophilic values. Contact angles in the range of 80 ± 10° indicate modest hydration repulsion incapable of competing with van der Waals attraction between carbon particles. Poor dispersibility in neat water is hence a common feature of studied materials.

Biochemistry & Molecular Biology↗

DOE OSTI2025

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗

DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗

DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Formation and Shape Changing of Conductive Helical Ribbons via Deposition of Highly Stressed Films on Mechanically Responsive Substrates

Abstract This work demonstrates that the electrodeposition of highly stressed films on compliant ribbons is a robust process to obtain helical structures with excellent mechanical stability and potentially high thermal and electrical conductance. Electrodeposition on end‐tethered ribbons alters their axial and bending stiffness while imparting mechanical stress to drive the formation of a helix with a microscale diameter and pitch in a controlled and scalable manner. The process generates helices with diameters and pitches between 80 and 200 µm and lengths as large as several millimeters. The approach is amenable to parallel processing a large number of 3D structures on any substrate, including large‐area semiconductor wafers. This phenomenon is explained in terms of the change of stress gradients as material is added. Applications of the fabricated helices include antennas, metamaterials, and slow‐wave structures in frequency ranges not previously attainable.

Chemistry↗

DOE OSTI2025

Photon–photon chemical thermodynamics of frequency conversion processes in highly multimode systems

Abstract Frequency generation in highly multimode nonlinear optical systems is inherently a complex process, giving rise to an exceedingly convoluted landscape of evolution dynamics. While predicting and controlling the global conversion efficiencies in such nonlinear environments has long been considered impossible, here, we formally address this challenge even in scenarios involving a very large number of spatial modes. By utilizing fundamental notions from optical statistical mechanics, we develop a universal theoretical framework that effectively treats all frequency components as chemical reactants/products, capable of undergoing optical thermodynamic reactions facilitated by a variety of multi-wave mixing effects. These photon–photon reactions are governed by conservation laws that directly determine the optical temperatures and chemical potentials of the ensued chemical equilibria for each frequency species. In this context, we develop a comprehensive stoichiometric model and formally derive an expression that relates the chemical potentials to the optical stoichiometric coefficients, in a manner akin to atomic/molecular chemical reactions. This advancement unlocks new predictive capabilities that can facilitate the optimization of frequency generation in highly multimode photonic arrangements, surpassing the limitations of conventional schemes that rely exclusively on nonlinear optical dynamics. Notably, we identify a universal regime of Rayleigh–Jeans thermalization where an optical reaction at near-zero optical temperatures can promote the complete and entropically irreversible conversion of light to the fundamental mode at a target frequency. Our theoretical results are corroborated by numerical simulations in settings where second-harmonic generation, sum-frequency generation and four-wave mixing processes can manifest.

Optics↗

DOE OSTI2026

AutoBEM: A scalable framework for nationwide building energy simulation and retrofit evaluation in the United States

This paper presents AutoBEM, an integrated, automated framework for nationwide building energy modeling and retrofit evaluation in the United States. Unlike prior UBEM platforms that either rely primarily on representative stock sampling or operate at city scale, AutoBEM automates the generation of building-resolved, physics-based EnergyPlus/OpenStudio simulation models at national scale using GIS-derived geometry, prototype-based assumptions, and standardized scalable workflows. Leveraging the Model America dataset and high-performance computing, AutoBEM generates and simulates energy models for 122.9 million buildings, representing 97.8% of the U.S. building stock. These models are being made publicly and freely available as the Model America v1.0 (MAv1) dataset. AutoBEM supports detailed, building-level assessments of energy consumption, CO2 emissions, and post-processed anthropogenic heat emissions (AHE), and evaluates 151 energy conservation measures (ECMs) using localized utility pricing and building characteristics. In addition, AutoBEM incorporates both typical and future climate conditions through integration with Typical Meteorological Year (TMY) and Future TMY (fTMY) weather data derived from IPCC scenarios. In a case study of Phoenix, Arizona, AutoBEM identified several high-efficiency HVAC upgrades and selected envelope measures with short modeled payback periods (1.5 years) for certain building types and standards. Simulations under future climate scenarios (SSP5–RCP8.5) project an 11.3% increase in electricity use and a 32% reduction in natural gas demand by 2100, underscoring the need for climate-adaptive retrofit planning. By enabling reproducible, bottom-up, and location-specific analysis at scale, AutoBEM provides a step toward a national digital twin of the built environment and supports data-driven screening and planning for decarbonization, resilience, and energy equity.

Li, Hang [ORNL] (ORCID:0000000306001920)↗

DOE OSTI2026

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Multi-Resolution UAV Path Replanning for Inspection of Tailings Dams

Autonomous inspection of large and complex structures with a commercial unmanned aerial vehicle (UAV) is a challenging problem that has been addressed in recent years. In this paper, we address the global motion planning problem of creating autonomous inspection missions for UAVs considering photogrammetry constraints. We focus on the inspection of large tailings dams, which are dam structures used to store waste byproducts of mining. Our method uses a prior sparse point cloud of the dam to generate a voxel grid, where paths satisfying photogrammetry constraints are tested for collisions. We then apply the A* algorithm as a local planner to avoid obstacles within the global mission. Moreover, we address the problem of changing routes online by using octree-based multi-resolution grids for efficient and fast pathfinding. Our results, obtained using tridimensional maps of an actual coal mine tailings dam, show that using octrees for multi-resolution motion planning is faster than using a fixed voxel grid in online missions while inspecting large structures.

42 ENGINEERING↗

DOE OSTI2026

The MuFusE large-volume diamond anvil cell for exploring muon-catalyzed fusion at higher pressures and temperatures

A new large-volume diamond anvil cell (DAC) has been developed for the Muon-catalyzed Fusion (μCF) Experiment (MuFusE), enabling the compression and heating of deuterium–tritium (d–t) mixtures to pressures and temperatures needed to advance μCF research. The MuFusE DAC achieves the large sample volumes necessary for high-precision fusion measurements while integrating cryogenic loading, all-metal sealing, and flexible bellows to maintain a secure environment during cell compression. Combined with remote pneumatic actuation and secondary containment, the DAC safely managed a 25 Ci tritium inventory while providing a clear optical path for in situ measurements of sample pressure and composition via laser spectroscopy. Utilizing 5 mm diameter diamond anvils oriented in the path of a high-intensity muon beam, the apparatus achieved a stable sample volume of 19.2 mm 3 at liquid density, pressures up to 933 MPa and temperatures up to 400 K—benchmarks that significantly exceed previously reported limits for static d–t targets.

Kalow, J. D. [Acceleron Fusion, Inc., Cambridge, M↗

DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗

DOE OSTI2026

Dark Energy Survey year 6 results: Magnification modeling and its impact on galaxy clustering and galaxy-galaxy lensing cosmology

Gravitational lensing magnification alters the observed spatial distribution of galaxies and must be accounted for to prevent biases in cosmological probes of the large-scale structure. We investigate its effects on the Dark Energy Survey Year 6 galaxy clustering and galaxy-galaxy lensing analyses using the fiducial lens (position tracer) sample M ag L im++. Magnification bias is parameterized by a coefficient that describes the response of the number of selected objects per unlensed area element to a change in the lensing convergence. We quantify this coefficient using the BALROG synthetic source injection catalog to account for the complexity of the selection function, and compare these results with simplified estimates. The resulting values of the magnification coefficients for each redshift bin are [3.16 ± 0.08, 2.76 ± 0.21, 4.09 ± 0.15, 4.42 ± 0.16, 4.90 ± 0.29, 4.83 ± 0.25]. Relative to Year 3, this analysis provides more precise and accurate magnification bias estimates through a larger BALROG area and reweighting to better match the data properties. Here, the cosmological results are robust when tested against various magnification parameter prior choices and also when adding cross-clustering between lens redshift bins. Neglecting magnification, however, introduces significant systematic shifts: relative to the fiducial analysis with Gaussian priors centered on the BALROG -derived estimates, we observe shifts of 1.37σ in S 8 and -0.84σ in Ω m (with cosmic shear included: -0.61σ in S 8 and -0.71σ in Ω m ), in agreement with findings from simulated data, demonstrating that magnification must be modeled to avoid biases. Freeing the magnification bias in lens bin 2 leads to unphysical negative values, further justifying its exclusion from the fiducial Year 6 analysis.

Cosmological parameters↗

DOE OSTI2025

Unveiling the lithium-ion transport mechanism in Li{sub 2}ZrCl{sub 6} Solid-State Electrolyte {ital via} deep learning-accelerated molecular dynamics simulations.

Lithium zirconium chlorides (LZCs) present a promising class of cost-effective solid electrolytes for next-generation all-solid-state batteries. The unique crystal structure of LZCs plays a crucial role in facilitating lithium-ion mobility, which further affects the electrochemical performance. To understand the underlying mechanism governing ion transport, we employed deep learning-accelerated molecular dynamics simulation on Li2ZrCl6 (trigonal alpha- and monoclinic beta-LZC), focusing specifically on the zirconium coordination environment. Our results reveal that disordered alpha-LZC exhibits the highest ionic conductivity, while beta-LZC demonstrates significantly lower conductivity, closely aligning with experimental findings. The study confirms that across all phases, lithium migration proceeds via the site-to-site hopping mechanism, where variations in site residence times critically impact the overall ionic conductivity. In alpha-LZCs, lithium ions prefer to anisotropically diffuse across interlayers as the result of a lower energy barrier, driven primarily by collective diffusion. In contrast, lithium ions in beta-LZC primarily isotropically diffuse within the intralayer, hindered by higher energy barriers and determined by individual diffusion. The variation in ZrCl6 2- octahedral unit softening, induced by the specific layered arrangement of zirconium atoms, emerges as a critical determinant of the energy barriers across the LZC phases. These atomic-scale insights into the transport processes provide valuable guidance for the rational design and optimization of LZCs-based electrolytes, accelerating their practical application in advanced energy storage technologies.

Guo, Hanzeng↗

DOE OSTI2026

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗