第3卷第4期
2004 年 12 月
热科学与技术
Jo urnal of Thermal Science and Tec hno lo gy
Vol.3No.4
Dec
.
2004
文章编号: 1671-8097(2004)04-0353-05
收稿日期: 2003-04-25; 修回日 期: 20 04-10-02 .
作者简介: 李海军(1971-), 女, 博 士生.
蒸汽喷射器内二维流场的数值模拟研究
李海军, 沈胜强, 张 博
( 大连理工大学 动力工程系, 辽宁 大连 116024 )
摘要:
通过求解二维 N -S 方程,对以水蒸气为工质的喷射器内复杂的流场进行数值模拟。分析比较了四种不
同的湍流模型,C hen-K im 修正的
k
-ε模型用于数值模拟,分析了工作参数及喷射器结构对喷射器内部的流场
及出口激波的形成的影响,得出了喷射器设计的优化方案。
关键词: 喷射器; 数值模拟; 激波
中图分类号:TB617 文献标识码
:A
0
引 言
喷 射制冷是一种绿色环保的制冷方式。它不
直接消耗机械能,可以利用太阳能、工业废汽、余
热等低品位能源;同时,其在工质选择上比较灵
活,可以避免使用氟利昂等造成温室效应的气体。
因此,随着人类对能源需求的飞速增长,以及对环
保意识的提高,喷射制冷作为一项开源节能的技
术受到越来越多的关注。
目 前,限制喷射制冷技术推广应用的主要原
因是喷射器的效率过低,从而影响整个制冷系统
的效率。以往对喷射器的研究主要是建立在一维
等熵基础上的计算模拟;其重点是分析工作参数
对喷射系数的影响,而不考虑其内部具体流动的
情况,因此,在指导喷射器结构设计上有一定局
限
[1, 2]
。本文借助 PH O E N IC S 软件平台,采用多种
湍流模型,通过求解二维 N -S 方程,模拟喷射器内
的流场。比较了结构、工作参数变化对喷射器效率
的影响,分析说明了壅塞现象出现的条件,为喷射
器的优化设计提供了依据。
1
喷射过程计算的控制方程
喷 射器内的流动具有强烈的湍动性,受计算
机计算能力及流动本身复杂性的限制,直接求解
三维 N-S 方程是很难实现的。本文利用在
B oussinesq 假设基础上的有效黏度概念的湍流模
型来进行数值求解。分别采用了 LV E L 零方程模
型、标准
k
-ε湍流模型、Chen-K im
k
-ε湍流模型和
RNG
k
-ε模型来模拟喷射器内的流动工况,计算
喷射器的性能参数,并将结果与实验值进行比较。
下面对各模型进行简要的介绍。
1
.
1
LVEL
模型
LVEL 湍流模型是零方程模型; 它是
P H O E N IC S 特有的,于 1993 年加入软件包的。它
引 入如下定义的无量纲距离
y
+
、无量纲速度
u
+
和无量纲有效黏度:
y
+
=
y
τ/ρ/υ
l
(1)
u
+
=
u
/τ/ρ (2)
υ
+
=d
y
+
/d
u
+
=υ/υ
l
(3)
式中:
y
为距最近壁的距离,
u
为平行于壁面的速
度,ν
l
为层流黏度。
y
+
与
u
+
的关系由 Spalding 壁
面定律给出,则流动中每一点的有效运动黏度 ν
+
可以由 N ew ton-R aphson 迭代计算得出。
1
.
2
标准
k
-ε两方程模型
k
-ε两方程模型是 H arlow 和 N akayam a 1968
年提出的,是目前应用最广泛的双方程涡黏度湍
流模型。这里使用的是 Launder 和 Spalding 1974
年提出的适于高雷诺数的
k
-ε模型:
∂
(ρ
k
)
∂
t
+
∂
ρ
u
k
-
ρ
v
t
Pr
k
·
∂
k
∂
x
i
∂
x
i
=ρ(
P
k
+
G
b
-ε)
(4)