LAUSR

In this paper, we consider a fully discrete modular grad-div stabilization algorithm for time-dependent thermally coupled magnetohydrodynamic (MHD) equations. The main idea of the proposed algorithm is to add an extra minimally intrusive module to penalize the divergence errors of velocity and improve the computational efficiency for increasing values of the Reynolds number and grad-div stabilization parameters. In addition, we provide the unconditional stability and optimal convergence analysis of this algorithm. Finally, several numerical experiments are performed and further indicated these advantages over the algorithm without grad-div stabilization.