Linear Full Decoupling, Velocity Correction Method for Unsteady Thermally Coupled Incompressible Magneto-Hydrodynamic Equations

We propose and analyze an effective decoupling algorithm for unsteady thermally coupled magneto-hydrodynamic equations in this paper. The proposed method is a first-order velocity correction projection algorithms in time marching, including standard velocity correction and rotation velocity correction, which can completely decouple all variables in the model. Meanwhile, the schemes are not only linear and only need to solve a series of linear partial differential equations with constant coefficients at each time step, but also the standard velocity correction algorithm can produce the Neumann boundary condition for the pressure field, but the rotational velocity correction algorithm can produce the consistent boundary which improve the accuracy of the pressure field. Thus, improving our computational efficiency. Then, we give the energy stability of the algorithms and give a detailed proofs. The key idea to establish the stability results of the rotation velocity correction algorithm is to transform the rotation term into a telescopic symmetric form by means of the Gauge–Uzawa formula. Finally, numerical experiments show that the rotation velocity correction projection algorithm is efficient to solve the thermally coupled magneto-hydrodynamic equations.