Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale Modeling”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Tethys Water Demand Data

U.S. water demand varies sharply by sector and region as land use, population, weather patterns, and economic activity co-evolve. High-resolution water demand data is required to capture these dynamics, support integrated energy-water-land modeling, and local-to-regional water scarcity assessments. This dataset contains gridded (1/8 degree), monthly, multi-sector water demand dataset for the contiguous United States (CONUS) covering 1980-2100 across eight future scenarios of human-Earth system change. The dataset covers irrigation, thermoelectric, municipal (public-supply and domestic), livestock, manufacturing, and mining demands, separately for withdrawals and consumption, and includes per-cell renewable vs. non-renewable water source attributions. The dataset is validated against the latest USGS 2010-2020 water-use data for the three largest water demand sectors (Domestic, Electricity, and Irrigation), with correlations ranging from 0.73-0.95 at the HUC6 scale. The two datasets largely agree on an aggregate basis with per-sector bias falling within +/-7%, but they disagree on the spatial allocation of water with individual HUC6 basins having normalized RMSE from 68-171% and median absolute percent difference from 37-86%. This dataset advances prior global products by combining state-resolved sectoral demands from GCAM-USA, future power-plant siting from the CERF model, and scenario-consistent high-resolution climate and population forcing data across the eight scenarios.

GCAM-USA↗
DOE OSTI2025

Investigating Radiation-Induced Actinide Species in Solution

Our fundamental understanding of actinide radiation-induced redox chemistry is crucial to nuclear fuel cycle development, due to the unavoidable exposure of these elements to ionizing radiation fields, both inherent and from in-process applications. Plutonium (Pu) and americium (Am) both possess multiple oxidation states, the careful manipulation of which are essential in the study and utilization of their rich chemistry, developing new nuclear technologies, and securing the long-term sustainability of nuclear power. However, knowledge in this area is far from complete. Consequently, we have studied the radiation-induced chemistry of both Pu and Am through a variety of techniques. Temperature-controlled electron pulse radiolysis has been used to study Am for the first time, determining the feasibility of Am redox reactions under used nuclear fuel reprocessing conditions, (e.g. nitric acid, non-ambient temperature). Additionally, we have developed an experimentally evaluated multi-scale computer model for the prediction of gamma radiation-induced Pu redox chemistry due to radiolysis and disproportionation reactions in concentrated nitric acid solutions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigating Radiation-Induced Actinide Species in Solution

Our fundamental understanding of actinide radiation-induced redox chemistry is crucial due to their unavoidable exposure to ionizing radiation fields, both inherent and from in-process applications. Plutonium (Pu) and americium (Am) both possess multiple oxidation states, the careful manipulation of which are essential in the study and utilization of their rich chemistry, developing new nuclear technologies, and securing the long-term sustainability of nuclear power. However, knowledge in this area is far from complete. Consequently, we have studied the radiation-induced chemistry of both Pu and Am through a variety of techniques. Temperature-controlled electron pulse radiolysis has been used to study Am for the first time, determining the feasibility of Am redox reactions under used nuclear fuel reprocessing conditions, (e.g. nitric acid, non-ambient temperature). Additionally, we developed an experimentally evaluated multi-scale computer model for the prediction of gamma radiation-induced Pu redox chemistry due to radiolysis and disproportionation reactions in concentrated nitric acid solutions.

actinide↗
DOE OSTI2024

A mean field homogenization model for the mechanical response of ceramic matrix composites

SiC/SiC composites offer exceptional mechanical stability at high temperatures and under irradiation. These ceramic matrix composites are therefore strong candidate materials for future nuclear energy applications. Their mechanical response, which exhibits pseudo-plasticity, is mediated by matrix cracking, fiber debonding, and fiber pull-out due to slip. Here, this study introduces a mechanistic model for the behavior of unidirectionally reinforced SiC/SiC composites. Specifically a mean field homogenization approach is proposed to account for all deformation and degradation modes during mechanical deformation. The homogenization scheme relies on a Mori Tanaka method that is extended to consider the effects of the coating’s elasto-plastic response on the development of micromechanical fields. Further, the model proposed introduces a method to effectively account for the role of localized damage (i.e., cracks) on mechanical fields within both the fiber and the matrix. Upon validating the model against experimental data, the roles of interface sliding, coating dimensions and intrinsic elastic response, as well as of microstructure (e.g. porosity, fiber volume fraction) are discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Computational Advances in Ionic Liquid Applications for Green Chemistry: A Critical Review of Lignin Processing and Machine Learning Approaches

The valorization and dissolution of lignin using ionic liquids (ILs) is critical for developing sustainable biorefineries and a circular bioeconomy. This review aims to critically assess the current state of computational and machine learning methods for understanding and optimizing IL-based lignin dissolution and valorization processes reported since 2022. The paper examines various computational approaches, from quantum chemistry to machine learning, highlighting their strengths, limitations, and recent advances in predicting and optimizing lignin-IL interactions. Key themes include the challenges in accurately modeling lignin’s complex structure, the development of efficient screening methodologies for ionic liquids to enhance lignin dissolution and valorization processes, and the integration of machine learning with quantum calculations. These computational advances will drive progress in IL-based lignin valorization by providing deeper molecular-level insights and facilitating the rapid screening of novel IL-lignin systems.

09 BIOMASS FUELS↗
DOE OSTI2025

Implications of AI and Data Center Service Demand on Energy, Water, and Economic Systems

This talk focuses on the challenges posed by the rapid growth of data centers (DCs), which rely heavily on electricity, water, land, and mineral resources, and have significant interactions with economic growth. The U.S. leads in DC deployment, but understanding how resource demands and constraints impact DC expansion requires strategic multisector analysis. To address this, the talk introduces a new conceptual framework and modeling capability within GCAM-USA, allowing integrated analysis of DC service demand, electricity, and direct and indirect water needs under varied scenarios. Insights include state-level water and electricity trends for DCs, uncertainty bounds, and the broader implications for resources and economics.

Data Centers↗
DOE OSTI2026

Active learning of a crystal plasticity flow rule from discrete dislocation dynamics simulations

Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Low-cost fin-tube heat exchanger design for building thermal energy storage using phase change material

Heat transfer in phase change materials (PCMs) is complex because the melting and freezing fronts change as functions of stored or released heat. In prior attempts to optimize heat exchangers (HXs) in one or two dimensions, complex geometry has often been used to maximize the melt and freeze front area. This complex geometry is difficult and hence expensive to construct. This paper proposes a multiple-scale 3D finite element modeling approach to design fin-tube HXs for low-cost latent thermal energy storage applications. Further, the optimal fin and tube designs were determined at three scales (unit-scale, medium-scale, and large-scale) by modeling the melt and freeze front in three dimensions and using measured bulk thermal properties. The finite element model was validated by comparing it with the experimental data for a referenced design of a similar type. The results indicate that commercially available organic PCMs with low conductivity (<0.3 W/m·K) can have charge and discharge times appropriate for building thermal energy storage (i.e., 4–5 h) with fin-tube HX designs at costs <$26/kWh, even when the temperature difference (5.56°C) between the heat transfer fluid and the PCM phase change temperature is small. However, as the HX increases in length, the temperature reduction along the tube limits some larger-scale designs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Development of Machine-Learned Interatomic Potentials to Predict Structure, Transport, and Reactivity in Platinum-Based Fuel Cells

Machine-learned interatomic potentials (MLIPs) have rapidly progressed in accuracy, speed, and data efficiency in recent years. However, training robust MLIPs in multicomponent systems remains a challenge. In this work, we train an MLIP to describe hydrated Nafion ionomers and platinum catalysts, which are important components of fuel cells, by constructing a diverse training set to describe the bulk polymer and interfacial catalyst–polymer interactions well. We use our trained MLIP to study the properties of the platinum–Nafion system, including polymer structure, proton mobility in a bulk Nafion polymer and near a platinum-Nafion interface, and reactions near and far from the interface, finding excellent results for structure and reactions contained within our training set. Transport seems to be well described, with both vehicular transport and Grotthuss hopping captured, although converged calculations of diffusivities were not computed because they require calculations of tens of nanoseconds that are challenging with current state-of-the-art MLIPs. The combined insights that this model provides can be leveraged to optimize fuel cell performance, and the approach can be applied to other chemical processes and devices where structure, transport, and reactivity all contribute to the overall observed performance.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2025

Computational Modeling to Advance Novel Medical Isotopes for Radiotheranostics: A DOE-NIH Joint Workshop Executive Summary

The DOE-NIH Joint Workshop on Computational Modeling to Advance Novel Medical Isotopes for Radiotheranostics, held on September 27, 2024, brought together experts from government, academia, and industry to address critical challenges in radionuclide production and clinical translation. Here, the workshop emphasized interdisciplinary collaboration, particularly between the Department of Energy (DOE) and the National Institutes of Health (NIH), to strengthen the domestic isotope supply, streamline regulatory pathways, and further integrate computational tools into radiopharmaceutical therapy (RPT). Key discussions explored the role of AI-driven modeling, machine learning, and digital twin technologies in optimizing dosimetry, dynamically personalizing treatments, and reducing time to clinical adoption. Advances in predictive computational modeling were highlighted as essential for improving radionuclide yield, purity, and synthesis efficiency. Regulatory considerations and equitable access were central themes, with participants advocating for harmonized global standards, adaptive trial designs, and expanded infrastructure for clinical implementation. DOE computational and production infrastructure was emphasized. Future priorities identified include increased investment in radionuclide production infrastructure, expanded workforce development in radiopharmaceutical sciences and computational modeling, and the creation of robust public-private partnerships. The workshop concluded that continued strategic collaboration and sustained resources will be vital for advancing next-generation radiotheranostics, ensuring safe and effective therapies accessible to all patients.

digital twins↗
DOE OSTI2025

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗
DOE OSTI2026

A scalable framework for efficient coupling of thermal and microstructural simulations in additive manufacturing

Predicting microstructure evolution in metal additive manufacturing (AM) is important for process optimization, but spatiotemporal scale disparities between thermal transport and microstructure evolution create significant challenges for efficient data transfer between simulation codes. To address this, we present Stork, a scalable framework for coupling thermal and microstructural simulations. Stork uses a sparse data representation to identify and store active solidification sub-volumes, enabling highly parallel quad-linear interpolation from coarse thermal grids to fine microstructure grids without large intermediate storage. We demonstrate the framework by coupling the semi-analytic heat transfer code 3DThesis with the time-parallel cellular automata code Toucan. This approach achieves over two orders of magnitude reduction in data generation time and file size compared to prior workflows. Numerical studies show that quad-linear interpolation preserves grain morphology and crystallographic texture in laser powder bed fusion (LPBF) simulations for coarsening ratios up to 16. Overall, Stork provides a scalable pathway for high-throughput, component-scale AM simulations on modern high-performance computing systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Multiscale analysis of large twist ferroelectricity and swirling dislocations in bilayer hexagonal boron nitride

With its atomically thin structure and intrinsic ferroelectric properties, heterodeformed bilayer hexagonal boron nitride (hBN) has gained prominence in next-generation non-volatile memory applications. However, studies to date have focused almost exclusively on small-twist bilayer hBN, leaving the question of whether ferroelectricity can persist under small heterostrain and large heterodeformation entirely unexplored. In this work, we establish the crystallographic origin of ferroelectricity in bilayer hBN configurations heterodeformed relative to high-symmetry configurations such as AA-stacking and 21.786789° twisted configurations (Σ7), using Smith normal form bicrystallography. We then demonstrate out-of-plane ferroelectricity in bilayer hBN across configurations vicinal to both the AA and Σ7 stackings. Atomistic simulations reveal that AA-vicinal systems support ferroelectricity under both small twist and small strain, with polarization switching in the latter governed by the deformation of swirling dislocations rather than the straight interface dislocations seen in the former. For Σ7-vicinal systems, where existing interatomic potentials underperform particularly under extreme out-of-plane compression, we develop a density-functional-theory-informed continuum framework—the bicrystallography-informed frame-invariant multiscale (BFIM) model, which captures out-of-plane ferroelectricity in heterodeformed configurations vicinal to Σ7 stacking. Interface dislocations in these large heterodeformed bilayer configurations exhibit markedly smaller Burgers vectors compared to interface dislocations in small-twist and small-strain bilayer hBN. The BFIM model reproduces experimental results and provides a powerful, computationally efficient framework for predicting ferroelectricity in large-unit-cell heterostructures where atomistic simulations are prohibitively expensive.

Ahmed, Md Tusher [Univ. of Illinois at Urbana-Cham↗