Electromagnetics

Type:section
Appearance:simple
Excludes:ContinuumMechanics, HeatConduction

Electromagnetic field problems address the computation of the electromagnetic field. This means the following impressed tensor fields are computed: the electric field strength \VField{E}, the electric flux density \VField{D}, the magnetic field strength \VField{H}, and the magnetic flux density \VField{B}.

Macroscopically, these fields are related by Maxwell’s equations

\begin{eqnarray*}
\curl \VField{E} & = & - \partial_t \VField{B}, \\
\curl \VField{H} & = &  \partial_t \VField{D}+\VField{J}, \\
\divo \VField{B} & = & 0,  \\
\divo \VField{D} & = & \SField{\rho},
\end{eqnarray*}

where \VField{J} and \SField{\rho} are the macroscopically free current and charge distributions, respectively. Maxwell’s equations are supplemented by constitutive relations characterizing the macroscopic medium. JCMsolve supports dispersive and anisotropic materials given by a linear relation in the frequency domain:

\begin{alignat*}{3}
\VField{D}(\pvec{x}, \omega) \phantom{x} &=& \phantom{x}  \TField{\varepsilon}(\pvec{x}, \omega) \VField{E}(\pvec{x}, \omega) & + \tilde{\TField{\varepsilon}}(\pvec{x}, \omega) \VField{H}(\pvec{x}, \omega)  \\
\VField{B}(\pvec{x}, \omega) \phantom{x} &=& \phantom{x} \TField{\mu}(\pvec{x}, \omega) \VField{H}(\pvec{x}, \omega) & + \tilde{\TField{\mu}}(\pvec{x}, \omega) \VField{E}(\pvec{x}, \omega)
\end{alignat*}

Here, \TField{\varepsilon} and \TField{\mu} are the permittivity and the permeability tensors, respectively. The tensors \tilde{\TField{\varepsilon}} and \tilde{\TField{\mu}} are used to define general bi-anisotropic, chiral materials , see bi-permittivity and bi-permeability. [1]

The material is further characterized by the electric conductivity \sigma, which relates the electric field intensity to the Ohmic electric current:

\begin{eqnarray*}
\VField{J}^{\mathrm{(Ohm)}}(\pvec{x}, \omega) & = & \TField{\sigma}(\pvec{x}, \omega) \VField{E}(\pvec{x}, \omega) \,.
\end{eqnarray*}

Transparent boundary conditions

Boundary conditions are an important topic when simulating Maxwell’s equations. In some cases it is justified to prescribe the tangential components of the electric field intensity or of the magnetic field intensity. More generally, an impedance boundary condition relates the tangential components of the electric and magnetic fields, c.f., boundary condition section Electromagnetic.

For high-frequency Maxwell-equations the electromagnetic field radiates over a huge distance of wavelengths into the surrounding of the device without a significant damping. In this case, the computational domain only encloses the device and JCMsolve has to take special care to guarantee field fidelity at the boundaries. In this case, so-called transparent boundaries are to be imposed.

However, within the FEM context of JCMsolve it is more appropriate to think of this transparency not as a boundary condition. Instead, JCMsolve contains numerical methods for solving the problem on the infinite domain. For doing this, a complete description of the physics in the exterior domain has to be passed. This includes the material distribution as well as the incoming electromagnetic field, which illuminates the device. Besides the near field solution, JCMsolve provides post-processes to evaluate the field in the exterior domain and to compute the far field with high precision.

Periodic boundary conditions

Many optical devices such as gratings, photonic crystals or photomasks contain periodic repetitions in one or several directions of a structure filling a unit cell \Omega. To put it in mathematical terms, the material tensors exhibit a translation invariance on a lattice, i.e.:

\begin{eqnarray*}
\TField{\varepsilon}(\pvec{x}+\pvec{a}_i, \omega) & = & \TField{\varepsilon}(\pvec{x}, \omega) \\
\TField{\mu}(\pvec{x}+\pvec{a}_i, \omega) & = & \TField{\mu}(\pvec{x}, \omega),
\end{eqnarray*}

where \pvec{a}_i is a lattice vector. Surely, a real world device is not infinitely periodic. However, the periodic arrangement often spreads over a region with a dimension much larger than the wavelength, so that simulating an infinite periodic lattice is a good model.

To exploit the periodicity in a numerical simulation it is required that the impressed tensor fields exhibit a compatible symmetry condition. Especially, this concerns the source fields \VField{J}, \SField{\rho} and the illuminating electromagnetic field given by the incoming electric and magnetic fields \VField{E}^{(\inc)}, \VField{H}^{(\inc)}:

\begin{eqnarray*}
\VField{E}^{(\inc)}(\pvec{x}+\pvec{a}_i, \omega) & = &  \VField{E}^{(\inc)}(\pvec{x}, \omega)e^{i\pvec{k_{\mathrm{B}}} \pvec{a}} \\
\VField{H}^{(\inc)}(\pvec{x}+\pvec{a}_i, \omega) & = &  \VField{H}^{(\inc)}(\pvec{x}, \omega)e^{i\pvec{k_{\mathrm{B}}} \pvec{a}} \\
\VField{J}(\pvec{x}+\pvec{a}_i, \omega) & = &  \VField{J}(\pvec{x}, \omega)e^{i\pvec{k_{\mathrm{B}}} \pvec{a}} \\
\SField{\rho}(\pvec{x}+\pvec{a}_i, \omega) & = &  \SField{\rho}(\pvec{x}, \omega)e^{i\pvec{k_{\mathrm{B}}} \pvec{a}} \\
\end{eqnarray*}

This condition is called Bloch periodicity. The fields are periodic up to a common phase factor e^{i\pvec{k_{\mathrm{B}}} \pvec{a}} with a common Bloch vector \pvec{k}_{\mathrm{B}}.

Lattice vectors \pvec{a}_i are accompanied by reciprocal lattice vectors \pvec{b}_i of the same multiplicity, which satisfy

\begin{eqnarray*}
\pvec{b}_i \cdot \pvec{a}_j & = & 2 \pi \delta_{ij} \,.
\end{eqnarray*}

For a one-fold periodicity, i.e., with one single lattice vector \pvec{a}_1, we demand that \pvec{b}_1 is a multiple of \pvec{a}_1. Precisely, we have \pvec{b}_1=2 \pi \pvec{a}_1/|\pvec{a}_1|^2. In this case the unit cell may extend infinitely in the directions perpendicular to \pvec{a}_1 or one may impose additional boundary conditions as discussed above.

For a two-fold periodicity, one demands that the lattice vectors \pvec{a}_1, \pvec{a}_2 and the reciprocal lattice vectors \pvec{b}_1, \pvec{b}_2 span the same plane with normal \pvec{n}. Again, the unit cell may be infinite in the directions -\pvec{n} or +\pvec{n}, or it may be bounded by an additional boundary condition.

Depending quantities

The table below lists field quantities, which are implicitly defined within the solution fieldbag. This way they can be used by subsequent post-processes.

Quantity Expression JCM Tag
electric flux density \VField{D}=
\varepsilon
\VField{E} +
\tilde{\varepsilon}
\VField{H} ElectricFluxDensity
magnetic flux density \VField{B}=
\mu
\VField{H} +
\tilde{\mu}
\VField{E} MagneticFluxDensity
electric field energy density \frac{1}{2}
\VField{D} \cdot
\VField{E} ElectricFieldEnergyDensity
magnetic field energy density \frac{1}{2}
\VField{B} \cdot
\VField{H} MagneticFieldEnergyDensity
electromagnetic field energy flux (Poynting vector) \VField{S}=
\VField{E} \times
\VField{H} ElectromagneticFieldEnergyFluxDensity
[1]We use our own notation \tilde{\TField{\varepsilon}} and \tilde{\TField{\mu}} for the bi-anistropic material tensors due to the lack of a standard notation in the literature and a comfortable readability (\tilde{\TField{\varepsilon}} == “magnetic permittivity”, \tilde{\TField{\mu}} == “electric permeability”).