Abstract:We develop a level set method for the computation of multi-valued physical observables (density, velocity, energy, etc.) for the high frequency limit of symmetric hyperbolic systems in any number of space dimensions. We take two approaches to derive the method.The first one starts with a weakly coupled system of an eikonal equation for phase S and a transport equation for density ρ:∂tS+H(x,∇S)=0,(t,x)∈R+×Rn,∂tρ+∇x·(ρ∇pH(x,∇xS))=0.The main idea is to evolve the density near the n-dimensional bi-characteristic manifold of the eikonal (Hamiltonian–Jacobi) equation, which is identified as the common zeros of n level set functions in phase space (x,k)∈R2n. These level set functions are generated from solving the Liouville equation with initial data chosen to embed the phase gradient. Simultaneously, we track a new quantity f=ρ(t,x,k)|det(∇kϕ)| by solving again the Liouville equation near the obtained zero level set ϕ=0 but with initial density as initial data. The multi-valued density and higher moments are thus resolved by integrating f along the bi-characteristic manifold in the phase directions.The second one uses the high frequency limit of symmetric hyperbolic systems derived by the Wigner transform. This gives rise to Liouville equations in the phase space with measure-valued solution in its initial data. Due to the linearity of the Liouville equation we can decompose the density distribution into products of function, each of which solves the Liouville equation with L∞ initial data on any bounded domain. It yields higher order moments such as energy and energy flux.The main advantages of these new approaches, in contrast to the standard kinetic equation approach using the Liouville equation with a Dirac measure initial data, include: (1) the Liouville equations are solved with L∞ initial data, and a singular integral involving the Dirac-δ function is evaluated only in the post-processing step, thus avoiding oscillations and excessive numerical smearing; (2) a local level set method can be utilized to significantly reduce the computation in the phase space. These methods can be used to compute all physical observables for multi-dimensional problems.Our method applies to the wave fields corresponding to simple eigenvalues of the dispersion matrix. One such example is the wave equation, which will be studied numerically in this paper.

Computing High Frequency Solutions of Symmetric Hyperbolic Systems with Polarized Waves

A Gaussian Beam Method for High Frequency Solution of Symmetric Hyperbolic Systems with Polarized Waves.

Computing Multi-Valued Physical Observables for the High Frequency Limit of Symmetric Hyperbolic Systems

A Hamiltonian-Preserving Scheme for High Frequency Elastic Waves in Heterogeneous Media

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

Polarized High-frequency Wave Propagation Beyond the Nonlinear Schrödinger Approximation

Numerical methods for hyperbolic systems with singular coe-cients: well-balanced scheme, Hamiltonian preservation, and beyond

Highly efficient Gauss's law-preserving spectral algorithms for Maxwell's double-curl source and eigenvalue problems based on eigen-decomposition

Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems

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

A high order discontinuous Galerkin method for the symmetric form of the anisotropic viscoelastic wave equation

Energy-Conserving Hermite Methods for Maxwell's Equations

An Accurate Method for Dispersion Characteristics of Surface Waves in Layered Anisotropic Semi-Infinite Spaces

Exact Solutions to the (1+1)-Dimensional Nonlinear Maxwell Equations in the Orthogonal Curvilinear Coordinates

Numerical Solution of the Helmholtz Equation with High Wavenumbers

Computational High Frequency Waves Through Curved Interfaces Via the Liouville Equation and Geometric Theory of Diffraction

Evanescent Plane Wave Approximation of Helmholtz Solutions in Spherical Domains

Exact Solutions to the Nonlinear Maxwell Equations in the Orthogonal Curvilinear Coordinates

Stable and Accurate Computation of Dispersion Relations for Layered Waveguides, Semi-Infinite Spaces and Infinite Spaces

Hamiltonian System and Algebro-Geometric Solution Associated with Dispersive Long Wave Equation

Hybrid Fourier-Continuation Method and Weighted Essentially Non-oscillatory Finite Difference Scheme for Hyperbolic Conservation Laws in a Single-Domain Framework