Volume 13, Issue 1
Numerical Analysis and Testing of a Fully Discrete, Decoupled Penalty-Projection Algorithm for MHD in Els?sser Variable

M. Akbas ,  S. Kaya ,  M. M. Jaman and L. G. Rebholz

Int. J. Numer. Anal. Mod., 13 (2016), pp. 90-113

Preview Full PDF BiBTex 305 661
  • Abstract

We consider a fully discrete, effcient algorithm for magnetohydrodynamic (MHD) ow that is based on the Elsässer variable formulation and a timestepping scheme that decouples the MHD system but still provides unconditional stability with respect to the timestep. We prove stability and optimal convergence of the scheme, and also connect the scheme to one based on handling each decoupled system with a penalty-projection method. Numerical experiments are given which verify all predicted convergence rates of our analysis on some analytical test problems, show the results of the scheme on a set of channel flow problems match well the results found when the computation is done with MHD in primitive variable, and finally show the scheme performs well on a channel flow over a step.

  • History

Published online: 2016-01

  • AMS Subject Headings

65M60, 76W05

  • Cited by