Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Monte Carlo Simulations of Crystal Defects in Open Ensembles

Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. In conclusion, the method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.

Grain boundaries↗
DOE OSTI2022

Novel Proposals for FAIR, Automated, Recommendable, and Robust Workflows

Lightning talks of the Workflows in Support of Large-Scale Science (WORKS) workshop are a venue where the workflow community (researchers, developers, and users) can discuss work in progress, emerging technologies and frameworks, and training and education materials. This paper summarizes the WORKS 2022 lightning talks, which cover five broad topics: data integrity of scientific workflows; a machine learning-based recommendation system; a Python toolkit for running dynamic ensembles of simulations; a cross-platform, high-performance computing utility for processing shell commands; and a meta(data) framework for reproducing hybrid workflows.

Abhinit, Ishan↗
DOE OSTI2020

Statistical Treatment of Inverse Problems Constrained by Differential Equations-based Models with Stochastic Terms

This paper introduces a statistical treatment of inverse problems constrained by models with stochastic terms. The solution of the forward problem is given by a distribution represented numerically by an ensemble of simulations. The goal is to formulate the inverse problem, in particular the objective function, to find the closest forward distribution (i.e., the output of the stochastic forward problem) that best explains the distribution of the observations in a certain metric. We use proper scoring rules, a concept employed in statistical forecast verification, namely energy, variogram, and hybrid (i.e., combination of the two) scores. We study the performance of the proposed formulation in the context of two applications: a coefficient field inversion for subsurface flow governed by an elliptic partial differential equation with a stochastic source and a parameter inversion for power grid governed by differential-algebraic equations. In both cases we show that the variogram and the hybrid scores produce better parameter inversion results than does the energy score, whereas the energy score leads to better probabilistic predictions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Influence of the Indian Ocean Dipole on the large-scale circulation in South America

The influence of each phase of the Indian Ocean Dipole (IOD) on the largescale circulation in South America is investigated using rainfall observations, fully-coupled, large-ensemble, historical simulations (LENS), and forced experiments using the coupled model’s atmospheric component. IOD events often occur when El Niño Southern Oscillation (ENSO), the largest source of interannual variability of precipitation in South America, is active. To distinguish from effects of ENSO, only cases during neutral ENSO conditions are analyzed in LENS and observations. During the positive IOD polarity, a perturbation in the localWalker circulation leads to increased convection over equatorial South America, resulting in wet anomalies in the Amazon basin. This signal is the opposite of what is expected during El Niño events. Tropical convection anomalies in the Indian Ocean also force an extratropical Rossby wave train that reaches subtropical South America. During positive IOD, the moisture flux from the Amazon to central and southeastern Brazil weakens, resulting in a drying of the area associated with the South Atlantic Convergence Zone. Meanwhile, the South Atlantic Subtropical High strengthens, contributing to a drying in southeastern Brazil. During negative IOD, the induced wave train from the Indian Ocean leads to increased moisture transport to the La Plata basin, leading to wet anomalies in the region.

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

Truchas Overview

Truchas and Truchas-PBF are two sister codes for part-scale multi-physics modeling of manufacturing processes. Both programs are open source and made publicly available. They’re designed for efficient use of HPC resources and can be programmatically driven from Python packages. This enables automatic execution and analysis of ensembles of simulations, in some cases allowing 1000s of simulations to be evaluated in a day on HPC. Beyond just giving engineers a window into the concealed internal state of a system, the goal of Truchas is to provide a framework for developing novel manufacturing processes by understanding how the entire space of engineering inputs affects thermal state. It often is used to explore combinations of capabilities uncommon in commercial software, or to scale up analyses beyond the capabilities of commercial software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

AI/ML-Enhanced Wind Forecasts for Reducing Uncertainty in Prescribed Fire Planning

Prescribed fire is a vital tool for ecosystem management and wildfire risk reduction but its escalation is constrained by overly conservative burn windows because of uncertainties, for instance, in wind forecasts. This review describes the state of the art in weather product use by fire/smoke models and identifies three priority research gaps that artificial intelligence/machine learning (AI/ML) is well positioned to address: (1) spatial and temporal downscaling to meter-scale, sub-hourly wind fields; (2) bias correction for systematic model errors in complex terrain; and (3) robust uncertainty quantification to inform ensemble-based simulations. Emerging AI/ML techniques offer promising frameworks to address all three challenges. By providing high-resolution, bias-corrected, and probabilistic wind fields, AI/ML-enhanced forecasts will allow for expanded burn windows, improved ignition strategy design and a reduced reliance on expert intuition, especially when a prescribed fire is introduced into new areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Vibrational Entropy of Crystalline Solids from Covariance of Atomic Displacements

The vibrational entropy of a solid at finite temperature is investigated from the perspective of information theory. Ab initio molecular dynamics (AIMD) simulations generate ensembles of atomic configurations at finite temperature from which we obtain the N-body distribution of atomic displacements, ρN. We calculate the information-theoretic entropy from the expectation value of lnρN. At a first level of approximation, treating individual atomic displacements independently, our method may be applied using Debye–Waller B-factors, allowing diffraction experiments to obtain an upper bound on the thermodynamic entropy. At the next level of approximation we correct the overestimation through inclusion of displacement covariances. We apply this approach to elemental body-centered cubic sodium and face-centered cubic aluminum, showing good agreement with experimental values above the Debye temperatures of the metals. Below the Debye temperatures, we extract an effective vibrational density of states from eigenvalues of the covariance matrix, and then evaluate the entropy quantum mechanically, again yielding good agreement with experiment down to low temperatures. Our method readily generalizes to complex solids, as we demonstrate for a high entropy alloy. Further, our method applies in cases where the quasiharmonic approximation fails, as we demonstrate by calculating the HCP/BCC transition in Ti.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Indirect climate impacts of the Hunga eruption

Injection of sulfur and water vapour by the Hunga volcanic eruption significantly altered chemical composition and radiative budget of the stratosphere. Yet, whether the eruption could also affect surface climate, especially via indirect pathways, remains poorly understood. Here we investigate these effects using large ensembles of simulations with the CESM2(WACCM6) Earth system model, incorporating interactive chemistry and aerosols in both coupled ocean and atmosphere-only configurations. We find some statistically significant extratropical regional climate responses to the eruption driven by circulation changes; these are partially linked to the modulation of El Niño–Southern Oscillation, and its associated teleconnections, and to perturbations of the stratospheric polar vortex in both hemispheres. The stratospheric anomalies affect surface climate through modulating the North Atlantic Oscillation in the Northern Hemisphere (up to three boreal winters following the eruption) and the Southern Annular Mode in the Southern Hemisphere in late 2023. The latter is partly related to a concurrent reduction in Antarctic ozone, as increased stratospheric aerosols and water vapor reach the polar vortex. Our study suggests that the eruption could have had a non-negligible influence on regional surface climate, and discusses the mechanisms via which such an influence could occur. However, the results also highlight that this forcing is relatively weak compared to interannual variability, and is subject to model uncertainties in the representation of key processes. More research is thus needed before definitive statements on the role of the eruption in contributing to surface climate and weather events in the following years are made.

Geosciences↗
DOE OSTI2021

Quantifying the role of moss in terrestrial ecosystem carbon dynamics in northern high latitudes

Mosses are ubiquitous in northern terrestrial ecosystems, and play an important role in regional carbon, water and energy cycling. Current global land surface models that do not consider mosses may bias the quantification of regional carbon dynamics. Here we incorporate mosses as a new plant functional type into the process-based Terrestrial Ecosystem Model (TEM 5.0), to develop a new model (TEM_Moss). The new model explicitly quantifies the interactions between vascular plants and mosses and their competition for energy, water, and nutrients. Compared to the estimates using TEM 5.0, the new model estimates that the regional terrestrial soils currently store 132.7 Pg more C and will store 157.5 and 179.1 Pg more C under the RCP8.5 and RCP2.6 scenarios, respectively, by the end of the 21st century. Ensemble regional simulations forced with different parameters for the 21st century with TEM_Moss predict that the region will accumulate 161.1±142.1 Pg C under the RCP2.6 scenario and 186.7±166.1 Pg C under the RCP8.5 scenario over the century. Our study highlights the necessity of coupling moss into Earth system models to adequately quantify terrestrial carbon–climate feedbacks in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Indices of extremes: geographic patterns of change in extremes and associated vegetation impacts under climate intervention

Extreme weather events have been demonstrated to be increasing in frequency and intensity across the globe and are anticipated to increase further with projected changes in climate. Solar climate intervention strategies, specifically stratospheric aerosol injection (SAI), have the potential to minimize some of the impacts of a changing climate while more robust reductions in greenhouse gas emissions take effect. However, to date little attention has been paid to the possible responses of extreme weather and climate events under climate intervention scenarios. We present an analysis of 16 extreme surface temperature and precipitation indices, as well as associated vegetation responses, applied to the Geoengineering Large Ensemble (GLENS). GLENS is an ensemble of simulations performed with the Community Earth System Model (CESM1) wherein SAI is simulated to offset the warming produced by a high-emission scenario throughout the 21st century, maintaining surface temperatures at 2020 levels. GLENS is generally successful at maintaining global mean temperature near 2020 levels; however, it does not completely offset some of the projected warming in northern latitudes. Some regions are also projected to cool substantially in comparison to the present day, with the greatest decreases in daytime temperatures. The differential warming–cooling also translates to fewer very hot days but more very hot nights during the summer and fewer very cold days or nights compared to the current day. Extreme precipitation patterns, for the most part, are projected to reduce in intensity in areas that are wet in the current climate and increase in intensity in dry areas. We also find that the distribution of daily precipitation becomes more consistent with more days with light rain and fewer very intense events than currently occur. In many regions there is a reduction in the persistence of long dry and wet spells compared to present day. However, asymmetry in the night and day temperatures, together with changes in cloud cover and vegetative responses, could exacerbate drying in regions that are already sensitive to drought. Overall, our results suggest that while SAI may ameliorate some of the extreme weather hazards produced by global warming, it would also present some significant differences in the distribution of climate extremes compared to the present day.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An ensemble of 48 physically perturbed model estimates of the 1/8° terrestrial water budget over the conterminous United States, 1980–2015

Terrestrial water budget (TWB) data over large domains are of high interest for various hydrological applications. Spatiotemporally continuous and physically consistent estimations of TWB rely on land surface models (LSMs). As an augmentation of the operational North American Land Data Assimilation System Phase 2 (NLDAS-2) four-LSM ensemble, this paper describes a dataset simulated from an ensemble of 48 physics configurations of the Noah LSM with multi-physics options (Noah-MP). The 48 Noah-MP physics configurations are selected to give a representative cross-section of commonly used LSMs for parameterizing runoff, atmospheric surface layer turbulence, soil moisture limitation on photosynthesis, and stomatal conductance. The dataset spans from 1980 to 2015 over the conterminous United States (CONUS) at a monthly temporal resolution and a 1/8° spatial resolution. The dataset variables include total evapotranspiration and its constituents (canopy evaporation, soil evaporation, and transpiration), runoff (the surface and subsurface components), as well as terrestrial water storage (snow water equivalent, four-layer soil water content from the surface down to 2 m, and the groundwater storage anomaly). The dataset is available at https://doi.org/10.5281/zenodo.7109816. Evaluations carried out in this study and previous investigations show that the ensemble per forms well in reproducing the observed terrestrial water storage, snow water equivalent, soil moisture, and runoff. Noah-MP complements the NLDAS models well, and adding Noah-MP consistently improves the NLDAS es timations of the above variables in most areas of CONUS. Besides, the perturbed-physics ensemble facilitates the identification of model deficiencies. The parameterizations of shallow snow, spatially varying groundwater dynamics, and near-surface atmospheric turbulence should be improved in future model versions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Downscaling atmospheric chemistry simulations with physically consistent deep learning

Abstract. Recent advances in deep convolutional neural network (CNN)-based super resolution can be used to downscale atmospheric chemistry simulations with substantially higher accuracy than conventional downscaling methods. This work both demonstrates the downscaling capabilities of modern CNN-based single image super resolution and video super-resolution schemes and develops modifications to these schemes to ensure they are appropriate for use with physical science data. The CNN-based video super-resolution schemes in particular incur only 39 % to 54 % of the grid-cell-level error of interpolation schemes and generate outputs with extremely realistic small-scale variability based on multiple perceptual quality metrics while performing a large (8×10) increase in resolution in the spatial dimensions. Methods are introduced to strictly enforce physical conservation laws within CNNs, perform large and asymmetric resolution changes between common model grid resolutions, account for non-uniform grid-cell areas, super-resolve lognormally distributed datasets, and leverage additional inputs such as high-resolution climatologies and model state variables. High-resolution chemistry simulations are critical for modeling regional air quality and for understanding future climate, and CNN-based downscaling has the potential to generate these high-resolution simulations and ensembles at a fraction of the computational cost.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Extreme windstorms in the Northeastern USA in the contemporary and future climate

Cold-season windstorms represent an important, and potentially changing, geophysical hazard in the Northeastern United States. Here we employ an integrated research methodology including both a storyline approach, where three intense windstorms from the current climate are subjected to pseudo-global warming (PGW) experiments, and a long-term transient simulation using the Weather Research and Forecasting (WRF) model. An ensemble of WRF simulations is built for each windstorm using different planetary boundary layer and microphysical parameterizations. The fidelity assessment suggests all ensemble members capture the windstorm evolution in contemporary climate. The configuration with highest fidelity is used in the PGW experiments performed with perturbed temperature fields, constant relative humidity, and deiced Great Lakes. These perturbation simulations indicate some evidence for a reduction of sea level pressure and increases in wind speed over and downwind of the Great Lakes and over the Atlantic Ocean plus an increase in precipitation accumulation but a reduction in snow coverage. These changes are spatially inhomogeneous in terms of magnitude and sign but are consistent with changes in potential vorticity. Alberta Clippers and Colorado Lows dominate the cyclones responsible for historical windstorms and thus are sampled in the PGW simulations. However, the transient simulation suggests an increasing role for tropical cyclones that undergo transition to extratropical cyclones. Furthermore, this reinforces the value of combining information from both PGW perturbation experiments within a storyline context and transient simulations when seeking to quantify the future risk associated with cold-season windstorms under changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Minimizing CGYRO HPC Communication Costs in Ensembles with XGYRO by Sharing the Collisional Constant Tensor Structure

First-principles fusion plasma simulations are both compute and memory intensive, and CGYRO is no exception. The use of many HPC nodes to fit the problem in the available memory thus results in significant communication overhead, which is hard to avoid for any single simulation. That said, most fusion studies are composed of ensembles of simulations, so we developed a new tool, named XGYRO, that executes a whole ensemble of CGYRO simulations as a single HPC job. By treating the ensemble as a unit, XGYRO can alter the global buffer distribution logic and apply optimizations that are not feasible on any single simulation, but only on the ensemble as a whole. The main saving comes from the sharing of the collisional constant tensor structure, since its values are typically identical between parameter-sweep simulations. This data structure dominates the memory consumption of CGYRO simulations, so distributing it among the whole ensemble results in drastic memory savings for each simulation, which in turn results in overall lower communication overhead.

CGYRO↗
DOE OSTI2021

Evaluation and Uncertainty of Radioxenon Transport with a Mesoscale Model after the February 2013 Underground Test in North Korea

Abstract The transport of radioxenon released from the February 2013 underground nuclear weapons test in North Korea was analyzed at two receptors—one at the Comprehensive Nuclear Test Ban Treaty site Rn58 in Russia (400 km downwind) and a second at Rn38 in Japan (1000 km downwind). Transport was modeled with two ensembles of mesoscale simulations, one generated with varying initial and lateral boundary conditions taken from the Global Forecasting System uncertainty ensemble, and a second created from different parameterizations and surface conditions. The wind variability was similar for the two ensembles and consistent with observations at 925 mb (1 mb = 1 hPa) but not at the surface. Biases in calculated surface winds and the radioxenon concentration in the ensembles were attributed mainly to poor simulation of the sea breeze at both locations and mountain lee affects at Rn38 in Japan. These wind regimes affected the timing of the surface radioxenon plume at Rn58 and its duration at Rn38. Surface wind variability induced by terrain and land–sea contrast (the sea breeze) also had a significant effect on the surface winds and plume dynamics, including blocking of flow approaching elevated terrain near Vladivostok and the west side of Japan. Increased plume uncertainty was seen at night because of surface wind variability. Measured surface chemical variability was larger than found in the first European Tracer Experiment in central Europe. The study found that horizontal model resolution contributes to uncertainty but not as much as vertical resolution, boundary layer parameterizations, and assimilation of surface meteorological data near the receptor.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

gRASPA

GPU Monte Carlo Simulation Code with a taste of RASPA We present enhancements in Monte Carlo simulation speed and functionality within an open-source code, gRASPA, which uses graphical processing units (GPUs) to achieve significant performance improvements compared to serial, CPU implementations of Monte Carlo. The code supports a wide range of Monte Carlo simulations, including canonical ensemble (NVT), grand canonical, NVT Gibbs, Widom test particle insertions, and continuous-fractional component Monte Carlo. Implementation of grand canonical transition matrix Monte Carlo (GC-TMMC) and a novel feature to allow different moves for the different components of metal-organic framework (MOF) structures exemplify the capabilities of gRASPA for precise free energy calculations and enhanced adsorption studies, respectively. The introduction of a High-Throughput Computing (HTC) mode permits many Monte Carlo simulations on a single GPU device for accelerated materials discovery. The code can incorporate machine learning (ML) potentials. The open-source nature of gRASPA promotes reproducibility and openness in science, and users may add features to the code and optimize it for their own purposes. The code is written in CUDA/C++ and SYCL/C++ to support different GPU vendors. The gRASPA code is publicly available at https://github.com/snurr-group/gRASPA.

Li, Zhao [Purdue/Northwestern/Notre Dame Universit↗
DOE OSTI2024

Changing effects of external forcing on Atlantic–Pacific interactions

Recent studies have highlighted the increasingly dominant role of external forcing in driving Atlantic and Pacific Ocean variability during the second half of the 20th century. This paper provides insights into the underlying mechanisms driving interactions between modes of variability over the two basins. We define a set of possible drivers of these interactions and apply causal discovery to reanalysis data, two ensembles of pacemaker simulations where sea surface temperatures in either the tropical Pacific or the North Atlantic are nudged to observations, and a pre-industrial control run. We also utilize large-ensemble means of historical simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6) to quantify the effect of external forcing and improve the understanding of its impact. A causal analysis of the historical time series between 1950 and 2014 identifies a regime switch in the interactions between major modes of Atlantic and Pacific climate variability in both reanalysis and pacemaker simulations. A sliding window causal analysis reveals a decaying El Niño–Southern Oscillation (ENSO) effect on the Atlantic as the North Atlantic fluctuates towards an anomalously warm state. The causal networks also demonstrate that external forcing contributed to strengthening the Atlantic's negative-sign effect on ENSO since the mid-1980s, where warming tropical Atlantic sea surface temperatures induce a La Niña-like cooling in the equatorial Pacific during the following season through an intensification of the Pacific Walker circulation. The strengthening of this effect is not detected when the historical external forcing signal is removed in the Pacific pacemaker ensemble. The analysis of the pre-industrial control run supports the notion that the Atlantic and Pacific modes of natural climate variability exert contrasting impacts on each other even in the absence of anthropogenic forcing. The interactions are shown to be modulated by the (multi)decadal states of temperature anomalies of both basins with stronger connections when these states are “out of phase”. We show that causal discovery can detect previously documented connections and provides important potential for a deeper understanding of the mechanisms driving changes in regional and global climate variability.

54 ENVIRONMENTAL SCIENCES↗