基于VOF方法的带相变的自由界面的计算.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号：066037 基于VOF方法的带相变的自由界面的计算 袁明豪杨燕华李天舒胡志华 上海交通大学核科学与工程学院200030 电话：021—mail=mhvuan(situ．edu．cn 摘要：对于气液自由界面上存在相变(燕发、沸腾或凝结)的情况，由于此时不仅物性分布是不连续 的，速度场也发生了不连续，这无论对于捕捉或追踪自由界面的算法本身还是流场的求解部带来了一 定的困难。本文提出一种基于VOF方法的计算带相变的自由界面的算法．对这种情况下的关键问题一 一速度场的不连续性给出了解决方法。并用编制的二维程序计算了水平壁面上的膜态沸腾，并将计算 结果与理论公式进行了比较。 关键词：自由界面；v0F：相变 1．引言 迄今为止，已经有各种基于连续介质假设的数值算法可用于计算带有气液自由界面 的问题[1】。由于界面附近流体的物性存在剧烈的变化，所以需要精确地捕捉或追踪自 由界面的位置，从而将此物性和压力(存在表面张力的情况)的不连续性反映到流场的 求解中去。各种自由界面算法，都通过各自的方式做到了这一点。例如，MAc方法通过示 of 踪粒子来追踪界面。VOF(VolumeFIUid)法通过流体体积函数，LevelSet方法通 过Level Set函数来表示自由界面的位置，一些动网格算法使网格线始终与自由界面的 位置相一致从而来准确处理界面上的运动学及动力学边界条件。但对于自由界面上存在 相变(蒸发或凝结)的情况。由于此时速度场也发生了不连续，这无论对于捕捉或追踪 自由界面的算法本身还是流场的求解都带来了一定的困难。对此类问题求解的尝试可见 文献[2】【3】。 本文提出一种基于VOF方法的计算带相变的自由界面的算法，对这种情况下的关键问 题——速度场的不连续性给出了解决方法。并用编制的二维程序计算了水平壁面上的膜 态沸腾，并将计算结果与理论公式进行了比较。 基金项目：国家自然科学基金资助项目 501 2物理模型 把液体和蒸汽视为不可压缩牛顿流体．动量方程为 1 ’ 兰㈣+Ⅳ啊。刁=却+厦Ⅳ桐+㈣rl+20嘶， c| 其中U是速度，p是密度，∥是动力粘度，P是压力，罗是体积力(重力)，口是表面 张力，。是当地的界面曲率．J，是集中在界面5上的狄拉克函数，刀是界面上的法向 矢量。能量方程如下； 蚂(詈+fi．vT卜(kvr)m 2 其中省略了粘性耗散。这里厂是温度。C是定压比热，k是导热系数。中是界面上相变 涉及到的潜热产生的能量源项。连续性方程如下： V-i=F． 3 其中．r是在两种密度不同的流体界面上发生传质而引起的体积的膨胀和收缩。以上方 程中的物性在每一相中保持常数．但相与相之间是不同的。 为了确定计算区域中的某个点是属于液相还是气相，引入一个特征函数x，在某一 相中取为1，而在另一相中取为0，描述z的运动方程为 VF方法中的体积份额函数C可以被看作特征函数x的离散形式。 当界面上存在蒸发或凝结时，r和通过自由界面的质量流量口之间的关系如下： 妒2孵一去声， 5 其中V是某个区域，S是包禽在V内的自由界面．如图1所示。P，和pt是相1和相 2的密度。在导出式5的过程中，使用了相1和相2的不可压缩性。 图1包含自由