Comparison of Two Numerical Models of Convection in the Earth’s Mantle

The three-dimensional convection in the Earth’s mantle is studied with a well-known mathematical model, which includes the Navier–Stokes equations in the Oberbeck–Boussinesq and geodynamic approximations. Two numerical models of convection are considered. The first is based on the implicit finite-difference schemes of splitting over spatial variables with correction of pressure. The second numerical model is based on the spectral difference method. The numerical models constructed were compared on model problems of convection in a rectangular parallelepiped in a liquid with constant viscosity, corresponding to the convection in the entire mantle of the Earth [1]. The calculation results are in good agreement with the test results.