论文部分内容阅读
反常扩散现象在自然界中普遍存在,尤其在某些复杂系统的扩散过程中该现象更为常见.为了更好地解释这一现象,不同的学者提出了不同的工具和理论.过去的二十年,许多研究人员发现分数阶微积分可以更为精确地描述这一现象,并且建立了大量的描述复杂系统中反常扩散输运过程的分数阶反常扩散方程.尽管这些分数阶反常扩散方程更精确地解释了所研究的实际问题,但求得这些方程的解析解是比较困难的.因此,寻求此类方程的数值解是十分有意义的.本文主要研究了描述反常扩散的时间分数阶(径向)扩散方程、分数阶Klein-Kramers方程、分数阶守恒律的数值求解问题,设计了这几类方程的稳定、有效的数值格式,建立了所给格式的误差估计,研究了这些方程的动力学行为.本文的主要工作包括以下三部分:第一部分考虑了描述次扩散的两类时间分数阶扩散方程的数值求解问题.首先,设计了时间分数阶扩散方程的两种高阶、容易实施、无条件稳定的正交样条配置格式.然后,给出了求解时间分数阶径向扩散方程的两种隐式差分格式,利用数学归纳法和离散极值原理证明了所给数值格式均为无条件稳定的.并用数值结果和数值模拟验证了理论分析的收敛阶和所给格式的有效性.第二部分讨论了两类分数阶Klein-Kramers方程的数值解.首先,给出了求解时间分数阶Klein-Kramers方程的有限差分格式,给出了数值格式的稳定性和收敛性的严格证明,随后的数值结果验证了理论分析的正确性.进一步,讨论了Levy分数阶Klein-Kramers方程的数值解,建立了Levy分数阶Klein-Kramers方程的隐式和显式的有限差分格式,借助广义的离散极值原理分析了数值格式的稳定性和收敛性,同时给出了几种提高精度的技巧,且数值结果表明所给方法是有效的.第三部分讨论了分数阶守恒律的数值算法.首先,给出了具有光滑解的周期分数阶非线性守恒律的半离散Fourier谱方法,详细讨论了半离散格式的误差分析,并利用四阶的积分因子-Runge-Kutta方法来求解半离散后的方程组.数值结果表明该方法空间上可以达到谱精度且时间上为四阶收敛的.接着,对于分数阶守恒律的非光滑初值问题,设计了该方程的分数步方法,数值结果表明所给数值算法对于光滑和间断的初值均是有效的.