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.

arXiv comment26 pages, 12 figures, accepted for publication in ApJ

References in corpus (25)

• Athena: A New Code for Astrophysical MHD
• Atmospheric Circulation and Tides of "51Peg b-like" Planets
• Atmospheric circulation of hot Jupiters: Coupled radiative-dynamical general circulation model simulations of HD 189733b and HD 209458b
• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• The 3.6-8.0 Micron Broadband Emission Spectrum of HD 209458b: Evidence for an Atmospheric Temperature Inversion
• Atmospheric circulation of tidally locked exoplanets: a suite of benchmark tests for dynamical solvers
• Atmospheric Circulation of Hot Jupiters: A Shallow Three-Dimensional Model
• Increased Heat Transport in Ultra-Hot Jupiter Atmospheres Through H${}_2$ Dissociation/Recombination
• Global Evolution of an Accretion Disk with Net Vertical Field: Coronal Accretion, Flux Transport, and Disk Winds
• Formation of Jets and Equatorial Superrotation on Jupiter
• The Atmospheric Circulation of Ultra-hot Jupiters
• Mechanisms of jet formation on the giant planets
• Generation of equatorial jets by large-scale latent heating on the giant planets
• Idealised simulations of the deep atmosphere of hot jupiters: Deep, hot, adiabats as a robust solution to the radius inflation problem
• Results from a set of three-dimensional numerical experiments of a hot Jupiter atmosphere
• 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
• Shear-driven instabilities and shocks in the atmospheres of hot Jupiters
• Global climate modeling of Saturn's atmosphere. Part II: multi-annual high-resolution dynamical simulations
• The Limits of the Primitive Equations of Dynamics for Warm, Slowly Rotating Small Neptunes and Super Earths
• Simulating Non-hydrostatic atmospheres on Planets (SNAP): formulation, validation and application to the Jovian atmosphere
• A High-performance Atmospheric Radiation Package: with applications to the radiative energy budgets of giant planets
• THOR 2.0: Major Improvements to the Open-Source General Circulation Model
• Extremely long convergence times in a 3D GCM simulation of the sub-Neptune Gliese 1214b
• Benchmark experiments with global climate models applicable to extra-solar gas giant planets in the shallow atmosphere approximation

Cited by in corpus (5)

• Aerosols in Exoplanet Atmospheres
• 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