Technical Reports

The ICICS/CS Reading Room

UBC CS TR-2010-07 Summary

A Mixed Finite Element Method with Exactly Divergence-free Velocities for Incompressible Magnetohydrodynamics, May 21, 2010 Chen Greif, Dan Li, Dominik Schoetzau and Xiaoxi Wei, 33 pages

We introduce and analyze a mixed finite element method for the numerical discretization of a stationary incompressible magnetohydrodynamics problem, in two and three dimensions. The velocity field is discretized using divergence-conforming Brezzi-Douglas-Marini (BDM) elements and the magnetic field is approximated by curl-conforming N\'{e}d\'{e}lec elements. The $H1$-continuity of the velocity field is enforced by a DG approach. A central feature of the method is that it produces exactly divergence-free velocity approximations, and captures the strongest magnetic singularities. We prove that the energy norm error is convergent in the mesh size in general Lipschitz polyhedra under minimal regularity assumptions, and derive nearly optimal a-priori error estimates for the two-dimensional case. We present a comprehensive set of numerical experiments, which indicate optimal convergence of the proposed method for two-dimensional as well as three-dimensional problems.

If you have any questions or comments regarding this page please send mail to