亚网格内两相流动界面位置的精确定位算法

孙红霞 1; 2; 苏军伟 3

doi:10.7652/xjtuxb201701013

您当前的位置：

首页 >

文章列表页 >

亚网格内两相流动界面位置的精确定位算法

更新时间：2025-07-09

- 亚网格内两相流动界面位置的精确定位算法
- A Sub-Grid Interface Locating Algorithm for Free Surface Two-Phase Flows
- 西安交通大学学报 2017年51卷第1期页码：79-87
- 作者机构：
  
  1. 中国石油大学(华东)石油工程学院,山东,青岛,266580
  2. 中国石化股份胜利油田分公司勘探开发研究院,山东,东营,257015
- 作者简介：
- 基金信息：
  
  : 2016-05-09。作者简介: 孙红霞(1983—),女,副研究员;苏军伟(通信作者),男,讲师。基金项目: 国家科技重大专项资助项目(2016ZX05011003-007);国家自然科学基金资助项目(21306145)
- DOI：10.7652/xjtuxb201701013
  中图分类号： O359
- 网络首发：2017-01-10，
  
  纸质出版：2017
- 稿件说明：
移动端阅览
孙红霞 1, 2, 苏军伟 3. 亚网格内两相流动界面位置的精确定位算法[J]. 西安交通大学学报, 2017,51(1):79-87.

A Sub-Grid Interface Locating Algorithm for Free Surface Two-Phase Flows[J]. 2017, 51(1): 79-87.
孙红霞 1, 2, 苏军伟 3. 亚网格内两相流动界面位置的精确定位算法[J]. 西安交通大学学报, 2017,51(1):79-87. DOI： 10.7652/xjtuxb201701013.

A Sub-Grid Interface Locating Algorithm for Free Surface Two-Phase Flows[J]. 2017, 51(1): 79-87. DOI： 10.7652/xjtuxb201701013.

摘要

针对流体体积(VOF)法处理两相流动界面张力的虚拟流体方法无法准确确定亚网格内界面位置点的问题

提出了一种高精度的界面位置定位算法。该方法根据界面两侧单元界面曲度在亚网格内构建变曲度弧线

实现了两相界面的高精度近似

完成了网格单元中心连线和界面交点位置以及该位置的曲度的准确求取

成功地考虑了界面曲度对亚网格内界面位置的影响

同时方法本身不受网格类型及空间维度的限制。针对毛细管上升、静态液滴、二维气泡上升以及三维气泡上升等物理过程进行了模拟

结果表明:界面曲度对亚网格界面位置可产生显著影响; 亚网格内界面位置的准确捕捉以及该位置点界面曲度的准确求取可以显著改善VOF方法的预测精度

模拟结果显示最好可提高3.5倍。该方法可为复杂结构内高表面张力驱动流提供新的预测手段。

Abstract

In ghost fluid method

for treating surface tension under the framework of volume of fluid(VOF)

the exact location of fluid interface cannot be easily determined. In order to resolve such a problem

a new sub-grid interface locating algorithm is proposed in this work. In the new algorithm

a variable curvature curve

used to approximate the interface in sub-grid

is constructed based on the curvature of the two grid cells on both sides of the interface. As a result

the intersection point of the connection line of centers of interfacial cells and the interface is calculated accurately. The influence of the curvature on the interface location is taken into account. Moreover

this algorithm is independent on mesh type and space dimension. Numerical tests

including capillary rise

static droplet

bubble rise in two-dimensions and three-dimensions

show that interface curvature has a significant influence on the sub-grid interface location. Calculating the sub-grid interface location and the curvature at such a location with high accuracy can improve the numerical prediction of VOF method significantly

and by 3.5 times at the best. This method can provide a new tool for predicting physical behaviors of high surface tension driven flow in complex geometrical structures.

关键词

Keywords

references

UNVERDI S O, TRYGGVASON G. A front-tracking method for viscous, incompressible, multi-fluid flows [J]. Journal of Computational Physics, 1992, 100(1): 25-37.

SUSSMAN M, SMEREKA P, OSHER S. A Level Set approach for computing solutions to incompressible two-phase flow [J]. Journal of Computational Physics, 1994, 114(1): 146-159.

HIRTC W, NICHOLSB D. Volume of fluid(VOF)method for the dynamics of free boundaries [J]. Journal of Computational Physics, 1981, 39(1): 201-225.

COQUERELLE M, GLOCKNER S. A fourth-order accurate curvature computation in a level set framework of two-phase flows subjected to surface tension forces [J]. Journal of Computational Physics, 2016, 305(1): 838-876.

YOUNGS D L. An interface tracking method for a 3D Eulerian hydrodynamics code: technical report 44/92/35 [R]. Berks, UK: AWRE, 1984: 3-9.

EDWARD P J J, GERRY E P. Second-order accurate volume-of-fluid algorithms for tracking material interfaces [J]. Journal of Computational Physics, 2004, 199(2): 465-502.

SCARDOVELLI R, ZALESKI S. Interface reconstruction with least-square fit and split Eulerian-Lagrangian advection [J]. Int J Numer Meth Fluids, 2003, 41(3): 251-274.

SUSSMAN M, GERRY E P. A coupled level set and volume-of-fluid method for computing 3d and axisymmetric incompressible two-phase flows [J]. Journal of Computational Physics, 2000, 162(2): 301-337.

BALCÁZAR N, ORIOL L, JOFRE L, et al. A coupled volume-of-fluid/level-set method for simulation of two-phase flows on unstructured meshes [J]. Computers Fluids, 2016, 124(1): 12-29.

SCARDOVELLI R, ZALESKI R. Direct numerical simulation of free-surface and interfacial flow [J]. Annual Review of Fluid Mechanics, 1999, 31(1): 567-603.

DIOT S, FRANCOIS M M. An interface reconstruction method based on an analytical formula for 3D arbitrary convex cells [J]. Journal of Computational Physics, 2016, 305(15): 63-74.

LPEZ J, HERNÁNDEZ J, GMEZ P, et al. A new volume conservation enforcement method for PLIC reconstruction in general convex grids [J]. Journal of Computational Physics, 2016, 316(1): 338-359.

FRANCOIS M M, CUMMINS S J, DENDY E D, et al. A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework [J]. Journal of Computational Physics, 2006, 213(1): 141-173.

ISSA R I. Solution of the implicitly discretized fluid flow equations by operator-splitting [J]. Journal of Computational Physics, 1985, 62(1): 40-65.

SU J, Gu Z L, CHEN C G, et al. A two-layer mesh method for discrete element simulation of gas-particle systems with arbitrarily polyhedral mesh [J]. International Journal for Numerical Methods in Engineering, 2015, 103(10): 759-780.

OpenCFD Ltd. OpenFOAM [CP/OL]. [2016-03-20]. www.openfoam.com.

ADAMI S, HU X Y, ADAMS N A. A new surface-tension formulation for multi-phase SPH using a reproducing divergence approximation [J]. Journal of Computational Physics, 2010, 229(13): 5011-5021.

孙鹏楠, 李云波, 明付仁. 自由上浮气泡运动特性的光滑粒子流体动力学模拟 [J]. 物理学报, 2015, 64(17): 1747011-1747015.

SUN Pengnan, LI Yunbo, MING Furen. Numerical simulation on the motion characteristics of freely rising bubbles using smoothed particle hydrodynamics method [J]. Acta Phys Sin, 2015, 64(17): 1747011-1747015.

WALTERS J K, DAVIDSON J F. The initial motion of a gas bubble formed in an inviscid liquid: Part 1 The two-dimensional bubble [J]. Journal of Fluid Mechanics, 1962, 12(3): 408-416.

BHAGA D, WEBER M E. Bubbles in viscous liquids: shapes, wakes, and velocities [J]. Journal of Fluid Mechanics, 1981, 105(1): 61-85.

王宁宁,刘海湖,张楚华.液桥重开过程的两相格子Boltzmann模型及计算.2016,50(9):101-105.[doi:10.7652/xjtuxb201603016]

范小军,李亮,李森,等.平板和静叶表面气流-水膜耦合流动特性的数值研究.2016,50(11):101-105.[doi:10.7652/xjtuxb201603016]

朱晓静,毕勤成.垂直上升内螺纹管中高压汽-水两相流截面含汽率的测量.2015,49(3):61-66.[doi:10.7652/xjtuxb 201510010]

曹文瑾,孙中国,席光.系泊型浮体运动的无网格法数值模型.2015,49(3):65-71.[doi:10.7652/xjtuxb201501011]

施东晓,毕勤成,周荣启.磁性液体两相界面演变特性的数值研究.2014,48(9):112-116.[doi:10.7652/xjtuxb201409019]

李涛,宗潇,杨小平,等.矩形通道内高速蒸汽与过冷水直接接触凝结换热流型的实验研究.2014,48(5):128-135.[doi:10.7652/xjtuxb201408022]

王太,李会雄,李阳.同轴两个气泡融合特性的数值研究.2013,47(7):99-103.[doi:10.7652/xjtuxb201305018]

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

液氢管路输运中静电积聚机理与电流密度分布特征研究

X射线辐照下聚酰亚胺材料空间电荷分布特性的数值模拟

异质裂隙内超临界CO2地热提取特性数值模拟研究

相对浸润性对双元液滴碰撞影响的数值研究

工质与叶轮传热对超临界二氧化碳向心透平气动性能与流动特性的作用机制