流体横掠圆柱相变融解过程的数值模拟研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083211 流体横流体横掠圆柱相变融解过程的数值模拟掠圆柱相变融解过程的数值模拟 流体横流体横掠掠圆柱相变融解过程的数值模拟圆柱相变融解过程的数值模拟 1 2 王华军 ，戴传山 （河北工业大学能源与环境工程学院，天津 300401） (Tel:022 Email:huajunwang@126.com) 摘要：对层流范围内流体横掠圆柱相变储热装置的融解过程进行了数值模拟。结果表明，圆柱相变融 解主要由固体导热和复杂的自然对流过程组成，其中顶部融解速度较快，底部融解较慢。当融解至一 定程度后，圆柱底部会出现明显的 Rayleigh- énard 对流现象。R 对流是影响相变融解过程的重要机 制，且融解时间与来流流速之间存在非线性变化关系。当流速处于某一范围时，R 对流表现比较强烈， 融解时间相对较短一些。 关键词：相变储热，融解特性，自然对流，数值模拟 0 前言 相变储热在建筑节能、工业废热回收、电子冷却以及新能源利用等领域具有十分重 [1] 要的应用前景 。前人研究表明，固液相界面经常存在导热－自然对流耦合现象，其复 杂的非线性特征使得精确解析解变得困难，而更多依赖于数值模拟或实验研究[2－3] 。目 前的 PCM 自然对流研究主要偏向于两种情形：1）温度梯度与重力方向垂直 （简称“自 然对流”）；2 ）温度梯度与重力方向平行 （简称“R 对流”）。对于流体横掠圆柱相变储 热装置的融解情况，温度梯度与重力的夹角范围为 0º－180º，因而上述两种自然对流可 能同时存在且相互影响。本文将对此问题进行初步数值模拟，旨在为类似相变储热装置 的优化设计提供一定的参考依据。 1 数值模型 1.11.1控制方程控制方程 1.11.1控制方程控制方程 图 1 给出了流体横掠圆柱相变储热装置的几何模型，其中流体流动方向与重力相互 垂直。计算控制方程如下： ∇⋅u = 0 (1) ∂u 2 ρ +ρ u⋅∇ u=−∇p+∇ u− ρ −ρ g+ S (2) ( ) ( o ) ∂τ ∂t ∂G ρc +ρc u⋅∇t =∇ λ∇t −ρL (3) p p ( ) ( ) ∂τ ∂τ G为液相体积百分比，定义如下： 资助项目：国家自然科学基金项目(No.50576065