A Finite Element Method for the Magnetostatic Problem in Terms of Scalar Potentials

The aim of this paper is to analyze a numerical method for solving the magnetostatic problem in a three-dimensional bounded domain containing prescribed currents and magnetic materials. The method discretizes a well-known formulation of this problem based on two scalar potentials: the total potential, defined in magnetic materials, and the reduced potential, defined in dielectric media and in nonmagnetic conductors carrying currents. The topology of the domain of each material is not assumed to be trivial. The resulting variational problem is proved to be well posed and is discretized by means of standard piecewise linear finite elements. Transmission conditions are imposed by means of a piecewise linear Lagrange multiplier on the surface separating the domains of both potentials. Error estimates for the numerical method are proved and the results of some numerical tests are reported to assess the performance of the method.

[1]  Ilaria Perugia,et al.  Vector potential formulation for Magnetostatics and modeling of permanent magnets , 1999 .

[2]  Stuart Crozier,et al.  Rapid computation of static fields produced by thick circular solenoids , 1997 .

[3]  A. Razek,et al.  Numerical modelling of 3D magnetostatic saturated structures with a hybrid FEM-BEM technique , 1992 .

[4]  C. Johnk,et al.  Engineering Electromagnetic Fields and Waves , 1975 .

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

[6]  A. Buffa,et al.  On traces for H(curl,Ω) in Lipschitz domains , 2002 .

[7]  A. Meir,et al.  Variational methods for stationary MHD flow under natural interface conditions , 1996 .

[8]  K. Preis,et al.  Comparison of different finite element formulations for 3D magnetostatic problems , 1988 .

[9]  Christian Magele,et al.  Boundary element analysis of 3-D magnetostatic problems using scalar potentials , 1992 .

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

[11]  Alfredo Bermúdez,et al.  A Finite Element Method with Lagrange Multipliers for Low-Frequency Harmonic Maxwell Equations , 2002, SIAM J. Numer. Anal..

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

[13]  Ilaria Perugia,et al.  A mixed formulation for 3D magnetostatic problems: theoretical analysis and face-edge finite element approximation , 1999, Numerische Mathematik.

[14]  K. Preis,et al.  On the use of the magnetic vector potential in the nodal and edge finite element analysis of 3D magnetostatic problems , 1996 .

[15]  Magnetostatics with scalar potentials in multiply connected regions (comments with reply) , 1989 .

[16]  K. Preis,et al.  Computation of 3-D Magnetostatic fields using a reduced scalar potential , 1992, Digest of the Fifth Biennial IEEE Conference on Electromagnetic Field Computation.

[17]  C. Magele,et al.  Numerical analysis of 3D magnetostatic fields , 1991 .

[18]  Peter Monk,et al.  Finite Element Methods for Maxwell's Equations , 2003 .

[19]  Ilaria Perugia,et al.  A field-based finite element method for magnetostatics derived from an error minimization approach , 2000 .

[20]  Fumio Kikuchi.  Mixed formulations for finite element analysis of magnetostatic and electrostatic problems , 1989 .

[21]  W. M. Rucker,et al.  Three-dimensional magnetostatic field calculation using boundary element method , 1988 .

[22]  L. R. Scott,et al.  Finite element interpolation of nonsmooth functions satisfying boundary conditions , 1990 .

[23]  Isaak D. Mayergoyz,et al.  A new scalar potential formulation for three-dimensional magnetostatic problems , 1987 .