Computation of Fields in Electromagnetic Systems with Hysteresis

A quasi-finite element method for calculation of static fields in media with hysteresis is presented in this paper. The proposed method utilizes vector Preisach hysteresis model to describe the material behavior. Time-stepping technique is employed in which the evolution of the electromagnetic system is traced to compute the field history as well as the current value of the field self-consistently. At every step in this technique a boundary value problem is solved by using coupled boundary-volume Galerkin's forms. In this formulation, the solution within the material medium is approximated in terms of finite element basis funtions, while the solution outside the material medium is approximated in terms of harmonic functions. An iterative technique for the solution of the resulting set of nonlinear equations is described. The proposed method is demonstrated on two dimensional electrostatic problems.