Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
1901000 | Wave Motion | 2006 | 22 Pages |
Abstract
We develop a numerical scheme for the wave front computation of complete transmissions and reflections in geometrical optics. Such a problem can be formulated by a reduced Liouville equation with a discontinuous local wave speed or index of refraction, arising in the high frequency limit of linear waves through inhomogeneous media. The key idea is to incorporate Snell's Law of Refraction into the numerical flux for the reduced Liouville equation. This scheme allows a hyperbolic CFL condition, under which positivity, and stabilities in both lâ and l1 norms, are established. Numerical experiments are carried out to demonstrate the validity and accuracy of this new scheme.
Keywords
Related Topics
Physical Sciences and Engineering
Earth and Planetary Sciences
Geology
Authors
Shi Jin, Xin Wen,