非线性变阶分数阶扩散方程的全隐差分格式

对于变阶的非线性分数阶扩散方程,提出了一种全隐的差分格式。然后,通过离散的能量方法证明了所提出的格式是无条件稳定的,其收敛阶为O(τ+h)。通过数值试验表明,全隐的差分格式是有效的和可靠的。     

二维空间分数阶变系数扩散方程的数值解法

研究二维有限域上的扩散系数与空间变量相关的空间分数阶扩散方程,通过移位的Grunwald公式对空间分数阶导数进行离散,得到交替差分格式,证明了格式的稳定性,最后给出了数值算例.     

时间分数阶扩散波动方程的二阶有限差分格式

基于时间分数阶扩散波动方程的等价积分形式,采用分数阶梯形法和Crank-Nicolson方法,对时间分数阶扩散波动方程初边值问题设计了一个计算稳定的有限差分格式,此格式在时间方向和空间方向都具有二阶精度。数值算例验证了该格式的精度和效果。     

一类分数阶对流扩散方程差分格式的理论分析

考虑一般的对流扩散方程,将一阶的时间导数用Caputo分数阶导数替换,二阶的空间导数用Riemann-Liouville分数阶导数替换,得到了一个Riemann-Liouville-Caputo分数阶对流扩散方程.给出了这个方程的一种计算有效的隐式差分格式,并证明了该差分格式是无条件稳定、无条件收敛的,其收敛阶为O(l+h).最后给出了数值例子.     

时间分数阶扩散方程的隐式差分近似

考虑时间分数阶扩散方程,它是从标准的扩散方程中用分数阶导数α(0α1)代替一阶时间导数而得到,提出了一个计算有效的隐式差分近似,并证明了这个隐式差分近似是无条件稳定和无条件收敛的。最后给出了数值例子。     

两边空间-时间分数阶扩散方程的加权有限差分格式

对于空间-时间分数阶扩散方程的初边值问题提出了一种加权差分格式. 利用能量估计, 得到了差分格式的稳定性. 然后使用数学归纳法证明了在相同的条件下, 所提出的的格式是收敛的. 最后通过一个例子说明了所提出的格式是可靠的、有效的.     

一种时间分数阶对流扩散方程的隐式差分近似

考虑时间分数阶对流扩散方程时,将一阶的时间导数用分数阶导数α(0<α<1)替换,给出一种计算有效的隐式差分格式,并证明这个隐式差分格式是无条件稳定、无条件收敛的.最后用数值例子说明差分格式是有效的.     

二维时间分数阶扩散方程Wilson元的收敛性分析

利用Wilson元提出了一类二维时间分数阶扩散方程的新的全离散逼近格式.基于单元的性质,在不需要外推和插值后处理技术的前提下,得到了u的比传统的H1-范数更大的模意义下相应的O(h~2+τ~(2-α/2))阶的误差分析结果,正好比通常的关于Wilson元的误差估计高出一阶.这里,h,τ分别表示空间剖分参数和时间步长.     

分数阶对流扩散方程的半加权有限差分格式（英文）

对于空间分数阶对流扩散方程的初边值问题提出了一系列半加权差分格式.可以证明此格式当分数阶导数属于[((17)~(1/2)-1)/2,2]时无条件稳定,且二阶收敛.最后给出数值算例验证了理论证明.     

时间分数阶对流-扩散方程的有限差分方法

研究时间分数阶常系数对流-扩散方程的数值解,提出了一种只需要存储部分历史数据的分数阶微分方程的数值计算方法,并给出了误差估计.     

分数阶对流-弥散方程的有限差分方法

本文对分数阶对流-弥散方程的初边值问题进行了数值研究.我们采用移位Grun-wald公式对空间分数阶导数进行离散,在此基础上建立Crank-Nichonlson(简称C-N)差分格式,并讨论了差分解的存在唯一性,然后分析了该方法的稳定性及收敛性,并利用外推法提高收敛阶.数值算例验证了格式的有效性.     

时间分数阶对流-扩散方程的有限差分法

时间分数阶对流-扩散方程可以用来模拟由传统的对流-扩散方程演变而来的反常扩散方程.本文针对一类时间分数阶对流-扩散方程提出了一个新的隐式差分格式,时间分数阶导数采用直接离散,空间导数采用中心差分格式离散,讨论了差分解的存在唯一性,并利用能量范数证明了该格式的无条件稳定性、收敛性,分析了收敛阶.数值试验验证了该格式的有效性.     

两类分数阶对流-扩散方程的有限差分方法

考虑两类分数阶偏微分方程,空间分数阶对流-扩散方程和时间-空间分数阶对流-扩散方程。基于移位的Grünwald公式,在第一类方程中,空间分数阶导数用加权平均有限差分法来近似,用特征值方法给出了稳定性分析,误差估计为O(τ+h);在第二类方程中,时间导数逼近用高阶近似,根据最大模估计方法证明了稳定性,其收敛阶为O(τ2-max{γ1,γ2}+h),这里γ1,γ2分别是方程中出现的两项Caputo时间分数阶导数的阶。数值实例验证了理论结果。     

时间-空间分数阶扩散方程求解的新方法

介绍了3种求解带有Caputo型导数的时间-空间分数阶扩散方程的方法.通过分离变量和级数展开求数值解,将Fourier变换和Laplace变换用于求解析解,并把时间和空间定义域上的分数阶导数分别限制在0γ≤1,0β≤2.     

分数阶随机微分方程的修正隐式数值格式

对H>1/2且随机积分为前向积分的分数阶布朗运动驱动的随机微分方程,为改善显式Euler格式和Milstein格式的稳定性,基于修正隐式技术构造了修正隐式Euler格式和Milstein格式,证明了修正隐式格式较显式格式有更大的稳定步长集,且在一定条件下修正隐式Euler格式是A稳定的.数值模拟显示,步长在稳定步长集内时数值格式稳定,步长在稳定步长集边界附近时误差几乎不改变,而步长在稳定步长集外时数值格式极度不稳定,从而验证了修正隐式格式在保持数值稳定性上的优越性和稳定步长集的合理性.     

二维扩散反应方程的三次样条高精度加权隐式差分格式及其多重网格方法

利用二阶微商的三次样条四阶紧致差分逼近公式,推导出两种数值求解二维扩散反应方程的两层9点加权隐式紧致差分格式.当θ=1/2时,该格式在时间和空间方向上分别达到二阶和四阶精度.通过Fourier方法讨论知,当1/2≤θ≤1时,格式是无条件稳定的;当0≤θ<1/2时,格式是条件稳定的.为了克服传统迭代法在求解隐格式方面的困难,差分方程采用多重网格方法进行求解并将本文格式的结果与P-R格式及C-N格式下的结果进行比较.数值实验结果验证本文方法的精确性和可靠性及多重网格方法的效率.