王昕 黄晨 黄武刚 杨建刚
摘要： 采用CFD数值模拟和现场实测的前期研究成果，针对具有上部侧墙开口和屋顶顶部开口两种上部开口形式的大空间建筑，使用PHOENICS数值模拟软件模拟了室内采用分层空调时各不同喷口高度、上部开口面积及其高度、下部开口面积等几何结构参数，以及室内负荷、上部开口背压、上部开口进风速度等运行参数多种组合工况下的室内热环境，讨论和分析了这两种上部开口形式工况下垂直温度分布、空调区平均温度、上部开口排风温度、上部开口排风量及室内排热量的差别。
关键词： 大空间建筑 室内热环境 数值模拟 上部侧墙开口 屋顶顶部开口
几乎所有大空间建筑因通风和结构的要求上部均设有开口。上部开口大致有上部侧墙开口和屋顶顶部开口两种形式。采用分层空调时，上部开口的形式和位置的不同对空调能耗和室内热环境特性的影响亦不尽相同，且差别较明显。从全年变化的室外气温看，除了冬季上部开口排风会增加室内负荷外，夏季或多或少地可以利用上部开口处的高温排风带走室内部分负荷，过渡季节则可关闭空调系统仅靠自然通风排走室内负荷，因此研究大空间建筑分层空调时上部开口等诸因素对室内热环境特性的影响尤为重要。
本文在开发和应用数值模拟预测大空间建筑室内温度场和速度场的研究基础上[1]~[3]，选用目前比较典型的侧喷送风方式，并将具有上部侧墙开口或具有顶部开口两种不同上部开口形式的大空间建筑作为研究对象，以夏季现场实测工况为分析基础[4]，重点讨论了这两种不同上部开口形式在不同工况下分层空调时的室内热环境特性的区别。其中在顶部开口的工况模拟中部分借用了侧墙开口工况的实测结果。
2 计算条件2.1 建筑模型
图1为数值模拟用某体育馆简化模型。建筑柱型部分直径68m，高18m，屋顶呈扁球体，信捷职称论文写作发表网，净高为8m，左右两侧为阶梯型观众席，室内采用中侧送风，由38个喷口组成环形对中喷射，其中28个短程喷口倾斜12°布置，负责观众席空调，10个长程喷口水平布置，负责场内中央的空调。回风采用以台阶均匀回风为主、侧墙回风为辅的方式。模拟主要基础参数见表1夏季游乐活动实测日工况[4]。上部开口形式分别模拟为上部侧墙开口或屋顶顶部开口。工作区入室大门及渗透缝隙则模拟为下部开口。
2.2 数学模型及其边界条件
表1 夏季现场实测日工况
送风量：34.57kg/s
回风量：31.31kg/s
送风温度：16.5℃
室外气温：36.5℃
环形外走廊平均温度:32℃
下开口进风温度：28℃
日射量：767W/m2
人体负荷：18.1kW
照明负荷：45kW
传热系数/W/（m2K）：
屋顶：2.75 外墙：2.48
内墙：2.59 楼板：2.21
选用Lam-Bremhorst低雷诺数K-ε模型[5]，采用第一类边界条件。边界条件以热平衡为基础，利用Gebhart吸收系数，借助现场实测数据，通过建立如下壁面i的壁温方程组求解[3]：
，
式中： ── i表面对流放热系数，W/(m2.K)；
qi，ti/Ti── 分别为i表面内侧壁温和相应的空气温度，℃/K；
Si，── i表面面积，m2；
QLfU，QLfD ── 分别为室内空中平面辐射热源向上和向下辐射分量，W；
Gij，GUj，GDj ── 分别为i表面、平面辐射热源上表面、平面辐射热源下表面对j表面GEBHART吸收系数。
eI ── i表面发射率；
s ── 玻尔兹曼常数，W/(m2.K4)；
qli ── i表面导热散热，W/m2。
2.3 计算工况与室内热环境参数
上部侧墙开口影响室内热环境的主要因素有喷口高度、上部开口面积及其高度、下部开口面积等几何结构参数，以及室内负荷、上部开口进风速度等运行参数；顶部开口影响室内热环境的主要因素有喷口高度，上部开口背压，上下开口面积等。通过模拟计算对室内垂直温度分布、空调区温度、上部开口排风温度、室内通风排风量及排热量随上述因素的变化进行分析与讨论，以得出上述因素对室内热环境特性的影响及其规律。上部侧墙开口与屋顶顶部开口的分析讨论分别在文献［6］和文献［7］中详述。本文仅针对这两种不同上部开口形式下的室内热环境特性参数变化的不同之处进行分析和阐述。表2列出了本文讨论的两种计算工况的变化参数。讨论中室内垂直温度分布以比赛内场为讨论对象。空调区温度为内场垂直方向上的空气平均温度，并定义不偏离平均温度1%时的区域为等温空调区，其高度为等温空调区高度。室内通风排热量为室内因上下开口引起的总通风排热量，定义通风排热为正，得热为负。室内通风排热量是衡量建筑在开口作用下，室内自然通风总排热的状况，它不仅反映了上部开口排热量的大小，也反映了下部开口进风带入室内热量的大小。无论开口结构形式如何变化，建筑开口所引起的室内通风排热量越大越好，它是开口节能性的一个标志。
表2 上部开口形式计算工况
上部侧墙开口
屋顶顶部开口
变化参数
参数计算工况
记号
参数计算工况
记号
喷口高度/m
13/15
N13/N15
11/13/15/17
N11/N13/N15/N17
上部开口高度/m
15/17/19
IN15IN17/IN19
——
——
下部开口面积/m2
0/0.93/1.86/3.72
0F/1F/2F/4F
0.32/0.644/1.365
FXQ、FXH、FX
上部开口面积/m2
56/28/10.45
SO/PO/WO
4/9/16
T2-O3/T3-O4/T4-O6
图2～图4中其他记号说明：
ALLEXH：指上部开口均处于排风状态
1Q：指室内热负荷为现场实测日测定值
-5P/0P/+5P：分别指顶部开口处背压为－5/0/5Pa
3 计算结果与分析3.1 垂直温度分布
在上部侧墙开口的计算工况中，上部开口面积对垂直温度分布的影响相对较大。开口面积越大，下部空调区温度则越低。喷嘴高度对屋顶附近温度影响较大，喷嘴高度越高，屋顶附近温度越低。见图2(a)、(b)。
在屋顶顶部开口的计算工况中，喷嘴高度对垂直温度分布的影响最大，其次是上部开口面积。喷嘴高度15m时，屋顶附近温度最低，等温空调区温度处于较低的水平。此外，上部开口面积较小时，非空调区温度较高。见图2(c)、(d)。
两种不同上部开口形式的计算工况下，垂直温度分布的规律大致是相同的。但上部侧墙开口时，垂直温度明显高于顶部开口工况5～10℃左右。
3.2 空调区平均温度
在上部侧墙开口的计算工况中，空调区平均温度随下开口面积的增加而增加，增幅不大；随上开口面积增加，空调区平均温度呈增加趋势，上开口面积大于28m2，增幅已不明显。不同喷嘴高度，温度变化趋势相同，高度增加，温度亦升高。见图3(b)、(d)。
在屋顶顶部开口的计算工况中，空调区平均温度受喷嘴高度影响最大，其次下部开口面积。喷嘴高度在11m~15m时，高度每增加1m, 等温空调区温度增加约为0.25℃，15m~19m时，喷嘴高度每增加1m，等温空调区温度增加约为0.5℃。而随下部开口面积变化，温度呈递增趋势。见图3(c)、(d)。
两种不同上部开口形式的计算工况下，空调区温度总体上在侧墙开口时的温度要高于顶部开口时，两者变化规律具有相似之处：随喷嘴高度增加，等温空调区温度都升高，同时随上部开口面积增加，空调区温度有增加趋势，但增幅不明显。随下部开口面积增加，空调区温度亦有增加趋势。
（a）不同喷嘴高度及上部开高度
（b）不同上部开口面积
（c）不同喷嘴高度
（d）不同上部开口面积
图2 侧开(a)(b)与顶开(c)(d)垂直温度分布比较
图3 侧开(a)(b)与顶开(c)(d)空调区温度/排风温度比较
3.3 上部开口排风温度
在上部侧墙开口的计算工况中，喷嘴高度升高，排风温度降低。随上部开口高度
