反常扩散方程与随机分数阶偏微分方程的数值方法研究

来源 :兰州大学 | 被引量 : 0次 | 上传用户:zj888666
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
反常扩散方程能够很好的刻画反常动力学的机制,包括空间幂律分布的扩散以及时间长程相关的扩散;因此吸引着各个领域的工作者去建立和研究反常扩散方程.确定性的方程能够呈现事物发展的主要规律,然而宇宙中的随机扰动无处不在,因此想要更全面的刻画事物的发展规律,学者们引入噪声项以刻画随机扰动.于是,随机微分方程的理论研究和数值研究也盛行起来.当方程中含有非局部算子和噪声项时,理论研究和数值研究会变得更具有挑战性;噪声的复杂程度也给研究带来了困难,如tempered高斯噪声.基于这些问题,本文将研究反常扩散方程的数值方法,以及随机分数阶偏微分方程的解正则性和数值逼近.本文由六章构成.第一章,简要地叙述反常扩散方程和随机分数阶偏微分方程研究意义及研究现状;详细说明本文的研究内容及创新之处.第二章,我们介绍了一些预备知识,包括分数阶拉普拉斯算子的定义、分数高斯过程、分数高斯过程的几种模拟方法并通过比较选出最合适的模拟方法.第三章,我们研究了控制tempered分数布朗运动概率密度函数的二维Fokker-Planck方程的数值格式.数值格式的主要挑战来自于时间在0时刻的奇异性.当0<H<1/2时,我们通过变量替换(?)(t2H)=2Ht2H-1(?)t消除了数值计算在0时刻的奇异性;这种形式的变量替换自然的导出了一个不均匀的时间离散格式,而且显著地提高了计算效率.对于H>1/2,为了保证计算的有效性和计算效率,我们引入了时间跨度相关的数值格式和非均匀时间离散化.通过傅立叶方法证明了数值格式的稳定性和收敛性.数值模拟相应的Fokker-Planck方程,我们获得了随机过程的均方位移,它符合tempered分数布朗运动的特征.第四章,非对称tempered分数阶拉普拉斯是各向异性的tempered Levy过程Xt的无穷小生成元,这一章我们研究了Xt的首次退出和非对称tempered分数阶Dirich-let问题.首次退出位置|XτD |和首次退出时间τD所有阶矩的上界被获得.我们发现| XτD|或τD的概率密度函数随|XτD|或τD的增加指数衰减;并且E[τD]~|E[XτD]|,E[τD]~E[|XτD-E[XτD]|2].因为 Δmα/2,λ是各向异性的 tempered Levy 过程Xt的无穷小生成元,因此我们导出了非对称tempered fractional Dirichlet问题的Feynman-Kac representation.此外,通过平均大量的随机过程的轨迹,我们获得了 Dirichlet问题的数值解.第五章,我们讨论了具有tempered分数高斯噪声的分数阶扩散方程.分数阶扩散方程控制subordinated killed布朗运动的概率密度函数.波动的外部源由tempered的分数高斯噪声表示,并且具有局部性.我们首先建立了无穷维tempered分数布朗运动的随机积分的正则性,然后建立了随机分数阶扩散方程的温和解的正则性.本章中,我们采用谱Galerkin方法进行空间逼近,之后将系统转化为一个等价形式,该等价形式比原系统在时间上具有更好正则性的.然后,我们使用半隐式Euler方法离散等价形式的时间导数.根据时空误差分裂技术,我们获得了均方L2-范数意义下的全离散格式的误差估计.大量的数值实验证实了理论估计.第六章,一个subordinated killed布朗运动的无穷生成元(分数阶拉普拉斯)被用来捕捉波传播的幂律衰减性质.这一章,我们研究了以分数阶拉普拉斯为空间算子的随机波方程的数值格式,其噪声项为无穷维布朗运动或分数阶布朗运动.首先,我们建立了随机分数阶波方程温和解的正则性.然后采用谱Galerkin方法进行空间半离散逼近,通过对无穷维高斯噪声的后处理,提高了空间收敛速度.在时间方向上,当温和解的时间导数在均方Lp-范数意义下有界时,我们提出了一种改进的随机trigonometric方法,得到了比现有结果更高的强收敛速度,即时间收敛速度大于1.特别地,我们所提供的方法的时间离散误差收敛速率可以达到2阶,但需要对温和解有一些额外的正则性要求.最后通过数值实验验证了理论误差估计.第七章,总结本文以及展望未来工作方向.
其他文献
学位
学位
电荷-宇称不守恒(CP破坏)是解释宇宙中物质-反物质不对称的基本要素。D介子是自然界含有重夸克的粒子中最轻的一组,其CP破坏在寻找新物理方面具有独特的优势。本文研究了 D介子衰变到中性K介子过程的CP破坏。我们发现在这类过程中存在一种新的CP破坏效应。它是由树图(包括Cabibbo-支持的和Cabibbo-双压低的)振幅和末态介子混合之间的干涉导致的。我们用因子化辅助的拓扑图方法系统地分析了 D→
花粉萌发和花粉管生长是被子植物完成有性生殖不可或缺的关键生理活动。花粉管生长是一种典型的顶端极性生长方式,依赖精准的“运输系统”保障其顶端快速地生长。Rab GTP酶是一类在真核生物中保守存在的囊泡运输相关蛋白,在囊泡形成、转运、栓系、融合、以及不同途径的内膜转运过程中都起着至关重要的作用。其中Rab A、Rab D等亚家族的少数成员被证明参与花粉发育和花粉管极性生长过程。尤其是GFP-Rab A
根的生长需要根尖分生组织细胞的连续分裂。ROOT MERISTEM GROWTH FACTORs(RGFs)是维持根尖干细胞微环境的一类关键多肽类激素。我们实验室和其他实验室前期的研究发现了五个进化关系相近的富含亮氨酸重复序列的类受体激酶(LRR-RLKs),将其命名为RGF1 INSENSITIVES(RGIs)或者RGF1RECEPTORS(RGFRs)。这五个蛋白在感知RGF1的信号并调控拟
强相互作用是人类已知的四种基本相互作用之一,量子色动力学是描述夸克层次上强相互作用的SU(3)局域规范理论。QCD的性质如渐进自由等特性使其呈现出丰富的物理内涵,并产生了复杂的相结构。QCD相图对早期宇宙演化、重离子碰撞以及致密星体结构等方面的研究具有重要意义。人们通过高能重离子碰撞实验产生高温高密的物质,即夸克胶子等离子体,来探索极端状态下物质的形态与特性。人们普遍认为在低化学势区域QCD相图上
本文主要研究了自由罗巴族代数,自由(三)叶型族代数和自由预-李族代数,全文共分为五章.第一章介绍了本文研究课题的背景及其进展,并给出了本文需要的基本概念和一些相关的记号,然后分析了本文的研究动机.第二章首先构造并证明了自由交换罗巴族代数.然后借助括号字,利用Gr(?)bnerShirshov基方法证明了自由非交换罗巴族代数.最后,引入了(三)叶型族代数的概念,证明了权为0(权为?=0)的罗巴族代数
通过介绍隧道工程的定义和对塑料材料的要求,指出运用于隧道工程的塑料材料应当具备良好的力学性能、耐化学腐蚀性能、防水性等特点。根据目前国内外塑料材料在隧道工程的最新研究成果,分别从聚乙烯塑料、聚氨酯塑料以及其他塑料材料方面,综述塑料材料在隧道防水、防火、应急救援、防冻害和堵水等方面的应用现状。指出目前关于塑料材料在隧道工程中应用主要包括材料的加工工艺、性能测试以及实地性能验证等方面,对于材料种类、材
随着无穷维动力系统理论的深入发展,许多由数学物理方程生成的耗散动力系统显现了一定的有限维属性.由此引发了一系列对无穷维动力系统进行有限维约化的研究.经典的惯性流形理论表明,如果一个偏微分方程存在一个N维惯性流形,则其长时间行为可以约化为一个N阶常微分方程组.这本质地简化了对原始偏微分方程动力学行为的理解.目前,惯性流形研究仍是无穷维动力系统中十分重要且具有挑战性的问题之一.本文研究惯性流形及其在耗
有限元方法作为数值求解偏微分方程的有效方法,其思想是把区域离散化,然后用分片多项式函数对解析解进行逼近,因此可以对不规则复杂区域的问题进行高效求解.随着有限元方法研究的深入,已经有很多成熟的有限元程序包在工程界广受欢迎,并在工程计算中大量应用,发挥着不可替代的作用.在科技飞速发展的当代,对实际问题进行模拟的偏微分方程也在快速发展之中,所以探索新有限元方法仍是数值研究领域的重要课题之一.本论文拟对几