×

耦合孔隙力学的增强松弛物理因子分解预处理。 (英语) Zbl 1524.74119号

摘要:松弛物理因子分解(RPF)预条件是一种新的算法,它可以有效且稳健地求解块线性系统,该系统由耦合孔隙力学的三场位移-速度-压力公式产生。然而,为了应用它,有必要反转代数形式的块,其中(C)是对称正定矩阵,(FF^T)是秩亏项,(β)是实非负系数。以不精确的方式执行的\(\hat{C}\)的反演对于\(\beta\)的大值可能会变得不稳定,因为它通常发生在全孔隙力学模拟的某些阶段。在这项工作中,我们提出了一系列代数技术来稳定不精确解。这种策略在可能出现此类问题的其他问题中也很有用,例如Navier-Stokes或不可压缩弹性的增强拉格朗日预处理技术。首先,我们引入了一个通过矩阵{C}的自然分裂得到的迭代格式。其次,我们开发了一种基于使用适当的投影算子来消除({C})的近核模式的技术。这两种方法都产生了一类新的预处理器,称为增强型RPF(ERPF)。所提算法的有效性和鲁棒性在理论基准测试和实际大型应用中都得到了验证,其性能优于本地RPF预条件器。

MSC公司:

74层10 流固相互作用(包括气动和水弹性、孔隙度等)
76S05号 多孔介质中的流动;过滤;渗流
76M10个 有限元方法在流体力学问题中的应用
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65层10 线性系统的迭代数值方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Biot,M.A.,《三维固结的一般理论》,J.Appl。物理。,12, 155-164 (1941)
[2] Coussy,O.,Poromechanics(2004),威利:威利奇切斯特,英国
[3] Frijns,A.J.H.,应用于软骨组织的四组分混合物理论:数值建模和实验(2000年),埃因霍温理工大学:荷兰埃因霍芬理工大学,博士论文·Zbl 0966.9202号
[4] Phillips,P.J。;Wheeler,M.F.,多孔弹性的混合和连续Galerkin有限元方法的耦合I:时间连续情况,计算。地质科学。,11, 2, 131-144 (2007) ·Zbl 1117.74015号
[5] Phillips,P.J。;Wheeler,M.F.,多孔弹性的混合和连续Galerkin有限元方法的耦合II:离散时间情况,计算。地质科学。,11, 2, 145-158 (2007) ·Zbl 1117.74016号
[6] Bause,M。;Radu,F.A。;Köcher,U.,具有迭代耦合的Biot孔隙弹性系统的时空有限元近似,计算。方法应用。机械。工程,320,745-768(2017)·Zbl 1439.74389号
[7] 胡,X。;罗德里戈,C。;Gaspar,F.J。;Zikatanov,L.T.,多孔弹性中Biot固结模型的非协调有限元方法,J.Compute。申请。数学。,310, 143-154 (2017) ·Zbl 1381.76175号
[8] 两者,J.W。;库马尔,K。;Nordbotten,J.M。;Radu,F.A.,Anderson非饱和多孔介质固结的加速固定应力劈裂方案,计算。数学。申请。,77, 6, 1479-1502 (2019) ·Zbl 1442.65251号
[9] Girault,V。;库马尔,K。;Wheeler,M.F.,断裂多孔弹性介质中地质力学与流动迭代耦合的收敛性,计算。地质科学。,20, 5, 997-1011 (2016) ·Zbl 1391.76650号
[10] Girault,V。;惠勒,M.F。;Almani,T。;Dana,S.,通过固定应力算法求解的离散多孔弹性系统的先验误差估计,石油天然气科学。Technol.-技术修订版IFP新能源,74,24(2019)
[11] Ferronato,M。;北卡罗来纳州卡斯特莱托。;Gambolia,G.,《Biot固结的全耦合三维混合有限元模型》,J.Compute。物理。,229, 12, 4813-4830 (2010) ·Zbl 1305.76055号
[12] 北卡罗来纳州卡斯特莱托。;怀特,J.A。;Ferronato,M.,三场混合有限元耦合孔隙力学的可缩放算法,J.Compute。物理。,327894-918(2016)·Zbl 1373.76312号
[13] 刘易斯,R.W。;Schrefler,B.A.,《多孔介质静态和动态变形与固结的有限元方法》(1998),英国威利:威利奇切斯特出版社·Zbl 0935.74004号
[14] 北卡罗来纳州卡斯特莱托。;Ferronato,M。;Gamboliti,G.,通过Godunov混合方法对流体地质储量进行的热力-水力-机械建模,Int.J.Numer。方法工程,90,8,988-1009(2012)·Zbl 1242.76157号
[15] Jha,B。;Juanes,R.,《多相流和孔隙力学耦合:孔隙压力对断层滑动和地震触发影响的计算模型》,水利部。研究,53776-3808(2014)
[16] Choo,J。;Lee,S.,局部质量守恒耦合多孔介质力学的强化Galerkin有限元,计算。方法应用。机械。工程,341,311-332(2018)·Zbl 1440.74120号
[17] Choo,J.,《具有局部质量守恒的大变形孔隙力学:一个丰富的Galerkin有限元框架》,国际期刊Numer。方法工程,116,1,66-90(2018)
[18] Lipnikov,K.,多孔弹性中Biot模型的数值方法(2002),休斯顿大学,博士论文
[19] 罗德里戈,C。;胡,X。;欧姆,P。;阿德勒,J.H。;Gaspar,F.J。;Zikatanov,L.T.,多孔弹性和Stokes方程的新稳定离散化,计算。方法应用。机械。工程,341467-484(2018)·Zbl 1440.76027号
[20] 牛,C。;鲁伊·H。;Hu,X.,多孔弹性的稳定杂交混合有限元方法,计算。地质科学。,757-774(2021)·Zbl 1460.65125号
[21] Honorio,H。;Maliska,C。;Ferronato,M。;Janna,C.,多孔弹性问题的基于稳定单元的有限体积法,J.Compute。物理。,364, 49-72 (2018) ·Zbl 1392.74090号
[22] 卡马戈,J.T。;怀特,J.A。;Borja,R.I.,《多相孔隙力学的宏观元素稳定性》,计算。地质科学。,25, 775-792 (2021) ·Zbl 1460.65118号
[23] 弗里戈,M。;北卡罗来纳州卡斯特莱托。;Ferronato,M。;White,J.A.,混合三场混合有限元耦合孔隙力学的高效求解器,计算。数学。申请。,91, 36-52 (2021) ·Zbl 1524.76451号
[24] Van der Vorst,H.A.,Bi-CGSTAB:用于非对称线性系统解的Bi-CG的快速平滑收敛变体,SIAM J.Sci。统计计算。,13, 2, 631-644 (1992) ·Zbl 0761.65023号
[25] 萨阿德,Y。;Schultz,M.H.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7, 3, 856-869 (1986) ·Zbl 0599.65018号
[26] Jha,B。;Juanes,R.,《流动与油藏地质力学耦合模拟的局部保守有限元框架》,《岩土学报》。,2, 3, 139-153 (2007)
[27] Almani,T。;库马尔,K。;Dogru,A。;辛格,G。;Wheeler,M.,流与地质力学耦合的多速率固定应力分裂迭代格式的收敛性分析,计算。方法应用。机械。工程,311180-207(2016)·Zbl 1439.74183号
[28] 两者,J.W。;博雷加莱斯,M。;Nordbotten,J.M。;库马尔,K。;Radu,F.A.,非均匀介质中Biot方程的稳健固定应力分裂,应用。数学。莱特。,68, 101-108 (2017) ·Zbl 1383.74025号
[29] 博雷加莱斯,M。;Radu,F.A。;库马尔,K。;Nordbotten,J.M.,非线性多孔力学的稳健迭代方案,计算。地质科学。,22, 4, 1021-1038 (2018) ·Zbl 1402.65109号
[30] 达纳,S。;Ganis,B。;Wheeler,M.F.,深层地下储层中耦合流动和孔隙力学的多尺度固定应力分裂迭代方案,J.Compute。物理。,352,1-22(2018)·Zbl 1375.76085号
[31] 达纳,S。;惠勒,M.F.,非均质多孔弹性介质中流动与变形耦合的两网格固定应力分裂迭代格式的收敛性分析,计算。方法应用。机械。工程师,341788-806(2018)·Zbl 1440.74122号
[32] 洪,Q。;克劳斯,J。;林伯里,M。;Wheeler,M.F.,多渗透孔隙弹性系统固定应力分裂迭代法的参数稳健收敛分析,多尺度模型。模拟。,18, 2, 916-941 (2020) ·Zbl 1447.65077号
[33] 库兹涅佐夫,Y。;Lipnikov,K。;里昂,S。;Maliassov,S.,《多孔弹性问题的数学建模和数值算法》,(Chen,Z.;Glowinski,R.;Li,K.,《科学计算的当前趋势》,《当代数学》,第329卷(2003),Amer。数学系:Amer。美国罗得岛普罗维登斯数学学会),191-202年·Zbl 1138.74393号
[34] 怀特,J.A。;Borja,R.,全耦合流动和地质力学的块-预处理Newton-Krylov解算器,计算。地质科学。,15464-659(2011年)·Zbl 1367.76034号
[35] 哈加,J.B。;Osnes,H。;Langtangen,H.P.,具有高度非均匀材料参数的大规模多孔弹性的并行块预处理程序,计算。地质科学。,16, 3, 723-734 (2012)
[36] 哈加,J.B。;Osnes,H。;Langtangen,H.P.,《关于低渗透和低压缩多孔介质中压力振荡的原因》,国际期刊数值。分析。方法地质力学。,36, 12, 1507-1522 (2012)
[37] O.阿克塞尔森。;布莱赫塔,R。;Byczanski,P.,多孔弹性问题的稳定离散化和产生鞍点型矩阵的有效预条件,计算。视觉。科学。,15, 4, 191-207 (2012) ·Zbl 1388.74035号
[38] Turan,E。;Arbenz,P.,三维骨孔隙弹性的大规模微观有限元分析,并行计算。,40, 7, 239-250 (2014)
[39] Gaspar,F.J。;Rodrigo,C.,关于流动与地质力学耦合的多重网格方法中更平滑的固定应力分裂方案,计算。方法应用。机械。工程,326526-540(2017)·兹比尔1439.74413
[40] 罗,P。;罗德里戈,C。;Gaspar,F.J。;Oosterlee,C.W.,关于多孔弹性方程多重网格中的Uzawa平滑器,数值。线性代数应用。,24,1,第2074条pp.(2017)·Zbl 1424.76041号
[41] Lee,J.J。;马尔达尔,K.-A。;Winther,R.,Biot固结模型的参数稳健离散化和预处理,SIAM J.Sci。计算。,39,1,A1-A24(2010)·Zbl 1381.76183号
[42] 阿德勒,F.J。;Gaspar,J.H。;胡,X。;罗德里戈,C。;Zikatanov,L.T.,《Biot模型的鲁棒块预处理器》,(Bjöstad,P.E.;Brenner,s.;Halpern,L.;Kornhuber,R.;Kim,H.H.;Rahman,T.;Widlund,O.B.,《科学与工程领域分解方法》第二十四卷(2018年),施普林格国际出版公司),3-16·Zbl 1442.65342号
[43] 洪,Q。;Kraus,J.,Biot固结模型经典三场公式的参数鲁棒稳定性,电子。变速器。数字。分析。,48, 202-226 (2018) ·Zbl 1398.65046号
[44] Ferronato,M。;弗朗西斯科尼,A。;詹纳,C。;北卡罗来纳州卡斯特莱托。;Tchelepi,H.,耦合多物理问题的一般预处理框架,J.Compute。物理。,398,第108887条pp.(2019)·Zbl 1453.65065号
[45] 阿德勒,J.H。;Gaspar,F.J。;胡,X。;欧姆,P。;罗德里戈,C。;Zikatanov,L.T.,多孔弹性方程新稳定离散化的鲁棒预条件,SIAM J.Sci。计算。,42、3、B761-B791(2020年)·Zbl 1448.65145号
[46] Bui,Q.M。;Osei-Kuffuor,D。;北卡罗来纳州卡斯特莱托。;White,J.A.,《非均质介质多相孔隙力学的可扩展多重网格简化框架》,SIAM J.Sci。计算。,42、2、B379-B396(2020)·Zbl 1435.65153号
[47] 弗里戈,M。;北卡罗来纳州卡斯特莱托。;Ferronato,M.,混合有限元耦合孔隙力学的松弛物理因子分解预处理,SIAM J.Sci。计算。,41,B694-B720(2019年)·兹比尔1420.65032
[48] Benzi,M。;Ng,M。;牛,Q。;Wang,Z.,不可压缩Navier-Stokes方程的松弛维因式分解预条件,J.Compute。物理。,230, 16, 6185-6202 (2011) ·Zbl 1419.76433号
[49] Benzi,M。;德帕里斯,S。;格兰佩林,G。;Quarteroni,A.,松弛维因式分解预条件的参数估计及其在血流动力学中的应用,计算。方法应用。机械。工程师,300,补遗C,129-145(2016)·Zbl 1423.76212号
[50] 布林克,美国。;Stein,E.,关于不可压缩和几乎不可压缩有限弹性的一些混合有限元方法,计算。机械。,19, 105-119 (1996) ·Zbl 0889.73066号
[51] Olshanskii,M。;Benzi,M.,《水动力稳定性线性化问题的增广拉格朗日方法》,SIAM J.Sci。计算。,30, 1459-1473 (2007) ·Zbl 1162.76031号
[52] Benzi,M。;Olshanskii,M.,线性化Navier-Stokes问题增广拉格朗日预条件的值域收敛性分析,SIAM J.Numer。分析。,49770-788(2011年)·Zbl 1245.76044号
[53] Benzi,M。;Olshanskii,M。;Wang,Z.,《不可压缩Navier-Stokes方程的修正增广拉格朗日预条件》,国际期刊数值。《液体方法》,66,486-508(2011)·Zbl 1421.76152号
[54] 法雷尔,体育。;Mitchell,L.公司。;Wechsung,F.,高雷诺数下三维定常不可压缩Navier-Stokes方程的增广拉格朗日预处理,SIAM J.Sci。计算。,41、5、A3073-A3096(2019)·Zbl 1448.65261号
[55] 北卡罗来纳州卡斯特莱托。;怀特,J。;Tchelepi,H.,双向耦合孔隙力学定应力迭代解的精度和收敛性,国际J。数值。分析。方法地质力学。,39, 1593-1618 (2015)
[56] 怀特,J.A。;北卡罗来纳州卡斯特莱托。;Tchelepi,H.A.,《耦合孔隙力学的分块求解器:统一框架》,计算。方法应用。机械。工程,303,55-74(2016)·Zbl 1425.74497号
[57] Bergamaschi,L。;Mantica,S。;Manzini,G.,黑油模型的有限元-有限元混合体积公式,SIAM J.Sci。计算。,20, 970-997 (1998) ·Zbl 0959.76039号
[58] Hackbusch,W.,大型稀疏方程组的迭代解(2016),施普林格:瑞士施普林格商学院·Zbl 1347.65063号
[59] Mandel,J.,《土壤综合》(Consolidation des sols),《岩土工程》,第3期,第7期,第287-299页(1953年)
[60] Saad,Y.,ILUT:双阈值不完全ILU因式分解,Numer。线性代数应用。,1, 387-402 (1994) ·Zbl 0838.65026号
[61] 林,C。;Moré,J.,记忆有限的不完全Cholesky因子分解,SIAM J.Sci。计算。,21, 1, 24-45 (1999) ·兹伯利0941.65033
[62] Ruge,J.Q。;Stüben,K.,代数多重网格,应用数学前沿,第3卷,73-130(1987),工业和应用数学学会:工业和应用算术学会,宾夕法尼亚州费城,美国,第4章·Zbl 0659.65094号
[63] Boyle,J。;Mihajlovic,医学博士。;Scott,J.A.,HSL_MI20:三维有限元问题的高效AMG预处理程序,国际期刊Numer。方法工程,82,64-98(2010)·Zbl 1183.76799号
[64] Ferronato,M。;Gazzola,L。;北卡罗来纳州卡斯特莱托。;蒂蒂尼,P。;Zhu,L.,北京地区地面沉降预测的耦合混合有限元Biot模型,(Vandamme,M.;Dangla,P.;Pereira,J.;Ghabezloo,S.,Poromechanics VI(2017),美国土木工程师学会),182-189
[65] 朱,L。;戴,Z。;龚,H。;盖博,C。;Teatini,P.,冲积扇含水层特征多区域过渡概率模型的统计反演,Stoch。环境。Res.风险评估。,30, 1005-1016 (2016)
[66] 朱,L。;弗朗西斯科尼,A。;龚,H。;Ferronato,M。;戴,Z。;Ke,Y。;潘,Y。;李,X。;王,R。;Teatini,P.,《朝拜平原地面沉降的三维相和地质力学模拟》,北京,水利。决议,56,3,条款e2019WR027026 pp.(2020)
[67] 克里斯蒂,M。;Blunt,M.,第十个SPE比较解决方案项目:升级技术的比较,SPE Reserve。评估。工程,4,4,308-316(2001)
[68] 陈,Z。;Huan,G.等人。;Ma,Y.,《多孔介质中多相流的计算方法》(2006),工业和应用数学学会:美国宾夕法尼亚州费城工业与应用数学学会·Zbl 1092.76001号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。