×

多孔介质中两相流和输运的基于聚合的非线性多重网格求解器。 (英语) Zbl 1504.65186号

摘要:提出了一种用于分数流量-压力-饱和度混合公式中两相流和输运的非线性多重网格求解器。该求解器在完全近似格式(FAS)的框架下,扩展了我们之前在非均匀扩散问题的非线性多重网格方面的工作。多重网格层次结构中的粗空间是通过首先聚合自由度,然后求解一些局部流问题来构造的。混合公式和粗空间的选择使我们能够在求解阶段组装粗问题,而无需访问更精细的层次,这对多重网格方法的可扩展性至关重要。具体来说,可以使用预先计算的粗略通量基向量在粗略水平上直接评估迎风通量的自然泛化。所得到的求解器适用于在一般非结构化网格上离散的问题。通过具有挑战性的数值算例,验证了所提出的非线性多重网格解算器与标准单层牛顿法的性能比较。结果表明,所提出的求解器对高度非线性问题具有鲁棒性,并且在高Courant-Freedrichs-Lewy(CFL)数的情况下明显优于Newton方法。

MSC公司:

6500万08 含偏微分方程初值和初边值问题的有限体积法
65M55型 多重网格方法;涉及偏微分方程初值和初边值问题的区域分解
76S05号 多孔介质中的流动;过滤;渗流
76T06型 液-液双组分流动
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 曹,H。;Tchelepi,H.A。;Wallis,J.R。;Yardumian,H.E.,用于油藏模拟的并行可扩展非结构化CPR型线性求解器,(SPE年度技术会议和展览。SPE年度科技会议和展览,OnePetro(2005)),pp.SPE-96809-MS
[2] Henson,V.E.,《非线性问题的多重网格方法:概述》(Bouman,C.A.;Stevenson,R.L.,《计算成像》,《SPIE学报》,第5016卷(2003)),第36-48页
[3] Deufhard,P.,《非线性问题的牛顿方法:仿射不变性和自适应算法》,第35卷(2011年),Springer Science&Business Media·Zbl 1226.65043号
[4] Younis,R.,《油藏模拟软件和求解算法的现代进展》(2011),斯坦福大学
[5] 詹妮·P。;Tchelepi,H.A。;Lee,S.H.,S型通量函数双曲守恒律的无条件收敛非线性解算器,J.Compute。物理。,228, 20, 7497-7512 (2009) ·Zbl 1391.76553号
[6] 王,X。;Tchelepi,H.A.,非均匀多孔介质非线性传输的基于信任域的求解器,J.Compute。物理。,253, 114-137 (2013) ·Zbl 1349.76406号
[7] 李,B。;Tchelepi,H.A.,粘性、浮力和毛细管力存在下多孔介质中多相传输的非线性分析,J.Compute。物理。,297, 104-131 (2015) ·Zbl 1349.76820号
[8] Möyner,O.,基于近似信赖域的三相运输问题非线性求解器,计算。地质科学。,21, 999-1021 (2017) ·Zbl 1396.76061号
[9] 尤尼斯,R。;Tchelepi,H.A。;阿齐兹,K.,《自适应局部连续化-纽顿方法-始终收敛的非线性解算器》,SPE J.,15,02,526-544(2009)
[10] 姜杰。;Tchelepi,H.A.,非均质多孔介质中多相流的基于耗散的延拓方法,J.Compute。物理。,375, 307-336 (2018) ·Zbl 1416.76152号
[11] 郭,F。;Tchelepi,H.A.,多孔介质中非线性多相流的基于势的约化牛顿算法,J.Compute。物理。,2271706-727(2007年)·Zbl 1388.76258号
[12] Natvig,J。;Lie,K.-A.,用具有最佳单元顺序的隐式间断Galerkin格式快速计算多孔介质中的多相流,J.Compute。物理。,227, 24, 10108-10124 (2008) ·Zbl 1218.76029号
[13] 哈蒙,F.P。;Tchelepi,H.A.,三相流全隐式模拟的基于排序的非线性求解器,计算。地质科学。,20, 5, 909-927 (2016) ·Zbl 1391.76805号
[14] Klemetsdal,Ø。美国。;拉斯穆森,A.F。;莫纳,O。;Lie,K.-A.,黑油模型的高阶高效重排序非线性Gauss-Seidel解算器,计算。地质科学。,24, 2, 593-607 (2020) ·Zbl 1434.86015号
[15] 蔡,X.C。;Keyes,D.E.,非线性预处理不精确牛顿算法,SIAM J.Sci。计算。,24, 1, 183-200 (2002) ·Zbl 1015.65058号
[16] 刘,L。;Keyes,D.E.,Field-split预处理不精确牛顿算法,SIAM J.Sci。计算。,37、3、A1388-A1409(2015)·Zbl 1328.65122号
[17] Dolean,V。;甘德,M.J。;Kheriji,W。;郭,F。;Masson,R.,《非线性预处理:如何使用非线性Schwarz方法预处理牛顿法》,SIAM J.Sci。计算。,38、6、A3357-A3380(2016)·Zbl 1352.65326号
[18] 斯科格斯塔德,J.O。;Keilegavlen,E。;Nordbotten,J.M.,多孔介质非线性流动问题的区域分解策略,J.Compute。物理。,234, 439-451 (2013)
[19] 斯科格斯塔德,J.O。;Keilegavlen,E。;Nordbotten,J.M.,多孔介质中两相非线性流动的双尺度预处理,Transp。多孔介质,114,2485-503(2016)
[20] Klemetsdal,Ø。美国。;Moncorgé,A。;莫纳,O。;Lie,K.-A.,Additive Schwarz预处理精确牛顿法作为多相多孔介质流动的非线性预处理剂,(ECMOR XVII会议录,第2020卷(2020),欧洲地球科学家和工程师协会),1-20
[21] Brandt,A.,边值问题的多级自适应解决方案,数学。计算。,31, 138, 333-390 (1977) ·Zbl 0373.65054号
[22] Christensen,M.l.C。;Eskildsen,K.L。;Engsig-Karup,A.P。;Wakefield,M.A.,油藏模拟非线性多重网格,SPE J.,21,3,888-898(2016)
[23] 克里斯滕森医学博士。;瓦西列夫斯基,P.S。;Villa,U.,《利用具有近似特性的AMGe粗空间的非线性多重网格解算器》,J.Compute。申请。数学。,340, 691-708 (2018) ·Zbl 1432.76160号
[24] 托夫特,R。;谎言,K.-A。;Möyner,O.,《储层模拟的完全近似方案》,(挪威信息技术协会会议录,挪威奥斯陆,挪威信息技术委员会会议录(2018))
[25] Lee,C.S。;哈蒙,F.P。;北卡罗来纳州卡斯特莱托。;瓦西列夫斯基,P.S。;White,J.A.,基于非均匀扩散问题局部谱粗化的非线性多重网格,计算。方法应用。机械。工程,372,第113432条pp.(2020)·兹比尔1506.76113
[26] 阿齐兹,K。;Settari,A.,《石油储层模拟》(1979),爱思唯尔应用科学出版社:英国伦敦
[27] Lie,K.-A.,《使用MATLAB/GNU Octave进行油藏模拟简介:MATLAB油藏模拟工具箱(MRST)用户指南》(2019年),剑桥大学出版社·Zbl 1425.76001号
[28] Eymard,R。;加洛特,t。;Herbin,R.,《有限体积方法》(Ciarlet,P.G.;Lions,J.L.,《R^n中方程的求解》(第3部分),科学计算技术(第三部分)。(R^n)中方程的求解(第3部分),科学计算技术(第3章),数值分析手册,第7卷(2000年),爱思唯尔出版社,713-1018·Zbl 0953.00016号
[29] Peaceman,D.W.,储层数值模拟中的井块压力解释,SPE J.,18,3,183-194(1978)
[30] 陈,Z。;Huan,G。;Ma,Y.,多孔介质中多相流的计算方法(2006),工业和应用数学学会·Zbl 1092.76001号
[31] Karypis,G。;Kumar,V.,一种用于划分不规则图的快速高质量多级方案,SIAM J.Sci。计算。,20, 1, 359-392 (1998) ·Zbl 0915.68129号
[32] Aarnes,J.E。;克罗格斯塔德,S。;Lie,K.-A.,基于混合有限元和非均匀粗网格的两相流分层多尺度方法,多尺度模型。模拟。,5, 2, 337-363 (2006) ·Zbl 1124.76022号
[33] Arbogast,T.,异质椭圆问题的混合多尺度方法,(Graham,I.G.;Hou,T.Y.;Lakkis,O.;Scheichl,R.,《多尺度问题的数值分析》(2012),斯普林格·柏林-海德堡:斯普林格尔·柏林-海德堡-柏林,海德堡),243-283·Zbl 1248.65119号
[34] 巴克,A.T。;Gelever,S.V。;Lee,C.S。;奥斯本,S。;Vassilevski,P.S.,图拉普拉斯问题的多级谱粗化及其在油藏模拟中的应用,SIAM J.Sci。计算。,43、4、A2737-A2765(2021)·Zbl 1482.65222号
[35] Wallis,J.R.,作为广义共轭梯度加速度预处理的不完全高斯消去,(SPE油藏模拟研讨会,SPE油藏仿真研讨会,OnePetro(1983))
[36] Wallis,J.R。;肯德尔,R.P。;Little,T.E.,共轭残差法的约束残差加速,(SPE油藏模拟研讨会。SPE油藏仿真研讨会,OnePetro(1985))
[37] Lacroix,S。;瓦西列夫斯基,Y。;Wheeler,M.F.,隐式并行精确油藏模拟器(IPARS)中的解耦预条件,数值。线性代数应用。,8, 8, 537-549 (2001) ·Zbl 1071.76583号
[38] Scheichl,R。;马森,R。;Wendebourg,J.,沉积盆地模拟的解耦和区块预处理,计算。地质科学。,7, 4, 295-318 (2003) ·Zbl 1076.76070号
[39] hypre:可伸缩线性解算器和多重网格方法
[40] Stüben,K。;克莱斯,T。;克利,H。;卢,B。;Wheeler,M.F.,《油藏模拟中高效求解全隐式公式的代数多重网格方法(AMG)》,(SPE油藏模拟研讨会。SPE油藏仿真研讨会,OnePetro(2007))
[41] 贝克,A.H。;R·D·法尔古特。;科列夫,Tz.V。;Yang,U.M.,超并行计算的多重网格平滑器,SIAM J.Sci。计算。,332864-2887(2011年)·Zbl 1237.65032号
[42] Christie,文学硕士。;Blunt,M.J.,第十个SPE比较解决方案项目:升级技术的比较,SPE Reserve。评估。工程,4,4,308-317(2001)
[43] 詹森,J.D。;丰塞卡,R.-M。;Kahrobaei,S。;Siraj,M.M。;范·埃森,G.M。;Van den Hof,P.M.J.,《卵形模型——油藏模拟的地质集成》,《地质学》。数据J.,1,2192-195(2014)
[44] Manzocchi,T。;卡特,J.N。;斯科斯塔德,A。;Fjellvoll,B。;斯蒂芬·K·D。;豪厄尔,J.A。;Matthews,J.D。;沃尔什·J·J。;Nepveu,M。;Bos,C.,地质不确定性对断层和未断层浅海油藏产量影响的敏感性:目标和方法,Pet。地质科学。,14, 1, 3-15 (2008)
[45] Cao,H.,《通用模拟器技术的发展》(2002),斯坦福大学
[46] 模块化有限元方法库
[47] smoothG:混合图Laplacian上尺度和解算器
[48] GLVis:openGL有限元可视化工具
[49] 克罗格斯塔德,S。;Lie,K.-A。;莫纳,O。;尼尔森,H.M。;雷诺德,X。;Skaflestad,B.,MRST-AD-油藏模拟问题快速原型制作和评估的开源框架,(SPE油藏模拟研讨会。SPE油藏仿真研讨会,OnePetro(2015))
[50] Aarnes,J.E.,《关于在油藏模拟中使用混合多尺度有限元方法以获得更大的灵活性、更快或更高的精度,多尺度模型》。模拟。,2, 3, 421-439 (2004) ·兹比尔1181.76125
[51] 基普,V。;阿尔内斯,J.E。;Lie,K.-A.,多孔介质流动椭圆问题多尺度方法的比较,计算。地质科学。,12, 3, 377-398 (2008) ·Zbl 1259.76047号
[52] Lee,S.H。;沃尔夫斯泰纳,C。;Tchelepi,H.A.,多孔介质中多相流的多尺度有限体积公式:可压缩、重力三相流的黑油公式,计算。地质科学。,12, 3, 351-366 (2008) ·Zbl 1259.76049号
[53] 卡宾斯基,S。;索林,美国。;Radu,F.A.,多孔介质中动态毛细效应两相流的内罚间断Galerkin方法线性化方案分析,国际期刊数值。方法工程,112,6,553-577(2017)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。