Split-step quintic B-spline collocation (SS5BC) methods are constructed for nonlinear Schrödinger equations in one, two and three dimensions in this paper. For high dimensions, new notations are introduced, which make the schemes more concise and achievable. The solvability, conservation and linear stability are discussed for the proposed methods. Numerical tests are carried out, and the present schemes are numerically verified to be convergent with second-order in time and fourth-order in space. The conserved quantity is also computed which agrees with the exact one. And solitary waves in one, two and three dimensions are simulated numerically which coincide with the exact ones. The SS5BC scheme is compared with the split-step cubic B-spline collocation (SS3BC) method in the numerical tests, and the former scheme is more efficient than the later one. Finally, the SS5BC scheme is also applied to compute Bose-Einstein condensates.
Citation: Shanshan Wang. Split-step quintic B-spline collocation methods for nonlinear Schrödinger equations[J]. AIMS Mathematics, 2023, 8(8): 19794-19815. doi: 10.3934/math.20231009
Split-step quintic B-spline collocation (SS5BC) methods are constructed for nonlinear Schrödinger equations in one, two and three dimensions in this paper. For high dimensions, new notations are introduced, which make the schemes more concise and achievable. The solvability, conservation and linear stability are discussed for the proposed methods. Numerical tests are carried out, and the present schemes are numerically verified to be convergent with second-order in time and fourth-order in space. The conserved quantity is also computed which agrees with the exact one. And solitary waves in one, two and three dimensions are simulated numerically which coincide with the exact ones. The SS5BC scheme is compared with the split-step cubic B-spline collocation (SS3BC) method in the numerical tests, and the former scheme is more efficient than the later one. Finally, the SS5BC scheme is also applied to compute Bose-Einstein condensates.
[1] | H. Wang, Numerical studies on the split-step finite difference method for nonlinear Schrödinger equations, Appl. Math. Comput., 170 (2005), 17–35. http://doi.org/10.1016/j.amc.2004.10.066 doi: 10.1016/j.amc.2004.10.066 |
[2] | H. Wang, X. Ma, J. Lu, W. Gao, An efficient time-splitting compact finite difference method for Gross-Pitaevskii equation, Appl. Math. Comput., 297 (2017), 131–144. http://dx.doi.org//10.1016/j.amc.2016.10.037 doi: 10.1016/j.amc.2016.10.037 |
[3] | L. Li, A split-step finite-element method for incompressible Navier-Stokes equations with high-order accuracy up-to the boundary, J. Comput. Phys., 408 (2020), 109274. http://doi.org/10.1016/j.jcp.2020.109274 doi: 10.1016/j.jcp.2020.109274 |
[4] | Y. Gao, L. Mei, Time-splitting Galerkin method for spin-orbit-coupled Bose-Einstein condensates, Comput. Math. Appl., 87 (2021), 77–90. http://doi.org/10.1016/j.camwa.2021.02.009 doi: 10.1016/j.camwa.2021.02.009 |
[5] | W. Bao, H. Li, J. Shen, A generalized-Laguerre-Fourier-Hermite pseudospectral method for computing the dynamics of rotating Bose-Einstein condensates, SIAM J. Sci. Comput., 31 (2009), 3685–3711. http://doi.org/10.1137/080739811 doi: 10.1137/080739811 |
[6] | W. Bao, Y. Cai, Y. Feng, Improved uniform error bounds on time-splitting methods for long-time dynamics of the nonlinear Klein-Gordon equation with weak nonlinearity, SIAM J. Numer. Anal., 60 (2022), 1962–1984. http://doi.org/10.1137/21M1449774 doi: 10.1137/21M1449774 |
[7] | S. Wang, L. Zhang, Split-step cubic B-spline collocation methods for nonlinear Schrödinger equations in one, two, and three dimensions with Neumann boundary conditions, Numer. Algorithms, 81 (2019), 1531–1546. http://doi.org/10.1007/s11075-019-00762-2 doi: 10.1007/s11075-019-00762-2 |
[8] | D. Irk, $\dot{I}$. Da$\check{g}$, Quintic B-spline collocation method for the generalized nonlinear Schrödinger equation, J. Frankl. Inst., 348 (2011), 378–392. http://doi.org/10.1016/j.jfranklin.2010.12.004 doi: 10.1016/j.jfranklin.2010.12.004 |
[9] | G. Strang, On the construction and comparison of difference schemes, SIAM J. Numer. Anal., 5 (1968), 506–517. https://doi.org/10.1137/0705041 doi: 10.1137/0705041 |
[10] | S. Yu, S. Zhao, G. W. Wei, Local spectral time splitting method for first- and second-order partial differential equations, J. Comput. Phys., 206 (2005), 727–780. http://doi.org/10.1016/j.jcp.2005.01.010 doi: 10.1016/j.jcp.2005.01.010 |
[11] | T. Wang, B. Guo, Q. Xu, Fourth-order compact and energy conservative difference schemes for the nonlinear Schrödinger equation in two dimensions, J. Comput. Phys., 243 (2013), 382–399. http://dx.doi.org/10.1016/j.jcp.2013.03.007 doi: 10.1016/j.jcp.2013.03.007 |
[12] | C. A. Hall, On error bounds for spline interpolation, J. Approx. Theory, 1 (1968), 209–218. https://doi.org/10.1016/0021-9045(68)90025-7 doi: 10.1016/0021-9045(68)90025-7 |
[13] | T. R. Taha, M. I. Ablowitz, Analytical and numerical aspects of certain nonlinear evolution equations. II: Numerical, nonlinear Schrödinger equations, J. Comput. Phys., 55 (1984), 203–230. http://doi.org/10.1016/0021-9991(84)90003-2 doi: 10.1016/0021-9991(84)90003-2 |
[14] | W. Bao, D. Jaksch, P. A. Markowich, Numerical solution of the Gross-Pitaevskii equation for Bose-Einstein condensation, J. Comput. Phys., 187 (2003), 318–342. http://doi.org/10.1016/S0021-9991(03)00102-5 doi: 10.1016/S0021-9991(03)00102-5 |