Computational high frequency waves through curved interfaces via the Liouville equation and geometric theory of diffraction
- Department of Mathematical Sciences, Tsinghua University, Beijing 100084 (China)
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.
- OSTI ID:
- 21159391
- Journal Information:
- Journal of Computational Physics, Vol. 227, Issue 12; Other Information: DOI: 10.1016/j.jcp.2008.02.029; PII: S0021-9991(08)00124-1; Copyright (c) 2008 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Tracker: A three-dimensional raytracing program for ionospheric radio propagation
Strong Stability Preserving Integrating Factor Two-Step Runge-Kutta Methods