变系数分数阶扩散方程的混合型Galerkin变分格式和快速算法

来源 :山东师范大学 | 被引量 : 0次 | 上传用户:lvy1234
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文我们考虑的是如下2-β(0<β<1)阶分数阶扩散方程,(此处公式省略)  这里未知函数c表示污染物的浓度,K(x)是扩散项系数,有正的上下界(此处公式省略)  本文通过引入通量函数u=-K(x)D c作为中间变量,对变系数双边分数阶扩散方程推导了一种混合型Galerkin变分形式,在空间上证明变分解的存在唯一性.在对解有一定正则性假设下,建立了(0.0.1)和变分形式的等价性.然后基于混合型Galerkin变分格式的好的性质的基础上,在通常所用的有限元空间上推导了混合型有限元离散格式,并且证明有限元解的存在唯一性以及给出未知函数和中间变量的误差估计.  由于分数阶微积分算子的非局部性,由分数阶扩散方程的数值方法产生的系数矩阵通常为满阵或稠密矩阵.传统的求解方法需要O(N3的计算量和O(N2)的存储量,其中N为网格节点.在对混合型有限元格式求解过程中,我们发现混合有限元格式所形成的系数矩阵具有特殊结构.系数矩阵由四个分块小矩阵构成,其中一个零阵、两个三对角的稀疏矩阵以及一个具有Toeplitz- l i k e结构的矩阵.通过引入快速傅里叶变换,结合共轭梯度法我们构造出了每次迭代计算量为 O(N logN)以及相应的存储量为O(N)的快速共轭梯度法.同样的,为了减少迭代次数,我们又对快速共轭梯度法做了进一步的改进,引入了预处理技术,构造了预条件快速共轭梯度算法,减少了迭代次数,同时保证了收敛精度.并且我们可以证明预条件的快速算法的总计算量仍然为O(N logN).本文提供的数值算例也验证了该快速算法的良好性质.
其他文献
摘 要:通过收集近些年来石化行业发生的事故以及自己日常从事安全管理工作的一些感悟,结合现代安全管理理论与心理行为分析理论,对人员习惯性违章的心理进行一些简单分析,为下一步保障油库安全、制订防范习惯性违章行为措施提供依据。  关键词:习惯性 违章 心理 透析  一直以来关注安全生产的人们不断在思考一个问题,就是如何杜绝因人员习惯性违章而造成的事故。笔者通过收集近些年来石化行业发生的事故以及自己日常从
移动IP是一种在全球Internet上提供移动功能的方案,它提供了一种IP路由机制使移动节点可以用一个永久IP地址连接到任何链路上.移动IP技术的目的就是无论接连在本地链路还是移动
本文研究了一类混合型数论函数的均值估计问题,这一方向一直是数论中备受关注的方向之一.本文在研究的问题中联系了尖形式的傅立叶系数A/(n),除数和函数a(n)和欧拉函数利n),有
Improving vehicle fuel consumption,performance and aerodynamic efficiency by drag reduction especially in heavy vehicles is one of the indispensable issues of a
非线性抛物方程的研究是偏微分方程的重要组成部分,研究其解的存在及爆破是非线性发展方程理论研究的一个非常重要的方向。本文主要研究两类非线性发展方程解的爆破行为,重点是
学位
本文主要研究了我国高校学费现状及定价的相关模型,分析了影响高校学费的几个重要因素,并结合BP(Back Propagation,误差反向传播神经网络)神经网络和平均影响值(MIV,Mean Impact V
在本论文中,我们研究了一类无限维Z2-分次但非有限Z-分次的Block型李代数B(q)的结构理论和表示理论,其中q是一个非零复参数.   第一章介绍了研究背景以及本论文的主要结果.Blo
近年来,随着离散系统在经济、物理和工程技术方面对广泛运用,使可以描述离散型变量的差分系统得到学者们的重视,并成为计算数学、应用数学和系统科学研究的热点问题.尽管已有许多
本文主要研究了关于M aass尖形式傅里叶系数在整变量三元二次型中的均值分布。论文主要运用经典圆法与指数和来估计该均值问题的上界,丰富了关于傅里叶系数性质的结果,并对其进