二维非线性四阶分数阶波动方程的BDF2-WSGI有限元算法.pdf
《二维非线性四阶分数阶波动方程的BDF2-WSGI有限元算法.pdf》由会员分享,可在线阅读,更多相关《二维非线性四阶分数阶波动方程的BDF2-WSGI有限元算法.pdf(9页珍藏版)》请在咨信网上搜索。
1、Advances in Applied Mathematics 应用数学进展应用数学进展,2024,13(4),1217-1225 Published Online April 2024 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2024.134112 文章引用文章引用:刘心愿.二维非线性四阶分数阶波动方程的 BDF2-WSGI 有限元算法J.应用数学进展,2024,13(4):1217-1225.DOI:10.12677/aam.2024.134112 二维非线性四阶分数阶波动方程的二维非线性
2、四阶分数阶波动方程的BDF2-WSGI有限元算法有限元算法 刘心愿刘心愿 内蒙古大学数学科学学院,内蒙古 呼和浩特 收稿日期:2024年3月9日;录用日期:2024年4月5日;发布日期:2024年4月12日 摘摘 要要 本文主要研究了二维非线性四阶分数阶波动方程的有效数值算法。通过结合二阶本文主要研究了二维非线性四阶分数阶波动方程的有效数值算法。通过结合二阶BDF2-WSGI时间离散格时间离散格式与有限元方法对二维非线性四阶分数阶方程进行求解。首先,引入辅助变量,将分数阶四阶波动问题式与有限元方法对二维非线性四阶分数阶方程进行求解。首先,引入辅助变量,将分数阶四阶波动问题转化为低阶耦合方程,然
3、后利用转化为低阶耦合方程,然后利用Riemann-Liouville分数阶积分对所得方程进行积分,最后使用分数阶积分对所得方程进行积分,最后使用WSGI逼逼近公式逼近分数阶积分,形成二阶近公式逼近分数阶积分,形成二阶BDF2有限元格式。本文给出了详细的数值算法,并通过一个二维算例有限元格式。本文给出了详细的数值算法,并通过一个二维算例进行了数值试验,验证了算法的有效性和收敛性。进行了数值试验,验证了算法的有效性和收敛性。关键词关键词 二维非线性四阶分数阶波动方程,有限元方法,二维非线性四阶分数阶波动方程,有限元方法,BDF2格式,格式,WSGI公式公式 BDF2-WSGI Finite Ele
4、ment Algorithm for a Two-Dimensional Nonlinear Fourth-Order Fractional Wave Equation Xinyuan Liu School of Mathematical Sciences,Inner Mongolia University,Hohhot Inner Mongolia Received:Mar.9th,2024;accepted:Apr.5th,2024;published:Apr.12th,2024 Abstract This article mainly studies effective numerica
5、l algorithms for two-dimensional nonlinear fourth-order fractional wave equations.We combine the second-order BDF2-WSGI time discretization scheme with the finite element method to solve two-dimensional nonlinear fourth-order fraction-刘心愿 DOI:10.12677/aam.2024.134112 1218 应用数学进展 al equations.Firstly
6、,introducing auxiliary variables transforms the fractional fourth-order wave problem into a low-order coupled equation.Then,the Riemann-Liouville fractional integration is used to integrate the resulting equation.Finally,the WSGI approximation formula approximates the fractional integration,forming
7、a second-order BDF2 finite element scheme.This article pro-vides a detailed numerical algorithm and conducts numerical experiments on a two-dimensional example to verify the effectiveness and convergence of the algorithm.Keywords Two-Dimensional Nonlinear Fourth-Order Fractional Wave Equation,Finite
8、 Element Method,BDF2 Scheme,WSGI Formula Copyright 2024 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 本文主要对以下二维非线性四阶分数阶波动方程进行了数值研究()()()()()()()2000w,w,0,0,w,0,w,0w,w,
9、0w,w.CttD uuf ugttTuutTuuuu+=(1)其中12,0T,二维空间2R。()f u是非线性项,()w,gt是源项,0CtD u是 Caputo 分数阶导数()()()21020w,1d.2tCtutD utsss=分数阶波动模型1-7作为分数阶偏微分方程的一个重要模型,在化学、生物和物理领域扮演着重要的角色,受到了专家学者们的广泛关注,越来越多的人开始研究分数阶波动模型的解析解和数值解。在4中,Dahaghin 和 Hassani 采用一种基于广义多项式的优化方法来求解分数阶波动模型。在7中,Ye等人研究了分数阶波动方程的紧致差分格式,并讨论了差分解的唯一可解性。近年来,
10、具有空间四阶导数的分数阶波动方程8 9 10 11 12由于其复杂的结构,精确解难以求解,很多学者使用不同的数值方法进行有关这类方程的数值研究,主要方法有有限元法、谱方法、有限体积法、有限差分法等。在11中,Huang 等人研究了具有空间四阶导数的时间分数阶非线性波动方程的两种线性化格式,并严格证明了这两种格式的收敛性和无条件稳定性。在12中,Li 等人在空间上采用参数五次样条的方法来求解四阶导数的分数阶波动方程,并验证了此方法的可解性,稳定性和收敛性。在13中,Liu 等人提出一种基于混合有限元方法的双网格算法,来求解具有 Caputo 型时间分数阶导数的非线性四阶反应扩散方程。在14中,W
11、ang 等人针对具有分数阶导数的四阶扩散波模型,提出了一种改进 L1 格式的混合单元算法,并分析了其误差估计。在15中,Wang 等人在空间方向上使用 Galerkin 混合有限元方法,在时间方向上使用具有 WSGD 公式的二阶 格式研究了二维非线性四阶时间分数阶波动模型。还有诸多研究,我们不能一一列举,但我们可以看到具有空间四阶导数的分数阶波动方程具有很好的研究空间。求解分数阶模型的重要步骤之一是对分数阶导数项进行处理,许多学者使用 WSGD 16 17逼近公Open AccessOpen Access刘心愿 DOI:10.12677/aam.2024.134112 1219 应用数学进展
12、式对分数阶导数项进行处理,在18中 Feng 研究了一类二维非线性时间分数阶耦合亚扩散方程组模型,在时间上采用 Crank-Nicolson 格式和 WSGD 逼近,在空间上采用有限元方法。近年来,WSGI 逼近公式作为一个求解分数阶偏微分方程的重要方法也受到了学者的广泛研究。在19中,Wang 等人在时间上采用 WSGI 近似的 Crank-Nicolson 格式,在空间上采用有限元方法研究了二维非线性四阶时间分数阶波动模型的算法。在20中,Cao 等人结合二阶 Crank-Nicolson 格式和 WSGI 逼近与有限元方法对多维时间分数阶波动方程进行求解。在本文中,我们将结合 BDF2-
13、WSGI 格式研究二维非线性四阶分数阶波动方程。为此,我们引入了一个非线性辅助变量,将分数阶四阶波动模型转换成一个低阶耦合系统,对其中一个方程进行积分,然后将时间上具有 WSGI 近似的 BDF2 格式与空间上的有限元方法相结合,得到完全离散格式。我们详细展示了数值算法以及如何执行数值计算。最后,通过数值算例验证了方法的可行性,并给出了收敛结果。本文的其余部分安排如下:在第 2 节中,我们构造了二维非线性四阶分数阶波动方程完全离散格式。在第 3 节中,我们提供了一个详细的数值算法。在第 4 节中,给出了一个数值例子来验证我们算法的有效性。最后,在第 5 节中,我们将对本文做一个简要的总结。2.
14、混合元全离散格式混合元全离散格式 在本节,我们构造式子(1)的全离散混合元格式。通过引入辅助变量()vuf u=,方程(1)可以改写为如下低阶耦合系统()()()()(0w,w,0,0,w,0,0,w,0,.CtD uvgttTvuf utTuvtT+=+=(2)为了构造全离散格式,对时间区间0,T进行均匀剖分使得010ntttT=,时间节点ntn=,这里TN=为时间不长,N 是一个正整数。为简便,我们记101tuuDu=,()nnuu t=。为了分析,我们需要引入 Riemann-Liouville 分数阶积分,定义如下()()()0101d,ttf sfsts=(3)其中01 1=,然后通
15、过此积分来积分(2)的第一个式子,以得到一下耦合系统()()()()()()(000w,w,w,w,0,0,w,0,.tttutvgtutTvuf utT +=+=(4)为了求解(4),我们需要引入 WSGI 近似来近似算子0t。()1111000,nnnkntktkuuEuE+=+(5)其中()2EO=,并且()()()()()0011,1,1,222kkkk =+(6)这里的系数()k定义如下()()()()()0111,1,1,1.kkkkkkk =+(7)刘心愿 DOI:10.12677/aam.2024.134112 1220 应用数学进展 我们使用 WSGI 近似和 BDF2 格式
16、得到以下时间方向的半离散格式:情况情况 1:1n ()111110001111134,2,nnnnnnttnnnnuuuvguEvuf uE+=+=(8)情况情况 2:0n=()111100011111,tttDuvguEvuf uE+=+=(9)其中()()()112nnnf uf uf u+=,()112nncEEEO+=+=,()()1121nncttnEDuu tO+=,()()()1111nnnEf uf u t+=。接下来我们给出(8)和(9)的弱格式,(8)和(9)的两个式子分别和10,H 做内积,然后对空间区域积分,我们得到:情况情况 1:1n ()()()()()()()()
17、()111110001111134,2,.nnnnnnttnnnnuuuvguEvuf uE+=+=(10)情况情况 2:0n=()()()()()()()()()()111100011111,.tttDuvguEvuf uE=+=(11)取有限元空间10hVH并得到如下全离散格式:情况情况 1:1n ()()()()()()()()()11111100011134,2,0,.nnnnnnnhhhhkhhkhhhhhkknnnhhhhhhhhuuuvguVvuf uV +=+=+=(12)情况情况 2:0n=()()()()()()()()()()()111111000111,0,.thhkh
18、hkhhhhhhhkkhhhhhhhhDuvguEVvuf uV =+=(13)3.算法实现算法实现 通过取有限元空间0|MhhhiiiVvvv=,其中 M 是空间单元的数量,我们有()()0Mhiiiutu t=,()()0Mhjjjvtvt=。将这些代入到(12)中并且令hm=,hl=,我们可以得到 刘心愿 DOI:10.12677/aam.2024.134112 1221 应用数学进展 ()()()()()()()()()()()()()00000000,0,1,2,0,0,1,2,.MMiimtjjmtmmtijMMMjjljiliiljiiu tvtg tumMvtutfu tlM
19、=+=+=(14)现在,考虑边界条件,我们得到(14)的半离散有限元格式的矩阵形式()()()()()()()000,0,htththhhtttttt=+=AuBvGuBuAvF u (15)其中()()()()11,imMM =A,()()()()11,imMM=B,()()()()T1100,MMhiiiiMiitfu tfu t=F u,()()T11,hMtvv=v,()()T11,hMtuu=u,()()()()()()T11,Mtg tg t=G,()()()T00101,Muu=u。接下来我们给出(15)的全离散格式 111100011134,20,nnnnnhhhthtnnnh
20、h+=+=uuuABvGuBuAvF (16)其中()()1111001nnnnkthhkhk+=+vvv,()T11,kkkhMuu=u,()T11,kkkhMvv=v,T1100,MMkkkiiiiMiifufu=F,(17)()()()T00101,Muu=u,()()()T11,kkkMgg=g,112nnn+=FFF,其中1,2,kn=。我们可以将(16)简化为 刘心愿 DOI:10.12677/aam.2024.134112 1222 应用数学进展 ()()()()1101111110000111113222,2,nnhhnnnnnnknkhhkknhnnnnnhhhgg +=+=
21、+=+=+=AuBvAAuuvuqBuAvFFw (18)所以当1n 时,矩阵形式为()10132nnhhnnhh +=AuqBvwBA (19)当0n=时,我们可以通过求解下列线性化的方程组来求解1hu和1hv()1000000000110,0,hhththh=+=+=uuABvGuGuBuAvF (20)然后我们将(20)写作()10000001100,hhhhhh=+=AAuBvGuuqBuAvFw (21)所以当0n=时,我们可以得到 1010hhhh=AuqOvwBA (22)4.数值实验数值实验 在本节中,我们将提供一个二维例子来验证算法的有效性和收敛阶。()()()()()()(
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 二维 非线性 分数 波动 方程 BDF2 WSGI 有限元 算法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。