×

含时定向空间分数平流扩散方程的快速迎风和欧拉-拉格朗日控制体积格式。 (英语) Zbl 1453.65247号

小结:我们发展了二维依赖时间的对流主导的定向空间-分形对流-扩散方程的控制体积方法,其中定向空间-分数导数通过单位圆上的概率测度在所有方向上加权,它们用于模拟地下水中溶质在地下非均质多孔介质中的各向异性超扩散运移。
我们为控制方程开发了一种快速迎风控制体积法,以消除对流项以空间为中心的数值离散中经常出现的虚假数值振荡,这些数值振荡相对容易实现。我们还为控制方程开发了一种欧拉-拉格朗日控制体积法,该方法通过将时滞项和平流项沿特征曲线组合为材料导数项来对称控制方程。这两种方法都是局部质量守恒的,这在这些应用中是必不可少的。
由于定向空分微分算子的非局部性质,相应的数值离散通常会生成全刚度矩阵。传统的直接求解器往往需要(O(N^2))内存需求,并且每个时间步长具有(O(N ^3))计算复杂性,其中(N)是空间未知数,在计算上比积分阶平流-扩散方程的数值近似要昂贵得多。在分析刚度矩阵结构的基础上,我们提出了一种快速的Krylov子空间迭代求解器来加速迎风和欧拉-拉格朗日控制体方法的数值逼近,通过直接求解器将计算复杂度从O(N^3)降低到O(Nlog N)\)每个时间步长的每个Krylov子空间迭代和从\(O(N^2)\)到\(O)(N)\)的内存需求。数值结果表明了该方法的有效性。

MSC公司:

6500万08 含偏微分方程初值和初边值问题的有限体积法
35卢比 积分-部分微分方程
26A33飞机 分数导数和积分
76S05号 多孔介质中的流动;过滤;渗流
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Arbogast,T。;Wheeler,M.F.,对流主导输运问题的特征混合有限元方法,SIAM J.Numer。分析。,32404-424(1995年)·Zbl 0823.76044号
[2] Bear,J.,《多孔材料中流体的动力学》(1972),爱思唯尔出版社:纽约·Zbl 1191.76001号
[3] Benson,D。;惠特克拉,西南部。;Meerschaert,M.M.,《勒维运动的分数阶控制方程》,《水资源》。决议,36,1413-1423(2000)
[4] Chan,R.H。;Ng,M.K.,Toeplitz系统的共轭梯度法,SIAM Rev.,38,427-482(1996)·Zbl 0863.65013号
[5] Cheng,A。;Wang,H。;Wang,K.,《反常扩散溶质运移的欧拉-拉格朗日控制体积法》,数值。方法部分差异。Equ.、。,31, 253-267 (2015) ·Zbl 1309.76136号
[6] 杜,N。;Wang,H.,(mathbb{R}^2)中有界区域上空间分数阶色散方程的快速有限元方法,SIAM J.Sci。计算。,37,A1614-A1635(2015)·Zbl 1331.65175号
[7] 欧文·V·J。;豪尔,N。;Roop,J.P.,时间相关、非线性、空间分数扩散方程的数值近似,SIAM J.Numer。分析。,45, 572-591 (2007) ·兹比尔1141.65089
[8] 欧文·V·J。;Roop,J.P.,(mathbb{R}^d)中有界区域上分数阶对流-弥散方程的变分解,数值。方法部分差异。Equ.、。,23, 256-281 (2007) ·Zbl 1117.65169号
[9] (尤因,R.E.,《油藏模拟数学》,油藏模拟数学,应用数学研究前沿,第1卷(1984年),SIAM:SIAM费城)·Zbl 0533.00031号
[10] 尤因·R·E。;罗素·T·F。;Wheeler,M.F.,用混合有限元和修正特征法对多孔介质中混溶驱替近似值的收敛性分析,计算。方法应用。机械。工程,47,73-92(1984)·Zbl 0545.76131号
[11] 尤因,R.E。;Wang,H。;大卫·斯隆;恩德雷·苏利;Vandewalle,Stefan,《含时对流偏微分方程数值方法概述》,《2000年20世纪的数值分析》,第7卷。《20世纪数值分析2000》,第7卷,《计算杂志》。申请。数学。,128, 423-445 (2001) ·Zbl 0983.65098号
[12] 郭,X。;李毅。;Wang,H.,空间分数阶扩散方程的四阶格式,计算。物理。,737, 410-424 (2018) ·Zbl 1416.65267号
[13] 贾,J。;Wang,H.,分数阶扩散方程局部精细网格上的快速有限体积法,东亚应用杂志。数学。,9, 755-779 (2019) ·Zbl 1465.65086号
[14] LeVeque,R.J.,《双曲问题的有限体积方法》,剑桥应用数学教材(2002),剑桥大学出版社:剑桥大学出版社·Zbl 1010.65040号
[15] 刘,F。;Anh,V。;特纳,I.,空间分数阶福克-普朗克方程的数值解,J.Compute。申请。数学。,166, 209-219 (2004) ·Zbl 1036.82019年
[16] Meerschaert,M.M。;Benson,D.A。;Baeumer,B.,《多维平流和分数弥散》,Phys。E版,595026-5028(1999)
[17] Meerschaert,M.M。;谢夫勒,H.P。;Tadjeran,C.,二维分数阶色散方程的有限差分方法,J.Compute。物理。,211, 249-261 (2006) ·Zbl 1085.65080号
[18] Meerschaert,M.M。;Sikorskii,A.,分数微积分的随机模型,德格鲁伊特数学研究(2011)·Zbl 1490.60004号
[19] 梅茨勒,R。;Klafter,J.,《随机行走结束时的餐厅:分数动力学描述异常运输的最新进展》,J.Phys。A、 数学。Gen.,37,R161-R208(2004)·2018年5月10日
[20] Roop,J.P.,《有界区域上分数阶对流-弥散方程有限元近似的计算方面》,(mathbb{R}^2),J.Compute。申请。数学。,193, 243-268 (2006) ·Zbl 1092.65122号
[21] Roos,H.-G。;苯乙烯,M。;Tobiska,L.,奇异摄动微分方程的稳健数值方法,计算数学中的Springer级数,第24卷(2008),Springer-Verlag:Springer-Verlag Berlin·Zbl 1155.65087号
[22] Wang,H。;Al-Lawatia,M.,瞬变平流扩散方程的局部保守欧拉-拉格朗日控制体积法,数值。方法部分差异。Equ.、。,22, 577-599 (2005) ·兹比尔1092.65076
[23] Wang,H。;Basu,T.S.,二维空间分数扩散方程的快速有限差分方法,SIAM J.Sci。计算。,34,A2444-A2458(2012)·Zbl 1256.35194号
[24] Wang,H。;香港大乐。;尤因·R·E。;Espedal,M.S。;夏普利,R.C。;Man,S.,《二维对流扩散方程的ELLAM格式》,SIAM J.Sci。计算。,20, 2160-2194 (1999) ·Zbl 0939.65109号
[25] Wang,H。;Du,N.,三维含时空间分数阶扩散方程的快速有限差分方法及其有效实现,J.Compute。物理。,253, 50-63 (2013) ·Zbl 1349.65341号
[26] Wang,H。;尤因·R·E。;秦刚。;Lyons,S.L。;Al-Lawatia,M。;Man,S.,多维对流-反应方程的欧拉-拉格朗日局部伴随方法家族,J.Compute。物理。,152, 120-163 (1999) ·Zbl 0956.76050号
[27] Wang,H。;Liang,D。;尤因·R·E。;Lyons,S.L。;Qin,G.,用欧拉-拉格朗日局部伴随法和混合有限元方法近似点源和点汇多孔介质中的混溶流体流动,SIAM J.Sci。计算。,22, 561-581 (2000) ·Zbl 0988.76054号
[28] Wang,H。;Liang,D。;尤因·R·E。;莱昂斯,S.L。;秦,G.,《点源和点汇多孔介质中可压缩流体流动的ELLAM-MFEM求解技术》,J.Compute。物理。,159, 344-376 (2000) ·Zbl 0979.76051号
[29] Wang,H。;王凯。;Sircar,T.,分数阶扩散方程的直接有限差分法,J.Compute。物理。,229, 8095-8104 (2010) ·Zbl 1198.65176号
[30] 王凯。;Wang,H.,分数阶对流扩散方程的快速特征有限差分法,Adv.Water Resour。,34, 810-816 (2011)
[31] Wu,S。;张,H。;周,T。通过对角化技术和多重网格求解时间周期分数阶扩散方程,数值。线性代数应用。,25,第2178条pp.(2018)·兹比尔1513.65347
[32] Wu,S。;Zhou,T.,分数阶扩散方程的快速准实迭代,J.Compute。物理。,329, 210-226 (2017) ·Zbl 1406.65077号
[33] Zhang,Y。;Benson,D.A。;Meerschaert,M.M。;LaBolle,E.M.,《可变参数的空间分数平流扩散方程:各种公式、数值解和MADE现场数据的应用》,水资源。研究,43(2007),W05439
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。