安德烈亚·博里奥；弗朗索瓦·哈蒙。；尼古拉·卡斯泰莱托；约书亚·A·怀特。；兰多夫·塞特卡斯特（Randolph R.Settgast）。 耦合孔隙力学的混合模拟有限差分和虚拟元公式。 （英语） Zbl 1506.74390号 计算。方法应用。机械。工程师。 383，文章ID 113917，26 p.（2021）. 摘要：我们提出了一种非结构网格上单相耦合多孔介质力学的混合模拟有限差分和虚拟元公式。该方案的主要优点是，它在包含具有任意形状的高度畸变单元的复杂网格上收敛。我们使用基于非结构化宏观元素的局部压力跳跃稳定方法，以防止接近不排水条件的不可压缩问题中出现虚假压力模式。使用一个专门为鞍点系统设计的块三角预条件器，可以获得一个可扩展的线性解策略。我们的方法在二维基准问题上的精度和效率得到了数值验证。 引用于6文件 MSC公司： 74S05号 有限元方法在固体力学问题中的应用 65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法 74层10 流固相互作用（包括空气弹性和水弹性、孔隙率等） 关键词：多孔弹性；模拟有限差分；虚元法；任意多边形网格 软件：Gms小时；HSL_MI20型；Matlab公司；MRST公司；PolyMesher公司 PDF格式BibTeX公司 XML格式引用 \textit{A.Borio}等人，《计算》。方法应用。机械。工程383，文章ID 113917，26 p.（2021；Zbl 1506.74390） 全文： 内政部 arXiv公司 参考文献： [1] Zoback，M.D.，《储层地质力学》（2010），剑桥大学出版社·Zbl 1189.74005号 [2] Rutqvist，J.，《深层沉积地层中CO（{}_2）储存的地质力学》，《岩土工程》。地质。工程师，30，3，525-551（2012） [3] Majer，E.L。；巴里亚·R。；斯塔克，M。；奥茨，S。；Bommer，J。；B.史密斯。；Asanuma，H.，与增强地热系统相关的诱发地震活动，地热学，36，3，185-222（2007） [4] Cowin，S.C.，《骨质疏松性》，J.Biomech。，32117-238（1999年） [5] 巴迪亚，S。；奎尼，A。；Quarteroni，A.，用于模拟流体-弹性介质相互作用的耦合Biot和Navier-Stokes方程，J.Compute。物理。，228, 21, 7986-8014 (2009) ·Zbl 1391.74234号 [6] Mallison，B。；剑，C。；Viard，T。；Milliken，W。；Cheng，A.，《模拟复杂储层的非结构化剖分网格》，SPE J.，19，02，340-352（2014） [7] Beirão da Veiga，L.公司。；Lipnikov，K。；Manzini，G.，《椭圆问题的模拟有限差分法》，第11卷（2014），Springer·Zbl 1286.65141号 [8] Beirão da Veiga，L.公司。；布雷齐，F。；马里尼，L.D。；Russo，A.，多边形网格上一般二阶椭圆问题的虚拟元方法，数学。模型方法应用。科学。，2015年4月26日，729-750·Zbl 1332.65162号 [9] Di Pietro，D.A。；Droniou，J.，《多边形网格的混合高阶方法》，第19卷（2019年），施普林格出版社 [10] Droniou，J.博士。；Eymard，R。；加洛特，t。；Guichard，C。；Herbin，R.，《梯度离散化方法》，第82卷（2018年），施普林格出版社·Zbl 1435.65005号 [11] Manzini，G。；Russo，A。；Sukumar，N.，《多边形和多面体有限元方法的新观点》，数学。模型方法应用。科学。，1665-1699年8月24日（2014年）·Zbl 1291.65322号 [12] 博蒂，M。；Di Pietro，D.A。；Sochala，P.，非线性多孔弹性的混合高阶离散方法，计算。方法应用。数学。，20, 2, 227-249 (2020) ·Zbl 1437.76027号 [13] Nordbotten，J.M.，《Biot方程的稳定单元中心有限体积离散化》，SIAM J.Numer。分析。，54, 2, 942-968 (2016) ·Zbl 1382.76187号 [14] Terekhov，K.M.，非均质各向异性孔隙力学问题的以细胞为中心的有限体积法，J.Compute。申请。数学。，365，第112357条pp.（2020）·Zbl 1423.76308号 [15] Beirão da Veiga，L.公司。；洛瓦迪纳，C。；Mora，D.，多面体网格上弹性和非弹性问题的虚拟单元法，计算。方法应用。机械。工程，295327-346（2015）·Zbl 1423.74120号 [16] Cangiani，A。；Georgoulis，E.H。；Pryer，T。；Sutton，O.J.，虚拟单元法的后验误差估计，数值。数学。，137, 4, 857-893 (2017) ·兹比尔1384.65079 [17] Berrone，S。；Borio，A.，虚拟元素法中形状不良多边形元素中的正交多项式，有限元。分析。设计。，129, 14-31 (2017) [18] 南卡罗来纳州布伦纳。；关，Q。；Sung，L.-Y.，虚拟元素方法的一些估计，计算。方法应用。数学。，17, 4, 553-574 (2017) ·Zbl 1434.65237号 [19] 贝内代托，M。；Berrone，S。；A.鲍里奥。；Pieraccini，S。；Scialó，S.，对流扩散问题虚元公式的保序SUPG稳定，计算。方法应用。机械。工程，311，18-40（2016）·Zbl 1439.76051号 [20] Berrone，S。；Borio，A.，虚拟元素法的残差后验误差估计，数学。模型方法应用。科学。，1423-1458年8月27日（2017）·Zbl 1376.65136号 [21] Andersen，O。；尼尔森，H.M。；Raynaud，X.，油藏模型地质力学模拟的虚拟单元法，计算。地质科学。，21, 5-6, 877-893 (2017) ·Zbl 1401.74208号 [22] 库莱特，J。；失败，I。；Girault，V。；盖伊，N。；Nataf，F.，一种使用虚拟单元法和有限体积的多孔弹性全耦合方案，计算。地质科学。，1-23 (2019) [23] 贝内代托，M。；Berrone，S。；A.鲍里奥。；Pieraccini，S。；Scialå，S.，用于离散裂缝网络模拟的混合砂浆虚拟单元方法，J.Comput。物理。，306, 148-166 (2016) ·Zbl 1351.76048号 [24] 贝内代托，M。；Berrone，S。；Scialó，S.，一种使用虚拟元素法求解离散裂隙网络流动的全局协调方法，有限元。分析。设计。，109, 23-36 (2016) [25] Lie，K.-A.，《使用MATLAB/GNU Octave进行油藏模拟简介：MATLAB油藏模拟工具箱（MRST）用户指南》（2019年），剑桥大学出版社·Zbl 1425.76001号 [26] 爱德华兹，M.G。；罗杰斯，C.F.，通用张量压力方程的有限体积离散化，计算。地质科学。，2, 4, 259-290 (1998) ·Zbl 0945.76049号 [27] Aavatsmark，I.，四边形网格多点通量近似介绍，计算。地质科学。，6, 3-4, 405-432 (2002) ·Zbl 1094.76550号 [28] Aarnes，J.E。；克罗格斯塔德，S。；Lie，K.-A.，角点网格上的多尺度混合/模拟方法，计算。地质科学。，12, 3, 297-315 (2008) ·Zbl 1259.76065号 [29] Droniou，J.博士。；Eymard，R。；加洛特，t。；Herbin，R.，《模拟有限差分、混合有限体积和混合有限体积方法的统一方法》，数学。模型方法应用。科学。，20, 02, 265-295 (2010) ·Zbl 1191.65142号 [30] Abushaikha，美国。；Terekhov，K.M.，《全张量渗透率通用地下油藏模拟的全隐式模拟有限差分格式》，J.Compute。物理。，406，第109194条pp.（2020）·Zbl 1453.76120号 [31] Eymard，R。；Guichard，C。；Herbin，R.，多孔介质中扩散流动的小型3D方案，ESAIM数学。模型。数字。分析。，46, 2, 265-290 (2012) ·Zbl 1271.76324号 [32] Klemetsdal，Ø。美国。；莫纳，O。；雷诺德，X。；Lie，K.-A.，多面体网格上椭圆问题一致离散化的比较，（复杂应用有限体积国际会议（2020），Springer），585-594·Zbl 1454.65147号 [33] 博纳尔迪，F。；Brenner，K。；Droniou，J.博士。；Masson，R.，裂隙多孔介质中两相流与机械变形耦合的梯度离散化（2020），arXiv:2004.09860 [34] Settari，A。；Mourits，F.，《储层与地质力学耦合模拟系统》，SPE J.，3，03，219-226（1998） [35] Babuška，I.，有限元方法的误差界限，数值。数学。，16, 4, 322-333 (1971) ·Zbl 0214.42001号 [36] Brezzi，F.，关于拉格朗日乘子引起的鞍点问题的存在性、唯一性和近似性，RAIRO Anal。编号。，8，R2，129-151（1974）·Zbl 0338.90047号 [37] Beirão da Veiga，L.公司。；Lipnikov，K.，带选定边气泡的Stokes问题的模拟离散化，SIAM J.Sci。计算。，32, 2, 875-893 (2010) ·Zbl 1352.76021号 [38] 卡马戈，J.T。；怀特，J.A。；Borja，R.I.，多相孔隙力学混合有限元/有限体积离散的宏元稳定性，计算。地质科学。，1-18 (2020) [39] M.Frigo，N.Castelletto，M.Ferronato，J.A.White，混合三场混合有限元耦合孔隙力学的高效求解器，计算。数学。申请。http://dx.doi.org/10.1016/j.camwa.2020.07.010。 ·兹比尔1524.76451 [40] Mandel，J.，《土壤综合》（Etude mathématique），《Géotechnique》，第3、7、287-299页（1953年） [41] Talischi，C。；佩雷拉，A。；Paulino，G.H。；梅内泽斯，I.F.M。；Carvalho，M.S.，《不可压缩流体流动的多边形有限元》，国际。J.数字。《液体方法》，74，2，134-151（2014）·Zbl 1455.76095号 [42] Biot，M.A.，《三维固结的一般理论》，J.Appl。物理。，12, 2, 155-164 (1941) ·肯尼迪67.0837.01 [43] Coussy，O.，《多孔力学》（2004），John Wiley&Sons [44] Frijns，A.J.H.，《应用于软骨组织的四组分混合物理论：数值建模和实验》（2000年），荷兰埃因霍温理工大学（博士论文）·Zbl 0966.9202号 [45] Lipnikov，K.，多孔弹性中Biot模型的数值方法（2002），休斯顿大学（博士论文） [46] Phillips，P.J。；Wheeler，M.F.，多孔弹性的混合和连续Galerkin有限元方法的耦合I：时间连续情况，计算。地质科学。，11, 2, 131-144 (2007) ·Zbl 1117.74015号 [47] Phillips，P.J。；Wheeler，M.F.，多孔弹性的混合和连续Galerkin有限元方法的耦合II：离散时间情况，计算。地质科学。，11, 2, 145-158 (2007) ·Zbl 1117.74016号 [48] Jha，B。；Juanes，R.，《流动与油藏地质力学耦合模拟的局部保守有限元框架》，《岩土学报》。，2, 3, 139-153 (2007) [49] Ferronato，M。；北卡罗来纳州卡斯特莱托。；Gambolia，G.，《Biot固结的全耦合三维混合有限元模型》，J.Compute。物理。，229, 12, 4813-4830 (2010) ·Zbl 1305.76055号 [50] 哈加，J.B。；Osnes，H。；Langtangen，H.P.，《关于低渗透和低压缩多孔介质中压力振荡的原因》，国际期刊数值。分析。方法地质力学。，36, 12, 1507-1522 (2012) [51] 两者，J.W。；博雷加莱斯，M。；Nordbotten，J.M。；库马尔，K。；Radu，F.A.，异质介质中Biot方程的鲁棒固定应力分裂，应用。数学。莱特。，68110-108（2017）·Zbl 1383.74025号 [52] Bause，M。；Radu，F.A。；Köcher，U.，具有迭代耦合的Biot孔隙弹性系统的时空有限元近似，计算。方法应用。机械。工程，320，745-768（2017）·Zbl 1439.74389号 [53] 胡，X。；罗德里戈，C。；Gaspar，F.J。；Zikatanov，L.T.，多孔弹性中Biot固结模型的非协调有限元方法，J.Compute。申请。数学。，310, 143-154 (2017) ·Zbl 1381.76175号 [54] 罗德里戈，C。；胡，X。；欧姆，P。；阿德勒，J.H。；Gaspar，F.J。；Zikatanov，L.T.，多孔弹性和Stokes方程的新稳定离散化，计算。方法应用。机械。工程，341467-484（2018）·Zbl 1440.76027号 [55] 两者，J.W。；库马尔，K。；Nordbotten，J.M。；Radu，F.A.，Anderson非饱和多孔介质固结的加速固定应力劈裂方案（2018） [56] 洪，Q。；Kraus，J.，Biot固结模型经典三场公式的参数鲁棒稳定性，Electron。事务处理。数字。分析。，48, 202-226 (2018) ·Zbl 1398.65046号 [57] 洪，Q。；克劳斯，J。；林伯里，M。；Wheeler，M.F.，多渗透孔隙弹性系统固定应力分裂迭代法的参数稳健收敛分析，多尺度模型。模拟。，18, 2, 916-941 (2020) ·Zbl 1447.65077号 [58] Eymard，R。；加洛特，t。；Herbin，R.，一般非协调网格上非均匀和各向异性扩散问题的离散化SUSHI：使用稳定化和混合界面的方案，IMA J.Numer。分析。，30, 4, 1009-1043 (2010) ·Zbl 1202.65144号 [59] Geuzaine，C。；Remacle，J.-F.，Gmsh：一个内置预处理和后处理设施的三维有限元网格生成器，国际。J.数字。方法工程，79，11，1309-1331（2009）·Zbl 1176.74181号 [60] Talischi，C。；Paulino，G.H。；佩雷拉，A。；Menezes，I.F.M.，PolyMesher:一个用Matlab编写的多边形元素通用网格生成器，Struct。多磁盘。最佳。，45, 3, 309-328 (2012) ·Zbl 1274.74401号 [61] 布雷齐，F。；Lipnikov，K。；Shashkov，M.，多面体网格上扩散问题的模拟有限差分法的收敛性，SIAM J.Numer。分析。，43, 5, 1872-1896 (2005) ·Zbl 1108.65102号 [62] 布雷齐，F。；Lipnikov，K。；Simoncini，V.，多边形和多面体网格上的模拟有限差分方法家族，数学。模型方法应用。科学。，15, 10, 1533-1551 (2005) ·Zbl 1083.65099号 [63] 博德，L。；马森，R。；洛佩兹，S。；Samier，P.，非等温两相Darcy流问题混合网格上基于面和基于节点的组合离散，ESAIM Math。模型。数字。分析。，53, 4, 1125-1156 (2019) ·Zbl 1435.65192号 [64] 达西，F。；Fumagalli，A。；Losapio，D。；Scialó，S。；斯科蒂，A。；Vacca，G.，《弯曲界面网格的混合虚拟元法》（2020），arXiv:2011.09332 [65] 达西，F。；Vacca，G.，《混合高阶虚元方法的砖：投影仪和微分算子》，应用。数字。数学。，155, 140-159 (2020) ·Zbl 1437.65184号 [66] Beirão da Veiga，L.公司。；Gyrya，V。；Lipnikov，K。；Manzini，G.，多边形网格上stokes问题的拟有限差分方法，J.Comput。物理。，228、19、7215-7232（2009），网址https://www.sciencedirect.com/science/article/pii/S002199910900343X ·Zbl 1172.76032号 [67] Beirão da Veiga，L.公司。；Gyrya，V。；Lipnikov，K。；Manzini，G.，《多边形网格上Stokes问题的模拟有限差分法》，J.Compute。物理。，228, 19, 7215-7232 (2009) ·Zbl 1172.76032号 [68] 北凯奇卡。；Silvester，D.J.，Stokes问题的局部稳定混合有限元方法分析，数学。公司。，58, 197, 1-10 (1992) ·Zbl 0738.76040号 [69] Silvester，D.J.，不可压缩流的最优低阶有限元方法，计算。方法应用。机械。工程，111，3-4，357-368（1994）·Zbl 0844.76059号 [70] Elman，H.C。；西尔维斯特·D·J。；Wathen，A.J.，《有限元和快速迭代求解器：在不可压缩流体动力学中的应用》（2014），牛津大学出版社：牛津大学出版社，美国·Zbl 1304.76002号 [71] Boffi，D。；布雷齐，F。；Fortin，M.，《混合有限元方法和应用》（2013），Springer·Zbl 1277.65092号 [72] 萨阿德，Y。；Schultz，M.H.，GMRES:求解非对称线性系统的广义最小残差算法，SIAM J.Sci。统计计算。，7, 3, 856-869 (1986) ·Zbl 0599.65018号 [73] Kim，J。；Tchelepi，H.A。；Juanes，R.，《耦合流和地质力学序贯方法的稳定性和收敛性：固定应力和固定应变分裂》，计算。方法应用。机械。工程，200,131591-1606（2011）·Zbl 1228.74101号 [74] 怀特，J.A。；北卡罗来纳州卡斯特莱托。；Tchelepi，H.A.，《耦合孔隙力学的分块求解器：统一框架》，计算。方法应用。机械。工程，303，55-74（2016）·Zbl 1425.74497号 [75] 北卡罗来纳州卡斯特莱托。；怀特，J.A。；Ferronato，M.，三场混合有限元耦合孔隙力学的可缩放算法，J.Compute。物理。，327, 894-918 (2016) ·兹比尔1373.76312 [76] 贝克，A.H。；R·D·法尔古特。；科列夫，Tz.V。；Yang，U.M.，超并行计算的多重网格平滑器，SIAM J.Sci。计算。，33, 5, 2864-2887 (2011) ·Zbl 1237.65032号 [77] 瓦西列夫斯基，P.S。；Yang，U.M.，使用加法变量减少代数多重网格中的通信，Numer。线性代数应用。，21, 2, 275-296 (2014) ·Zbl 1340.65304号 [78] 里斯，M。；美国特罗滕贝格。；Winter，G.，关于MGR方法的注释，线性代数应用。，49, 1-26 (1983) ·Zbl 0515.65070号 [79] Bui，Q.M。；Osei-Kuffuor，D。；北卡罗来纳州卡斯特莱托。；White，J.A.，《非均质介质多相孔隙力学的可扩展多重网格简化框架》，SIAM J.Sci。计算。，42、2、B379-B396（2020）·Zbl 1435.65153号 [80] Bui，Q.M。；哈蒙，F。；北卡罗来纳州卡斯特莱托。；Osei-Kuffuor，D。；塞特卡斯特，R。；White，J.，多孔和断裂介质中耦合过程的多重网格还原预处理框架（2021），arXiv:2101.11649 [81] Boffi，D。；博蒂，M。；Di Pietro，D.A.，一般网格上Biot问题的非协调高阶方法，SIAM J.Sci。计算。，38、3、A1508-A1537（2016）·Zbl 1337.76042号 [82] O.阿克塞尔森。；Gustafsson，I.，求解Navier弹性方程的迭代方法，计算。方法应用。机械。工程，15，2，241-258（1978）·Zbl 0402.73028号 [83] Blaheta，R.，线性弹性问题的位移分解-完全因式分解预处理技术，数值。线性代数应用。，1, 2, 107-128 (1994) ·Zbl 0837.65021号 [84] Ruge，J.W。；Stüben，K.，代数多重网格，（多重网格方法（1987），SIAM），73-130 [85] 波义耳，J.W。；Mihajlović，医学博士。；Scott，J.A.，HSL_MI20：三维有限元问题的有效AMG预处理器，国际。J.数字。方法工程，82，64-98（2010）·Zbl 1183.76799号 [86] 北卡罗来纳州卡斯特莱托。；怀特，J.A。；Tchelepi，H.A.，双向耦合孔隙力学固定应力迭代解的精度和收敛性，国际数值杂志。分析。方法地质力学。，1593-1618（2015年） [87] Yi，S.-Y.，Biot固结模型新混合有限元法的收敛性分析，Numer。方法偏微分方程，30，4，1189-1210（2014）·Zbl 1350.74024号 [88] Phillips，P.J。；Wheeler，M.F.，《克服线弹性和多孔弹性中的锁定问题：启发式方法》，计算。地质科学。，13, 1, 5-12 (2009) ·Zbl 1172.74017号 [89] 布雷齐，F。；Lipnikov，K。；Shashkov，M.，曲面多面体网格上扩散问题的模拟有限差分法的收敛性，数学。模型方法应用。科学。，16, 02, 275-297 (2006) ·Zbl 1094.65111号 [90] Beirão da Veiga，L.公司。；Russo，A。；Vacca，G.，带曲边的虚拟单元法，ESAIM数学。模型。数字。分析。，53, 2, 375-404 (2019) ·Zbl 1426.65163号 [91] 达西，F。；Scacchi，S.，鞍点问题三维虚拟元离散化的并行块预条件，计算。方法应用。机械。工程，372，第113424条pp.（2020）·Zbl 1506.65204号 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。