Stability and convergence analysis of A linear, fully decoupled and unconditionally energy stable scheme for magneto-hydrodynamic equations
This work provides a provably stable and efficient numerical method for solving coupled Navier-Stokes and Maxwell equations, which is important for computational fluid dynamics and plasma physics.
The authors developed a linear, fully decoupled, first-order time-stepping scheme for magneto-hydrodynamic equations that is unconditionally energy stable and proved optimal error estimates. Numerical experiments demonstrated stability and accuracy in simulating Kelvin-Helmholtz instability and magnetic-frozen phenomenon.
In this paper, we consider numerical approximations for solving the nonlinear magneto-hydrodynamical system, that couples the Navier-Stokes equations and Maxwell equations together. A challenging issue to solve this model numerically is about the time marching problem, i.e., how to develop suitable temporal discretizations for the nonlinear terms in order to preserve the energy stability at the discrete level. We solve this issue in this paper by developing a linear, fully decoupled, first order time-stepping scheme, by combining the projection method and some subtle implicit-explicit treatments for nonlinear coupling terms. We further prove that the scheme is unconditional energy stable and derive the optimal error estimates rigorously. Various numerical experiments are implemented to demonstrate the stability and the accuracy in simulating some benchmark simulations, including the Kelvin-Helmholtz shear instability and the magnetic-frozen phenomenon in the lid-driven cavity.