时间分数阶奇异摄动对流扩散方程的自适应网格算法.pdf
《时间分数阶奇异摄动对流扩散方程的自适应网格算法.pdf》由会员分享,可在线阅读,更多相关《时间分数阶奇异摄动对流扩散方程的自适应网格算法.pdf(6页珍藏版)》请在咨信网上搜索。
1、 年月南宁师范大学学报(自然科学版)J u n 第 卷 第期J o u r n a l o fN a n n i n gN o r m a lU n i v e r s i t y(N a t u r a l S c i e n c eE d i t i o n)V o l N o D O I:/j c n k i i s s n 文章编号:()时间分数阶奇异摄动对流扩散方程的自适应网格算法王伟,姜雨君,刘利斌(广西开放大学 文理学院,广西 南宁 ;南宁师范大学 数学与统计学院,广西 南宁 )摘要:该文提出一个求解时间分数阶奇异摄动对流扩散方程的自适应网格算法首先在均匀网格下利用L 格式对时间
2、导数进行离散然后针对空间方向,在任意非均匀网格上建立相应的迎风有限差分方法,并利用网格等分布原理设计一个适合所有时间层的自适应网格算法最后用数值实验验证算法的有效性关键词:奇异摄动;C a p u t o分数阶导数;自适应网格算法中图分类号:O 文献标志码:A本文考虑如下的时间分数阶奇异摄动对流扩散问题:L u(x,t)DtuLx,u(x,t)f(x,t),(x,t)(,)(,T,u(x,)u(x),x(,),u(,T)u(,t),t(,T,()其中,Lx,u ux xa(x)ux xb(x)u;为摄动参数;Dt为C a p u t o分数阶导数,其定义为Dtu(x,t)()tu(x,s)s(
3、ts)ds;a(x),b(x)与f(x,t)为充分光滑函数,且存在常数和使得a(x),b(x),x,()由条件()可知问题()存在唯一解,且该解当时在x处存在边界层时间分数阶偏微分方程在粘弹性力学、混沌学、生物学、金融财务等领域有着广泛的应用 大多数分数阶微分方程都无法精确求解,只有少数问题的精确解可由M i t t a g L e f f l e r等复杂函数来表示,但这些函数的计算比较困难因此分数阶偏微分方程的数值方法是当前计算数学领域中的一个热点课题S u n和Wu研究了C a p u t o时间分数阶扩散波方程的全离散有限差分方法,并在最大范数下利用能量法分析了离散格式的稳定性和收敛性
4、L i n和X u针对一类()阶C a p u t o时间分数阶扩散方程,在时间和空间方向分别构造了均匀网格下的L 有限差分格式和高精度的G a l e r k i n谱方法 G a o和S u n研究了一类带R i e m a n n L i o u v i l l e分数阶导数的次扩散方程的高阶紧致差分方法,并证明了离散格式的精度在时间方向为阶,空间方向为四阶针对阶C a p u t o分数阶算子,A l i k h a n o v 提出了一个具有阶精度的有限差分格式,并将其应用在时间分数阶扩散方程上L v和X u 提出了一个阶的有限差分格式,并给出了严格的理论分析L i和X u 提出了数
5、值求解时间分数阶扩散方程的时空谱方法,并给出了先验误差估计 M u s t a p h a和M c l e a n,针对时间分数阶扩散方程和扩散波方程,提出了间断有限元方法对于变时间分数阶次扩散方程,C h e n等 发展了一类高精度的有限差分格收稿日期:基金项目:广西自然科学基金(G X N S F AA )第一作者简介:王伟(),男,江西丰城人,讲师,研究方向:微分方程数值解通信作者简介:刘利斌(),男,湖南娄底人,教授,博士,硕士生导师,研究方向:微分方程数值解及其应用E m a i l:l i u l i b i n c o m 南 宁 师 范 大 学 学 报(自 然 科 学 版)第
6、卷式,并证明了离散格式的稳定性和收敛性当时间分数阶偏微分方程的二阶导数项包含摄动参数时,这类问题就称为时间分数阶奇异摄动问题最近,S a h o o和G u p t a 研究了问题()的参数一致的有限差分方法首先针对时间分数阶导数,在均匀网格下构造了一个L 有限差分格式;然后针对空间导数,在S h i s h k i n网格下构造了一个迎风有限差分格式;最后证明数值方法的一致收敛性姜雨君等 针对问题()的空间方向,在V u l a n o v i c B a k h v a l o v网格下严格分析了离散格式的一致收敛性本文在,的基础上进一步讨论问题()的数值算法首先,针对时间分数阶导数,在均
7、匀网格给出L 离散格式然后,在空间方向,我们在任意网格下,给出一个迎风有限差分格式,并采用网格等分布原理,设计出空间方向的一个自适应网格生成算法最后,在该自适应网格下给出相应的数值实验时间方向半离散格式的构造令t,t,并将其均匀分成N个小区间,得均匀网格Nttii t,i,N,其中网格步长 tT/N为书写方便,令(x,t)Dtu(x,t)在该网格Nt下,采用经典L 格式对C a p u t o分数阶算子Dt进行离散,得(x,tj)()j ku(x,tj)u(x,tj)ttk stk(tjs)dsdj,()u(x,tj)dj,j()u(x,t)()j ku(x,tjk)(dj,k dj,k),(
8、)其中(x,tj)为(x,t)的近似值,u(x,tj)表示问题()的半离散解,且dj,k(tjtjk)(tjtjk)t,k,j,j,N因此,我们基于()构造时间半离散格式如下:Lu(x,tj)(x,tj)Lx,u(x,tj)f(x,tj),x(,),u(x,)u(x),u(,tj)u(,tj),()其中L表示时间方向的离散算子,j,N空间离散格式及自适应网格的构造差分格式的构造令MxxxxM 为空间上任意非均匀网格,其步长为hixixi,i,M对于给定的任意网格函数viMi,定义差分算子Dvivivi hi,Dvivi vihi,DDvii(vi vihi vivi hi),其中ihihi 在
9、任意非均匀网格MxNt上,可得全离散格式LM,NUi DDUjia(xi)DUjib(xi)Ujif(xi,tj)(xi,tj),UiU(xi),xM,UjUjM,tN()其中Uji为u(x,t)在节点(xi,tj)处的数值解第期王伟,等:时间分数阶奇异摄动对流扩散方程的自适应网格算法 自适应网格算法奇异摄动问题自适应网格算法的关键技术是构造一个合适的网格控制函数M(x),使得存在一组网格xiMi满足如下网格等分布原理:xixi M(s)dsMM(s)ds,i,M()近年来,针对整数阶奇异摄动抛物方程的自适应网格算法,文献,在每个时间层上构造了一个自适应网格生成算法随后,C h e n和L i
10、 u 提出了一个适合所有时间层的自适应网格生成算法,其控制函数为M(x)m a xx,U(x),UN(x),()其中Uj(x)为经过点(xi,Uji)的分段线性插值函数,j,N,i,M基于此,本文在文 的基础上,选择()作为网格控制函数,构造出空间方向的自适应网格算法,其算法迭代步骤如下设k,设均匀网格xii/MMi为初始网格,并且选取常数C用于判断网格生成过程的终止对任意给定的k,计算出在网格点(x(k)i,tj)上方程()的数值解Uj,(k)i,i,M,j,N并令h(k)ix(k)ix(k)i,且l(k)ih(k)i m a xjN(DUj,(k)i)h(k)i m a xjN(Uj,(k
11、)iUj,(k)i)为点(x(k)i,Uj,(k)i)与点(x(k)i,Uj,(k)i)之间的最大弧长,那么总弧长为L(k)Mil(k)iMih(k)i m a xjN(DUj,(k)i)Mih(k)i m a xjN(Uj,(k)iUj,(k)i)对现有的网格进行测试,如果有m a xiMl(k)iL(k)CM,那么转到第步,否则进行第步由弧长等分布原理生成新网格xk xk xk M,它满足x(k)ix(k)i m a xx,(U(k)(x),(U(k)N(x)L(k)M令kk,并基于第步进行迭代令x(k)i为最终网格,并计算出相应的解Uj,(k)i,迭代结束数值实验为了验证本文自适应网格算
12、法的有效性,我们考虑如下时间分数阶奇异摄动对流扩散微分方程:Dtu ux x(x)ux(x)(t)u tetx(x),(x,t)(,)(,T),u(x,),x,u(,t)u(,t),t(,T()由于该问题的精确解没有给出,我们采用以下公式来计算数值解的绝对误差:EM,N m a xiMjN|UM,Ni,jUM,Ni,j|,其中UM,Ni,j表示问题()在网格MxNt下的数值解,UM,Ni,j是将网格MxNt下进行二等分加密后,在网格节点(xi,tj)上的的数值解相应的收敛阶的计算公式为 南 宁 师 范 大 学 学 报(自 然 科 学 版)第 卷M l o gEM,NEM,N当 k(k,),M
13、j(j,)时,表给出了N 时,由本文的自适应网格算法得到的数值结果从表中可以看出,本文的自适应网格算法在空间方向的收敛阶为一阶表自适应网格方法计算得到的数值结果(N ,)M M M M M M eM E E E E E E pM eM E E E E E E pM eM E E E E E E pM eM E E E E E E pM eM E E E E E E pM eM E E E E E E pM 为了进一步展示出自适应网格生成算法的迭代过程,当 ,N ,M 时,图画出了每一次迭代过程中网格点的移动过程,图画出了相应的数值解的图像从图和图可以看出问题()的解在x处存在边界层图自适应网格
14、算法的网格移动过程(,N,M)第期王伟,等:时间分数阶奇异摄动对流扩散方程的自适应网格算法 图问题()的数值解(,N,M)结束语本文给出了一类时间分数阶奇异摄动对流扩散方程的自适应网格算法首先在均匀网格上利用L 格式对时间分数阶导数进行离散然后基于任意非均匀网格,对空间导数采用迎风有限差分方法进行离散,得到相应的全离散格式接着利用网格等分布原理设计自适应网格网格算法值得指出的是本文的自适应网格算法能进一步推广到其他分数阶奇异摄动问题,甚至带弱奇性的时间分数阶偏微分方程参考文献:K i l b a sAA,S r i v a s t a v aH M,T r u j l l oJJ T h e
15、o r ya n da p p l i c a t i o n so f f r a c t i o n a l d i f f e r e n t i a l e q u a t i o n sM E l s e v i e rS c i e n c eBV,Am s t e r d a m,H a j i p o u rM,J a j a r m iA,B a l e a n uD A ne f f i c i e n t n o n s t a n d a r d f i n i t ed i f f e r e n c e s c h e m e f o r a c l a s s o
16、 f f r a c t i o n a l c h a o t i cs y s t e m sJ J o u r n a l o fC o m p u t a t i o n a l a n dN o n l i n e a rD y n a m i c s,():H a j i p o u rA,H a j i p o u rM,B a l e a n uD O nt h ea d a p t i v es l i d i n gm o d e c o n t r o l l e r f o r ah y p e r c h a o t i c f r a c t i o n a l o
17、 r d e r f i n a n c i a l s y s t e mJ P h y s i c a,A:S t a t i s t i c a lm e c h a n i c sa n d i t sa p p l i c a t i o n s,:L o n e s c uC,L o p e sA,C o p o tD,e ta l T h er o l eo f f r a c t i o n a l c a l c u l u s i nm o d e l l i n gb i o l o g i c a lp h e n o m e n a:ar e v i e wJC o
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 时间 分数 奇异 摄动 对流 扩散 方程 自适应 网格 算法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。