loading page

A pointwise conservative method for thermochemical convection under the compressible anelastic liquid approximation
  • Nathan Sime,
  • Cian R. Wilson,
  • Peter van Keken
Nathan Sime
Carnegie Institution for Science

Corresponding Author:[email protected]

Author Profile
Cian R. Wilson
Carnegie Institution for Science
Author Profile
Peter van Keken
Carnegie Institution for Science
Author Profile

Abstract

In prior work we found that precise approximation of the continuity constraint is crucial for accurate propagation of tracer data when advected through a background incompressible velocity field (Sime et al., 2021). Here we extend this investigation to compressible flows using the anelastic liquid approximation (ALA) and address four related issues: 1. exact conservation of tracer discretized fields through a background compressible velocity; 2. exact mass conservation; 3. addition and removal of tracers without affecting (exact) conservation to preserve a consistent number of tracers per cell; and 4. the diffusion of tracer data, for example, as induced by thermal or chemical effects. In this process we also present an abstract formulation of the interior penalty hybrid discontinuous Galerkin (HDG) finite element formulation for diffusion problems, and apply it to the advection-diffusion and compressible Stokes systems. Finally we present numerical experiments exhibiting the HDG compressible Stokes momentum formulation’s superconvergent compressibility approximation and reproduce community numerical benchmarks from the literature for the ALA.
Feb 2022Published in Geochemistry, Geophysics, Geosystems volume 23 issue 2. 10.1029/2021GC009922