时空分数阶扩散方程的扩展混合有限元方法

发布时间：2020-07-29 12:04

【摘要】：本文第一部分讨论下列由双边Riemann-Liouville分数阶导数刻画的守恒型分数阶扩散方程其中,Ω =(0,1),0β1,0 ≤θ ≤ 1,K是常扩散系数,f ∈ L2(Ω)表示源或汇项;(?)和D=dd/dx分别为关于时间、空间的一阶导数算子,0Ixβ和xI1β是由(2.2.1)式定义的β阶左、右Riemann-Liouville分数阶积分算子.大量的实验表明,上述分数阶扩散方程较二阶扩散方程能更准确地刻画诸如湍流,混沌动力学,粘弹性力学,地下污染渗流传输等反常扩散或非菲克现象,具有鲜明的物理背景与广泛的应用.由于难于利用傅里叶变换,拉普拉斯变换等方法求得一般分数阶扩散问题的解析解,人们通常采用有限元,有限差分方法求其数值解.为了满足工程上的需要,理想的数值方法不仅应同时关注未知函数及其通量,而且需满足某种意义下的守恒律以反映物质扩散输运过程的数学物理特征.常规的思路是,通过引入分数阶通量p=-K(0Ixβ+(1-θ)xI1β)作为中间变量,建立恰当的鞍点变分原理,以此构造标准的混合有限元方法.但由于算子oIxβDu,xI1βDu和刀不是对称的,我们无法构造出恰当的双线性形式b(·,·)来满足inf-sup条件.对此,陈焕贞和王宏[10]对空间分数阶扩散方程,通过引入q = Du作为第二个中间变量,建立了相应的鞍点原理,据此构造了能同时高精度逼近未知函数与分数阶通量,且保持单元守恒的扩展混合有限元方法.在本文中,我们借鉴[10]的思路,通过引入分数阶通量和未知函数的导数作为中间变量,对上述依赖时间的分数阶扩散方程建立了相应的鞍点变分形式,构造了局部守恒,可同时高精度逼近未知函数,未知函数导数和分数阶通量的扩展混合有限元格式,并利用向后欧拉方法对时间导数进行离散,建立了全离散扩展混合有限元格式.进一步,我们利用所导出的混合元方程组系数矩阵的可逆性,证明了全离散格式解的存在唯一性;借助离散的Growall不等式,证明了该离散格式的无条件稳定性.在数值分析中,利用混合元投影算子的投影误差估计,得到全离散扩展混合有限元格式的收敛阶为O(τ十hmin{k+1,s-1 +β/2),这里τ,h,k分别表示时间步长、空间步长与有限元空间指数.数值实验表明,文中提出的扩展混合元数值格式可有效地模拟由上述守恒型分数阶扩散方程所刻画的反常扩散输运模型.但大量的反常扩散或非菲克输运过程都表现为超扩散现象,即关于时间与空间的导数均为分数阶情形.因此,在本文的第二部分,我们主要研究下列时空分数阶扩散问题其中,(?)αu(x,t)表示α,0α1阶Caputo 导数.我们仍通过引入未知函数的通量p =-K(θ0Ixβ+(1-θ)xI1β)Du和导数q =Du作为中间变量,以此建立相应的鞍点变分形式与相应的扩展混合有限元格式,并对α阶的Caputo导数采用L1格式离散,构造出了全离散扩展混合有限元方法.在数值分析中,我们采用数学归纳的思想以代替离散的Gronwall不等式,通过复杂的分析论证,建立了全离散扩展混合元方法的收敛性理论,关于未知函数u收敛阶为O(τ2α十hmin{k+1,s-1+β/2}),关于函数导数与分数阶数值通量p的收敛阶为O(τ23α/2+hmrn-k+1,s-1+β/2}).文中数值实验表明,所提出的L1全离散扩展混合有限元格式具有理想的数值逼近效果.
【学位授予单位】：山东师范大学
【学位级别】：硕士
【学位授予年份】：2018
【分类号】：O241.82

【相似文献】

相关期刊论文 前10条

1 卢曰海;丘东元;张波;陈艳峰;;大功率分数阶电感的电路实现[J];电源学报;年期

2 王娴;李东;;分数阶混沌系统的间歇控制同步[J];重庆工商大学学报(自然科学版);2018年04期

3 李明;陈旭;郑永爱;;一类分数阶混沌系统的自适应滑模同步[J];扬州大学学报(自然科学版);2016年03期

4 程春蕊;朱军辉;毛北行;;一类分数阶混沌系统的投影同步[J];河南科学;2016年11期

5 郝孟丽;任勤;;标量控制下的分数阶Lü系统的参数辨识和自适应同步[J];河南理工大学学报(自然科学版);2017年01期

6 王志兰;;分数阶电路阶跃响应特性研究[J];电子测试;2016年24期

7 李亮;李庆宾;毛北行;;分数阶同步发电机系统的混沌同步[J];河南科学;2017年03期

8 李特;袁建宝;吴莹;;一类不确定分数阶混沌系统同步的自适应滑模控制方法[J];动力学与控制学报;2017年02期

9 赖满丰;金玲玉;;分数阶Klein-Gordon-Schr?dinger方程弱解的存在性[J];佛山科学技术学院学报(自然科学版);2017年03期

10 田魏巍;;非线性分数阶动力系统的控制研究[J];教育现代化;2017年22期

相关会议论文 前10条

1 罗绍凯;;分数阶动力学基本理论与方法的研究进展[A];第十二届全国分析力学学术会议摘要集[C];2016年

2 孟瑞繁;殷德顺;;分数阶本构模型描述流变现象时间效应的物理意义[A];中国力学大会-2015论文摘要集[C];2015年

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

4 顾葆华;单梁;李军;王执铨;;一种新分数阶混沌系统及其复合快速同步控制[A];2009年中国智能自动化会议论文集（第七分册）[南京理工大学学报（增刊）][C];2009年

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

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

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

8 王淑英;赵建峰;常迎香;李险峰;;基于假分数阶统一混沌系统同步的图像加密[A];第十五届全国非线性振动暨第十二届全国非线性动力学和运动稳定性学术会议摘要集[C];2015年

9 刘杰;董鹏真;尚钢;;分数阶非线性系统动力学分析中数值算法可靠性及其诱导的复杂现象[A];中国力学学会学术大会'2009论文摘要集[C];2009年

10 杨海天;赵潇;;蚁群算法求解二维分数阶黏弹性参数反问题[A];中国力学学会学术大会'2009论文摘要集[C];2009年

相关博士学位论文 前10条

1 赵猛;时空分数阶偏微分方程的快速算法及其应用[D];山东大学;2018年

2 李政;若干非局部和非线性问题的数值方法[D];华东师范大学;2018年

3 Azmat Ullah Khan Niazi;分数阶中立型微分方程的稳定性和可控性[D];安徽大学;2018年

4 郝朝鹏;双边分数阶扩散方程的高精度和高效数值方法[D];东南大学;2017年

5 练婷婷;Banach空间中分数阶发展系统的能控性与优化控制问题[D];扬州大学;2018年

6 姜翠美;分数阶复动力系统的定性分析与同步研究[D];山东大学;2017年

7 孙红;基于多项式插值逼近的分数阶偏微分方程高精度差分方法[D];东南大学;2016年

8 吴艳萍;分数阶非线性系统同步与非线性电路理论若干问题研究[D];西北农林科技大学;2015年

9 王乔;分数阶混沌系统控制与同步理论研究[D];浙江大学;2015年

10 丁策;机载光电稳定平台的分数阶控制研究[D];中国科学院研究生院(长春光学精密机械与物理研究所);2013年

相关硕士学位论文 前10条

1 袁琼;时空分数阶扩散方程的扩展混合有限元方法[D];山东师范大学;2018年

2 刘娅芳;时间变分数阶扩散方程的紧致差分格式[D];山东大学;2018年

3 洪丹;分数阶微分方程约束优化问题的预处理方法[D];华东师范大学;2018年

4 闫书庆;带有一般外力的时间分数阶福克—普朗克方程的有限差分格式[D];山东大学;2018年

5 乐琪;基于分数阶方程的MRI模型研究[D];华东师范大学;2018年

6 魏龙飞;分数阶粘弹性模型参数识别及其有限元实现[D];西南交通大学;2018年

7 董佳华;隐式分数阶微分方程及耦合系统的研究[D];武汉科技大学;2018年

8 张薇;分数阶积分微分方程数值算法研究[D];电子科技大学;2018年

9 刘娜;分数阶长短波方程组的吸引子和近似惯性流形[D];鲁东大学;2018年

10 蒋伟;几类带边值条件的分数阶微分方程解的存在性研究[D];安徽大学;2018年

本文编号：2773941