分数阶偏微分方程的高阶差分方法及其应用研究

来源 :东南大学 | 被引量 : 0次 | 上传用户:amoyzhu
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
在过去的十几年里,随着分数阶偏微分方程的广泛应用,证实了带分数阶导数的模型比传统的整数阶导数模型能更精确的描述科学与工程领域中的系统现象.然而,分数阶模型的解析解很难得到.所以,分数阶偏微分方程的数值方法研究已成为当前应用数学领域一个重要的分支.本文主要研究几类时间分数阶偏微分方程、空间分数阶非线性Schr(o)dinger方程、变时间分数阶偏微分方程及相关应用。  首先考虑两类空间带变系数的时间分数阶偏微分方程的差分方法.对带Dirichlet边界条件的空间变系数反常扩散方程,结合时间分数阶导数的L1离散和空间导数的四阶逼近.建立了紧差分格式.用能量方法分析了紧差分格式的无条件稳定性和收敛性,对带Neumann边界条件的空间变系数反常扩散方程,利用降阶法构造了Box型格式,并相应的给出了理论分析结果.通过数值试验验证了理论分析的精度。  对一维分数阶变系数Cattaneo方程.结合了对一阶时间导数和分数阶时间导数的Crank-Nicolson型离散,构造了Crank-Nicolson型紧格式.通过能量方法,分析了紧格式在H1半范下的稳定性和收敛性.收敛阶为O(τ2-α+ h4),其中τ和h分别是时间和空间方向的步长,α是时间分数阶导数的阶数.对二维问题,结合对时间导数的离散,选择适当的小量项建立了不降低时间方向精度的Crank-Nicolson型紧交替方向差分格式.一维问题和二维问题的数值算例证明了所提差分格式的有效性。  其次,对α∈(0,1)阶Caputo导数,在半网格点提出了中点公式.并分析了公式的截断误差.结合中点公式对分数阶导数的离散及二阶中心差商对拉普拉斯算子的离散,对反常扩散问题构造了时空均是二阶的分数阶Crank-Nicolson格式(FCN).对二维问题,利用添加小量项的方法,建立了分数阶Crank-Nicolson交替方向隐格式,并通过数值方法验证了格式在时间方向的精度是1+α。通过数值试验,与隐格式(INM)(Zhuang et al.,SIAM J.Numer.Anal.,46(2008),pp.1079-1095)与Crank-Nicolson型(CN)格式(Zhang et al.,SIAMJ.Numer.Anal.,49(2011),pp.2302-2322)进行比较,结果显示FCN格式对任意的α在时间方向都是二阶收敛的.而当α接近1时,已有的格式的精度较低.对α∈(1,2)阶Caputo导数,构造二阶的逼近公式并分析其截断误差.结合二阶逼近公式与空间二阶中心差商,对扩散波问题提出了时空均是二阶收敛的差分格式.此外.对空间导数运用紧算子,建立了空间四阶的紧格式.数值算例验证了所提差分格式在时间方向和空间方向的精度。  研究非线性空间分数阶Schr(o)dinger方程的数值解法.利用已有的二阶中心差分格式(Celik和M.Duman,J.Comput.Phys.,231(2012),pp.1743-1750)对α∈(1,2]阶Riesz导数构造紧算子并证明算子具有四阶精度.对一维的非线性空间分数阶Schr(o)dinger方程,建立了三层线性化紧格式并给出了格式的唯一可解性与收敛性的理论分析.对二维问题,结合紧算子对空间导数的逼近,提出了紧交替方向隐格式,并证明了格式的唯一可解性与收敛性,收敛阶为O(τ2+ h4),其中τ是时间步长,h=max{h1,h2},h1,h2分别是空间x方向和y方向的网格步长.数值试验验证了理论分析的精度同时利用CPU时间说明了新的紧算子能大大提高格式的精度而不带来额外的计算量。  最后研究时间变分数阶偏微分方程的应用.借助于时间变分数阶导数,揭示反常扩散及扩散波对应不同的阶数在不同的时间段内的现象.通过调节依赖于时间和空间变量的时间分数阶导数的阶数,将波传播问题在边界附近逐渐转变为扩散问题,使得波的能量在边界附近逐渐耗散从而有效地防止波的反弹.把变分数阶导数算子应用于Burgers方程的数值求解,利用导数的扩散性质,有效的消除了数值解的震荡。
其他文献
在古典风险过程模型中,对破产概率求解的方法可谓多种多样,尤其当索赔额及索赔时间间隔满足一定的分布时,我们会发现一些非常规的方法比一般方法在解决破产概率问题方面显得
该文首先讨论了算子正则性的一般结果.着重考查了算子正则性与绝对值的关系,即绝对值的存在性问题,并得到相应的结论:当空间具有某种序完备性时,两者等价;而当空间不具备序完
2010年9月4日,由冀中能源峰峰集团组织的“冀峰2×300 MW煤矸石综合利用发电项目低热值燃料评估报告评审会”在北京召开,煤炭洗选加工、资源综合利用、电厂运行、电厂设计等
森林是地球上最大的陆地生态系统,具有重要经济价值、生态价值和社会价值。它不但为人类提供了大量的林产品,还对调节气候、水分平衡、净化大气、土壤肥力改善,以及居民的游
该文尝试给出有限群G的极大子群M的h-完备的定义,即称C为M的h-完备,若C是M的θ-完备,且有C=G或存在G的一个子群D使得D不是M的θ-完备,但有C
2010年12月6-9日,“第十届中国国际现代化铁路技术装备展览会”在北京国家会议中心举办。来自世界十多个国家和地区的180多家铁路企业参展。本次展览会重点围绕高速铁路工程
该论文由六节组成.第一节给出了小波子空间上抽样定理的发展状况;第二节呈现出已知结果和准备情况.第三节至第六节是本论文的主要内容.具体为:第三节构造了一类正交基插值尺
随机加权方法,也称之为Bayesian bootstrap,它有别于Bootstrap方法.在第一章中,我们介绍了Randomization Tests,Jackknife和Bootstrap等三种再抽样方法.在第二章中,我们发现
在这篇文章里,我们首先考察了一类p-Laplace方程的Dirichlet问题(公式略)其中Ω( )R是边界光滑的有界区域,1
教学中,教师起着重要的作用。教学的方向、内容、方式方法、进程、结果等都离不开教师的引导。学生的思想觉悟、学习动机、学习方法、以及学习效果,都要受到教师的影响。可见,教