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.展开更多
The solution of the Riemann Problem (RP) for the one-dimensional (1D) non-linear Shallow Water Equations (SWEs) is known to produce four potential wave patterns for the scenario where the water depth is always positiv...The solution of the Riemann Problem (RP) for the one-dimensional (1D) non-linear Shallow Water Equations (SWEs) is known to produce four potential wave patterns for the scenario where the water depth is always positive. In this paper, we choose four test problems with exact solutions for the 1D SWEs. Each test problem is a RP with one of the four possible wave patterns as its solution. These problems are numerically solved using schemes from the family of Weighted Essentially Non-Oscillatory (WENO) methods. For comparison purposes, we also include results obtained from the Random Choice Method (RCM). This study has three main objectives. Firstly, we outline the procedures for the implementation of the methods employed in this paper. Secondly, we assess the performance of the schemes in conjunction with a second-order Total Variation Diminishing (TVD) flux on a variety of RPs for the 1D SWEs (for both short- and long-time simulations). Thirdly, we investigate if a single method yields optimal outcomes for all test problems. Optimal outcomes refer to numerical solutions devoid of spurious oscillations, exhibiting high resolution of discontinuities, and attaining high-order accuracy in the smooth parts of the solution.展开更多
文摘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.
文摘The solution of the Riemann Problem (RP) for the one-dimensional (1D) non-linear Shallow Water Equations (SWEs) is known to produce four potential wave patterns for the scenario where the water depth is always positive. In this paper, we choose four test problems with exact solutions for the 1D SWEs. Each test problem is a RP with one of the four possible wave patterns as its solution. These problems are numerically solved using schemes from the family of Weighted Essentially Non-Oscillatory (WENO) methods. For comparison purposes, we also include results obtained from the Random Choice Method (RCM). This study has three main objectives. Firstly, we outline the procedures for the implementation of the methods employed in this paper. Secondly, we assess the performance of the schemes in conjunction with a second-order Total Variation Diminishing (TVD) flux on a variety of RPs for the 1D SWEs (for both short- and long-time simulations). Thirdly, we investigate if a single method yields optimal outcomes for all test problems. Optimal outcomes refer to numerical solutions devoid of spurious oscillations, exhibiting high resolution of discontinuities, and attaining high-order accuracy in the smooth parts of the solution.
