Divergence-Free $\boldsymbol{\mathcal{H}}(\mathbf{div})$-Conforming Hierarchical Bases for Magnetohydrodynamics (MHD)

Abstract

In order to solve the magnetohydrodynamics (MHD) equations with a $\boldsymbol{\mathcal{H}}(\mathbf{div})$-conforming element, a novel approach is proposed to ensure the exact divergence-free condition on the magnetic field. The idea is to add on each element an extra interior bubble function from a higher order hierarchical $\boldsymbol{\mathcal{H}}(\mathbf{div})$-conforming basis. Four such hierarchical bases for the $\boldsymbol{\mathcal{H}} (\mathbf{div})$-conforming quadrilateral, triangular, hexahedral, and tetrahedral elements are either proposed (in the case of tetrahedral) or reviewed. Numerical results have been presented to show the linear independence of the basis functions for the two simplicial elements. Good matrix conditioning has been confirmed numerically up to the fourth order for the triangular element and up to the third order for the tetrahedral element.

Keywords

Hierarchical bases / $\boldsymbol{\mathcal{H}}(\mathbf{div})$-conforming elements')">$\boldsymbol{\mathcal{H}}(\mathbf{div})$-conforming elements / Divergence-free condition

References

[1] Brackbill J.U., Barnes D.C. The effect of nonzero product of magnetic gradient and B on the numerical solution of the magnetohydrodynamic equations. J. Comput. Phys.. 1980, 35 426-430

[2] Tóth G. The ∇⋅B=0 constraint in shock-capturing magnetohydrodynamics codes. J. Comput. Phys.. 2000, 161 605-652

[3] Evans C.R., Hawley J.F. Simulation of magnetohydrodynamic flows: a constrained transport method. Astrophys. J.. 1988, 332 659-677

[4] Monk P. Finite Element Methods for Maxwell Equations. 2003 Oxford: Oxford University Press

[5] Nédélec J.C. Mixed finite elements in $\mathbb{R}^{3}$. Numer. Math.. 1980, 35 315-341

[6] Nédélec J.C. A new family of mixed finite elements in $\mathbb{R}^{3}$. Numer. Math.. 1986, 50 57-81

[7] Ainsworth M., Coyle J. Hierarchic finite element bases on unstructured tetrahedral meshes. Int. J. Numer. Methods Eng.. 2003, 58 2103-2130

[8] Falk R., Gatto P., Monk P. Hexahedral H (div) and H (curl) finite elements. Modél. Math. Anal. Numér.. 2011, 45 1 115-143

[9] Kwak D., Pyo H.C. Mixed finite element methods for general quadrilateral grids. Appl. Math. Comput.. 2011, 217 14 6556-6565

[10] Wheeler M., Xue G., Yotov I. A multipoint flux mixed finite element method on distorted quadrilaterals and hexahedra. Numer. Math.. 2012, 121 1 165-204

[11] Zaglmayr, S.: High order finite element methods for electromagnetic field computation. Ph.D. Dissertation, Johannes Kepler Universität, Linz (2006)

[12] Xin J., Cai W., Guo N. On the construction of well-conditioned hierarchical bases for $\boldsymbol{\mathcal{H}}(\mathbf{div})$-conforming $\mathbb{R}^{n}$ simplicial elements. Commun. Comput. Phys.. 2013, 14 3 621-638

[13] Magnus W., Oberhettinger F., Soni R.P. Formulas and Theorems for the Special Functions of Mathematical Physics. 1966 3 New York: Springer

[14] Cockburn B., Kanschat G., Schötzau D. A note on discontinuous Galerkin divergence-free solutions for the Navier–Stokes equations. J. Sci. Comput.. 2007, 31 61-73

[15] Greif C., Li D., Schötzau D., Wei X. A mixed finite element method with exactly divergence-free velocities for incompressible magnetohydrodynamics. Comput. Methods Appl. Mech. Eng.. 2010, 199 2840-2855