Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system 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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dissecting Anvil Cloud Response to Sea Surface Warming

Abstract We derive an anvil cloud diagnostic from the continuity equation of cloud ice and apply it to the output of convection‐permitting Energy Exascale Earth System Model (E3SM) simulations run in radiative‐convective equilibrium mode. This diagnostic shows that anvil cloud fraction can be reliably diagnosed as a product of cloud detrainment and lifetime. Detrainment is found to be approximated well by a product of clear sky convergence and cloud ice mixing ratio, while cloud lifetime is dominated by sedimentation. Taken together, this diagnostic expresses anvil cloud fraction as a function of five physically measurable quantities. Of these, clear‐sky convergence changes drive the anvil cloud reduction with warming while an increase in cloud ice mixing ratio buffers the decrease. Accordingly, this study provides a theoretical foundation upon which the Stability‐Iris hypothesis can be tested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Model Evaluation and Intercomparison of Marine Warm Low Cloud Fractions With Neural Network Ensembles

Abstract Low cloud fractions (LCFs) and meteorological factors (MFs) over an oceanic region containing multiple cloud regimes are examined for three data sets: one Energy Exascale Earth System Model (E3SM) simulation with the default 72‐layer vertical grid (E3SM72), another one with 8‐times vertical resolution via the Framework for Improvement by Vertical Enhancement (E3SM 8), and one with MFs from ERA5 reanalysis and LCFs from the CERES SSF product (ERA5‐SSF). Neural networks (NNs) are trained to capture the relationship between MFs and LCF and to select the best‐performing MF subsets for predicting LCF. NN ensembles are used to (a) confirm the performance of selected MF subsets, (b) to serve as proxy models for each data set to predict LCFs for MFs from all data sets, and (c) to classify MFs into those in shared and uniquely occupied MF subspaces. Overall, E3SM72 and E3SM 8 have large fractions of MFs in shared MF subspace, but less so near the Californian and Peruvian stratocumulus decks. E3SM 8 and ERA5 have small fractions of MFs in shared MF subspace but greater than E3SM72 and ERA5, especially in the Southeast Pacific. The differences in LCFs between three pairs of data sets are decomposed into those associated with the differences in the LCF‐MF relationship and those involving different MFs. Given the same MFs, LCFs produced by E3SM 8 are greater than those produced by E3SM72 but are still different from those in ERA5‐SSF. In general, the shift in MFs dominates the difference in the LCFs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Maximum Supersaturation in the Marine Boundary Layer Clouds Over the North Atlantic

Abstract The maximum supersaturation ( S x ) in clouds is a key parameter affecting the cloud's microphysical and radiative properties. We investigate the S x of the marine boundary layer clouds by combining airborne and surface observations in the Eastern North Atlantic. The cloud droplet number concentration ( N c ) in the least diluted cloud cores agrees well with the number concentration of particles larger than the Hoppel Minimum (HM) ( N >HM ) below clouds, indicating that the HM represents the average size threshold above which particles are activated to form cloud droplets. The S x values derived from surface observations vary from 0.10% to 0.50% from June 2017 to June 2018, with a clear seasonal variation exhibiting higher values during winter. Most of the S x variance (∼60%) can be explained by the cloud condensation nuclei (CCN) concentration and updraft velocity ( w ), with the CCN concentration playing a more important role than w in explaining the variation of S x . The influence of CCN concentration on S x leads to a buffered response of N c to aerosol perturbations. The response of N c to low aerosol concentration during winter is further buffered by the high w . The global Community Earth System Model (CESM) simulated S x values in the Azores have a positive bias compared to measured S x , likely due to overestimated w and underestimated CCN concentration. The CESM simulated S x exhibits higher values further north over the North Atlantic Ocean, which is attributed to stronger w . The suppression of S x by aerosol is also evident in regions with high CCN concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

High‐Latitude Stratospheric Aerosol Injection to Preserve the Arctic

Abstract Stratospheric aerosol injection (SAI) has been shown in climate models to reduce some impacts of global warming in the Arctic, including the loss of sea ice, permafrost thaw, and reduction of Greenland Ice Sheet (GrIS) mass; SAI at high latitudes could preferentially target these impacts. In this study, we use the Community Earth System Model to simulate two Arctic‐focused SAI strategies, which inject at 60°N latitude each spring with injection rates adjusted to either maintain September Arctic sea ice at 2030 levels (“Arctic Low”) or restore it to 2010 levels (“Arctic High”). Both simulations maintain or restore September sea ice to within 10% of their respective targets, reduce permafrost thaw, and increase GrIS surface mass balance by reducing runoff. Arctic High reduces these impacts more effectively than a globally focused SAI strategy that injects similar quantities of SO 2 at lower latitudes. However, Arctic‐focused SAI is not merely a “reset button” for the Arctic climate, but brings about a novel climate state, including changes to the seasonal cycles of Northern Hemisphere temperature and sea ice and less high‐latitude carbon uptake relative to SSP2‐4.5. Additionally, while Arctic‐focused SAI produces the most cooling near the pole, its effects are not confined to the Arctic, including detectable cooling throughout most of the northern hemisphere for both simulations, increased mid‐latitude sulfur deposition, and a southward shift of the location of the Intertropical Convergence Zone. For these reasons, it would be incorrect to consider Arctic‐focused SAI as “local” geoengineering, even when compared to a globally focused strategy.

Lee, Walker Raymond↗
DOE OSTI2022

Machine Learning–Based Detection of Weather Fronts and Associated Extreme Precipitation in Historical and Future Climates

Extreme precipitation events, including those associated with weather fronts, have wide-ranging impacts across the world. Here we use a deep learning algorithm to identify weather fronts in high resolution Community Earth System Model (CESM) simulations over the contiguous United States (CONUS), and evaluate the results using observational and reanalysis products. We further compare results between CESM simulations using present-day and future climate forcing, to study how these features might change with climate change. We find that detected front frequencies in CESM have seasonally varying spatial patterns and responses to climate change and are found to be associated with modeled changes in large scale circulation such as the jet stream. We also associate the detected fronts with precipitation and find that total and extreme frontal precipitation mostly decreases with climate change, with some seasonal and regional differences. Here, decreases in Northern Hemisphere summer frontal precipitation are largely driven by changes in the frequency of different front types, especially cold and stationary fronts. On the other hand, Northern Hemisphere winter exhibits some regional increases in frontal precipitation that are largely driven by changes in frontal precipitation intensity. While CONUS mean and extreme precipitation generally increase during all seasons in these climate change simulations, the likelihood of frontal extreme precipitation decreases, demonstrating that extreme precipitation has seasonally varying sources and mechanisms that will continue to evolve with climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Different Strategies of Stratospheric Aerosol Injection Would Significantly Affect Climate Extreme Mitigation

Abstract Stratospheric aerosol injection (SAI) has been proposed as a potential supplement to mitigate some climate impacts of anthropogenic warming. Using Community Earth System Model ensemble simulation results, we analyze the response of temperature and precipitation extremes to two different SAI strategies: one injects SO 2 at the equator to stabilize global mean temperature and the other injects SO 2 at multiple locations to stabilize global mean temperature as well as the interhemispheric and equator‐to‐pole temperature gradients. Our analysis shows that in the late 21st century, compared with the present‐day climate, both equatorial and multi‐location injection lead to reduced hot extremes in the tropics, corresponding to overcooling of the mean climate state. In mid‐to‐high latitude regions, in comparison to the present‐day climate, substantial decreases in cold extremes are observed under both equatorial and multi‐location injection, corresponding to residual winter warming of the mean climate state. Both equatorial and multi‐location injection reduce precipitation extremes in the tropics below the present‐day level, associated with the decrease in mean precipitation. Overall, for most regions, temperature and precipitation extremes show reduced change in response to multi‐location injection than to equatorial injection, corresponding to reduced mean climate change for multi‐location injection. In comparison with equatorial injection, in response to multi‐location injection, most land regions experience fewer years with significant change in cold extremes from the present‐day level, and most tropical regions experience fewer years with significant change in hot extremes. The design of SAI strategies to mitigate anthropogenic climate extremes merits further study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Climate Conundrum: A Wet or Dry European and Northern African Climate During the Middle Miocene

Abstract End of 21st‐century hydroclimate projections suggest an expansion of subtropical dry zones, with Mediterranean and Sahel regions becoming much drier. However, paleobotanical assemblage evidence from the middle Miocene (17‐12 Ma), suggests both regions were instead humid environments. Here we show that by modifying regional sea surface temperatures (SST) in an Earth System Model (CESM1.2) simulation of the middle Miocene, the increased ocean evaporation and integrated water vapor flux overrides any drying effects associated with warming‐induced land‐surface evaporation driven by atmospheric CO 2 concentrations. These modifications markedly reduce the bias in the model‐data comparison for this period. A vegetation model (BIOME4) forced with simulated climatologies predicts both regions were dominated by mixed forest, which is largely consistent with the paleobotanical record. This study unveils the potential for wetter subtropical Mediterranean climates associated with warming, presenting an alternative scenario from future drying projections with localized SST warming governing regional climate change.

Acosta, R. P.↗
DOE OSTI2025

Environmental Influences on Deep Convective Upscale Growth Rate in Central Argentina From a Convection‐Permitting Simulation

This study uses a convection‐permitting model simulation to describe the environmental conditions under which convective upscale growth occurs in central Argentina, particularly examining environmental parameters when deep convection initially forms that could differentiate the rate of initial upscale growth. Simulated mesoscale convective systems (MCSs) are separated into slow and rapid growth by the rate of spatial growth from convection initiation until reaching the MCS scale. A low‐level jet (LLJ) is found more frequently near the deep convection that experiences rapid growth to an MCS, but its presence alone is not predictive of rapid growth. Using spatially‐averaged parameters, we find that rapid growth to MCSs also occurs in environments that are significantly more thermodynamically favorable with greater low‐level moisture and instability. Fewer significant differences are found in the kinematic environment with only the 0–2 km vertical wind shear magnitude being significantly larger for rapid growth MCSs compared to slow growth MCSs, potentially related to LLJs often peaking near this height. When focusing only on MCSs with the slowest and fastest growth rates, elevated‐layer shear is significantly smaller for very rapid growth MCSs, suggesting elevated‐layer shear may help discriminate between the upper and lower bounds of growth rate. Finally, when upscale growth occurs near the Sierras de Córdoba (SDC) with a LLJ present, rapid growth is also supported by favorable wind shear orientation. However, this does not hold for upscale growth occurring away from the SDC, highlighting the importance of interpreting shear direction relative to the orientation of features initiating deep convection.

Sasaki, Clayton R. S. [Univ. of Washington, Seattl↗
DOE OSTI2026

Vegetation responses to air dryness amplify future land surface warming

Temperature exerts a first-order control on vegetation photosynthesis and transpiration. Yet most studies investigating temperature impacts on plants rely on near-surface air temperature, rather than canopy temperature—the temperature plants actually experience. Because canopy temperature directly regulates ecosystem function, it provides a more accurate measure of vegetation–climate interactions. Combining Earth System Model (ESM) simulations and satellite observations in a dual emergent constraint, here we show that canopy temperature is projected to increase substantially more than air temperature (~0.11-degrees more or a 16% increase in their difference) over the 21st century. The ESM ensemble median fails to capture these stronger increases in the majority of vegetated regions. We find that the largest projected increases in the difference between canopy and air temperature are predicted to occur in regions where elevating moisture stress—particularly rising vapor pressure deficit—increasingly constrains vegetation growth and transpiration. This implies that future warming will impose stronger constraints on plant function than currently estimated. Relying on air temperature alone will therefore lead to systematic underestimation of temperature effects on photosynthesis, vegetation growth, and the land carbon sink. Accurate representation of canopy temperature in ESMs is thus essential to improve projections of ecosystem responses and feedbacks to climate change.

Green, Julia K↗
DOE OSTI2022

Sampling lattices in semi-grand canonical ensemble with autoregressive machine learning

Calculating thermodynamic potentials and observables efficiently and accurately is key for the application of statistical mechanics simulations to materials science. However, naive Monte Carlo approaches, on which such calculations are often dependent, struggle to scale to complex materials in many state-of-the-art disciplines such as the design of high entropy alloys or multi-component catalysts. To address this issue, we adapt sampling tools built upon machine learning-based generative modeling to the materials space by transforming them into the semi-grand canonical ensemble. Furthermore, we show that the resulting models are transferable across wide ranges of thermodynamic conditions and can be implemented with any internal energy model U, allowing integration into many existing materials workflows. We demonstrate the applicability of this approach to the simulation of benchmark systems (AgPd, CuAu) that exhibit diverse thermodynamic behavior in their phase diagrams. Finally, we discuss remaining challenges in model development and promising research directions for future improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Continuous recoil-driven lasing and cavity frequency pinning with laser-cooled atoms

Abstract Laser-cooled gases of atoms interacting with the field of an optical cavity are a versatile tool for quantum sensing and the simulation of quantum systems. These systems can exhibit phenomena such as self-organization phase transitions, lasing mechanisms, squeezed states and protection of quantum coherence. However, investigations of these phenomena typically occur in a discontinuous manner due to the need to reload atomic ensembles. Here we demonstrate hours-long continuous lasing from laser-cooled 88 Sr atoms loaded into a ring cavity. The required inversion to produce lasing arises from inversion in the atomic-momentum degrees of freedom, which is linked to the self-organization phase transitions and collective atomic recoil lasing observed previously only in a cyclic fashion. We find that over a broad parameter range, the sensitivity of the lasing frequency to changes in cavity frequency is significantly reduced due to an atomic loss mechanism, suggesting a potential approach for mitigating low-frequency cavity noise. Our findings open opportunities for continuous cavity quantum electrodynamics experiments and robust and continuous super-radiant lasers.

Schäfer, Vera M. (ORCID:0000000210003737)↗
DOE OSTI2021

AI-based design of a nuclear reactor core

The authors developed an artificial intelligence (AI)-based algorithm for the design and optimization of a nuclear reactor core based on a flexible geometry and demonstrated a 3× improvement in the selected performance metric: temperature peaking factor. The rapid development of advanced, and specifically, additive manufacturing (3-D printing) and its introduction into advanced nuclear core design through the Transformational Challenge Reactor program have presented the opportunity to explore the arbitrary geometry design of nuclear-heated structures. The primary challenge is that the arbitrary geometry design space is vast and requires the computational evaluation of many candidate designs, and the multiphysics simulation of nuclear systems is very time-intensive. Therefore, the authors developed a machine learning-based multiphysics emulator and evaluated thousands of candidate geometries on Summit, Oak Ridge National Laboratory’s leadership class supercomputer. The results presented in this work demonstrate temperature distribution smoothing in a nuclear reactor core through the manipulation of the geometry, which is traditionally achieved in light water reactors through variable assembly loading in the axial direction and fuel shuffling during refueling in the radial direction. The conclusions discuss the future implications for nuclear systems design with arbitrary geometry and the potential for AI-based autonomous design algorithms.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Relaxation of stationary states on a quantum computer yields a unique spectroscopic fingerprint of the computer’s noise

Abstract Quantum computing has the potential to revolutionize computing, but its significant sensitivity to noise requires sophisticated error correction and mitigation. Traditionally, noise on the quantum device is characterized directly through qubit and gate measurements, but this approach has drawbacks in that it does not adequately capture the effect of noise on realistic multi-qubit applications. In this paper, we simulate the relaxation of stationary quantum states on a quantum computer to obtain a unique spectroscopic fingerprint of the computer’s noise. In contrast to traditional approaches, we obtain the frequency profile of the noise as it is experienced by the simulated stationary quantum states. Data from multiple superconducting-qubit IBM processors show that noise generates a bath within the simulation that exhibits both colored noise and non-Markovian behavior. Our results provide a direction for noise mitigation but also suggest how to use noise for quantum simulations of open systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Unraveling membrane electrode assembly design for electrochemical conversion of carbon dioxide to formate/formic acid

This work presents a one-dimensional continuum modeling approach to investigate various cell architectures used for electrochemical conversion of CO 2 to formate/formic acid. Ion transport is simulated by a system of generalized modified Poisson–Nernst–Planck (GMPNP) equations that reflect the reactive transport phenomena including steric effects as the electrolyte solutions become concentrated. In the cathode catalyst layer, ionic current contributions from both the supporting electrolyte and solid-state ionomer are considered. Voltage and CO 2 utilization breakdowns are utilized to deconvolute the impacts of the cell architecture. The origins of (bi)carbonate formation in the cathode are explored, as the subsequent decrease in CO 2 availability is a key reason for low faradaic efficiencies to formate/formic acid. In addition, the role of a supporting electrolyte (KOH) is investigated to understand its tradeoffs: while the K + ions can improve both conductivity and electrochemically active surface area in the cathode, the presence of OH − ions raises the pH and leads to deleterious formation of (bi)carbonates. To this end, we also present parametric studies on the concentration and flow rate of supplied KOH to the cell, to establish a path towards eliminating the need for a supporting electrolyte.

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

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counterpropagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scattering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. Finally, we compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Writable spin wave nanochannels in an artificial-spin-ice-mediated ferromagnetic thin film

We report magnonics, which employs spin-waves to transmit and process information, is a promising venue for low-power data processing. One of the major challenges is the local control of the spin-wave propagation path. Here, we introduce the concept of writable magnonics by taking advantage of the highly flexible reconfigurability and rewritability of artificial spin ice systems. Using micromagnetic simulations, we show that globally switchable spin-wave propagation and locally writable spin-wave nanochannels can be realized in a ferromagnetic thin film underlying an artificial pinwheel spin ice. The rewritable magnonics enabled by reconfigurable spin wave nanochannels provides a unique setting to design programmable magnonic circuits and logic devices for ultra-low power applications.

42 ENGINEERING↗
DOE OSTI2023

A parallel, distributed memory implementation of the adaptive sampling configuration interaction method

The many-body simulation of quantum systems is an active field of research that involves several different methods targeting various computing platforms. Many methods commonly employed, particularly coupled cluster methods, have been adapted to leverage the latest advances in modern high-performance computing. Selected configuration interaction (sCI) methods have seen extensive usage and development in recent years. However, the development of sCI methods targeting massively parallel resources has been explored only in a few research works. Here, we present a parallel, distributed memory implementation of the adaptive sampling configuration interaction approach (ASCI) for sCI. In particular, we will address the key concerns pertaining to the parallelization of the determinant search and selection, Hamiltonian formation, and the variational eigenvalue calculation for the ASCI method. Load balancing in the search step is achieved through the application of memory-efficient determinant constraints originally developed for the ASCI-PT2 method. The presented benchmarks demonstrate near optimal speedup for ASCI calculations of Cr 2 (24e, 30o) with 10 6 , 10 7 , and 3 × 10 8 variational determinants on up to 16 384 CPUs. Importantly, to the best of the authors’ knowledge, this is the largest variational ASCI calculation to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗