第三类边界条件数值模拟室内热环境及其实验验证【工程建筑论文】.doc

工程建筑论文-第三类边界条件数值模拟室内热环境及其实验验证贾玉凤 邹志军 黄晨 罗行 李俊红摘要： 本文应用数值模拟软件，利用第三类边界条件对某实验房的室内热环境进行数值模拟，并通过实验进行验证验证结果表明模拟值与实际测量值基本吻合在数值模拟验证的基础上，论文通过设置不同围护结构热工特性、室外空气温度、以及送风参数的模拟，得到了相应室内热环境随围护结构热工特性、室外温度、送风参数变化的特性与规律，进一步扩大了实验范围，充实了实验手段 关键词： 数值模拟 实验验证 变参数模拟 0.引言 随着计算机的大容量化和高速度化以及计算流体力学的发展，在室内热环境方面，特别是大空间建筑室内热环境设计中已逐渐普及采用 CFD 来解决室内气流组织、热环境等问题的研究[1]，从而使室内热环境特性研究及其全面评价成为可能本文应用软件 Airpak，利用第三类边界条件对某实验室室内热环境进行数值模拟，并通过实验予以验证，进而利用数值模拟对室内热环境特性进行分析1.环境实验室简介 如图 1 所示，环境实验室内尺寸为 4.9m×3.5m×2.5m，墙体均采用保温材料气流组织采用顶送下回，送风口尺寸为 16cm×69cm，距东墙中侧设有一30cm×30cm 的回风口。

室内东西墙附近各有一个散热器，图 1 中 Z 向为北向2.数值模拟计算与结果 2.1 物理模型及数学模拟概况模拟用物理模型如图 1 所示，其墙体传热系数为 0.383W/(K×m2)，墙外侧温度28℃送风速度为 2.35m/s，送风温度 17.8℃，靠近东、西墙处的散热器散热量分别为 840W、2410 W，且室内日光灯关闭数值模拟用数学模型为 K-ε 紊流模型，利用第三类边界条件对房间进行热环境模拟对送、回风口及回风管处、散热器等采用了网格加密的处理，总网格数18655 个2.2 数值模拟结果2.2.1 温度场分布如图 2（a）、（b）所示，沿着风口自上而下，温度逐渐变化近风口处等温线密集，温度分布存在明显的扩散现象在图 2（a）中，由于右侧存在一个散热器，导致了两边温度分布并不对称在图 2（b）中，水平方向温度梯度明显变小，存在衰减现象，回风口处等温线相对稀疏，房间居住域温度变化相对缓慢图 2（c）为南墙表面的温度分布，从图中可以看出，墙面自下而上温度逐渐升高，离风口较远处的温度相对较高，等温线较密集图中所标数字单位均为℃2.2.2 速度场的分布图 3 为室内速度场模拟结果模拟结果表明，射流断面速度从射流中心开始逐渐向边界衰减并沿射程有所变化，导致流量沿程增加，射流直径略有增大。

回风口的气流近似于流体力学中所述的汇流离开汇点距离越大，流速衰减越大，呈二次方衰减[2]从图中可以看出，风口下方速度较大，自上而下存在衰减现象其余区域速度较小图 3（a）中，气流在左右两侧各形成一个较小的涡流图 3（b）中，除送风口与回风口处速度较大，整个房间的速度较小，且分布比较均匀3.实验验证 3.1 实验布点与测量方法实验中共布置九个速度测点，在宽度方向上取中间截面布置七个点，两个散热器附近各布置一个测点空气速度采用万向风速仪，其输出信号通过 Fluke 采集器进行集中采集布点位置如图 4（a）所示采用垂直方向上均匀布点的原则，实验中布置二十个温度测点，采用带防辐射屏蔽罩的 T 型热电偶进行测试，数据采集通过 Anjelun 采集器集中采集，每分钟采集一次，布点位置如图 4（b）所示3.2 实验结果与模拟值的对比分析表 1、表 2 分别为图 4（a）、图 4（b）各测点实验值定义系列测定误差为：其中 xs——实测值；xm——模拟值；n——测点总数计算 σ 时剔除最大偏差值经计算，速度系列误差 σv=0.15m/s，温度系列误差 σt=1.66℃速度误差相对较大，这是由于在速度均匀区域测点较少，某些点实测值与模拟值相差较大造成的。

温度误差相对较小对比表中的各个数值，说明模拟热环境与实际热环境基本一致，数值模拟结果可靠表 1 速度模拟值与实测值比较 测点序号1 2 3 4 5 6 7 8 9实测值(m/s)2.131.731.690.110.170.09 0.70.181.16模拟值(m/s)2.041.711.350.140.110.160.520.120.08表 2 温度的模拟值与实测值 测点序号1 2 3 4 5 6 7 8 9 10实测值（℃）25.225.9925.2618.4920.220.8220.5724.3823.9423.75模拟值（℃）22.9622.7622.4618.6919.5520.1120.1623.8623.4823.2测点序号11 12 13 14 15 16 17 18 19 20实测值（℃）23.7522.8924.32 2531.4225.1524.1223.0322.9623.35模拟值（℃）22.8222.1424.6923.5522.6722.4522.3723.7423.226.914.室内环境特性模拟 对围护结构传热系数、室外空气温度、以及送风温度、速度等参数进行了变参数模拟。

选取在房间中间位置点 10 和靠近出风口处点 5 作为观察对象参看图 4（b））4.1 变送风参数模拟结果分别设置送风温度 14、16、17.8、20、22℃，由图 5（a）可知，随着送风温度的增加，点 10，5 的温度都在增加，室内温度也随之升高，点 5 温度增加的趋势要高于点 10分别设置送风速度 1.8、2.1、2.35、2.7、3.0m/s，由图 5（b）可知，随着送风速度的增加，点 5 的速度增加趋势略为明显，这是与点 5 位于风口附近，受送风速度影响较大有关点 10 的风速变化并不明显，速度较均匀4.2 变热工参数的模拟结果分别设置墙体的传热系数为 0.383、2.5、4.5、6.5、8.5W/(K×m2)，其它参数不变，由图 6 可知，随着传热系数的增加，室内温度略有升高这是因为随着传热系数的增加，材料的保温性能降低，比较容易受到室外参数的影响，点 5影响较小5.3 变室外温度的模拟结果分别设置送风温度 20、24、28、32 和 36℃，由图 7 可知，随着室外温度的升高，点 10，5 的温度略有增加，室内热环境受室外温度影响较小6.结论 采用 Airpak 软件对某实验室热环境数值模拟，经实验验证结果表明基本吻合，模拟结果可靠。

利用经验证后的数值模拟体系进行一些列变参数模拟结果表明，随传热系数增加，室温提高，当传热系数增加到 2.5 W/(K×m2)以上后，室温影响减弱，这是由于室内热源较大，墙体热工参数影响相对减弱所致此外室内温度受送风参数影响较大通过论文研究表明，借助一定的实验，利用数值模拟研究室内热环境是一种比较有效、可靠的研究方法，其研究成果可为空调设计提供参考和指导参考书目： [1] 黄晨等.大空间建筑室内垂直温度分布的研究.暖通空调.1999,No.5.[2] 赵荣义等.空气调节.北京: 中国建筑工业出版社,2002.151-156.。