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.

N A ] 5 F eb 2 01 8 A Novel Sixth Order Energy-Conserved Method for Three-Dimensional Time-Domain Maxwell ’ s Equations

A Novel High‐order Time‐domain Scheme for Three‐dimensional Maxwell's Equations

The Symplectiness Of Maxwell'S Equations

Decomposition Methods for Time‐domain Maxwell's Equations

Time exponential integrator Fourier pseudospectral methods with high accuracy and multiple conservation laws for three-dimensional Maxwell's equations

An exact in time Fourier pseudospectral method with multiple conservation laws for three-dimensional Maxwell’s equations

Energy Conserving Higher Order Mixed Finite Element Discretizations of Maxwell's Equations

Simulating Maxwell–Schrödinger Equations by High-Order Symplectic FDTD Algorithm

A novel class of explicit divergence-free time-domain methods for efficiently solving Maxwell's equations

The High-Order Symplectic Finite-Difference Time-Domain Scheme

A Novel Matrix-Free Finite Element Method for Time-Harmonic Maxwell’s Equations

High-oder Symplectic FDTD Scheme for Solving Time-Dependent Schrdinger Equation

Energy-Conserving Hermite Methods for Maxwell's Equations

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

Time-domain methods based on Exponential Evolution Operator approximation for the solution of Maxwell's equations

High-order accurate schemes for Maxwell's equations with nonlinear active media and material interfaces

Optimal error estimates and energy conservation identities of the ADI-FDTD scheme on staggered grids for 3D Maxwell's equations

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

A 3-D Unconditionally Stable Precise Integration Time Domain Method in Circular Cylindrical Coordinates

Optimal Error Estimate for Energy-Preserving Splitting Schemes for Maxwell's Equations.