Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Simulations”

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

Tardigrade-examples V0.1.0

Tardigrade-examples is a repository of computational workflows that exercise the Tardigrade software package. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling. These capabilities include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This repository investigates continuum upscaling of various direct numerical simulations (DNS) conducted in Abaqus finite element (FE), Ratel FE, and GEOS material point method (MPM) software. Verification of the upscaling workflow is first investigated by considering DNS of trivial stress states for homogeneous materials, results of which indicate that classical continuum behavior is recovered as expected. DNS of heterogeneous materials are then considered.

Allard, Thomas↗
DOE OSTI2026

Towards robust laser beam propagation in atmospheric turbulence

High-fidelity optical propagation through the atmosphere is essential for free-space optical technologies, including laser-based remote sensing and optical communication. However, atmospheric turbulence severely distorts beams and compromises system performance. In this work, we employ hypergeometric-Gaussian (HyGG) vortex beams as probes to characterize and mitigate atmospheric turbulence. Using over 250,000 experimental and simulated frames, we show that refining the power spectrum density (PSD) can reduce numerical prediction errors by up to 79.8%. Concurrently, experimental observations supported by numerical simulations demonstrate that HyGG beams exhibit superior turbulence resilience across multiple metrics compared to conventional Gaussian beams, particularly in their ability to withstand over 5 times stronger turbulence while maintaining similar intensity fluctuations. These dual investigations, on both turbulence mitigation and robust beam solutions, converge to form a unified strategy for enhancing free-space optical system performance. Collectively, our findings provide new insights into light–turbulence interactions and highlight the practical utility of vortex beams under atmospheric conditions.

Zhang, Boyu↗
DOE OSTI2026

Final DOE-ASR Report for the Project “Using LASSO to bridge the gap between model and observations and to learn about atmospheric convection”

Atmospheric convection spans a wide range of spatial and temporal scales and involves complex interactions with the surrounding dynamic and thermodynamic environment, particularly over tropical continental regions. These processes remain a major source of uncertainty in weather and climate models, including persistent biases in the diurnal cycle of convective precipitation that directly affect estimates of climate sensitivity. Addressing these challenges requires the combined use of high-resolution observations and cloud-resolving modeling frameworks. In this context, the DOE Atmospheric Radiation Measurement (ARM) program’s Large-Eddy Simulation ARM Symbiotic Simulation and Observation (LASSO) activity provides a powerful platform that pairs comprehensive observations with numerical simulations to enable process-level understanding of atmospheric convection. Within this context, this Research and Development Partnership Pilot (RDPP) project was designed to initiate and expand DOE ARM/ASR research capacity at minority-serving institutions, while advancing scientific understanding of convective processes over the Amazon rainforest. Consistent with the RDPP mission, the project emphasized partnership development, training, and workforce capacity building alongside exploratory research activities. On the scientific side, the project produced two peer-reviewed journal articles, and one manuscript currently under review (see list in section 3.1). Together, these studies combine long-term ARM observations and cloud-resolving and convection-permitting modeling to investigate the environmental controls on the shallow-to-deep convective transition during the Amazon wet season. The results demonstrate the central role of early-day moisture preconditioning and large-scale dynamical forcing in regulating isolated deep convection, provide mechanistic insight into convective evolution, and establish physically informed modeling frameworks for future sensitivity experiments. These scientific outcomes are described in sections 2.1 to 2.3 and were disseminated in 8 conference presentations (see section 3.2) and 5 invited talks (see section 3.3), reflecting broad engagement with our community. Equally important, the project achieved its RDPP capacity-building objectives (see section 2.4). A sustained research partnership was established among the University of Maryland, Baltimore County (UMBC), Morgan State University (MSU), and Howard University (HU), and extended to include collaboration with Pacific Northwest National Laboratory (PNNL). The project organized multiple multi-day training events focused on ARM data, LASSO simulations, and quantitative analysis methods, directly engaging students, postdoctoral researchers, and faculty across institutions. These activities broadened participation in ASR research and led to independent adoption of LASSO workflows by students beyond the immediate project team. Finally, the project successfully positioned the participating institutions to pursue future DOE research. Preliminary scientific results, coupled with strengthened partnerships and technical capacity, enabled the submission of follow-on proposals to DOE ASR funding opportunities. In this way, the project fulfilled the RDPP goal of seeding durable research capacity and laying the foundation for larger-scale, sustained engagement with DOE ARM and ASR programs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Formation and Evolution Simulations of Saturn, Including Composition Gradients and Helium Immiscibility

In this paper, the formation of Saturn is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. Previous models are enhanced to include the dissolution of accreting planetesimals, composed of water ice, rock, and iron, in the gaseous envelope of the planet, leading to a nonuniform composition with depth. The immiscibility of helium in metallic hydrogen layers is also considered. The calculations start at a mass of 0.5 Earth masses and are extended to the present day. At 4.57 Gyr, the model, proceeding outward, has the following structure: (i) a central core composed of 100% heavy elements and molecules, (ii) a region with a decreasing heavy-element mass fraction, down to a value of 0.1, (iii) a layer of uniform composition with the helium mass fraction Y enhanced over the primordial value, (iv) a helium rain region with a gradient in Y, (v) an outer convective, adiabatic region with uniform composition in which Y is reduced from the primordial value, and (vi) the very outer layers where cloud condensation of the heavy elements occurs. Models of the distribution of heavy elements as a function of radius are compared with those derived to fit the observations of the Cassini mission, with rough qualitative agreement. The helium mass fraction in Saturn’s outer layers is estimated to be around 20%. Models are found that provide good agreement with Saturn’s intrinsic luminosity and radius.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Direct statistical simulation of the Lorenz96 system in model reduction approaches

Direct statistical simulation (DSS) of nonlinear dynamical systems bypasses the traditional route of accumulating statistics by lengthy direct numerical simulations by solving the equations that govern the statistics themselves. DSS suffers, however, from the curse of dimensionality as the statistics (such as correlations) generally have higher dimensions than the underlying dynamical variables. Here we investigate two approaches to reduce the dimensionality of DSS, illustrating each method with numerical experiments with the Lorenz96 dynamical system. The forms of DSS chosen here involve approximate closures at second and third order in the equal-time cumulants. We demonstrate significant reduction in computational effort that can be achieved without sacrificing the accuracy of DSS. The methods developed here can be applied to turbulent fluid and magnetohydrodynamical systems. Published by the American Physical Society 2025

Li, Kuan↗
DOE OSTI2026

Modeling glass degradation and release of radionuclides from vitrified waste for performance assessment simulations

The release of radionuclides initially encapsulated in a slowly degrading solid waste form and contained in an eventually corroding canister defines the source term for numerical simulations for the assessment of a geologic repository for high-level radioactive waste. While the details of waste degradation, canister corrosion, and dissolution and mobilization of the radionuclides in pore water include complex chemical reaction and transport processes that are coupled to the thermal, hydrological, microbiological, and mechanical conditions in the repository, the source-term model suitable for use in a numerical performance assessment model should be a defensible abstraction of these mechanisms. We developed a radiological source-term model and implemented it into a non-isothermal flow and transport simulator. While the proposed source-term model is applicable to various waste forms, canister systems, and disposal concepts, we specifically considered radionuclide releases from vitrified high-level waste placed in a cylindrical canister disposed in a deep vertical borehole repository. In this model, waste degradation is a function of temperature, and it can be adjusted to evaluate the influence of and propagate uncertainties in pH, passivation reactions, and chemical conditions as well as geometrical factors. The time-dependent, congruent release of safety-relevant radionuclides present in the decaying inventory is then calculated. Finally, the radionuclides are mobilized by diffusive and advective transport according to the thermo-hydraulic conditions prevailing in the near field of the repository, from where they migrate through the geosphere to the accessible environment. We examine the influence of the source-term model’s parameters on performance assessment calculations through sensitivity and uncertainty propagation analyses, identifying influential factors and confirming the upper bound of their impact. These considerations align with the overarching goal of repository design, which is to demonstrate that engineered and natural barriers can collectively delay radionuclide migration for timescales far exceeding human planning, thereby providing multiple, redundant barriers against environmental contamination.

iTOUGH2↗
DOE OSTI2026

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗
DOE OSTI2024

Quantifying Memory in Spin Glasses

Rejuvenation and memory, long considered the distinguishing features of spin glasses, have recently been proven to result from the growth of multiple length scales. This insight, enabled by simulations on the Janus II supercomputer, has opened the door to a quantitative analysis. We combine numerical simulations with comparable experiments to introduce two coefficients that quantify memory. A third coefficient has been recently presented by Freedberg et al. Here, we show that these coefficients are physically equivalent by studying their temperature and waiting-time dependence.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantitative Voltammetry Measurements of High-Concentration Actinides in Molten Chloride Fuel Salts

Molten chloride mixtures containing high concentrations of UCl 3 (i.e., >70 wt%) are candidates for next-generation molten salt reactor (MSR) fuel salts that enable new plant optimization and fuel cycle strategies. MSR licensing, however, is contingent on the development of effective material control and accounting (MC&A) practices for these salts. To close this gap, we have developed an electroanalytical approach that enables accurate in situ measurements of uranium concentrations in these highly loaded fuel salts. This methodology uses a multielectrode array voltammetry sensor combined with digital simulations that account for non-idealities such as uncompensated resistance effects. The theory-based numerical simulations of voltammetry responses yielded correction factors that are agnostic of electrochemical cell geometry and produced consistent voltammetry peak current results when correcting for varying amounts of ohmic resistance. In doing so, we have demonstrated accurate concentration measurements for UCl 3 in LiCl-KCl-UCl 3 across a range of mass fractions between 55 and 72.3 wt% at 550 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Sensitivity of the Shallow‐To‐Deep Convective Transition to Moisture and Wind Shear in the Amazon

Deep convection is the primary influence on weather and climate in tropical regions. However, understanding and simulating the shallow‐to‐deep (STD) convective transition has long been challenging. Here, we conduct high‐resolution numerical simulations to assess the environmental controls on the evolution of isolated convection in the Amazon during the wet season. The large‐scale forcing derived through a constrained variational analysis approach for the GoAmazon2014/5 Experiment is used in the simulations. Through sensitivity experiments, we examine the relative importance of moisture and wind shear in controlling the shallow‐to‐deep convective transition for isolated convective events. Convection exhibits the greatest sensitivity to humidity within the lowest 1.5 km, where a 4 mm reduction in column water vapor nearly suppresses ice water formation on deep convective days. In contrast, a reduction in column water vapor in the free troposphere by a factor of two or more is necessary to produce a comparable impact on convection. Increasing low‐level wind speed from 6 to 9 m s -1 enhances afternoon deep convection, raising the cloud ice mixing ratio by approximately 25%. Conversely, upper‐level wind shear reveals the weakest correlation with daytime convection in our simulations. Our results help characterize the role of moisture and wind shear on the STD transition and our understanding of the underlying mechanisms.

environmental controls↗
DOE OSTI2025

Observation and Simulation of Methane Plumes During the Morning Boundary Layer Transition

Abstract Methane (CH 4 ) contributes significantly to global warming. However, accurate identification of CH 4 sources for reducing CH 4 emissions is often hampered by inadequate accuracy and spatiotemporal coverage of CH 4 detection, and lack of accurate CH 4 forward modeling used in top‐down inversion systems. In this study, a field experiment was conducted in Pampa, Texas using two CH 4 sensors (LI‐COR and OGI camera) to detect CH 4 releases. We investigated whether high‐resolution simulations using the Weather Research and Forecasting (WRF) model with greenhouse gases (WRF‐GHG) could accurately simulate the CH 4 plumes in the presence of evolving atmospheric boundary layer from sunrise to noon. CH 4 plumes showed substantial variation in time. At a release rate of ∼17.5 kg hr −1 , the maximum enhancement of CH 4 measured by LI‐COR was 2.6 ppm at sunrise (7:36 a.m.), 250 m from the release location. Within half an hour after sunrise, this enhancement decreased to 0.3–0.4 ppm. The enhancement was 0.2 ppm by 10:00 a.m. and further dropped to less than 0.1 ppm after 11:30 a.m. Due to the low temperature at sunrise, the OGI camera failed to detect the CH 4 plume. The WRF‐GHG large‐eddy simulation (LES) with 32 m grid spacing successfully reproduced these CH 4 enhancements. In situ measurements together with numerical simulations illustrate the impact of the transition from a stable boundary layer in the early morning to a convective boundary layer at noon on the dispersion of CH 4 plumes. Additionally, CH 4 plumes from a cattle farm in Oklahoma are briefly examined using the same modeling approach.

Hu, Xiao‐Ming [Center for Analysis and Prediction ↗
DOE OSTI2024

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗
DOE OSTI2025

Multiphysics simulation of recent experiments on alkali‐silica reaction expansion in reinforced concrete members

Alkali‐silica reaction (ASR) is an important degradation process that causes volumetric expansion and damage in concrete, and is affected significantly by the local temperature, moisture and stress conditions that often vary across the regions of a structure. Numerical simulation is essential to predict the progression and effects of ASR on the performance of structures. Because of the interactions between thermal and moisture transport and mechanical deformation, it is important for numerical models to represent all these physical phenomena and the coupling between them. Simulations of ASR in reinforced concrete (RC) structures are further complicated by the need to capture interactions between concrete and embedded reinforcing bars. Here, this paper describes the implementation of a scalable, coupled‐physics ASR model for simulating RC structures and assesses the ability of that model to predict ASR‐induced expansion in recent laboratory tests on RC block and beam specimens. These laboratory tests and the simulation approach were selected because of their applicability to RC structural‐scale simulations. This validation study helps builds confidence the ability of this approach to model ASR expansion in large, complex RC structures, which is a current high‐priority need.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Problem-tailored Simulation of Energy Transport on Noisy Quantum Computers

The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L = 12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Understanding Recombination Pathways of CdSeTe/CdTe Devices Using Temperature Dependence of Fill Factor

Understanding the temperature dependence of photovoltaic parameters is critical for identifying performance limiting mechanisms in CdTe-based solar cells. In this study, we investigate the temperature dependence of fill factor (FF(T)) in CdSeTe/CdTe solar cells using a combination of experimental current density-voltage (JV) measurements and numerical simulations. The measured FF(T) exhibits a nonlinear response with temperature, peaking at moderate temperature values (10–30°C). Through numerical modeling, we reveal that this peaking behavior in FF(T) arises when the dominant recombination mechanism transitions from bulk to interface with increasing temperature. This transition is primarily driven by low hole mobility (0.6–1.2 cm2/V·s), while electron mobility, interface recombination velocities, and carrier lifetime also contribute to the effect. Our findings demonstrate that FF(T) analysis is a promising diagnostic tool for distinguishing recombination pathways and guiding the optimization of thin-film CdTe solar cells under real-world conditions.

CdTe photovoltaics↗
DOE OSTI2025

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗