The Jiles-Atherton Hysteresis Model
Magnetic hysteresis originates from the alignment of electron spins into magnetic domains and the reorganization of such domains under the influence of an externally applied magnetizing field. In the Jiles‑Atherton model, for each domain there is a reversible rotation and an irreversible rotation, which corresponds to a reversible and irreversible magnetization. In a magnetic material without hysteresis losses, there is only a reversible part of the magnetization, usually modeled with a relative permeability or with a reversible BH curve. The Jiles-Atherton hysteresis model includes a nonlinear relationship between the change in magnetization M and the change in the magnetizing field H. The magnetic flux density B follows from the relation.
The Jiles-Atherton model is a phenomenological model based on several parameters. The version available in COMSOL Multiphysics is full vectorial and anisotropic [4, 5], based upon the following five parameters:
•
Interdomain coupling: α
•
Maximum (saturation) magnetization: Ms
•
Langevin slope: a
•
Reversibility: c
•
Pinning: k
These quantities are second-order tensors, something that is not reflected by the equation notation here and in the graphical user interface.
The interdomain coupling parameter is used for defining the effective magnetic field He, which is the field that is assumed to affect the magnetic moments on a microscopical level.
The maximum magnetization and Langevin slope determines the anhysteretic part of the material behavior, which is supposed to follow the Langevin function
through the formula for anhysteretic magnetization Man.
Ms and α enter the expression with the nonlinear Langevin function and consistent with [4], they are available as diagonal matrices. The use of the Langevin function should be interpreted as applied component wise.
The reversibility parameter c, affects the degree of anhysteretic versus hysteretic behavior. If it is set to unity, the model is purely anhysteretic.
The remaining parameter, k, is controlling the pinning and hysteresis of the magnetic moments.
Given the parameters above, the key equation in the Jiles-Atherton model, determining the magnetization M, is
where the auxiliary vector χ is defined as:
The Jiles‑Atherton hysteresis model is available mainly for time-dependent studies, but it also supports stationary simulation. The latter is especially useful if the material is starting from a magnetized configuration that is not the solution of a previous transient study. In the graphical user interface, it is possible to select the initial value of magnetization, in which case the initial magnetic field is set accordingly.
It is implemented for the following physics interfaces: Magnetic Fields (based on vector potential); Magnetic Fields, No Currents (based on scalar potential); and Rotating Machinery, Magnetic (mixed vector and scalar potentials).
The internal state variables are discretized directly in the Gauss points, and an auxiliary variable is added to the same domain also for aiding support to simulations with moving domains.
The order of Gauss point discretization can be selected in the Discretization section under the Magnetic Field section, which is shown when Show>Discretization is selected in the Model Builder. For the Gauss point discretization order the default is Automatic. This default sets the Gauss points discretization order depending on the order of the parent physics. If the parent physics variable has order n, the Gauss point order is 2*(n–1). It is good practice to test that the results do not change significantly when using a Gauss point discretization of order 2*n since a good discretization order depends on the kind of magnetic potential, mesh element shape, and space dimension, rather than the pure mesh size. For an example, see.
Vector Hysteresis Modeling: Application Library path ACDC_Module/Other_Industrial_Applications/vector_hysteresis_modeling