×

具有纹理平衡孔隙的材料的水平集方法。 (英语) Zbl 1349.76816号

摘要:结构平衡控制着许多天然多孔材料中液相的分布,例如部分熔融的岩石和合金、盐水和冰水系统。在这些材料中,孔隙几何形状的演变使固液界面能最小化,同时在固液接触线处保持恒定的二面角(θ)。我们提出了一种水平集方法来计算在三维多个固体颗粒的空间填充镶嵌下处于结构平衡的液固界面的隐式表示。每个晶粒由一个单独的水平集函数表示,通过将表面扩散下的固液界面演化为恒定平均曲率表面,实现界面能最小化。使用虚拟对流和法向速度项将液体体积和二面角约束添加到公式中。这导致了非线性耦合PDE系统的初值问题,该系统使用固体颗粒的隐式表示作为初始条件,控制每个颗粒水平集的演化。设计了一种区域分解方案,将每个晶粒的计算域限制在晶粒周围的几个网格点。界面之间的耦合是在原始计算域的更高层次上实现的。通过高阶空间微分方案和通过区域合成实现计算本地化,提高了离散化的空间分辨率。还获得了不同晶粒分布网络的三维解示例,说明了该方法的几何灵活性。

MSC公司:

76S05号 多孔介质中的流动;过滤;渗流
76米25 其他数值方法(流体力学)(MSC2010)
80平方米 其他数值方法(热力学)(MSC2010)
80A10号 经典热力学和相对论热力学
80A22型 Stefan问题、相位变化等。
2008年6月 地球物理学相关问题的计算方法

软件:

工具箱LS
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] von Bargen,N。;Waff,H.S.,《部分熔融体系的渗透率、界面面积和曲率:平衡微观结构的数值计算结果》,J.Geophys。决议,91,B9,9261-9276(1986)
[2] 雷佩尔(A.W.Rempel)。;E.D.沃丁顿。;Wettlaufer,J.S。;Worster,M.G.,《通过预融化和异常扩散可能取代古冰中的气候信号》,《自然》,411,6837,568-571(2001)
[3] 刘易斯,S。;Holness,M.,平衡石盐-(H_2O)二面角:浅层地壳中的高岩盐渗透性?,地质学,24,5,431-434(1996)
[4] Smith,C.S.,《晶粒、相和界面:微观结构的解释》,Trans。金属。美国医学会学会,175,15-51(1948)
[5] 甘巴扎德,S。;普罗达诺维奇,M。;Hesse,M.A.,《各向异性结构平衡孔隙网络中的渗流和晶界润湿》,Phys。修订稿。,113, 4, 048001 (2014)
[6] 孔B,M。;福斯特,硕士。;约翰·D·菲茨·杰拉尔德(John D.Fitz Gerald),《解读熔融岩石和固化岩石中的二面角》,《虚拟探险家杂志》,第35、2页(2010年)
[7] Smith,C.S.,《多晶微观结构的一些基本原理》,《金属》。修订版,9、1、1-48(1964年)
[8] de Gennes,P.G.,《润湿:静力学和动力学》,修订版。物理。,57, 3, 827-863 (1985)
[9] Bulau,J.R。;瓦夫,H.S。;Tyburczy,J.A.,《部分熔体中流体分布的机械和热力学约束》,J.Geophys。Res.,84,6102-6108(1979)
[10] van der Marck,S.C.,多孔介质中非零传输阈值的证据,水资源。决议,35,2,595-599(1999)
[11] 朱伟。;Gaetani,G.A。;Fusseis,F。;Montesi,L.G.J。;Carlo,F.D.,部分熔融岩石的显微断层扫描:地幔橄榄岩中的三维熔体分布,科学,3326025,88-91(2011)
[12] Schoenherr,J。;Urai,J.L。;Kukla,P.A。;Littke,R。;施莱德,Z。;拉罗克,J.-M。;M.J.纽瓦尔。;Al-Abry,N。;Al-Siyabi,H.A。;Rawahi,Z.,《岩盐封闭能力的限制:来自南阿曼盐盆地的卡布里亚-阿拉盐案例研究》,Am.Assoc.Pet。地质。公牛。,91, 11, 1541-1557 (2007)
[13] Sobolev,A.V。;Shimizu,N.,超亏损熔体和海洋地幔渗透性,横贯。多克。俄罗斯科学院。科学。地球科学。第节。,328182(1993年)
[14] Lundstrom,C.C。;吉尔·J。;威廉姆斯,Q。;Perfit,M.R.,地幔熔融和通过平衡多孔流提取玄武岩,《科学》,2701958-1961(1995)
[15] 施皮格曼,M。;Elliott,T.,《年轻熔岩中铀系不平衡的熔体输送后果》,《地球行星》。科学。莱特。,118, 1-20 (1993)
[16] 霍尔内斯,M.B。;Lewis,S.,根据岩盐-(H_2O-CO_2)系统平衡二面角的压力和温度变化推断的岩盐-海水界面结构,Geochim。科斯莫辛。《学报》,61,4795-804(1997)
[17] Downey,M.W.,《油气聚集密封性评估》,美国协会宠物。地质。公牛。,68, 11, 1752-1763 (1984)
[18] Beere,W.,《渗透液相和烧结孔稳定性的统一理论》,《金属学报》。,23, 1, 131-138 (1975)
[19] Takei,Y.,《孔隙几何形状对VP/VS的影响:从平衡几何形状到裂纹》,J.Geophys。固体地球研究,107,B2,ECV61-ECV612(2002)
[20] 佩尔武基纳,M。;Kuwahara,Y.,固液复合材料的电性能和弹性性能与界面能量控制的平衡微观结构之间的相关性,《地球行星》。科学。莱特。,265, 3-4, 410-422 (2008)
[21] Wimert,J。;Hier-Mayumder,S.,地球深部内部熔体几何形状的三维微观地球动力学模型,J.Geophys。固体地球研究,117,B4,B04203(2012)
[22] Oprea,J.,《肥皂电影的数学:与Maple的探索》,学生数学图书馆,第10卷(2000年),Amer。数学学会·Zbl 0972.53003号
[23] Nye,J.F.,《多晶冰中水脉和节点的几何形状》,J.Glaciol。,35, 119, 17-22 (1989)
[24] Cheadle,M.J.,《织构平衡两相聚集体的特性》(1989年),剑桥大学博士。
[25] Hirth,G。;Kohlstedt,D.L.,《部分熔融上地幔动力学的实验约束:扩散蠕变状态下的变形》,J.Geophys。固体地球研究,100,B21981-2001(1995)
[26] Endres,A.L。;Murray,T。;A.D.布斯。;West,L.J.,《利用雷达和地震波速联合估算地层含水量和孔隙几何形状的新框架》,《地球物理》。Res.Lett.公司。,36、4、L04501(2009)
[27] Osher,S。;Sethian,J.A.,曲率相关速度的前沿传播:基于Hamilton-Jacobi公式的算法,J.Comput。物理。,79, 1, 12-49 (1988) ·Zbl 0659.65132号
[28] Osher,S.J。;Fedkiw,R.P.,《水平集方法和动态隐式曲面》(2002),施普林格出版社
[29] Sethian,J.A.,《水平集方法和快速行进方法》(1999),剑桥大学出版社:剑桥大学出版社·Zbl 0929.65066号
[30] 彭,D。;梅里曼,B。;Osher,S。;赵,H。;Kang,M.,一种基于PDE的快速局部水平集方法,J.Compute。物理。,155, 2, 410-438 (1999) ·Zbl 0964.76069号
[31] Smereka,P.,曲率和表面扩散运动的半隐式水平集方法,科学杂志。计算。,19, 1-3, 439-456 (2003) ·Zbl 1035.65098号
[32] Spelt,P.D.M.,模拟具有多个滞后移动接触线的流动的水平集方法,J.Compute。物理。,207, 2, 389-404 (2005) ·Zbl 1213.76127号
[33] 刘,H。;Krishnan,S。;Marella,S。;Udaykumar,H.S.,夏普界面笛卡尔网格法II:模拟液滴与任意形状表面相互作用的技术,J.Compute。物理。,210, 1, 32-54 (2005) ·Zbl 1154.76358号
[34] Li,Z。;Lai,M.-C。;He,G。;Zhao,H.,移动接触线自由边界问题的一种增广方法,计算。流体,39,6,1033-1040(2010)·Zbl 1242.76047号
[35] Jettestuen,E。;Helland,J.O。;Prodanovic,M.,用非零接触角在孔隙尺度上模拟毛细控制位移的水平集方法,Water Resour。研究,49,8,4645-4661(2013)
[36] 普罗达诺维奇,M。;Bryant,S.L.,《确定排水和渗吸临界曲率的水平集方法》,《胶体界面科学杂志》。,304, 2, 442-458 (2006)
[37] Lee,W。;儿子,G。;Jeong,J.J.,气泡生长和微腔脱离的数值分析,数值。热传输。,B部分,Fundam。,58323-342(2010年)
[38] 肖普,D.L。;Sethian,J.A.,曲率内禀拉普拉斯运动,界面自由边界。,1, 1-18 (1999) ·Zbl 0938.65144号
[39] Mitchell,I.M.,《水平集方法的灵活、可扩展和高效工具箱》,J.Sci。计算。,35, 2-3, 300-329 (2008) ·Zbl 1203.65295号
[40] 洛伦森,W.E。;Cline,H.E.,《行进立方体:高分辨率三维曲面构造算法》(SIGGRAPH 1987)。SIGGRAPH 1987,美国加利福尼亚州阿纳海姆,第21卷,共4卷(1987)),163-169
[41] Armstrong,R.T。;波特,M.L。;Wildenschild,D.,《将孔隙尺度界面曲率与柱状毛细压力联系起来》,《高级水资源》。,46, 55-62 (2012)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。