Time-Harmonic Magnetic Model¶
Introduction¶
The time-harmonic magnetic model describes electromagnetic fields in terms of their complex amplitudes at a single angular frequency \(\omega = 2\pi f\), under the low-frequency (eddy-current) approximation where displacement currents are neglected. This approximation is applicable when wave-propagation effects are negligible compared to diffusion/induction effects, which is typical for many low-frequency magnetic devices.
The time-harmonic magnetic model solves
\begin{align}
\nabla \times \left( \nu \nabla \times \tilde{\mathbf{A}} \right)
+ j \omega \sigma \tilde{\mathbf{A}}
= \tilde{\mathbf{J}}_\mathrm{s},
\end{align}
where - \(\nu = 1/\mu\) is the magnetic reluctivity - \(\sigma\) is the electric conductivity - \(\tilde{\mathbf{J}}_\mathrm{s}\) is the impressed (source) current density.
In a conductive material, the induced (eddy) current density is given by \(\tilde{\mathbf{J}} = \sigma \tilde{\mathbf{E}}\).
Model¶
The model can be created and added to the simulation using
time_harmonic_magnetic_model = TimeHarmonicMagneticModel(
marker=magnetic_domain,
frequency=60,
order=2
)
sim.get_model_manager().add_model(time_harmonic_magnetic_model)
marker determines where the domain on which the model is solved
- frequency \(\left[ \rm{Hz} \right]\) the frequency of the fields
- order the model order
Materials¶
In the most general case, the material parameters \(\nu\) and \(\sigma\) in Eq. (3) can be complex, frequency-dependent tensors that vary spatially and may also depend on other physical properties, such as temperature. A material is specified through the Time-Harmonic General Magnetic Material. For details refer to Time-Harmonic General Magnetic Material.
Most materials can be setup using the General Material with an example provided in Table 1.
![]() Example from the Compumag TEAM-7 case. The simulation contains three materials: Air, Copper and Aluminum. |
Air For non-magnetic and non-conductive regions such as air with $$ \mu = \mu_0,\qquad \sigma = 0,\qquad $$ we can create the material using:
CopperFor non-magnetic but electrically conductive regions such as copper with $$ \mu = \mu_0,\qquad \sigma = \sigma_{\text{Cu}},\qquad $$ we can create the material using:
(note that we ignore eddy currents here due to the stranded coil).
AluminumFor electrically conductive regions such as aluminum we use $$ \mu = \mu_0,\qquad \sigma = \sigma_{\text{al}},\qquad $$ we can create the material using:
|
Note that once created the materials need to be added to the model
Conditions¶
Equation (3) is solved within a specific domain and must be accompanied by appropriate boundary conditions and excitations.
| Name | Supported Entities | Description |
|---|---|---|
| Tangential Magnetic Flux | Boundary | Enforces the magnetic flux density to be tangential to the boundary. |
| Normal Magnetic Field | Boundary | Enforces the magnetic field to be normal to the boundary. |
| Tangential Magnetic Field | Boundary | Sets the tangential component of the magnetic field on the boundary. |
Reports¶
Coefficients¶
The following functions are available in the time-harmonic magnetic model for visualization or querying:
| Name | Field Type | Description |
|---|---|---|
| Magnetic Vector Potential-Real | Vector | The real part of the complex amplitude of magnetic vector potential: $$ \mathrm{Re}\left\{\tilde{\mathbf{A}}\right\} $$ |
| Magnetic Vector Potential-Imag | Vector | The imaginary part of the complex amplitude of magnetic vector potential: $$ \mathrm{Im}\left\{\tilde{\mathbf{A}}\right\} $$ |
| Magnetic Flux Density-Real | Vector | The real part of the complex amplitude of magnetic flux density: $$ \mathrm{Re}\left\{\tilde{\mathbf{B}}\right\} = \mathrm{Re}\left\{\nabla \times \tilde{\mathbf{A}}\right\} $$ |
| Magnetic Flux Density-Imag | Vector | The imaginary part of the complex amplitude of magnetic flux density: $$ \mathrm{Im}\left\{\tilde{\mathbf{B}}\right\} = \mathrm{Im}\left\{\nabla \times \tilde{\mathbf{A}}\right\} $$ |
| Magnetic Field-Real | Vector | The real part of the magnetic field vector, calculated as $$ \mathrm{Re}\left\{\tilde{\mathbf{H}}\right\} = \mathrm{Re}\left\{\nu \, \tilde{\mathbf{B}}\right\}, $$ where, $\nu$ is the magnetic reluctivity, and $\tilde{\mathbf{B}}$ is the magnetic flux density vector. |
| Magnetic Field-Imag | Vector | The imaginary part of the magnetic field vector, calculated as $$ \mathrm{Im}\left\{\tilde{\mathbf{H}}\right\} = \mathrm{Im}\left\{\nu \, \tilde{\mathbf{B}}\right\}, $$ where, $\nu$ is the magnetic reluctivity, and $\tilde{\mathbf{B}}$ is the magnetic flux density vector. |
| Electric Current Density-Real | Vector | The real part of the electric current density: $$ \mathrm{Re}\left\{\tilde{\mathbf{J}}\right\} = -\omega \, \mathrm{Re}\left\{j \sigma \tilde{\mathbf{A}}\right\}, $$ $\omega$ is the angular frequency, $\sigma$ is the electric conductivity, and $\mathbf{A}$ is the magnetic vector potential. |
| Electric Current Density-Imag | Vector | The imaginary part of the electric current density: $$ \mathrm{Im}\left\{\tilde{\mathbf{J}}\right\} = -\omega \, \mathrm{Im}\left\{j \sigma \tilde{\mathbf{A}}\right\}, $$ $\omega$ is the angular frequency, $\sigma$ is the electric conductivity, and $\mathbf{A}$ is the magnetic vector potential. |
| Ohmic Heating | Scalar | Time-averaged ohmic loss density: $$ P_\Omega = \frac{1}{2} \mathrm{Re}\left\{ \tilde{\mathbf{J}} \cdot \tilde{\mathbf{E}}^* \right\}, $$ where $^*$ denotes the complex conjugate. |
| Relative Magnetic Permeability | Scalar | The relative magnetic permeability is the inverse of the magnetic reluctivity: $\mu = 1/\nu$. |
