Abstract:In some applications, it is reasonable to assume that geodesics (rays) have a consistent orientationso that Maxwell's equations may be viewed as an evolution equation in one of the spatial directions.With such applications in mind, we propose a new Eulerian geometrical-optics method,dubbed the fast Huygens sweeping method, for computing Green's functions of Maxwell's equationsin inhomogeneous media in the high-frequency regime and in the presence of caustics.The first novelty of the fast Huygens sweeping method is that a new dyadic-tensor-type geometrical-optics ansatz isproposed for Green's functions which is able to utilize some unique features of Maxwell's equations.The second novelty is that the Huygens--Kirchhoff secondary source principle is usedto integrate many locally valid asymptotic solutions to yield a globally valid asymptotic solution sothat caustics associated with the usual geometrical-optics ansatz can be treated automatically.The third novelty is that a butterfly algorithm is adapted to carry out the matrix-vector productsinduced by the Huygens--Kirchhoff integration in $O(N\\log N)$ operations, where $N$ is the total numberof mesh points, and the proportionality constant depends on the desired accuracy and is independentof the frequency parameter. To reduce the storage of the resulting traveltime and amplitude tables, wecompress each table into a linear combination of tensor-product based multivariate Chebyshev polynomialsso that the information of each table is encoded into a small number of Chebyshev coefficients.The new method enjoys the following desired features:(1) it precomputes a set of local traveltime and amplitude tables;(2) it automatically takes care of caustics;(3) it constructs Green's functions of Maxwell's equations for arbitrary frequencies and for many point sources;(4) for a specified number of points per wavelength it constructs each Green's function in nearly optimalcomplexity $O(N\\log N)$ in terms of the total number of mesh points $N$, where the prefactor of the complexitydepends only on the specified accuracy and is independent of the frequency parameter.Three-dimensional numerical experiments are presented to demonstrate the performance and accuracy of the new method.

A Simple Green Function for Diffraction-Radiation of Time-Harmonic Waves with Forward Speed

Eigenfunction expansion of dyadic Green's function in biisotropic medium and its applications

A ‘‘compact Green function’’ approach to the time‐domain direct and inverse problems for a stratified dissipative slab

Spurious Mode Free Broadband Green’s Function Technique for Periodic Scatterers Using the Combined Field Integral Equation Formulation

Generalized Reflection and Transmission of Beams

Green's function of periodic scatterers applied to scattering by finite periodic arrays of scatterers

Wave analysis in the complex Fourier transform domain: A new method to obtain the Green's functions of dispersive linear partial differential equations

Scattering of Waves by a Half-Space of Periodic Scatterers Using Broadband Green's Function.

Propagation of waves from an arbitrary shaped surface -- a generalization of the Fresnel diffraction integral

Time Domain Green Functions Technique for a Point Source over a Dissipative Stratified Half-Space with a Phase Velocity Mismatch at the Surface

Exact boundary-value solution for an electromagnetic wave propagating in a linearly-varying index of refraction

Frequency and Time Domain Green-Function Technique for Nonuniform Lcrg Transmission-Lines with Frequency-Dependent Parameters

Computational High Frequency Wave Diffraction by a Corner Via the Liouville Equation and Geometric Theory of Diffraction

Green'S Functions, Including Scatterers, for Photonic Crystals and Metamaterials

Mathematical Analysis of Wave Radiation by a Step-Like Surface

A time‐harmonic Green’s function technique and wave propagation in a stratified nonreciprocal chiral slab with multiple discontinuities

Rapidly convergent quasi-periodic Green functions for scattering by arrays of cylinders---including Wood anomalies

Truncated Hadamard-Babich Ansatz and Fast Huygens Sweeping Methods for Time-Harmonic Elastic Wave Equations in Inhomogeneous Media in the Asymptotic Regime

Eulerian Geometrical Optics and Fast Huygens Sweeping Methods for Three-Dimensional Time-Harmonic High-Frequency Maxwell's Equations in Inhomogeneous Media.

Three-dimensional quasi-periodic shifted Green function throughout the spectrum--including Wood anomalies

Non-Reflecting Artificial Boundaries For Transient Scalar Wave Propagation In A Two-Dimensional Infinite Homogeneous Layer