SEARCH · Search NASA

Results for “Structural modeling and simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2026

MARVEL Reactor Fuel Performance Report

The Microreactor Applications Research Validation and EvaLuation (MARVEL) project is producing a high temperature liquid metal-cooled nuclear test bed at Idaho National Laboratory (INL) to ultimately improve the integration of microreactors to end-user applications. This ambitious effort seeks to design, authorize, construct, test, and operate the reactor within five years. In order to construct and operate the MARVEL reactor in a timely manner, the system will utilize materials and component designs which have already been used, qualified, or licensed from previous reactors. The MARVEL reactor will be located at the INL Transient Reactor Test (TREAT) facility in the north high-bay equipment pit and will use the existing 304 stainless steel-clad U-ZrH1.6 pin-type fuel system developed by General Atomics and purchased from TRIGA International. This fuel has been previously qualified under the United States Department of Energy’s (US DOE) Reduced Enrichment for Research and Test Reactors (RERTR) Program. Even though the regulator of the MARVEL reactor is the US DOE, the standards and overall approach recommended by the Nuclear Regulatory Commission is well-defined and utilized here. Following NUREG-1537 regulatory guidance, this report documents the authorization case for the MARVEL fuel system’s application to MARVEL and establishes stable and predictable fuel performance during the most thermophysically unfavorable conditions achievable in the MARVEL reactor. To that end, this report provides a comprehensive survey of the known thermophysical properties, performance, and quantitative relationships associated with the MARVEL reactor fuel element and uses this information to determine its mechanical integrity and risk of reaching unacceptable conditions during the most extreme accident scenarios predicted for the reactor using the most conservative assumptions available. The information contained herein is compiled from a combination of historical reports and peer reviewed scientific publication manuscripts. Known mechanisms under which the fuel is susceptible to failure are highlighted and compared to conditions that could exist in the MARVEL reactor during an unanticipated transient or accident scenario. The two scenarios considered for analysis in this report are (1) an unprotected loss of flow accident and (2) a hypothetical unprotected loss of coolant accident during the loss of flow accident. Preliminary 2D steady-state analyses herein indicate that both fuel-cladding chemical interactions and fuel-cladding mechanical interactions are negligible throughout the fuel’s operational cycle under both normal and high temperature accident scenario conditions. Although higher fidelity 3D time-dependent modeling and simulations are planned, the following may be concluded presently. The MARVEL fuel element maintains its geometric stability and structural integrity during the most extreme accident scenarios predicted for the MARVEL reactor. The hoop stress during the unprotected loss of flow accident reaches about -1.27 MPa; this negative stress indicates that it is compressive rather than tensile. The compressive stress is a result of the NaK pressure on the outside of the fuel element, caused by the restricted thermal expansion of the NaK coolant, exceeding the internal pressures generated inside of the fuel element. The hoop stress generated in the cladding during the unprotected loss of coolant accident reaches a maximum of about approximately 10 MPa, which is nearly an order of magnitude less than the predicted yield strength of the cladding under high-temperature accident scenario conditions. Calculations were compared with results from high performance computational simulations using BISON and are in very close agreement. A conservative MARVEL fuel meat peak temperature limit of 900 °C is recommended presently, which is about 180 °C higher than the peak fuel temperature predicted to occur during the most extreme accident. Based on the known properties and behavior of the MARVEL fuel element, the fuel successfully meets its design and safety requirements under normal and most extreme accident conditions with a large safety margin.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

A method for assessing economic, environmental, and reliability tradeoffs of interregional transmission connecting ERCOT (the Texas grid) to the eastern and western grids

Reliable development of the power grid is an evolving concern for humanity due to extreme weather that frequently threatens power sector infrastructure. The state of Texas is a uniquely structured testbed for grid planners to study when looking for solutions to development, innovation, and overcoming such challenges. Because of its size and islanded structure, Texas is small enough to model, but big enough to matter. Texas is a global leader in energy production, energy consumption, and maintains an unusually diverse fuel mix. In addition, the state has experienced winter freezes, heat waves, wind storms, droughts and floods that have threatened power sector infrastructure or caused recent blackouts and calls for demand side conservation. One of the most devastating of these events was the North American winter storm, dubbed “Winter Storm Uri” by the Weather Channel, that froze the region in February 2021 and led to an extended power outage event that put the majority of Texan residents in darkness for days. While preparing to avoid such outage events in the future, various tools have been proposed to improve grid reliability, including energy efficiency, demand response, and distributed energy resources. An additional option would be to develop interregional transmission that connects the Texas grid to other national grids. To assess the merits of this idea, we developed a novel, universally-applicable and internationally-relevant framework to study how the Texas grid would evolve alongside access to various interregional ties. This method allows us to stress the synthetic grid structure and analyze how it would respond to the shock of a simulated winter storm event. Our method leverages open-source modeling tools, such as PowerGenome, pyGRETA, and GenX to synthesize unique zonal grid data, construct a consolidated network of model regions, and simulate different developmental pathways of capacity expansion and operational dispatch. We demonstrate our method with an analysis connecting the Electric Reliability Council of Texas (ERCOT), the grid that serves most of Texas, the Western Electricity Coordinating Council (WECC), the grid that serves the western half of the contiguous U.S., and the Eastern Interconnect, the grid that serves the eastern half of the contiguous U.S. Our results indicate that the cost-optimal capacity of interregional transmission connecting the ERCOT grid to other grids lies between 9–13 GW assuming baseline conditions. Building this amount of connecting capacity in one or multiple directions lowers the costs and emissions of development and operation by up to $16 billion and 257 million metric tonnes (MMT) respectively. Additionally, our results show that the interregional connections between ERCOT and other national grids reduce the amount of total load shed required through mild winter storm events. However, our results also show that there is a threshold of very extreme winter storm conditions, spanning multiple service areas, above which the connections exacerbate resource adequacy problems. Therefore, the results indicate that the connections need to be carefully planned alongside the rest of the grid infrastructure to avoid over-reliance on specific resources or technology options.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

OC7 phase I: CFD investigation of viscous forces on rectangular members of semisubmersibles

In Phase I of the Offshore Code Comparison Collaboration 7 (OC7) project, we aim to address key questions related to the modeling of viscous hydrodynamic loads on floating platforms. Here, this paper presents results from a collaborative effort within Phase I of the OC7 project. We investigate viscous forces acting on the rectangular members of semisubmersible platforms using high-fidelity computational fluid dynamics (CFD) simulations. The parametric CFD analyses include two-and-a-half-dimensional (2.5D) simulations of pontoon sections of two different semisubmersible platform designs (the VolturnUS-S and the INO WINDMOOR) and full three-dimensional (3D) structure simulations of the VolturnUS-S platform. We extract both the added mass and drag coefficients that are typically used in engineering-level hydrodynamic models from the 2.5D forced oscillation simulations with various Keulegan-Carpenter numbers and compare them to the coefficients from the 3D forced oscillation simulations. We observe consistent variation of the coefficients along the pontoon legs due to 3D platform and free-surface effects. With the breadth of the simulations and component-level detailed analyses, this work establishes a useful database and several empirical formulas for selecting hydrodynamic coefficients when modeling semisubmersible platforms in engineering tools.

13 HYDRO ENERGY↗

DOE OSTI2024

Exploring dielectric properties in atomistic models of amorphous boron nitride

Abstract We report a theoretical study of dielectric properties of models of amorphous Boron Nitride, using interatomic potentials generated by machine learning. We first perform first-principles simulations on small (about 100 atoms in the periodic cell) sample sizes to explore the emergence of mid-gap states and its correlation with structural features. Next, by using a simplified tight-binding electronic model, we analyse the dielectric functions for complex three dimensional models (containing about 10.000 atoms) embedding varying concentrations of sp 1 , sp 2 and sp 3 bonds between B and N atoms. Within the limits of these methodologies, the resulting value of the zero-frequency dielectric constant is shown to be influenced by the population density of such mid-gap states and their localization characteristics. We observe nontrivial correlations between the structure-induced electronic fluctuations and the resulting dielectric constant values. Our findings are however just a first step in the quest of accessing fully accurate dielectric properties of as-grown amorphous BN of relevance for interconnect technologies and beyond.

Materials Science↗

DOE OSTI2025

Cross slip of extended dislocations in face-centered cubic metals through phase-field modeling

Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Atmospheric Modeling and Denoising for Millimeter-Wave Line Intensity Mapping

Line-intensity mapping (LIM) offers a promising approach to mapping large-scale cosmic structure, and the greatest obstacle for ground-based observations at millimeter wavelengths is foreground contamination from atmospheric emission. In this work, we present a simulation and denoising framework designed to isolate and subtract atmospheric fluctuations from LIM data, modeled after the instrument parameters of the South Pole Telescope Summertime Line Intensity Mapper (SPT-SLIM). We generate mock observations spanning 125-175 GHz containing cosmic signals, precipitable water vapor screens, ice crystal fluctuations, and photon noise. We then implement a spatial-spectral atmospheric removal pipeline combining per-pixel linear template regression with a two-dimensional Fourier-domain filter. The framework is evaluated under simulated conditions in the South Pole and the Atacama Desert across three key metrics: cosmic signal preservation, foreground subtraction efficiency, and instrument noise injection. Our pipeline achieves atmospheric suppression at large spatial scales, and these results establish a physically grounded foundation for atmosphere removal in ground-based LIM data collection.

Saye, Laney [UC, Berkeley (main)] (ORCID:000900078↗

DOE OSTI2026

Atmospheric Modeling and Denoising for Millimeter-Wave Line Intensity Mapping

Saye, L. K. [UC, Berkeley (main)] (ORCID:000900078↗

DOE OSTI2025

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗

DOE OSTI2024

Platform for Integrated Land use And Transportation Experiments and Simulation (PILATES) v1.0

PILATES allows for flexibly and at-scale coupling of multiple models to allow for multi-scale and multi-resolution simulation of regional-scale transport networks. In particular, it couples the MATSim-derived transportation modeling framework for Behavior, Energy, Autonomy and Mobility (BEAM) with other models operating at different time scales. Rather than tightly coupling supply and demand models using shared agents and memory within the same software process, PILATES orchestrates different model runs in a containerized framework. This structure requires passing information from the demand models to BEAM in the format of a synthetic population and agent plans, and from BEAM to the demand models in terms or origin/destination tables (also known as "skims"). This allows it to take advantage of the behavioral sophistication of existing activity-based models as well as the reinforcement learning structure of MATSim replanning and adopted by BEAM, in a way that requires minimal changes to existing models. It also takes advantage of the computational performance of BEAM, which allows for simulations with millions of agents to complete in reasonable time as well as allowing for detailed mechanistic simulation of the operation of on-demand modes.

Needell, Zachary↗

DOE OSTI2026

A generative machine learning model for designing metal hydrides applied to hydrogen storage

Developing new metal hydrides is a critical step toward efficient hydrogen storage in carbon-neutral energy systems. However, existing materials databases, such as the Materials Project, contain a limited number of well-characterized hydrides, which constrains the discovery of optimal candidates. This work presents a framework that integrates causal discovery with a lightweight generative machine learning model to generate novel metal hydride candidates that may not exist in current databases. Using a dataset of 450 samples (270 training, 90 validation, and 90 testing), the model generates 1000 candidates. After ranking and filtering, six previously unreported chemical formulas and crystal structures are identified, four of which are validated by density functional theory simulations and show strong potential for future experimental investigation. Overall, the proposed framework provides a scalable and time-efficient approach for expanding hydrogen storage datasets and accelerating materials discovery.

generative model↗

DOE OSTI2025

Micropolar deep material network

This study extends the Deep Material Network (DMN), a physics-informed machine learning framework, to predict the homogenized mechanical response of composite materials with micropolar (Cosserat-type) constitutive behavior. This extension incorporates microstructure-dependent size effects, enabling accurate, efficient, and size-aware predictions for composites with complex internal architectures. While traditional, direct numerical simulation micropolar models effectively capture size effects by introducing extra local degrees of freedom, they bring significant computational challenges, particularly for multiscale analyses relevant to engineering applications. The micropolar DMN developed in this paper achieves high accuracy while significantly reducing computation time compared to micropolar direct numerical simulations. This advancement enables multiscale analyses and parameter studies that were previously impractical, such as high-cycle fatigue simulations and comprehensive investigations of internal length scale effects notably in size-dependent plastic response and the optimization of lattice structures. By uniting microstructure-sensitive modeling, physics-driven learning, and scalable surrogate modeling, the micropolar DMN paves the way for accelerated material design, large-scale parametric studies, and the reliable incorporation of size-dependent effects across a wide range of engineering applications, including optimization and next-generation composite design.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Electron-Only Magnetic Reconnection and Inverse Magnetic-Energy Transfer at Subion Scales

We derive, and validate numerically, an analytical model for electron-only magnetic reconnection applicable to strongly magnetized plasmas. Our model predicts subion-scale reconnection rates significantly higher than those pertaining to large-scale reconnection, aligning with recent observations and simulations. Here, we apply this reconnection model to the problem of inverse magnetic energy transfer at subion scales. We derive time-dependent scaling laws for the magnetic energy decay and the typical magnetic structure dimensions that differ from those previously found in the magnetohydrodynamics regime. These scaling laws are validated via two- and three-dimensional simulations, demonstrating that subion-scale magnetic fields can reach large, system-size scales via successive coalescence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Cation Vacancies Enable Anion Redox in Li Cathodes

Conventional Li-ion battery intercalation cathodes leverage charge compensation that is formally associated with redox on the transition metal. Employing the anions in the charge compensation mechanism, so-called anion redox, can yield higher capacities beyond the traditional limitations of intercalation chemistry. Here, we aim to understand the structural considerations that enable anion oxidation and focus on processes that result in structural changes, such as the formation of persulfide bonds. Using a Li-rich metal sulfide as a model system, we present both first-principles simulations and experimental data that show that cation vacancies are required for anion oxidation. First-principles simulations show that the oxidation of sulfide to persulfide only occurs when a neighboring vacancy is present. To experimentally probe the role of vacancies in anion redox processes, we introduce vacancies into the Li 2 TiS 3 phase while maintaining a high valency of Ti. When the cation sublattice is fully occupied and no vacancies can be formed through transition metal oxidation, the material is electrochemically inert. Upon introduction of vacancies, the material can support high degrees of anion redox, even in the absence of transition metal oxidation. The model system offers fundamental insights to deepen our understanding of structure–property relationships that govern reversible anion redox in sulfides and demonstrates that cation vacancies are required for anion oxidation, in which persulfides are formed.

25 ENERGY STORAGE↗

DOE OSTI2025

Active-noise-induced dynamic clustering of passive colloidal particles

Active fluids generate spontaneous, often chaotic mesoscale flows. Harnessing these flows to drive soft materials embedded within an active fluid into structures with controlled length scales and lifetimes is a key challenge at the interface between the fields of active matter and nonequilibrium self-assembly. Here, we present a simple and efficient computational approach to model soft materials advected by active fluids, by simulating particles moving in a spatiotemporally correlated noise field. To illustrate our approach, we simulate the dynamical self-organization of repulsive colloids within such an active noise field. The colloids form structures whose sizes and dynamics can be tuned by the correlation time and length of the active fluid, and range from small rotating droplets to clusters with internal flows and system-spanning sizes that vastly exceed the active correlation length. Our results explain how the interplay between active fluid time and length scales and emergent driven assembly can be used to rationally design functional assemblies. More broadly, our approach can be used to efficiently simulate diverse active fluids and other systems with spatiotemporally correlated noise.

Brownian dynamics↗

DOE OSTI2025

Perspectives on Systematic Cloud Microphysics Scheme Development With Machine Learning

Cloud microphysics—the collection of processes that govern the small‐scale formation, evolution, and interactions of liquid droplets and ice crystals in clouds and precipitation—remains a major source of uncertainty in weather and climate models. Although too small in scale to be explicitly resolved in any large‐eddy simulation, weather, or climate model, the representation of cloud microphysical processes has significant impact at the climate scale. Current microphysical schemes are limited by both parametric uncertainty, linked to uncertainty in physical parameter values, and structural uncertainty, arising from incomplete physical understanding of the processes at play or approximations made for computational efficiency. Recent advances in the application of machine learning (ML) to the physical sciences show significant potential for minimizing these limitations by leveraging high‐fidelity simulations and observations. Here we outline the challenges that must be addressed to apply ML toward cloud microphysics scheme development. This perspectives paper synthesizes recent progress in using data‐driven methods, including ML, to improve cloud microphysics parameterizations and highlights opportunities to address key uncertainties. We discuss the roles of aleatoric (irreducible, or statistical) and epistemic (reducible, or systematic) errors in contributing to microphysics parameterization uncertainty. ML can leverage observations to improve microphysical schemes via bottom‐up and top‐down constraints. Methods such as differentiable programming and ML‐enhanced sampling strategies and the creation of large scale benchmark data sets promise to bridge the gap between observations and models and to improve the consistency of cloud microphysical representation across temporal and spatial scales.

Lamb, Kara D. [Columbia Univ., New York, NY (Unite↗

DOE OSTI2026

Achieving high rate performance in hybrid pristine-recycled cathodes using model-informed electrode designs

Direct recycling lithium-ion battery cathodes, a process that retains the engineered oxide structures from end-of-life materials, presents a cost-effective and energy-efficient alternative to other battery recycling methods. However, while direct-recycled cathodes have demonstrated performance comparable to that of pristine materials at low cycling rates, their high-rate performance remains uncertain. Morphology changes in cathode particles, a main mode of degradation, directly impact rate performance by limiting surface kinetics and solid-phase diffusion. If direct recycling processes do not sufficiently restore pristine-like morphologies, the recycled materials may retain structural defects that hinder high-rate performance. The present work uses a physics-based pseudo-2D model to simulate hybrid electrodes with pristine and artificially “aged/recycled” NMC materials to investigate potential impacts of incorporating performance-limited aged cathode materials into cells. The study highlights how differences in transport and kinetic properties can influence rate capabilities in mixed electrodes — particularly in high-loading cells in high-demand applications. However, model results also reveal a possible mitigation strategy via dual-layer electrode architectures with lower-performing materials positioned near the current collector. Simulations of 4.0 mAh cm −2 cells cycled at 4C using a dual-layer architecture provided approximately 5%–30% more capacity in constant-current protocols compared to homogeneously blended electrode architectures with the same loadings and mixed-material compositions. These findings highlight the importance of strategic electrode design in minimizing potential performance losses and facilitating the integration of recycled materials into high-performance batteries, advancing sustainable and cost-effective battery manufacturing.

25 ENERGY STORAGE↗

DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Subseasonal Tropical Convection Characteristics in the Energy Exascale Earth System Model Version 2

Accurate simulation of subseasonal tropical moist convection remains a key challenge for Earth system models. The difficulties stem from the reliance of cumulus cloud processes on model parameterizations and the need to represent the multiscale nature of interactions among clouds, radiation, moisture, circulation, and surface energy fluxes. Equatorial convection drives circulation anomalies that can affect weather patterns and extremes globally, motivating efforts to better understand and simulate these tropical disturbances. Here, a detailed review of subseasonal tropical convective behavior as simulated in the Energy Exascale Earth System Model version 2 (E3SMv2) is presented, with comparison to its predecessor version 1 (E3SMv1) and reference data sets. Model structural changes to the deep convective trigger function and surface fluxes, along with parametric tuning of the cloud and microphysics schemes, together result in an improved depiction of organized tropical convection across scales. In particular, E3SMv2 exhibits a more realistic Madden‐Julian oscillation (MJO) and low‐frequency Kelvin waves—owing to a sharper time mean equatorial meridional moisture gradient and improved convection‐circulation coupling —as well as a better depiction of MJO Northern Hemisphere teleconnections. Despite these improvements, subseasonal precipitation variance continues to be strongly underestimated in E3SMv2. Use of a cloud plume model also reveals that the coupling between daily averaged tropical precipitation and lower tropospheric instability in E3SM is inconsistent with observations, a bias that could potentially impact the simulation of intraseasonal disturbances.

54 ENVIRONMENTAL SCIENCES↗