Testing the magnetohydrodynamic analysis software with natural convection and geodynamo problems

Using the control volume method we developed the software for the numerical solution of viscous incompressible fluid resistive magnetohydrodynamics problems on structured staggered meshes in spherical coordinates. The constrained transport algorithm and the QUICK method with delayed correction for the approximation of the convective terms were used for the discretization of the magnetic field induction equation. The SIMPLER algorithm was applied to solving the hydrodynamic equations. We developed software for modeling natural convection and the hydromagnetic dynamo in a rotating sphere or spherical shell. We proposed an algorithm for the numerical solution of the geodynamo problem with vacuum boundary conditions. The results of solving natural convection and geodynamo benchmark problems with vacuum boundary conditions are presented; they demonstrate a fairly accurate agreement with the reference calculations. The software supports CUDA-enabled accelerators and uses a set of extensions to the Fortran programming language.