An efficient method is proposed for numerical solutions of nonlinear SchrÃ¶dinger equations on an unbounded domain. Through approximating the kinetic energy term by a one-way equation and uniting it with the potential energy equation, absorbing boundary conditions are designed to truncate the unbounded domain, which are in nonlinear form and can perfectly absorb waves outgoing from the boundaries of the truncated computational domain. The stability of the induced initial boundary value problem defined on the computational domain is examined by a normal mode analysis. Numerical examples are given to illustrate the stable and tractable advantages of the method.
|Journal||Physical Review E|
|Publication status||Published - 29 Aug 2008|
Scopus Subject Areas
- Statistical and Nonlinear Physics
- Statistics and Probability
- Condensed Matter Physics