Skip to main content Accessibility help
×
Home

High Order Hierarchical Divergence-Free Constrained Transport H(div) Finite Element Method for Magnetic Induction Equation

  • Wei Cai (a1), Jun Hu (a2) and Shangyou Zhang (a3)

Abstract

In this paper, we propose to use the interior functions of an hierarchical basis for high order BDMp elements to enforce the divergence-free condition of a magnetic field B approximated by the H(div) BDMp basis. The resulting constrained finite element method can be used to solve magnetic induction equation in MHD equations. The proposed procedure is based on the fact that the scalar (p–1)-th order polynomial space on each element can be decomposed as an orthogonal sum of the subspace defined by the divergence of the interior functions of the p-th order BDMp basis and the constant function. Therefore, the interior functions can be used to remove element-wise all higher order terms except the constant in the divergence error of the finite element solution of the B-field. The constant terms from each element can be then easily corrected using a first order H(div) basis globally. Numerical results for a 3-D magnetic induction equation show the effectiveness of the proposed method in enforcing divergence-free condition of the magnetic field.

Copyright

Corresponding author

*Corresponding author. Email addresses: wcai@uncc.edu (W. Cai), hujun@math.pku.edu.cn (J. Hu), szhang@udel.edu (S. Y. Zhang)

References

Hide All
[1] Ainsworth, M. and Coyle, J., Hierarchic finite element bases on unstructured tetrahedral meshes, Int. J. Numer. Methods Eng., 58 (2003), pp. 21032130.
[2] Cai, W., Wu, J. and Xin, J. G., Divergence-free H(div)-conforming hierarchical bases for magnetohydrodynamics (MHD), Commun. Math. Stat., 1 (2013), pp. 1935.
[3] Brackbill, J. U. and Barnes, D. C., The effect of nonzero product of magnetic gradient and B on the numerical solution of the magnetohydrodynamic equations, J. Comput. Phys., 35 (1980), pp. 426430.
[4] Brezzi, F., Douglas, J., Duran, R. and Fortin M, M., Mixed finite elements for second order elliptic problems in three variables, Numerische Mathematik, 51(2) (1987), pp. 237–50.
[5] Evans, C. R. and Hawley, J. F., Simulation of magnetohydrodynamic flows: a constrained transport method, Astrophys. J., 332 (1988), pp. 659677.
[6] Tóth, G., The ∇ · B = 0 constraint in shock-capturing magnetohydrodynamics codes, J. Comput. Phys., 161 (2000), pp. 605652.
[7] Hu, J., Finite element approximations of symmetric tensors on simplicial grids in Rn: the higher order case, J. Comput. Math., 33 (2015), pp. 283296.
[8] Hu, J., A new family of efficient conforming mixed finite elements on both rectangular and cuboid meshes for linear elasticity in the symmetric formulation, SIAM J. Numer. Anal., 53 (2015), pp. 14381463.
[9] Hu, J. and Zhang, S., A family of symmetric mixed finite elements for linear elasticity on tetredral grids, Sci. China Math., 58 (2015), pp. 297307.

Keywords

MSC classification

High Order Hierarchical Divergence-Free Constrained Transport H(div) Finite Element Method for Magnetic Induction Equation

  • Wei Cai (a1), Jun Hu (a2) and Shangyou Zhang (a3)

Metrics

Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed