A hybrid discrete exterior calculus and finite difference method for anelastic convection in spherical shells

Hamid Hassan Khan,Pankaj Jagad,Matteo Parsani

2023-11-12

Abstract:The present work develops, verifies, and benchmarks a hybrid discrete exterior calculus and finite difference (DEC-FD) method for density-stratified thermal convection in spherical shells. Discrete exterior calculus (DEC) is notable for its coordinate independence and structure preservation properties. The hybrid DEC-FD method for Boussinesq convection has been developed by Mantravadi et al. (Mantravadi, B., Jagad, P., & Samtaney, R. (2023). A hybrid discrete exterior calculus and finite difference method for Boussinesq convection in spherical shells. Journal of Computational Physics, 491, 112397). Motivated by astrophysics problems, we extend this method assuming anelastic convection, which retains density stratification; this has been widely used for decades to understand thermal convection in stars and giant planets. In the present work, the governing equations are splitted into surface and radial components and discrete anelastic equations are derived by replacing spherical surface operators with DEC and radial operators with FD operators. The novel feature of this work is the discretization of anelastic equations with the DEC-FD method and the assessment of a hybrid solver for density-stratified thermal convection in spherical shells. The discretized anelastic equations are verified using the method of manufactured solution (MMS). We performed a series of three-dimensional convection simulations in a spherical shell geometry and examined the effect of density ratio on convective flow structures and energy dynamics. The present observations are in agreement with the benchmark models.

Fluid Dynamics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to develop, validate, and evaluate a hybrid discrete exterior calculus and finite - difference method (DEC - FD) for density - stratified thermal convection in a spherical shell. Specifically, the paper attempts to solve the following key problems:

1. **Numerical simulation of density - stratified thermal convection**:
- Develop a new hybrid discrete exterior calculus and finite - difference method (DEC - FD) for simulating density - stratified thermal convection in a spherical shell.
- This method combines the coordinate - independence and structure - preserving properties of discrete exterior calculus (DEC) with the efficiency of the finite - difference method (FD) in the radial direction.

2. **Validation and benchmark testing of the method**:
- Use the Method of Manufactured Solutions (MMS) to validate the discretized anelastic equations.
- Examine the influence of the density ratio on the convective flow structure and energy dynamics through a series of convection simulations in three - dimensional spherical - shell geometries.
- Compare the simulation results with existing benchmark models to ensure the accuracy and reliability of the method.

3. **Application of the anelastic approximation**:
- Extend the existing hybrid DEC - FD method to make it applicable to anelastic convection, retaining the density - stratified effect while filtering out sound waves in fully compressible flows.
- The anelastic approximation is widely used in understanding thermal convection in stars and giant planets, so this method has important applications in astrophysics.

### Formula summary

1. **Continuity equation**:
\[
\nabla\cdot(\rho_a\mathbf{u}) = 0
\]

2. **Momentum equation**:
\[
\frac{\partial\mathbf{u}}{\partial t}+(\mathbf{u}\cdot\nabla)\mathbf{u}=-\frac{\nabla p}{\rho_a}+Ra Pr\mathbf{g}(r)s+\frac{Pr}{\rho_a}\tau
\]
where \(\tau\) is the viscous stress tensor:
\[
\tau=\nabla^2\mathbf{u}+\frac{1}{3}\nabla(\nabla\cdot\mathbf{u})
\]

3. **Energy equation**:
\[
\frac{\partial s}{\partial t}+(\mathbf{u}\cdot\nabla)s=\chi(r)u_r+\frac{1}{\rho_a}\nabla^2s+\Phi
\]
where the viscous heating term \(\Phi\) is expressed as:
\[
\Phi=\frac{2\epsilon}{Ra T_a\rho_a}\left(e_{ij}-\frac{1}{3}(\nabla\cdot\mathbf{u})\delta_{ij}\right)^2
\]
here \(e_{ij}=\frac{1}{2}(\partial_j u_i+\partial_i u_j)\) is the strain - rate tensor, and \(\delta_{ij}\) is the Kronecker delta.

### Conclusion

This paper successfully solves the problem of numerical simulation of density - stratified thermal convection in a spherical shell by developing and validating the hybrid DEC - FD method. This method not only retains the coordinate - independence and structure - preserving properties of discrete exterior calculus but also achieves efficient computation in the radial direction through the finite - difference method. Through a series of simulations and benchmark tests, the application potential of this method in astrophysics has been demonstrated.