摘要
在提出一种基于时间分裂格式的纯无网格有限点集(split-step finite pointset method, SS-FPM)法的基础上,数值模拟了含孤立波的二维非线性薛定谔(nonlinear Schr?dinger, NLS)/(Gross-Pitaevskii, GP)方程.SS-FPM的构造过程为:1)基于时间分裂的思想将非线性薛定谔方程分成线性导数项和非线性项; 2)采用基于Taylor展开和加权最小二乘法的有限点集法,借助Wendland权函数,对线性导数项进行数值离散.随后,模拟了带有Dirichlet和周期性边界条件的NLS方程,将所得结果与解析解做对比.数值结果表明:给出的SS-FPM粒子法的优点是在粒子分布非均匀情况下仍具有近似二阶精度,且较网格类有限差分算法实施容易,较已有改进的光滑粒子动力学方法计算误差小.最后,运用SS-FPM对无解析解的二维周期性边界NLS方程和Dirichlet边界玻色-爱因斯坦凝聚二分量GP方程进行了数值预测,并与其他数值结果进行对比,准确展现了非线性孤立波奇异性现象和量子化涡旋过程.
In this paper, a split-step finite pointset method(SS-FPM) is proposed and applied to the simulation of the nonlinear Schr?dinger/Gross-Pitaevskii equation(NLSE/GPE) with solitary wave solution. The motivation and main idea of SS-FPMisas follows. 1) The nonlinear Schr?dinger equation is first divided into the linear derivative term and the nonlinear term based on the time-splitting method. 2) The finite pointset method(FPM) based on Taylor expansion and weighted least square method is adopted, and the linear derivative term is numerically discretized with the help of Wendland weight function. Then the two-dimensional(2D) nonlinear Schr?dinger equation with Dirichlet and periodic boundary conditions is simulated, and the numerical solution is compared with the analytical one. The numerical results show that the presented SS-FPM has second-order accuracy even if in the case of non-uniform particle distribution, and is easily implemented compared with the FDM, and its computational error is smaller than those in the existed corrected SPH methods. Finally, the 2D NLS equation with periodic boundary and the two-component GP equation with Dirichlet boundary and outer rotation BEC, neither of which has an analytical solution, are numerically predicted by the proposed SS-FPM.Compared with other numerical results, our numerical results show that the SS-FPM can accurately display the nonlinear solitary wave singularity phenomenon and quantized vortex process.
引文
[1]Bao W Z,Chern I L,Lim F Y 2006 J.Comput. Phys.219 836
[2]Qu C,Sun K,Zhang C 2015 Phys.Rev.A 91 053630
[3]Mason P,Aftalion A 2011 Phys.Rev.A 84 033611
[4]Antoine X,Bao W,Besse C 2013 Comput.Phys.Commun.184 2621
[5]Wang D S,Xue Y S,Zhang Z F 2016 Rom.J.Phys.61 827
[6]Wang D S,Shi Y R,Feng W X,Wen L 2017 Physica D351-352 30
[7]Wang H 2005 Appl. Math.Comput.170 17
[8]Gao Y L,Mei L Q 2016 Appl. Num.Math.109 41
[9]Blanes S,Casas F,Murua A 2015 J.Comput.Phys.303 396
[10]Dehghan M,Taleei A 2010 Comput.Phys.Commun.181 43
[11]Wang T C,Guo B L,Xu Q B 2013 J.Comput.Phys.243 382
[12]Chen R Y,Pan W L,Zhang J Q,Nie L R 2016 Chaos 26093113
[13]Chen R Y,Tong L M,Nie L R,Wang C I,Pan W 2017Physica A:Statist. Mech.Appl. 468 532
[14]Chen R Y,Nie L R,Chen C Y 2018 Chaos 28 053115
[15]Gong Y Z,Wang Q,Wang Y S,Cai J X 2017 J.Com.put.Phys.328 354
[16]Cheng R J,Cheng Y M 2016 Chin.Phys.B 25 020203
[17]Dehghan M,Mirzaei D 2008 Int.J.Numer.Meth.76 501
[18]Abbasbandy S,Roohani Ghehsareh H,Hashim I 2013 Eng.Anal. Bound.Elem.37 885
[19]Liu M B,Liu G R 2010 Arch.Comput. Meth.Eng.17 25
[20]Liu M B,Chang J Z 2010 Acta Phys.Sin.59 3654(in Chinese)[刘谋斌,常建忠2010物理学报59 3654]
[21]Huang C,Lei J M,Liu M B,Peng X Y 2015 In.J.Num.Meth.Flu.78 691
[22]Jiang T,Chen Z C,Ren J L,Li G 2017 Acta Phys.Sin.66130201(in Chinese)[蒋涛,陈振超,任金莲,李刚2017物理学报66 130201]
[23]Jiang T,Chen Z C,Lu W G,Yuan J Y,Wang D S 2018Comput.Phys.Commun.231 19
[24]Kuhnert J,Tiwari S 2001 Beri.chte des Fraunhofer ITWMNr.25
[25]Kuhnert J,Tiwari S 2001 Berichte des Fraunhofer ITWMNr.30
[26]Resendiz-Flores E O,Garcia-Calvillo I D 2014 Int. J.Heat Mass Trans.71 720
[27]Wendland H 1995 Adv.Comput.Math.4 389