We present a steepest descent energy minimization scheme for micromagnetics. The method searches on a curve that lies on the sphere which keeps the magnitude of the magnetization vector constant. The step size is selected according to a modified Barzilai-Borwein method. Standard linear tetrahedral finite elements are used for space discretization. For the computation of static hysteresis loops the steepest descent minimizer is faster than a Landau-Lifshitz micromagnetic solver by more than a factor of two. The speed up on a graphic processor is 4.8 as compared to the fastest single-core CPU implementation.

The development of permanent magnets containing less or no rareearth elements is linked to profound knowledge of the coercivity mechanism. Prerequisites for a promising permanent magnet material are a high spontaneous magnetization and a sufficiently high magnetic anisotropy. In addition to the intrinsic magnetic properties the microstructure of the magnet plays a significant role in establishing coercivity. The influence of the microstructure on coercivity, remanence, and energy density product can be understood by using micromagnetic simulations. With advances in computer hardware and numerical methods, hysteresis curves of magnets can be computed quickly so that the simulations can readily provide guidance for the development of permanent magnets. The potential of rare-earth reduced and free permanent magnets is investigated using micromagnetic simulations. The results show excellent hard magnetic properties can be achieved in grain boundary engineered NdFeB, rare-earth magnets with a ThMn 12 structure, Co-based nano-wires, and L1 0 -FeNi provided that the magnet's microstructure is optimized.

International audienceThe coercive field and angular dependence of the coercive field of single-grain Nd$_{2}$Fe$_{14}$B permanent magnets are computed using finite element micromagnetics. It is shown that the thickness of surface defects plays a critical role in determining the reversal process. For small defect thicknesses reversal is heavily driven by nucleation, whereas with increasing defect thickness domain wall de-pinning becomes more important. This change results in an observable shift between two well-known behavioral models. A similar trend is observed in experimental measurements of bulk samples, where a Nd-Cu infiltration process has been used to enhance coercivity by modifying the grain boundaries. When account is taken of the imperfect grain alignment of real magnets, the single-grain computed results appears to closely match experimental behaviour

The coercive field of permanent magnets decreases with increasing grain size. The grain size dependence of coercivity is explained by a size dependent demagnetizing factor. In Dy free Nd$_2$Fe$_{14}$B magnets the size dependent demagnetizing factor ranges from 0.2 for a grain size of 55 nm to 1.22 for a grain size of 8300 nm. The comparison of experimental data with micromagnetic simulations suggests that the grain size dependence of the coercive field in hard magnets is due to the non-uniform magnetostatic field in polyhedral grains

A direct integration algorithm is described to compute the magnetostatic field and energy for given magnetization distributions on not necessarily uniform tensor grids. We use an analytically-based tensor approximation approach for function-related tensors, which reduces calculations to multilinear algebra operations. The algorithm scales with N4/3 for N computational cells used and with N2/3 (sublinear) when magnetization is given in canonical tensor format. In the final section we confirm our theoretical results concerning computing times and accuracy by means of numerical examples.