| Popis: |
This paper develops a weak Galerkin (WG) finite element method of arbitrary order for the steady thermally coupled incompressible Magnetohydrodynamics (MHD) flow. The WG scheme uses piecewise polynomials of degrees $k(k\geq 1),k,k-1,k-1$ and $k$ respectively for the approximations of the velocity, the magnetic field, the pressure, the magnetic pseudo-pressure, and the temperature in the interior of elements, and uses piecewise polynomials of degrees $l(l =k-1,k),k,k,k$ and $l$ for the numerical traces of velocity, the magnetic fields, the pressure, the magnetic pseudo-pressure, and the temperature on the interfaces of elements. The method is shown to yield globally divergence-free approximations of the velocity and magnetic fields. We give existence and uniqueness results for the discrete scheme and derive optimal a priori error estimates. We also present a convergent linearized iterative algorithm. Numerical experiments are provided to verify the obtained theoretical results. |
