立式陈列柜风幕数值模拟的双流体模型研究.pdfVIP

立式陈列柜风幕数值模拟的双流体模型 The Two-Fluid Model of Air Curtain Numerical Simulation in a Vertical Display Case 1, 2 1 2 余克志 丁国良 陈天及 1 上海交通大学 制冷与低温工程研究所 上海200030 2 上海水产大学 制冷工程系 上海200090 Tel: 021 E-mail: kzyu@shfu.edu.cn 摘要 考虑到陈列柜风幕的流动特点 采用双流体模型对立式陈列柜风幕进行数值模拟 其中内外风幕采用传统的K 紊流模型 陈列柜外的环境空气采用非紊流模型 两种流体之间的质量 动量和能量交换采用经验关系式处理 通过对风幕 温度的测试结果比较发现 双流体模型的模拟结果相当吻合 是很有价值的的风幕计算模型 同时指出了双流体模型的不足 和需要改进的地方 关键词 陈列柜 风幕 双流体模型 紊流 风幕对陈列柜的保温性能有十分重大的影响 果符合得很好[9-14] 但目前没有应用于陈列柜风幕 国内外许多研究者对其进行过研究[1 — 8] 考虑到陈 的实例 从物理机理上分析 相对于紊流的单流体 列柜的开口处的边界条件难于确定 为计算方便起 模型 双流体模型更能够准确地反应紊流的流动状 见 很多研究者将开口处计算区域加以扩充 并认 况 但由于同时求解两套方程组 同时两种流体的 为扩充区域的边界不受陈列柜风幕的影响 其速度 方程之间存在强烈的耦和关系 其求解的难度更 按无滑移的边界条件考虑 温度按实际测试的结果 大 需要一定的求解技巧 确定 由于陈列柜外由于风幕卷吸的环境空气流速 2风幕的双流体模型 很小 不会象风幕一样处于强烈的紊流状态 但他 们的计算模型往往采用单一的紊流模型 K 模 2 1控制方程 型 雷诺应力模型或者大涡模拟 这显然与实际 情况不符 风幕双流体模型的通用方程是[13] f r 考虑到陈列柜外环境空气的非紊流特点 采用 (r r f ) (r r v f ) (r k ) (f D k ) S I t k k k x k k kj k x k fk x x k fk x fk fk 双流体模型对陈列柜风幕进行数值模拟 其中内外 j j j j j 风幕采用传统的K — 紊流模型 陈列柜外的环境空 ( ) ( ) ( ) ( ) ( ) ( ) 气采用非紊流模型 同时考虑到两种流体之间的质 其中下标k 表示流体种类 可取为 1 或2 为 量 动量和能量交换 方便起见 流体1 表示紊流流体 流体2 表示非紊 流流体 下标j 表示空间坐标 r 表示该种流体占 1双流体模型的思路和特点 有的体积分数 或存在几率 f 表示因变量 取 紊流的双流体模型是Spalding 提出的 其基本 为1 时对应连续性方程 如取速度 温度分别对应 思路是 认为紊流流动可以看作两种流体各自的运