We study the quasi-random choice method (QRCM) for the Liouville equation of ge- ometrical optics with discontinuous locM wave speed. This equation arises in the phase space computation of high frequency waves throu...We study the quasi-random choice method (QRCM) for the Liouville equation of ge- ometrical optics with discontinuous locM wave speed. This equation arises in the phase space computation of high frequency waves through interfaces, where waves undergo partial transmissions and reflections. The numerical challenges include interface, contact discon- tinuities, and measure-valued solutions. The so-called QRCM is a random choice method based on quasi-random sampling (a deterministic alternative to random sampling). The method not only is viscosity-free but also provides faster convergence rate. Therefore, it is appealing for the prob!em under study which is indeed a Hamiltonian flow. Our analy- sis and computational results show that the QRCM 1) is almost first-order accurate even with the aforementioned discontinuities; 2) gives sharp resolutions for all discontinuities encountered in the problem; and 3) for measure-valued solutions, does not need the level set decomposition for finite difference/volume methods with numerical viscosities.展开更多
文摘We study the quasi-random choice method (QRCM) for the Liouville equation of ge- ometrical optics with discontinuous locM wave speed. This equation arises in the phase space computation of high frequency waves through interfaces, where waves undergo partial transmissions and reflections. The numerical challenges include interface, contact discon- tinuities, and measure-valued solutions. The so-called QRCM is a random choice method based on quasi-random sampling (a deterministic alternative to random sampling). The method not only is viscosity-free but also provides faster convergence rate. Therefore, it is appealing for the prob!em under study which is indeed a Hamiltonian flow. Our analy- sis and computational results show that the QRCM 1) is almost first-order accurate even with the aforementioned discontinuities; 2) gives sharp resolutions for all discontinuities encountered in the problem; and 3) for measure-valued solutions, does not need the level set decomposition for finite difference/volume methods with numerical viscosities.