An interior penalty Galerkin method for the MHD equations in heterogeneous domains

The Maxwell equations in the magnetohydrodynamic (MHD) limit in heterogeneous domains composed of conducting and non-conducting regions are solved by using Lagrange finite elements and by enforcing continuities across interfaces using an Interior Penalty technique a la Baker [Finite element methods for elliptic equations using non-conforming elements, Math. Comp. 31 (137) (1977) 45-59]. The method is shown to be stable and convergent and is validated by convergence tests. It is used to compute Ohmic decay in various compact conducting domains and to simulate the kinematic dynamo action in two different geometries.

[1]  Xiaoge Wang,et al.  SPLIB: A Library of Iterative Methods for Sparse Linear Systems , 1997 .

[2]  Ilaria Perugia,et al.  The hp-local discontinuous Galerkin method for low-frequency time-harmonic Maxwell equations , 2003, Math. Comput..

[3]  A. Bossavit.  Computational Electromagnetism: Variational Formulations, Complementarity, Edge Elements , 1997 .

[4]  정원식.  Non-linear 편집시스템 , 1996 .

[5]  Emmanuel Dormy,et al.  An integro-differential formulation for magnetic induction in bounded domains: boundary element-finite volume method , 2004 .

[6]  G. A. Baker.  Finite element methods for elliptic equations using nonconforming elements , 1977 .

[7]  H. K. Mo Att,et al.  Magnetic field generation in electrically conducting fluids , 1978 .

[8]  M. Costabel.  A coercive bilinear form for Maxwell's equations , 1991 .

[9]  Peter Frick,et al.  Screw dynamo in a time-dependent pipe flow. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  A. Bossavit.  The computation of eddy-currents, in dimension 3, by using mixed finite elements and boundary elements in association , 1991 .

[11]  Christophe Geuzaine,et al.  Modeling of thin insulating layers with dual 3-D magnetodynamic formulations , 2003 .

[12]  Jean-Luc Guermond,et al.  Mixed finite element approximation of an MHD problem involving conducting and insulating regions: The 3D case , 2003 .

[13]  Jean-Luc Guermond,et al.  A new Finite Element Method for magneto-dynamical problems: two-dimensional results , 2003 .

[14]  Jun Zou,et al.  A non-linear, 3-D spherical α2 dynamo using a finite element method , 2001 .

[15]  Jean-Luc Guermond,et al.  Effects of conductivity jumps in the envelope of a kinematic dynamo flow , 2006 .

[16]  S. Rebay.  Efficient Unstructured Mesh Generation by Means of Delaunay Triangulation and Bowyer-Watson Algorithm , 1993 .

[17]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[18]  Peter D. Minev,et al.  Mixed Finite Element approximation of an MHD problem involving conducting and insulating regions: the 2D case , 2002 .

[19]  Peter Frick,et al.  Non-stationary screw flow in a toroidal channel: way to a laboratory dynamo experiment , 2002 .

[20]  J. Lions,et al.  Problèmes aux limites non homogènes et applications , 1968 .

[21]  Ana Alonso,et al.  A MATHEMATICAL JUSTIFICATION OF THE LOW-FREQUENCY HETEROGENEOUS TIME-HARMONIC MAXWELL EQUATIONS , 1999 .

[22]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[23]  A. Bossavit,et al.  Magnetostatic problems in multiply connected regions: some properties of the curl operator , 1988 .