时间分布阶偏微分方程的几类有限元方法研究

来源 :内蒙古大学 | 被引量 : 0次 | 上传用户:whenhm
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
时间分布阶微分方程常用于描述一些扩散指数随时间变化的复杂过程,如加速亚扩散过程等,目前已经在诸多领域发挥着重要作用,成为国际学术界的热门研究课题.但是分布阶算子所具有的复杂性和非局部性使得求解分布阶微分方程的精确解困难重重,因此学者们转而求其数值解,并取得了重要进展.在众多算法中,有限元方法以其较强的区域适应性、更灵活的网格剖分、更低的光滑性要求以及较强的通用性等显著优势,备受学者们的青睐.本文研究了时间分布阶偏微分方程的几类有限元方法,重点讨论了H1-Galerkin混合有限元(GMFE)方法、两层网格有限元方法和交替方向隐式(ADI)有限元方法.具体的研究内容概括为以下几个部分:第一部分,使用基于二阶σ格式的H1-Galerkin混合有限元方法数值求解带有非线性项的时间分布阶反应扩散模型.时间方向上使用二阶σ格式逼近,空间方向上采用混合元算法进行离散,进而形成全离散格式.文中详细给出了格式的稳定性证明,并推导了未知函数p及中间辅助函数u在L2范数下的最优误差估计.最后用数值算例给出格式在时间和空间上的误差和收敛阶,可以看到时间和空间方向均达到最优阶,与理论分析结果一致,验证了理论分析结果的正确性.同时通过数值图像可以看到数值解与精确解的图像吻合,这说明我们给出的H1-GMFE方法对于求解非线性时间分布阶偏微分方程是有效的.第二部分,通过建立基于加权移位Grunwald差分(WSGD)算子的两层网格ADI有限元格式数值求解二维非线性时间分布阶反应扩散方程.使用WSGD算子结合数值求积公式对时间分布阶导数进行逼近,进一步形成非线性时间分布阶反应扩散模型的两层网格ADI有限元全离散格式.文中给出了数值格式的稳定性分析和误差估计的推导过程,得到了空间方向的二阶收敛精度.该方法降低了分布阶微分方程的维度、减少了存储空间并提高了计算效率.根据数值算例可以看出数值格式在空间方向上为二阶收敛,符合理论分析结果.同时,根据数值解与精确解的对比图,可以看到图像是相吻合的,这表明该算法可以有效的数值求解非线性时间分布阶反应扩散方程.第三部分,首次讨论了使用基于移位分数阶梯形公式(SFTR)的ADI有限元方法求解非线性时间分布阶反应扩散耦合系统.我们用SFTR结合数值求积公式逼近时间分布阶导数,进一步形成非线性分布阶反应扩散耦合系统的ADI有限元全离散格式.文中详细证明了格式的稳定性并借助投影算子得到了未知函数u和v的误差估计结果.理论分析和数值计算结果均显示,基于SFTR的ADI有限元算法可以在空间方向上达到最优收敛阶.该方法可以有效地降低对存储量的巨大需求,提高计算效率.
其他文献
在大尺度运动中,大气与海洋运动表现为非线性Rossby孤立波,孤立波因其能保持振幅结构不发生变化,在解释实际问题时有重要意义。本文研究了大尺度大气与海洋运动的非线性正压–斜压相互作用的演化机制,特别是对非线性Rossby孤立波振幅相干结构的激发、传播到消散进行深入讨论。从经典的两层准地转模型出发,引入正压流函数与斜压流函数,并考虑底部地形和–平面近似对其的影响,分别采用多重尺度法和摄动展开法,得到
本文首次采用QCD求和规则计算了Λb→Λc/(?)b→(?)c跃迁的半轻形状因子,考虑了包括微扰(dim-0)、夸克凝聚(dim-3)、夸克-胶子混合凝聚(dim-5)的贡献。由于通常认为胶子-胶子凝聚贡献较小,因此本文不予考虑。本文首先给出初末态的内插流,并定义关联函数。然后分别在强子层次和夸克层次上计算关联函数。最后令强子层次和夸克层次的关联函数表达式相等,并利用夸克-强子对偶减除激发态和连续
小鼠多能干细胞(mouse pluripotent stem cell,PSC)的多能性状态根据其多能性调控特征的可以分为原始态(na?ve)和启发态(primed)。随着研究人员对胚胎发育过程诱导机制的理解不断加深,人们陆续通过诱导获得了一些具有独特多能性特征的新型多能干细胞。本实验室分别利用在化学成分明确的培养基中使用Activin A+b FGF(AF)和Activin A+BMP4+CHI
体细胞重编程技术在干细胞生物学,疾病建模,治疗性克隆等方面具有巨大的潜力,但是影响体细胞重编程效率的障碍以及不同体细胞重编程技术之间的差异仍未被完全阐明。本研究以牛和小鼠的体细胞核移植(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的精准检测变得尤为关键,逐渐成为了当下化学生物学、医学等前沿学科研究的热点。为了实现单碱基突变的检测,近年测序型检测和非测序型检测发展迅猛。二代测序成本随着测序技术的高速发展而大幅度降低,但测序技术仍存在一些