收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

分数阶扩散方程的几种数值解法

苏丽娟  
【摘要】: 尽管分数阶微积分的历史几乎和整数阶的一样长,但是由于缺少相关的实际应用背景,分数阶微积分在其初期发展十分缓慢.众所周知,对于解释和模拟许多应用科学领域的动力学过程,经典微积分都是一个强有力的工具.但是,越来越多的实验和现实告诉我们,在自然界的反常动力学中有许多复杂系统,不能用经典的导数模型来描述.因此,在最近的十几年里,分数阶微积分已经被应用于几乎所有科学、工程和数学的领域中去. 物理中,反常扩散或许是一种最常研究的复杂问题.我们利用分数阶导数,可以将经典的整数阶扩散与波的偏微分方程,推广到时间和空间的分数阶上去.进而再扩展到各类非线性方程并给出其初边值问题的解,是近几年来分数阶微积分应用的一个主要领域.一般来讲,这些问题大都具非常重要的实际应用背景,如在分形和多孔介质中的弥散、半导体物理、湍流及凝聚态物理等. 本文主要研究一些分数阶扩散方程及其数值解法,共由四个彼此相关而又相互独立的章节构成.第一章简要介绍了分数阶微积分的历史、理论及其应用,以及文中将用到的一些基本知识,和相关数值解的现有研究成果;第二章和第三章研究的都是双边空间分数阶对流扩散方程,在这两章中我们分别给出了此类方程的几类不同有限差分法,主要有分数阶权平均法、改进型权平均法和特征有限差分法等;而在最后一章中,我们则是给出了时间分数阶扩散方程的一种高精度隐式数值解法. 第一章为序言.首先了介绍分数阶微积分的历史及其发展情况,并给出了几种常用的分数阶算子定义以及它们的一些基本性质,例如:Riemann-Liouville分数阶算子,Caputo分数阶算子和Grunwald-Letnikov分数阶算子等,同时还列出了几个相关的运算性质.然后,在§1.3中,我们对Mittag-Leffler型特殊函数和它的基本性质也进行了一定的叙述,这类特殊函数主要包括单参数的Mittag-Leffler函数和两参数的广义Mittag-Leffler函数.这类特殊函数常常是很多分数阶微分方程的基本解,其它具有类似性质的特殊函数还有Wright函数和H-fox函数,等等. 此外,在本章中的§1.4,我们还归纳叙述了目前为止,几类常见的分数阶微分方程的一些数值解法.例如,有限差分法,有限元法,微分变换法,Adomian区域分解法,变分迭代法,同伦摄动法,等等.同时对每种方法分别列举出了一些相关研究成果.最后,在本章的最后一节中,我们较详细的介绍了分数阶微积分在当前非线性物理复杂系统的各个领域中的应用. 在接下来的章节中,我们将研究两种不同的反常扩散模型.在第二章中,主要研究1维空间分数阶对流扩散方程.我们根据移位Grunwald公式离散Riemann-Liouville分数阶导数,从而提出了方程的分数阶权平均法.通过理论研究和算例分析,可以得知以前出现过的一些相关数值算法,它们大都是此方法的某些特例.在§2.3中,我们利用圆盘定理和矩阵法证明了分数阶权平均法的稳定性,具体理论结果由定理2.1详细给出. 然后,在§2.4中,我们又讨论了分数阶权平均法的一种新的改进格式,并再次给出了相关稳定性分析.最后,则是用数值例子来验证理论的正确性,同时又计算了分数阶权平均法的特例,分数阶Crank-Nicholson (FCN)法.显然,无条件稳定又拥有2阶时间精确度的FCN法更好一些.本章部分内容已经公开发表在Physics Letters A. 在第三章中,关于双边空间分数阶对流扩散方程,据我们所知,目前为止,它的数值解法全都是Eulerian法.结果,这些方法都具有和2阶对流扩散方程相同的数值局限性.在本章中,结合移位Griinwald-Letnikov有限差分过程以及Lagrangian法,我们在§3.3中首次提出了一种分数阶特征有限差分法(CFDM).此法保留了2阶对流扩散方程特征法和分数阶对流扩散方程有限差分法的所有数值优点.在§3.4中,我们证明了这种方法是无条件稳定、相容和收敛的,并且给出了本方法误差估计的最大值. 在§3.5中,我们给出了一个实际算例的数值模拟,并把分数阶特征有限差分法和其它的分数阶标准差分法相比较.算例结果表明,这种分数阶新CFDM在精度和稳定性上都大大优于其它已知方法,例如显式迎风差分法和隐式迎风差分法等.并且,此法对于对流占优问题,显得尤为高效、优越.本章内容已投到Journal of Computational Physics. 在第四章中,我们主要讨论的是一类时间分数阶扩散方程.反常次扩散运动是复杂系统中一个特别重要的内容,如在一些有机和无序材料中,它的运动路径被一些几何或能量因子约束着.对于反常次扩散随机游走过程的数学模型,一般扩散方程则会被Riemann-Liouville分数阶时间扩散方程所替代.分析表明,这些分数阶模型显然比经典的整数阶模型更加符合实际背景. 在本章中,首先我们我们利用移位Grunwald公式来逼近时间分数阶导数,并且使用中心差分格式去逼近1阶时间导数和2阶空间导数,从而提出了此类扩散方程的一种新的隐式差分法.它是一种三层差分格式,其中第一时间层的数值解可以由全隐式格式或Crank-Nicholson格式给出,这两种格式都是无条件稳定的.接着,我们利用广义化Fourier-Von Neumann分析法,证明了这种新方法的无条件稳定性,并导出此法关于时间的2阶精确度.最后,则是用数值试验和对比法,来验证和观察本算法的性质和特征.本章内容已投到Applied Mathematics and Computation.


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 吕龙进;肖建斌;;分数阶反常扩散方程及其精确解[J];杭州电子科技大学学报;2008年04期
2 孙建强;秦孟兆;戴桂冬;;解扩散方程的指数时间差分方法[J];数值计算与计算机应用;2008年04期
3 陈泮勤;非定常两维平流扩散方程的数值实验[J];气象学报;1985年01期
4 张忠诚;;Zn-MnO_2干电池扩散方程的求解[J];山东大学学报(工学版);1989年02期
5 陈军;李秋姝;;一类非线性扩散方程解的估计及非负解的存在性[J];贵州大学学报(自然科学版);1990年01期
6 张庆红;郁国瑞;;扩散方程的具有振动边值的Dirichlet问题[J];河北师范大学学报(自然科学版);2006年03期
7 陈昊;钟尔杰;;扩散方程的高精度恒稳定差分格式[J];电子科技大学学报;2007年S1期
8 杨传富;;扩散方程系数的半逆问题[J];数学年刊A辑(中文版);2011年01期
9 杨志峰;周同明;;含源扩散方程的一种高精度差分方法[J];北京师范大学学报(自然科学版);1992年01期
10 范学刚;一类扩散方程逆问题的稳定性[J];焦作师范高等专科学校学报;2003年01期
11 邓滨,施保昌,王广超;求解含源项扩散方程的两种格子BGK模型比较[J];华中科技大学学报(自然科学版);2004年12期
12 张新东;胡月宏;;同伦分析方法求具有边界条件的扩散方程的精确解[J];山东大学学报(理学版);2008年12期
13 马昌前;扩散方程的岩浆动力学意义[J];地质科技情报;1988年02期
14 苏秀娟;用K理论研究南京大厂地区污染物扩散[J];南京气象学院学报;1989年03期
15 张海亮,于鸣歧;一类非线性反应扩散方程解的Blow-up问题[J];数学杂志;1997年04期
16 梅家斌;分形在电镀中的应用[J];中南民族大学学报(自然科学版);2002年01期
17 赵明登,徐林春,严佳庆,槐文信;求解二维扩散方程的有限近似解法[J];水动力学研究与进展A辑;2005年05期
18 张锴;;点源迭加求解扩散方程[J];湖南文理学院学报(自然科学版);2007年04期
19 马菊香;;扩散方程的一种跳点格式[J];宁夏师范学院学报;2010年06期
20 马菊香;;扩散方程的一种跳点格式[J];石河子大学学报(自然科学版);2011年01期
中国重要会议论文全文数据库 前10条
1 苏娟;沙宪政;;复杂毛细管分布组织中氧及葡萄糖传输数学模型及计算机模拟的研究[A];中国生物医学工程学会成立30周年纪念大会暨2010中国生物医学工程学会学术大会报告论文[C];2010年
2 吴振翔;缪柏其;;MCMC方法在金融中的一点应用[A];2003中国现场统计研究会第十一届学术年会论文集(下)[C];2003年
3 张庆红;谭亦丽;韩伟;;基于扩散方程的上解与解的研究[A];数学·力学·物理学·高新技术研究进展——2006(11)卷——中国数学力学物理学高新技术交叉研究会第11届学术研讨会论文集[C];2006年
4 金蔚青;;材料制备的界面质量输运现象[A];全国第三届溶胶—凝胶科学技术学术会议论文摘要集[C];2004年
5 杨水龙;;一个化学反应扩散方程奇异行波摄动解[A];数学·力学·物理学·高新技术研究进展——2002(9)卷——中国数学力学物理学高新技术交叉研究会第9届学术研讨会论文集[C];2002年
6 袁光伟;盛志强;;大变形网格上扩散方程计算格式的设计[A];全国计算物理学会第六届年会和学术交流会论文摘要集[C];2007年
7 李西成;;经皮吸收的分数阶药物动力学模型[A];中国力学学会学术大会'2009论文摘要集[C];2009年
8 谭文长;;心肌细胞内钙火花的反常扩散[A];中国力学学会学术大会'2009论文摘要集[C];2009年
9 袁光伟;;扩散方程计算方法研究进展[A];中国工程物理研究院科技年报(2009年版)[C];2010年
10 金蔚青;;材料制备中的界面质量输运现象[A];中国空间科学学会空间材料专业委员会’2004学术交流会论文集[C];2004年
中国博士学位论文全文数据库 前10条
1 苏丽娟;分数阶扩散方程的几种数值解法[D];山东大学;2010年
2 张瀛;一维分数阶扩散方程的若干反问题研究[D];复旦大学;2012年
3 李灿;几类分数阶反常扩散方程的数值分析[D];兰州大学;2012年
4 许秋燕;二维泊松方程和扩散方程的一类显式并行算法[D];山东大学;2010年
5 李西成;分数阶可动边界问题及其在药物控释系统中的某些应用[D];山东大学;2009年
6 靳朝辉;离子交换动力学的研究[D];天津大学;2004年
7 赵强;扩散方程有限体积格式及守恒型并行差分格式研究[D];中国工程物理研究院;2011年
8 王少伟;分数阶微积分理论在粘弹性流体力学及量子力学中的某些应用[D];山东大学;2007年
9 李颖花;具周期位势或周期源的高阶扩散方程[D];吉林大学;2008年
10 金春花;具非线性源的非散度型扩散方程[D];吉林大学;2008年
中国硕士学位论文全文数据库 前10条
1 肖玲玲;非线性分数阶反应子扩散方程数值解法[D];湘潭大学;2010年
2 张方方;柱坐标系下时间分数阶波扩散方程及分析[D];山东大学;2011年
3 徐倩倩;扩散方程高阶格式的分组迭代法[D];山东大学;2011年
4 郭宾;二维分数阶反应子扩散方程的有限差方法[D];湘潭大学;2010年
5 郑生富;有损均匀传输线的时域有限差分法研究[D];重庆大学;2007年
6 刘晓蕾;局部相平衡模型(LPEM)方程数值解法的实现与分析[D];中国地质大学(北京);2007年
7 赵小红;一类非线性偏微分方程的若干有限差分格式[D];燕山大学;2008年
8 史蓉;扩散方程的势对称和Degasperis-Procesi方程的精确解[D];西北大学;2010年
9 郑营;非线性矩阵方程X+A~*X~(-n)A=Q的理论与数值解法[D];大连理工大学;2009年
10 杜萍;延迟微分方程的再生核数值解法[D];哈尔滨工业大学;2006年
中国重要报纸全文数据库 前4条
1 ;维系生命之源[N];科技日报;2006年
2 韩琦;前人所未及 后世不可无[N];中华读书报;2004年
3 贾朋群;用数字计算天气的先行者[N];中国气象报;2003年
4 ;我国古代数学成就[N];西藏日报;2002年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978