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.

GEOMETRIC OPTICS FOR 3D-HARTREE-TYPE EQUATION WITH COULOMB POTENTIAL

Geometric Optics for One-Dimensional Schrödinger-Poisson System

HMO in geometrical optics

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

GLOBAL GEOMETRICAL OPTICS APPROXIMATION TO THE SCHR(O)DINGER EQUATION

Global Geometrical Optics Approximation to the High Frequency Helmholtz Equation with Discontinuous Media

Global Geometrical Optics Method

Global Geometrical Optics Method for Vector-Valued Schrödinger Problems

Weakly Nonlinear Geometric Optics for Hyperbolic Systems of Conservation Laws

Global Solutions for Nonlinear Klein–Gordon Equations in Infinite Homogeneous Waveguides

Semilinear Geometric Optics for Generalized Solutions

ACTA PHOTONICA SINICA EXACT BOUND STATE SOLUTIONS OF N-DIMENSIONAL KLEIN-GORDON EQUATION FOR COULOMB POTENTIAL

Decay estimates for defocusing energy-critical Hartree equation

A Hamiltonian-Preserving Scheme for the Liouville Equation of Geometrical Optics with Partial Transmissions and Reflections.

ON GLOBAL SMOOTH SMALL DATA SOLUTIONS OF 3-D QUASILINEAR KLEIN-GORDON EQUATIONS ON R2 x T

Planck’s Theory About"geometrical Optics of Non-Homogeneous Bodies Relationships to Quantum Mechanics"

The energy-critical inhomogeneous generalized Hartree equation in 3D

Recovering a Potential from Cauchy Data via Complex Geometrical Optics Solutions

Generalized Hunter–Saxton equation and the geometry of the group of circle diffeomorphisms

Babich-Like Ansatz for Three-Dimensional Point-Source Maxwell's Equations in an Inhomogeneous Medium at High Frequencies.

Numerical analysis of eikonal equation