Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Dispersion and Aging of Volcanic Aerosols after the La Soufriere Eruption in April 2021

Volcanic aerosols change the atmospheric composition and thereby affect weather and climate. Aerosol dynamic processes such as nucleation, condensation, and coagulation modify the shape, size, and mass of aerosol particles, which influence their atmospheric lifetime and radiative properties. Nevertheless, most models omit these processes for ash particles. In this work, we explore the ash aerosol aging and sulfate production during the first 4 days following the 2021 La Soufrière (St. Vincent) eruption with the ICON-ART model (ICOsahedral Nonhydrostatic model with Aerosol and Reactive Trace gases). Online coupling of ICON-ART with a one-dimensional volcanic plume model calculates volcanic emission, which makes it possible to resolve the different eruption phases of the noncontinuous La Soufrière eruption. We compared our simulated aerosol distribution and composition with observations from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument, the Multiangle Imaging SpectroRadiometer (MISR) Research Aerosol (RA) Algorithm, and the Barbados Cloud Observatory (BCO). We show that online coupling is essential to adequately model the emissions and plume development close to the volcano. The modeled aerosol aging is in very good agreement with observations from MISR near the emission source and with CALIOP at larger distances. Furthermore, particle aging occurs faster in the troposphere than in the stratosphere due to the availability of water vapor and OH, but a layer of coated ash appears at the plume top due to faster oxidation of SO2 and lofting by aerosol-radiation interaction. This paper gives the first direct comparison of aerosol aging in volcanic eruption plumes between simulations and observations

volcanic aerosols↗
NASA NTRS2017

Ion Velocity Distributions in Dipolarization Events: Beams in the Vicinity of the Plasma Sheet Boundary

Using combined MHD/test particle simulations, we further explore characteristic ion velocity distributions in relation to magnetotail reconnection and dipolarization events, focusing on distributions at and near the plasma sheet boundary layer (PSBL). Simulated distributions right at the boundary are characterized by a single earthward beam, as discussed earlier. However, farther inside, the distributions consist of multiple beams parallel and antiparallel to the magnetic field, remarkably similar to recent Magnetospheric Multiscale observations. The simulations provide insight into the mechanisms: the lowest earthward beam results from direct acceleration at an earthward propagating dipolarization front (DF), with a return beam at somewhat higher energy. A higher-energy earthward beam results from dual acceleration, first near the reconnection site and then at the DF, again with a corresponding return beam resulting from mirroring closer to Earth. Multiple acceleration at the X line or the propagating DF with intermediate bounces may produce even higher-energy beams. Particles contributing to the lower energy beams are found to originate from the PSBL with thermal source energies, increasing with increasing beam energy. In contrast, the highest-energy beams consist mostly of particles that have entered the acceleration region via cross-tail drift with source energies in the suprathermal range.

Birn, J.↗
NASA NTRS2024

Modeling Impacts of Dust Mineralogy on Fast Climate Response

Mineralogical composition drives dust impacts on Earth's climate systems. However, most climate models still use homogeneous dust, without accounting for the temporal and spatial variation in mineralogy. To quantify the radiative impact of resolving dust mineralogy on Earth's climate, we implement and simulate the distribution of dust minerals (i.e., illite, kaolinite, smectite, hematite, calcite, feldspar, quartz, and gypsum) from Claquin et al. (1999) (C1999) and activate their interaction with radiation in the GFDL AM4.0 model. Resolving mineralogy reduces dust absorption compared to the homogeneous dust used in the standard GFDL AM4.0 model that assumes a globally uniform hematite volume content of 2.7 % (HD27). The reduction in dust absorption results in improved agreement with observation-based single-scattering albedo (SSA), radiative fluxes from CERES (the Clouds and the Earth's Radiant Energy System), and land surface temperature from the CRU (Climatic Research Unit) compared to the baseline HD27 model version. It also results in distinct radiative impacts on Earth's climate over North Africa. Over the 19-year (from 2001 to 2019) modeled period during JJA (June–July–August), the reduction in dust absorption in AM4.0 leads to a reduction of over 50 % in net downward radiation across the Sahara and approximately 20 % over the Sahel at the top of the atmosphere (TOA) compared to the baseline HD27 model version. The reduced dust absorption weakens the atmospheric warming effect of dust aerosols and leads to an alteration in land surface temperature, resulting in a decrease of 0.66 K over the Sahara and an increase of 0.7 K over the Sahel. The less warming in the atmosphere suppresses ascent and weakens the monsoon inflow from the Gulf of Guinea. This brings less moisture to the Sahel, which combined with decreased ascent induces a reduction of precipitation. To isolate the effect of reduced absorption compared to resolving spatial and temporal mineralogy, we carry out a simulation where the hematite volume content of homogeneous dust is reduced from 2.7 % to 0.9 % (HD09). The dust absorption (e.g., single-scattering albedo) of HD09 is comparable to that of the mineralogically speciated model on a global mean scale, albeit with a lower spatial variation that arises solely from particle size. Comparison of the two models indicates that the spatial inhomogeneity in dust absorption resulting from resolving mineralogy does not have significant impacts on Earth's radiation and climate, provided there is a similar level of dust absorption on a global mean scale before and after resolving dust mineralogy. However, uncertainties related to emission and distribution of minerals may blur the advantages of resolving minerals to study their impact on radiation, cloud properties, ocean biogeochemistry, air quality, and photochemistry. On the other hand, lumping together clay minerals (i.e., illite, kaolinite, and smectite), but excluding externally mixed hematite and gypsum, appears to provide both computational efficiency and relative accuracy. Nevertheless, for specific research, it may be necessary to fully resolve mineralogy to achieve accuracy.

Soil dust↗
NASA NTRSDate not supplied

Evolving the SpaceFOM: Lessons Learned and Future Development

For over fifty years, simulation has been a cornerstone technology for space missions. In many cases, one single simulator, developed by one team, is not enough to meet the simulation objectives. Astronaut crew training, docking, vehicle testing and complex mission planning all require several simulations to be connected and to interoperate, in order to achieve the goal of the simulation. To facilitate the development of interoperable simulations, an open standard, the SISO SpaceFOM was developed in collaboration between government, industry and academia. The standard was released in 2020. The initial release of SpaceFOM focuses on the core interoperability topics for space simulation: (i) specifying spatial data using several, well defined reference frames, (ii) managing time (including real-time, hard real-time, and as-fast-as-possible execution), (iii) composition of systems with subsystems such as space vehicles, and (iv) execution control, such as initializing, running, pausing, and terminating simulations. A number of requirements for the next version of SpaceFOM have been identified and proposed solutions are already under development. One key area is the object classes used for data exchange, where a richer set of classes have been proposed, based on practical use cases. For time management, the use of variable time resolutions and variable time steps is also proposed. For execution control, save and restore of simulations is considered as well as improved robustness, to handle failing components. A number of clarifications are also under consideration. This paper provides both a forward-looking view of the technical developments underway and a retrospective on the standard's evolution. The authors emphasize that prioritizing core challenges and delivering a practical first version were crucial steps. While not every idea was implemented in the initial release, the timely delivery of a functional standard enabled real-world applications, the growth of a user community, and the collection of valuable insights to shape future versions. The authors look forward to working with the next version of SpaceFOM together with the growing user community.

FOMs↗
NASA NTRS2001

Models of Small-Scale Patchiness

Patchiness is perhaps the most salient characteristic of plankton populations in the ocean. The scale of this heterogeneity spans many orders of magnitude in its spatial extent, ranging from planetary down to microscale. It has been argued that patchiness plays a fundamental role in the functioning of marine ecosystems, insofar as the mean conditions may not reflect the environment to which organisms are adapted. Understanding the nature of this patchiness is thus one of the major challenges of oceanographic ecology. The patchiness problem is fundamentally one of physical-biological-chemical interactions. This interconnection arises from three basic sources: (1) ocean currents continually redistribute dissolved and suspended constituents by advection; (2) space-time fluctuations in the flows themselves impact biological and chemical processes, and (3) organisms are capable of directed motion through the water. This tripartite linkage poses a difficult challenge to understanding oceanic ecosystems: differentiation between the three sources of variability requires accurate assessment of property distributions in space and time, in addition to detailed knowledge of organismal repertoires and the processes by which ambient conditions control the rates of biological and chemical reactions. Various methods of observing the ocean tend to lie parallel to the axes of the space/time domain in which these physical-biological-chemical interactions take place. Given that a purely observational approach to the patchiness problem is not tractable with finite resources, the coupling of models with observations offers an alternative which provides a context for synthesis of sparse data with articulations of fundamental principles assumed to govern functionality of the system. In a sense, models can be used to fill the gaps in the space/time domain, yielding a framework for exploring the controls on spatially and temporally intermittent processes. The following discussion highlights only a few of the multitude of models which have yielded insight into the dynamics of plankton patchiness. In addition, this particular collection of examples is intended to furnish some exposure to the diversity of modeling approaches which can be brought to bear on the problem. These approaches range from abstract theoretical models intended to elucidate specific processes, to complex numerical formulations which can be used to actually simulate observed distributions in detail.

McGillicuddy, D. J.↗
DOE OSTI2026

Superstructure Optimization of Waste Plastic Pyrolysis, Integrating Thermal, Catalytic, and Plasma Technologies with Machine Learning

Global plastic waste generation exceeds 430 million tonnes per year, yet fewer than 9% are recycled in the United States. Pyrolysis offers a chemical recycling route at scale, but existing techno-economic and life cycle assessments fix product yields to single pure polymers, producing economic and environmental outputs that break down when the feed composition changes. Here, we present a superstructure optimization framework that addresses this by embedding a composition-aware random forest yield predictor, trained on 566 pyrolysis experiments, within a full-scale process simulation. Product distributions update automatically as feed allocation shifts across four reactor chemistries: conventional thermal, catalytic (HZSM-5), thermal oxo-degradation, and nonequilibrium CO2 plasma. The optimal superstructure achieves minimum selling prices of −0.56 to −0.76/kg feed and global warming potentials of −0.276 to −0.322 kg CO2-eq/kg feed across four commodity price scenarios, confirming profitable, carbon-negative operation without tipping fees. Carbon abatement costs of $\$$0.46 to $\$$1.25/kg CO2-eq are competitive with direct air capture. Sensitivity analysis shows that the catalytic-plasma split fraction is the single largest driver of both economic and climate performance, while hydrocracking allocation in the wax upgrading stage is emission-neutral across the full variable range. Mixed plastic waste streams, evaluated as composition-variable feedstocks rather than pure resins, are profitable and carbon-negative across realistic market conditions. These results give a quantitative basis for reactor selection, circular economy investment, and policy design targeting chemical recycling on a large scale.

Life cycle assessment↗
DOE OSTI2025

Global Simulations of Phase State and Equilibration Time Scales of Secondary Organic Aerosols with GEOS-Chem

The phase state of secondary organic aerosols (SOA) can range from liquid through amorphous semisolid to glassy solid, which is important to consider as it influences various multiphase processes including SOA formation and partitioning, multiphase chemistry, and cloud activation. In this study, we simulate the glass transition temperature and viscosity of SOA over the globe using the global chemical transport model, GEOS-Chem. The simulated spatial distributions show that SOA at the surface exist as liquid over equatorial regions and oceans, semisolid in the midlatitude continental regions, and glassy solid over lands with low relative humidity. The predicted SOA viscosities are mostly consistent with the available measurements. In the free troposphere, SOA particles are mostly predicted to be semisolid at 850 hPa and glassy solid at 500 hPa, except over tropical regions including Amazonia, where SOA are predicted to be low viscous. Phase state also exhibits seasonal variation with a higher frequency of semisolid and solid particles in winter compared to warmer seasons. We calculate equilibration time scales of SOA partitioning (τ eq ) and effective mass accommodation coefficient (α eff ), indicating that τ eq is shorter than the chemical time step of GEOS-Chem of 20 min and αeff is close to unity for most locations at the surface level, supporting the application of equilibrium SOA partitioning. However, τ eq is prolonged and α eff is lowered over drylands and most regions in the upper troposphere, suggesting that kinetically limited growth would need to be considered for these regions in future large-scale model studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Stable Simulation of the Community Atmosphere Model Using Machine‐Learning Physical Parameterization Trained With Experience Replay

In recent years, machine learning (ML) models have been used to improve physical parameterizations of general circulation models (GCMs). A significant challenge of integrating ML models into GCMs is the online instability when they are coupled for long‐term simulation. We present a new strategy that demonstrates robust online stability when the physical parameterization package of an atmospheric GCM is replaced by a deep ML model. The method uses experience replay with a multistep training scheme of the ML model in which the model's own output at the previous time step is used in the training. Predicted physics tendencies in the replay buffer with the most recent errors in the training iterations are reused, making the ML model learn from its own errors. The training method reduces the gap between the offline and online environments of the ML model. The method is used to train the ML model as the physical parameterization of the Community Atmosphere Model (CAM5) with training data from the Multi‐scale Modeling Framework high resolution simulations. Three 6‐year online simulations of the CAM5 are carried out by using the ML physics package. The simulated spatial distributions of precipitation, surface temperature and zonally averaged atmospheric fields demonstrate overall better accuracy than that of the standard CAM5 and benchmark model even without the use of additional physical constraints or tuning. This work is the first to demonstrate a solution to address the online instability problem in climate modeling with ML physics by using experience replay.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Cosmological constraints using Minkowski functionals from the first year data of the Hyper Suprime-Cam

We use Minkowski functionals to analyse weak lensing convergence maps from the first-year data release of the Subaru Hyper Suprime-Cam (HSC-Y1) survey. Minkowski functionals provide a description of the morphological properties of a field, capturing the non-Gaussian features of the Universe matter-density distribution. Using simulated catalogues that reproduce survey conditions and encode cosmological information, we emulate Minkowski functionals predictions across a range of cosmological parameters to derive the best-fit from the data. By applying multiple scales cuts, we rigorously mitigate systematic effects, including baryonic feedback and intrinsic alignments. From the analysis, combining constraints of the angular power spectrum and Minkowski functionals, we obtain S8≡σ8Ωm/0.3=0.808−0.046+0.033 and Ωm=0.293−0.043+0.157⁠. These results represent a 40 per cent improvement on the S8 constraints compared to using power spectrum only. Minkowski functionals results are consistent with other two-point, and higher order statistics constraints using the same data, being in agreement with CMB results from the Planck S8 measurements. Our study demonstrates the power of Minkowski functionals beyond two-point statistics to constrain and break the degeneracy between Ωm and σ8⁠.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Adaptive Computing (AC) [SWR-24-106]

The Adaptive Computing (AC) software stack supports goal-based computing, for which a simulation workload is created on the fly adapting to the results of calculations. Application-specific code defines an objective, which may be to solve an optimization problem or to train a surrogate model with minimal uncertainty. Then, the AC driver decides where in the design parameter space to run simulations to best achieve that objective. This process is iterative and online; as new data is returned from simulations, the AC driver chooses new simulations to run. The AC driver can strategically run simulations on distributed hardware resources (including high performance computing machines, cloud resources, and edge devices) to maximize throughput and obey resource constraints.

Griffin, Kevin [National Renewable Energy Laborato↗
DOE OSTI2025

DISTRI: Distributed Multi-Facility HPC Simulator (DISTRI) v2.1

DISTRI is an advanced network simulator designed for multi-facility computational infrastructures with agentic behavior. It simulates HPC facilities where computational resources act as autonomous agents, making intelligent decisions about job scheduling, load balancing, and resource allocation. The simulator focuses on developing and testing decentralized algorithms that promote resilience and efficiency in multi-facility environments. Key Features: - Agentic Resource Behavior: Processors and DTNs act as autonomous agents with decision-making capabilities - Pheromone-Based Load Balancing: Decentralized load balancing inspired by ant colony optimization - Dual Topology Support: Mesh (normal operations) and Dumbell (network testing) topologies - Comprehensive TCP Simulation: Realistic TCP implementations with multiple congestion control algorithms - Failure Resilience Testing: Processor failure simulation with automatic job reassignment - Extensive Visualization: Detailed performance analysis and metrics collection - Research-Ready: Designed for algorithm development and benchmarking

Bez, Jean Luca [Lawrence Berkeley National Laborat↗
DOE OSTI2026

Adaptive Computing (AC) (Open Source) [SWR-24-106]

The Adaptive Computing (AC) software stack supports goal-based computing, for which a simulation workload is created on the fly, adapting to the results of calculations. Application-specific code defines an objective, which may be to solve an optimization problem or to train a surrogate model with minimal uncertainty. Then, the AC driver decides where in the design parameter space to run simulations to best achieve that objective. This process is iterative and online; as new data is returned from simulations, the AC driver chooses new simulations to run. The AC driver can strategically run simulations on distributed hardware resources (including high performance computing machines, cloud resources, and edge devices) to maximize throughput and obey resource constraints.

Griffin, Kevin [National Laboratory of the Rockies↗
DOE OSTI2025

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗
DOE OSTI2025

AI Applications to Physics Experiments at Jefferson Lab

We survey how AI/ML is being deployed across Jefferson Lab's experimental and accelerator programs. In EPSCI, Hydra applies computer vision to automate real-time data-quality monitoring across all four experimental halls, replacing manual inspection of hundreds to thousands of histograms per shift. AIEC (AI Experiment Controls) uses ML to stabilize drift chamber gains and is now part of standard CEBAF production running, while AI Optimized Polarization (AIOP) targets autonomous control of polarized targets and photon beam angular alignment. In CASA, cavity fault classification models identify faulted cavities and trip types from waveform data with ~85% and ~78% agreement to labeled data, respectively, and are deployed in production; a separate effort applies LLMs and hybrid search to make the CEBAF operations logbook AI-ready. QCD-focused work includes transformer- and GAN-based generative models for particle-level event simulation, with distributed GAN training scaling studies on Polaris. Additional efforts span ML-on-FPGA for the EIC and a new Data Science Department coordinating anomaly detection, uncertainty quantification, and HPC-scalable ML lab-wide. Collectively, these projects illustrate AI's growing role in improving efficiency across JLab's nuclear physics mission.

Mei, Xinxin [Thomas Jefferson National Accelerator↗
DOE OSTI2025

Thermo-Mechanical Phase-Field Modeling of Fracture in High-Burnup UO2 Fuels Under Transient Conditions

This study presents a novel multiphysics phase-field fracture model to analyze high-burnup uranium dioxide (UO2) fuel behavior under transient reactor conditions. Fracture is treated as a stochastic phase transition, which inherently accounts for the random microstructural effects that lead to variations in the value of fracture strength. Moreover, the model takes into consideration the effects of temperature and burnup on thermal conductivity. Therefore, the model is able to predict crack initiation, propagation, and complex morphologies in response to thermal gradients and stress distributions. Several simulations were conducted to investigate the effects of operational and transient conditions on fracture behavior and the resulting cracking patterns. High-burnup fuels exhibit reduced thermal conductivity, elevating temperature gradients and resulting in extensive radial and circumferential cracks. Transient heating rates and temperatures significantly affect fracture patterns, with higher heating rates generating steeper gradients and more irregular crack trajectories. This approach provides critical insights into fuel integrity during accident scenarios and supports the safety evaluation of extended burnup limits.

Chemistry↗