Mixed finite element methods for stationary incompressible magneto-hydrodynamics

Summary.: A new mixed variational formulation of the equations of stationary incompressible magneto-hydrodynamics is introduced and analyzed. The formulation is based on curl-conforming Sobolev spaces for the magnetic variables and is shown to be well-posed in (possibly non-convex) Lipschitz polyhedra. A finite element approximation is proposed where the hydrodynamic unknowns are discretized by standard inf-sup stable velocity-pressure space pairs and the magnetic ones by a mixed approach using Nédélec's elements of the first kind. An error analysis is carried out that shows that the proposed finite element approximation leads to quasi-optimal error bounds in the mesh-siz

Finite element eigenvalue enclosures for the Maxwell operator

We propose employing the extension of the Lehmann-Maehly-Goerisch method developed by Zimmermann and Mertins, as a highly effective tool for the pollution-free finite element computation of the eigenfrequencies of the resonant cavity problem on a bounded region. This method gives complementary bounds for the eigenfrequencies which are adjacent to a given real parameter. We present a concrete numerical scheme which provides certified enclosures in a suitable asymptotic regime. We illustrate the applicability of this scheme by means of some numerical experiments on benchmark data using Lagrange elements and unstructured meshes.Comment: arXiv admin note: substantial text overlap with arXiv:1306.535

Error Analysis of a Fully Discrete Projection Method for Magnetohydrodynamic System

In this paper, we develop and analyze a finite element projection method for magnetohydrodynamics equations in Lipschitz domain. A fully discrete scheme based on Euler semi-implicit method is proposed, in which continuous elements are used to approximate the Navier–Stokes equations and H(curl) conforming Nédélec edge elements are used to approximate the magnetic equation. One key point of the projection method is to be compatible with two different spaces for calculating velocity, which leads one to obtain the pressure by solving a Poisson equation. The results show that the proposed projection scheme meets a discrete energy stability. In addition, with the help of a proper regularity hypothesis for the exact solution, this paper provides a rigorous optimal error analysis of velocity, pressure and magnetic induction. Finally, several numerical examples are performed to demonstrate both accuracy and efficiency of our proposed scheme

A nodal-based finite element approximation of the Maxwell problem suitable for singular solutions

A new mixed finite element approximation of Maxwell’s problem is proposed, its main features being that it is based on a novel augmented formulation of the continuous problem and the introduction of a mesh dependent stabilizing term, which yields a very weak control on the divergence of the unknown. The method is shown to be stable and convergent in the natural H(curl; ) norm for this unknown. In particular, convergence also applies to singular solutions, for which classical nodal based interpolations are known to suffer from spurious convergence upon mesh refinement.Postprint (published version