×

多孔介质非线性流动问题的非协调时间网格迭代方法。 (英语) Zbl 1516.65082号

小结:多孔介质中的部分饱和流动通常由Richards方程模拟,该方程是非线性抛物线方程,并且可能退化。本文提出了基于区域分解的Richards方程数值格式,其中不同的时间步长可以用于不同的子域。在混合公式中导出了两种时域全局分解方法:第一种方法基于物理传输条件,第二种方法基于等效Robin传输条件。对于每种方法,我们使用子结构技术将原始问题重写为定义在子域之间的时空接口上的非线性问题。这种时空界面问题用牛顿法线性化,然后用GMRES迭代求解;每个GMRES迭代都涉及子域中时间相关问题的并行求解。通过二维数值实验,验证并比较了所提方法与局部时间步长方法的收敛性和准确性。

MSC公司:

65M55型 多重网格方法;涉及偏微分方程的初值和初边值问题的区域分解
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
6500万06 含偏微分方程初值和初边值问题的有限差分方法
65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65时10分 方程组解的数值计算
65层10 线性系统的迭代数值方法
2005年5月 并行数值计算
65M50型 涉及偏微分方程初值和初边值问题数值解的网格生成、细化和自适应方法
76S05号 多孔介质中的流动;过滤;渗流
35K55型 非线性抛物方程
35问题35 与流体力学相关的PDE

软件:

自由Fem++
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 艾哈迈德·E。;Ali Hassan,S。;贾普特,C。;科恩,M。;Vohralík,M.,不同岩石类型之间两相流时空域分解的后验误差估计和停止标准,SMAI J.Compute。数学。,5, 195-227 (2019) ·Zbl 1432.65130号 ·doi:10.5802/smai-jcm.47
[2] Alt,硬件;Luckhaus,S.,拟线性椭圆-抛物微分方程,数学。Z.,183,3,311-341(1983)·Zbl 0497.35049号 ·doi:10.1007/BF01176474
[3] Arbogast,T。;Wheeler,MF,多孔介质流动中退化抛物方程的非线性混合有限元法,SIAM J.Numer。分析。,33, 4, 1669-1687 (1996) ·Zbl 0856.76033号 ·doi:10.1137/S0036142994266728
[4] 贝尔·J。;Bachmat,Y.,《多孔介质中传输现象建模导论》(1991年),诺威尔:克鲁沃,诺威尔·Zbl 0780.76002号
[5] 本杰明·D·。;MJ甘德;Gouarin,L。;Halpern,L.,二维对流反应扩散方程的优化Schwarz波形松弛,数值。数学。,134, 3, 513-567 (2016) ·Zbl 1356.65225号 ·doi:10.1007/s00211-015-0784-8
[6] 本杰明·D·。;MJ甘德;Halpern,L.,应用于优化Schwarz波形弛豫的同形最佳逼近问题,数学。公司。,78, 265, 185-223 (2009) ·Zbl 1198.65179号 ·doi:10.1090/S0025-5718-08-02145-5
[7] Bergamaschi,L。;Putti,M.,Richards方程解的混合有限元和牛顿型线性化,Int.J.Numer。方法。工程师。,45, 1025-1046 (1999) ·Zbl 0943.76047号 ·doi:10.1002/(SICI)1097-0207(19990720)45:8<1025::AID-NME615>3.0.CO;2克
[8] Berninger,H.、Kornhuber,R.、Sander,O.:关于跳跃非线性的非线性Dirichlet-Neumann算法。收录于:Widlund,O.B.,Keyes,D.E.(eds.)《科学与工程领域分解方法》第十六卷,LNCSE第55卷,第483-490页。施普林格,柏林(2007)
[9] Berninger,H.,Kornhuber,R.,Sander,O.:非线性传输问题的Dirichlet Neumann和Robin方法的收敛性。收录于:Huang,Y.,Kornhuber,R.,Widlund,O.,Xu,J.(编辑)《科学与工程领域分解方法十九》,LNCSE第78卷。柏林施普林格出版社(2011)·Zbl 1217.65200号
[10] Berninger,H。;Sander,O.,Signorini型问题的子结构和非均质土壤中Richards方程的Robin方法,计算。视觉。科学。,13, 5, 187-205 (2010) ·Zbl 1337.74035号 ·doi:10.1007/s00791-010-0141-5
[11] Berninger,H。;Kornhuber,R。;Sander,O.,分层土壤中Richards方程的多域离散化,计算。地质科学。,19, 1, 213-232 (2015) ·Zbl 1326.76055号 ·doi:10.1007/s10596-014-9461-8
[12] 布莱奥,E。;德布鲁,L。;Lemarié,F.,关于具有不连续和空间可变系数的扩散方程的优化全局实时Schwarz算法。第1部分:恒系数情况,电子。事务处理。数字。分析。,40, 170-186 (2013) ·Zbl 1288.65138号
[13] Blayo,E.,Halpern,L.,Japhet,C.:不连续系数耦合对流扩散问题的优化Schwarz波形松弛算法,具有非协调时间离散化。收录于:科学与工程领域分解方法十六。莱克特。注释计算。科学。《工程》,第55卷,第267-274页。施普林格,柏林(2007)
[14] Brenner,K。;Cancès,C.,《通过参数化提高牛顿法的性能:理查兹方程的情况》,SIAM J.Numer。分析。,55, 4, 1760-1785 (2017) ·Zbl 1369.65106号 ·doi:10.137/16M1083414
[15] 马萨诸塞州西莉亚;东部布卢塔斯;Zarba,RL,非饱和流动方程的一般质量守恒数值解,水资源。Res.,26,7,1483-1496(1990)·doi:10.1029/WR026i007p01483文件
[16] Dolean,V。;MJ甘德;Kheriji,W。;郭,F。;Masson,R.,《非线性预处理:如何使用非线性Schwarz方法预处理牛顿法》,SIAM J.Sci。计算。,38,A3357-3380(2016)·Zbl 1352.65326号 ·数字对象标识码:10.1137/15M102887X
[17] Dryja,M。;Hackbusch,W.,关于非线性区域分解方法,BIT,37,2,296-311(1997)·Zbl 0891.65126号 ·doi:10.1007/BF02510214
[18] van Duijn,CJ;Pelecier,LA,部分饱和多孔介质中的非稳态过滤,Arch。理性力学。分析。,78, 173-198 (1982) ·Zbl 0502.76101号 ·doi:10.1007/BF00250838
[19] Eymard,R。;Hilhorst,D。;Vohralík,M.,退化抛物问题的组合有限体积非一致/混合有限元格式,Numer。数学。,105, 1, 73-131 (2006) ·Zbl 1108.65099号 ·数字对象标识代码:10.1007/s00211-006-0036-z
[20] Eymard,R。;古特尼克,M。;Hillhorst,D.,理查兹方程的有限体积法,计算。地质科学。,3, 259-294 (1999) ·Zbl 0953.76060号 ·doi:10.1023/A:1011547513583
[21] MJ甘德;Halpern,L。;Nataf,F.,一维波动方程的最佳Schwarz波形松弛,SIAM J.Numer。分析。,41, 1643-1681 (2003) ·Zbl 1085.65077号 ·doi:10.137/S003614290139559X
[22] MJ甘德;Halpern,L.,对流反应扩散问题的优化schwarz波形松弛,SIAM J.Numer。分析。,45, 2, 666-697 (2007) ·Zbl 1140.65063号 ·doi:10.1137/050642137
[23] Gander,M.J.,Halpern,L.,Kern,M.:用于不连续系数和非匹配网格的对流-扩散-反应问题的Schwarz波形松弛方法。收录于:科学与工程领域分解方法十六,Lect。注释计算。科学。《工程》,第55卷,第283-290页。施普林格,柏林(2007)
[24] MJ甘德;郭,F。;Mandal,BC,Dirichlet Neumann和Neumann Neumann抛物型问题的波形松弛算法,电子。事务处理。数字。分析。,45, 424-456 (2016) ·Zbl 1355.65128号
[25] MJ甘德;郭,F。;Mandal,BC,Dirichlet-Neumann多子域抛物线和双曲线问题的波形松弛方法,BIT-Numer。数学。,61, 173-207 (2021) ·Zbl 1466.65111号 ·doi:10.1007/s10543-020-00823-2
[26] Gander,M.J.,Lunowa,S.B.,Rohde,C.:非线性平流扩散方程的非重叠Schwarz波形松弛。U资产计算数学预印本编号UP-21-03(2021)·Zbl 1516.65081号
[27] 格洛温斯基,R。;丁,QV;Periaux,J.,流体动力学非线性问题的区域分解方法,计算。方法应用。机械。工程,40,1,27-109(1983)·Zbl 0505.76068号 ·doi:10.1016/0045-7825(83)90045-2
[28] Haeberlein,F.,Halpern,L.,Michel,A.:Newton-Schwarz优化了非线性反应输运的波形松弛Krylov加速器。收录于:科学与工程领域分解方法XX。莱克特。注释计算。科学。《工程》,第91卷,第387-394页。斯普林格,海德堡(2013)
[29] Halpern,L.,Japhet,C.,Omnes,P.:多孔介质应用的非一致时域分解方法。In:Pereira,J.C.F.,Sequeira,A.(eds.)《第五届欧洲计算流体动力学会议论文集》,葡萄牙里斯本(2010)
[30] Halpern,L。;贾普特,C。;Szeftel,J.,针对异质性问题的优化Schwarz波形松弛和间断Galerkin时间步进,SIAM J.Numer。分析。,50, 5, 2588-2611 (2012) ·Zbl 1262.65127号 ·doi:10.137/120865033
[31] Hecht,F.,《FreeFem++的新发展》,J.Numer。数学。,20, 251-265 (2012) ·Zbl 1266.68090号 ·doi:10.1515/jnum-2012-0013
[32] Helmig,R.,《地下多相流和输运过程:对水文系统建模的贡献》(1997),柏林:斯普林格出版社,柏林·doi:10.1007/978-3-642-60763-9
[33] Hoang,TTP,混合公式中平流-扩散问题的完全隐式局部时间步长方法,计算。数学。申请。,118248-264(2022)·Zbl 1524.65925号
[34] Hoang,TTP;贾夫雷,J。;贾普特,C。;科恩,M。;Roberts,JE,混合公式中扩散问题的时空域分解方法,SIAM J.Numer。分析。,51, 6, 3532-3559 (2013) ·Zbl 1295.65095号 ·数字对象标识代码:10.1137/130914401
[35] Hoang,TTP;贾普特,C。;科恩,M。;Roberts,JE,混合公式中平流-扩散问题的时空域分解,数学。计算。仿真。,137, 366-389 (2017) ·Zbl 07313831号 ·doi:10.1016/j.matcom.2016.11.002
[36] Hoang,TTP;Lee,H.,耦合非线性Stokes和Darcy流的全局时域分解方法,J.Sci。计算。,87, 1, 1-22 (2021) ·Zbl 1466.65112号 ·doi:10.1007/s10915-021-01422-1
[37] Jäger,W。;Kačur,J.,用松弛格式求解双重非线性和退化抛物问题,RAIRO模型。数学。分析。数字。,29, 605-627 (1995) ·Zbl 0837.65103号 ·doi:10.1051/m2安/1995290506051
[38] 克拉森,RA;半径,FA;Eigestad,GT,MPFA在三角剖分和Richards方程中的收敛性,国际。J.数字。《液体方法》,58,12,1327-1351(2008)·Zbl 1391.76341号 ·doi:10.1002/fld.1787
[39] 列表,F。;Radu,FA,求解理查兹方程的迭代方法研究,计算。地质科学。,20, 2, 341-353 (2016) ·Zbl 1396.65143号 ·doi:10.1007/s10596-016-9566-3
[40] 洛特,P。;Walker,H。;伍德沃德,C。;Yang,U.,与可变饱和流量相关的非线性系统的加速Picard方法,Adv.Water Resour。,38, 92-101 (2012) ·doi:10.1016/j.advwatres.2011.12.013
[41] Martin,V.,二维非定常对流扩散方程的优化Schwarz波形松弛方法,Appl。数字。数学。,52, 401-428 (2005) ·Zbl 1070.65088号 ·doi:10.1016/j.apnum.2004.08.022
[42] Mitra,K。;Pop,IS,解决非线性扩散问题的改进L格式,计算。数学。申请。,77, 1722-1738 (2019) ·Zbl 1442.65215号
[43] Mitra,K.,Vohralik,M.:理查兹方程的后验误差估计。arXiv:2108.12507(2021)
[44] Otto,F.,拟线性椭圆-抛物方程的L1-压缩和唯一性,《微分方程》,131,1,20-38(1996)·Zbl 0862.35078号 ·doi:10.1006/jdeq.1996.0155
[45] Pop,IS,Richards方程时间离散化方法的误差估计,计算。地质科学。,6, 141-160 (2002) ·Zbl 1079.65542号 ·doi:10.1023/A:1019936917350
[46] 波普,IS;半径,FA;Knabner,P.,《Richards方程的混合有限元:线性化程序》,J.Compute。申请。数学。,168, 365-373 (2004) ·Zbl 1057.76034号 ·doi:10.1016/j.cam.2003.04.008
[47] 波普,IS;Schweizer,B.,退化理查兹方程和流出条件的正则化方案,数学。模型方法应用。科学。(M3AS),21,1685-1712(2011)·Zbl 1231.35104号 ·doi:10.1142/S0218202511005532
[48] 半径,FA;波普,IS;Knabner,P.,一些退化抛物方程的混合有限元离散化的误差估计,Numer。数学。,109, 285-311 (2008) ·兹比尔1141.65071 ·doi:10.1007/s00211-008-0139-9
[49] 理查兹,洛杉矶,液体通过多孔介质的毛细传导,物理学,1,5,318-333(1931)·Zbl 0003.28403号 ·数字对象标识代码:10.1063/1.1745010
[50] Richardson,LF,《数值过程天气预报》(1922),剑桥:剑桥大学出版社,剑桥
[51] 施奈德,E。;Knabner,P。;Radu,FA,Richards方程混合有限元离散化的先验误差估计,Numer。数学。,98, 353-370 (2004) ·Zbl 1075.76042号 ·doi:10.1007/s00211-003-0509-2
[52] Schweizer,B.,具有干燥区域的非饱和多孔介质中流出问题的正则化,J.Differ。等于。,237, 2, 278-306 (2007) ·Zbl 1137.35043号 ·doi:10.1016/j.jde.2007.03.011
[53] 塞乌斯,D。;Mitra,K。;波普,IS;半径,FA;Rohde,C.,多孔介质中部分饱和流动的线性区域分解方法,计算。方法应用。机械。工程,33331-355(2018)·Zbl 1440.76149号 ·doi:10.1016/j.cma.2018.01.029
[54] JO斯科格斯塔德;Keilegavlen,E。;Nordbotten,JM,多孔介质非线性流动问题的区域分解策略,J.Compute。物理。,234, 439-451 (2013) ·doi:10.1016/j.jcp.2012.10.001
[55] 泰语,X-C;Espedal,M.,线性和非线性问题的一些空间分解方法的收敛速度,SIAM J.Numer。分析。,35, 4, 1558-1570 (1998) ·Zbl 0915.65063号 ·doi:10.137/S0036142996297461
[56] 伍德沃德,CS;Dawson,CN,《非线性抛物方程模拟变饱和多孔介质流动的扩展混合有限元方法分析》,SIAM J.Numer。分析。,37, 3, 701-724 (2000) ·Zbl 0948.65096号 ·doi:10.137/S0036142996311040
[57] Yotov,I.,多孔介质流动中退化抛物方程在非匹配多块网格上的混合有限元离散,东西方J.Numer。数学。,5, 211-230 (1997) ·Zbl 0897.76057号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。