基于高阶逼近公式的时间分数阶偏微分方程的有限元方法

来源 :内蒙古大学 | 被引量 : 0次 | 上传用户:alex851123
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文使用高阶逼近公式对分数阶导数进行了离散,并结合有限元方法对三类时间分数阶偏微分方程进行数值求解.第二章研究非线性分布阶亚扩散模型的基于时间两网格(TT-M)的1-Galerkin混合元方法.该算法相比1-Galerkin混合元方法能够提升计算效率,节省计算时间.时间方向采用快速的TT-M Crank-Nicolson格式逼近,其中时间分布阶导数采用FBN-结合复化梯形公式离散.空间方向采用1-Galerkin混合方法逼近.证明数值格式解的存在唯一性,给出稳定性和先验误差估计的详细分析.通过带有光滑解的算例验证算法的有效性和快速算法的计算效率,进一步通过非光滑解算例的数值计算结果说明修正格式能够有效地解决由弱奇异导致的掉精度问题.第三章采用一类全离散混合有限元方法数值求解非线性分布阶四阶亚扩散模型,其中时间方向采用分数阶导数的高阶FBT-逼近公式逼近,并采用TT-M快速算法加速计算,空间方向采用一种混合元方法逼近.给出理论误差估计的详细分析,通过带有光滑或非光滑解的数值算例验证理论结果的正确性,说明快速算法的计算效率和校正格式在提升计算精度方面的有效性.第四章提出变阶分数阶导数的一类高阶逼近公式,证明其具有2阶逼近精度.基于该逼近公式,给出变分数阶Allen-Cahn模型全离散有限元格式,通过算例检验数值方法的可行性和高阶逼近公式的效果.
其他文献
体细胞重编程技术在干细胞生物学,疾病建模,治疗性克隆等方面具有巨大的潜力,但是影响体细胞重编程效率的障碍以及不同体细胞重编程技术之间的差异仍未被完全阐明。本研究以牛和小鼠的体细胞核移植(somatic cell nuclear transfer,SCNT)和诱导多能性干细胞(induced pluripotent stem cells,iPSC)的重编程过程为研究对象,通过单细胞转录组测序和整合公
拓扑半金属在费米能级处的非平庸线性交叉可产生独特的拓扑表面态,有望在未来拓扑量子计算方面得以广泛应用。因此,近年来拓扑半金属材料的研究成为人们关注的热点。传统的拓扑半金属材料(topological semimetallic materials,TSM)通常由具有高强度自旋轨道耦合(spin-orbit coupling,SOC)作用的重元素构成。这些重元素往往价格昂贵、不环保,对材料的实际应用极
作为新型量子态,拓扑半金属具有诸多奇特的物理性质,含超高速、低能耗量子传输,费米弧表面态,高载流子迁移率,负磁阻等。当前,探索这种新奇的拓扑量子态并寻找相应的真实稳定材料已成为凝聚态物理的一个热点。随着新算法的发展,人们在已有的3D晶体数据库中发现了数千种拓扑半金属。然而,它们中的大多数由有毒重元素构成,这对实际应用非常不利。鉴于此,本文以环境友好的Ge元素为研究对象,采用第一性原理方法,探究了具
本文运用手征SU(3)夸克模型和扩展的手征SU(3)夸克模型,采用共振群方法分别对S波重子Σ(*)与介子DΣ(*)、B(*)和重子Σ(*)与反介子(?)(*)、(?)(*)组合而成五夸克态系统进行了动力学研究。首先,我们对S波重子Σ(*)与介子D(*)或B(*)组成的五夸克系统做了研究。我们先后计算了所有可能量子数的S波系统的束缚能、系统内部的哈密顿量矩阵元和系统内部各部分相互作用势的贡献。由数值
本文针对一维对流扩散反应方程和Burgers方程,根据方程各自的特点构造了相应的三次有限体积元数值格式,给出其收敛性分析,并利用数值算例验证理论分析结果和所提数值格式的有效性.具体研究内容如下.第一部分针对一维对流扩散反应方程,用最佳应力节点构建了对偶网格剖分,采用基于分片三次Lagrange插值的试探函数空间和分片常数检验函数空间,建立了Crank-Nicolson三次有限体积元格式,获得了数值
B介子衰变,是唯一一类涉及三代夸克弱相互作用的强子衰变,研究B介子的衰变对于精确检验标准模型有重要意义,特别的,末态没有轻子的B介子含粲三体衰变过程可以用于测量CKM矩阵的γ角和β角。本文计算了B(s)到D(s)P1P2的三体衰变过程,其中P1P2可以是π介子或K介子。这个过程在实验上被认为是由中间共振态主导的准两体过程,所以我们可以合理使用准两体的框架对此过程进行研究。准两体过程就是B介子先衰变
单核苷酸多态性(SNPs)是最常见的人类可遗传变异,包括单碱基转换、颠换等,在各类多态性中占80%左右。SNP位点数在人类基因组中可能达到了300万,并且与人类各类疾病密切相关。因此实现SNP的精准检测变得尤为关键,逐渐成为了当下化学生物学、医学等前沿学科研究的热点。为了实现单碱基突变的检测,近年测序型检测和非测序型检测发展迅猛。二代测序成本随着测序技术的高速发展而大幅度降低,但测序技术仍存在一些
时间分布阶微分方程常用于描述一些扩散指数随时间变化的复杂过程,如加速亚扩散过程等,目前已经在诸多领域发挥着重要作用,成为国际学术界的热门研究课题.但是分布阶算子所具有的复杂性和非局部性使得求解分布阶微分方程的精确解困难重重,因此学者们转而求其数值解,并取得了重要进展.在众多算法中,有限元方法以其较强的区域适应性、更灵活的网格剖分、更低的光滑性要求以及较强的通用性等显著优势,备受学者们的青睐.本文研
近年来,二维材料固有磁性的实验发现引起了人们对二维磁性材料的极大兴趣。其中,具有高自旋极化的二维铁磁材料(如,铁磁半导体,铁磁半金属)被认为是下一代自旋电子器件的关键。然而,目前实验发现的铁磁材料存在成本高、居里温度低等缺陷。这极大的限制了二维材料在自旋电子学领域的应用。随着计算硬件水平的提高及理论算法的发展,第一性原理计算成为探寻新型室温铁磁材料的有效途径。基于此,本文采用第一性原理计算方法,设
本文研究了 一类不连续的分数阶Sturm-Liouville边值问题,其边界条件中带有谱参数,且转移条件是分数阶型的.利用算子理论,我们将边界条件和转移条件中的参数相结合,定义了一个新的内积,然后将分数阶边值问题转化为一个新的Hilbert空间中的算子问题,从而得到主问题的特征值和本征函数与该算子的特征值和本征函数的关系.另外文中给出了基本解,且说明了特征函数的解是我们研究的问题的特征值.文章的最