EM: Time-Harmonic Electromagnetic Fields

Understanding how electromagnetic waves behave is essential for designing everything from radio antennas to fiber-optic cables. However, directly solving Maxwell’s time-dependent equations for sinusoidal sources can be algebraically messy and time-consuming. Time-harmonic analysis provides a powerful shortcut by using complex number representations, or phasors, which transform differential calculus into simpler algebra. This method is the cornerstone of analyzing steady-state wave propagation, power transfer, and the resonant behavior of structures like waveguides and cavities.

The Phasor Transformation: From Time Derivatives to Algebra

The entire framework rests on a key assumption: all fields and sources are varying sinusoidally at a single angular frequency, denoted by $\omega$. This is the time-harmonic condition. In the real-world time domain, an electric field component might look like $E(x,y,z,t)=E_0(x,y,z)\cos (\omega t+\varphi )$. Working with such cosine functions and their derivatives is cumbersome.

The phasor technique simplifies this by using Euler’s formula, $e^{j\theta }=\cos \theta +j\sin \theta$. We represent the sinusoidal field by a complex phasor, $\tilde{\mathbf{E}}(x,y,z)$, which contains the amplitude and phase information but not the explicit time dependence. The rule for conversion is:

$$E(x,y,z,t)=\text{Re}\left\{\tilde{\mathbf{E}}(x,y,z)e^{j\omega t}\right\}$$

Here, $\tilde{\mathbf{E}}$ is a complex vector. The most powerful consequence of this representation is its effect on time derivatives. Taking the partial time derivative $\frac{\partial }{\partial t}$ of the phasor representation introduces a factor of $j\omega$. Therefore, in the phasor domain, the operation $\frac{\partial }{\partial t}$ is replaced by multiplication by $j\omega$.

Maxwell’s Equations in Algebraic Phasor Form

This transformation is applied directly to Maxwell’s equations. Consider Ampère’s Law in differential form for a region with conductivity $\sigma$, permittivity $ϵ$, and permeability $\mu$:

$$\nabla \times \mathbf{H}=\mathbf{J}+\frac{\partial \mathbf{D}}{\partial t}$$

The current density $\mathbf{J}$ has two parts: a source current ${\mathbf{J}}_s$ and a conduction current $\sigma \mathbf{E}$. Using the phasor forms $\tilde{\mathbf{H}}$, $\tilde{\mathbf{E}}$, and ${\tilde{\mathbf{J}}}_s$, and replacing $\frac{\partial }{\partial t}$ with $j\omega$, the equation becomes:

$$\nabla \times \tilde{\mathbf{H}}={\tilde{\mathbf{J}}}_s+\sigma \tilde{\mathbf{E}}+j\omega ϵ\tilde{\mathbf{E}}={\tilde{\mathbf{J}}}_s+j\omega {ϵ}_c\tilde{\mathbf{E}}$$

Here, we have combined terms into a complex permittivity, ${ϵ}_c=ϵ-j\frac{\sigma }{\omega }$, which neatly accounts for both dielectric polarization and conductive loss. The complete set of time-harmonic (phasor domain) Maxwell’s equations is:

$$\begin{array}{rl}\nabla \times \tilde{\mathbf{E}}& =-j\omega \mu \tilde{\mathbf{H}}\\ \nabla \times \tilde{\mathbf{H}}& ={\tilde{\mathbf{J}}}_s+j\omega {ϵ}_c\tilde{\mathbf{E}}\\ \nabla \cdot \tilde{\mathbf{D}}& ={\tilde{\rho }}_v\\ \nabla \cdot \tilde{\mathbf{B}}& =0\end{array}$$

Notice that the time derivatives are gone, replaced by multiplicative factors of $j\omega$. This turns the problem of solving coupled partial differential equations into a more manageable problem of solving coupled algebraic-differential equations in space.

The Helmholtz Wave Equation and Propagation

To find how waves propagate in a source-free region (${\tilde{\mathbf{J}}}_s=0,{\tilde{\rho }}_v=0$), we combine the phasor form of Maxwell’s curl equations. Taking the curl of $\nabla \times \tilde{\mathbf{E}}$ and substituting, we eliminate $\tilde{\mathbf{H}}$ and obtain the Helmholtz equation for the electric field:

$${\nabla }^2\tilde{\mathbf{E}}+k^2\tilde{\mathbf{E}}=0$$

The constant $k$ is the complex wavenumber, defined by $k^2={\omega }^2\mu {ϵ}_c$. In a lossless medium ($\sigma =0$), ${ϵ}_c=ϵ$ and $k=\omega \sqrt{\mu ϵ}=\frac{\omega }{v_p}$ is a real number, where $v_p$ is the phase velocity. The Helmholtz equation’s solutions describe propagating plane waves, such as $\tilde{\mathbf{E}}(z)={\tilde{\mathbf{E}}}_0{e}^{-jkz}$, and form the basis for analyzing waves in more complex structures. Solving this equation with appropriate boundary conditions is how you determine the modal field patterns inside waveguides and cavity resonators.

Power Flow: The Complex Poynting Vector

In circuit analysis, you use $V{I}^{\ast }$ to find complex power. The electromagnetic analogue is the complex Poynting vector, $\tilde{\mathbf{S}}$, defined as:

$$\tilde{\mathbf{S}}=\frac{1}{2}\tilde{\mathbf{E}}\times {\tilde{\mathbf{H}}}^{\ast }$$

The factor of $\frac{1}{2}$ appears because $\tilde{\mathbf{E}}$ and $\tilde{\mathbf{H}}$ are peak-valued phasors, and we want to average over time. The real part of $\tilde{\mathbf{S}}$ gives the time-averaged power density (in watts per square meter) flowing through a surface:

$${\mathbf{S}}_{avg}=\text{Re}\{\tilde{\mathbf{S}}\}=\frac{1}{2}\text{Re}\{\tilde{\mathbf{E}}\times {\tilde{\mathbf{H}}}^{\ast }\}$$

The imaginary part of $\tilde{\mathbf{S}}$ represents the time-averaged reactive power density, which oscillates between the field and the source without net energy transfer. Calculating the total time-averaged power through a surface, like the cross-section of a waveguide, involves integrating ${\mathbf{S}}_{avg}$ over that area. This is critical for determining the power-handling capability and efficiency of any radiating or guiding structure.

Key Applications: Antennas, Waveguides, and Cavities

The phasor methodology unlocks analysis of fundamental electromagnetic components. For antenna analysis, you solve for the radiated fields $\tilde{\mathbf{E}}$ and $\tilde{\mathbf{H}}$ from a known current distribution ${\tilde{\mathbf{J}}}_s$ on the antenna. You then use the complex Poynting vector to find the total radiated power and the radiation pattern.

In waveguide analysis, you solve the source-free Helmholtz equation subject to the boundary condition that the tangential electric field must be zero at the perfectly conducting walls. This leads to discrete, quantized solutions called modes (e.g., TE${}_{10}$, TM${}_{11}$), each with a unique cutoff frequency and field pattern. The propagation constant for each mode is derived from $k$.

A cavity resonator is essentially a waveguide closed at both ends. The boundaries force standing waves, and resonance occurs at specific frequencies where the cavity length is an integer multiple of half-wavelengths. Phasor analysis allows you to find these resonant frequencies and the corresponding three-dimensional field modes inside the cavity, which are vital for applications like filters and particle accelerators.

Common Pitfalls

1. Dropping the Complex Conjugate in the Poynting Vector: A frequent error is defining $\tilde{\mathbf{S}}$ as $\tilde{\mathbf{E}}\times \tilde{\mathbf{H}}$ without the complex conjugate. This is incorrect. The conjugate is essential to correctly compute time-averaged power. The correct formula is $\tilde{\mathbf{S}}=\frac{1}{2}\tilde{\mathbf{E}}\times {\tilde{\mathbf{H}}}^{\ast }$.
2. Misinterpreting Phasor Magnitudes: Remember that the magnitude of a field phasor, $|\tilde{\mathbf{E}}|$, represents the peak amplitude of the sinusoidal oscillation, not the RMS value. For power calculations, this is why the factor of $\frac{1}{2}$ appears, analogous to the difference between $\frac{1}{2}{V}_m{I}_m$ and $V_{rms}{I}_{rms}$ in AC circuit theory.
3. Applying Time-Harmonic Analysis to Transients: The phasor method only describes the steady-state response to a sinusoidal source. It cannot capture the initial transient behavior when a source is first switched on. Attempting to use it for transient analysis will yield incomplete or incorrect results.
4. Ignoring the Complex Nature of $k$ in Lossy Media: In a lossy medium, the wavenumber $k$ becomes complex ($k=\beta -j\alpha$). The real part $\beta$ is the phase constant, while the imaginary part $\alpha$ is the attenuation constant. Forgetting this and treating $k$ as purely real will lead to incorrect predictions of wave propagation and attenuation.

Summary

• Time-harmonic analysis uses phasors to represent sinusoidal fields, converting time derivatives in Maxwell's equations into multiplicative factors of $j\omega$, thus simplifying them to an algebraic form.
• In source-free regions, the phasor fields obey the Helmholtz wave equation (${\nabla }^2\tilde{\mathbf{E}}+k^2\tilde{\mathbf{E}}=0$), where the complex wavenumber $k$ governs wave propagation and attenuation.
• Time-averaged power flow is calculated using the complex Poynting vector, $\tilde{\mathbf{S}}=\frac{1}{2}\tilde{\mathbf{E}}\times {\tilde{\mathbf{H}}}^{\ast }$, where its real part gives the average power density.
• This formalism is the essential tool for solving the modal fields in waveguides, the radiation patterns of antennas, and the resonant frequencies of cavity resonators.
• Always remember the method's limitations: it is strictly a steady-state technique for single-frequency sinusoidal excitation.