分数阶微分方程的谱方法和间断Galerkin方法研究
【摘要】摘要：分数阶微分方程以其带有全局性和记忆性特征而被应用于许多工程问题的数值模拟中,应用发展的需要使得分数阶微分方程的求解成为紧迫和重要的研究课题。多区域谱方法和间断Galerkin法由于仅在局部单元独立构造方程的逼近格式,不同的单元通过惩罚边界或者边界上的数值流函数来进行信息交换,在基函数的选取和网格的剖分方面具有很好的灵活性和较高的精度,已在传统微分方程的求解中发挥了重要的作用,本文将主要基于多区域谱方法和间断Galerkin法对分数阶微分方程的数值求解进行研究,以提出适应于任意网格剖分的高精度、高效率算法。主要研究成果如下：(1)用间断Galerkin法方程对空间进行离散、有限差分法对时间进行离散,构造了时间分数阶扩散方程的有限差分-间断Galerkin法,证明了空间半离散格式的L2稳定性和最优收敛性,用数值算例验证了理论结果。(2)基于多区域谱方法的思想,为了达到较高的求解精度而同时尽量不增加计算量,提出了一种计算分数阶导数记忆项的混合方法,构造了求解时间分数阶微分方程的高阶多区域谱方法。根据分数阶导数的“短记忆原理”,将具有“无限步数”的多步法简化为具有固定步数的线性多步法,利用特征值分析法给出了方法的稳定性区域,通过数值算例分析了方法的收敛性。数值结果显示,这种方法具有N+1-α(N为基函数的最高次数)阶的精度,弥补了有限差分-间断Galerkin法时间离散精度不高的弱点。(3)用多项式逼近方法构造了单区域内分数阶导数的逼近、用分片连续多项式逼近方法构造了分数阶导数在多区域内的逼近,利用正交多项式的性质给出逼近方法的误差分析,证明了当N-[α]为奇数时,所给方法具有N+1-[α]阶精度；当N-[α]为偶数时,所给方法具有N+1-α阶精度。通过添加罚函数项的方式为方程附加单元边界条件,分别构造了分数阶对流方程和扩散方程的罚函数谱方法,用矩阵特征值方法分析了所构造的罚函数谱方法的稳定性,并用数值算例验证了方法的收敛性。(4)研究了带有分数阶Laplace算子的对流扩散方程的间断Galerkin方法。为了保证间断Galerkin法的相容性和高精度,将α(1<α<2)阶Laplace算子分解为两次一阶导数和2-α次分数阶积分的复合形式,通过辅助变量得到了低阶偏微分方程组,并对方程的扩散项采用交替方向数值流函数,对非线性项采用单调数值流函数,推导出了分数阶对流扩散方程的局部间断Galerkin法数值格式,证明了所给算法的L2稳定性；针对分数阶扩散方程,证明了所给算法具有最优收敛阶；对带非线性对流项的分数阶对流扩散方程证明了所给算法具有N+1/2阶精度。以分数阶扩散方程和分数阶Burgers方程为算例,验证了方法的稳定性和收敛阶。(5)利用有限差分方法对时间域进行离散、配置方法对空间域进行离散,针对一类带广义分数阶导数的二维时间分数阶扩散方程构造了有限差分-配置法,并基于ADI给出了快速计算方法；证明了全离散格式对于初值和源项的无条件稳定性,通过数值算例验证了方法的稳定性、分析了方法的收敛性和精度。
【作者】徐勤武；
【导师】郑洲顺；JanS.Hesthaven；
【作者基本信息】中南大学，数学，2014，博士
【关键词】分数阶微分方程；谱方法；间断Galerkin法；最优收敛阶；
【参考文献】
[1]李艳红.半夏白术天麻汤加味联合针刺治疗老年原发性高血压（痰瘀阻络型）的临床研究[D].云南中医学院,中西医结合临床,2012,硕士.
[2]陈乔.涡流缓速器电磁优化设计研究[D].华中科技大学,电气工程,2013,硕士.
[3]王偕蓉.00Cr17Ni14Mo2超低碳不锈钢热穿孔管裂纹缺陷分析[J].四川冶金.1989(03)
[4]王正,张喻,何盛.适张木工带锯条金相组织与硬度的试验分析[J].木工机床.2007(01)
[5]杨琛.再生材料在装置艺术中的应用研究[D].云南艺术学院,设计艺术学,2013,硕士.
[6]杜沙沙.外啮合齿轮泵模型建立以及振动和噪声研究[D].辽宁工程技术大学,机械电子工程,2012,硕士.
[7]魏旭旭.CYP2C19基因分型方法的建立及其在氯吡咯雷个体化用药的初步应用[D].西北大学,生物化学与分子生物学,2013,硕士.
[8]张跃.蜂胶降血糖有效部位化学成分研究[D].吉林农业大学,野生动植物保护与利用,2012,硕士.
[9]张智强.提升庄河市志愿服务水平研究[D].大连理工大学,MPA（专业学位）,2012,硕士.
[10]陆紫微.非负重状态下3.0T MR T2mapping及IDEAL成像在评价椎间盘退变中的应用[D].苏州大学,影像医学与核医学,2013,硕士.
[11]王剑平,张云生,张果,张晶.实时系统时间并行模型确定性分析[J].控制与决策,2012,12:1859-1863+1869.
[12]罗媛.以环境行为学理念探讨医院病房楼设计[D].湖南大学,建筑学,2011,硕士.
[13]孙宝民.运动想象疗法对脑卒中偏瘫患者日常生活能力的影响[D].青岛大学,康复医学与理疗学,2011,硕士.
[14]张洁.鲨鱼肝脏中与肝再生相关的microRNAs研究[D].浙江理工大学,生物化学与分子生物学,2013,硕士.
[15]王双.我国公司董事勤勉义务研究[D].西南大学,民商法学,2014,硕士.
[16]田胜利.城市轨道交通电力监控系统结构的标准化问题思考[J].城市轨道交通研究.2004(04)
[17]李倩.严羽“以禅喻诗”研究[D].华东师范大学,文艺学,2013,硕士.
[18]李燕.艺术歌曲《春晓》声乐教学中的审美特征探析[J].音乐大观,2014,03:67-68.
[19]章骁程.用反向有限元法从小冲杆试验中获得材料的断裂韧性参数[D].华东理工大学,化工过程机械,2014,硕士.
[20]石华.与时俱进:税收信息化明确战略方向[J].中国税务,2002,03:10-14.
[21]刘文干.三株红壤高效溶磷菌的分离、鉴定、溶磷特性及其对花生促生效应的研究[D].南京农业大学,微生物学,2012,硕士.
[22]王科.基于TiO_2薄膜染料敏化太阳能电池的制备及性能优化[D].中央民族大学,2013.
[23]吕行.以表演为重心：徐渭南戏理论研究[D].湖北民族学院,文艺学,2014,硕士.
[24]韩明正.采煤塌陷矿区土地整理模式研究[D].中国农业大学,2004.
[25]张旭.基于BP神经网络和模态曲率法的拱桥损伤识别研究[D].兰州理工大学,桥梁与隧道工程,2013,硕士.
[26]常乐.郑州市中学篮球队课余训练现状及对策研究[D].哈尔滨体育学院,体育教育训练学,2013,硕士.
[27]赵少磊.唐山矿西南区深部瓦斯赋存特征及聚集规律研究[D].中国矿业大学,矿产普查与勘探,2014,硕士.
[28]谢佳.译者主体性下《庄子》核心词“道” 的英译研究[D].湖南师范大学,英语语言文学,2014,硕士.
[29]童锟.dBASE—Ⅲ在总调度计算机管理中的应用[J].冶金自动化,1987,02:52-54.
[30]魏连启.酚醛树脂纳米复合材料及其在摩阻材料中的应用研究[D].武汉理工大学,材料学,2004,硕士.
[31]陈玉波.甘蓝型油菜A1上主效含油量QTL OilAl与千粒重QTL SWA1的精细定位[D].杭州师范大学,植物学,2013,硕士.
[32]程度海.变粗放为精益[J].医药产业资讯,2005,22:56-57.
[33]陈显坤.急倾斜松软煤层综合机械化开采技术研究[D].西安科技大学,矿业工程,2013,硕士.
[34]邓聪迩.6-烷氨基/烷氧基-2-丙硫基-8-氮杂嘌呤核苷的合成及抗血小板凝聚活性评价[D].北京化工大学,化学,2013,硕士.
[35]丁洋洋.联合国安理会第1973号决议的国际法探析[D].贵州大学,国际法,2014,硕士.
[36]姚永玲.网络新闻评论及其功能[D].中央民族大学,新闻学,2004,硕士.
[37]余品能.离散Y变换及其快速算法[J].石油地球物理勘探,1994,05:558-565+668.
[38]沈瑜平.冷轧厂天车(行车)定位及智能导航系统[J].冶金自动化,2015,01:80-84.
[39]戴青原.内脏脂肪素对氧化应激HUVECs的作用及其分子机制研究[D].昆明医科大学,外科学,2014,博士.
[40]王迟宇.柴油机热平衡数值仿真与试验研究[D].浙江大学,2007.
[41]陈楠纬,孙水裕,任杰,曾佳俊,刘敬勇.咖啡渣燃烧特性及动力学研究[J].环境科学学报.
[42]张德锋,郭玉霞,宋志刚.数据挖掘技术[J].航空计算技术,2005,03:76-79.
[43]姚传好.燃料电池用Ti_4O_7的合成及其载Pt催化剂的性能研究[D].北京工业大学,2012.
[44]江丽.出口、集聚与产业效率[D].浙江大学,国际商务（专业学位）,2013,硕士.
[45]江杰.内蒙古居民收入增长与经济增长关系的研究[D].内蒙古农业大学,技术经济及管理,2013,硕士.
[46]闫雪菲.城市社区公共服务机制创新研究[D].华中师范大学,宪政与法治,2014,硕士.
[47]杨涛,李静,潘泉,张艳宁.基于场景模型与统计学习的鲁棒行人检测算法[J].自动化学报,2010,04:499-508.
[48]王磊.新课程促进中学化学教师专业发展研究[D].河北师范大学,学科教学,2014,硕士.
[49]马自伟,葛华才.氧化石墨烯/壳聚糖复合材料的制备及吸附研究进展[J].材料导报,2014,S1:60-63.
[50]樊国强.纳米结构金属的表面等离效应在有机太阳能电池中的应用[D].苏州大学,2013.
