守恒型分数阶扩散方程的混合有限元数值方法

发布时间：2017-09-09 14:21

  本文关键词：守恒型分数阶扩散方程的混合有限元数值方法

  更多相关文章： 守恒型分数阶扩散方程 双边分数阶导数 鞍点理论 混合变分原理 扩展混合有限元方法 最优阶误差估计 共轭梯度法 快速傅里叶变换 预处理快速共轭梯度法 数值实验

【摘要】：本文讨论由双边Riemann-Liouville导数刻画的守恒型分数阶扩散方程,这里K是扩散系数,f ∈L2(Ω)是源项或汇项；D=(?)是一阶导数算子,0(?)和(?)是由(2.2.1)和(2.2.3)式定义的β阶左、右Riemann-Liouville分数阶积分算子.(?)分别对应左、右Riemann-Liouville导数刻画的单边分数阶扩散方程.上述分数阶扩散方程刻画了依赖于全局性态的反常或非菲克扩散现象.为了满足工程应用中对扩散通量的辨识需求,无论是传统的二阶扩散方程还是分数阶扩散方程,所设计的数值方法都应同时关注未知函数及其通量,且满足质量守恒律以反映扩散问题的原始数学物理特征.与二阶扩散方程比较,分数阶微积分算子的非局部性质导致其数值离散格式的系数矩阵为非稀疏矩阵,计算复杂性显著增大.因此,构造快速算法也已成为高性能数值模拟分数阶扩散问题的极富挑战性的内容之一对此,本文中我们基于鞍点理论与负指数分数阶导数空间,引入分数阶通量p=-K(θ0Ixβ+(1-θ)xI1β)q与导数q=DU作为中间变量,在H1(Ω)×H-β/2)(Ω)× L2(Ω)上建立了与分数阶扩散问题等价的混合变分原理.据此,构造了能同时高精度数值逼近未知函数、扩散通量以及未知函数导数的扩展混合有限元方法.理论分析表明,文中所构造的扩展混合有限元离散格式逐单元保持质量守恒,这一点对实际的工程计算是至关重要的；离散格式的解存在唯一,且具有对未知函数、分数阶扩散通量及其函数导数具有在某种意义的最优L2或H-β/2-模收敛精度.在收敛性分析中,我们摒弃了文献[16]中基于未经证明正则性假定的对偶论证方法,代之以利用真解在相应空间上投影的良好逼近性质,建立了不依赖于对偶问题正则性假定的最优阶收敛性理论.我们还给出一系列数值实验结果,表明文中所提离散格式具有与收敛性理论分析相匹配的数值逼近精度,从而说明了离散格式的有效性.注意到分数阶微分算子的非局部性质所导致的离散问题系数矩阵复杂性,若使用传统的高斯迭代法求解,将会产生O(N3)和O(N2)的计算量与存储量,致使计算时间过长甚至无法运行.因此,设计快速求解上述离散格式的有效算法是十分必要的.我们发现,离散格式的刚度矩阵可分块为四个零阵、四个稀疏矩阵和一个Toeplitz矩阵,而快速傅里叶变换在求解具Toeplitz型矩阵的矩阵-向量积时,具有理想的计算量O(Nlog N)鉴于此,我们将快速傅里叶变换与共轭梯度法结合,构造出了每次迭代计算量为O(N log N),存储量为O(N)的扩展混合有限元快速算法(FCG),较高斯消去法有了明显的改善.遗憾的是,在设计的FCG中,由于混合有限元离散格式系数矩阵的病态性质,造成了迭代次数强烈的依赖于未知量个数N,使FCG的计算量未达到理想的O(N log N).对此,我们将预处理技术与FCG相结合,设计了扩展混合有限元预处理快速算法(PFCG).数值实验表明,PFCG的迭代次数不依赖于未知量个数,计算量与存储量分别达到了理想的0(Nlog N)和O(N).
【关键词】：守恒型分数阶扩散方程 双边分数阶导数 鞍点理论 混合变分原理 扩展混合有限元方法 最优阶误差估计 共轭梯度法 快速傅里叶变换 预处理快速共轭梯度法 数值实验
【学位授予单位】：山东师范大学
【学位级别】：硕士
【学位授予年份】：2016
【分类号】：O241.82
【目录】：

• 中文摘要5-7
• 英文摘要7-10
• 第一章 绪论10-13
• 第二章 守恒分数阶扩散方程的扩展混合有限元方法13-33
• 2.1 引言13-14
• 2.2 预备知识14-16
• 2.3 鞍点变分格式16-21
• 2.4 扩展混合有限元方法21-24
• 2.5 收敛性分析24-30
• 2.6 数值算例30-33
• 第三章 守恒分数阶扩散方程的扩展混合有限元快速算法33-47
• 3.1 引言33
• 3.2 共轭梯度法33-35
• 3.3 快速共轭梯度法35-39
• 3.4 预处理快速共轭梯度法39-44
• 3.5 数值算例44-47
• 第四章 评注47-48
• 参考文献48-52
• 攻读学位期间发表的学术论文52-53
• 致谢53

【相似文献】

中国期刊全文数据库 前10条

1 王德金;郑永爱;;分数阶混沌系统的延迟同步[J];动力学与控制学报;2010年04期

2 杨晨航,刘发旺;分数阶Relaxation-Oscillation方程的一种分数阶预估-校正方法[J];厦门大学学报(自然科学版);2005年06期

3 王发强;刘崇新;;分数阶临界混沌系统及电路实验的研究[J];物理学报;2006年08期

4 夏源;吴吉春;;分数阶对流——弥散方程的数值求解[J];南京大学学报(自然科学版);2007年04期

5 张隆阁;;一类参数不确定混沌系统的分数阶自适应同步[J];中国科技信息;2009年15期

6 陈世平;刘发旺;;一维分数阶渗透方程的数值模拟[J];高等学校计算数学学报;2010年04期

7 辛宝贵;陈通;刘艳芹;;一类分数阶混沌金融系统的复杂性演化研究[J];物理学报;2011年04期

8 黄睿晖;;分数阶微方程的迭代方法研究[J];长春理工大学学报;2011年06期

9 蒋晓芸,徐明瑜;分形介质分数阶反常守恒扩散模型及其解析解[J];山东大学学报(理学版);2003年05期

10 陈玉霞;高金峰;;一个新的分数阶混沌系统[J];郑州大学学报(理学版);2009年04期

中国重要会议论文全文数据库 前10条

1 李西成;;经皮吸收的分数阶药物动力学模型[A];中国力学学会学术大会'2009论文摘要集[C];2009年

2 谢勇;;分数阶模型神经元的动力学行为及其同步[A];第四届全国动力学与控制青年学者研讨会论文摘要集[C];2010年

3 张硕;于永光;王亚;;带有时滞和随机扰动的不确定分数阶混沌系统准同步[A];中国力学大会——2013论文摘要集[C];2013年

4 李常品;;分数阶动力学的若干关键问题及研究进展[A];中国力学大会——2013论文摘要集[C];2013年

5 李常品;;分数阶动力学简介[A];第三届海峡两岸动力学、振动与控制学术会议论文摘要集[C];2013年

6 蒋晓芸;徐明瑜;;时间依靠分数阶Schr銉dinger方程中的可动边界问题[A];中国力学学会学术大会'2009论文摘要集[C];2009年

7 王花;;分数阶混沌系统的同步在图像加密中的应用[A];第二届全国随机动力学学术会议摘要集与会议议程[C];2013年

8 王在华;;分数阶动力系统的若干问题[A];第三届全国动力学与控制青年学者研讨会论文摘要集[C];2009年

9 张硕;于永光;王莎;;带有时滞和随机扰动的分数阶混沌系统同步[A];第十四届全国非线性振动暨第十一届全国非线性动力学和运动稳定性学术会议摘要集与会议议程[C];2013年

10 李西成;;一个具有糊状区的分数阶可动边界问题的相似解研究[A];中国力学大会——2013论文摘要集[C];2013年

中国博士学位论文全文数据库 前10条

1 陈善镇;两类空间分数阶偏微分方程模型有限差分逼近的若干研究[D];山东大学;2015年

2 任永强;油藏与二氧化碳埋存问题的数值模拟与不确定性量化分析以及分数阶微分方程的数值方法[D];山东大学;2015年

3 蒋敏;分数阶微分方程理论分析与应用问题的研究[D];电子科技大学;2015年

4 卜红霞;基于分数阶傅里叶域稀疏表征的CS-SAR成像理论与算法研究[D];北京理工大学;2015年

5 杨变霞;分数阶Laplace算子的谱理论及其在微分方程中的应用[D];兰州大学;2015年

6 邵晶;几类微分系统的定性理论及其应用[D];曲阜师范大学;2015年

7 方益;分数阶Yamabe问题的一些紧性结果[D];中国科学技术大学;2015年

8 王国涛;几类分数阶非线性微分方程解的存在理论及应用[D];西安电子科技大学;2014年

9 陈明华;分数阶微分方程的高阶算法及理论分析[D];兰州大学;2015年

10 孟伟;基于分数阶拓展算子的灰色预测模型[D];南京航空航天大学;2015年

中国硕士学位论文全文数据库 前10条

1 黄志颖;非线性时间分数阶微分方程的数值解法[D];华南理工大学;2015年

2 赵九龙;基于分数阶微积分的三维图像去噪增强算法研究[D];宁夏大学;2015年

3 楚彩虹;单载波分数阶傅里叶域均衡系统及关键技术研究[D];郑州大学;2015年

4 全晓静;非线性分数阶积分方程的Adomian解法[D];宁夏大学;2015年

5 黄洁;非线性分数阶Volterra积分微分方程的小波数值解法[D];宁夏大学;2015年

6 庄峤;复合介质中时间分数阶热传导正逆问题及其应用研究[D];山东大学;2015年

7 高素娟;分数阶延迟偏微分方程的紧致有限差分方法[D];山东大学;2015年

8 赵珊珊;时—空分数阶扩散方程的快速算法以及MT-TSCR-FDE的快速数值解法[D];山东大学;2015年

9 王珍;分数阶奇异边值问题的研究[D];山东师范大学;2015年

10 冯静;一类分数阶奇异脉冲边值问题正解的存在性研究[D];山东师范大学;2015年

，

本文编号：820943