Numerical solution of the high frequency asymptotic expansion for the scalar wave equation
Journal of Computational Physics
Sticky Particles and Scalar Conservation Laws
SIAM Journal on Numerical Analysis
High-frequency wave propagation by the segment projection method
Journal of Computational Physics
Geometric optics in a phase-space-based level set and Eulerian framework
Journal of Computational Physics
Discretization of Dirac delta functions in level set methods
Journal of Computational Physics
Journal of Computational Physics
The numerical approximation of a delta function with application to level set methods
Journal of Computational Physics
Journal of Computational Physics
Journal of Computational Physics
SIAM Journal on Numerical Analysis
A fast phase space method for computing creeping rays
Journal of Computational Physics
High order numerical methods to a type of delta function integrals
Journal of Computational Physics
Journal of Computational Physics
High Order Numerical Methods to Three Dimensional Delta Function Integrals in Level Set Methods
SIAM Journal on Scientific Computing
Journal of Computational Physics
Hi-index | 31.46 |
We construct a class of numerical schemes for the Liouville equation of geometric optics coupled with the Geometric Theory of Diffractions to simulate the high frequency linear waves with a discontinuous index of refraction. In this work [S. Jin, X. Wen, A Hamiltonian-preserving scheme for the Liouville equation of geometric optics with partial transmissions and reflections, SIAM J. Numer. Anal. 44 (2006) 1801-1828], a Hamiltonian-preserving scheme for the Liouville equation was constructed to capture partial transmissions and reflections at the interfaces. This scheme is extended by incorporating diffraction terms derived from Geometric Theory of Diffraction into the numerical flux in order to capture diffraction at the interface. We give such a scheme for curved interfaces. This scheme is proved to be positive under a suitable time step constraint. Numerical experiments show that it can capture diffraction phenomena without fully resolving the wave length of the original wave equation.