×

三场混合有限元耦合多孔介质力学的可缩放算法。 (英语) Zbl 1373.76312号

摘要:我们引入了一类块预条件,用于加速基于三场公式的耦合孔隙力学方程的迭代解。将位移/速度/压力混合有限元方法与一阶后向差分公式相结合,用于时间导数的近似,产生了一系列具有3乘3非对称不定块矩阵的线性系统。预处理器是通过在基于物理的自变量的帮助下近似两级Schur补码来获得的,这些自变量也可以在纯代数方法中推广。通过理论和实验分析,证明了该算法的鲁棒性、有效性和可扩展性。该性能还针对现实世界中具有挑战性的浅层地层固结实验进行了评估。

MSC公司:

76S05号 多孔介质中的流动;过滤;渗流
74层10 流固相互作用(包括气动和水弹性、孔隙度等)
65平方米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65F08个 迭代方法的前置条件
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] Hettema,M.H.H。;P.M.T.M.舒特延斯。;Verboom,B.J.M。;Gussinklo,H.J.,砂岩储层的生产诱导压实:应力路径的强烈影响,SPE Reserve。评估。工程,3,4,342-347(2000)
[2] 加里波夫,T.A。;Karimi-Fard,M。;Tchelepi,H.A.,耦合流和地质力学的离散断裂模型,计算。地质科学。,20, 1, 149-160 (2016) ·Zbl 1392.76079号
[3] 蒂蒂尼,P。;Ferronato,M。;甘布拉蒂,G。;Gonella,M.,《意大利Emilia-Romagna海岸的地下水抽水和地面沉降:模拟过去的发生和未来的趋势》,《水资源》。研究,42,1,文章W01406 pp.(2006)
[4] 罗志杰。;曾凤,基于粘塑性Biot固结理论的地面沉降和地下水开采有限元数值模拟,J.Hydrodyn。,23, 5, 615-624 (2011)
[5] 马哈茂德普尔,M。;Khamehchiyan,M。;M.R.Nikedul。;Ghassemi,M.R.,伊朗德黑兰西南平原地下水开采引起的区域地面沉降数值模拟与预测,工程地质学。,201, 6-28 (2016)
[6] 邓普西,D。;Kelkar,S。;达瓦茨,N。;希克曼,S。;Moos,D.,沙漠峰值增强地热系统剪切刺激期间注入、应力和渗透率增强的数值模拟,国际岩石力学杂志。最小科学。,78, 190-206 (2015)
[7] 潘·L。;弗赖费尔德,B。;道蒂,C。;扎肯,S。;Sheu,M。;卡特里特,B。;Terrall,T.,使用(CO_2)作为工作流体的地热开采的全耦合井筒-储层建模,地热,53,100-113(2015)
[8] Ferronato,M。;甘布拉蒂,G。;贾纳,C。;Teatini,P.,已开采气田中人为二氧化碳封存的地质力学问题,《能源转换》。管理。,51, 10, 1918-1928 (2010)
[9] 蒂蒂尼,P。;北卡罗来纳州卡斯特莱托。;Gambolia,G.,《断层地层中CO_2地质封存的三维地质力学建模:意大利亚得里亚海北部近海储层的案例研究》,国际J.Greenh。《气体控制》,22,63-76(2014)
[10] 怀特,J.A。;Chiaramonte,L。;埃泽丁,S。;Foxall,W。;Hao,Y。;A.拉米雷斯。;McNab,W.,in Salah(CO_2)储存项目储层和盖层系统的地质力学行为,Proc。国家。阿卡德。科学。,111, 24, 8747-8752 (2014)
[11] Cowin,S.C.,《骨质疏松性》,J.Biomech。,32, 3, 217-238 (1999)
[12] Frijns,A.J.H.,《应用于软骨组织的四组分混合物理论:数值建模和实验》(2000年),埃因霍温理工大学:荷兰埃因霍芬理工大学,博士论文·Zbl 0966.9202号
[13] Geertsma,J.,《石油开采工程中的岩石力学问题》,(国际岩石力学学会第一届会议,第1卷(1966年),国际岩石力学协会),585-594
[14] Biot,M.A.,《三维固结的一般理论》,J.Appl。物理。,12, 155-164 (1941) ·肯尼迪67.0837.01
[15] Coussy,O.,《多孔力学》(2004),Wiley:Wiley Chichester,英国
[16] 布雷齐,F。;Fortin,M.,《混合和混合有限元方法》(1991),Springer-Verlag:Springer-Verlag纽约,纽约,美国·Zbl 0788.7302号
[17] 刘易斯,R.W。;Schrefler,B.A.,《多孔介质静态和动态变形与固结的有限元方法》(1998),英国威利:威利奇切斯特出版社·Zbl 0935.74004号
[18] Murad,医学硕士。;Loula,A.F.D.,《关于Biot固结问题有限元近似的稳定性和收敛性》,国际J·数值。方法工程,37,4,645-667(1994)·Zbl 0791.76047号
[19] Wan,J.,地质力学与多相流耦合的稳定有限元方法(2002),斯坦福大学,博士论文
[20] 怀特,J.A。;Borja,R.I.,用于耦合固体变形/流体扩散的稳定低阶有限元及其在断层带瞬变中的应用,计算。方法应用。机械。工程,197,49-50,4353-4366(2008)·Zbl 1194.74480号
[21] 罗德里戈,C。;Gaspar,F。;胡,X。;Zikatanov,L.,Biot固结模型某些离散化的稳定性和单调性,计算。方法应用。机械。工程,298183-204(2016)·Zbl 1425.74164号
[22] Jha,B。;Juanes,R.,《流动与油藏地质力学耦合模拟的局部保守有限元框架》,《岩土学报》。,2, 3, 139-153 (2007)
[23] Phillips,P.J。;Wheeler,M.F.,多孔弹性的混合和连续Galerkin有限元方法的耦合I:时间连续情况,计算。地质科学。,131-144年11月2日(2007年)·Zbl 1117.74015号
[24] Phillips,P.J。;Wheeler,M.F.,孔隙弹性混合和连续Galerkin有限元方法的耦合II:离散时间情况,Comput。地质科学。,1145-158年11月2日(2007年)·Zbl 1117.74016号
[25] Ferronato,M。;北卡罗来纳州卡斯特莱托。;Gambolia,G.,《Biot固结的全耦合三维混合有限元模型》,J.Compute。物理。,229, 12, 4813-4830 (2010) ·Zbl 1305.76055号
[26] 哈加,J.B。;Osnes,H。;Langtangen,H.P.,《关于低渗透和低压缩多孔介质中压力振荡的原因》,国际期刊数值。分析。方法地质力学。,36, 12, 1507-1522 (2012)
[27] Turan,E。;Arbenz,P.,三维骨孔隙弹性的大规模微观有限元分析,并行计算。,40, 7, 239-250 (2014)
[28] Yi,S.-Y.,Biot固结模型新混合有限元法的收敛性分析,Numer。方法部分差异。Equ.、。,30, 4, 1189-1210 (2014) ·Zbl 1350.74024号
[29] Ferronato,M。;皮尼,G。;Gambolia,G.,《预处理在用Krylov子空间方法求解有限元耦合固结方程中的作用》,Int.J.Numer。分析。方法地质力学。,33, 3, 405-423 (2009) ·Zbl 1272.74610号
[30] 怀特,J.A。;Borja,R.,全耦合流动和地质力学的块-预处理Newton-Krylov解算器,计算。地质科学。,15, 4, 647-659 (2011) ·Zbl 1367.76034号
[31] Bergamaschi,L。;阿拉巴马州马丁内斯。,RMCP:地质力学中出现的鞍点线性系统的松弛混合约束预条件,计算。方法应用。机械。工程,221-222,54-62(2012)·Zbl 1253.74032号
[32] Kim,J。;Tchelepi,H.A。;Juanes,R.,耦合流和地质力学顺序方法的稳定性、准确性和效率,SPE J.,16,2,249-262(2011)·Zbl 1228.74106号
[33] 怀特,J.A。;北卡罗来纳州卡斯特莱托。;Tchelepi,H.A.,《耦合孔隙力学的分块求解器:统一框架》,计算。方法应用。机械。工程,303,55-74(2016)·Zbl 1425.74497号
[34] Jha,B。;Juanes,R.,《多相流和孔隙力学耦合:孔隙压力对断层滑动和地震触发影响的计算模型》,水利部。研究,53776-3808(2014)
[35] 北卡罗来纳州卡斯特莱托。;Ferronato,M。;Gambolia,G.,《利用Godunov混合方法对流体地质储量进行热-水-机械模拟》,国际期刊Numer。方法工程,90,8,988-1009(2012)·Zbl 1242.76157号
[36] 王海峰,《线性多孔弹性理论》(2000),普林斯顿大学出版社:普林斯顿大学出版,美国新泽西州普林斯顿
[37] 布克·J·R。;Small,J.C.,《关于Biot固结方程数值解稳定性的研究》,国际固体结构杂志。,907-917年7月11日(1975年)·兹比尔0311.73047
[38] 夸特罗尼,A。;Sacco,R。;Saleri,F.,《数值数学》(2007),施普林格出版社:施普林格出版社,纽约,美国·Zbl 0913.65002号
[39] O.阿克塞尔森。;布莱赫塔,R。;Byczanski,P.,多孔弹性问题的稳定离散化和产生鞍点型矩阵的有效预条件,计算。视觉。科学。,15, 4, 191-207 (2012) ·Zbl 1388.74035号
[40] Van der Vorst,H.A.,Bi-CGSTAB:非对称线性系统解的Bi-CG的一种快速平滑收敛变体,SIAM J.Sci。统计计算。,13, 2, 631-644 (1992) ·Zbl 0761.65023号
[41] Bramble,J.H。;Pasciak,J.E.,由椭圆问题的混合近似产生的不定系统的预处理技术,数学。计算。,50, 181, 1-17 (1988) ·Zbl 0643.65017号
[42] 西尔维斯特·D·J。;Elman,H.C。;Wathen,A.J.,不可压缩流线性化Navier-Stokes方程的有效预处理,J.计算。申请。数学。,128, 1-2, 261-279 (2001) ·Zbl 0983.76051号
[43] Elman,H.C。;西尔维斯特·D·J。;Wathen,A.J.,离散稳态Navier-Stokes方程鞍点预条件的性能和分析,数值。数学。,90, 4, 665-688 (2002) ·Zbl 1143.76531号
[44] Simoncini,V.,对称鞍点问题的块三角预条件,应用。数字。数学。,49, 1, 63-80 (2004) ·Zbl 1053.65033号
[45] 本兹,M。;Golub,G.H。;Liesen,J.,鞍点问题的数值解,数值学报。,14, 1-137 (2005) ·Zbl 1115.65034号
[46] 萨阿德,Y。;Schultz,M.H.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7, 3, 856-869 (1986) ·Zbl 0599.65018号
[47] 格林鲍姆,A。;普塔克,V。;Strakoš,Z.,对于GMRES,任何非增量收敛曲线都是可能的,SIAM J.Matrix Anal。申请。,17, 3, 465-469 (1996) ·Zbl 0857.65029号
[48] Bergamaschi,L。;Mantica,S。;Saleri,F.,多孔介质中达西定律的混合有限元近似(1994),技术代表CRS4-ApplMath-94-17,CRS4,意大利
[49] 赖斯,J.R。;Cleary,M.P.,《含可压缩成分的饱和流体弹性多孔介质的一些基本应力扩散解》,《地球物理学评论》。,14, 2, 227-241 (1976)
[50] Verruijt,A.,《含水层的弹性储存》(De Wiest,R.J.M.,《流经多孔介质》(1969),学术出版社:美国纽约州纽约市学术出版社),929-936
[51] 甘布拉蒂,G。;蒂蒂尼,P。;巴奥,D。;Ferronato,M.,意大利波河流域动态活动含水层中孔隙弹性耦合的重要性,水资源。研究,36,9,2443-2459(2000)
[52] Kim,J。;Tchelepi,H.A。;Juanes,R.,《耦合流和地质力学序列方法的稳定性和收敛性:固定应力和固定应变分裂》,计算。方法应用。机械。工程,200,131591-1606(2011)·兹比尔1228.74101
[53] Mikelić,A。;惠勒,M.F.,耦合流和地质力学迭代耦合的收敛性,计算。地质科学。,17, 3, 455-461 (2013) ·Zbl 1392.35235号
[54] Girault,V。;库马尔,K。;Wheeler,M.F.,断裂多孔弹性介质中地质力学与流动迭代耦合的收敛性,计算。地质科学。,20, 5, 997-1011 (2016) ·Zbl 1391.76650号
[55] Almani,T。;库马尔,K。;Dogru,A。;辛格,G。;Wheeler,M.F.,流与地质力学耦合的多速率固定应力分裂迭代格式的收敛性分析,计算。方法应用。机械。工程,311180-207(2016)·Zbl 1439.74183号
[56] 加布西,J。;Wilson,E.L.,多孔弹性介质中可压缩流体的流动,国际期刊数值。方法工程,5,3,419-442(1973)·Zbl 0248.76037号
[57] 维米尔,P.A。;Verruijt,A.,《有限元固结的精度条件》,国际期刊数字。分析。方法地质力学。,5, 1, 1-14 (1981) ·兹比尔0456.73060
[58] Ferronato,M。;甘布拉蒂,G。;Teatini,P.,有限元孔隙弹性方程的Ill条件,国际固体结构杂志。,第38、34、5995-6014页(2001年)·Zbl 1075.74643号
[59] Mandel,J.,《土壤综合》(Consolidation des sols),《岩土工程》,第3期,第7期,第287-299页(1953年)
[60] Abousleiman,Y。;Cheng,A.H.-D。;崔,L。;Detournay,E。;Roegiers,J.-C.,《重新审视曼德尔的问题》,《岩土工程》,46,2,187-195(1996)
[61] 北卡罗来纳州卡斯特莱托。;怀特,J.A。;Tchelepi,H.A.,双向耦合孔隙力学定应力迭代解的精度和收敛性,国际期刊数值。分析。方法地质力学。,39, 14, 1593-1618 (2015)
[62] HSL,用于大规模科学计算的Fortran代码集(2013)
[63] 贾纳,C。;Ferronato,M。;Gambolia,G.,对称正定线性系统的块FSAI-ILU并行预处理器,SIAM J.Sci。计算。,32, 5, 2468-2484 (2010) ·Zbl 1220.65037号
[64] 贾纳,C。;Ferronato,M。;Gambolia,G.,使用区域分解技术的增强块FSAI预处理,SIAM J.Sci。计算。,35、5、S229-S249(2013)·Zbl 1288.65034号
[65] 贾纳,C。;Ferronato,M。;萨托雷托,F。;Gambolia,G.,FSAIPACK:一个用于高性能因子稀疏近似逆预处理的软件包,ACM-Trans。数学。软质。,41, 2, 10 (2015) ·Zbl 1369.65052号
[66] 林,C。;Moré,C.C.,记忆有限的不完全Cholesky因子分解,SIAM J.Sci。计算。,21, 1, 24-45 (1999) ·Zbl 0941.65033号
[67] 甘布拉蒂,G。;皮尼,G。;Ferronato,M.,有限元固结模型中投影方法的数值性能,国际期刊数值。分析。方法地质力学。,25, 14, 1429-1447 (2001) ·Zbl 1112.74513号
[68] O.阿克塞尔森。;Gustafsson,I.,求解Navier弹性方程的迭代方法,计算。方法应用。机械。工程师,15,2,241-258(1978)·Zbl 0402.73028号
[69] Blaheta,R.,线性弹性问题的位移分解-完全因式分解预处理技术,数值。线性代数应用。,1, 2, 107-128 (1994) ·Zbl 0837.65021号
[70] 古斯塔夫森,I。;Lindskog,G.,《线性弹性问题的并行求解》,第一部分:理论,数值。线性代数应用。,5, 2, 123-139 (1998) ·Zbl 0935.74079号
[71] Gee,M.W。;西弗特,C.M。;胡建杰。;杜米纳罗,R.S。;Sala,M.G.,ML 5.0平滑聚合用户指南(2006),桑迪亚国家实验室,技术代表SAND2006-2649
[72] 萨拉,M。;Heroux,M.,《使用IFPACK 3.0的稳健代数预条件》(2005),Sandia国家实验室,技术代表SAND2005-0662
[73] 班格思,W。;哈特曼,R。;Kanschat,G.,Deal,II一个通用的面向对象有限元库,ACM Trans。数学。软质。,33、4、24(2007年)·Zbl 1365.65248号
[74] Heroux,医学硕士。;Bartlett,R.A。;豪尔,V.E。;霍克斯特拉,R.J。;胡建杰。;科尔达,T.G。;Lehoucq,R.B。;Long,K.R。;Pawlowski,R.P。;Phipps,E.T。;塞林格,A.G。;香港桑奎斯特。;杜米纳罗,R.S。;Willenbring,J.M。;威廉姆斯。;Stanley,K.S.,Trilinos项目概述,ACM Trans。数学。软质。,31, 3, 397-423 (2005) ·兹比尔1136.65354
[75] 贝伦戈,V。;Leoni,M。;Simonini,P.,威尼斯泻湖淤泥随时间变化行为的数值模拟,(Singh,D.,《国际计算机方法与地质力学进展协会第十二届国际会议论文集》,国际计算机方法和地质力学进展学会第十二届会议论文集,IACMAG(2008)),929-936
[76] 北卡罗来纳州卡斯特莱托。;甘布拉蒂,G。;Teatini,P.,威尼斯海岸大规模荷载试验的耦合MFE孔隙力学模型,计算。地质科学。,19, 1, 17-29 (2015)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。