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 HIGH ACCURACY NUMERICAL ALGORITHM FOR WAVEFRONT TIME COMPUTATION IN 3 D COMPLEX MEDIA

Performance Comparison of Optical Path Difference Calculation Algorithms for Wavelength Scanning Interferometry

A New Scheme for Fast Calculation of Seismic Traveltimes and Ray Paths in Heterogeneous Media

SIMULATION OF ACOUSTIC WAVE PROPAGATION IN COMPLEX MEDIA USING MRFD METHOD

Fast Wavefront Propagation (FWP) for Computing Exact Geodesic Distances on Meshes

Wavelet-Galerkin Method for Computations of Electromagnetic Fields

Calculation of the coupling coefficients for spatial waveguides using the field matching method

Numerical simulation of elastic wave based on the staggered grid finite difference method

A time-space domain stereo finite difference method for 3D scalar wave propagation

Finite difference frequency domain method with QR-decomposition-based complex-valued adaptive coefficients for 3D diffusive viscous wave modelling

An Efficient and high accuracy finite-difference scheme for the acoustic wave equation in 3D heterogeneous media

3-D Irregular Grids Finite Difference Method for Elastic Wave Propagation in Anisotropic Media

Seismic Propagation Simulation in Complex Media with Non-Rectangular Irregular-Grid Finite-Difference

Variable-order optimal implicit finite-difference schemes for explicit time-marching solutions to wave equations

A Three-dimensional Immersed Boundary Method for Accurate Simulation of Acoustic Wavefields with Complex Surface Topography

A Central Difference Method with Low Numerical Dispersion for Solving the Scalar Wave Equation

An Efficient Method for Estimating the Electromagnetic Wave Propagation in Three Dimensional Optical Waveguide Structures

An novel finite difference dispersion error elimination mechanism in the Lax–Wendroff high‐order time discretization

A Parallel Weighted Runge-Kutta Discontinuous Galerkin Method For Solving Acousitc Wave Equations In 3d D'Alembert Media On Unstructured Meshes

An Optimized Time-Space-domain Finite Difference Method with Piecewise Constant Interpolation Coefficients for Scalar Wave Propagation

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