首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
介绍求解多项四阶时间分数阶慢扩散方程的有限差分方法.利用L1公式逼近时间分数阶导数,用降阶法处理空间四阶导数项,再借助离散能量方法证明差分格式是无条件稳定的且在无穷范数下其收敛阶为O(τ2-β+h2),其中τ和h分别为时间方向和空间方向的步长,β是时间分数导数的最大阶.最后用数值实验验证所提出差分格式的精度和有效性.  相似文献   

2.
介绍求解多项四阶时间分数阶慢扩散方程的有限差分方法.利用L1公式逼近时间分数阶导数,用降阶法处理空间四阶导数项,再借助离散能量方法证明差分格式是无条件稳定的且在无穷范数下其收敛阶为O(t^(2-B)+h^2),其中t和h分别为时间方向和空间方向的步长,B是时间分数导数的最大阶.最后用数值实验验证所提出差分格式的精度和有效性.  相似文献   

3.
变系数分数阶反应-扩散方程的数值解法   总被引:1,自引:0,他引:1  
考虑了变系数分数阶反应一扩散方程,将一阶的时间偏导数和二阶的空间偏导数分别用Caputo分数阶导数和Riemann-Liouville分数阶导数替换,利用L1算法和G算法对方程的变系数分数阶导数进行适当的离散,给出了该方程的一种计算有效的隐式差分格式,并证明了这个差分格式是无条件稳定和无条件收敛的,且具有o(τ+h)收敛阶.最后用数值例子说明差分格式是有效的.  相似文献   

4.
针对一维波动方程提出了一种有限差分方法.首先,采用泰勒级数展开公式和原方程代入的方法推导出了第一个时间层未知函数值的四阶紧致差分格式.然后,用四阶紧致差分公式近似空间导数项,采用中心差分格式截断误差余项修正的方法处理时间导数项,推导出了第二个时间层以后未知函数的四阶紧致差分格式.该方法时间和空间具有整体四阶精度.利用Fourier方法分析了所提格式的稳定性.由于本文格式在未知时间层仅涉及3个网格点,因此可采用追赶法求解离散化后所得到的线性方程组.最后,用数值算例验证了本文格式的精确性和稳定性.  相似文献   

5.
针对非线性变阶空间-时间分数阶对流-扩散方程的初边值问题,提出一种全隐式有限差分格式.首先,分别对Riemann-Liouville型变时间分数阶导数算子和Riemann-Liouville型变空间分数阶导数算子和广义Riesz分数阶导数算子进行离散化处理;然后,通过离散的能量方法证明全隐式有限差分格式的稳定性和收敛性,并验证其收敛阶为O(τ+h);最后,通过数值算例检验该方法.试验结果表明:全隐式有限差分格式求解非线性变阶空间-时间分数阶对流-扩散方程初边值问题是可行和有效的.  相似文献   

6.
作者在已有的的具有非局部边值条件的抛物型方程的一个高精度差分格式的基础上,通过改变初值的近似公式,利用不同的能量估计方法,去掉了对边值积分核函数的限制条件,证明了所给差分格式的无条件可解性和L∞范数收敛性.此差分格式的收敛速率关于时间步长是二阶的、关于空间步长是四阶的.数值例子用以说明收敛性.  相似文献   

7.
提出了求解有限区域上的一维时空分数阶变系数对流扩散方程的两种隐式有限差分格式,就格式的精度和收敛阶比较这两种差分格式的优劣.当使用Caputo分数阶导数对a阶时间导数项进行离散时,在两个不同的点上分别采用中心差分,而对β阶空间导数项均使用转化的Grünwald公式进行离散.对得到的两种格式进行稳定性和收敛性分析.用几个已知精确解的数值例子验证和比较这两种有限差分格式的精确性和有效性.  相似文献   

8.
该文提出了在周期和Dirichlet边界条件下的1维对流扩散方程的紧致差分格式.在这2种边界条件下对空间变量使用4阶紧致差分格式,对时间变量利用3次Hermite插值公式构造空间和时间同时具有4阶精度的数值格式,并证明了格式的绝对稳定性,最后通过对2种边界条件下的算例进行数值实验和比较,验证了格式的精确性和可靠性.  相似文献   

9.
郭非凡  张新东  王硕 《河南科学》2019,37(6):878-886
近些年,越来越多的研究表明,随着时间或者空间变化,方程的扩散系数也会改变.主要研究了变系数分数阶扩散方程的有限差分解法.首先,引入半整数点,在空间网格上进行对偶剖分,再通过差分方法离散空间二阶偏导数.其次,利用两种分数阶导数,即Grünwald-Letnikov导数与Caputo导数的关系,近似替代时间分数阶导数,从而得到了收敛精度为o(t+h~2)的有限差分格式,并且该有限差分格式的解是存在且唯一的.最后,通过利用数学归纳法和最大模方法,证明出差分格式的稳定性和收敛性,并用一个一维时间分数阶变系数扩散方程的数值算例来验证差分格式的收敛阶.  相似文献   

10.
考虑了一类时间分数阶变系数时滞扩散微分方程,方程中对时间的导数由传统一阶导数变为α(0α1)阶导数,对此类方程利用差分法构造了有效的差分格式,并对此格式的收敛性和稳定性进行证明,数值算例检验该格式解决此类方程是有效的.  相似文献   

11.
对流扩散方程的研究大多局限于常扩散系数或整数阶的范围,为了能更加精确的描述溶质的运动特征,将它拓广到变扩散系数的情形,用Caputo型分数阶导数取代时间上的整数阶导数.对这种变系数时间分数阶对流扩散方程建立了一种隐式的有限差分格式,证明了该格式差分解的存在唯一性,分析了差分解的收敛性和稳定性,并用数值实验验证了此差分格式的有效性.  相似文献   

12.
提出一种带有Caputo导数的时间分数阶变系数扩散方程的数值解法.方程的解在初始时刻附近通常具有弱正则性,采用非一致网格上的L1公式离散时间分数阶导数,并使用局部间断Galerkin (local discontinuous Galerkin, LDG)方法离散空间导数,给出方程的全离散格式.基于离散的分数阶Gronwall不等式,证明了格式的数值稳定性和收敛性,且所得结果关于α是鲁棒的,即当α→1-时不会发生爆破.最后,通过数值算例验证理论分析的结果.  相似文献   

13.
刘明鼎  张艳敏 《河南科学》2014,(9):1688-1691
给出了求解一类时间分数阶时滞微分方程的数值解法,将传统对时间的一阶导数利用分数阶导数α(0α1)阶导数代替,给出了求解微分方程的差分格式,并对差分格式证明了收敛性和稳定性,数值算例检验该格式解决此类方程是有效的.  相似文献   

14.
研究时间分数阶扩散方程中分数导数阶的估计问题.首先,定义了一个带自然对数核的Caputo分数阶导数算子,推导出了时间分数阶扩散方程的分数阶导数所满足的方程,称之为时间分数阶扩散的伴随方程.其次,我们分别对两个方程进行时间离散构造有限差分格式和弱形式,再对弱形式中的半离散解进行Legendre多项式逼近得到全离散格式.然...  相似文献   

15.
二维变系数空间分数阶电报方程数值解   总被引:1,自引:0,他引:1  
针对二维变系数空间分数阶电报方程,利用Grünwald-Letnikov分数阶导数的定义,在交替方向法的基础上提出了一种分数阶Peaceman-Rachford差分格式.通过Gerschgorin定理和Lax等价定理证明了所提出的分数阶Peaceman-Rachford差分格式是无条件稳定和收敛的.数值试验表明:分数阶Peaceman-Rachford差分格式是有效和可靠的.  相似文献   

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

17.
提出一种求解Riesz空间分布阶的分数阶扩散方程的数值方法.利用辛普森数值求积公式,将分布阶微分方程离散为一个多项分数阶导数的微分方程;利用四阶差分格式求解此具有多项分数阶导数的微分方程,并运用能量法分析数值格式的稳定性和收敛性.同时,给出数值例子,说明所建立的数值离散格式的有效性.  相似文献   

18.
将一阶的时间偏导数用Coimbra变时间分数阶导数算子进行替换,提出了一种新隐式差分解法.首先,对Coimbra型变时间分数阶导数算子和二阶空间导数进行离散化处理,将Coimbra变时间分数阶扩散-波动方程转化为代数方程组求解;然后,借助于数学归纳法给出了新隐式差分方法的收敛性分析,并证明了新隐式差分方法是无条件收敛的;最后,通过数值例子检验该方法,计算结果表明新隐式差分方法的理论分析是正确的,所构造的离散格式是可行的和有效的.  相似文献   

19.
提出一个纳米尺度的分数阶抛物两步模型,得到金属纳米尺度热传导的精确数值格式.该模型是通过引入Caputo-Hadamard时间分数阶导数到抛物型两步能量输运方程中,并将其温度跃变边界条件耦合得到.数值格式基于空间四阶紧格式和Caputo-Hadamard时间分数阶导数的L1逼近格式而建立.通过2个算例验证模型和数值方法的准确性和适用性.  相似文献   

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

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号