Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmospheric 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

The ICON-1.2 Hydrostatic Atmospheric Dynamical Core on Triangular Grids - Part 1: Formulation and Performance of the Baseline Version

A hydrostatic atmospheric dynamical core is developed for the purpose of global climate modelling. The model applies a finite difference method to discretize the primitive equations on spherical icosahedral grids, using C-type staggering with triangles as control volume for mass. This paper documents the numerical method employed in the baseline version of this model, discusses its properties, and presents results from idealized test cases. The evaluation shows that the new dynamical core is able to correctly represent the evolution of baroclinic eddies in the atmosphere and their role in meridional heat and momentum transport. The simulations compare well with the reference solutions, and converge as the horizontal resolution increases. First results from two aqua planet experiments are also presented, in which the equatorial wave spectra derived from tropical precipitation agree well with those simulated by a spectral transform model. The new dynamical core thus provides a good basis for further model development. Certain aspects of the model formulation that need further investigation and improvement are also pointed out.

atmospheric↗

HOMMEXX 1.0: a performance-portable atmospheric dynamical core for the Energy Exascale Earth System Model

Abstract. We present an architecture-portable and performant implementation of the atmospheric dynamical core (High-Order Methods Modeling Environment, HOMME) of the Energy Exascale Earth System Model (E3SM). The original Fortran implementation is highly performant and scalable on conventional architectures using the Message Passing Interface(MPI) and Open MultiProcessor (OpenMP) programming models.We rewrite the model in C++ and use the Kokkos library to express on-node parallelism in a largely architecture-independent implementation. Kokkos provides an abstraction of a compute node or device,layout-polymorphic multidimensional arrays, and parallel execution constructs. The new implementation achieves the same or better performance on conventional multicore computers and is portable to GPUs. We present performance data for the original and new implementations on multiple platforms, on up to 5400 compute nodes, and study several aspects of the single- and multi-node performance characteristics of the new implementation on conventional CPU (e.g., Intel Xeon), many core CPU (e.g., Intel Xeon Phi Knights Landing),and Nvidia V100 GPU.

58 GEOSCIENCES↗

HOMMEXX 1.0: A Performance Portable Atmospheric Dynamical Core for the Energy Exascale Earth System Model

Abstract. We present an architecture-portable and performant implementation of the atmospheric dynamical core (HOMME) of the Energy Exascale Earth System Model (E3SM). The original Fortran implementation is highly performant and scalable on conventional architectures using MPI and OpenMP. We rewrite the model in C++ and use the Kokkos library to express on-node parallelism in a largely architecture-independent implementation. Kokkos provides an abstraction of a compute node or device, layout-polymorphic multidimensional arrays, and parallel execution constructs. Here, the new implementation achieves the same or better performance on conventional multicore computers and is portable to GPUs. We present performance data for the original and new implementations on multiple platforms, on up to 5400 compute nodes, and study several aspects of the single- and multi-node performance characteristics of the new implementation on conventional CPU, Intel Xeon Phi Knights Landing, and Nvidia V100 GPU.

54 ENVIRONMENTAL SCIENCES↗

A Petascale Non-Hydrostatic Atmospheric Dynamical Core in the HOMME Framework

The High-Order Method Modeling Environment (HOMME) is a framework for building scalable, conserva- tive atmospheric models for climate simulation and general atmospheric-modeling applications. Its spatial discretizations are based on Spectral-Element (SE) and Discontinuous Galerkin (DG) methods. These are local methods employing high-order accurate spectral basis-functions that have been shown to perform well on massively parallel supercomputers at any resolution and scale particularly well at high resolutions. HOMME provides the framework upon which the CAM-SE community atmosphere model dynamical-core is constructed. In its current incarnation, CAM-SE employs the hydrostatic primitive-equations (PE) of motion, which limits its resolution to simulations coarser than 0.1 per grid cell. The primary objective of this project is to remove this resolution limitation by providing HOMME with the capabilities needed to build nonhydrostatic models that solve the compressible Euler/Navier-Stokes equations.

58 GEOSCIENCES↗

A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models

A benchmark calculation is proposed for evaluating the dynamical cores of atmospheric general circulation models (GCMs) independently of the physical parameterizations. The test focuses on the long-term statistical properties of a fully developed general circulation; thus, it is particularly appropriate for intercomparing the dynamics used in climate models. To illustrate the use of this benchmark, two very different atmospheric dynamical cores--one spectral, one finite difference--are compared. It is found that the long-term statistics produced by the two models are very similar. Selected results from these calculations are presented to initiate the intercomparison.

Held, Isaac M.↗

A moist aquaplanet variant of the Held–Suarez test for atmospheric model dynamical cores

Abstract. A moist idealized test case (MITC) for atmospheric model dynamical cores is presented. The MITC is based on the Held–Suarez (HS) test that was developed for dry simulations on “a flat Earth” and replaces the full physical parameterization package with a Newtonian temperature relaxation and Rayleigh damping of the low-level winds. This new variant of the HS test includes moisture and thereby sheds light on the nonlinear dynamics–physics moisture feedbacks without the complexity of full-physics parameterization packages. In particular, it adds simplified moist processes to the HS forcing to model large-scale condensation, boundary-layer mixing, and the exchange of latent and sensible heat between the atmospheric surface and an ocean-covered planet. Using a variety of dynamical cores of the National Center for Atmospheric Research (NCAR)'s Community Atmosphere Model (CAM), this paper demonstrates that the inclusion of the moist idealized physics package leads to climatic states that closely resemble aquaplanet simulations with complex physical parameterizations. This establishes that the MITC approach generates reasonable atmospheric circulations and can be used for a broad range of scientific investigations. This paper provides examples of two application areas. First, the test case reveals the characteristics of the physics–dynamics coupling technique and reproduces coupling issues seen in full-physics simulations. In particular, it is shown that sudden adjustments of the prognostic fields due to moist physics tendencies can trigger undesirable large-scale gravity waves, which can be remedied by a more gradual application of the physical forcing. Second, the moist idealized test case can be used to intercompare dynamical cores. These examples demonstrate the versatility of the MITC approach and suggestions are made for further application areas. The new moist variant of the HS test can be considered a test case of intermediate complexity.

54 ENVIRONMENTAL SCIENCES↗

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↗

The Flux-Differencing Discontinuous Galerkin Method Applied to an Idealized Fully Compressible Nonhydrostatic Dry Atmosphere

Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite-volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux-Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high-order accuracy and stability for large-eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge-layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez

54 ENVIRONMENTAL SCIENCES↗

Dynamical Core Model Intercomparison Project (DCMIP) tracer transport test results for CAM‐SE

The Dynamical Core Model Intercomparison Project (DCMIP) provides a set of tests and procedures designed to facilitate development and intercomparison of atmospheric dynamical cores in general circulation models (GCMs). Test category 1 examines the advective transport of passive tracers by three‐dimensional prescribed wind velocity fields, on the sphere. These tests are applied to the Spectral Element (SE) dynamical core of the Community Atmosphere Model (CAM), the default for high‐resolution simulations in the Community Earth System Model (CESM). Test case results are compared with results from the CAM‐FV (Finite Volume) and MCore models where possible. This analysis serves both to evaluate the performance of CAM‐SE's spectral‐element tracer transport routines as well as to provide a baseline for comparison with other atmospheric dynamical cores and f or future improvements to CAM‐SE itself.

Hall, David M.↗

A Finite Volume Scheme on the Cubed Sphere Grid

The performance of a multidimensional finite-volume scheme for global atmospheric dynamics is evaluated on the cubed-sphere geometry. We will explore the properties of the finite volume scheme through traditional advection and shallow water test cases. Baroclinic evaluations performed via a recently developed deterministic initial value baroclinic test case from Jablonowski and Williamson that assesses the evolution of an idealized baroclinic wave in the Northern Hemisphere for a global 3-dimensional atmospheric dynamical core. Comparisons will be made when available to the traditional latitude longitude discretization of the finite-volume dynamical core, as well as other traditional gridpoint and spectral formulations for atmospheric dynamical cores.

Putman, William M.↗

DCMIP2016: a review of non-hydrostatic dynamical core design and intercomparison of participating models

Abstract. Atmospheric dynamical cores are a fundamental component of global atmospheric modeling systems and are responsible for capturing the dynamical behavior of the Earth's atmosphere via numerical integration of the Navier–Stokes equations. These systems have existed in one form or another for over half of a century, with the earliest discretizations having now evolved into a complex ecosystem of algorithms and computational strategies. In essence, no two dynamical cores are alike, and their individual successes suggest that no perfect model exists. To better understand modern dynamical cores, this paper aims to provide a comprehensive review of 11 non-hydrostatic dynamical cores, drawn from modeling centers and groups that participated in the 2016 Dynamical Core Model Intercomparison Project (DCMIP) workshop and summer school. This review includes a choice of model grid, variable placement, vertical coordinate, prognostic equations, temporal discretization, and the diffusion, stabilization, filters, and fixers employed by each system.

54 ENVIRONMENTAL SCIENCES↗

Introducing Enabling Computational Tools to the Climate Sciences: Multi-Resolution Climate Modeling with Adaptive Cubed-Sphere Grids

The research investigates and advances strategies how to bridge the scale discrepancies between local, regional and global phenomena in climate models without the prohibitive computational costs of global cloud-resolving simulations. In particular, the research explores new frontiers in computational geoscience by introducing high-order Adaptive Mesh Refinement (AMR) techniques into climate research. AMR and statically-adapted variable-resolution approaches represent an emerging trend for atmospheric models and are likely to become the new norm in future-generation weather and climate models. The research advances the understanding of multi-scale interactions in the climate system and showcases a pathway how to model these interactions effectively with advanced computational tools, like the Chombo AMR library developed at the Lawrence Berkeley National Laboratory. The research is interdisciplinary and combines applied mathematics, scientific computing and the atmospheric sciences. In this research project, a hierarchy of high-order atmospheric models on cubed-sphere computational grids have been developed that serve as an algorithmic prototype for the finite-volume solution-adaptive Chombo-AMR approach. The foci of the investigations have lied on the characteristics of both static mesh adaptations and dynamically-adaptive grids that can capture flow fields of interest like tropical cyclones. Six research themes have been chosen. These are (1) the introduction of adaptive mesh refinement techniques into the climate sciences, (2) advanced algorithms for nonhydrostatic atmospheric dynamical cores, (3) an assessment of the interplay between resolved-scale dynamical motions and subgrid-scale physical parameterizations, (4) evaluation techniques for atmospheric model hierarchies, (5) the comparison of AMR refinement strategies and (6) tropical cyclone studies with a focus on multi-scale interactions and variable-resolution modeling. The results of this research project demonstrate significant advances in all six research areas. The major conclusions are that statically-adaptive variable-resolution modeling is currently becoming mature in the climate sciences, and that AMR holds outstanding promise for future-generation weather and climate models on high-performance computing architectures.

58 GEOSCIENCES↗

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↗

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↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Subgrid-Scale Mixing in Climate Models: A Novel Look at Diffusion, Accuracy, Stability and Climate Sensitivity

This project focuses on evaluating the role of subgrid-scale dissipation in the dynamical core of atmospheric models. All dynamical cores of atmospheric general circulation models (GCMs) employ some form of subgrid-scale dissipation, either explicitly specified or inherent in the chosen numerical schemes. The dissipation processes are needed to keep the simulation stable or to satisfy important physical properties, and the hope is that they capture and mimic in some poorly understood way the true processes at the unresolved subgrid scale. There is no physical basis that such dissipation can accomplish this. We originally posed a set of numerical test cases chosen or designed to isolate the role of the filters and fixers on both the dynamical variables (pressure, temperature, velocity, vorticity) and trace constituents. From these test case results, we synthesize the information to determine the impact of the subgrid-scale assumptions on weather and climate models.

54 ENVIRONMENTAL SCIENCES↗

AMIP Simulation with the CAM4 Spectral Element Dynamical Core

We evaluate the climate produced by the Community Earth System Model, version 1, running with the new spectral-element atmospheric dynamical core option. The spectral-element method is congured to use a cubed-sphere grid, providing quasi-uniform resolution over the sphere, increased parallel scalability and removing the need for polar filters. It uses a fourth order accurate spatial discretization which locally conserves mass and moist total energy. Using the Atmosphere Model Intercomparison Project protocol, we compare the results from the spectral-element dy- namical core with those produced by the default nite-volume dynamical core and with observations.

CAM-SE↗