A conservative spectral collocation method for the nonlinear Schrödinger equation in two dimensions
Introduction
The nonlinear Schrödinger (NLS) equation plays an important role in many fields of physics, such as plasma physics, quantum physics and nonlinear optics [4], [7]. In this paper, we consider the following cubic NLS equation on two-dimensional domain with periodic boundary conditions and initial condition Here is the complex unit, is a complex-valued function, Δ is the Laplace operator, and β is a given real constant. It is easy to show that the periodic-initial value problem (1) and (2) admits mass conservation law and the energy conservation law for t > 0.
Extensive numerical methods have been studied for the NLS equation in the literature. These methods include the spectral (pseudospectral) method [3], [23], finite difference method (FDM) [6], [8], [9], [25], [26], [28], finite element method (FEM) [1], [10], [12], discontinuous Galerkin (DG) method [15], [29], [30]. Recently, Li et al. [16] studied a sixth-order alternating direction implicit method for two-dimensional Schrödinger equation. Shi et al. [21] studied the superconvergence analysis of conforming finite element method for NLS equation. Kong et al. [14] proposed a compact and efficient conservative scheme for coupled NLS equations. Taleei and Dehghan [22] studied the pseudo-spectral domain decomposition method in space discretization and time-splitting method in time discretization for NLS equation. Katsaounis and Mitsotakisc [13] used an implicit-explicit type Crank–Nicolson finite element scheme to study the phenomenon of soliton reflection. For more numerical methods the reader may consult the Ref. [2].
It has been well known that spectral methods can provide a very useful tool for the solution NLS equation because of their spectral accuracy, when the geometry of the problem is smooth and regular [5]. The Fourier pseudo-spectral method is one of the spectral methods which is applied for various nonlinear Schrödinger equations by many authors such as Muruganandam and Adhikari [18], and Bao et al. [3]. Recently, Fourier [17] and Chebyshev [11], [22] spectral collocation method were proposed for solving NLS equation. These collocation methods have spectral accuracy and can be proved to keep conservation. However, these methods are limited in one-dimensional case. In this paper, we apply Fourier spectral collocation method for solving two-dimensional NLS equation. After the space discretization by the spectral differentiation matrices, the two-dimensional NLS equation is discretized into a system of nonlinear ordinary differential equations (ODEs) in matrix formulation. We will prove that the proposed method preserve the discrete mass and energy conservation laws.
In the time integration, the alternating direction implicit (ADI) schemes are usually applied in FDM discretization which consists of a number of tridiagonal matrix equations [9], [28]. However, the spectral collocation method is not appropriate for the ADI method because the spectral differentiation matrices are full. Here we consider the compact implicit integration factor (cIIF) method for time discretization [19], [20], [31]. In the space discretization, the diffusion terms are approximated by the spectral differentiation matrix. Then the cIIF method applies matrix exponential operations sequentially in x- and y-direction. As the results, we can calculate and store the exponential in small sizes. Another novel property of the methods is that the exact evaluation of the diffusion terms is decoupled from the implicit treatment of the nonlinear terms. We only solve a local nonlinear system at each spatial grid point.
The rest of the paper is organized as follows: in Section 2 we present the Fourier spectral collocation method combined with compact integration factor method to solve NLS equation in 2D. We also prove the mass and energy conservation laws in semi-discrete formulation. Numerical experiments are reported in Section 3. Finally, we summarize our conclusion in Section 4.
Section snippets
Fourier spectral collocation method
In the process of spectral collocation method, the essential part is the generation of the spectral differentiation matrix. We first give the Fourier spectral differentiation matrix on the interval [0, 2π]. Other intervals can be easily handled by a scale factor.
Let N be positive even integer. The spacing of the grid is and the collocation points are . We assume that f(x) is a function on [0, 2π]. Then f(x) can be interpolated by a sum of periodic sinc functions [24]
Numerical tests
In this section, we present numerical experiments on some examples to support conservation analysis and demonstrate convergence of our numerical schemes. The Fourier differentiation matrices can be generated using the algorithms discussed in [24]. Our cIIF scheme (27) relies on matrix exponentials of the spectral differentiation matrices which are small in every direction. In the calculation, the matrix exponentials are computed using a scaling and squaring algorithm with a Pade approximation
Conclusion
In this paper we have combined the Fourier spectral collocation method with the compact implicit integration factor method for computing the two-dimensional NLS equation. We have proved that the method can keep the conservation law in mass and energy. Moreover this method has spectral accuracy in space. The numerical examples confirmed the excellent qualities of our method. We introduce the cIIF method in time discretization. This method can reduce the computational cost significantly in both
Acknowledgments
R. Zhang's work was also supported by Brazilian Young Talent Attraction Program via National Council for Scientific and Technological Development (CNPq), the Foundation of Liaoning Educational Committee (Grant No. L201604). J. Zhu and A. Loula's work was partially supported by CNPq. X. Yu's work was supported by National Natural Science Foundation of China (Grant No. 11571002), the Science Foundation of China Academy of Engineering Physics (Grant No. 2013A0202011, 2015B0101021) and the Defense
References (31)
- et al.
Computational methods for the dynamics of the nonlinear Schrödinger/Gross–Pitaevskii equations
Comput. Phys. Commun.
(2013) - et al.
Numerical solution of the Gross–Pitaevskii equation for Bose–Einstein condensation
J. Comput. Phys.
(2003) - et al.
Difference schemes for solving the generalized nonlinear Schrödinger equation
J. Comput. Phys.
(1999) - et al.
A compact split-step finite difference method for solving the nonlinear Schrödinger equations with constant and variable coefficients
Comput. Phys. Commun.
(2010) - et al.
Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations
Appl. Numer. Math.
(2011) - et al.
B-spline finite element studies of the non-linear Schrödinger equation
Comput. Methods Appl. Mech. Eng.
(1993) - et al.
Numerical studies on nonlinear Schrödinger equations by spectral collocation method with preconditioning
J. Math. Anal. Appl.
(2007) - et al.
A combined discontinuous Galerkin method for the dipolar Bose–Einstein condensation
J. Comput. Phys.
(2014) - et al.
A spatial sixth-order alternating direction implicit method for two-dimensional cubic nonlinear Schrödinger equations
Comput. Phys. Commun.
(2015) - et al.
Efficient semi-implicit schemes for stiff systems
J. Comput. Phy.
(2006)
Compact integration factor methods in high spatial dimensions
J. Comput. Phys.
Superconvergence analysis of conforming finite element method for nonlinear Schrödinger equation
Appl. Math. Comput.
Time-splitting pseudo-spectral domain decomposition method for the soliton solutions of the one- and multi-dimensional nonlinear Schrödinger equations
Comput. Phys. Commun.
High-order compact ADI (HOC-ADI) method for solving unsteady 2D Schrödinger equation
Comput. Phys. Commun.
Compact finite difference schemes with high accuracy for one-dimensional nonlinear Schrödnger equation
Comput. Methods Appl. Mech. Eng.
Cited by (10)
Numerical analysis of fractional viscoelastic fluid problem solved by finite difference scheme
2022, Computers and Mathematics with ApplicationsCitation Excerpt :They can provide relatively accurate simulation results for engineering practice and save a lot of calculating time. So far, finite difference algorithm [31–33], finite element algorithm [34–36], finite volume algorithm [37–40], spectral method [41,42] and other numerical techniques are adopted to a great degree. Among them, finite difference method is widely used in the calculation of FPDEs because of its simple operation, high efficiency and accuracy.
Smooth Quintic spline approximation for nonlinear Schrödinger equations with variable coefficients in one and two dimensions
2018, Chaos, Solitons and FractalsCitation Excerpt :Two-dimensional Schrödinger equations are solved via differential quadrature method by Golbabai and Nikpour [49]. Rongpei Zhang et al. [51] presented a conservative Fourier spectral collocation (FSC) method to solve the two-dimensional nonlinear Schrödinger (NLS) equation. An effective differential quadrature method (DQM) which is based on modified cubic B- spline (MCB) has been implemented to obtain the numerical solutions for the nonlinear Schrödinger (NLS) equation by Bashan et al. [53].
An Efficient Parallel Laplace Transform Spectral Element Method for the Two-Dimensional Schrödinger Equation
2022, Mathematical Problems in EngineeringNumerical simulation of a generalized nonlinear derivative Schrödinger equation
2022, Electronic Research ArchiveA new high-order compact finite difference scheme based on precise integration method for the numerical simulation of parabolic equations
2020, Advances in Difference EquationsA conservative spectral method for the coupled Schrödinger-KdV equations
2020, International Journal of Modern Physics C