A Hybrid Phase Flow Method for Solving the Liouville Equation in a Bounded Domain
Hao Wu,Xu Yang
DOI: https://doi.org/10.1137/090773726
IF: 3.039
2011-01-01
SIAM Journal on Numerical Analysis
Abstract:The phase flow method, originally introduced in [L. X. Ying and E. J. Candes, J. Comput. Phys., 220 (2006), pp. 184-215], can efficiently solve autonomous ordinary differential equations. In [S. Jin, H. Wu, and Z. Y. Huang, SIAM J. Sci. Comput., 31 (2008), pp. 13031321], the method was generalized to solve Hamiltonian system where the Hamiltonian function was discontinuous. However, both of these methods require a phase flow map constructed on an invariant manifold. This can increase computational cost when the invariant domain is big or unbounded. Following the idea of [S. Jin, H. Wu, and Z. Y. Huang, SIAM J. Sci. Comput., 31 (2008), pp. 1303-1321], we propose a hybrid phase flow method for solving the Liouville equation in a bounded domain, which is smaller than the invariant manifold of a phase flow map. By using some proper boundary conditions, this method can help solve the problem where the invariant manifold of a phase flow map determined by the Liouville equation is unbounded. We verify numerical accuracy and efficiency by several examples of the semiclassical limit of the Schrodinger equation. Analysis of numerical stability and convergence is given for the semiclassical limit equation with inflow boundary condition.