Linear Algebra Methods in a Mixed Approximation of Magnetostatic Problems
Authors:
- I. PerugiaIPI. Perugia
- V. SimonciniVSV. Simoncini
- M. ArioliMAM. Arioli
- 1 January 1999
- journal article
- Published by Society for Industrial & Applied Mathematics (SIAM) in SIAM Journal on Scientific Computing
- Vol. 21 (3) , 1085-1101
Abstract
In the solution of magnetostatic problems the use of a mixed formulation based on both the magnetic and magnetic displacement fields is particularly appropriate as it allows us to impose the physical conditions exactly and to maintain the continuity properties of the two fields, together with an efficient treatment of boundary conditions. The discretization by means of a proper finite element method yields a strongly structured algebraic linear system. This paper is concerned with the solution of this large, very sparse indefinite linear system. In particular, we present the implementation of some known direct and preconditioned iterative methods and discuss their performance on two-dimensional (2-D) and three-dimensional (3-D) specific models. We show that the 2-D system can be efficiently handled by appropriate variants of these schemes, while preliminary tests on the 3-D system give some insight in the understanding of the analysis that needs to be done.Keywords
This publication has 31 references indexed in Scilit:
- The Use ofQRFactorization in Sparse Quadratic Programming and Backward Error IssuesSIAM Journal on Matrix Analysis and Applications, 2000
- A mixed face-edge finite element formulation for 3D magnetostatic problemsIEEE Transactions on Magnetics, 1998
- Preconditioning in H(^cp) and applicationsMathematics of Computation, 1997
- A novel field-based mixed formulation of magnetostaticsIEEE Transactions on Magnetics, 1996
- Solution of augmented linear systems using orthogonal factorizationsBIT Numerical Mathematics, 1994
- On the use of the magnetic vector potential in the finite-element analysis of three-dimensional eddy currentsIEEE Transactions on Magnetics, 1989
- A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problemsMathematics of Computation, 1988
- A rationale for 'edge-elements' in 3-D fields computationsIEEE Transactions on Magnetics, 1988
- Robust pole-placement algorithm for adaptive controlIEE Proceedings D Control Theory and Applications, 1988
- New vector finite elements for three-dimensional magnetic field computationJournal of Applied Physics, 1987