Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large 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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Flower‐Type Organized Trade‐Wind Cumulus: A Multi‐Day Lagrangian Large Eddy Simulation Intercomparison Study

Shallow cumulus cloud fields in subtropical marine trade wind environments, particularly over the tropical Atlantic Ocean, show distinct organizational patterns. Among these, Flower‐type clouds are characterized by expansive stratiform cloud patches surrounded by regions of scattered convection. The objectives of this study were (a) to construct a case study of a time period during the EUREC 4 A/ATOMIC field campaign when Flower‐type organization was observed, (b) to evaluate the fidelity of a multi‐model ensemble of large eddy simulations of that case, and (c) to analyze the interaction between cloud and precipitation processes and mesoscale organization in the simulations. The simulations follow a quasi‐Lagrangian trajectory, allowing mesoscale features to develop over time in a domain that follows the boundary‐layer airmass. The results show a broad agreement in simulated thermodynamic properties across different LES codes, with Flower‐type cloud patches appearing within hours of each other. The consensus among models is consistent with observations made during the EUREC 4 A/ATOMIC field campaign on the specific day of interest. The cloud structure reveals three distinct peaks in the joint probability densities of cloud base and cloud top height, with the dominant peak at any given time influenced by the stage of cloud organization. The simulated cloud system evolution reveals consistent occurrence of maxima in liquid water path and rain rate before Flower reaches its maximum length scale. Targeted sensitivity tests reveal a weak relationship between Cloud Droplet Number concentration and the extent/degree/type of organization.

EUREC4A↗
DOE OSTI2026

Quantum-Inspired Bayesian Sampling for Uncertainty Quantification and Machine Learning (Final Technical Report)

With increasing simulation and measurement data, machine learning and artificial intelligence have been widely used in computational decision-making of complex engineering systems. The resulting tools, such as uncertainty quantification solvers, reinforcement learning, and physics-informed machine learning, have achieved great success in critical DOE tasks such as material discovery and design, energy system modeling and control, and numerical weather and climate prediction. A core topic in scientific machine learning and artificial intelligence is Bayesian inference: given an observed data set, people want to estimate the posterior distribution of a (possibly large) number of hidden parameters. Due to the flexibility and weak assumptions, Bayesian sampling has been the mainstream Bayesian inference solvers despite the rapid progress of approximate Bayesian inference. Classical Bayesian sampling methods such as Markov-chain Monte Carlo suffer from a low-acceptance rate due to the random walk nature, therefore state-of-the-art techniques use Hamiltonian Monte Carlo and its variants to efficiently draw posterior samples in a high dimension. The key idea of Hamiltonian Monte Carlo and its variants is to simulate the Hamiltonian dynamics of a classical particle with a fixed mass, and their performance significantly degrades when the posterior distribution is highly spiky or has multiple modes. Leveraging the idea of quantum physics, this project has investigated new theory, algorithms and applications of Bayesian inference (especially Bayesian sampling). The main results include: (1) novel quantum-inspired Bayesian sampling methods that can lead to better accuracy for challenging multi-modal or spiky distributions, (2) more scalable machine learning framework leveraging tensor-compressed Bayesian inference, and (3) Bayesian and sampling approaches for verifying the robustness of continuous and binary neural networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Interfacial Atomistic Mechanisms of Lithium Metal Stripping and Plating in Solid–State Batteries

All-solid-state batteries based on a Li metal anode represent a promising next-generation energy storage system, but are currently limited by low current density and short cycle life. Further research to improve the Li metal anode is impeded by the lack of understanding in its failure mechanisms at lithium–solid interfaces, in particular, the fundamental atomistic processes responsible for interface failure. Here, using large-scale molecular dynamics simulations, the first atomistic modeling study of lithium stripping and plating on a solid electrolyte is performed by explicitly considering key fundamental atomistic processes and interface atomistic structures. In the simulations, the interface failure initiated with the formation of nano-sized pores, and how interface structures, lithium diffusion, adhesion energy, and applied pressure affect interface failure during Li cycling are observed. By systematically varying the parameters of solid-state lithium cells in the simulations, the parameter space of applied pressures and interfacial adhesion energies that inhibit interface failure during cycling are mapped to guide selection of solid-state cells. Furthermore, this study establishes the atomistic modeling for Li stripping and plating, and predicts optimal solid interfaces and new strategies for the future research and development of solid-state Li-metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Parameterization of generic positive sequence models to represent behavior of inverter based resources in low short circuit scenarios

The connection of large inverter based resources (IBR) in transmission systems is often located geographically and electrically far away from load centers. This, coupled with the displacement of synchronous machine plants, results in a reduction of the network short circuit strength at the point of connection. Under these conditions, state-of-the-art positive sequence simulation platforms and models can have difficulties maintaining numerical stability and/or providing an accurate representation of IBR plant dynamic behavior. As a result, computationally heavy time domain electromagnetic transient (EMT) simulations may be required to evaluate these systems. In this paper, a recently developed improved generic positive sequence model is parameterized to provide a representation of IBR behavior in low short circuit networks. Further, comparisons against generic and detailed EMT models demonstrate the suitability of the improved positive sequence model to study practical stability issues experienced with presently in-service plants. Such a model can provide some of the accuracy of an EMT representation with a much lower computational burden. The performance of the positive sequence model is validated against the behavior shown by both open white box and closed black box EMT domain models from around the world.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

High-fidelity forced convection simulations of the University of Wisconsin–Madison air-cooled reactor cavity cooling system

Among the various types of advanced nuclear technologies that exist, high-temperature gas-cooled reactors (HTGRs) stand out for their co-generation capabilities and exceptional passive safety systems. HTGRs rely on a reactor cavity cooling system (RCCS) to remove decay heat from the reactor pressure vessel (RPV) during the course of long-term transients. Several studies have focused on experimental and numerical modeling of the RCCS; however, the literature reveals a notable gap when it comes to high-fidelity simulations such as direct numerical simulations (DNS) or large eddy simulation (LES) models of the RCCS. An LES model was developed for forced convection inside the RCCS so as to obtain a high-fidelity representation of the RCCS at a scale never before attempted. Furthermore, the LES model was compared against three different Reynolds-averaged Navier–Stokes (RANS) models (Realizable k − ϵ, k − ω SST, and Reynolds stress model [RSM]). This comparison verified the robustness and accuracy of the turbulence modeling approach of the RANS models for the RCCS. Finally, based on the comparison between the LES and the RANS models, the k − ω SST RANS model of the experimental facility was developed and compared against the available experimental data.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Neutrino fast flavor instability in three dimensions for a neutron star merger

The flavor evolution of neutrinos in core collapse supernovae and neutron star mergers is a critically important unsolved problem in astrophysics. Following the electron flavor evolution of the neutrino system is essential for calculating the thermodynamics of compact objects as well as the chemical elements they produce. Accurately accounting for flavor transformation in these environments is challenging for a number of reasons, including the large number of neutrinos involved, the small spatial scale of the oscillation, and the nonlinearity of the system. We take a step in addressing these issues by presenting a method which describes the neutrino fields in terms of angular moments. We apply our moment method to neutron star merger conditions and show it simulates fast flavor neutrino transformation in a region where this phenomenon is expected to occur. By comparing with particle-in-cell calculations we show that the moment method is able to capture the three phases of growth, saturation, and decoherence, and correctly predicts the lengthscale of the fastest growing fluctuations in the neutrino field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Physical properties of the quasi-two-dimensional square lattice antiferromagnet Ba 2 FeSi 2 O 7

In this work, we report magnetization (χ, M ), magnetic specific heat ( C M ), and neutron powder diffraction results on a quasi-two-dimensional (2D) S = 2 square lattice antiferromagnet Ba 2 FeSi 2 O 7 consisting of FeO 4 tetrahedrons with highly compressive tetragonal distortion (27%). Despite of the quasi-2D lattice structure, both χ and C M present three-dimensional magnetic long-range ordering below the Néel temperature T N = 5.2 K . Neutron diffraction data show a collinear Q m = (1,0,1/2) antiferromagnetic (AFM) structure below T N but the ordered moment aligned in the a b plane is suppressed by 26% from the ionic spin S = 2 value ( 4 μ B ). Both the AFM structure and the suppressed moments are well explained by using Monte Carlo simulations with a large single-ion in-plane anisotropy D = 1 .4 meV and a rather small Heisenberg exchange J intra = 0 .15 meV in the plane. The characteristic 2D spin fluctuations are recognized in the magnetic entropy release and diffuse scattering above T N . This new quasi-2D magnetic system also displays unusual nonmonotonic dependence of T N as a function of magnetic field H .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Improving collaborations between empiricists and modelers to advance grassland community dynamics in ecosystem models

Climate change, increasing atmospheric CO 2 , and land use change have altered biogeochemical and hydrologic cycles world-wide, with grassland systems being particularly vulnerable to resulting vegetation shifts (Komatsu et al., 2019). Therefore, incorporating plant community dynamics into ecosystem models is critical for accurate forecasting of ecosystem responses to global change (Levine, 2016). Process-based ecosystem models, which simulate the biogeochemical transfers of mass and energy among biota, the subsurface, and atmosphere, require representation of dynamic composition of organisms within ecosystems. For example, these models simulate leaf and plant-level characteristics, such as electron transport rate and allometry of carbon (C) allocation, to predict how net primary productivity and other ecosystem processes respond to abiotic drivers. These models are particularly useful in scaling from organismal to ecosystem levels but are still underdeveloped in their ability to capture community change, especially in grassland ecosystems. To represent compositional changes, these models must simulate competition, mortality, establishment, and reproduction of plant populations within communities. Yet, current ecosystem modeling approaches to forecast plant community change have derived from studies of forested systems and are either too coarse to capture fine-scale community dynamics (e.g. dynamic global vegetation models (DGVMs)) or too complex to be used at large spatial scales (e.g. forest gap models).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

CONTROLLED END GAS AUTO IGNITION WITH EXHAUST GAS RECIRCULATION ON A STOICHIOMETRIC, SPARK IGNITED, NATURAL GAS ENGINE

The goal of this study is to address fundamental limitations to achieving diesel-like efficiencies in heavy duty on-highway natural gas (NG) engines. Engine knock and misfire are barriers to pathways leading to higher efficiency engines. This study explores enabling technologies for development of high efficiency stoichiometric, spark ignited, natural gas engines. These include design strategies for fast and stable combustion and higher dilution tolerance. Additionally, advanced control methodologies are implemented to maintain stable operation between knock and misfire limits. To implement controlled end- gas autoignition (C-EGAI) strategies a Combustion Intensity Metric (CIM) is used for ignition control with the use of a Woodward large engine control module (LECM). Tests were conducted using a single cylinder, variable compression ratio, cooperative fuel research (CFR) engine with baseline conditions of 900 RPM, engine load of 800 kPa indicated mean effective pressure (IMEP), and stoichiometric air/fuel ratio. Exhaust gas recirculation (EGR) tests were performed using a custom EGR system that simulates a high pressure EGR loop and can provide a range of EGR rates from 0 to 40%. The experimental measurements included the variance of EGR rate, compression ratio, engine speed, IMEP, and CIM. These five variables were optimized through a Modified Box- Benken design Surface Response Method (RSM), with brake efficiency as the merit function. A positive linear correlation between CIM and f-EGAI was identified. Consequently, CIM was used as the feedback control parameter for C-EGAI. As such, implementation of C-EGAI effectively allowed for the utilization of high EGR rates and CRs, controlling combustion between a narrower gap between knock and lean limits. The change from fixed to parametric ignition timing with CIM targeted select values of f-EGAI with an average coefficient of variance (COV) of peak pressure of 5.4. The RSM efficiency optimization concluded with operational conditions of 1080 RPM, 1150 kPa IMEP, 10.55:1 compression ratio, and 17.8% EGR rate with a brake efficiency of 21.3%. At this optimized point of peak performance, a f-EGAI for C-EGAI was observed at 34.1% heat release due to auto ignition, a knock onset crank angle value of 10.3° aTDC and ignition timing of - 24.7° aTDC. This work has demonstrated that combustion at a fixed f-EGAI can be maintained through advanced ignition control of CIM without experiencing heavy knocking events.

Cooperative Fuel Research Engine, Auto Ignition, E↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers part I: Computational models and workflow

Computational simulations have become central to the seismic analysis and design of major infrastructure over the past several decades. Most major structures are now “proof tested” virtually through representative simulations of earthquake-induced response. More recently, with the advancement of high-performance computing (HPC) platforms and the associated massively parallel computational ecosystems, simulation is beginning to play a role in increased understanding and prediction of ground motions for earthquake hazard assessments. However, the computational requirements for regional-scale geophysics-based ground motion simulations are extreme, which has restricted the frequency resolution of direct simulations and limited the ability to perform the large number of simulations required to numerically explore the problem parametric space. In this article, recent developments toward an integrated, multidisciplinary earth science-engineering computational framework for the regional-scale simulation of both ground motions and resulting structural response are described with a particular emphasis on advancing simulations to frequencies relevant to engineered systems. This multidisciplinary computational development is being carried out as part of the US Department of Energy (DOE) Exascale Computing Project with the goal of achieving a computational framework poised to exploit emerging DOE exaflop computer platforms scheduled for the 2022–2023 timeframe.

58 GEOSCIENCES↗
DOE OSTI2024

Uncertainty Propagation from Experiment Measurements to Modeling Approaches: A Case for SMR Steam Entrainment Testing

To license new and advanced reactor designs, regulators must be convinced that their unique safety cases—relative to existing large scale reactors—have been adequately addressed by the designed reactor protection systems. In water cooled small modular reactors (SMRs), droplet entrainment in steam flow has significant implications on the progression of accident scenarios due to its compact design features, which requires representative test data applicable to SMR designs. Computer code, modeling and simulation (M&S) tools and models require adequate verification, assessment, and qualification. This includes M&S results validation against scaled empirical data within allowable uncertainty bands to gain regulatory approvals during the various stages of reactor system design, demonstration, and commercialization. However, measurement uncertainty within the empirical datasets and test data applicability ranges requires careful consideration of M&S inputs (i.e., boundary conditions, and initial conditions), and verification and validation efforts. This study focuses on uncertainty quantification in designing scaled test facilities for SMR applications with appropriate measurements and a standard data-reduction method to estimate thermal hydraulics characteristics parameters that incorporate physics phenomena of interest. In addition, this study supports the evaluation model development and assessment process using M&S that interfaces with advanced computing tools and digital twin capabilities. This will allow synchronization between experiment and modeling approaches for droplet entrainment testing and analysis, improving diagnostics, prognostics, and decision-making to accelerate regulatory approval.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

SCUC-DER Integration Report: Integrating Distributed Energy Resources (DER) Using Advanced Unit Commitment Models and DER Aggregation Methodologies

Distributed energy resources (DERs) are continuing to grow due to regulatory, policy, and market shifts, and it needs to be ensured that small DERs are given a level playing field with traditional resources. Legacy market processes such as unit-commitment problems were designed for a power grid consisting largely of centralized power plants. In contrast, DERs consist of many small devices with distinct operating characteristics that may or may not be connected at the transmission interconnection points, limiting their visibility to the independent system operators (ISOs) who operate wholesale electricity markets. This report details the development and initial results from a simulation platform that integrates state-of-the-art security-constrained unit commitment software, detailed feeder models, and a DER aggregator model to quantify potential DER integration issues. Quantitative results to date illustrate potential infeasible scheduling solutions from SCUC when the DERA includes aggregations of energy-limited energy storage resources. Likewise, if aggregations are not penalized for dispatch deviations, they may have incentives to deviate from the SCUC-determined resource schedules. Assumptions about the amount of aggregated demand response resources (DRRs) and the ability of DERAs to follow profit incentives have an important impact – DRRs have significant flexibility and can typically feasibly meet their SCUC schedule, but on the other hand, their profit incentives can cause unscheduled increases in load before and after DRR dispatch. Computation time results on the SCUC solver and simulation platform only show a modest increase in SCUC solution time as the number of DERAs is increased, but results to date only reflect the RTS-GMLC test system; results may show more significant solver slowdown in larger transmission systems. The largest contribution to simulation time is attributed to the DERA offer generation method, which is suggested for future improvements.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Compare low-gain 1D and 3D models toward a storage ring-based X-ray free-electron-laser oscillator

Fourth-generation diffraction-limited synchrotron light sources (DL-SLS) present a compelling platform for advanced light source applications due to their ultra-low emittance (around tens of picometers), sub-micrometer beam stability, and high repetition rates (greater than 1 MHz). Despite these advantages, fundamental limitations, such as long bunch durations (tens of picoseconds), large relative energy spread (~10⁻³), and low peak current (<300 A), have hindered their effectiveness in free-electron laser (FEL) applications. Using a low-gain three-dimensional (3D) model based on NSLS-II upgrade lattice parameters, we demonstrate that a 6% single-pass gain is achievable with an optimized transverse gradient undulator (TGU). Concurrently, recent advances in X-ray optics—specifically, a four-mirror bow-tie resonator utilizing diamond mirrors and beryllium compound refractive lenses (CRLs)—enable total round-trip losses below 4% and an output coupling efficiency near 1%, rendering a Storage Ring based X-ray FEL oscillator (XFELO) a realistic possibility. In the 3D framework, the TGU effectively mitigates the detrimental effects of large energy spread by spatially dispersing the beam, leading to enhanced gain relative to conventional one-dimensional (1D) models. We establish both analytical and numerical connections between the 1D and 3D descriptions, leveraging GENESIS simulations to map equivalent gain conditions. While 1D resonance is defined in the energy domain, the 3D TGU system shifts the resonance dependence to transverse gradient and dispersion parameters. Once optimized, this framework allows predictive modeling using 1D formulas, which can then be validated through full 3D simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Seminar 13 - Advanced Methods for Grid Integration of High-Performance Residential Communities

High-performance residential communities such as net-zero energy (NZE) and NZE-ready communities are emerging across the United States as many cities and states have committed to achieve ambitious energy efficiency and renewable energy goals between now and 2050. The achievement of these goals will transform the nation's communities into the holistic paradigms consisting of pervasive solar photovoltaics (PVs) and other renewable distributed energy resources, energy storage and flexible loads. However, integrating residential communities with high PV penetration may bring many challenges to the grid. For example, during many hours of the year, electric power is flowing out of PV-enabled homes and back into the grid. When this happens simultaneously in many nearby homes, distribution feeder voltage can be raised significantly. Intermittent solar power can also cause voltage flickers, leading to potential damages on both distribution feeder and home appliances. Existing solutions, such as PV curtailment, battery storage, and direct load control, are not able to address all the challenges in a cost-effective way. Better understanding is needed on how to deploy rooftop PV, how to predict building loads, and how to coordinate the homes in a community to enhance grid reliability. This seminar consists of four integral presentations that cover different aspects of a solution to these challenges. First, we will present an automated method to identify the best surfaces to deploy PV and size PV according to the PV penetration goals in a community. Second, we will present a streamlined way to create gray-box building models that are suitable for large-scale building-grid co-simulation. Next, we will present a data-driven method to estimate the end-use energy consumption to better coordinate the behind-the-meter resources to reduce utility cost, provide grid services, and improve renewable energy integration. Lastly, we will present a hierarchical control system that integrates the building models and the estimated non-dispatchable loads into the decision-making process to enhance the reliability and resilience of the residential community. Simulation results from an NZE-ready community in Fort Collins, Colorado will be presented in all four presentations and the guidelines on better grid integration of residential communities will be provided.

14 SOLAR ENERGY↗
DOE OSTI2024

SIMULATION CLONING FOR DIGITAL TWINS: A SCALABLE APPROACH

Digital Twin (DT) methods represent an important technology in which a simulated model of the operations of a physical system uses real-time sensor data to simulate, monitor, and consequently, improve its operations. One of the primary objectives of such a DT is to inform the physical system of measures to be taken in response to one or multiple intervening events that can change the state of the physical system. As such, a capability that is able to carry out multiple scenario assessments in real time in readiness for such events is a very effective tool in the use of simulations as DTs. However, continuous evaluation with highly probable event simulation scenarios are challenging due to the constraints of finite memory and a large exploration space. This paper reports a novel methodology for the continuous evaluation of $k$ probabilistic \textit{what-if} event scenarios under finite resource constraints and demonstrates its use as a digital-twin for a real-world application.

Yoginath, Srikanth↗
DOE OSTI2026

Evaluation of Best Practices in Mitigating Startup Costs on Leadership-Class Supercomputers

Supercomputers at Department of Energy (DOE) National Laboratories face a widening range of workloads, from traditional modeling and simulation to Artificial Intelligence model training or complex multi-stage workflows, and beyond. At DOE Leadership Computing Facilities like the Oak Ridge Leadership Computing Facility (OLCF), these workloads demand concurrent access to large portions of the supercomputer’s resources. Launching a job across massive supercomputers is challenging from the start; the file system struggles with a large backlog of metadata requests as tens of thousands of processes read thousands of the same files, and the compute job cannot start until this is completed. There are multiple existing approaches to calm this metadata storm, ranging from vendor-developed tools like sbcast to National Laboratory-developed tools like Spindle and Copper. In this paper, we benchmark and discuss three common approaches to improving compute job launch latencies on Frontier: Slurm’s sbcast tool, Spindle, and Copper. We evaluate these tools by measuring the launch latencies of four workloads: OSU Microbenchmark’s osu_init, Pynamic, Python import mpi4py, and Python import torch. We provide discussion of the results, highlighting data that meet expectations and that do not meet expectations.

Hagerty, Nick [ORNL] (ORCID:0000000330014414)↗
DOE OSTI2021

A Co-design Framework for Online Data Analysis and Reduction

Science applications preparing for the exascale era are increasingly exploring in situ computations comprising of simulation-analysis-reduction pipelines coupled in-memory. Efficient composition and execution of such complex pipelines for a target platform is a codesign process that evaluates the impact and tradeoffs of various application- and system-specific parameters. In this article, we describe a toolset for automating performance studies of composed HPC applications that perform online data reduction and analysis. We describe Cheetah, a new framework for composing parametric studies on coupled applications, and Savanna, a runtime engine for orchestrating and executing campaigns of codesign experiments. Furthermore, this toolset facilitates understanding the impact of various factors such as process placement, synchronicity of algorithms, and storage versus compute requirements for online analysis of large data. Ultimately, we aim to create a catalog of performance results that can help scientists understand tradeoffs when designing next-generation simulations that make use of online processing techniques. We illustrate the design of Cheetah and Savanna, and present application examples that use this framework to conduct codesign studies on small clusters as well as leadership class supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantifying CanESM5 and EAMv1 sensitivities to Mt. Pinatubo volcanic forcing for the CMIP6 historical experiment

Large volcanic eruptions reaching the stratosphere have caused marked perturbations to the global climate including cooling at the Earth's surface, changes in large-scale circulation and precipitation patterns and marked temporary reductions in global ocean heat content. Many studies have investigated these effects using climate models; however, uncertainties remain in the modelled response to these eruptions. This is due in part to the diversity of forcing datasets that are used to prescribe the distribution of stratospheric aerosols resulting from these volcanic eruptions, as well as uncertainties in optical property derivations from these datasets. To improve this situation for the sixth phase of the Coupled Model Intercomparison Project (CMIP6), a two-step process was undertaken. First, a combined stratospheric aerosol dataset, the Global Space-based Stratospheric Aerosol Climatology (GloSSAC; 1979–2016), was constructed. Next, GloSSAC, along with information from ice cores and Sun photometers, was used to generate aerosol distributions, characteristics and optical properties to construct a more consistent stratospheric aerosol forcing dataset for models participating in CMIP6. This “version 3” of the stratospheric aerosol forcing has been endorsed for use in all contributing CMIP6 simulations. Recent updates to the underlying GloSSAC from version 1 to version 1.1 affected the 1991–1994 period and necessitated an update to the stratospheric aerosol forcing from version 3 to version 4. As version 3 remains the official CMIP6 input, quantification of the impact on radiative forcing and climate is both relevant and timely for interpreting results from experiments such as the CMIP6 historical simulations. This study uses two models, the Canadian Earth System Model version 5 (CanESM5) and the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1), to estimate the difference in instantaneous radiative forcing in simulated post-Pinatubo climate response when using version 4 instead of version 3. Differences in temperature, precipitation and radiative forcings are generally found to be small compared to internal variability. An exception to this is differences in monthly temperature anomalies near 24 km altitude in the tropics, which can be as large as 3 °C following the eruption of Mt. Pinatubo.

58 GEOSCIENCES↗