Abstract:A new non-hydrostatic and cloud-resolving atmospheric model is developed for studying moist convection and cloud formation in planetary atmospheres. It is built on top of the Athena++ framework, utilizing its static/adaptive mesh-refinement, parallelization, curvilinear geometry, and dynamic task scheduling. We extend the original hydrodynamic solver to vapors, clouds, and precipitation. Microphysics is formulated generically so that it can be applied to both Earth and Jovian planets. We implemented the Low Mach number Approximate Riemann Solver (LMARS) for simulating low speed atmospheric flows in addition to the usual Roe and HLLC Riemann solvers. Coupled with a fifth-order Weighted Essentially Nonoscillatory (WENO) subgrid-reconstruction method, the sharpness of critical fields such as clouds is well-preserved, and no extra hyperviscosity or spatial filter is needed to stabilize the model. Unlike many atmospheric models, total energy is used as the prognostic variable of the thermodynamic equation. One significant advantage of using total energy as a prognostic variable is that the entropy production due to irreversible mixing process can be properly captured. The model is designed to provide a unified framework for exploring planetary atmospheres across various conditions, both terrestrial and Jovian. First, a series of standard numerical tests for Earth's atmosphere is carried out to demonstrate the performance and robustness of the new model. Second, simulation of an idealized Jovian atmosphere in radiative-convective equilibrium shows that 1) the temperature gradient is superadiabatic near the water condensation level because of the changing of the mean molecular weight, and 2) the mean profile of ammonia gas shows a depletion in the subcloud layer down to nearly 10 bars. Relevance to the recent Juno observations is discussed.

A Global Non-Hydrostatic Atmospheric Model with a Mass and Energy Conserving Vertically-Implicit-Correction (VIC) Scheme

Simulating Non-hydrostatic atmospheres on Planets (SNAP): formulation, validation and application to the Jovian atmosphere

A Three-Dimensional Finite Volume Method for Conservation Laws in Conjunction with Modified Solution for Nonlinear Coupled Constitutive Relations

An Undecomposed Hybrid Algorithm for Nonlinear Coupled Constitutive Relations of Rarefied Gas Dynamics

Numerical Simulation of Hypersonic Reaction Flows with Nonlinear Coupled Constitutive Relations

Numerical Investigation on Local Non-equilibrium Flows Using a Diatomic Nonlinear Constitutive Model

A New Coupled Computational Method in Conjunction with Three-Dimensional Finite Volume Schemes for Nonlinear Coupled Constitutive Relations

Modified Nonlinear Coupled Constitutive Relations Model for Hypersonic Nonequilibrium Flows

ExoCubed: A Riemann-Solver based Cubed-Sphere Dynamic Core for Planetary Atmospheres

Review on development of a scalable high-order nonhydrostatic multi-moment constrained finite volume dynamical core

THOR: A New and Flexible Global Circulation Model to Explore Planetary Atmospheres

A global nonhydrostatic dynamical core on cubed sphere using multi-moment finite volume method: formulation and preliminary test

A semi-implicit compressible model for atmospheric flows with seamless access to soundproof and hydrostatic dynamics

ExoCubed: A Riemann-solver-based Cubed-sphere Dynamic Core for Planetary Atmospheres

Algorithm Development for Vertical Lagrangian Coordinate Based Non-Hydrostatic and Hydrostatic Models

A Performance Study of Horizontally Explicit Vertically Implicit (HEVI) Time-Integrators for Non-Hydrostatic Atmospheric Models

An energetically balanced, quasi-Newton integrator for non-hydrostatic vertical atmospheric dynamics

Supercells and Tornado‐Like Vortices in an Idealized Global Atmosphere Model

Improving the parametrization of transport and mixing processes in planetary atmospheres: the importance of implementing the full Coriolis acceleration

Eigenvectors, Circulation and Linear Instabilities for Planetary Science in 3 Dimensions (ECLIPS3D)

Exploring the effects of numerical methods and slope limiters in heliospheric modeling