An unconventional algorithm is presented to compute quasi-static magnetic fields. It aims to be as close to the physics as possible for the class of strongly heterogeneous media. Using edge-element expansion functions for the magnetic field strength and face-element expansion functions for the magnetic flux density, a system of linear algebraic equations in the expansion coefficients is constructed from the exact satisfaction, for each element, of the domain-integrated field equations and of the compatibility relations, combined with the least-square satisfaction, for each element, of the constitutive relation. The resulting system of equations is over-determined and is solved by minimizing the L/sup 2/-norm of the residual. Accuracy and convergence are tested by applying the method to some two-dimensional problems whose solution is known analytically. A novel method to encompass an unbounded exterior domain is included in the method; it performs remarkably well.
[1]
Richard Perry,et al.
McGraw-Hill Book Company
,
1973
.
[2]
A. Bossavit.
Whitney forms: a class of finite elements for three-dimensional computations in electromagnetism
,
1988
.
[3]
A. T. Hoop.
Handbook of radiation and scattering of waves
,
1995
.
[4]
Ilaria Perugia,et al.
A novel field-based mixed formulation of magnetostatics
,
1996
.
[5]
A. Bossavit.
Complementary formulations in steady-state eddy-current theory
,
1992
.
[6]
A. Bossavit,et al.
Constitutive inconsistency: rigorous solution of Maxwell equations based on a dual approach
,
1994
.
[7]
On Non-Linear Magnetostatics: Dual-Complementary Models and “Mixed” Numerical Methods
,
1990
.