时间(缓增)分数阶扩散方程的局部间断Galerkin方法

来源 :西安理工大学 | 被引量 : 0次 | 上传用户:lovejr622
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
分数阶微积分在众多领域有着广泛的应用,而用分数阶扩散方程描述粒子的反常扩散现象则是分数阶微积分的一个典型应用。因分数阶微积分的历史依赖性和非局部性质,增加了分数阶微分方程数值求解时的计算量,并给计算机的存储带来不便,一般通过构造高阶算法和快速算法克服此困难。局部间断Galerkin有限元方法因其极强的灵活性和高精度性能备受关注。本文结合有限差分法和局部间断Galerkin法,研究时间分数阶扩散方程和时间缓增分数阶扩散方程的高阶数值格式,并进行相应的理论分析。主要内容如下:  (1)建立时间缓增分数阶扩散方程的高阶数值格式。首先用局部间断Galerkin有限元方法离散空间导数,得到空间方向具有(k+1)阶精度的半离散格式,其中k为有限元空间中多项式基函数的最高次数。然后用q(q=1,2,3,4)阶缓增加权位移Lubich算子离散时间导数,得到收敛阶为O(τq+hk+l)的全离散格式,并在L2范数意义下证明了半离散格式和全离散格式的稳定性和收敛性。最后数值算例和数值模拟结果表明了该格式的有效性。  (2)建立时间分数阶扩散方程的具有时间二阶精度的两种数值格式。采用局部间断Galerkin有限元方法得到空间半离散格式,并进行理论分析。时间方向采用二阶加权位移Grunwald差分算子和Crank-Nicolson差分格式及二阶向后差分格式逼近,得到收敛阶为O(τ2+hk+1)的全离散格式Ⅰ和格式Ⅱ,并在L2范数意义下进行了详细的理论分析。最后通过算例验证了所构造的两种数值格式是有效的。  (3)建立时间分数阶扩散方程的高阶快速算法。结合离散Caputo时间分数阶导数的快速算法和具有高精度的局部间断Galerkin有限元方法,建立有效求解时间分数阶扩散方程的高阶全离散格式,并在L2范数意义下证明了数值格式是无条件稳定的且收敛阶为O(△t2-α+hk+1+ε)。  时间(缓增)分数阶扩散方程;局部间断Galerkin方法;加权位移Grünwald差分算子;稳定性;收敛性
其他文献
本文主要应用上下解方法、摄动方法和最大值原理, 结合二阶线性椭圆型方程的估计理论和正则性理论, 讨论非线性椭圆型方程正整体解的存在性。 本文分四章进行论述。第二章
本文考虑的是谱负的Lévy过程,也即没有正跳的Lévy过程。把开始于u(u≥0)的谱负Lévy过程看作是推广的风险模型,文章得出了破产时刻和破产瞬间前后余额的三者联合密度函数。然
本文研究组合弹性结构的数学理论、数值分析方法及其在建筑结构分析等领域的应用,该研究将涉及组合弹性结构的本构关系、变分模型、微分方程模型和数值分析方法。本文选择弹性
规范化管理是所有企业管理者都在努力实现的管理目标,在目标实现过程中不乏有一些理论体系与方法论做支撑,但如何通过信息技术手段与规范化管理的理论体系与方法论融合,实现
我们知道,当随机微分方程的系数函数足够光滑时其Euler逼近的弱收敛阶数为1.   但当系数函数不光滑且仅仅满足Lipschitz条件时,Euler逼近的弱收敛阶数似乎不再为1.直到2003
新媒体技术的发展和应用,使我们进入一个“互联网+”的全新时代,大学生社会主义核心价值观教育工作也随之面临到全新的挑战、发展和机遇.在这一背景下,高校大学生核心价值观
谈话,是教师与学生沟通感情,对学生施行教育的重要手段.艺术性的谈话,能使学生如沐春风,开启学生的心灵之窗,教会他们做人的道理,从而建立和谐的师生关系,让学生在善意的教诲
本篇文章主要考虑的是如何最大化全部折现分红和最终的固定资产P的和。其中最终固定资产P表示保险公司在破产时所拥有的固定的可折算成资金的资产。本篇文章中,保险公司的资金
序集抽样方法是上世纪50年代McIntyre在寻求能较好估计牧场草的产量时提出来的.它是假定在一个无穷总体中,对每一个体进行准确测量的费用很高或时间很长,而对他们或其相关的量
董事长作为股东利益代表行使决策权,总经理作为公司经营者行使执行权,两者本应在《公司法》和公司章程等法律或制度构筑好的权力框架下各司其职、各负其责,相互合作、相互支