论文部分内容阅读
本文主要利用两类二阶时间逼近格式结合有限元方法数值求解含有时间分数阶导数的非线性Cable方程和含有时空分数阶导数的一类非线性偏微分方程.在此前的研究中,对于含有分数阶导数的Cable方程,通常采用L1逼近离散Riemann-Liouville分数阶导数,并得到2-α阶的近似精度.在本文中,我们利用具有2阶近似精度的一类加权离散公式逼近Riemann-Liouville分数阶导数,进而形成二阶离散格式.另外我们利用有限元方法结合WSGD算子的离散技巧数值求解一类时空分数阶偏微分方程,并得到时间二阶逼近结果. 首先,对于非线性时间分数阶Cable方程,在tk-α/2时间点处,我们用一个带有参数α的两步格式来近似整数阶时间导数,采用具有2阶精度的加权离散格式来近似Riemann-Liouville分数阶导数,并利用Galerkin有限元法进行空间方向离散,形成时间二阶逼近的有限元数值方法.进一步,我们给出了一些数值理论,证明了格式的稳定性,并得到基于L2模意义下的带有O(△t2+(1+△t-α)hm+1)精度的误差估计结果.为了说明理论结果的正确性,我们通过几个算例给出数值验证. 其次,利用二阶WSGD逼近结合有限元数值方法数值求解含有时空分数阶导数的一类偏微分方程,给出了详细的数值计算过程,并通过三类数值例子对格式进行有效验证.数值结果显示我们的数值格式具有二阶近似精度.