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

来源 :东南大学 | 被引量 : 0次 | 上传用户:JK0803_tangkai
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
分数阶微积分最近十多年来越来越多地引起人们的关注.它可以成功地描述工程、物理、化学、生物以及经济学等领域的许多现象.由于分数阶导数固有的非局部属性,在最近几年,它已被证明可以用来描述一些与记忆和遗传有关的现象或过程.这些重要的应用促使我们努力的寻求高效、稳定并且易于执行的算法来求解分数阶微分方程.本文用有限差分方法研究了时间空间分数阶Bloch-Torrey方程的高精度的数值解法,并证明所建立的算法的稳定性和收敛性.文章分为两大部分.第一部分,研究了一维时间空间分数阶Bloch-Torrey方程的两个高阶数值算法,并给出相应的先验估计式.首先,利用加权位移的Grunwald-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格式虽然可以使问题得到简化,但时间方向的精度分别达到α阶或2a阶,比高阶差分格式C和D的精度降低了.因此,我们只给出了ADI求解格式,在理论上没有进一步的研究.
其他文献
今年以来,内蒙古自治区赤峰市将全国禁毒示范城市创建活动作为整体提升禁毒工作的抓手,坚持建机制、固基础,按照预防宣传要“新”、动态管控要“细”、禁吸收戒要“准”、戒
报纸
风湿寒性关节痛是原因未明,以风、寒、湿环境因素及体力负荷为主要危险因子,以膝、腰椎等多关节疼痛为主要临床表现的一种良性风湿性综合征.本文从风关痛的流行病学、基础研
<正>高血压是全球主要的公共卫生问题之一。我国成人高血压患者约有1.6亿[1],而其中8%~15%为顽固性高血压。顽固性高血压的定义:经使用包括利尿剂在内的、足量且合理的3种或以
<正>抗利尿激素分泌异常综合征(syndrome of inappropriate antidiuretic hormone secretion,SIADH)由Schwartz、Bartter等[1]于1957年首次报道,故又称Schwartz-Bartter综合
目的分析宁夏回族人群线粒体DNA编码区10394、5176、4831、7598、13259、663等6个位点基因型及等位基因频率的分布特征。方法采用PCR-RFLP技术,分析线粒体DNA编码区6个位点及
由于非洲经济的特殊性,我国对外承包工程所需要.的大部分物资严重依赖进口。因此,做好物资抵港后的}快速清关工作对于节约时间与财力成本颇为重要。
伴随着我国广播电视以及广播电台技术的普及以及节目的多样化。我国的广播电台已经成为我国信息传播的重要途径,作为我国媒体传播途径的重要组成部分之一,我国的电台广播有着
动态时间规整(DynamicTimeWarping)是语音识别中的一种经典算法,由于算法简单又有效,因此在实现孤立词识别系统中获得了广泛的应用,但仍然存在一些不足,本文从提高识别率和识
判断反腐败斗争的历史转型要兼顾反腐败要素上的变化和结构上的调整两个方面。从中心要素考察,反腐败的基本类型有自发反腐、自觉反腐和自律反腐。接续不断的中国反腐败史经
目的分析研究慢性肺心病并发心衰患者采取急诊护理措施的意义。方法特选我院2018年诊治并确诊为慢性肺心病并发心衰患者82例,根据护理方法以及随机数表法将这82例分别分为对