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

We propose an adaptive approach in picking the wave-number parameter of absorbing boundary conditions for Schr\"{o}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\"{o}dinger equations in one- and two dimensions are presented to demonstrate the properties of the discussed absorbing boundary conditions.Comment: 18 pages; 12 figures. A short movie for the 2D NLS equation with absorbing boundary conditions can be downloaded at http://home.ustc.edu.cn/~xuzl/movie.avi. To appear in Journal of Computational Physic

An approximate solution of the MHD Falkner-Skan flow by Hermite functions pseudospectral method

Based on a new approximation method, namely pseudospectral method, a solution for the three order nonlinear ordinary differential laminar boundary layer Falkner-Skan equation has been obtained on the semi-infinite domain. The proposed approach is equipped by the orthogonal Hermite functions that have perfect properties to achieve this goal. This method solves the problem on the semi-infinite domain without truncating it to a finite domain and transforming domain of the problem to a finite domain. In addition, this method reduces solution of the problem to solution of a system of algebraic equations. We also present the comparison of this work with numerical results and show that the present method is applicable.Comment: 15 pages, 4 figures; Published online in the journal of "Communications in Nonlinear Science and Numerical Simulation

Structural stability of finite dispersion-relation preserving schemes

The goal of this work is to determine classes of travelling solitary wave solutions for a differential approximation of a finite difference scheme by means of a hyperbolic ansatz. It is shown that spurious solitary waves can occur in finite-difference solutions of nonlinear wave equation. The occurance of such a spurious solitary wave, which exhibits a very long life time, results in a non-vanishing numerical error for arbitrary time in unbounded numerical domain. Such a behavior is referred here to has a structural instability of the scheme, since the space of solutions spanned by the numerical scheme encompasses types of solutions (solitary waves in the present case) that are not solution of the original continuous equations. This paper extends our previous work about classical schemes to dispersion-relation preserving schemes

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters