Electromagnetic waves

This page describes the available boundary conditions, port types, and interactions for electromagnetic wave propagation in Allsolve. The weak formulation is stated up front; the full derivation is at the end of this page.

Weak formulation

The weak form solved by Allsolve is

\begin{align} -\int_{\Gamma}\ (\boldsymbol{n} \times \frac{\partial \boldsymbol{H}}{\partial t}) \cdot \boldsymbol{E}^{\prime}\ d\Gamma + \int_{\Omega}\ (\boldsymbol{\mu}^{-1} \nabla \times \boldsymbol{E}) \cdot \ (\nabla \times \boldsymbol{E}^{\prime})\ d\Omega \\[10pt] + \int_{\Omega}\ \boldsymbol{\sigma} \frac{\partial \boldsymbol{E}}{\partial t} \cdot \boldsymbol{E}^{\prime}\ d\Omega + \int_{\Omega}\ \boldsymbol{\epsilon} \frac{\partial^2 \boldsymbol{E}}{\partial t^2} \cdot \boldsymbol{E}^{\prime}\ d\Omega = 0. \end{align}

where $\boldsymbol{E}$ is the electric field, $\boldsymbol{H}$ the magnetic field intensity, $\mu$ the magnetic permeability, $\sigma$ the conductivity, and $\epsilon$ the permittivity.

See Formulation derivation for the full strong-to-weak derivation.

Interactions

Constraint

Fixes the $\boldsymbol{E}$ field value at the defined region, enforcing idealized behavior such as perfect reflection.


How to use	You can use either the matrix editor or the expression editor. With the matrix editor, provide a `3×1` matrix for a 3D problem (`2×1` for 2D, etc.) with electric field strength vector values for each dimension. In the expression editor, write the direct expression.
Example	`[1; 5; 0]` applies a 3D electric field strength constraint of $E_x = 1\:\frac{V}{m}$ , $E_y = 5\:\frac{V}{m}$ and $E_z = 0\:\frac{V}{m}$ to the specified region. Same principles apply for 2D.
Unit	Electric field strength in Volts per meter (V/m)

Perfectly matched layer

Implements an absorbing boundary region to simulate infinite space and prevent reflections from truncated boundaries. Placed at the outer boundary of the air domain, it absorbs outgoing waves without reflection.

Key properties:

Used to simulate open boundaries for radiating structures.

Two types available:

AML type: Suitable for smooth boundaries.
Box PML: Suitable for rectangular or cornered boundaries.
The number of PML layers can be tuned in the Shared PML settings.

Rectangular port

A special case of the Eigenmode port, used when the port geometry and material are well-defined and homogeneous. Provides analytical solutions for the field distribution under certain geometric assumptions.

Assumptions:

The port surface must be rectangular.
The dielectric across the port must be uniform.
All four port edges must be bounded by perfect conductors.

Key properties:

Perfectly absorbing (impedance matched, no reflection).
Provides electric field $\boldsymbol{E}$ and magnetic field $\boldsymbol{H}$ access, but not voltage $V$ or current $I$ .

Perfect conductor

Imposes a Dirichlet boundary condition on the electric field. Sets $\boldsymbol{E} = 0$ on the target boundary or region.

Key properties:

Models highly conductive materials as ideal conductors.
Reduces the number of unknowns and computational cost.

Eigenmode port

Defines a port that perfectly absorbs the excited mode with matched impedance (no reflection). This port type is typically used when the cut-off frequency or mode shapes are unknown. It supports both low and high-frequency signals.

Key properties:

No voltage $V$ or current $I$ output — only electric field $\boldsymbol{E}$ and magnetic field $\boldsymbol{H}$ are available.

Dielectric loss

As the EM wave propagates through a dielectric, some of the energy is absorbed by the dielectric and is dissipated as heat. This loss of EM energy is simulated using the Dielectric loss interaction.

The target region for this interaction is the dielectric region. In Allsolve, the dielectric loss is parameterized using the loss tangent ( $\tan \delta$ ):

$\tan \delta = \varepsilon^{\prime\prime} / \varepsilon^{\prime}$

$\varepsilon^{\prime}$ is the real part of the complex permittivity.
$\varepsilon^{\prime\prime}$ is the negative of the imaginary part of the complex permittivity.
The complex electric permittivity is given by $\varepsilon = \varepsilon^{\prime} - j\varepsilon^{\prime\prime}$

For script users:

The predefinedemwave function in Allsolve takes as arguments, among many, the real part and imaginary part of the complex permittivity.
The real part ( $\varepsilon^{\prime}$ ) is given by the permittivity defined in the material properties of the dielectric.
The imaginary part ( $-\varepsilon^{\prime\prime}$ ) is defined as $-\tan \delta * \varepsilon^{\prime}$

Lump V/I

Defines a lumped voltage or current source for feeding compact structures where the port dimensions are much smaller than the wavelength. This type of port assumes that the electric potential $v$ is physically meaningful and that the local electric field is derived from $\boldsymbol{E} = –\nabla v$ .

Key properties:

Provides direct access to voltage $V$ and current $I$ at the port.
Allows connecting external circuit elements (RLC components).
Suitable for compact feed structures such as microstrip, CPW, or GCPW.

Periodicity

Imposes periodic boundary conditions on the electromagnetic $\boldsymbol{E}$ and $\boldsymbol{H}$ fields between two boundaries. Reduces the computational domain size for geometrically symmetric or antisymmetric problems, avoiding the need to model the full geometry.


Example	Periodicity follows the same principles, regardless of which physics module it belongs to. See how periodicity is used in magnetism as a reference: Periodicity in electric motors

Absorbing boundary

Applies a first-order absorbing boundary condition to the outer boundary of the electromagnetic domain. Allows outgoing waves to leave the computational domain with minimal reflection. Lighter-weight than a Perfectly matched layer — suitable for cases where some residual reflection is acceptable in exchange for simpler setup and fewer mesh elements.


How to use	Select the outer boundary where waves should exit. No additional parameters are needed — the condition approximates a matched impedance termination.

Boundary admittance

Applies a surface admittance condition that characterizes how easily current passes through a conductive boundary. Admittance ( $Y$ ) is the reciprocal of impedance ( $Z = 1/Y$ ), expressing the same physics in terms of conductance rather than resistance.


How to use	Select the boundary surface of the conductor and set the boundary admittance.
Type options	Admittance — specify the real (Yr) and imaginary (Yi) parts of the boundary admittance directly. Good conductor — specify the conductivity and magnetic permeability (defaults to `mu0`) of the conductor material.
Unit	Boundary admittance Yr and Yi in Siemens (S)

Compatibilities with analysis types

Electromagnetic waves supports Harmonic and Eigenmode analysis. Static analysis is not applicable. Multiharmonic and Transient analysis are not supported.

See the full compatibility table in Simulations overview.

Use cases

Formulation derivation

Strong formulation

We start from Maxwell’s equations along with the material relations

\begin{align*} \boldsymbol{B} &= \boldsymbol{\mu} \boldsymbol{H} + \boldsymbol{B_r} \tag{1} \\[5pt] \boldsymbol{D} &= \boldsymbol{\varepsilon} \boldsymbol{E} + \boldsymbol{D_r} \tag{2} \\[5pt] \boldsymbol{J} &= \boldsymbol{\sigma} \boldsymbol{E} + \boldsymbol{J_r}\ . \tag{3} \end{align*}

Where $\boldsymbol{\varepsilon}$ is the electric permittivity tensor, $\boldsymbol{\mu}$ is the magnetic permeability tensor and $\boldsymbol{\sigma}$ is the electric conductivity tensor. Quantities $\boldsymbol{D_r}$ , $\boldsymbol{B_r}$ and $\boldsymbol{J_r}$ are typically associated with remanent effects. We assume that $\boldsymbol{\epsilon}$ , $\boldsymbol{\mu}$ , $\boldsymbol{\sigma}$ , $\boldsymbol{D_r}$ , $\boldsymbol{B_r}$ , and $\boldsymbol{J_r}$ change slowly compared to the electromagnetic wave frequency, so their time derivatives can be neglected.

Under these assumptions, Faraday’s law becomes

\begin{align*} \nabla \times \boldsymbol{E} &= -\boldsymbol{\mu} \frac{\partial \boldsymbol{H}}{\partial t} \tag{4} \\[10pt] \boldsymbol{\mu}^{-1} (\nabla \times \boldsymbol{E}) &= -\frac{\partial \boldsymbol{H}}{\partial t}. \tag{5} \end{align*}

We apply the curl operator to both sides and assume the curl and time derivative can be interchanged using the time derivative of a curl identity

\begin{align*} \nabla \times (\boldsymbol{\mu}^{-1} (\nabla \times \boldsymbol{E})) &= -\frac{\partial}{\partial t} (\nabla \times \boldsymbol{H}). \tag{6} \end{align*}

Substituting in the Ampère-Maxwell, material relations $\text{(2)}$ and $\text{(3)}$ and by rearranging we obtain

\begin{align*} \nabla \times (\boldsymbol{\mu}^{-1} (\nabla \times \boldsymbol{E})) + \boldsymbol{\sigma} \frac{\partial \boldsymbol{E}}{\partial t} + \boldsymbol{\varepsilon} \frac{\partial^2 \boldsymbol{E}}{\partial t^2} &= 0. \tag{7} \end{align*}

If we assume isotropic material, all tensor quantities become scalars, the equation results in

\begin{align*} \nabla \times (\frac{1}{\mu} (\nabla \times \boldsymbol{E})) + \sigma \frac{\partial \boldsymbol{E}}{\partial t} + \varepsilon \frac{\partial^2 \boldsymbol{E}}{\partial t^2} &= 0. \tag{8} \end{align*}

Deriving the weak form

To obtain the weak formulation of $\text{(7)}$ , multiply by a test function $\boldsymbol{E}'$ and integrate over the domain $\Omega$ :

\begin{align*} \int_\Omega\ \nabla \times (\boldsymbol{\mu}^{-1} \, \nabla \times \boldsymbol{E}) \cdot \boldsymbol{E}' \, d\Omega + \int_\Omega\ \boldsymbol{\sigma} \frac{\partial \, \boldsymbol{E}}{\partial t} \cdot \boldsymbol{E}' \, d\Omega + \int_\Omega\ \boldsymbol{\epsilon} \frac{\partial^2 \, \boldsymbol{E}}{\partial t^2} \cdot \boldsymbol{E}' \, d\Omega &= 0. \tag{9} \end{align*}

We can rewrite the first term using divergence of cross product

\begin{align*} \int_{\Omega}\ \nabla \cdot ((\boldsymbol{\mu}^{-1} \nabla \times \boldsymbol{E}) \times \boldsymbol{E}^{\prime})\ d\Omega + \int_{\Omega}\ (\boldsymbol{\mu}^{-1} \nabla \times \boldsymbol{E}) \cdot \ (\nabla \times \boldsymbol{E}^{\prime}), \tag{10} \end{align*}

Applying divergence theorem on the divergence term:

\begin{align*} \int_{\Gamma}\ ((\boldsymbol{\mu}^{-1} \nabla \times \boldsymbol{E}) \times \boldsymbol{E}^{\prime}) \cdot \boldsymbol{n}\ d\Gamma + \int_{\Omega}\ (\boldsymbol{\mu}^{-1} \nabla \times \boldsymbol{E}) \cdot \ (\nabla \times \boldsymbol{E}^{\prime})\ d\Omega. \tag{11} \end{align*}

Substituting in the Faraday’s law and using the material relation $\text{(1)}$ we obtain

\begin{align*} -\int_{\Gamma}\ (\frac{\partial \boldsymbol{H}}{\partial t} \times \boldsymbol{E}^{\prime}) \cdot \boldsymbol{n}\ d\Gamma + \int_{\Omega}\ (\boldsymbol{\mu}^{-1} \nabla \times \boldsymbol{E}) \cdot \ (\nabla \times \boldsymbol{E}^{\prime})\ d\Omega. \tag{12} \end{align*}

Finally, by using the scalar triple product identity we derive the final weak formulation stated at the top of this page.