1、置换通风热力分层高度的数值研究同济大学 刘猛 龙惟定 张改景摘要：本文根据置换通风系统的工作原理，介 绍了置换通 风热力分层高度的定义及数值计算方法，并 对送风温度、送风速度等因素对热 力分层高度的影响作了定性的分析，在此基础上拟合出无量纲热力分层高度关于送风温度 T、送风速度 V、热源间距离 L 和围护结构传热 Q 的经验公式。本文的工作对今后置换通风系统的优化设计提供参考。关键词：置换通风 热力分层高度 数值模拟 热舒适性0 引言随着人们对可持续发展认识的不断深入，节能和环保已成为当今空调发展的两大主题。相对于传统的混合通风方式，置换通风具有较高的室内空气品质和较好的节能效果，因而备受人们
2、的关注 1。置换通风两大优点的发挥依赖于流场的分布，而热力分层高度是置换通风流场分布的标志，也是置换通风系统设计中的重要参数，不同影响因素对流场的影响将通过热力分层的高低反映出来，事实上，也只有掌握了分层高度随影响因素的变化关系，才可能设计自控线路，进而调节送入室内气体的相关参数，使节能和高的空气品质同步实现。因此，本文主要工作将围绕建立热力分层高度与流场影响因素之间定量关系来开展，重点导出热力分层高度与流场影响因素之间的定量关系。1 热力分层高度的定义假定置换通风系统通送风量为 QS，热源产生的自然对流射流量即热烟羽流量为 QP，在热源周围自然对流射流形成的初期阶段，热烟羽流量主要是靠置换通
3、风的送风量来补充，因而卷吸周围较少的空气量，此时 QPQS。根据连续性原理，在热烟羽的上升过程中，必有某高度 Z 处，此时自然对流射流吸卷周围空气的量与置换通风系统送风量相等，即 QP=QS，如下图1 所示：图 1 热力分层高度示意图刘猛，男，1979年9月生，在读博士研究生，地址：上海市曹安公路4800号同济大学嘉定校区13-115信箱，邮政编码：201804，电话：（021）69587882， Email：随着空气在下部区域的不断送入以及顶棚处热污浊空气的不断排出，室内热烟羽状态稳定下来，形成以 Z 高度处为分界的上、下两个区域。Z 高度以下为单向流动的清洁区，Z 高度以上为混合紊流的污浊
4、区，Z 高度即为热力分层高度。2 置换通风房间模型的建立在本文的研究中，采用 Srebric 的办公室置换通风测试作为算例 2，这也是 ASHRAE（美国供暖、制冷空调工程师协会）用于验证数值计算结果的报告中的算例之一，如下图 2 所示。图 2 置换通风房间模型该办公室的几何尺寸为：长 宽 高5.16m3.65m2.43m，置换通风送风口为扁平型置换通风散流器，送风口的中心线与墙体中心线相重合，距离地面 0.03m，送风温度是 21.2，送风速度是0.12m/s。回风口位于房顶顶板的中心，尺寸为 0.43m0.43m，日光灯尺寸为 0.2m0.2m0.8m，均匀分布在房间顶部。室内有计算机、人
5、员、灯光等热源，还有桌子和壁柜等家具，室内热源功率如表1 所示。表 1 室内热源功率热源项目 热源大小 （W）人体 752计算机 1 108计算机 2 173日光灯 344围护结构 180总计 7473 置换通风房间流场的计算置换通风系统的气流组织形式与其它通风系统的不同之处在于浮力的影响很明显，室内气流方式由送风口送入的冷空气流和室内热源上方产生的浮力烟羽所控制。本文数值模拟中将流体假定为不可压缩流体，使用 Boussinesq 近似 3来处理浮力的影响，采用涡团扩散湍流模型。S.J.Ress 等人的研究指出，置换通风系统中，当室内负荷在 4572W/m 2时，室内气流会产生类似周期性的流动
6、，热源上方的烟羽到达吊顶撞击顶棚，会有复杂的侧向摆动，产生明显回流 4。此房间热负荷为：747W/ (5.16m3.65m)39.6 W/m2，不超过 40 W/m2，根据 S.J.Ress 的理论，数值模拟可以产生稳定的流场。经过多次调试，得到如下所示的温度场和速度场，迭代计算的收敛过程大约需要8h。图 3 Z1.825 截面上的温度分布 图 4 Z1.825 截面上的速度分布图 3 和图 4 与文献2的模拟结果进行比较，发现二者基本一致，只是图中所显示的温度区间有所不同，但两者温度梯度是相符的。这说明本文关于流场的模拟正确，在此基础上可进行热力分层高度的数值研究。4 热力分层高度的数值计算
7、方法在计算收敛的流场中取图 5 所示的模拟测点，分别取 A、B、C、D、E、F 六个不同位置，由于 A 点在室内流场中的位置较有代表性，A 点处的温度分布能很好的反映室内的流场分布（计算表明，B 点、C 点、E 点三处的平均温度和 A 点温度值相差甚微），因此，选取 A 点数据进行计算热力分层高度的计算。将 A 点在不同高度的温度值在温度高度坐标系中标出，再拟合通过这些点的曲线，由温跃层理论可知 5，曲线拐点所对应的高度 Z 就是置换通风房间的热力分层高度。由图 6 中曲线看出，在 0.75m1.25m 的范围内，空气温度发生了较大的变化，因此可以判断热力分层高度就在此高度区间范围内。曲线前半
8、部分向上凸后半部分向下凹，由拐点的性质可知，在此高度区间内存在曲线的拐点。图 5 数值模拟测点的位置分布 图 6 A 点温度的拟合曲线为求出拐点的坐标值，将 对 x 求二阶导数 ，并令60()myxa()yx，解方程得 ，而 ，此值即为该置换通风房间的热力()0yx18.29*K*()1.98分层高度。5 不同工况下计算结果的分析本文就不同送风温度、送风速度、热源间距离和不同围护传热 14 种工况下的热力分层高度分别进行了计算，其数值如表 2 所示。表 2 不同工况下的热力分层高度值工况 无量纲热力分层高度 z 送风温度 T（） 送风速度V（m/s） 热源间距离 L（m） 围护结构传 热 P（
9、W）1 0.404115 21.2 0.08 4.2 1802 0.422634 21.2 0.10 4.2 1803 0.451852 21.2 0.12 4.2 1804 0.474074 21.2 0.14 4.2 1805 0.494239 21.2 0.16 4.2 1806 0.383951 20.2 0.12 4.2 1807 0.479012 22.2 0.12 4.2 1808 0.513169 23.2 0.12 4.2 1809 0.616461 21.2 0.12 0.4 18010 0.558436 21.2 0.12 0.8 18011 0.483128 21.2
10、0.12 1.6 18012 0.44321 21.2 0.12 4.2 013 0.446914 21.2 0.12 4.2 9014 0.457613 21.2 0.12 4.2 240对比表 2 的各组数据，可以得到如下的结论：（1）送风温度增大，热力分层高度增高，工作区域空气温度与送风温度差减小。因此，增大送风温度，有利于热力分层高度的增高，可避免因送风温度低对人体产生不舒适的吹风感。（2）室内热负荷、送风温度等其它条件相同的情况下,随着送风速度的增大即送风量的增大,热力分层高度也在增高；反之降低。这是因为送风速度增大时，加强了上升热气流的卷吸作用，使其扩散加快,促使热气流与周围空气掺
11、混。（3）热力分层高度随着热源之间距离的增大而减小，分散的热源不利于热力分层高度的提高。（4）围护结构在不同的传热功率下，房间下部区域的高度温度曲线基本上是重合的，这说明其热力分层高度基本上是一致的，热力分层高度值都在 1.1m 左右，这说明围护结构热损失的大小,不影响置换通风的分层特性,不影响热力分层的高度,而仅仅影响了室内空气的温度值。6 无量纲热力分层高度经验公式工程设计过程中需要快速准确地确定热力分层的高度值，因此必须建立送风速度、送风温度及热源间距离等参数关于热力分层高度的关联式。在设计参数确定的情况下，通过关联式即可快速估算出热力分层的高度，为系统的优化设计提供参考。考虑到设计过程
12、中房间会有不同的高度，为使关联式在不同高度房间设计中适用，建立的关联式须采用无量纲热力分层高度，即 z=Z/H，此处 H 为房间的高度。假设所要建立的关联式形式为： 3124nnnnzaTVaLQ将表 2 中的数据代入上式中将得到 14 个方程的方程组，该方程组中含有a1、n 1，a 4、n 4八个未知数，该方程组方程的个数超过了未知数的个数，因此该方程组属于超定方程组，而超定方程组在极特殊的情况下才有精确解,一般来说超定方程组在普通意义下是无解的，只能在新设定的准则下定义它的解。因此，为求解的方便，先给定n1、n 2、n 3、n 4的值如表 3 所示，使方程组转化为线性超定方程组，用最小二乘
13、法求解矛盾方程组 6得到 a1、a 2、a 3、a 4的值，如表 4 所示，这样就建立了无量纲热力分层高度 z 关于送风温度 T、送风速度 V、热源间距离 L 及围护结构传热 Q 的关联方程式。表 3 不同 n 的取值n1 n2 n3 n40.2 0.08 0.25 0.1表 4 系数 a 的取值a1 a2 a3 a40.07564 0.75556 0.23206 0.00354于是得无量纲热力分层高度 z 的表达式为： 08.2.075664. VTz 1.025.06QL对于不同高度的房间，为得出实际的热力分层高度，用关联式算出的无量纲热力分层高度乘以设计房间的高度 H，即 Z=zH 就是
14、设计房间的热力分层高度值。本文将将经验公式和数值方法计算得到的热力分层高度值对比如表 5 所示。表 5 经验公式计算值和数值计算值的对比工况 1 2 3 4 5 6 7数值计算 0.982 1.027 1.098 1.152 1.201 0.933 1.164经验公式 1.0457 1.0728 1.0953 1.1145 1.1314 1.0185 1.0984相对误差（） 6.492 4.464 0.244 3.254 5.798 9.169 5.631工况 8 9 10 11 12 13 14数值计算 1.247 1.498 1.357 1.174 1.077 1.086 1.112经验
15、公式 1.1014 1.4541 1.3693 1.2684 1.0809 1.0944 1.09574相对误差（） 11.670 2.929 0.905 8.039 0.358 0.769 1.462从上表看出，除工况 8 的相对误差是 11.670以外，其余工况的相对误差都在 10以下，有些工况的相对误差不到 1。因为在热力分层高度 Z 处，存在一个厚度为 的过渡层，该层是由于对流紊流和热力扩散的平衡作用而形成的，该空气层极不稳定，易受各种外界因素而扰动。对于典型办公室一般 为 0.15m0.25m 厚，这说明热力分层高度 Z 本身就是不稳定的，会有一定的浮动，因此采用该经验公式可以达到工
16、程设计的精度要求，在工程设计中可以直接用此公式估算热力分层的高度。将文献7中的实验结果与该公式的计算结果列于表 6 进行对比，文献7中使用的是三热源模型，且三个热源间距离相当，都在 0.25m 左右，故经验公式中取热源间距离 L 为0.25m。工况送风温度T（）送风速度V（m/s）热源间距离L（m）围护结构传热P（W）实验结果（m）经验公式计算值（m）相对误差（）1 20 0.15 0.25 0 0.9188 0.9342 1.67612 21 0.15 0.25 0 0.9200 0.9362 1.76093 22 0.15 0.25 0 0.9388 0.9381 0.7456表 6 经验
17、公式计算值与实验结果的对比从表 6 看出实验值与经验公式计算所得的热力分层高度吻合得较好，由此可见，如果工况的参数在样本区间之内，用此经验公式计算所得的 Z 与实验值吻合得较好，可靠性较高；如果工况的参数落在样本区间之外，用经验公式计算所得的 Z 值与实际热力分层高度值可能会有较大的偏差，其解决的办法是扩大样本的容量，例如送风温度 T 取值区间扩大至 1826，其它送风参数取值区间也相应扩大，用数值方法计算出样本扩大后各种工况下的热力分层高度，在此基础上再拟合经验公式，就能保证经验公式有更大的适用范围。因此，本文的工作有待于进一步深入。7 结束语本文结合置换通风系统的工作原理，介绍了热力分层高
18、度的定义，对所建模型的置换通风房间进行了流场的数值模拟计算与分析，在此基础上对置换通风房间的热力分层高度进行了数值计算，并分析其影响因素，如：送风温度、送风速度、围护结构传热和热源分布等，最后拟合了无量纲热力分层高度相应于送风温度 T、送风速度 V、热源间距离 L 和围护结构传热 Q 的经验公式，为设计者提供有价值的设计参数和设计依据。参 考 文 献1 李龙宇,李强民.置换通风的原理及应用.通风除尘,1996(1)2 Srebric JSimplified methodology for indoor environment design.Ph.d.thesis.USA:Massachuset
19、ts Institute of Technology,20003 陶文铨.数值传热学，西安交通大学出版社,20014 S.J.Rees, J.J.McGuirk, P.Haves. Numerical investigation of transient buoyant flow in a room with a displacement ventilation and chilled ceiling system，International Journal of Heat and Mass Transfer 44（2001）306730805 倪波.置换通风的实验研究J. 暖通空调,2000,30(5):246 李信真,车刚明,欧阳洁等.计算方法.西安:西北工业大学出版社,2000 （2）7 陈俊俊,王晓彤,武文斐.三污染热源置换通风热力分层高度实验研究,内蒙古科技与经济,2004（2）4 23 0.15 0.25 0 0.9588 0.9401 1.9504
