基于耦合纯无网格方法时间分数阶下孤立子波碰撞过程的数值模拟研究

Numerical Simulation of the Solitary Wave Collision Process in Time Fractional Orders Based on the Coupled Pure Meshless Method

为数值预测时间分数阶耦合非线性Schrödinger(TF-CNLS)方程描述的孤立子波非弹性碰撞过程,首次发展了一种耦合纯无网格有限点集法(coupled finite pointset method,CFPM).其构造过程为:1)对时间分数阶Caputo导数项采用一种高精度的差分格式;2)对空间导数采用基于Taylor展开和加权最小二乘法的有限粒子法(FPM)离散格式;3)对区域进行局部加密和采用稳定性好的双曲余弦核函数以提高数值精度.数值研究中,首先,运用CFPM对有解析解的一维TF-CNLS方程进行求解,分析了节点均匀分布或局部加密情况下的误差和收敛阶,表明给出的耦合无网格法具有近似二阶精度和易局部加密求解的灵活性;其次,运用CFPM对无解析解一维TF-CNLS方程描述的孤立子波非弹性碰撞过程进行了数值预测,其出现的波塌缩现象与整数阶下出现的多波现象截然不同;最后,与有限差分结果作对比,表明CFPM数值预测时间分数阶下孤立子波非弹性碰撞过程的复杂传播现象是可靠的.

A coupled pure meshless finite pointset method(CFPM)was developed for the first time to numerically predict the inelastic collision process of solitary waves described with the time fractional coupled nonlinear Schrödinger(TF-CNLS)equation.Its construction process was formulated as:1)a high-precision difference scheme was used for the Caputo time fractional derivative;2)the FPM discrete scheme based on the Taylor expansion and the weighted least square method was adopted for spatial derivatives;3)the region was locally refined and the double cosine kernel function with good stability was used to improve the numerical accuracy.In the numerical study,the 1D TF-CNLS equations were analytically solved with the CFPM,and the errors and convergence rates were analyzed with the nodes uniformly distributed or locally refined,which shows that the proposed method has the approximate 2nd-order accuracy and the flexibility of easy local refinement.Secondly,the inelastic collision process of solitary waves,which was described with the 1D TF-CNLS equation without analytical solutions,was numerically predicted with the CFPM,and the wave collapse phenomenon is completely different from the multi-wave phenomenon in the integer order.Meanwhile,the comparison of the results with those from the finite difference method shows that,the CFPM is reliable to predict the complex propagation of the inelastic collision process of the solitary waves in the time fractional order.