Applications of Mathematics, Vol. 66, No. 2, pp. 213-232, 2021


A spatially sixth-order hybrid $L1$-CCD method for solving time fractional Schrödinger equations

Chun-Hua Zhang, Jun-Wei Jin, Hai-Wei Sun, Qin Sheng

Received December 4, 2019.   Published online December 16, 2020.

Abstract:  We consider highly accurate schemes for nonlinear time fractional Schrödinger equations (NTFSEs). While an $L1$ strategy is employed for approximating the Caputo fractional derivative in the temporal direction, compact CCD finite difference approaches are incorporated in the space. A highly effective hybrid $L1$-CCD method is implemented successfully. The accuracy of this linearized scheme is order six in space, and order $2-\gamma$ in time, where $0<\gamma<1$ is the order of the Caputo fractional derivative involved. It is proved rigorously that the hybrid numerical method accomplished is unconditionally stable in the Fourier sense. Numerical experiments are carried out with typical testing problems to validate the effectiveness of the new algorithms.
Keywords:  nonlinear time fractional Schrödinger equations; $L1$ formula; hybrid compact difference method; linearization; unconditional stability
Classification MSC:  65M06, 65M20, 65M60


References:
[1] D. Baleanu, K. Diethelm, E. Scalas, J. J. Trujillo: Fractional Calculus: Models and Numerical Methods. Series on Complexity, Nonlinearity and Chaos 5. World Scientific, Hackensack (2012). DOI 10.1142/8180 | MR 2894576 | Zbl 1248.26011
[2] D. A. Benson, R. Schumer, M. M. Meerschaert, S. W. Wheatcraft: Fractional dispersion, Lévy motion, and the MADE tracer tests. Transp. Porous Media 42 (2001), 211-240. DOI 10.1023/A:1006733002131 | MR 1948593
[3] A. H. Bhrawy, E. H. Doha, S. S. Ezz-Eldien, R. A. Van Gorder: A new Jacobi spectral collocation method for solving 1+1 fractional Schrödinger equations and fractional coupled Schrödinger systems. Eur. Phys. J. Plus 129 (2014), Article ID 260, 21 pages. DOI 10.1140/epjp/i2014-14260-6
[4] A. G. Bratsos: A linearized finite-difference scheme for the numerical solution of the nonlinear cubic Schrödinger equation. Korean J. Comput. Appl. Math. 8 (2001), 459-467. DOI 10.1007/BF02941979 | MR 1848910 | Zbl 1015.65042
[5] Q. Chang, E. Jia, W. Sun: Difference schemes for solving the generalized nonlinear Schrödinger equation. J. Comput. Phys. 148 (1999), 397-415. DOI 10.1006/jcph.1998.6120 | MR 1669707 | Zbl 0923.65059
[6] X. Chen, Y. Di, J. Duan, D. Li: Linearized compact ADI schemes for nonlinear time-fractional Schrödinger equations. Appl. Math. Lett. 84 (2018), 160-167. DOI 10.1016/j.aml.2018.05.007 | MR 3808513 | Zbl 06892656
[7] B. Chen, D. He, K. Pan: A linearized high-order combined compact difference scheme for multi-dimensional coupled Burgers' equations. Numer. Math., Theory Methods Appl. 11 (2018), 299-320. DOI 10.4208/nmtma.OA-2017-0090 | MR 3849456 | Zbl 1424.65124
[8] B. Chen, D. He, K. Pan: A CCD-ADI method for two-dimensional linear and nonlinear hyperbolic telegraph equations with variable coefficients. Int. J. Comput. Math. 96 (2019), 992-1004. DOI 10.1080/00207160.2018.1478415 | MR 3911287
[9] P. C. Chu, C. Fan: A three-point combined compact difference scheme. J. Comput. Phys. 140 (1998), 370-399. DOI 10.1006/jcph.1998.5899 | MR 1616146 | Zbl 0923.65071
[10] M. Cui: Compact finite difference method for the fractional diffusion equation. J. Comput. Phys. 228 (2009), 7792-7804. DOI 10.1016/j.jcp.2009.07.021 | MR 2561843 | Zbl 1179.65107
[11] M. Dehghan, A. Taleei: A compact split-step finite difference method for solving the nonlinear Schrödinger equations with constant and variable coefficients. Comput. Phys. Commun. 181 (2010), 43-51. DOI 10.1016/j.cpc.2009.08.015 | MR 2575184 | Zbl 1206.65207
[12] R. P. Feynman, A. R. Hibbs: Quantum Mechanics and Path Integrals. Dover Publications, New York (2010). MR 2797644 | Zbl 1220.81156
[13] G.-H. Gao, H.-W. Sun: Three-point combined compact alternating direction implicit difference schemes for two-dimensional time-fractional advection-diffusion equations. Commun. Comput. Phys. 17 (2015), 487-509. DOI 10.4208/cicp.180314.010914a | MR 3371511 | Zbl 1388.65052
[14] G.-H. Gao, H.-W. Sun: Three-point combined compact difference schemes for time-fractional advection-diffusion equations with smooth solutions. J. Comput. Phys. 298 (2015), 520-538. DOI 10.1016/j.jcp.2015.05.052 | MR 3374563 | Zbl 1349.65294
[15] G.-H. Gao, H.-W. Sun, Z.-Z. Sun: Stability and convergence of finite difference schemes for a class of time-fractional sub-diffusion equations based on certain superconvergence. J. Comput. Phys. 280 (2015), 510-528. DOI 10.1016/j.jcp.2014.09.033 | MR 3273149 | Zbl 1349.65295
[16] G.-H. Gao, Z.-Z. Sun: A compact finite difference scheme for the fractional sub-diffusion equations. J. Comput. Phys. 230 (2011), 586-595. DOI 10.1016/j.jcp.2010.10.007 | MR 2745445 | Zbl 1211.65112
[17] G.-H. Gao, Z.-Z. Sun, H.-W. Zhang: A new fractional numerical differentiation formula to approximate the Caputo fractional derivative and its applications. J. Comput. Phys. 259 (2014), 33-50. DOI 10.1016/j.jcp.2013.11.017 | MR 3148558 | Zbl 1349.65088
[18] A. K. Golmankhaneh, D. Baleanu: Calculus on fractals. Fractional Dynamics. De Gruyter, Berlin (2015), 307-332. DOI 10.1515/9783110472097-019 | MR 3468175 | Zbl 1333.00024
[19] D. He: An unconditionally stable spatial sixth-order CCD-ADI method for the two-dimensional linear telegraph equation. Numer. Algorithms 72 (2016), 1103-1117. DOI 10.1007/s11075-015-0082-7 | MR 3529835 | Zbl 1350.65088
[20] D. He, K. Pan: A fifth-order combined compact difference scheme for the Stokes flow on polar geometries. East Asian J. Appl. Math. 7 (2017), 714-727. DOI 10.4208/eajam.200816.300517a | MR 3766403 | Zbl 1426.76471
[21] D. He, K. Pan: An unconditionally stable linearized CCD-ADI method for generalized nonlinear Schrödinger equations with variable coefficients in two and three dimensions. Comput. Math. Appl. 73 (2017), 2360-2374. DOI 10.1016/j.camwa.2017.04.009 | MR 3648018 | Zbl 1373.65056
[22] T. N. Jones, Q. Sheng: Asymptotic stability of a dual-scale compact method for approximating highly oscillatory Helmholtz solutions. J. Comput. Phys. 392 (2019), 403-418. DOI 10.1016/j.jcp.2019.04.046 | MR 3948727
[23] N. A. Khan, M. Jamil, A. Ara: Approximate solutions to time-fractional Schrödinger equation via homotopy analysis method. Int. Sch. Res. Not. 2012 (2012), Article ID 197068, 11 pages. DOI 10.5402/2012/197068
[24] R. Klages, G. Radons, I. M. Sokolov (eds.): Anomalous Transport: Foundations and Applications. Wiley, Weinheim (2008). DOI 10.1002/9783527622979
[25] N. Laskin: Fractional quantum mechanics and Lévy path integrals. Phys. Lett., A 268 (2000), 298-305. DOI 10.1016/S0375-9601(00)00201-2 | MR 1755089 | Zbl 0948.81595
[26] S. T. Lee, J. Liu, H.-W. Sun: Combined compact difference scheme for linear second-order partial differential equations with mixed derivative. J. Comput. Appl. Math. 264 (2014), 23-37. DOI 10.1016/j.cam.2014.01.004 | MR 3164100 | Zbl 1294.65098
[27] L. Z. Li, H.-W. Sun, S.-C. Tam: A spatial sixth-order alternating direction implicit method for two-dimensional cubic nonlinear Schrödinger equations. Comput. Phys. Commun. 187 (2015), 38-48. DOI 10.1016/j.cpc.2014.10.008 | MR 3283134 | Zbl 1348.35238
[28] D. Li, J. Wang, J. Zhang: Unconditionally convergent $L1$-Galerkin FEMs for nonlinear time-fractional Schrödinger equations. SIAM. J. Sci. Comput. 39 (2017), A3067-A3088. DOI 10.1137/16M1105700 | MR 3738320 | Zbl 1379.65079
[29] H.-L. Liao, H.-S. Shi, Y. Zhao: Numerical study of fourth-order linearized compact schemes for generalized NLS equations. Comput. Phys. Commun. 185 (2014), 2240-2249. DOI 10.1016/j.cpc.2014.05.002 | MR 3211950 | Zbl 1344.35137
[30] A. Mohebbi, M. Abbaszadeh, M. Dehghan: The use of a meshless technique based on collocation and radial basis functions for solving the time fractional nonlinear Schrödinger equation arising in quantum mechanics. Eng. Anal. Bound. Elem. 37 (2013), 475-485. DOI 10.1016/j.enganabound.2012.12.002 | MR 3018826 | Zbl 1352.65397
[31] D. A. Murio: Implicit finite difference approximation for time fractional diffusion equations. Comput. Math. Appl. 56 (2008), 1138-1145. DOI 10.1016/j.camwa.2008.02.015 | MR 2437884 | Zbl 1155.65372
[32] M. Naber: Time fractional Schrödinger equation. J. Math. Phys. 45 (2004), 3339-3352. DOI 10.1063/1.1769611 | MR 2077514 | Zbl 1071.81035
[33] S. Z. Rida, H. M. El-Sherbiny, A. A. M. Arafa: On the solution of the fractional nonlinear Schrödinger equation. Phys. Lett., A 372 (2008), 553-558. DOI 10.1016/j.physleta.2007.06.071 | MR 2378723 | Zbl 1217.81068
[34] E. Scalas, R. Gorenflo, F. Mainardi: Fractional calculus and continuous-time finance. Phys. A 284 (2000), 376-384. DOI 10.1016/S0378-4371(00)00255-7 | MR 1773804
[35] E. Shivanian, A. Jafarabadi: Error and stability analysis of numerical solution for the time fractional nonlinear Schrödinger equation on scattered data of general-shaped domains. Numer. Methods Partial Differ. Equations 33 (2017), 1043-1069. DOI 10.1002/num.22126 | MR 3652177 | Zbl 1379.65082
[36] H.-W. Sun, L. Z. Li: A CCD-ADI method for unsteady convection-diffusion equations. Comput. Phys. Commun. 185 (2014), 790-797. DOI 10.1016/j.cpc.2013.11.009 | MR 3161142 | Zbl 1360.35193
[37] Z.-Z. Sun, X. Wu: A fully discrete difference scheme for a diffusion-wave system. Appl. Numer. Math. 56 (2006), 193-209. DOI 10.1016/j.apnum.2005.03.003 | MR 2200938 | Zbl 1094.65083
[38] Y.-M. Wang, L. Ren: Efficient compact finite difference methods for a class of time-fractional convection-reaction-diffusion equations with variable coefficients. Int. J. Comput. Math. 96 (2019), 264-297. DOI 10.1080/00207160.2018.1437262 | MR 3893479
[39] Z. Wang, S. Vong: Compact difference schemes for the modified anomalous fractional sub-diffusion equation and the fractional diffusion-wave equation. J. Comput. Phys. 277 (2014), 1-15. DOI 10.1016/j.jcp.2014.08.012 | MR 3254221 | Zbl 1349.65348
[40] L. Wei, Y. He, X. Zhang, S. Wang: Analysis of an implicit fully discrete local discontinuous Galerkin method for the time-fractional Schrödinger equation. Finite Elem. Anal. Des. 59 (2012), 28-34. DOI 10.1016/j.finel.2012.03.008 | MR 2949417
[41] Y. Xu, L. Zhang: Alternating direction implicit method for solving two-dimensional cubic nonlinear Schrödinger equation. Comput. Phys. Commun. 183 (2012), 1082-1093. DOI 10.1016/j.cpc.2012.01.006 | MR 2888993 | Zbl 1277.65073
[42] L. Zhu, Q. Sheng: A note on the adaptive numerical solution of a Riemann-Liouville space-fractional Kawarada problem. J. Comput. Appl. Math. 374 (2020), Article ID 112714, 14 pages. DOI 10.1016/j.cam.2020.112714 | MR 4066721 | Zbl 1435.65136

Affiliations:   Chun-Hua Zhang, Department of Mathematics, University of Macau, Avenida da Universidade, Taipa, Macau, China, e-mail: chzringlang@163.com; Jun-Wei Jin, Key Laboratory of Grain Information Processing and Control, Ministry of Education of China, and College of Information Science and Engineering, Henan University of Technology, Zhengzhou 450001, China, e-mail: jinjunwei24@163.com; Hai-Wei Sun, Department of Mathematics, University of Macau, Avenida da Universidade, Taipa, Macau, China, e-mail: hsun@um.edu.mo; Qin Sheng (corresponding author), Department of Mathematics and Center for Astrophysics, Space Physics and Engineering Research, Baylor University, Waco, TX 76798-7328, USA, e-mail: Qin_Sheng@baylor.edu


 
PDF available at: