A direction splitting scheme for Navier–Stokes–Boussinesq system in spherical shell geometries

Abstract

AbstractThis article introduces a second‐order direction splitting method for solving the incompressible Navier–Stokes–Boussinesq system in a spherical shell region. The equations are solved on overset Yin–Yang grids, combined with spherical coordinate transforms. This approach allows to avoid the singularities at the poles and keeps the grid size relatively uniform. The downside is that the spherical shell is subdivided into two equally sized, overlapping subdomains that requires the use of Schwarz‐type iterations. The temporal second‐order accuracy is achieved via an artificial compressibility scheme with bootstrapping. The spatial discretization is based on second‐order finite differences on the Marker‐And‐Cell stencil. The entire scheme is implemented in parallel using a domain decomposition iteration and a direction splitting approach for the local solves. The stability, accuracy, and weak scalability of the method is verified on a manufactured solution of the Navier–Stokes–Boussinesq system while its practicality is demonstrated on the natural convection problem in the gap between two concentric spheres.