A Finite Element Method with Lagrange Multipliers for Low-Frequency Harmonic Maxwell Equations

The aim of this paper is to analyze a finite element method to solve the low-frequency harmonic Maxwell equations in a bounded domain containing conductors and dielectrics. This system of partial differential equations is a model for the so-called eddy currents problem. After writing this problem in terms of the magnetic field, it is discretized by Nedelec edge finite elements on a tetrahedral mesh. Error estimates are easily obtained if the curl-free condition is imposed on the elements in the dielectric domain. Then, the curl-free condition is imposed, at a discrete level, by introducing a piecewise linear multivalued potential. The resulting problem is shown to be a discrete version of other continuous formulation in which the magnetic field in the dielectric part of the domain has been replaced by a magnetic potential. Moreover, this approach leads to an important saving in computational effort. Problems related to the topology are also considered in that the possibility of having a nonsimply connected dielectric domain is taken into account. Implementation issues are discussed, including an amenable procedure to impose the boundary conditions by means of a Lagrange multiplier. Finally, the method is applied to solve a three-dimensional model problem: a cylindrical electrode surrounded by dielectric.

[1]  Habib Ammari,et al.  A Justification of Eddy Currents Model for the Maxwell Equations , 2000, SIAM J. Appl. Math..

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

[3]  Alberto Valli,et al.  Some remarks on the characterization of the space of tangential traces ofH(rot;Ω) and the construction of an extension operator , 1996 .

[4]  Rainer Picard,et al.  On the boundary value problems of electro- and magnetostatics , 1982, Proceedings of the Royal Society of Edinburgh: Section A Mathematics.

[5]  Peter Monk,et al.  A finite element method for approximating the time-harmonic Maxwell equations , 1992 .

[6]  Alberto Valli,et al.  An optimal domain decomposition preconditioner for low-frequency time-harmonic Maxwell equations , 1999, Math. Comput..

[7]  Peter P. Silvester,et al.  Finite elements for electrical engineers: Finite Elements for Electrical Engineers , 1996 .

[8]  A. Bossavit.  Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism , 1988 .

[9]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[10]  Lauri Kettunen,et al.  Formulation of the eddy current problem in multiply connected regions in terms of h , 1998 .

[11]  Alain Bossavit.  “Hybrid” electric–magnetic methods in eddy-current problems , 1999 .

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

[13]  Alfredo Bermúdez,et al.  A finite element method for the thermoelectrical modelling of electrodes , 1998 .

[14]  R. C. MacCamy,et al.  A skin effect approximation for eddy current problems , 1985 .

[15]  J.-C. Verite,et al.  A mixed fem-biem method to solve 3-D eddy-current problems , 1982 .

[16]  Patrick Ciarlet,et al.  On traces for functional spaces related to Maxwell's equations Part II: Hodge decompositions on the boundary of Lipschitz polyhedra and applications , 2001 .

[17]  Peter P. Silvester,et al.  Finite elements for electrical engineers: Frontmatter , 1996 .

[18]  V. Girault,et al.  Vector potentials in three-dimensional non-smooth domains , 1998 .

[19]  Martin Costabel,et al.  Symmetric coupling of finite elements and boundary elements for a parabolic-elliptic interface problem , 1990 .

[20]  A. Bu.  On traces for functional spaces related to Maxwell's equations Part I: An integration by parts formula in Lipschitz polyhedra , 2001 .

[21]  Manil Suri,et al.  A time-dependent interface problem for two-dimensional eddy currents , 1987 .

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

[23]  Alberto Valli,et al.  A domain decomposition approach for heterogeneous time-harmonic Maxwell equations , 1997 .

[24]  A. Bermúdez,et al.  Numerical computation of the electromagnetic field in the electrodes of a three‐phase arc furnace , 1999 .

[25]  Vivette Girault,et al.  Finite Element Methods for Navier-Stokes Equations - Theory and Algorithms , 1986, Springer Series in Computational Mathematics.

[26]  Patrick Ciarlet,et al.  On traces for functional spaces related to Maxwell's equations Part I: An integration by parts formula in Lipschitz polyhedra , 2001 .

[27]  Peter Monk,et al.  A mixed method for approximating Maxwell's equations , 1991 .

[28]  Gianni Gilardi,et al.  Magnetostatic and Electrostatic Problems in Inhomogeneous Anisotropic Media with Irregular Boundary and Mixed Boundary Conditions , 1997 .

[29]  S. I. Hariharan,et al.  Integral equation procedures for eddy current problems , 1982 .

[30]  Branko D. Popović,et al.  Introductory Engineering Electromagnetics , 1971 .

[31]  Peter Monk,et al.  Analysis of a finite element method for Maxwell's equations , 1992 .