Micromagnetic simulations with periodic strayfield calculation of soft magnetic composite-materials

TL;DR

This paper introduces a numerical framework to simulate and analyze soft magnetic composite materials, predicting their magnetic properties and energy losses.

Contribution

A novel numeric framework combining two methods to decouple and analyze distinct loss contributions in soft magnetic materials.

Findings

Simulated trends align well with experimentally measured improvements in magnetic properties.

Distinct loss contributions are analyzed, offering insights into material behavior under different conditions.

Numerical simulations are validated against measurements of iron powders.

Abstract

This work provides a complete numeric framework, with which soft magnetic composite materials can be numerically analyzed through a combination of two distinct numerical methods. Key soft magnetic properties connected to energy efficiency, namely the permeability and the energy-loss can be predicted. The latter consists of different contributions. These contributions are treated separately in our work and therefore 2 decoupled numerical models are required to fully encompass these materials. Experimentally established concepts with the aim to improve magnetic properties for soft magnetic materials have been simulated, resulting in trends that agree well with the experimentally measured improvements. Moreover, we analyze the distinct loss contributions-an uncommon practice in experimental physics-providing valuable insights into the working principles of these materials, and therefore understanding which contributions become dominant under which conditions. Furthermore measurements of iron powders were conducted, in order to compare the numerical simulations with experimental data.

Figures

The comparison between the measured and the simulated powerloss for a Carbonyl iron powder^[35](#CR35)^ is shown as a function of the magnetic flux density *B*.

A schematic of the simulation cell and how it relates to the simulated lattice and the magnetic core (a) A toroidal magnetic core, as it would be produced from SMCs. (b) Zooming into the core, reveals the, oversimplified model, that the iron based particles form a perfect fcc-lattice, where the particles are not exchanged coupled but have an interparticle layer between them (c) The unit-cell of a fcc-structure, as it is utilized in the micromagnetic simulation with PBCs.Fig. 2The geometry that was chosen for the Maxwell-solver model is the primitive cell of the fcc-structure repeated in a $\do

Two different implementations for the grain-distribution (a) Regular cuboid grains with equal volumes (b) Voronoi grains with unequal volumes.

Two different implementations for the grain-distribution (a) Regular cuboid grains with equal volumes (b) Voronoi grains with unequal volumes.

A hysteresis loop of a soft magnetic composite simulation is shown, along with the magnetization images at different field magnitudes.Fig. 5The results of a Maxwell-solver simulation, with only the center fcc unit cell cut out. (a) Magnetic flux density distribution $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{B}$$\end{document}$ for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysy

The magnetic powerloss $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document}$ , and the permeability $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document}$ of the micromagnetic SMC-simulation models util

The magnetic powerloss $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document}$ , and the permeability $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document}$ of the micromagnetic SMC-simulation models util

The comparison between the measured and the simulated powerloss for a Carbonyl iron powder^[35](#CR35)^ is shown as a function of the magnetic flux density *B*.

Full text

Introduction

2022 showed a worldwide increase in energy consumption by 1.11%^1^, if compared to the prior year, this yields a total energy usage of 178.389 TWh, which marks a record in yearly global energy consumption. 2022 is not an outlier in this regard, with the exception of 2019 and 2009, since 1990 the trend of the global power consumption is growing with each year. As one way to combat this the EU has committed itself to reduce the energy consumption by 32.5%^2^ until 2030, as compared to 2007.

Electrical energy has accounted for roughly 20%^3^ of the worldwide energy consumption in 2022 and is forecast to increase to 50% until 2050^3^. Nearly all means of generating electrical energy, whether they are based on fossil fuels or renewable energies are based on generators powered by turbines, with the only exception being solar power. Soft magnetic materials are vital to generators and essentially every electric device that is currently being used, therefore we believe strongly that a simulation model which can predict the performance of soft magnets would be of great value. This work is focused on modeling a specific kind of soft magnet, however our methods could be used to study and improve different kinds of soft magnets as well.

A subcategory within soft magnetic materials are soft magnetic composites(SMCs). The market of SMCs in 2019 was approximated at 19.23 billion USD^4^ and is estimated to grow by 2029 to a value of 80 billion USD^5^. These materials have emerged in the 1990s and consist of a ferromagnetic powder, with particle sizes which range from  1 to 100 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\upmu }\hbox {m}$$\end{document} , which is coated by an electrical insulating material. The advantage of using an insulated powder, as opposed to a bulk material, is the decreased eddy-current loss at medium and high frequencies. The obvious disadvantage is that the permeability decreases, due to the airgap between the ferromagnetic particles. Until now there have been numerous studies^6–8^, which attempt at modeling SMCs utilizing numeric methods. The work of Arzbacher et al.^6^ is focused on modeling the demagnetizing field through Maxwells equations based on tomography data of the SMCs microstructure, while Weir et al.^7^ establish an analytical model in order to predict the permeability of SMCs. Bruckner et al.^8^ come closest to the study at hand, since it is also calculating SMCs utilizing micromagnetics. However the study in question, is a proof of concept, for the calculation of the demagnetization fields with PBCs and does not attempt to realistically depict the shape, the size, the material parameters or results which are comparable with commercially available SMCs. Within this work we perform efficient micromagnetic simulations utilizing true periodic boundary conditions proposed by Bruckner et al.^8^, along with numerically solving Maxwells equations, in order to realistically model SMCs and to optimize the SMC microstructure for low energy losses and high permeability.

We define the relative magnetic permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mu _r = \mu /\mu _0 = \frac{1}{\mu _0}\frac{dB}{dH}, \end{aligned}$$\end{document}

with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document} being the magnetic permeability, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _0$$\end{document} the magnetic permeability in vacuum, B the magnetic flux density, and H the applied external magnetic field. In this work \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} is always evaluated at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H=0$$\end{document} . In the following \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} will be referred to as the permeability.

The energy loss of SMCs can be split up into different contributions, the majority of literature^9–12^ available agrees on the definition as follows

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} W(f)=W_{\text {hyst}}+W_{\text {eddy}}(f)+W_{\text {res}}(f), \end{aligned}$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_{\text {hyst}}$$\end{document} is the static hysteresis loss, which arises from irreversible domain-wall motion. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_{\text {eddy}}$$\end{document} is the classic eddy current loss, the cause of which are electrical resistance losses within the material caused by a magnetic field which alternates with the frequency f, one can further separate this into inter- and intra-particle losses, meaning losses which arise from eddy currents within the particle and losses attributed to defects in the coating. Lastly \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_{\text {res}}$$\end{document} is the residual loss, which is usually defined as all other losses of unclear origin.

In this work, we have chosen to examine SMCs utilizing computational micromagnetism and Maxwell solvers. The micromagnetic model can resolve the structure of magnetic domains and domain walls, this yields accurate static hysteresis computation, with which \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_{\text {hyst}}$$\end{document} can be predicted. Another feature of micromagnetics is its ability to resolve the time-dependent magnetic switching processes and domain-wall movements. This is crucial due to the fact that applied frequencies of the external field are in the high-frequency regime( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 1$$\end{document} MHz), where the residual loss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_{\text {res}}$$\end{document} becomes relevant. A part of the residual losses originate in the increase of the hysteresis, which results from the fact that the magnetization lags behind the applied magnetic field due to its high frequency, we will refer to this as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_{\text {dyn}}$$\end{document} . This means that if one were to calculate the time-dependent hysteresis with the chosen micromagnetic-code, part of the total power loss is given, by calculating the area of the hysteresis. We call this the micromagnetic power-loss and define it as follows

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} P_m&= f\cdot W_{m}=f\cdot (W_{\text {hyst}}+W_{\text {dyn}})\nonumber \\&= f\cdot \oint {B(H)dH} \le f \cdot (W_{\text {hyst}}+W_{\text {res}}). \end{aligned}$$\end{document}

From here on out we will refer to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_{m}$$\end{document} as micromagnetic power-loss. Although our micromagnetic model does not predict eddy current losses, they can be calculated using Maxwell solvers, which effectively solve Maxwells equations. These solvers enable us to accurately compute the electric currents induced by an alternating magnetic field, allowing us to determine the eddy-current loss with the following equation

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} P_e=\frac{1}{2\sigma }\int |\mathbf {J(\textbf{x})}|^{2}dA, \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma$$\end{document} is the conductivity and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{J}$$\end{document} is the current density.

Micromagnetic theory

This section aims to give a short overview over both the micromagnetic theory, which lies at the core of the simulation work at hand. Firstly the individual micromagnetic energy contributions will be shortly discussed, with emphasis on how the inclusion of periodic boundary conditions(PBCs) alter the strayfield-calculation, as was already outlined in Ref.^8^. PBCs are essential in order to accurately simulate SMCs, due to the fact that the magnetic samples are very large, typically in the millimeter-regime, however only a small part of the structure can be simulated, since the application range of micromagnetism is in the nano- to micrometer regime. Simply truncating the magnetic domain would lead to a strong shape anisotropy and therefore to a much smaller permeability, and since the shape anisotropy is scale invariant, the actual size of the truncated geometry does not matter. For the same reason also pseudo periodic boundary conditions, as commonly used by many micromagnetic codes, will not produce correct results. All micromagnetic simulations in this work will be performed utilizing a GPU-accelerated finite difference code called \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {magnum.np}$$\end{document} ^13^, which will dramatically reduce the computation time.

The energy of a ferromagnet highly depends on the external field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{H}_{\text {zee}}$$\end{document} , therefore the simulations will employ the Zeeman-Energy. A SMC which is pressed into a transformer core, as can be seen in Fig. 1a, is subject to an external field provided by the primary coil. Therefore the external field within the simulation will have the form

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \textbf{H}_{\text {zee}}=\textbf{H}_{\text {ext}} \cdot sin(\omega t + \phi ), \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {H_{\text {ext}}}=(H_{\text {ext}},0,0)$$\end{document} is the amplitude, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega =2\pi f$$\end{document} is the angular frequency and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi$$\end{document} is the phase-difference. Since this is a local energy contribution, the calculation and numerical implementation doesn’t change with PBCs, when compared to open-boundary conditions.

Since iron has a cubic lattice symmetry, with a body-centered cubic structure, and most SMCs are iron based powders the simulations will utilize cubic anisotropy. Similarly to the Zeeman energy, the anisotropy energy is a purely local energy contribution, therefore the computation remains unaltered.

Every simulation cell is exchange coupled to its nearest neighboring cells through the exchange field, therefore the calculation changes slightly when PBCs are taken into account. Periodic boundary condition take the initial simulation box and replicate it infinitely in all directions, this means that the simulation cell on an edge surface is exchange coupled with the corresponding cell on the opposing edge surface, when evaluating discretized differential operators.

The demagnetizing field is notoriously the most cost-heavy energy-contribution to calculate, due to the fact that it is a long range interaction. Usually one solves

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \Delta u = \nabla \cdot \textbf{M}, \end{aligned}$$\end{document}

with u being the scalar potential of the demagnetizing field and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{M}$$\end{document} the magnetization. A standard method to solve (6) is to use the FFT and to apply zero padding in order to solve for the open boundary problem^14^ and for the scalar potential^15^. However since the demagnetizing energy is a long range interaction, we again have to alter the computation to fit the PBCs. Some standard micromagnetic codes use pseudo PBCs where the summation is truncated after a finite number of periodic images^16^. We will employ the methods of Ref.^8^ in order to calculate the demagnetizing field with true PBCs.

Maxwell solver theory

To calculate the eddy-current-losses we have to utilize a different model than our micromagnetic framework. In order to calculate the current density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{J}$$\end{document} we are going to solve Maxwells-equations numerically. However we will not discuss here in detail how we are bringing the equations into the weak formulation to solve them utilizing finite-elements. We will focus here on explaining the simulation model, and how we achieve periodic boundary conditions with the Maxwell–Solver. We have to utilize a different method for PBCs, because we were not able to unify the Fourier-based^8^ approach with the vectorfield-formulation. Therefore we want to implement periodic boundary conditions by using symmetry boundary conditions inspired by Ref.^17^. Therein a matrix of magnetic cubes is calculated, but instead of modeling all of the cubes, only one 1/8-cube is considered with symmetry conditions. Meaning that effectively the calculated space dependent, magnetic flux density of the 1/8-cube should be equivalent to the result for that part of a whole cube which is surrounded on all sides by cubes of the same size. The symmetry boundary conditions are

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \begin{aligned} &\textbf{B}\cdot \textbf{n}=0 \text { : On the 4 side surfaces} \\&\textbf{B}\times \textbf{n}=0 \text { : On the top \& bottom surfaces}, \end{aligned} $$\end{document}

with n being the orientation vector of the surface in question. The boundary conditions in (7) mean that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{B}$$\end{document} has to be tangential to the side planes and normal to the bottom and top plan. The simplest way to implement this is to utilize the reduced form of the Maxwell equations within Comsol.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \textbf{J}&=\nabla \times \textbf{H}, \\ \textbf{B}&=\nabla \times (\textbf{A}_b+\textbf{A}_r), \\ \textbf{J}&=\sigma \textbf{E}+j\omega \textbf{D}+\textbf{J}_e, \\ \textbf{E}&=-j\omega (\textbf{A}_b+\textbf{A}_r), \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma$$\end{document} is the electric conductivity, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{D}$$\end{document} is the displacement current field, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{E}$$\end{document} is the electric field and the magnetic vector potential \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{A}=\textbf{A}_b+\textbf{A}_r$$\end{document} is separated into a background field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{A}_b$$\end{document} and the reduced field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{A}_r$$\end{document} . Due to this set of equations the symmetry conditions from (7) change to

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \begin{aligned} &\textbf{n}\times \mathbf {A_r}=0 \Rightarrow \textbf{B}\cdot \textbf{n} = 0\text { : On the 4 side surfaces}, \\&\textbf{B}\times \textbf{n}=0 \text { : On the top \& bottom surfaces}. \end{aligned} $$\end{document}

The obvious limit of these boundary conditions is that symmetry is required, i.e. in this model we would not be able to apply the magnetic field in a direction that isn’t parallel to a coordinate axes. However for the study at hand this limitation is of no consequence.

The simulation model

Geometry

There are a wide variety of shapes when it comes to the magnetic powder, e.g. spheres^18–22^ , flakes^23–26^ or ribbons^27^. However the vast majority of ferromagnetic powders employed in SMCs are spherical, therefore we will also employ spherical particles. The ferromagnetic spheres were arranged in an fcc-structure. The particles are all separated by an interparticle layer, which is non-magnetic and non-conductive, therefore the particles are decoupled with regard to the exchange interaction. This gives us an upper threshold for the filling-factor of 74%, with the filling factor being the (volume-)fraction of magnetic material. The actual filling-factor in the simulation will be slightly reduced, due to the fact that we have to introduce a small gap between the individual particles. The filling factor used throughout all simulations in this work is 70%. Since we utilize PBCs we do not need to model an entire fcc structure, as is shown in Fig. 1b but can simply model a unit-cell with half and eighth spheres. This gives a cubic simulation box with six half-spheres centered at the 6 edge surfaces and eight eighth-spheres situated at each corner of the cube, which if repeated an infinite number of times, in all directions, gives a complete fcc structure. The micromagnetic model uses the simulation box as it is shown in Fig. 1c.

To ensure that the results of the eddy-current calculations are comparable with the micromagnetic model, it would be advantageous to simulate just the primitive cell of the fcc structure. However, this approach proved impractical because the simulation box in the eddy-current model cannot intersect magnetic material. If we cut through a region where current density flows-such as a half-sphere-the current density would form closed loops around the half-sphere. However, since the magnetic field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{H}$$\end{document} is always applied in the z-direction, the current paths are constrained to the xy-plane. Therefore, cutting through the xy-plane does not interfere with the current paths and is permissible.

The solution, shown in Fig. 2, involves repeating the primitive cell with whole spheres in a 3 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document} 3 grid, which is then sliced along the xy-plane. Additionally, a cube the size of the primitive cell is placed in the center of the model. This central cube does not represent any physical structure but is used for averaging purposes, to make the results comparable to the micromagnetic simulations. For example, when we report that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B = 100\,\text {mT}$$\end{document} in an eddy-current simulation, this means that the z-component of the magnetic flux density, averaged over the central cube, equals 100 mT.Fig. 1A schematic of the simulation cell and how it relates to the simulated lattice and the magnetic core (a) A toroidal magnetic core, as it would be produced from SMCs. (b) Zooming into the core, reveals the, oversimplified model, that the iron based particles form a perfect fcc-lattice, where the particles are not exchanged coupled but have an interparticle layer between them (c) The unit-cell of a fcc-structure, as it is utilized in the micromagnetic simulation with PBCs.Fig. 2. The geometry that was chosen for the Maxwell-solver model is the primitive cell of the fcc-structure repeated in a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3\times 3$$\end{document} grid with the spheres being cut along the xy-plane. A cube is added in the center of the mesh, this cube has no physical meaning, but is solely there so that one could average over a central primitive fcc-cell.

Grains and particles

Not only does the spacial distribution of the spheres matter, but the spheres also have a microstructure, which can be resolved by micromagnetics. The spherical iron particles consist of smaller grains which are exchange coupled, and have altering random easy axes in each grain. This only needs to be considered in the micromagnetic model, since Maxwells equations do not account for crystal anisotropy. In this work the axes will be chosen randomly, so that they are evenly distributed on the unit sphere, the three orthogonal easy axes are \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[\textbf{e}_i,\textbf{e}_j,\textbf{e}_k]$$\end{document} . The entire simulation box is then divided into \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\times n\times n$$\end{document} equally sized cubic grains, as can be seen in Fig. 3A. This gives the advantage that the dependence on grain-size is very easily controllable, since one can simply increase n in order to decrease the size of the grain, but the disadvantage is that it is not necessarily realistic. We will refer to a set of directions which is needed for one simulation as a seed of directions. Another possibility is to generate Voronoi cells as grains, as is shown in Fig. 3B.Fig. 3. Two different implementations for the grain-distribution (a) Regular cuboid grains with equal volumes (b) Voronoi grains with unequal volumes.

Simulation setup

In Tables 1 and 2 the parameters of the particle and the filler material are shown alongside the references for the respective values, for the two different simulation models. For the saturation polarization of the iron alloys we use \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _0 M_s=1.5\text {T}$$\end{document} , which is within the range of experimentally reported values^28^. Note that when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_s$$\end{document} is discussed, it is always supposed to mean the magnetic saturation of the magnetic material the spherical particles consist of. Since the micromagnetic simulation averages the magnetization over the entire simulation box, the hysteresis would not saturate at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_s$$\end{document} , but at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_{s,SMC}=\rho M_s$$\end{document} , where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document} is the filling factor. The anisotropy constants \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_2$$\end{document} chosen are those of pure iron^29^. SMCs come in a wide variety of available particle sizes, they approximately range from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1-150\,\upmu \text {m}$$\end{document} . The following simulations employs particles with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D=1\,\upmu \text {m}$$\end{document} , which places the simulated particles at the lower end of the spectrum. This is necessary in order for the simulation to be numerically feasible since the size of the particle one can simulate is limited by the memory of the GPU at hand. The size of the grains are varied in the study between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10-300\text { nm}$$\end{document} , if not stated explicitly the grain size is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=140\text { nm}$$\end{document} . Choosing d much larger would result in particles consisting of an unrealistically low number of grains, meaning that the seed of directions would influence the result significantly, however choosing the grainsize to small would mean that the random anisotropy model, which will be discussed later, would become relevant. Due to the cubic simulation box, finite difference methods are well suited, which are computationally less expensive compared to finite element methods. Hence, we use our in house developed finite difference code magnum.np^13^ for the micromagnetic model. One micromagnetic simulation consists of integrating the Landau-Lifschitz-Gilbert equation for a given geometry and parameters for a duration of three periods \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{3}{f}$$\end{document} . Obviously only one period is needed for the evaluation of the powerloss, however since the damping is relatively small, we calculate two periods prior to evaluating so that the system is in a reproducible state. The result is a hysteresis curve \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_x(H_x)$$\end{document} , from which one can evaluate the micromagnetic powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} and the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} . Figure 4 shows a standard hysteresis for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=1\text { MHz}$$\end{document} and material parameters as defined in table 1.

An eddy-current simulation involves solving Maxwell’s equations, as shown in Eq. (8), using the symmetry boundary conditions from Eq. (7) in the frequency domain. We used Comsol to carry out these simulations. Figure 5 presents the outcome of a particular simulation with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=1\text { MHz}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B=100\text { mT}$$\end{document} and the material parameters as shown in Table 2. Specifically, Fig. 5a,b display the norm of the magnetic flux density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\textbf{B}|$$\end{document} , and the norm of the electric current density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\textbf{J}|$$\end{document} as a function of space, respectively. The region of the simulation box shown in the figures represents only the parts which were evaluated to calculate the power loss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_e$$\end{document} .Fig. 4A hysteresis loop of a soft magnetic composite simulation is shown, along with the magnetization images at different field magnitudes.Fig. 5. The results of a Maxwell-solver simulation, with only the center fcc unit cell cut out. (a) Magnetic flux density distribution \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{B}$$\end{document} for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=10\text { MHz}$$\end{document} & \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B=100\text { mT}$$\end{document} . (b) Electric current density distribution \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textbf{J}$$\end{document} for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=10\text { MHz}$$\end{document} & \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B=100\text { mT}$$\end{document} .Table 1. Unless specifically stated otherwise, this table shows the parameters regarding the magnetic and geometric properties of the magnetic particle and the filler material which together forms the entire micromagnetic simulation box.ConstantSymbolValue magnetic particleValue fillerUnitSaturation magnetization \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_s$$\end{document} 1.5/ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _0$$\end{document} ^28^0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {A/m}$$\end{document} Cubical anisotropy 1st order \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} 48.0^29^0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {kJ/m}^{3}$$\end{document} Cubical anisotropy 2nd order \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_2$$\end{document} 0.5^29^0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {kJ/m}^{3}$$\end{document} Exchange couplingA10.00 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {pJ/m}$$\end{document} Alpha damping \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} 0.10 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {a.u.}$$\end{document} Particle sizeD1.0/ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu \text {m}$$\end{document} Grain sized140.0/ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {nm}$$\end{document} Filling factor \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document} 70.030.0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\%$$\end{document} Table 2. Unless specifically stated otherwise, this table shows the parameters regarding the magnetic and geometric properties of the magnetic particles and the filler material which together forms the entire Maxwell-solver simulation box.ConstantSymbolValue magnetic particleValue fillerUnitConductivity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma$$\end{document} 1e61 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {S/m}$$\end{document} Permittivity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon _r$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\infty$$\end{document} 1[a.u.]Magnetic bulk permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _b$$\end{document} 1e31 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {[a.u.]}$$\end{document} Particle sizeD1.0/ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu \text {m}$$\end{document} Filling factor \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho$$\end{document} 70.030.0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\%$$\end{document}

Results and discussion

The effect of random grains

The micromagnetic calculations at hand have a random contribution, namely the direction of the anisotropy within the grains, as well as the grain size and shape. Hence, we first show that our statistics is sufficiently good in order to obtain reliable results. This is only necessary for the micromagnetic simulations, since Maxwell equations do not account for anisotropy.

For this matter, we perform 100 simulations with different directional seed for each calculation, while keeping the total number of cubic grains constant. The frequency was chosen to be \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=100\text { MHz}$$\end{document} , since we care only about the statistics, we increased the frequency, which results in significantly reduced simulation time. Figure 6a along with Fig. 6b shows the results of the aforementioned study. The permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} seems to be completely independent with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r=7.39 \pm 0.06$$\end{document} , however the powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} does have a dependence on the directional seed with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m=(166.8 \pm 30.0)\text { kW/cm}^{3}$$\end{document} . Since a standard deviation of 18% is accurate enough for a first study, the following simulations are performed only with a single seed. However, if very small deviations are investigated the statistic has to be improved for more accurate results. A similar test was carried out with the Voronoi grains where we again performed 100 simulations where only the grain structure, was changed while the seed direction remained the same with each simulation, each time with a different seed of new centers which always gives a different Voronoi-structure, as is demonstrated in Fig. 6b. Again we observe that the permeability is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r=7.47 \pm 0.06$$\end{document} , meaning that the permeability is rather independent on the Voronoi seed, and also that it does not matter for the permeability wether regular or Voronoi grains are used. Examining the powerloss with Voronoi grains shows \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m=(130.1 \pm 19.9)\text { kW/cm}^{3}$$\end{document} , this means that different structures result in a smaller deviation than different directional seeds. However since both models show qualitatively similar results, in the following study the regular grain model is used. This has the additional advantage that they are easier to scale in terms of grainsize d.Fig. 6. The magnetic powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} , and the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} of the micromagnetic SMC-simulation models utilizing an fcc-particle-structure with differing (a) directional seeds in cubic grains. (b) Voronoi structures in Voronoi grains.Fig. 7. The magnetic powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} , the eddy-current powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_e$$\end{document} and the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} of the two decoupled SMC-simulation models utilizing an fcc-particle-structure and cubic grains for a sweep over (a) the magnetic flux density B at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=1\text { MHz}$$\end{document} . (b) the frequency f at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B=100\text { mT}$$\end{document} .

Sweep over magnetic flux density and frequency

The classic sweeps by which the magnetic losses of SMCs are characterized are the dependencies on the frequency f and the magnetic flux density B. This is represented by a phenomenological equation which is known as the Steinmetz equation

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} P_m=kB^{\alpha }f^{\beta }, \end{aligned}$$\end{document}

where k, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document} are the Steinmetz coefficients, which are material parameters, generally calculated from the hysteresis by fitting. It is well known that the eddy-current losses scale with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f^2$$\end{document} . Figure 7 shows the sweeps over magnetic flux density B and frequency f for the magnetic powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} , the eddy current powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_e$$\end{document} and the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} . \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} remains almost completely unchanged with increasing B or f, which is also what is to be expected according to literature^9,30^. The trends for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} show a linear increase in the log-log plot which means that they increase exponentially, through fitting we determined that the Steinmetz coefficient are \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =1.93$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta =1.71$$\end{document} , as is expected from the literature^31^. This gives validity to our model, since the simulations obeys Steinmetz equation, along with the Steinmetz coefficients being in the correct range.

From the log-log plot we can easily determine that the eddy-current losses increase with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f^2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^2$$\end{document} , moreover we can see \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_e$$\end{document} is multiple magnitudes smaller than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} . Therefore according to our model the magnetic losses, meaning the static hysteresis losses as well as the dynamic losses, which add on to the former, are significantly larger than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_e$$\end{document} . This is not entirely surprising since SMCs are designed toward having very low \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_e$$\end{document} . Therefore there is no reason to further optimize these materials to lower \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_e$$\end{document} , but only to optimize towards decreasing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} .

Sweep over anisotropy constant

A parameter which is crucial to the performance of SMCs is the crystal anisotropy constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} . However the influence of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} on the power-losses and the permeability has up to now hardly been investigated at all. This is is in part due to the fact that the measurement of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} is rather difficult for polycrystalline samples. This is where simulations as a tool can become helpful, since the anisotropy constant is one of the input parameters of the micromagnetic simulations. Figure 8a shows the losses \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} and the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} as function of the anisotropy constant \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} , for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B=100\text {mT}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f=1\text { MHz}$$\end{document} . It is evident that an increase in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} causes an increase in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} while it slightly decreases the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} . The reason \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} only slightly changes the permeability is, that the air gaps between the magnetic sphere dominate the effective permeability of the entire structure. The increased power-loss is intuitively clear since higher anisotropies result in more irreversible switching processes, which leads to increased hysteresis and therefore higher power-loss.Fig. 8(a) Sweep over the anisotropy constant Kc1 is shown, for the micromagnetic powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} and the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} . (b)Simulating the random anisotropy model, by changing the grain-size d for the micromagnetic powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} and the relative permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} , along with (c) the magnetization images at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B=100\text {mT}$$\end{document} .

Design concepts aimed toward improving magnetic performance

The random anisotropy model

Based on the results of the previous section one might conclude that to produce a better softmagnet one simply has to use a material with a small or vanishing crystal anisotropy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_1$$\end{document} . While this would certainly be a solution, an easier one might be to decrease the grainsize d, therefore utilizing the random anisotropy model. The random-anisotropy model was first proposed by Herzer^32^. To understand it one needs to understand that the critical scale at which the exchange energy balances the anisotropy energy is called the ferromagnetic exchange length \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_0$$\end{document} , which, in the case of uniaxial anisotropy, can be written as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} L_0=\sqrt{A/K_{u1}}, \end{aligned}$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_0$$\end{document} is the minimum length at which the magnetization can change in direction and therefore also determines the domainwall-width, A is the exchange constant and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_{u1}$$\end{document} is the uniaxial anisotropy constant. If one were to reduce the grainsize d below the exchange length \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_0$$\end{document} , the magnetization will not follow the direction of the individual grains, but instead is forced to align parallel by the exchange interaction. Therefore the effective anisotropy is averaged over several grains and thus considerably reduced in magnitude. The difference this makes for the coercive field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_c$$\end{document} and therefore also for the powerloss in the experiment can be observed in Ref.^32^. We have performed the same parameter-study with our simulations by reducing the grainsize in our model, as shown in Fig. 3A. The grains in our model are all equally sized cubes, meaning that decreasing d decreases the cubes until they get to the extreme where every single simulation cell is a different grain. The calculated results are illustrated in Fig. 8b, along with the magnetization configuration at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B=100\,mT$$\end{document} for different grainsizes in Fig. 8c. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} increases with increasing d, due to the fact that the effective anisotropy is decreased, until it plateaus at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\approx 70\text { nm}$$\end{document} , and then decreases again, reflecting the more conventional rule that good soft magnetic properties require large grains^33^. This behavior was already observed experimentally^32,34^. The permeability remains completely unaffected by the change in grainsize, which is not what one would expect from Ref.^32^. This discrepancy stems from the difference in hysteresis shapes between^32^ and our calculations. The SMCs we simulate can be characterized as airgap-dominated magnetic composites, since as already shown in Ref.^8^ the hysteresis depends heavily on the airgap, since it gets more and more sheared the bigger the airgap is or in other words, the smaller the fill-factor becomes. For small airgaps the hysteresis looks rectangular and the coercive field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_c$$\end{document} is directly related to the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r \propto H_c^{-1}$$\end{document} . This does not hold for hysteresis which are more sheared, as they are in our calculations due to the smaller fill-factor.

Reducing hysteresis loss under transverse magnetic field

Feng et al.^30^ have proposed a method with which the power loss of SMCs decreases at the cost of reducing the permeability to a lesser extent. It was found that if a transverse magnetic field is applied to the SMC core the power loss reduced by about 1/3 at the expense of less than 10% decrease in permeability. The effect is explained by a reduction of the hysteresis loss which is attributed to spin rotations and domain wall displacements which is completely locked in the direction of the applied field. Since this method effects the hysteresis loss, the proposed micromagnetic model should be able to replicate this effect.

To replicate the condition of the experiment we have repeated the simulations from Fig. 7a but this time we have applied a constant magnetic field in the z-direction with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\mu _0 \textbf{H}_{trans}|=175\text { mT}$$\end{document} while the AC-field is applied in the x-direction. Figure 9 shows the comparison between the losses of the performance of the model with and without a transversal field. The results show that the overall trend is the same in the simulation and in Ref.^30^. However in detail the trends differ, the simulations show a decrease in the powerloss by factor of roughly 4 while the permeability decreases by 20%. This discrepancy can be explained by the fact that in contrast to our simulations the experiments in Ref.^30^ were conducted in the kHz regime, therefore one could call them quasistatic and the losses given are simply the hysteresis losses, in contrast, we simulate the dynamic loss.Fig. 9(a)The powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} and (b) the permeability \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _r$$\end{document} are shown for different magnetic flux densities B with and without applying a transversal field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _0 H_{trans}=175\text { mT}$$\end{document} .

Comparison with experiments

In order to compare our findings to experimental results, we have performed measurements of Carbonyl iron powder^35^. SMC samples were prepared by mixing insulated soft magnetic carbonyl iron powder of grade SWS-S from BASF (median particle diameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_{50} = 3.8\, {\upmu }\text {m}$$\end{document} ) with an epoxy resin. The mixture was then molded into toroidal geometries under uniaxial compaction at a pressure below \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$30 \text { MPa}$$\end{document} . This low-pressure approach was chosen to avoid mechanical deformation of the iron particles and to avoid the introduction of internal stresses that could influence the magnetic performance. After molding, the samples were thermally cured and demolded without further post-processing. The magnetic properties of the toroidal samples were analyzed using an IWATSU SY-8218 B-H analyzer. All measurements were carried out under purely alternating current excitation at a fixed frequency of 1 MHz. To evaluate the field-dependent magnetic behavior and associated losses, magnetic flux density sweeps were performed. The use of toroidal ring cores ensured a closed magnetic path, minimizing demagnetizing effects and stray fields during measurement. For the purpose of comparison, this means that certain simulation parameters need to be altered from Tables 1 and 2 to that of pure iron, as can be seen in Table 3, all parameters not explicitly mentioned therein were taken from Tables 1 and 2. Furthermore we have seen that the directional seed influences the outcome of the micromagnetic simulation. Therefore the experimental results will be compared to micromagnetic simulations which have been carried out for ten different directional seeds and averaged afterwards, this is not necessary for the eddy-current simulation, since it doesn’t have a random contribution. Figure 10 shows the results of this comparison. All results have been averaged over the multiple measurements or simulations which were carried out. Since there are only few measuring points available, a fit was drawn through the data so that it is more easily comparable. We see that the simulated data is below the experimental data. This discrepancy decreases at higher field and is at its highest at lower fields. Adding the two simulated energy contribution would hardly change the result since the eddy-current losses are orders of magnitude below the micromagnetic losses, which is to be expected. One reason for this discrepancy could be an additional loss-source, which is not in the scope of this work. Anomalous eddy current losses are distinct from classical eddy current losses and are caused by the movement of magnetic domain walls within the material^36^. The losses would add onto our calculated loss and might explain the gap between the simulated and the experimental results. Calculating these losses would require to couple the micromagnetic model with the Maxwell model, which is outside the scope of this work. However this is not necessarily an issue. The simulations show good enough agreement and correctly predict the slope as a function of the magnetic flux density B.Table 3. This table shows the parameters used for the simulation of iron-powder, which were performed with the aim to compare the results to the measurements.SymbolValue magnetic particleUnit \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_s$$\end{document} 2.2/ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu _0$$\end{document} ^28^A/mA22.0pJ/m \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} 0.01a.u.D4.0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document} md140.0nm

Fig. 10. The comparison between the measured and the simulated powerloss for a Carbonyl iron powder^35^ is shown as a function of the magnetic flux density B.

Conclusion

We have presented a numerical model, which allows to simulate soft magnetic materials without finite size effects. This was accomplished through the implementation of periodic boundary conditions within the strayfield computation. All relevant energy contributions were demonstrated to be calculable using two distinct simulation models: the first based on micromagnetism and the second solving Maxwell’s equations. We have chosen to verify our results by reproducing experiments which aim to alter the powerloss or the permeability. The simulations obey Steinmetz-law when it comes to the B and f dependent powerloss \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_m$$\end{document} , with the Steinmetz-coefficients being \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =1.93$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta =1.71$$\end{document} . It was shown that the powerloss can be reduced by utilizing the random anisotropy model through decreasing the grainsize d. Alternatively the powerloss can be decreased by applying a transversal field \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{trans}$$\end{document} , while at the same time decreasing the permeability by a lesser degree. Additionally, we compared our simulations to experimental findings. While the results did not match perfectly, the fact that they are within the same order of magnitude is highly encouraging. Furthermore, we have demonstrated how power loss is distributed among various loss contributions-an analysis rarely conducted, particularly for frequencies in the MHz range. Notably, our findings reveal that micromagnetic losses are the dominant factor, highlighting the importance of focusing future development efforts on minimizing these losses. The financial support of Huawei is acknowledged.