High-order Simulation of Nearly Hydrostatic Atmospheric Flows
Abstract
Accurate simulation of atmospheric flows is challenging. Low-Mach-number dynamics often need to be resolved in the presence of a dominant, nearly hydrostatic background state. Conventional upwind numerical fluxes can introduce excessive acoustic-speed-based dissipation in the low-Mach regime. Discretization errors in the balance law can also generate spurious motion. This work investigates these two issues within a high-order Direct Flux Reconstruction (DFR) framework coupled with a third-order exponential time integrator for the 3D compressible Euler equations with gravity. Low-Mach-number accuracy is addressed using a modified AUSM+-up flux, denoted AUSM+-upH. Preservation of the background equilibrium is improved through a deviation-based formulation. The resulting approach is assessed using a hierarchy of idealized atmospheric test cases. The results show that AUSM+-upH preserves low-Mach flow structures more effectively than the more dissipative Rusanov flux. However, reduced numerical dissipation also makes AUSM+-upH more sensitive to residual errors in the discrete hydrostatic balance. The deviation formulation substantially suppresses these errors by several orders of magnitude, with the strongest improvement observed for flows over orography. Overall, AUSM+-upH exhibits greater robustness when combined with the deviation formulation across the atmospheric test cases considered.
Related
- Hierarchical Fractional-Order Disturbance-Rejection Control of an Electro-Hydrostatic Actuator for Marine Steering Gear
- Amortized Probabilistic Retrieval of Atmospheric CO2 from OCO-2 Spectra Using Deep Learning with Laplace Approximations and Normalizing Flows
- Well-balanced second-order approximation of the compressible atmospheric Euler equations
- Conditional Mean Flows: Few-Step Densities for Simulation-Based Inference
- A Massively Parallel Hybridizable Discontinuous Galerkin Solver for Direct Numerical Simulation of Compressible Flows on GPUs
BibTeX
@misc{panday2026highorder,
title = {{High-order Simulation of Nearly Hydrostatic Atmospheric Flows}},
author = {Shoyon Panday and Stéphane Gaudreault and Carlos Pereira and Philipp Birken},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.31836},
eprint = {2610.31836},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.31836},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - High-order Simulation of Nearly Hydrostatic Atmospheric Flows
AU - Shoyon Panday
AU - Stéphane Gaudreault
AU - Carlos Pereira
AU - Philipp Birken
PY - 2026
DA - 2026/10/22
PB - arXiv
DO - 10.48550/arxiv.2610.31836
UR - https://arxiv.org/abs/2610.31836
AB - Accurate simulation of atmospheric flows is challenging. Low-Mach-number dynamics often need to be resolved in the presence of a dominant, nearly hydrostatic background state. Conventional upwind numerical fluxes can introduce excessive acoustic-speed-based dissipation in the low-Mach regime. Discretization errors in the balance law can also generate spurious motion. This work investigates these two issues within a high-order Direct Flux Reconstruction (DFR) framework coupled with a third-order exponential time integrator for the 3D compressible Euler equations with gravity. Low-Mach-number accuracy is addressed using a modified AUSM^+-up flux, denoted AUSM^+-upH. Preservation of the background equilibrium is improved through a deviation-based formulation. The resulting approach is assessed using a hierarchy of idealized atmospheric test cases. The results show that AUSM^+-upH preserves low-Mach flow structures more effectively than the more dissipative Rusanov flux. However, reduced numerical dissipation also makes AUSM^+-upH more sensitive to residual errors in the discrete hydrostatic balance. The deviation formulation substantially suppresses these errors by several orders of magnitude, with the strongest improvement observed for flows over orography. Overall, AUSM^+-upH exhibits greater robustness when combined with the deviation formulation across the atmospheric test cases considered.
ER -