时间空间分数阶BlocH-Torrey方程的高阶数值解法

来源 :东南大学 | 被引量 : 0次 | 上传用户:ya4516623
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
分数阶微积分最近十多年来越来越多地引起人们的关注.它可以成功地描述工程、物理、化学、生物以及经济学等领域的许多现象.由于分数阶导数固有的非局部属性,在最近几年,它己被证明可以用来描述一些与记忆和遗传有关的现象或过程.这些重要的应用促使我们努力的寻求高效、稳定并且易于执行的算法来求解分数阶微分方程.本文用有限差分方法研究了时间空间分数阶Bloch-Torrey方程的高精度的数值解法,并证明所建立的算法的稳定性和收敛性.  文章分为两大部分.  第一部分,研究了一维时间空间分数阶Bloch-Torrey方程的两个高阶数值算法,并给出相应的先验估计式.首先,利用加权位移的Grünwald-Letnikov算子来离散Caputo分数阶导数,得到一个关于时间分数阶导数的三阶逼近;对于空间方向的分数阶导数,利用中心差商公式来离散Riesz分数阶导数,得到了空间方向的二阶逼近.这样对方程构造了高阶差分格式A,使其时间方向和空间方向在L1(L2)范数下分别达到三阶和二阶精度.利用离散的能量分析方法,严格证明了该差分格式的唯一可解性、无条件稳定性和收敛性,并且给出了数值算例验证了数值解的精确性和差分格式的有效性.然后,对于时间方向的分数阶导数仍采用刚才提到的离散方法来离散;对于空间方向的分数阶导数,用中心差商公式来离散Riesz分数阶导数在三点的加权值,得到空间方向的一个四阶逼近.这样对方程构造了高阶差分格式B,其在时间和空间方向达到三阶和四阶精度.同样,我们利用离散的能量分析方法,严格证明了高阶差分格式B的唯一可解性、无条件稳定性和收敛性,并且给出了数值算例验证了数值解的精确性和差分格式的有效性.第二部分,研究了二维时空分数阶Bloch-Torrey方程的两个高阶的数值解法,并给出相应的先验估计式;然后给出四个交替方向的差分格式.首先,对于时间方向和空间方向的分数阶导数,采用一维情形时得到差分格式A时用到的离散技巧,分别利用三阶逼近公式来离散Caputo分数阶导数、利用二阶逼近公式来离散Riesz分数阶导数;这样对二维问题构造了高阶差分格式C,使其在时间方向和空间方向在L1(L2)范数下分别达到三阶精度和二阶精度.利用离散能量分析方法,严格证明了所建立的高阶差分格式C的唯一可解性、无条件稳定性和收敛性,并且给出了数值算例验证了数值解的精确性和差分格式的有效性.然后,对空间方向的分数阶导数的离散,采用一维情形时得到差分格式B时用到的离散技巧,得到了空间方向的四阶逼近.这样对二维问题构造了高阶差分格式D,使其在时间方向和空间方向在L1(L2)范数下分别达到三阶精度和四阶精度.利用离散能量分析方法,严格证明了高阶差分格式D的唯一可解性、无条件稳定性和收敛性,并且给出了数值算例验证了数值解的精确性和差分格式的有效性.最后对二维时空分数阶Bloch-Torrey方程建立了4个ADI求解格式.ADI格式虽然可以使问题得到简化,但时间方向的精度分别达到α阶或2α阶,比高阶差分格式C和D的精度降低了.因此,我们只给出了ADI求解格式,在理论上没有进一步的研究.
其他文献
本文主要讨论了在求解自然电位测井的正问题时所用到的一个数值方法。在第一章中,我们首先对自然电位测井的原理进行了阐述,并对地层信息进行合理的假设,进而建立相应的数学
介绍了某核电厂制氯站循环水处理系统的工艺流程,PLC控制系统及过程检测仪表的设计要求,PLC硬件配置及仪表选型安装中应注意的问题等,系统采用AB 1756系列PLC,CPU及通讯接口
角点检测在计算机视觉与模式识别学科领域有着重要的意义与广泛的应用。其中,多尺度角点检测既能够完整描述信号的整体特征又不至于忽视信号的细节特征,因此受到人们的广泛关
一、注意力和注意力营销 1.什么是注意力和注意力营销注意力营销,就是以获取消费者注意为目的的营销。要搞清楚什么是注意力营销,就要先弄明白什么是注意力。注意力是表达注