第2卷第4期
2003 年 12 月
热科学与技术
Jo urnal of Thermal Science and Tec hno lo gy
V ol.2,N o.4
Dec
.
2003
文章编号: 1671-8097(2003)04-0318-06
收稿日期: 2003-06-24; 修回日 期: 20 03-11-07 .
基金项目: 留学归国人员科研启动基 金资助项目(2002).
作者简介: 赵大勇(1978-), 男, 硕 士生.
VOF
方法中几种界面重构技术的比较
赵大勇, 李维仲
( 大连理工大学 动力工程系, 辽宁 大连 116024 )
摘要:
通过理论分析和数值模拟,比较了界面跟踪技术中 V O F 方法的三种界面重构技术。利用三种流场,即
平移场、旋转场和剪切场对它们进行了数值模拟和比较。结果表明,除有些尖角会被少许抹平外,PLIC 方法无
论在界面的精细度上还是在锐利性上都要好于前两者,而 FLA IR 方法次之。
关键词: V O F; 界面重构; 施主 -受主; FL A IR ; PL IC
中图分类号:O359 文献标识码
:A
0
引 言
具有运动界面的流体流动现象广泛存在于机
械、能源、电子、动力、水利和建筑等多种工程领域
中。运动界面的追踪是了解这类现象的关键,特别
是在多相、多组分的气体、液体和固体的内部速
度、压力、密度等的数值模拟中起着非常重要的作
用。在理论上这种流动也是流体力学中一个具有
挑战性的研究工作,因此长期以来受到各国学者
的重视。曾提出各种不同的追踪技术对不同问题
进行了较为广泛的研究,取得了较丰富的成果。归
结起来界面追踪技术可以分为三类:Euler 法、
L agrange 法和 E uler-Lagrange 混合法。发展比较
成熟且用得较多的是 E uler 法。在 E uler 法中,根
据 网格变化情况,又有三种方法,即固定网格法、
自适应网格法和保角变换法。其中固定网格法最
直接并且计算量小,因而被广泛采用。
在 固定网格法中,自由面形状求解的方法分
为自由界面追踪(surface tracking) 技术和流量追
踪(volum e tracking) 技术。自由面追踪技术用一
系列插值点和曲线段来表示自由面形状,然后按
流动发展方程和自由面演化方程计算出每一个时
间步内,这些离散点和曲线线段的新位置,从而构
造出新的自由面形状。它具有方法简单、存储量小
等优点,但是在处理自由面翻卷和合并时,对出现
的多值问题显得无能为力而造成逻辑上的混
淆
[1]
。流量追踪是根据每个网格单元的流量流动
情况来构造自由面形状。在流量追踪中的
V O F(volum e of fluid) 方法以其存储小、应用简
单、追踪的界面锐利性好、精细,而深受广大 C FD
工作者的喜爱。而且随之发展起来的各种界面重
构技术在 V O F 中越来越占据主流地位。从最初的
施主 - 受主格式到完全的重构技术。界面变得越
来越精细,精度越来越高。为了选择较好的重构技
术,本文研究和比较了流行的 H irt 和 N ichols、
F LA IR 和 PL IC 三种重构技术。利用速度反向还
原方法比较这三种方法的精度和界面的锐利性,
得到了明显的结果;这为今后的研究和 V O F 法的
使用奠定了基础。
1
VOF
方法的基本原理
VOF方法是Hirt和Nichols
[2]
在 1981 年首先
提出的;他通过定义了一个流体体积函数
F
,使
F
的值等于一个单元内流体体积与该单元体积之
比。如果
F
= 1,说明该单元内充满流体;如果
F
= 0,说明该单元内不含流体;如果 0 <
F
<1,说
明该单元内含有自由界面。
F
的守恒型传输方程
(二维) 如下:
∂
F
∂
t
+
∂
(
uF
)
∂
x
+
∂
(
vF
)
∂
y
=0 (1)
式中:
u
、
v
分别为主场沿
x
和
y
方向的速度;
F
为
步进函数,用通常的差分格式进行离散求解会抹
平
F
的间断性,从而失去函数的原有定义,所以必
须采用特殊的方法进行计算。由于
F
值沿自由面