高压高过冷度下圆管流动沸腾流动传热特性的数值模拟.docxVIP

高压高过冷度下圆管流动沸腾流动传热特性的数值模拟 高压水过冷思想实验研究 在石油、化工、能源、能源等领域，过度冷、流动和沸腾作为一种非常重要的传热方式，得到了广泛应用。其中石油化工管式加热炉炉管内介质的过冷沸腾流动, 较为常见的工况是管内介质的压力较高且初始温度较低, 目前针对此领域内过冷沸腾传热的研究相对较少。所以深入研究高压高过冷度下沸腾及气液两相流的传热特性, 进一步了解沸腾传热机理, 对于提高系统的传热性能及优化系统设计都具有十分重要的意义。 Bartolomej等较早对垂直管中高压水过冷沸腾进行了实验研究。Kurul等采用双流体模型分析了高压过冷沸腾中气液两相流局部参数分布特性。随后较多研究者对过冷沸腾进行了数值模拟及实验的研究, 多是针对一些低压低过冷度工况。通过相关理论及实验的研究, 逐步认识到气泡直径作为过冷流动沸腾模型的关键参数, 包括气泡脱离直径及主流区气泡直径。气泡直径大小不同, 两相相界面浓度必然发生变化, 形成两相间输运特性的差异, 也就对模型的预测精度起到显著影响。孙奇等对低压沸腾管内气液两相流气泡尺寸的径向分布特性进行了测量与研究。Tu分析了低压过冷沸腾流动中气泡尺寸模型对于气泡相分布的影响。上述研究揭示了更适合实验工况的相间传输模型, 对于过冷沸腾模型的完善起到了一定发展, 然而这些模型多是适合一些低压工况的模型, 且在理论上都未深入分析气泡直径对壁面传热方式的影响。本文借助Fluent流体计算平台, 首先采用高压下的不同气泡直径模型并结合相关实验数据进行了对比分析计算, 揭示了气泡直径对壁面传热方式的影响并确定所应用模型。在此基础上, 对高压高过冷度工况下水的过冷流动沸腾进行了数值模拟, 得到气相体积分数及表面传热系数的分布规律, 并着重分析了压力及壁面热通量对高压高过冷度工况流动及传热特性的影响。 1 两个流模型和过度冷却模型 1.1 量、动能及能量守恒方程 双流体模型假定气泡相与流体相为相互贯穿的连续性介质, 建立两相各自的控制方程, 两相之间的相互作用通过控制方程中的相间传输项来描述。每一相的基本控制方程包括质量、动量及能量守恒方程, 具有如下的形式: 质量守恒方程 动量守恒方程 能量守恒方程 式中下角标k和j为相标记, 取L或g时分别代表液相或气相;α、ρ、v、H、μe、k′和S分别表示体积分数、密度、速度矢量、焓、有效黏性系数、有效热导率及源项;Γjk、Fkj和Qkj分别为相间的质量、动量和能量传输项;Γjkvk和ΓjkHk分别表示由于相间的质量传递引起的相间动量传递及能量传递。为使方程能够封闭, 需要补充相间传输项的关系式, 以描述过冷沸腾过程中气液两相之间质量、动量及能量的交换。 1.2 过度冷和沸腾的模型 1.2.1 汽化核心区域的热传热 在过冷流动沸腾过程中, 气泡产生以前, 固体壁面完全被液体覆盖, 传热机理为单相对流传热。随着气泡不断在汽化核心位置处形成、长大并脱离壁面, 周围的冷液体补充到气泡所占据的区域, 瞬间被加热升温, 存在蒸发传热和激冷传热两种不同的机理。所以固体壁面上的热通量分为3部分:蒸发热通量Qe, 激冷热通量Qq和液体与壁面间的单相对流热通量Qc。热通量的计算式分别为 壁面总热通量为 1.2.2 共同常数的传递 气液两相间的动量传递可以通过相间作用力描述, 主要考虑曳力Fd、升力FL、湍流耗散力FTD及壁面润滑力FW。 相间总作用力 1.2.3 蒸发传质速率ts 相间质量传递过程来自于两个区域:加热壁面附近液体的汽化及主流区域液体的汽化或冷凝。相间质量传递速率计算式为 式 (9) 右端第一项为壁面附近液体的蒸发传质速率, 第二项为主流区的传质速率。主流区两相之间的质量传递速率依赖于液体的温度, 当液体过冷时, 传质方向为从气相到液相, 气体冷凝。当液体过热时, 传质方向为从液相到气相, 液体汽化。式中, hLg为相间传热系数, 根据Ranz-Marshall关联式获得;Ts为液体饱和温度, 作为当地压力的函数;ALg为界面面积浓度 1.2.4 壁面加热对热速率的计算 气液两相间的能量传递为经过相界面面积ALg的传热速率, 计算式为 气相温度设置为饱和温度, 假定在相变过程中保持不变, 壁面加热引起的源项只在液体相能量方程中考虑。液体相能量方程的源项为 1.2.5 气泡游离模型 正确预测气泡脱离直径dbW非常重要, 关系到壁面蒸发传热速率的大小。 Tolubinsky等的气泡脱离直径模型关系式为 式 (12) 的气泡脱离直径模型源自于高压水沸腾实验数据, 只依赖于液体的过冷度。而Unal的模型假定气泡在过热的微液层内蒸发产生, 同时气泡顶部受冷凝的影响。在加热壁面气泡可以滑动或破碎, 当生长到最大尺寸时脱离壁面。Unal模型中气泡脱离直径作为压力、液体过