Preprint 2005-024

Hamiltonian-Preserving Schemes for the Liouville Equation of Geometrical Optics with Transmissions and Reflections

Shi Jin and Xin Wen

Abstract: We construct a class of Hamiltonian-preserving numerical schemes for a Liouville equation of geometrical optics, with transmissions and reflections. This equation arises in the high frequency limit of the linear wave equation, with discontinuous local wave speeds. In our previous work we introduced the Hamiltonian-preserving schemes for the same equation when only complete transmissions or reflections occur at the interfaces. These schemes are extended in this paper to the general case of coexistence of both transmission and reflection satisfying the Snell Law of Refraction, with the correct transmission and reflection coefficients. This scheme allows a hyperbolic stability condition, under which positivity, and stabilities in both $l^1$ and $l^\infty$ norms, are established. Numerical experiments are carried out to study the numerical accuracy.

Available as PDF (1080 Kbytes)
Shi Jin, <>
Xin Wen, <>
Publishing information:
Submitted by:
<> May 13 2005.

[ 1996 | 1997 | 1998 | 1999 | 2000 | 2001 | 2002 | 2003 | 2004 | 2005 | All Preprints | Preprint Server Homepage ]
© The copyright for the following documents lies with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use.

Conservation Laws Preprint Server <>
Last modified: Sat May 14 11:13:17 MEST 2005