VOF方法中自由液面重构的数值模拟.pdfVIP

科技情报开发与经济 SCI～TECHINFORMATIONDEVELOPMENTECONOMY 2011年 第21卷 第5期 文章编号：1005—6033(2011)05—0196-02 收稿 日期：0211—O1—13 VOF方法中自由液面重构的数值模拟 范石伟 ，陈义胜 ，庞赘佶 ，罗剑慧 ，蒋大伟 (1．内蒙古科技大学材料与冶金学院，内蒙古包头，014010；2内蒙古科技大学能源与环境学院，内蒙古包头，014010) 摘 要：介绍了VOF方法中的斜线法，详细写出了重构的具体过程，并与经典的计算方 法——HIRT—Nichols方法做 了对比，发现斜线法易于编程，结果更接近真实的液面情 况 。 关键词 ：自由液面；VOF方法；自由液面重构；斜线法 中图分类号：TV131．4 文献标识码 ：A VOF(流体体积函数)法是固定网格法中追踪 自由表面的一 道每个单元格的F值，就可以编制程序迭代求解出直线与每个 种方法。自由面的重构是VOF方法中研究的一个重要部分，即利 单元格的交点，即可模拟出自由表面的形状，并能求出液面与单 用单元格的流体体积计算 自由面的位置和形状。F为 目标流体的 元格交点的具体位置。 体积与网格体积的比值。只要知道每个网格的，值，再计算出自 1．2 斜线重构法中流量的计算 由表面的方向，就可以实现对 自由界面的追踪 ¨]。流量追踪是根 法向可以定出自由面的大致方向，自由面可分为下述8种 据每个网格单元的流量来构造 自由面形状。流量追踪中的VOF 情况 ： =O， O；以 ：O， O； O， =O； O，J--O；n, (volumeoffluid)方法存储小、应用简单、模拟结果接近真实值，而 0，以J0；nJO，以√0； J0， JO；n~jO，峨oo。 且随之发展起来 的各种界面重构技术在 VOF中越来越占据主流 在上述 8种情况中，每一种情况根据速度不同又对应下面 8 地位。VOF自由面追踪方法包括两部分研究内容：一是根据给定 种情况 ：UO，=O；UO，=0；U=O，O；U=O，o；UO，o； 的F值确定 自由液面位置及方向；二是根据重构好的自由液面 ， 0，V0；UO，VO；UO，VO。 计算时间 内每个单元格的流量并确定下一时刻F值。 图1为自由面的法向图，图2为网格液体体积流量图。 事实已经证明不能采用F方程的离散形式来计算流体的体 积流量，因为结果会产生重叠现象，发生数值耗散 [2]。要根据具 体情况，结合速度来计算单元格液体的体积流量。为了选择较好 的重构技术 ，本文详细介绍了斜线法的具体过程，并与其 HIRT-Nichols【方法进行了比较。发现斜线法 [s精度更高，方法 也更简单。 1 VOF方法中的斜线重构法 图1 自由面的法向图 1．1 斜线重构法中法向的确定 确定自由面要有两个因素：一是 自由面的法向n( ， )；二 是每个网格上 自由面的交点。但是当网格F=I和F=0时，可认为 固 流向(，+1，J+1) 此网格的不存在法向，即n(o，O)。因为当网格F=I和F=O时已 暖 流向(，，^1) ×D