×

变密度两相流相场法的压力修正和保界离散化。 (英语) Zbl 07524769号

小结:在本文中,我们提出了一种有效的数值算法来求解模拟不同密度两相流的含时Cahn-Hilliard-Navier-Stokes方程。投影法中的压力校正步骤包括一个具有修正右侧的泊松问题。空间离散化基于分段线性或分段二次多项式的间断Galerkin方法。这项工作的一个重要贡献是制定了通量和斜率限制技术,成功地消除了顺序参数中的整体偏移、超调和欠调。证明了离散序参数的保有界性。多个数值结果表明,所提出的数值算法对于模拟多孔结构和数字岩石中的两组分非混溶流动是有效和稳健的。

MSC公司:

7.6亿 流体力学基本方法
76天xx 不可压缩粘性流体
6500万 偏微分方程、初值和含时初边值问题的数值方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 莱克莱尔,S。;Parmigiani,A。;马拉菠菜。;肖伯德,B。;Latt,J.,多孔介质中非混溶两相孔隙尺度吸渗和排水的广义三维晶格Boltzmann色粒法,Phys。E版,95,3,第033306条,pp.(2017)
[2] Raeini,A.Q。;杨,J。;邦迪诺,I。;Bultreys,T。;布伦特,M.J。;Bijeljic,B.,使用两相流显微CT图像验证广义孔隙网络模型,Transp。多孔介质,130,2,405-424(2019)
[3] Kim,J.,多组分流体流动的相场模型,Commun。计算。物理。,12, 3, 613-661 (2012) ·Zbl 1373.76030号
[4] 沈,J.,用相场方法对两相不可压缩流动进行建模和数值近似,(材料模拟的多尺度建模和分析(2012),《世界科学》,147-195
[5] Feng,X.,两相流体流动Navier-Stokes-Cahn-Hilliard扩散界面模型的全离散有限元近似,SIAM J.Numer。分析。,44, 1049-1072 (2006) ·Zbl 1344.76052号
[6] Bao,K。;Shi,Y。;Sun,S。;Wang,X.-P.,移动接触线问题Cahn-Hilliard和Navier-Stokes耦合系统数值解的有限元方法,J.Compute。物理。,231, 24, 8083-8099 (2012) ·Zbl 1284.65131号
[7] 迪格尔,A.E。;王,C。;王,X。;Wise,S.M.,Cahn-Hilliard-Navier-Stokes系统二阶精确有限元方法的收敛分析和误差估计,Numer。数学。,1373495-534(2017)·Zbl 1523.65081号
[8] 刘,C。;Riviere,B.,Cahn-Hilliard-Navier-Stokes方程间断Galerkin方法的先验误差分析,CSIAM Trans。申请。数学。,1, 1, 104-141 (2020)
[9] 丁·H。;Spelt,P.,接触线运动扩散界面模拟中的润湿条件,Phys。E版,75,第046708条,pp.(2007)
[10] Abels,H。;Garcke,H。;Grün,G.,《不同密度不可压缩两相流的热力学一致、框架无关扩散界面模型》,数学。模型方法应用。科学。,22,03,第1150013条pp.(2012)·兹比尔1242.76342
[11] Aland,S。;Voigt,A.,《二维气泡动力学扩散界面模型的基准计算》,国际期刊Numer。《液体方法》,69,3,747-761(2012)
[12] Shokrpour-Roudbari,M。;van Brummelen,E.H。;Verhoosel,C.V.,多孔介质中两相流的多尺度扩散界面模型,计算。流体,141212-222(2016)·Zbl 1390.76355号
[13] Van Brummelen,E。;Shokrpour-Roudbari,M。;Van Zwieten,G.,基于Navier-Stokes-Cahn-Hilliard方程的弹性毛细模拟,(计算流体-结构相互作用和流动模拟进展(2016)),451-462·Zbl 1356.74068号
[14] 吉蒙德,J.-L。;Minev,P。;Shen,J.,不可压缩流投影方法概述,计算。方法应用。机械。工程,195,44-47,6011-6045(2006)·Zbl 1122.76072号
[15] 刘建国。;Shu,C.-W.,二维不可压缩流动的高阶间断Galerkin方法,J.Compute。物理。,160, 577-596 (2000) ·Zbl 0963.76069号
[16] 沙巴齐,K。;菲舍尔,P。;Ethier,C.,非定常不可压缩Navier-Stokes方程的高阶间断Galerkin方法,J.Compute。物理。,222, 391-407 (2007) ·Zbl 1216.76034号
[17] 博蒂,L。;Di Pietro,D.A.,具有间断速度和连续压力的对流主导的不可压缩流的压力校正方案,J.Compute。物理。,230572-585(2021)·Zbl 1283.76030号
[18] 董,S。;Shen,J.,《大密度比两相不可压缩流相场模拟中涉及常数矩阵的时间步进格式》,J.Compute。物理。,231, 17, 5788-5804 (2012) ·Zbl 1277.76118号
[19] 多德,M。;Ferrante,A.,《不可压缩双流体流动的快速压力校正方法》,J.Compute。物理。,273, 416-434 (2014) ·Zbl 1351.76161号
[20] Lee,A。;Münch,A。;Süli,E.,具有退化迁移率的Cahn-Hilliard方程的Sharp界面极限,SIAM J.Appl。数学。,76, 433-456 (2016) ·Zbl 1343.35129号
[21] 弗兰克·F。;刘,C。;Alpak,F.O。;Riviere,B.,微孔域上具有简并迁移率的对流Cahn-Hilliard方程的有限体积/间断Galerkin方法,源于显微CT成像,计算。地质科学。,22, 2, 543-563 (2018) ·Zbl 1405.65104号
[22] 弗兰克·F。;拉普,A。;Kuzmin,D.,守恒定律DG离散的束缚守恒通量限制方案及其在Cahn-Hilliard方程中的应用,计算。方法应用。机械。工程,359,第112665条pp.(2019)·Zbl 1441.76059号
[23] 库兹明,D。;Möller,M.,代数通量校正。I.标量守恒定律,(通量校正输运。原理、算法和应用。科学计算(2005),施普林格),155-206·兹比尔1094.76040
[24] Kuzmin,D.,p-自适应间断Galerkin方法的基于顶点的分层斜率限制器,J.Compute。申请。数学。,233, 12, 3077-3085 (2010) ·Zbl 1252.76045号
[25] Kuzmin,D.,具有可能非正交Taylor基的间断Galerkin近似的斜率限制,国际期刊Numer。《液体方法》,71,9,1178-1190(2013)·Zbl 1431.65173号
[26] 刘,C。;弗兰克·F。;蒂勒,C。;Alpak,F.O。;伯格,S。;查普曼,W。;Riviere,B.,求解多孔介质中粘性对比度Cahn-Hilliard-Navier-Stokes系统的高效数值算法,J.Compute。物理。,400,第108948条pp.(2020)·Zbl 1453.76078号
[27] 刘,C。;Shen,J.,两种不可压缩流体混合物的相场模型及其傅立叶谱方法的近似,Phys。D: 非线性现象。,179, 3, 211-228 (2003) ·Zbl 1092.76069号
[28] 巴达拉西,V。;塞尼塞罗斯,H。;Banerjee,S.,用相场模型计算多相系统,J.Compute。物理。,190, 317-397 (2003) ·Zbl 1076.76517号
[29] Feng,X.,两相流体流动Navier-Stokes-Cahn-Hilliard扩散界面模型的全离散有限元近似,SIAM J.Numer。分析。,44, 3, 1049-1072 (2006) ·兹比尔1344.76052
[30] 弗兰克·F。;刘,C。;Scanziani,A。;Alpak,F.O。;Riviere,B.,相场法锯齿表面上基于能量的平衡接触角边界条件,J.胶体界面科学。,523, 282-291 (2018)
[31] D.J.Eyre,梯度系统的无条件稳定一步格式,1998年,第1-15页,未发表文章。
[32] Girault,V。;里维埃尔,B。;Wheeler,M.,Stokes和Navier-Stokes问题的非重叠区域分解间断Galerkin方法,数学。计算。,74, 249, 53-84 (2005) ·Zbl 1057.35029号
[33] Riviere,B.,《解椭圆型和抛物型方程的间断Galerkin方法:理论与实现》,应用数学前沿(2008),工业和应用数学学会·Zbl 1153.65112号
[34] 库兹明,D。;Gorb,Y.,一种用于处理稠密悬浮液中封闭极限的通量校正传输算法,J.Compute。申请。数学。,236、184944-4951(2012年)·Zbl 1426.76280号
[35] 刘,C。;弗兰克·F。;Alpak,F.O。;Riviere,B.,三维不可压缩Navier-Stokes方程的内罚间断Galerkin方法,用于多孔介质渗透率估算,J.Compute。物理。,396, 669-686 (2019) ·Zbl 1452.76095号
[36] 马加莱蒂,F。;皮卡诺,F。;Chinappi,M。;马里诺,L。;Casciola,C.M.,《二元流体Cahn-Hilliard/Navier-Stokes模型的尖锐界面极限》,《流体力学杂志》。,714, 95 (2013) ·Zbl 1284.76116号
[37] Lake,L.W.,《提高石油采收率》(1989年),Prentice Hall:新泽西州Prentice Hall。
[38] Hilfer,R。;Armstrong,R。;伯格,S。;乔治亚迪斯,A。;Ott,H.,《毛细管饱和度和去饱和度》,Phys。E版,92,第062023条,pp.(2015)
[39] Alpak,F.O。;扎卡鲁迪奥,I。;伯格,S。;Dieterich,J。;Saxena,N.,在通用图形处理单元上使用相场晶格Boltzmann方法在大型数字岩石图像上直接模拟孔隙尺度两相粘毛细流,计算。地质科学。,第23、5、849-880页(2019年)·Zbl 1425.76196号
[40] Stegemeier,G.,《圈闭油饱和度与多孔介质岩石物理性质的关系》(1974年),SPE 4754
[41] Yeganeh,M。;Hegner,J。;Lewandowski,E。;莫汉,A。;湖泊,L.W。;Cherney,D。;Jusufi,A。;Jaishankar,A.,《毛细管去饱和曲线基本原理》(2016),SPE-179574-MS
[42] Armstrong,R.T。;乔治亚迪斯,A。;奥特·H。;Klemin,D。;Berg,S.,《临界毛细管数:用快速X射线计算机显微层析研究的去饱和作用》,地球物理学。Res.Lett.公司。,41,1-6(2014)
[43] Hilfer,R。;Oren,P.,孔隙尺度和场尺度非混溶驱替的量纲分析,Transp。多孔介质,22,1,53-72(1996)
[44] Nystrom,N.A。;莱文,M.J。;罗斯基斯,R.Z。;Scott,J.R.,Bridges:新社区和数据分析的独特灵活的hpc资源,(《2015年XSEDE会议论文集:增强网络基础设施带来的科学进步》,XSEDE’15(2015)),第30条pp。
[45] 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
[46] 蒂勒,C。;Araya-Polo,M。;Alpak,F.O。;里维埃尔,B。;Frank,F.,《非精确层次尺度分离:基于非连续Galerkin离散的线性系统的双尺度方法》,计算。数学。申请。,74, 8, 1769-1778 (2017) ·Zbl 1398.65088号
[47] 海伦布鲁克,B。;Mavrilis,D。;Atkins,H.,连续和不连续有限元离散的“p”多重网格分析,(第16届AIAA计算流体动力学会议(2003),美国航空航天研究所)
[48] Tamstorf,R。;琼斯,T。;McCormick,S.F.,用于布料模拟的平滑聚合多重网格,ACM Trans。图表。,34, 6, 1-13 (2015)
[49] 蒂勒,C。;Araya-Polo,M。;Alpak,F.O。;Riviere,B.,用于孔隙尺度两相流加速高精度模拟的相场方法的分布式并行混合CPU-GGPGPU实现,(SPE储层模拟会议(2019),石油工程师学会)
[50] Thiele,C.,可组合p-层次解算器
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。