喷口射流流场的数值模拟 转载

喷口射流流场的数值模拟
摘要 本文在对数学模型及物理模型进行理论分析的基础上，基于ε-k 方程双方程模型与SIMPLE 算法，运用模拟软件6.1，对不同喷口流场进行了三维模拟，证实了射流流场与喷口内部流线有密切关系，并对结果进行了分析比较，从而为喷口的优化设计提供了依据。

关键词 喷口；计算流体力学；数值模拟 1 引言
目前，高大空间的空调一般采用分层空调的方式。

该系统一般采用喷口送风，喷口的结构，位置对室内气流组织影响较大。

国内的生产厂家较少涉足该喷口的制作，有的厂家尽管生产，一般都套用国外产品的性能参数.由于制造工艺等方面的原因，其性能参数与国外的产品不同，按此设计其空调效果较差。

实践证明，一些进口的喷口，空调效果比较理想。

从物理角度而言，空调通风房间的空气流动是由风口射流引起的，射流动量流量和质量流量对室内空气分布情况起关键作用。

为了能准确的描述风口，产生了很多风口模型的模拟方法。

就目前的研究情况来说，比较广泛使用的风口模型模拟方法[1]
主要包括传统模型法、基本模型法、动量方法、盒子方法和指定速度方法。

通常这些模型对末端装置的特性不予考虑或考虑得不够，模拟主要是针对室内的气流运动而进行的，而没有模拟送风装置内部的气流[2]。

本文通过计算机直接模拟喷口形成入口边界条件对不同结构的喷口进行了模拟。

2 数学物理模型 2.1 风口模型
计算机直接模拟喷口形成入口边界条件是先从管道到喷口出口处的三维区域进行模拟，得出出口断面处的各变量参数，然后将此结果作为对室内空气流动模拟的入口边界条件。

2.2 基本方程
［4］
Fluent6.1 把连续方程、动量方程、能量方程写成如式2.1所示的通用方程形式：
()φφφφρρφS x u x t i i i
=⎪⎪⎭⎫
⎝
⎛∂∂Γ-∂∂
+∂∂ （2.1） 式中：φ _____通用变量，如u ,v ,w ,T ,,K ε等；
ρ _____流体的密度， kg/m 3；
i u _____速度矢量，m/s ；
φΓ _____通用变量φ的有效热交换系数； φS _____变量的源项。

如空调室内的气流属于常温、常压及低速下的流动，可以视为不可压缩流动，但空气的温差产生的自然对流对室内的流场和温度场存在着显著的影响。

对这类流动现象的描述引入著名的Boussinesq 假设，可以大大简化所讨论的问题。

Boussinesq 假设为：流体密度的变化并不显著地改变流体的性质，即流体的其它物性不变：密度的变化对惯性力项、压力差项和粘性力项的影响可忽略不计，仅考虑密度的变化对质量力的影响。

考虑质量力时，流体的运动方程为：
()()⎪⎪⎭
⎫
⎝⎛∂∂∂∂+∂∂-=∂∂
+
∂∂j i
j
i i j i i i x
x x p f u u x t u μμρρρ （2.2） 当考虑密度的变化仅是由温度的变化引起的，并且质量力仅是重力的情况时，由热膨胀系数β的定义式：⎪⎭
⎫
⎝⎛∂∂-=∂∂=
T T V V ρρβ11（2.3） 由（2.3）式可得：
)(00T T --=-ρβρρ (2.4)
()()i i i i i i g T T g g g g f 00000ρβρρρρρρ+--=+-== （2.5）
上式中()00T T g i --βρ即为流体温差导致的所谓Boussinesq 假设，代入运动方程（2.2）中得：
()()()⎪⎪⎭
⎫
⎝⎛∂∂∂∂+∂∂-+--=∂∂
+
∂∂j i
j
i i i j i i i x
x x p g T T g u u x t u μμρβρρρ000 （2.6） 表2.1 广义扩散系数与源项取值
表2.2 ε-K 模型常数取值
在直角坐标系下，φ，φΓ，φS 的值如表2.1所示。

其中u ,v ,w 分别表示x,y,z 三个方向的速度（m/s ）；T ,ρ,0T ,0ρ，,K ε，P 分别表示送风温度(℃)及密度（kg/ m 3
），室内气流温度(℃)及密
度（m 3
/h ），紊流脉动动能及紊流脉动动能耗散率和压力（P a ）；x g ,y g ,z g 分别表示x,y,z 三个方向的
重力加速度（m/s 2
）；l μ，t μ，eff μ为层流粘性系数，紊流粘性系数和有效粘性系数；T σ，K σ，ε
σ分别为T ,,K ε的当量普朗特数。

基于以上模型及湍流模型，对式2.1所示的通用微分方程采用有限容积法进行离散，数值差分格式采用幂函数格式，求解法为动量方程在交错网格进行求解的SIMIPLE 计算法。

下面采用以上的数学物理模型及数值计算法对几种喷口流场进行数值模拟。

3 数值模拟研究 3.1研究对象
选取四个不同形状的喷口如图3.1，研究其喷口内部流线对射流流场的气体动力特性参数的影响。

()
(
)
2
/12
222
22
2
1
1/31/3111-⎪⎭
⎪
⎬⎫⎪⎩⎪⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡+-⎪⎪⎭
⎫
⎝⎛--=l x l x r
r r r （3.1）
式中0r 、1r 是喷口进口和出口截面的半径，l 为喷口进口到喷口出口的距离，x 为壁面上任一点至喷口进口的轴向距离。

喷口d 为某公司生产的喷口。

3.2初始条件
五个喷口的计算条件完全相同，外流场的计算域为20d ×5d ，计算域的入口采用速度入口边界条件，入口 速度设为10m/s ，出口采用出流条件，给定环境压力；固壁处采用无滑移边界条件，计算域及网格的划分如图3.2所示。

a 喷口a 的部分网格
b 喷口b 的部分网格
c 喷口c 的部分网格
d 喷口d 的部分网格
图3.2喷口a 、b 、c 、d 的部分网格示意图
3.3结果分析
从图3.3,3.4喷口速度矢量图看（颜色深或浅表示速度大或小），四个喷口的轴心速度都较高。

喷口b 较喷口a 轴心速度大，说明喷口b 较喷口a 对气流的阻滞作用大。

喷口c 与喷口d 速度分布基本一
直一致。

从速度颜色深度可以看出喷口a、b速度变化从接近出口突然开始，而喷口c、d速度变化从喷口内部沿流道曲线较缓变化。

从图3.5,3.6喷口速度矢量图看，对于收缩喷口a、b在出口附近的速度场差异较大，若喷口内部设计为维多辛斯基及样品喷口就可以获得更好的流场参数，有利于提高喷口的工作效率。

从图3.7喷口轴心速度分布看，喷a、b轴线上的速度最大值基本相同，喷口c、d的速度最大值基本相同。

喷口a、b较喷口c、d的速度最大值大，但喷口a、b轴心速度衰减较喷口c、d快，喷口d
衰减最慢。

从该图上可以看出喷口内部气流速度变化情况：喷口a、b内部轴线上的速度经过速度增大不变减小增大减小；喷口c、d轴线上的速度则经过增大不变
增大减小。

喷口a速度最大值在（5，2.3）处，喷口b在（7，2.2），喷口c、d在（2，3.7），由此可以看出在相同条件下，锥角越大，对射流的阻滞效果越明显，内部流道为维多辛斯基曲线的喷口c及样品喷口d对射流的阻滞最小。

通过图3.8可以看出，出口截面r/r0在0.8以下时，喷口c、d无因次速度u/u0基本一致，均小于喷口a、b的速度且喷口b的速度最大。

出口截面在0.8~1.0时，速度由大至小的喷口依次为喷口d、c、a、b。

从图3.9,3.10中可以看出，喷口a、b出口截面静压变化大，总压衰减快。

内部流道为维多辛斯基曲线的喷口c及样品喷口d其出口截面产生的横向压力梯度及分速度逐渐减小，从而得到较均匀的速度场，样品喷口d较喷口c流场更均匀。

(m/s)
图3.3喷口a、b的速度分布图
(m/s)
图3.4喷口c、d的速度分布
(m/s)
图3.5喷口a、b的速度矢量图
(m/s) 图3.6 喷口c、d的速度矢量图
u
/
u
x/d
u
/
u
r/r
图3.7 喷口轴线的无因次速度分布图3.8 喷口出口截面的无因次速度分布P
S
(
1
5
P
a
)
r/r
图 3.9 喷口出口截面的静压分布图3.10 喷口出口截面的总压分布
4 结语
CFD模拟越来越广泛地应用于解决暖通空调领域的空气流动和换热问题。

合理地给定入口边界条件对于得到正确的模拟结果是很重要的。

本文通过计算机模拟喷口形成入流边界条件对喷口流场数值模拟是可行的。

这种方法入口速度的设置特征尺度都与实际的物理模型保持一致，保证了射流扩散和衰减的正确性，较好地描述了射流过程。

而且它还有其它的优势，不需要测量风速，入口速度分布是根据风量通过模拟计算出的，根据模拟结果可得出风口的特性。

样本喷口与标准喷口在分析的几方面，特性基本一致，而实际工程运用中对样本喷口反映良好，制造商可以根据这一特性研究喷口的制造模具，制造出适合实际工程的喷口，因此这种对喷口流场数值模拟的方法对开发新型产品是重要且必要的。

参考文献：
[1] 赵彬，李先庭，彦启森．室内空气流动数值模拟的风口模型综述．暖通空调2000，30（5）：33～
36
[2] 杜国付，端木琳，舒海文．工位空调送风气流数值模拟风口模拟比较．热科学与技术，2003，2
（2）：162～167
[3] 何枫，谢峻石，杨京龙．喷嘴内部流道型线对射流流场的影响．应用力学学报．2001，18（4）：115～
119
[4] 梅启元．空调房间气流组织的数值计算与模拟．南京：南京理工大学硕士学位论文，2002，１：26～
28。