多项时间分数阶扩散方程类Carey非协调元的误差分析.pdf
《多项时间分数阶扩散方程类Carey非协调元的误差分析.pdf》由会员分享,可在线阅读,更多相关《多项时间分数阶扩散方程类Carey非协调元的误差分析.pdf(5页珍藏版)》请在咨信网上搜索。
1、第 43 卷 第 2 期许昌学院学报Vol.43.No.2 2024 年 3 月JOURNAL OF XUCHANG UNIVERSITYMar.2024收稿日期:2023-09-28作者简介:马国锋(1971),男,河南许昌人,副教授,硕士,研究方向:有限元方法及应用.文章编号:1671-9824(2024)02-0007-05多项时间分数阶扩散方程类 Carey 非协调元的误差分析马国锋(许昌学院 数理学院,河南 许昌 461000)摘 要:基于 L1 全离散格式,针对具有 Caputo 导数的二维多项时间分数阶扩散方程给出了类 Carey 非协调有限元方法.利用该单元的特殊性质和分数阶导
2、数巧妙的处理技巧导出了 L2模和 H1模意义下的最优误差估计.关键词:多项时间分数阶扩散方程;类 Carey 元;全离散格式;最优误差估计中图分类号:O242.21 文献标识码:A 考虑如下的多项时间分数阶扩散方程P,1,2,m(Dt)u-u=f(X,t),(X,t)(0,T,u(X,t)=0,(X,t)(0,T,u(X,0)=u0(X),X ,(1)其中 R2是有界的凸多边形区域,是 的边界,X=(x,y),u0(X),f(X,t)为已知适当光滑的函数,算子 P,1,2,m(Dt)为P,1,2,m(Dt)=Dt+mi=1liDit,li 0,m N+,0 1 2 m 1,Dt是关于 t 的左
3、 阶 Caputo 导数算子,其定义为Dtu(X,t)=1(1-)t0u(x,s)s1(t-s)ds,这里()是 Gamma 函数.分数阶偏微分方程可以用来模拟工程、物理、生物等科学领域的许多现象.研究者发现分数阶微分算子具有空间全域相关性和历史记忆性,因此分数阶微分方程非常适合用来描述现实生活中具有记忆和遗传特性的问题,如电容理论、电解化学、半导体物理、粘弹性系统、生物数学及统计力学等等.由于它的解析解不易求得,因而分数阶微分方程的数值算法研究备受关注.例如,文1考虑了多项时间分数阶扩散方程的谱方法,给出了全离散格式下的误差分析.文2对多项时间分数阶对流扩散方程构造了一类显-隐和隐-显差分格
4、式,并证明了格式的无条件稳定性和收敛性.文3对多项时间分数阶扩散方程提出了一种高阶数值求解方法,并对其稳定性和收敛性进行了严格的证明.文4对单项时间分数阶扩散方程进行了非协调新混合元方法研究,得到了原始变量和中间变量的最优误差估计.文5对问题(1)利用非协调EQ1rot元进行了收敛性分析.文6和7分别针对单项和多项时间分数阶扩散方程进行类 Wilson 非协调逼近并导出了超收敛结果.类 Carey 元是非协调三角形元,关于整数阶偏微分方程的类 Carey 元方法已有一些讨论,如文8-10将类 Carey 元分别应用于双曲积分微分方程、非线性 Sobolev-Galpern 型湿气迁移方程和 P
5、oisson 方程,得到了有限元解的超逼近性质和整体超收敛结果.但据我们所知关于问题(1)的三角形非协调元逼近目前许昌学院学报2024 年 3 月还未见报道.本文将类 Carey 非协调有限元应用到二维多项时间分数阶扩散方程进行数值逼近.基于分数阶导数巧妙的处理技巧及文3中给出的该单元性质导出了 L2模和 H1模意义下的最优误差估计.1 类 Carey 元构造和重要引理设 Jh是 的一族直角三角形剖分,K Jh,设其两条边分别平行于 x 轴和 y 轴,并且所有水平边和竖直边分别相等(记为 GATM 三角形网格).与文5 类似定义类 Carey 非协调元空间 Vh及空间上相应的插值算子.设Vh
6、Vh为 Jh上线性三角形有限元空间,K Jh,记插值算子 Ih:H2()Vh.引理 1 3假设 H10()H3(),对任意 vh Vh,有KKnvhds Ch23vhh,(2)(h(-Ih),hvh)Ch23vhh Ch 3vh0,(3)其中,h=(K21,K)12,(w1(x,y),w2(x,y)=KKw1(x,y)w2(x,y)dxdy,h是分片求导算子,C是不依赖空间步长 h 和时间步长 的正常数,不同的地方可以不同,n 是 K 上的外法线方向向量.问题(1)的变分形式是:求 u(X,t):(0,T H10(),满足(P,1,2,m(Dt)u,v)+(u,v)=(f,v),v H10()
7、,u(X,0)=u0(X),X .(4)2全离散逼近格式和稳定性分析假设 0=t0 t1 t2 tN-1 B1 BN-1 0.当 1 k n-1 时有 Bnk 0.问题(4)的全离散逼近格式为:给出 Un-1 Vh,求 Un Vh,使得(P,1,2,m(Dt)Un,vh)+(hUn,hvh)=(f n,vh),vh Vh,U0=Ihu0(X),X .(6)为方便起见,令 Rn=P,1,2,m(Dt)un-P,1,2,m(Dt)un.若 utt(X,t)L2(),t (0,T,根据文10 有Rn0 Cmax0tTutt(X,t)02-.(7)为了给出有限元解的稳定性,利用文6 中的思想,先给出如
8、下引理.引理 2设 nNn=0是 上的函数列,有8第 43 卷第 2 期马国锋:多项时间分数阶扩散方程类 Carey 非协调元的误差分析(n,nk=0Bnkk)=12(B0n20+n-1k=0Bnkk20-n-1k=0Bnkk-n20).引理 3按照 b,N-1,BN-1的定义,则b,N-1(1-)T-,BN-1,这里 =(1-)T-+mi=1(2-)(2-i)Ti(1-i)li.基于文5 中的方法可得逼近格式的稳定性如下:定理 1设 Un 是方程(6)的解,则有Un2h U020+C maxt0,Tf(X,t)20.3最优误差估计定理 2设 un,Un分别是(4)和(6)在 t=tn的解.对
9、任意 t 0,T,若 u H3(),ut H2(),utt L2(),那么对任意 1 n N,有un-Un0=O(h2+2-).证明 对于 vh Vh,由(4)和(6)可得误差方程:(P,1,2,m(Dt)(un-Un),vh)+(h(un-Un),hvh)=KKunvhds-(Rn,vh).(8)为方便起见,定义 un-Un=(un-Ihun)+(Ihun-Un)=n+n.在(8)中令 vh=n,有 (P,1,2,m(Dt)n,n)+n2h =-(P,1,2,m(Dt)n,n)-(n,hn)+KKunnnds-(Rn,n).(9)根据引理 2,得(P,1,2,m(Dt)n,n)=-2(2-)
10、(B0n20+n-1k=0Bnkk20-n-1k=0Bnkk-n20).(10)注意到 0=0,由式(9)和(10),我们得到 B0n20+Bn-1n20-nk=1Bnkk-n20+2(2-)n2h =-n-1k=0Bnkk20-2nk=0Bnk(k,n)-2(2-)(n,hn)-2(2-)(Rn,n)+2(2-)KKunnnds=5j=1Wj.(11)接下来估计(11)式的右端,定义(t)L(Hm()=maxt0,T(t)m,并由tn-k0=1tn-ktn-k-1tdt0 Ch2utL(H2(),得到W2=-2n-1k=0Bk(tn-k,n)2n-1k=0Bktn-k0n0 Ch2utL(H
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 多项 时间 分数 扩散 方程 Carey 协调 误差 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。