Adaptive absorbing boundary conditions for Schrodinger-type equations: application to nonlinear and multi-dimensional problems

Zhenli Xu,Houde Han,Xiaonan Wu
DOI: https://doi.org/10.1016/j.jcp.2007.02.004
2007-02-11
Abstract:We propose an adaptive approach in picking the wave-number parameter of absorbing boundary conditions for Schrödinger-type equations. Based on the Gabor transform which captures local frequency information in the vicinity of artificial boundaries, the parameter is determined by an energy-weighted method and yields a quasi-optimal absorbing boundary conditions. It is shown that this approach can minimize reflected waves even when the wave function is composed of waves with different group velocities. We also extend the split local absorbing boundary (SLAB) method [Z. Xu and H. Han, {\it Phys. Rev. E}, 74(2006), pp. 037704] to problems in multidimensional nonlinear cases by coupling the adaptive approach. Numerical examples of nonlinear Schrödinger equations in one- and two dimensions are presented to demonstrate the properties of the discussed absorbing boundary conditions.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: when solving Schrödinger - type equations (especially nonlinear and multi - dimensional problems), how to effectively handle artificial boundary conditions to reduce the interference of reflected waves. Specifically, the author proposes a method for adaptively selecting wavenumber parameters to construct absorbing boundary conditions, thereby minimizing the influence of reflected waves. ### Problem Background In scientific and engineering applications, many physical phenomena can be described by partial differential equations, especially on unbounded domains. For example, the nonlinear Schrödinger equation can describe gravity waves in fluid mechanics, pulse propagation in optical fibers, and Bose - Einstein condensation at low temperatures. However, due to the unboundedness of the physical domain, one of the main difficulties in numerically solving these equations is how to handle artificial boundary conditions to ensure the accuracy of the calculation results. ### Research Objectives 1. **Reduce Reflected Waves**: For problems involving wave propagation, it is hoped that the artificial boundary conditions can eliminate all incident waves as much as possible, so that almost no or only very small reflected waves enter the internal region. 2. **Applicable to Nonlinear and Multi - Dimensional Problems**: Traditional absorbing boundary conditions have made a lot of progress in dealing with linear problems, but for nonlinear and multi - dimensional problems, effective solutions are still lacking. ### Solutions The author proposes an adaptive method based on the Gabor transform to select wavenumber parameters in the absorbing boundary conditions. The specific steps are as follows: 1. **Capture Local Frequency Information**: Use the Gabor transform to capture local frequency information near the artificial boundary. 2. **Energy - Weighted Method**: Determine the wavenumber parameters through the energy - weighted method, thereby constructing quasi - optimal absorbing boundary conditions. 3. **Extension to Multi - Dimensional and Nonlinear Cases**: Extend this method to multi - dimensional nonlinear Schrödinger equations and combine it with the time - splitting method (SLAB method) for numerical solution. ### Mathematical Formula Representation - Gabor transform: \[ \hat{\psi}(k, t)=\int_{x_{r}-b}^{x_{r}} \psi(x, t) e^{-ikx} \, dx \] where the window function is: \[ W(x)=\begin{cases} 1 & \text{if } x\in[x_{r}-b, x_{r}] \\ 0 & \text{otherwise} \end{cases} \] - Wavenumber Parameter Selection: \[ k_{0}=\frac{\int_{0}^{\infty}|\hat{\psi}(k, t)|^{p} k \, dk}{\int_{0}^{\infty}|\hat{\psi}(k, t)|^{p} \, dk} \] where \(p\) is a positive real number. ### Conclusion Through the above method, the author has successfully reduced the influence of reflected waves, especially when dealing with nonlinear and multi - dimensional Schrödinger equations. Numerical experiments have verified the effectiveness of this method and demonstrated its superior performance in different application scenarios.