Electromagnetic Theory and Simulation Computational Electromagnetics Informational

What is the PML absorbing boundary condition and how does it work in FDTD simulation?

The Perfectly Matched Layer (PML) is an artificial absorbing material surrounding the FDTD simulation domain that absorbs outgoing electromagnetic waves without reflection, simulating an infinite open space. PML works by introducing a frequency-dependent conductivity that increases smoothly from zero at the inner boundary to a maximum at the outer boundary, attenuating waves as they propagate through the PML region. The key property is impedance matching: the PML is designed so that its wave impedance exactly matches free space (or the adjacent medium) at the interface, producing zero reflection for all frequencies, polarizations, and angles of incidence in the continuous case. In practice, the discretized PML produces small reflections that depend on: (1) PML thickness: typically 8-16 cells, thicker provides lower reflection. (2) Grading profile: polynomial or geometric conductivity increase. A cubic Profile (sigma proportional to d^3) provides a good balance between absorption and reflection. (3) Distance from scatterer: PML should be at least lambda/4 to lambda/2 from the nearest object to avoid evanescent field coupling. (4) Incidence angle: PML performance degrades for waves arriving at grazing angles. Typical PML reflection: -40 to -80 dB for well-configured setups. The stretched-coordinate PML (SC-PML or CPML) is the current standard implementation, handling both propagating and evanescent waves and being stable for long simulations.

Category: Electromagnetic Theory and Simulation

Updated: April 2026

Product Tie-In: Simulation Software, PCB Materials

PML Theory and Configuration

Before PML's invention by Jean-Pierre Berenger in 1994, FDTD simulations of open-boundary problems suffered from significant artificial reflections from the domain boundary. PML revolutionized computational electromagnetics by enabling accurate simulation of antennas, scattering, and radiation problems in finite computational domains.

Technical Considerations

PML introduces an anisotropic absorbing medium where the conductivity is applied independently in each coordinate direction. A wave propagating in the x-direction encounters x-directed conductivity sigma_x, which attenuates it exponentially: E = E_0 × exp(-sigma_x × x / (2 × epsilon_0 × c)). The mathematical trick is that the wave impedance remains equal to the adjacent medium (eta = sqrt(mu/epsilon) is unchanged because both mu and epsilon are modified by the same complex stretch factor), so there is zero reflection at the interface in the continuous formulation. The complex coordinate stretching formulation: x → x_tilde = x × (1 + sigma_x / (j × omega × epsilon_0)), transforms the real coordinate into a complex coordinate that causes exponential decay. This formulation is more general and numerically stable than the original split-field approach.

Performance verification: confirm specifications against the application requirements before finalizing the design
Environmental factors: temperature range, humidity, and vibration affect long-term reliability and parameter drift
Cost vs. performance: evaluate whether the application demands premium components or standard commercial grades
Interface compatibility: verify impedance, connector type, and mechanical form factor match the system architecture
Margin allocation: include sufficient design margin to account for manufacturing tolerances and aging effects

Performance Analysis

For FDTD simulation: (1) PML thickness: 8 cells for quick simulations, 12-16 cells for high accuracy. Each additional cell reduces reflection by approximately 6-10 dB. (2) Grading order: cubic polynomial (m=3) for the conductivity profile: sigma(d) = sigma_max × (d/d_total)^3, where d is distance from the inner PML boundary. (3) Maximum conductivity: sigma_max = -(m+1) × ln(R_0) / (2 × eta × d_total), where R_0 is the desired reflection coefficient at normal incidence (typically 10^-6 to 10^-8) and eta is the wave impedance. (4) Distance to objects: maintain at least lambda/4 (and preferably lambda/2) between any radiating or scattering object and the PML inner boundary. Objects too close to PML couple evanescent near-fields into the absorber, causing numerical artifacts. (5) Cornersand edges: PML regions overlap at simulation domain corners and edges, where two or three PML directions are simultaneously active. Modern implementations handle this automatically, but older codes may require explicit corner treatment.

Common Questions

Frequently Asked Questions

How many PML cells do I need?

8 cells minimum for engineering accuracy (-40 to -50 dB reflection). 12 cells for high accuracy (-60 to -70 dB). 16 cells for precision (-70 to -80 dB). Each doubling of PML cells reduces reflection by approximately 20 dB but doubles the PML memory and increases total simulation domain size. For most antenna and scattering problems, 10-12 cells with cubic grading provides adequate accuracy. For near-field calculations very close to PML (e.g., evanescent wave analysis), use 16+ cells and CFS-PML formulation.

What is the difference between PML and ABC?

Absorbing Boundary Conditions (ABCs) like Mur or Liao are single-layer boundary conditions applied at the domain edge that approximate the radiation condition. They are simpler (no additional cells needed) but provide limited absorption (-20 to -30 dB) and degrade significantly at non-normal incidence angles. PML is a multi-cell absorbing region that provides much better absorption (-40 to -80 dB) at all angles. PML has completely replaced ABCs in modern FDTD codes. The only remaining use for ABCs is in educational codes or extremely memory-constrained situations where the 8-16 cell PML overhead is unacceptable.

Can PML absorb DC or static fields?

Standard PML does not absorb DC (zero frequency) fields because the conductivity-based absorption requires propagating wave behavior. The CFS-PML formulation includes a real-valued frequency-shift parameter alpha that provides absorption even at DC. This is important for simulations involving charged particles, static field problems, or very low-frequency transient responses. Without CFS-PML, DC or very-low-frequency components reflect from the PML boundary and can cause non-physical late-time fields that corrupt the simulation results.

Keep Reading