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.

Computational high-frequency wave propogation using the level-set method with applications to the semi-classical limit of the Schrödinger equations

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

Computation of the Semiclassical Limit of the Schrödinger Equation with Phase Shift by a Level Set Method

A level set method for the semiclassical limit of the Schrödinger equation with discontinuous potentials

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

Computing High Frequency Solutions of Symmetric Hyperbolic Systems with Polarized Waves

Mathematical and Computational Methods for Semiclassical Schrödinger Equations

Computing Multivalued Physical Observables for the Semiclassical Limit of the Schrödinger Equation

Hamiltonian-preserving Schemes for the Liouville Equation of Geometrical Optics with Discontinuous Local Wave Speeds

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

Numerical Simulation for the Wave of the Variable Coefficient Nonlinear Schrödinger Equation Based on the Lattice Boltzmann Method

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

A compact finite difference scheme for the nonlinear Schrödinger equation with wave operator.

Systematic exploration of solitary wave characteristics for the high-order dispersive extended nonlinear Schrödinger model

Numerical L Methods for Schrödinger Equations

Symplectic and Multi-Symplectic Wavelet Collocation Methods for Two-Dimensional Schrödinger Equations

A LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS

Gaussian Beam Methods for the Schrodinger Equation in the Semi-classical Regime

Numerical Study of Time-Splitting Spectral Discretizations of Nonlinear Schrödinger Equations in the Semiclassical Regimes

Multi-phase Computations of the Semiclassical Limit of the Schrödinger Equation and Related Problems: Whitham Vs Wigner

Stationary Schrödinger equation in the semi-classical limit: numerical coupling of oscillatory and evanescent regions