×
罗，李;刘璐璐;蔡晓川;戴维·凯斯（David E.Keyes）。

三维非结构网格上多孔介质两相流的全隐式混合两层区域分解算法。 （英语） Zbl 1435.76038号

J.计算。物理学。 409，文章ID 109312，19 p.（2020）.
摘要：由于算子的非线性和材料系数的高度不均匀性，多孔介质中地下流动的模拟很困难。本文基于重叠区域分解方法，提出了一种可扩展的不可压缩两相流全隐式求解器。具体地说，使用带有解析雅可比矩阵的不精确Newton-Krylov算法求解三维非结构网格上控制方程的间断Galerkin离散化所产生的非线性系统。线性雅可比系统通过加法Schwarz算法进行预处理，该算法自然适合并行计算。我们提出了由嵌套粗空间组成的加性Schwarz预条件器的混合两级版本，以提高经典一级版本的鲁棒性和可扩展性。在粗水平上，通过使用带有一级预条件的GMRES求解由粗网格上问题的相同离散化产生的较小线性系统，直到达到相对公差。数值实验证明了所提出的求解器在三维非均匀介质问题中的有效性和效率。我们还报告了所提算法在具有8192个处理器核的超级计算机上的并行可伸缩性。

引用于11文件

MSC公司：

76M10个 有限元方法在流体力学问题中的应用
76S05号 多孔介质中的流动；过滤；渗流
76T06型 液-液双组分流动
2005年5月 并行数值计算

关键词：

多孔介质中的两相流;二级区域分解;间断Galerkin;非结构化网格;不精确牛顿;并行可扩展性

软件：

libMesh（libMesh）;MRST公司;METIS公司;PETSc公司;Matlab公司
PDF格式BibTeX公司 XML格式引用
\textit{L.Luo}等人，J.Compute。物理学。409，文章ID 109312，19 p.（2020；Zbl 1435.76038）
全文： 内政部 链接

参考文献：

[1] Arbogast，T。；Juntunen，M。；Pool，J。；Wheeler，M.F.，多孔介质中两相流的间断Galerkin方法，增强速度和连续毛细管压力，计算。岩土工程。，17, 1055-1078 (2013) ·Zbl 1393.76059号
[2] 巴格格，A。；阿特金斯，H。；Keyes，D.E.，非连续Galerkin方法的并行实现，（并行CFD’99（1999）论文集），115-122
[3] 巴莱，S。；Abhyankar，S。；M.F.亚当斯。；Brown，J。；布鲁纳，P。；Buschelman，K。；艾伊霍特，V。；格罗普，W.D。；考希克，D。；Knepley，M.G。；McInnes，L.C.公司。；鲁普，K。；B.F.史密斯。；Zhang，H.，《PETSc用户手册》（2019年），阿贡国家实验室
[4] Bastian，P。；Rivière，B.，多孔介质中两相流的间断Galerkin方法（2004），海德堡大学国际水利研究所，技术报告2004-28
[5] Bastian，P.，多孔介质中不连续毛细压力两相流的完全耦合不连续Galerkin方法，计算。地质科学。，18, 779-796 (2014) ·兹比尔1392.76072
[6] 巴克，A.T。；Cai，X.C.，流体-结构相互作用的两级牛顿和混合Schwarz预条件，SIAM J.Sci。计算。，32, 2390-2417 (2010) ·Zbl 1214.92014年
[7] 巴克，A.T。；南卡罗来纳州布伦纳。；Sung，L.-Y.，椭圆问题局部间断Galerkin方法的重叠Schwarz区域分解预条件，J.Numer。数学。，19, 165-187 (2011) ·Zbl 1232.65164号
[8] 伯纳基，M。；费佐伊，L。；Lanteri，S。；Piperno，S.，用于计算三维波传播问题的并行间断Galerkin非结构化网格解算器，应用。数学。型号。，30, 744-763 (2006) ·Zbl 1101.78009号
[9] 蔡，X.C。；格罗普，W.D。；Keyes，D.E。；梅尔文·R·G。；Young，D.P.，跨声速全势方程的并行Newton-Krylov-Schwarz算法，SIAM J.Sci。计算。，192, 46-65 (1998)
[10] 蔡，X.C。；Sarkis，M.，一般稀疏线性系统的限制加性Schwarz预条件，SIAM J.Sci。计算。，21, 792-797 (1999) ·Zbl 0944.65031号
[11] Chen，Z。；Ewing，R.，《从单相流到成分流：混合有限元的适用性》，Transp。多孔介质，57，225-242（1997）
[12] Chen，Z。；Huan，G。；Li，B.，多孔介质中两相流的改进IMPES方法，Transp。多孔介质，54，361-376（2004）
[13] Chen，Z。；Huan，G。；Ma，Y.，多孔介质中多相流的计算方法（2006），SIAM：美国宾夕法尼亚州费城SIAM·Zbl 1092.76001号
[14] 科尔曼，T.F。；Moré，J.J.，稀疏雅可比矩阵的估计和图着色问题，SIAM J.Numer。分析。，20, 187-209 (1983) ·Zbl 0527.65033号
[15] 库西尼，M。；卢基亚诺夫，A.A。；Natvig，J。；Hajibeygi，H.，多孔介质中多相流全隐式模拟的约束压力剩余多尺度（CPR-MS）方法，J.Compute。物理。，299, 472-486 (2015) ·Zbl 1351.76055号
[16] Dawson，C.N。；克利，H。；惠勒，M.F。；Woodward，C.S.，用预处理Newton-Krylov解算器求解两相流的并行、隐式、单元中心方法，Comput。地质科学。，1, 215-249 (1997) ·Zbl 0941.76062号
[17] 丹尼斯，J.E。；Schnabel，R.B.，《无约束优化和非线性方程的数值方法》（1996），SIAM:SIAM Philadelphia·Zbl 0847.65038号
[18] Dogru，A.H。；Fung，L.S.K。；美国米迪亚。；Al Shaalan，T.M.，技术经理。；Pita，J.A。；HemanthKumar，K。；苏海杰。；Tan，J.C.T。；霍伊，H。；Dreiman，W.T。；哈恩，W.A。；阿尔·哈比（Al-Harbi，R.）。；Al-Youbi，A。；新墨西哥州Al-Zamel。；梅兹加尼，M。；Al-Mani，T.，下一代大型油藏并行油藏模拟器，（SPE 119272在SPE油藏模拟研讨会上发表。SPE 119722在德克萨斯州SPE油藏仿真研讨会上发表，2009年2月2日至4日）
[19] 艾森斯塔特，S.C。；Walker，H.F.，《用不精确牛顿法选择强迫项》，SIAM J.Sci。计算。，17, 16-32 (1996) ·Zbl 0845.65021号
[20] Epshteyn，Y。；Rivière，B.，两相流的全隐式间断有限元方法，应用。数字。数学。，57, 383-401 (2007) ·Zbl 1370.76085号
[21] Epshteyn，Y。；Rivière，B.，《不可压缩两相流hp间断Galerkin方法分析》，J.Compute。申请。数学。，225, 487-509 (2009) ·Zbl 1157.76024号
[22] Ern，A。；莫佐列夫斯基，I。；Schuh，L.，非连续毛细压力非均匀多孔介质中两相流的间断Galerkin近似，计算。方法应用。机械。工程，1991491-1501（2010）·Zbl 1231.76143号
[23] Eslinger，O.，《应用于两相气水流动问题的间断Galerkin有限元方法》（2005），德克萨斯大学奥斯汀分校博士论文
[24] 尤因，R。；Heinemann，R.，混合有限元方法在成分模拟中的应用，以减少数值弥散（油藏模拟研讨会，编号SPE12267（1983））
[25] 格里斯，S。；Stüben，K。；Brown，G.L。；陈，D。；Collins，D.A.，在完全隐式油藏模拟中有效应用代数多重网格的预处理，SPE J.，19726-736（2014）
[26] Griewank，A.，评估导数：算法微分的原理和技术（2000），SIAM:SIAM Philadelphia·Zbl 0958.65028号
[27] 格罗普，W.D。；Keyes，D.E。；麦克因斯，L。；Tidriri，M.，并行隐式CFD的全球化Newton-Krylov-Schwarz算法和软件，Int.J.高性能计算。申请。，14, 102-136 (2000)
[28] Gunasekera，D。；Childs，P。；赫林，J。；Cox，J.，通用多面体网格的多点通量离散格式，（国际石油和天然气会议和展览，SPE48855（1998））
[29] Hoteit，H。；Firoozabadi，A.，具有不同毛细管压力的非均质渗透介质中两相流的数值模拟，高级水资源。，31, 56-73 (2008)
[30] Hou，J。；陈，J。；Sun，S。；Chen，Z.，非均匀介质中两相流的自适应混合和惩罚间断Galerkin方法，J.Compute。申请。数学。，307, 262-283 (2016) ·Zbl 1382.76166号
[31] Karypis，G。；Kumar，V.，MeTis:《非结构化图形划分和稀疏矩阵排序系统》，4.0版（2009），明尼苏达大学：明尼阿波利斯大学
[32] 柯克，B.S。；彼得森，J.W。；Stogner，R.H。；Carey，G.F.，libMesh：一个用于并行自适应网格细化/粗化模拟的C++库，工程计算。，22, 237-254 (2006)
[33] Klieber，W。；Rivière，B.，用间断Galerkin方法自适应模拟两相流，计算。方法应用。机械。工程，196，404-419（2006）·Zbl 1120.76327号
[34] Knoll，D.A。；Keyes，D.E.，《无雅可比牛顿-克利洛夫方法：方法和应用调查》，J.Compute。物理。，193, 357-397 (2004) ·Zbl 1036.65045号
[35] 寇，J。；Sun，S.，非均匀介质中不可压缩两相流的间断Galerkin方法的收敛性，SIAM J.Numer。分析。，51, 3280-3306 (2013) ·Zbl 1282.76124号
[36] 寇，J。；Sun，S.，多孔介质中不可压缩两相流的迎风间断Galerkin质量守恒方法，数值。方法部分差异。Equ.、。，30, 1674-1699 (2014) ·Zbl 1308.76169号
[37] Lie，K.A.，《使用MATLAB/GNU Octave进行油藏模拟简介：MATLAB油藏模拟工具箱（MRST）用户指南》（2019），剑桥大学出版社：剑桥大学出版社·Zbl 1425.76001号
[38] 刘，H。；王凯。；Chen，Z.，用于并行油藏模拟的一系列约束压力剩余预条件，数值。线性代数应用。，23, 120-146 (2016) ·Zbl 1413.65045号
[39] 刘，L。；Keyes，D.E。；Sun，S.，用加性Schwarz预处理不精确牛顿法进行全隐式两相油藏模拟，（SPE Reserve.Charact.Simul.Conf.Exhib.（2013），石油工程师学会）
[40] 罗，L。；Shiu，W.S。；陈，R。；Cai，X.C.，《人类动脉狭窄血流问题的非线性消除预处理不精确牛顿法》，J.Compute。物理。，399，第108926条pp.（2019）·Zbl 1453.65109号
[41] Michel，A.，多孔介质中两相不可压缩流动模拟的有限体积格式，SIAM J.Numer。分析。，41, 1301-1317 (2003) ·Zbl 1049.35018号
[42] 蒙特亚古多，J.E.P。；Firoozabadi，A.，《裂缝性和非裂缝性介质注水模拟的完全隐式和IMPES公式的比较》，国际期刊数值。方法工程，69698-728（2007）·Zbl 1194.76160号
[43] 奈尔·R·D。；Choi，H.W。；Tufo，H.M.，可扩展高阶间断Galerkin大气动力核心的计算方面，计算。流体，30309-319（2009）·Zbl 1237.76129号
[44] Nayagum，D。；谢弗，G。；Mosé，R.，使用混合和不连续有限元模拟多孔介质中的两相不可压缩流动，计算。地质科学。，8, 49-73 (2004) ·Zbl 1221.76119号
[45] Peaceman，D.W.，《数值油藏模拟基础》（1977年），爱思唯尔：爱思唯尔阿姆斯特丹
[46] Reichenberger，V。；雅各布斯，H。；Bastian，P。；Helmig，R.，裂隙多孔介质中两相流的混合维有限体积法，高级水资源。，29, 1020-1036 (2006)
[47] 里维埃，B。；惠勒，M.F。；Banas，K.，第二部分。间断Galerkin方法应用于多孔介质中的单相流动，计算。地质科学。，4, 337-349 (2000) ·Zbl 1049.76565号
[48] 里维埃，B。；Wheeler，M.F.，多孔介质中的混合驱替，（Hassanizadeh，S.M.；Schotting，R.J.，《第十四届水资源计算方法国际会议论文集》（2002），Elsevier:Elsevier Amsterdam）·Zbl 1060.76072号
[49] Saad，Y.，《一种灵活的内外预处理GMRES算法》，SIAM J.Sci。计算。，14, 461-469 (1993) ·Zbl 0780.65022号
[50] Saad，Y.，《稀疏线性系统的迭代方法》（2003），SIAM:SIAM Philadelphia，PA·Zbl 1002.65042号
[51] 沙迪德，J.N。；杜米纳罗，R.S。；Devine，K.D。；Hennigan，G.L。；Lin，P.T.，有限元传输/反应模拟中全耦合区域分解预处理器的性能，J.Compute。物理。，205, 24-47 (2005) ·Zbl 1087.76069号
[52] B.史密斯。；比约尔斯塔德，P。；Gropp，W.，《区域分解：椭圆偏微分方程的并行多层方法》（1996），剑桥大学出版社·Zbl 0857.65126号
[53] Sun，S。；Wheeler，M.F.，耦合流和反应输运问题的间断Galerkin方法，应用。数字。数学。，52, 273-298 (2005) ·Zbl 1079.76584号
[54] Sun，S。；Wheeler，M.F.，多孔介质中反应输运的对称和非对称间断Galerkin方法，SIAM J.Numer。分析。，43, 195-219 (2005) ·Zbl 1086.76043号
[55] 托塞利，A。；Widlund，O.，《区域分解方法-算法和理论》（2005年），施普林格出版社：施普林格-柏林·Zbl 1069.65138号
[56] 王凯。；刘，H。；Chen，Z.，分布式存储并行计算机上的可扩展并行黑油模拟器，J.Compute。物理。，301, 19-34 (2015) ·Zbl 1349.76833号
[57] 杨，C。；Cai，X.-C.，立方体上非光滑地形浅水方程隐式解的并行多级方法，J.Comput。物理。，230, 2523-2539 (2011) ·Zbl 1316.76016号
[58] 杨，C。；Cai，X.C.，大气流动中尺度非静力模拟的可缩放全隐式可压缩Euler解算器，SIAM J.Sci。计算。，36，S23-S47（2014）·Zbl 1305.86005号
[59] Yang，H。；Cai，X.C.，非定常不可压缩流分布式控制的并行全隐式双网格Lagrange-Newton-Krylov-Schwarz方法，国际数值杂志。方法流体，72，1-21（2013）·Zbl 1455.76050号
[60] Yang，H。；杨，C。；Sun，S.，多孔介质中两相流问题的带非线性消除的Active-set缩减空间方法，SIAM J.Sci。计算。，38，B593-B618（2016）·兹比尔1383.76384
[61] Yang，H。；Sun，S。；杨，C.，多孔介质两相流变分不等式解的非线性预处理半光滑牛顿方法，J.Compute。物理。，332, 1-20 (2017) ·兹比尔1378.76115
[62] Yang，H。；Sun，S。；李毅。；Yang，C.，并行计算机上油藏模拟的可扩展全隐式框架，计算。方法应用。机械。工程，330，334-350（2018）·Zbl 1439.76103号
[63] Yang，H。；Sun，S。；李毅。；Yang，C.，《油藏黑油模型的完全隐式约束保持模拟器》，J.Compute。物理。，396, 347-363 (2019) ·Zbl 1452.76138号
[64] Yang，H。；Sun，S。；李毅。；Yang，C.，多组分可压缩流完全隐式互补公式的并行油藏模拟器，计算。物理学。社区。，244, 2-12 (2019) ·Zbl 07674826号
[65] 齐达内，A。；Firoozabadi，A.，多孔介质中多组分可压缩两相流的隐式数值模型，Adv.Water Resour。，85，64-78（2015）
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配，并且可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。