Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamical core”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A mountain-induced moist baroclinic wave test case for the dynamical cores of atmospheric general circulation models

Idealized test cases for the dynamical cores of atmospheric general circulation models are informative tools to assess the accuracy of the numerical designs and investigate the general characteristics of atmospheric motions. A new test case is introduced that is built upon a baroclinically unstable base state with an added orographic barrier. The topography is analytically prescribed and acts as a trigger of both baroclinic Rossby waves and inertia–gravity waves on a rotating, regular-sized planet. Both dry and idealized moist configurations are suggested. The latter utilizes the Kessler warm-rain precipitation scheme. The test case enhances the complexity of the existing test suite hierarchy and focuses on the impacts of two midlatitudinal mountain ridges on the circulation. Selected simulation examples from four dynamical cores are shown. These are the Spectral Element and Finite Volume dynamical cores, which are part of the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM), versions 2.1.3 and 2.2, and the Cubed-Sphere Finite Volume dynamical cores, which is new to CESM version 2.2. In addition, the Model for Prediction Across Scales (MPAS) is tested. The overall flow patterns agree well in the four dynamical cores, but the details can vary greatly. The examples highlight the broad palette of use cases for the test case and reveal physics–dynamics coupling issues.

58 GEOSCIENCES↗
DOE OSTI2019

Evaluation of Implicit‐Explicit Additive Runge‐Kutta Integrators for the HOMME‐NH Dynamical Core

Abstract The nonhydrostatic High‐Order Method Modeling Environment (HOMME‐NH) atmospheric dynamical core supports acoustic waves that propagate significantly faster than the advective wind speed, thus greatly limiting the time step size that can be used with standard explicit time integration methods. Resolving acoustic waves is unnecessary for accurate climate and weather prediction. This numerical stiffness is addressed herein by considering implicit‐explicit additive Runge‐Kutta (ARK IMEX) methods that can treat the acoustic waves in a stable manner without requiring implicit treatment of nonstiff modes. Various ARK IMEX methods are evaluated for their efficiency in producing accurate solutions, ability to take large time step sizes, and sensitivity to grid cell length ratio. Both the gravity wave test and baroclinic instability test from the 2012 Dynamical Core Model Intercomparison Project are used to recommend 5 of the 27 ARK IMEX methods tested for use in HOMME‐NH.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Comparative Hydrophobic Core Dynamics Between Wild‐Type Amyloid‐β Fibrils, Glutamate‐3 Truncation, and Serine‐8 Phosphorylation

Abstract Post‐translational modifications (PTMs) of amyloid‐β (Aβ) species are implicated in the modulation of overall toxicities and aggregation propensities. We investigated the internal dynamics in the hydrophobic core of the truncated ΔE3 mutant fibrils of Aβ 1–40 and compared them with prior and new data for wild‐type fibrils as well as with phosphorylated S8 fibrils. Deuteron static solid‐state NMR techniques, spanning line‐shape analysis, longitudinal relaxation, and chemical exchange saturation transfer methods, were employed to assess the rotameric jumps of several methyl‐bearing and aromatic groups in the core of the fibrils. Taken together, the results indicate the rather significant influence of the PTMs on the hydrophobic core dynamics, which propagates far beyond the local site of the chemical modification. The phosphorylated S8 fibrils display an overall rigidifying of the core based on the higher activation barriers of motions than the wild‐type fibrils, whereas the ΔE3 fibrils induce a broader variety of changes, some of which are thermodynamic in nature rather than the kinetic ones.

Vugmeyster, Liliya↗
DOE OSTI2025

A GPU‐Based Ocean Dynamical Core for Routine Mesoscale‐Resolving Climate Simulations

Abstract We describe an ocean hydrostatic dynamical core implemented in Oceananigans optimized for Graphical Processing Unit (GPU) architectures. On 64 A100 GPUs, equivalent to 16 computational nodes in current state‐of‐the‐art supercomputers, our dynamical core can simulate a decade of near‐global ocean dynamics per wall‐clock day at an 8‐km horizontal resolution; a resolution adequate to resolve the ocean's mesoscale eddy field. Such efficiency, achieved with relatively modest hardware resources, suggests that climate simulations on GPUs can incorporate fully eddy‐resolving ocean models. This removes a major source of systematic bias in current IPCC coupled model projections, the parameterization of ocean eddies, and represents a major advance in climate modeling. We discuss the computational strategies, focusing on GPU‐specific optimization and numerical implementation details that enable such high performance.

Silvestri, Simone [Massachusetts Institute of Tech↗
DOE OSTI2022

An Assessment of Nonhydrostatic and Hydrostatic Dynamical Cores at Seasonal Time Scales in the Energy Exascale Earth System Model (E3SM)

Abstract In global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central‐southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

58 GEOSCIENCES↗
DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗
DOE OSTI2021

Coupled valence carrier and core-exciton dynamics in WS 2 probed by few-femtosecond extreme ultraviolet transient absorption spectroscopy

Few-femtosecond extreme ultraviolet (XUV) transient absorption spectroscopy, performed with optical 500-1000 nm supercontinuum and broadband XUV pulses (30-50 eV), simultaneously probes dynamics of photoexcited carriers in WS$_{2}$ at the W O$_3$ edge (37-45 eV) and carrier-induced modifications of core-exciton absorption at the W N$_{6,7}$ edge (32-37 eV). Access to continuous core-to-conduction band absorption features and discrete core-exciton transitions in the same XUV spectral region in a semiconductor provides a novel means to investigate the effect of carrier excitation on core-exciton dynamics. The core-level transient absorption spectra, measured with either pulse arriving first to explore both core-level and valence carrier dynamics, reveal that core-exciton transitions are strongly influenced by the photoexcited carriers. A $1.2\pm0.3$ ps hole-phonon relaxation time and a $3.1\pm0.4$ ps carrier recombination time are extracted from the XUV transient absorption spectra from the core-to-conduction band transitions at the W O$_{3}$ edge. Global fitting of the transient absorption signal at the W N$_{6,7}$ edge yields $\sim 10$ fs coherence lifetimes of core-exciton states and reveals that the photoexcited carriers, which alter the electronic screening and band filling, are the dominant contributor to the spectral modifications of core-excitons and direct field-induced changes play a minor role. This work provides a first look at the modulations of core-exciton states by photoexcited carriers and advances our understanding of carrier dynamics in metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Polarization-resolved core exciton dynamics in LiF using attosecond transient absorption spectroscopy

The ability to control absorption by modifying the polarization of light presents an exciting opportunity to experimentally determine the orbital alignment of absorption features. Here, attosecond extreme ultraviolet (XUV) transient absorption spectroscopy is used to investigate the polarization dependence of core exciton dynamics in LiF thin films at the Li + K edge. XUV pulses excite electrons from the Li 1s core level into the conduction band, allowing for the formation of a p-orbital-like core exciton, aligned along the XUV light polarization axis. A sub-5 fs near-infrared (NIR) probe pulse then arrives at variable time delays, perturbing the XUV-excited states and allowing the coherence decay of the core exciton to be mapped. The coherence lifetimes are found to be ≈ 2.4 ± 0.4 fs, which is attributed to a phonon-mediated dephasing mechanism as in previous core exciton studies. The differential absorption features are also shown to be sensitive to the relative polarization of the XUV and NIR fields. The parallel NIR probe induces couplings between the initial XUV-excited p-like bright exciton and s-like dark excitons. When crossed pump and probe polarizations are used, the coupling between the bright and dark states is no longer dipole-allowed, and the transient absorption signal associated with the coupling is suppressed by approximately 90%. This interpretation is supported by simulations of a few-level model system, as well as analysis of the calculated band structure. Here, the results indicate that laser polarization can serve as a powerful experimental tool for exploring the orbital alignment of core excitonic states in solid-state materials.

Attosecond laser spectroscopy↗
DOE OSTI2020

Attosecond spectroscopy reveals alignment dependent core-hole dynamics in the ICl molecule

The removal of electrons located in the core shells of molecules creates transient states that live between a few femtoseconds to attoseconds. Owing to these short lifetimes, time-resolved studies of these states are challenging and complex molecular dynamics driven solely by electronic correlation are difficult to observe. Here, few-femtosecond core-excited state lifetimes of iodine monochloride are obtained by attosecond transient absorption on iodine 4d-16p transitions around 55 eV. Core-level ligand field splitting allows direct access of excited states aligned along and perpendicular to the ICl molecular axis. Lifetimes of 3.5 fs and 4.3 fs are obtained for core-hole states parallel to the bond and 6.5 fs and 6.9 fs for perpendicular states, while nuclear motion is essentially frozen on this timescale. Theory shows that the dramatic decrease of lifetime for core-vacancies parallel to the covalent bond is a manifestation of non-local interactions with the neighboring Cl atom of ICl.

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

DCMIP2016: the tropical cyclone test case

This paper describes and analyzes the Reed–Jablonowski (RJ) tropical cyclone (TC) test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). This intermediate-complexity test case analyzes the evolution of a weak vortex into a TC in an idealized tropical environment. Reference solutions from nine general circulation models (GCMs) with identical simplified physics parameterization packages that participated in DCMIP2016 are analyzed in this study at 50 km horizontal grid spacing, with five of these models also providing solutions at 25 km grid spacing. Evolution of minimum surface pressure (MSP) and maximum 1 km azimuthally averaged wind speed (MWS), the wind–pressure relationship, radial profiles of wind speed and surface pressure, and wind composites are presented for all participating GCMs at both horizontal grid spacings. While all TCs undergo a similar evolution process, some reach significantly higher intensities than others, ultimately impacting their horizontal and vertical structures. TCs simulated at 25 km grid spacings retain these differences but reach higher intensities and are more compact than their 50 km counterparts. These results indicate that dynamical core choice is an essential factor in GCM development, and future work should be conducted to explore how specific differences within the dynamical core affect TC behavior in GCMs.

Willson, Justin L.↗
DOE OSTI2021

Solid state core-exciton dynamics in NaCl observed by tabletop attosecond four-wave mixing spectroscopy

Nonlinear wave mixing in solids with ultrafast x-rays can provide insight into complex electronic dynamics of materials. Here, tabletop-based attosecond noncollinear four-wave mixing (FWM) spectroscopy using one extreme ultraviolet (XUV) pulse from high harmonic generation and two separately timed few-cycle near-infrared (NIR) pulses characterizes the dynamics of the Na + L 2,3 edge core-excitons in NaCl around 33.5 eV. An inhomogeneous distribution of core-excitons underlying the well-known doublet absorption of the Na+ Γ point core-exciton spectrum is deconvoluted by the resonance-enhanced nonlinear wave mixing spectroscopy. In addition, other dark excitonic states that are coupled to the XUV-allowed levels by the NIR pulses are characterized spectrally and temporally. Approximately <10fs coherence lifetimes of the core-exciton states are observed. Furthermore, the core-excitonic properties are discussed in the context of strong electron-hole exchange interactions, electron-electron correlation, and electron-phonon broadening. This investigation successfully indicates that tabletop attosecond FWM spectroscopies represent a viable technique for time-resolved solid state measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Seasonal dynamics of core fungi in the switchgrass phyllosphere, and co-occurrence with leaf bacteria

Plant leaves harbor complex microbial communities that influence plant health and productivity. Nevertheless, a detailed understanding of phyllosphere community assembly and drivers is needed, particularly for phyllosphere fungi. Here, we investigated seasonal dynamics of epiphytic phyllosphere fungal communities in switchgrass ( Panicum virgatum L.), a focal bioenergy crop. We also leverage previously published data on switchgrass phyllosphere bacterial communities from the same experimental plants, allowing us to compare fungal and bacterial dynamics and explore inter-Domain network associations in the switchgrass phyllosphere. Overall, we found a strong impact of sampling date on fungal community composition, with multiple taxonomic levels exhibiting clear temporal patterns in relative abundance. In addition, leaf nitrogen concentration, leaf dry matter content, plant height, and minimum daily air temperature explained significant variation in phyllosphere fungal communities, likely due to their correlation with sampling date. Finally, among the core taxa, fungal-bacterial network associations were much more common than bacteria-bacteria associations, suggesting the importance of inter-Domain phylogenetic diversity in microbiome assembly. Although our findings highlight the complexity of phyllosphere microbiome assembly, the clear temporal patterns in lineage-specific fungal abundances give promise to the potential for accurately predicting shifts in fungal phyllosphere communities throughout the growing season, a key research priority for sustainable agriculture.

16S rRNA↗
DOE OSTI2022

Islet: interpolation semi-Lagrangian element-based transport

Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.

58 GEOSCIENCES↗
DOE OSTI2020

Improving Time Step Convergence in an Atmosphere Model With Simplified Physics: The Impacts of Closure Assumption and Process Coupling

Convergence testing is a common practice in the development of dynamical cores of atmospheric models but is not as often exercised for the parameterization of sub-grid physics. An earlier study revealed that the stratiform cloud parameterizations in several predecessors of the Energy Exascale Earth System Model (E3SM) showed strong time-step sensitivity and slower-than-expected convergence when the model's time step was systematically refined. In this work, a simplified atmosphere model is configured that consists of the spectral-element dynamical core of the E3SM atmosphere model coupled with a large-scale condensation parameterization based on commonly used assumptions. This simplified model also resembles E3SM and its predecessors in the numerical implementation of process coupling and shows poor time-step convergence in short ensemble tests. We present a formal error analysis to reveal the expected time-step convergence rate and the conditions for obtaining such convergence. Numerical experiments are conducted to investigate the root causes of convergence problems. We show that revisions in the process coupling and closure assumption help to improve convergence in short simulations using the simplified model; the same revisions applied to a full atmosphere model lead to significant changes in the simulated long-term climate. This work demonstrates that causes of convergence issues in atmospheric simulations can be understood by combining analyses from physical and mathematical perspectives. Addressing convergence issues can help to obtain a discrete model that is more consistent with the intended representation of the physical phenomena.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sensitivity of the latitude of the westerly jet stream to climate forcing

The latitude of the westerly jet stream is influenced by a variety of climate forcings. Their effects often manifest as a tug-of-war between tropical forcing, characterized by tropical upper-tropospheric warming, and polar forcing, with a signature of Antarctic stratospheric cooling or Arctic amplification. Here we present a unified forcing-feedback frame work relating different climate forcings to their forced jet changes, in which the interactions between the westerly jet and synoptic storms are synthesized by a zonal advection feedback, analogous to the feedback framework for assessing climate sensitivity. This framework is supported by a prototype feedback analysis in the atmospheric dynamical core of a climate model with diverse thermal and mechanical forcings. Our analysis indicates that the latitude of a westerly jet is most sensitive to the thermally-induced jet speed changes near the tropopause. The jet shift also displays a larger sensitivity than linearity for the equatorward shift.

atmospheric dynamics, atmospheric variability, cli↗