工业厂房自然通风的数值模拟.docxVIP

工业厂房自然通风的数值模拟 0 热分布系数的确定 随着工业产品质量要求的提高，对工业工厂环境质量的要求越来越高。很多地区的工业厂房在夏季和过渡季节都采用自然通风形式, 以减少建筑运行能耗。而提高厂房内环境质量的重要环节, 则是加强和改善厂房的自然通风设计, 组织好厂房内部的进、排气流, 使厂房内获得应有的新鲜空气, 并带走影响产品质量和工人健康的大量余热和有害气体。因此, 如何作好厂房的自然通风设计, 已成为厂房设计急需解决的关键课题。 在国内针对带内热源的工业厂房自然通风的研究还不多。在对自然通风进行设计或者计算时, 仍按照文献所提供的方法, 排风温度通过计算热分布系数m(如果室内的散热量为Q, 直接进入到工作区的那部分热量为mQ(称为有效余热量) ,m定义为m= (tn-to) / (tp-to) , 其中tn,tp,to分别为工作区温度、排风口温度、进风口温度 (可看作室外计算温度) ) 得到, 从而把求解排风温度转化成求解热分布系数m。热车间自然通风计算中应用热分布系数确定工作区温度及室内外温差、压力的方法是前苏联学者巴图林对自然通风研究的一个重要成果, 但此法确定的m值并不全面。文献利用下送风房间模型, 通过大量的实验及计算机仿真模拟, 对m值进行进一步研究, 得出其应用方程式。由于热分布系数本质上是一个描述受限空间内的热对流扩散问题的参数, 因此它与建筑、热源和通风条件有关, 是多个因素的函数, 所以m值很难用理论公式来精确计算得到。 随着数值计算和模拟技术的发展, 利用计算流体力学CFD方法对自然通风过程进行系统全面模拟研究逐渐引起了很多研究人员的注意。Linden等人研究了单空间建筑置换式自然通风, 提出了单空间Emptying Filling Box自然通风模型。Li等人则提出了Emptying Air-Filling Box, 并对风压和热压作用下的自然通风作了大量研究。Fracastoro等人采用CFD方法对单区和两区建筑热压驱动的单侧通风进行了研究。Kindangen用实验方法和CFD方法研究了自然通风建筑中不同自然通风条件对室内气流分布和流速的影响。此外, 居住区室内外环境的自然通风问题也成为近年来研究的热点。 笔者利用CFD方法对具有内热源的工业厂房进行了自然通风模拟研究。在热源参数一定的情况下, 针对工业厂房的建筑结构对自然通风的影响, 分别改变进风口面积、进风口离地高度、矩形天窗喉部开口面积及排风口面积等, 对具有典型特征的工业厂房自然通风流场进行了数值模拟分析, 得出了不同计算工况下的通风量、排风温度、工作区的平均温度、热分布系数m值及其变化规律, 为工业厂房自然通风系统的设计提供了理论依据。 1 物理模型和数值方法 1.1 压力和速度关系 以上海某工业区的厂房为研究对象, 如图1所示, 该厂房的基本尺寸为50 m×30 m×10 m (长×宽×高) 。一般建筑内人员的活动高度在2 m以下, 所以将工作区的高度设为2 m。参照与模型面积相同的轧钢厂中棒材厂房的炉体表面散热量及根据文献中工业炉的散热量计算公式确定出模型中热源的散热量为500 kW, 热源表面平均温度为300 ℃, 热源尺寸为40 m×7 m×1.5 m (长×宽×高) 。室外大气压为标准大气压101 325 Pa;根据《工业企业设计卫生标准》 (GBZ 1—2002) , 车间作业地点夏季空气温度应根据实际的本地区夏季通风室外计算温度确定, 因此取上海夏季通风室外计算温度32 ℃ (305 K) 为室外环境温度。 城市风环境是个复杂的系统, 风压的作用不是一个可靠的稳定因素, 在一般工程设计中不考虑。然而在实际建筑中的自然通风是热压和风压共同作用的结果, 且两者共同作用时并不是简单的线性叠加。如果单一地研究热压而忽视风压, 就无法研究风压和热压共同作用对自然通风的影响。因此, 在理论上对室内热环境所进行的分析中, 必须考虑风压的作用。通过分析, 城市内部风环境中风速参数在0.5～3 m/s之间变化。因此取不利气候条件下室外平均风速0.5 m/s为进口风速;迎风面采用速度进口的边界条件, 由于背风面风向预先不确定, 因此定义为压力进口边界条件;迎风面的空气动力系数为0.8, 背风面的空气动力系数为-0.2;由于一般工业厂房的矩形天窗两边都设有挡风板 (为了消除外部风环境的影响, 不至于使气流从矩形天窗倒灌入厂房) , 因此排风口设置为自由出口。 1.2 模型网格划分 本文研究的流场属于自然对流湍流流场。通过大量的试算和模拟结果对比分析, 选用RealizableK-ε模型为湍流模型, 壁面附近采用标准壁面函数方法进行处理。另外考虑到壁面和热源之间存在辐射换热, 在模型中还需要加上DO辐射模型。采用二阶迎风格式对方程进行离散, 采