Mixed finite element approximation of incompressible MHD problems based on weighted regularization

Authors:
Urs Hasler;Anna Schneebeli;Dominik Schötzau
Affiliations:
Department of Mathematics, University of Basel, Rheinsprung 21, 4051 Basel, Switzerland;Department of Mathematics, University of Basel, Rheinsprung 21, 4051 Basel, Switzerland;Mathematics Department, University of British Columbia, 1984 Mathematics Road, Vancouver, BC V6T 1Z2, Canada
Venue:
Applied Numerical Mathematics
Year:
2004

Citing 5
Cited 2

Stationary stokes and Navier-Stokes systems on two-or three-dimensional domains with corners. Part I: linearized equations

SIAM Journal on Mathematical Analysis
Stability of higher-order Hood-Taylor methods

SIAM Journal on Numerical Analysis
Mixed and hybrid finite element methods

Mixed and hybrid finite element methods
Thehp-local discontinuous Galerkin method for low-frequency time-harmonic Maxwell equations

Mathematics of Computation
Mixed finite element methods for stationary incompressible magneto–hydrodynamics

Numerische Mathematik

Approximation of the thermally coupled MHD problem using a stabilized finite element method

Journal of Computational Physics
On an unconditionally convergent stabilized finite element approximation of resistive magnetohydrodynamics

Journal of Computational Physics

Quantified Score

Hi-index	0.01

Visualization

Abstract

We introduce and analyze a new mixed finite element method for the numerical approximation of stationary incompressible magneto-hydrodynamics (MHD) problems in polygonal and polyhedral domains. The method is based on standard inf-sup stable elements for the discretization of the hydrodynamic unknowns and on nodal elements for the discretization of the magnetic variables. In order to achieve convergence in non-convex domains, the magnetic bilinear form is suitably modified using the weighted regularization technique recently developed in [Numer. Math. 93 (2002) 239]. We first discuss the well-posedness of this approach and establish a novel existence and uniqueness result for non-linear MHD problems with small data. We then derive quasi-optimal error bounds for the proposed finite element method and show the convergence of the approximate solutions in non-convex domains. The theoretical results are confirmed in a series of numerical experiments for a linear two-dimensional Oseen-type MHD problem, demonstrating that weighted regularization is indispensable for the resolution of the strongest magnetic singularities.