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

Abstract

Global non-hydrostatic atmospheric models are becoming increasingly important for studying the climates of planets and exoplanets. However, such models suffer from computational difficulties due to the large aspect ratio between the horizontal and vertical directions. To overcome this problem, we developed a global model using a vertically-implicit-correction (VIC) scheme in which the integration time step is no longer limited by the propagation of acoustic waves in the vertical. We proved that our model, based on the $\mathrm{Athen}{\mathrm{a}}^{++}$ framework and its extension for planetary atmospheres - SNAP (Simulating Non-hydrostatic Atmosphere on Planets), rigorously conserves mass and energy in finite volume simulations. We found that traditional numerical stabilizers such as hyper-viscosity and divergence damping are not needed when using the VIC scheme, which greatly simplifies the numerical implementation and improves stability. We present simulation results ranging from 1D linear waves to 3D global circulations with and without the VIC scheme. These tests demonstrate that our formulation correctly tracks local turbulent motions, produces Kelvin-Helmholtz instability, and generates a super-rotating jet on hot Jupiters. Employing this VIC scheme improves the computational efficiency of global simulations by more than two orders of magnitude compared to an explicit model and facilitates the capability of simulating a wide range of planetary atmospheres both regionally and globally.

26 pages, 12 figures, accepted for publication in ApJ

References in corpus (9)

• Athena: A New Code for Astrophysical MHD
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Atmospheric Circulation of Hot Jupiters: A Shallow Three-Dimensional Model
• The Atmospheric Circulation of Ultra-hot Jupiters
• Idealised simulations of the deep atmosphere of hot jupiters: Deep, hot, adiabats as a robust solution to the radius inflation problem
• Atmospheric convection plays a key role in the climate of tidally-locked terrestrial exoplanets: insights from high-resolution simulations
• THOR: A New and Flexible Global Circulation Model to Explore Planetary Atmospheres
• Simulating Non-hydrostatic atmospheres on Planets (SNAP): formulation, validation and application to the Jovian atmosphere
• Extremely long convergence times in a 3D GCM simulation of the sub-Neptune Gliese 1214b

Cited by in corpus (4)

• The THOR+HELIOS general circulation model: multi-wavelength radiative transfer with accurate scattering by clouds/hazes
• A new convective parameterization applied to Jupiter: implications for water abundance near the $24^{\circ }$ N region
• Nonuniform Water Distribution in Jupiter's Mid Latitudes: Influence of Precipitation and Planetary Rotation
• Cloudy and Cloud-free Thermal Phase Curves with PICASO: Applications to WASP-43b