Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
467861 | Computers & Mathematics with Applications | 2015 | 21 Pages |
In the paper we consider the penalty finite element method for the stationary incompressible magnetohydrodynamics (MHD) problem with a factor of penalty parameter. Stability and convergence of numerical solutions are established. Furthermore, two level penalty methods are also developed for the MHD problem. Our methods consist of solving a nonlinear MHD problem by the usual penalty method on a coarse mesh with mesh size HH, and then a linearized MHD problem based on the Stokes, Newton and Oseen iterations on a fine mesh is solved by the penalty method with mesh size h(h≪H), respectively. Stability and error estimates of numerical solutions in two level penalty methods are presented. Finally, some numerical tests are provided to demonstrate the effectiveness of the developed algorithms.