求解两相蒸发和冷凝问题的气液相变模型.pdf

第 卷第期 年月 西安交通大学学报 收稿日期：作者简介：孙东亮（），男，博士，副教授；徐进良（通信作者），男，教授基金项目：国家自然科学基金青年基金资助项目（）；国家重点基础研究发展计划资助项目（）；北京市自然科学基金资助项目（）；河北省自然科学基金资助项目（）网络出版时间：网络出版地址：求解两相蒸发和冷凝问题的气液相变模型孙东亮，徐进良，王丽（华北电力大学新能源电力系统国家重点实验室，北京；华北电力大学低品位能源多相流与传热北京市重点实验室，北京）摘要：为了更精确简单地求解两相流中的蒸发和冷凝问题，基于 中的流体体积（）方法提出了一种气液相变模型，该模型适应于两相中一相为非饱和相，即处于过热或过冷状态，另一相为饱和相，即处于饱和状态 该气液相变模型中：非饱和相导热系数和比热容为真实的物性参数；假设饱和相导热系数等于，饱和相比热容等于非饱和相比热容；界面处相变率仅由非饱和相决定 最后，通过一维 问题相界面位置的分析解和二维膜态沸腾气相体积比的精确解，验证了该相变模型的精确性和可行性 文中工作为该相变模型的推广应用奠定了基础关键词：两相流；气液相变模型；流体体积方法中图分类号：文献标志码：文章编号：（），（，；，）：（），；，：；模拟气液两相流动与传热的数值方法主要分为颗粒轨道模型、双流体模型和运动界面追踪技术，其中运动界面追踪技术最能准确、全面地反映气液两相间的信息目前，是最广泛使用的数值模拟软件，它采用流体体积（）方法追踪相间的运动界面，但该方法仍缺乏准确的气液相变模型，无法准确模拟相间的蒸发和冷凝现象在已有的相变模型中，提出的模型得到了 ：最广泛的应用，其具体表达式为蒸发过程 ，（）冷凝过程 ，（）式中：和分别为气相和液相的相变率，（）；为控制相变强度的因子，该变量的设定依据为使包含相界面的计算网格的温度值维持在饱和温度，针对不同问题该变量具有不同的取值 等采用该模型对换热管内的沸腾过程进行了模拟，设为 ；等对热管内的蒸发冷凝现象进行了计算，也设为 ；文献 中均设为 在蒸发和冷凝问题中，相界面温度通常设定为饱和温度 根据傅里叶定律，可计算出相界面两侧热流密度的界跃值式中：为相界面单位法线矢量，方向指向气体侧文献 根据式（）推导出了以下相变模型（）（）（）由于推导过程的假设不合理，导致在部分问题中该模型的计算结果与实际物理过程具有较大的偏差例如在一维 问题中，界面位置随时间的变化规律与液体的导热系数无关，具体见其分析解表达式（），而在式（）中却包含了的信息，导致模拟结果与具有相关性目前，许多作者基于 方法开发了求解相变问题的程序代码，这些方法的核心内容就是如何精确求解相界面两侧的热流密度，即式（）中的和 虽然这些方法可以精确地模拟蒸发和冷凝问题，但求解过程较为复杂，需要获取相界面位置和相界面附近网格的信息，很难在 中推广本文基于 中的 方法，提出了一种求解两相蒸发和冷凝问题的精确简单的气液相变模型，该模型适应于两相中一相为非饱和相，即处于过热或过冷状态，另一相为饱和相，即处于饱和状态 本文还通过两个经典的相变问题验证了该相变模型的精确性和可行性，表明该模型可应用于微通道内的蒸发和冷凝问题控制方程 方法是 等在 年首先提出的，它引入了流体体积分数的概念，各相流体体积分数之和等于，即 （）式中：下标、分别代表饱和相、非饱和相饱和相与非饱和相体积分数的控制方程为（）（）（）（）式中：动量方程为（）（）（）（）其中 （）（）表面张力为表面力，通过 （）模型 可将表面力转化为体积力 （）（）式中：为表面张力系数；界面曲率 （）（）能量方程为（其中 （），（），（）（）（）饱和相为气相 饱和相为液相（）式中：为汽化潜热相变模型的构建该相变模型主要包括以下部分（）非饱和相导热系数 和比热容，为真实的物性参数（）假设饱和相导热系数，比热容，西安交通大学学报 第 卷 ：，饱和相各点温度相等，不存在导热，因此饱和相导热系数的大小对结果没有影响 本模型从数值计算角度出发，强制设定，从而可以强制保证饱和相内温度保持不变 同时，饱和相温度保持不变，的大小对计算结果不产生影响，因此为了计算方便，本模型设定，（）求解界面处相变率 该模型针对一相为非饱和相、另一相为饱和相的情况，故式（）可写成 （）（）式中：为相界面上热流密度界跃值，；指向饱和相侧与将表面张力由表面力转化为体积力的 模型相似，将表面上的转化为相界面两侧有限宽度内单位体积的热量值，即（）式中：表示相界面两侧有限宽度值，；表示单位体积的热量值，假设相界面附近非饱和相温度呈线性分布，则可构建出 （）（）以下基于图证明式（）的合理性 根据式（），得出在网格点处，（，）（）（）（）其中 （），（）（）（）将式（）、（）代入式（），得出，（）（）同理，可得出在网格点处，（）（）将式（）、（）代入式（）的右侧，得出（，）（）（）图相界面、饱和相、非饱和相及温度分布示意图通过式（）（），即可验证式（）的合理性最后，基于式（），可以得出界面处的相变率 （）（）相变模型的应用及验证将本文提出的气液相变模型应用于一维 问题和二维膜态沸腾问题，以验证模型的精确性和可行性，其中采用 的 方法捕捉相界面，采 用 算 法 对 速 度 和 压 力 进 行 耦 合求解 一维 问题一维 问题是带有相变的两相流动的经典问题，如图所示 当给气体侧的固体壁面上施加一个较高的温度边界条件时，气体侧就会产生温度梯度，热量向液体部分传递，由于液体部分已经处于饱和状态，因此传递到液体部分的热量会使得液体部分开始气化，在相界面处产生相变 该问题相界面位置（）的分析解 为（）（，）（）式中：是一个超越方程的解，该超越方程为 （）（），（）（）图一维 问题示意图该问题中非饱和相为过热气体，（），（）；饱和相为饱和液体，根据第期孙东亮，等：求解两相蒸发和冷凝问题的气液相变模型 ：本文提出的模型，假设，（）；，壁面温度与饱和温度的差值 通过网格独立性考核，选取宽度为 的网格进行计算，时间步长为 图给出了一维 问题中不同密度比率下气液界面位置随时间的变化曲线，从图中可以看出，本文相变模型的模拟结果与分析解相一致，验证了模型的精确性和可行性图相界面位置随时间的变化 二维膜态沸腾图给出了二维膜态沸腾示意图 该问题中非饱和相为过热气体，动力黏度 ，（），（）；饱和相为饱和液体，根据本文提出的模型，假设，（）；，图中为 不稳定性波长（）（）图二维膜态沸腾示意图由于流动结构的对称性，本文选取一半区域作为计算区域，宽度为，高度为 边界条件：底部壁面速度为无滑移边界条件，温度为定壁温边界条件，壁面温度比饱和温度高，即 ；左右均为对称边界条件；上部为压力出口边界条件 初始条件：速度初始条件为；液体侧初始温度为饱和温度，气体侧初始温度分布为沿着气液界面到底部壁面呈线性分布；初始气液界面位置为 （）（）（）通过网格独立性考核，计算区域网格数选为 ，同时相界面附近进行局部加密，最小网格为初始网格的，时间步长为 图显示了初始气泡形成过程中不同时刻的气（）（）（）（）（）（）图不同时刻相界面分布图西安交通大学学报 第 卷 ：液界面分布 图给出了该过程中气相体积比（当前时刻的气相体积初始时刻的气相体积）随时间变化的曲线图气相体积比随时间的变化图中精确解为 等基于自己开发的程序代码得出的计算结果，从图中可以看出，本文相变模型的模拟结果与精确解完全吻合，进一步验证了模型的精确性和可行性结论本文基于 中的 方法，提出了一种求解两相蒸发和冷凝问题的精确简单的气液相变模型，并通过一维 问题和二维膜态沸腾问题验证了该模型的精确性和可行性 该气液相变模型主要针对两相中一相为非饱和相、另一相为饱和相的情况，主要包括以下部分：（）非饱和相导热系数 和比热容，为真实的物性参数；（）假设饱和相导热系数，比热容，；（）界面处相变率 （）在该气液相变模型的基础上，作者正在开发可以应用于两相同时处于非饱和状态的相变模型，以进一步扩大模型的应用范围参考文献：，：，（）：，（）：，（）：，（）：毛文彬微流控系统中混合及汽液相变传热的数值模拟 广州：中国科学院广州能源研究所，（）：，：，（）：，（），（）：，（）：，：，（编辑荆树蓉）第期孙东亮，等：求解两相蒸发和冷凝问题的气液相变模型