# Efficient high order accurate staggered semi-implicit discontinuous Galerkin methods for natural convection problems

## TL;DR

This paper introduces high-order staggered semi-implicit discontinuous Galerkin methods for natural convection, combining explicit and semi-Lagrangian schemes to improve efficiency and stability in simulating buoyancy-driven flows.

## Contribution

It develops a novel class of high-order staggered semi-implicit DG schemes for natural convection, extending previous work to include buoyancy effects and compressible flows.

## Key findings

- Schemes show good agreement with benchmark data.
- Semi-Lagrangian approach enhances computational efficiency.
- Method effectively handles both incompressible and compressible flows.

## Abstract

We propose a new family of high order staggered semi-implicit discontinuous Galerkin (DG) methods for the simulation of natural convection problems. Assuming small temperature fluctuations, the Boussinesq approximation is valid and the flow can simply be modeled by the incompressible Navier-Stokes equations coupled with a transport equation for the temperature and a buoyancy source term in the momentum equation. Our numerical scheme is developed starting from the work presented in [TD14], in which the spatial domain is discretized using a face-based staggered unstructured mesh. For the computation of the advection and diffusion terms, two different algorithms are presented: i) a purely Eulerian explicit upwind-type scheme and ii) a semi-Lagrangian approach. The first methodology leads to a conservative scheme whose major drawback is the time step restriction imposed by the CFL stability condition. On the contrary, computational efficiency can be notably improved relying on a semi-Lagrangian approach. This method leads to an unconditionally stable scheme if the diffusive terms are discretized implicitly. Once the advection and diffusion contributions have been computed, the pressure Poisson equation is solved and the velocity is updated. As a second model for the computation of buoyancy-driven flows, we also consider the full compressible Navier-Stokes equations. The staggered semi-implicit DG method first proposed in [TD17] for all Mach number flows is properly extended to account for the gravity source terms arising in the momentum and energy conservation laws. The validity and the robustness of our novel class of staggered semi-implicit DG schemes is assessed at the aid of several classical benchmark problems, showing in all cases a good agreement with available numerical reference data. Finally, a detailed comparison between the incompressible and the compressible solver is presented.