×

两相流的三维移动网格有限元方法。 (英语) 兹比尔1349.76166

总结:发展了一种三维ALE有限元方法,用于研究两相流现象,使用一种新的离散化方法计算表面张力。该计算方法基于任意拉格朗日-欧拉公式(ALE)和有限元方法(FEM),创建了液气界面改进模型的两相方法。由于计算网格节点根据流移动,因此还提出了一种自适应网格更新程序,以有效地管理网格以删除、添加和修复元素。ALE描述通过一组相互连接的节点明确定义了两相界面位置,确保了边界的清晰表示,包括表面张力的作用。所提出的计算曲率的方法能够以适度的编程工作量和计算成本获得准确的结果。为了验证该方法,进行了静态和动态测试,结果与文献中的分析解和实验结果进行了很好的比较,表明新提出的方法能够很好地描述界面力和气泡动力学。本文重点描述了所提出的方法,特别强调了表面张力的离散化、新的网格重划分技术和验证结果。此外,还对两个细长气泡进行了复杂几何形状的微通道模拟。

MSC公司:

76M10个 有限元方法在流体力学问题中的应用
65平方米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
76T10型 液气两相流,气泡流
82C26型 统计力学中的动态和非平衡相变(一般)
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 胡,H。;北帕坦卡。;Zhu,M.,使用任意拉格朗日-欧拉技术对流体-固体系统进行直接数值模拟,J.Compute。物理。,169, 427-462 (2001) ·Zbl 1047.76571号
[2] 希尔特,C。;Nichols,B.,自由边界动力学的流体体积(VOF)方法,J.Compute。物理。,39, 201-225 (1981) ·Zbl 0462.76020号
[3] 萨斯曼,M。;斯梅雷卡,P。;Osher,S.,《计算不可压缩两相流解的水平集方法》,J.Compute。物理。,114, 146-159 (1994) ·Zbl 0808.76077号
[4] Bourlioux,A.,用于追踪材料界面的耦合水平集流体体积算法,(第六届计算物理国际研讨会论文集,加州塔霍湖(1995))
[5] 萨斯曼,M。;Puckett,E.,用于计算三维和轴对称不可压缩两相流的耦合水平集和流体体积法,J.Compute。物理。,162, 301-337 (2000) ·Zbl 0977.76071号
[6] 哈洛,F.H。;Welch,J.E.,含自由表面流体随时间变化的粘性不可压缩流动的数值计算,物理。流体,82182-2189(1965)·Zbl 1180.76043号
[7] Amsden,A。;Harlow,F.H.,《SMAC方法:计算不可压缩流体流动的数值技术》(1970年),技术代表,洛斯·阿拉莫斯·Zbl 0206.55002号
[8] Glimm,J。;格罗夫,J。;Lindquist,W。;O.麦克布莱恩。;Tryggvadson,G.,跟踪标量波的分岔,SIAM J.Compute。,9, 1, 61-79 (1988) ·Zbl 0636.65132号
[9] 佩罗,B。;Nallapati,R.,用于模拟不可压缩自由表面流动的移动非结构化交错网格方法,J.Compute。物理。,184, 1, 192-214 (2002) ·Zbl 1118.76307号
[10] Quan,S。;Schmidt,D.,三维不可压缩两相流的移动网格界面跟踪方法,J.Compute。物理。,221, 761-780 (2006) ·Zbl 1216.76050号
[11] 戈伊斯,J。;Nakano,A。;诺纳托,L。;Buscaglia,G.,移动最小二乘曲面的前方跟踪,J.Comput。物理。,227, 9643-9669 (2008) ·Zbl 1153.65017号
[12] Braess,H。;Wriggers,P.,自由表面流动的任意拉格朗日-欧拉有限元分析,计算。方法应用。机械。工程,190,95-109(2000)·Zbl 0967.76053号
[13] Ganesan,S。;Tobiska,L.,计算含可溶性表面活性剂两相流的任意拉格朗日-欧拉有限元法,J.Compute。物理。,231, 3685-3702 (2012) ·Zbl 1402.76072号
[14] Babuska,I.,有限元方法的误差界限,数值。数学。,16, 322-333 (1971) ·Zbl 0214.42001号
[15] Brezzi,F.,《拉格朗日乘子鞍点问题的存在性、唯一性和逼近》,RAIRO,Ana。编号。,R2,129-151(1974)·Zbl 0338.90047号
[16] Souza,F。;Mangiavacchi,N.,《模拟不可压缩双流体流动的拉格朗日水平集方法》,《国际数值杂志》。方法流体,471393-1401(2004)·Zbl 1155.76360号
[17] 德斯布鲁姆,M。;梅耶,M。;施罗德,P。;Barr,A.,使用扩散和曲率流对任意网格进行隐式光顺,(Siggraph论文集(1999)),317-324
[18] Geuzaine,C。;Remacle,J.-F.,Gmsh:一个具有内置预处理和后处理设施的三维有限元网格生成器,国际期刊Numer。方法工程,79,11,1309-1331(2009)·Zbl 1176.74181号
[19] Si,H.,通过约束Delaunay精化的自适应四面体网格生成,Int.J.Numer。方法工程,46,7,856-880(2008)·Zbl 1195.65129号
[20] Gordon,W。;Thiel,L.,超有限映射及其在网格生成中的应用,应用。数学。计算。,10-11, 171-233 (1982) ·Zbl 0495.65053号
[21] Anjos,G.,《含相变两相流的三维ALE有限元方法》(2012年),埃科尔理工学院-法国洛桑分校博士论文
[22] Chorin,A.J.,Navier-Stokes方程的数值解,数学。计算。,22, 745-762 (1968) ·Zbl 0198.50103号
[23] Tomiyama,A.,与计算气泡动力学的斗争,(第三届多相流会议论文集-ICMF(1998))
[24] 华,J。;Stene,J。;Lin,P.,使用前沿跟踪方法对粘性液体中三维气泡上升的数值模拟,J.Compute。物理。,227, 6, 3358-3382 (2008) ·Zbl 1329.76262号
[25] 杜,J。;固定,B。;Glimm,J。;贾,X。;李,X。;李毅。;Wu,L.,一个简单的前端跟踪软件包,J.Compute。物理。,213, 613-628 (2006) ·兹比尔1089.65128
[26] 巴加,D。;韦伯,M.,《粘性液体中的气泡:形状、尾迹和速度》,《流体力学杂志》。,105, 61-85 (1981)
[27] Hua,J。;Lou,J.,粘性液体中气泡上升的数值模拟,J.Compute。物理。,222, 769-795 (2007) ·Zbl 1158.76404号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。