Numerical Solution of Steady-State Electromagnetic Scattering Problems Using the Time-Dependent Maxwell's Equations

Allen Taflove, M.E. Brodwin

IEEE Transactions on Microwave Theory and Techniques · 1975 · 1.3K citations · 6 references

Concepts

TL;DR

The paper proposes a numerical method to compute electromagnetic fields in dielectric scatterers roughly one wavelength across, with a view toward extending it to three‑dimensional geometries. The approach treats the scatterer’s illumination as an initial‑value problem, launching a plane‑wave at t = 0 and repeatedly solving a finite‑difference time‑domain discretization of Maxwell’s equations until sinusoidal steady‑state is reached, from which the standing‑wave envelope is taken as the scattered field. Applied to a dielectric cylinder, the method locates and evaluates standing‑wave peaks with less than 10 % error in one minute of computation.

Abstract

A numerical method is described for the solution of the electromagnetic fields within an arbitrary dielectric scatterer of the order of one wavelength in diameter. The method treats the irradiation of the scatterer as an initial value problem. At t = 0, a plane-wave source of frequency f is assumed to be turned on. The diffraction of waves from this source is modeled by repeatedly solving a finite-difference analog of the time-dependent Maxwell's equations. Time stepping is continued until sinusoidual steady-state field values are observed at all points within the scatterer. The envelope of the standing wave is taken as the steady-state scattered field. As an example of this method, the computed results for a dielectric cylinder scatterer are presented. An error of less than /spl plusmn/10 percent in locating and evaluating the standing-wave peaks within the cylinder is achieved for a program execution time of 1 min. The extension of this method to the solution of the fields within three-dimensional dielectric scatterers is outlined.

References

6