摘要：由于ISO5级(百级)洁净室运行能耗较大，考虑到节能的必要性，本文利用计算流体动力学(CFD)方法对拟采用风机过滤器单元(FFU)洁净空调方案的ISO5级电子工业洁净室进行模拟，得出室内气流速度场，分析其性能，并通过理论公式计算所能达到的洁净度。我们认为结合洁净室的具体用途，通过合理布置末端FFU送风口位置及选择回风形式，以及选择较高级别的末端过滤器，可以在满布率较低时达到较高的洁净度级别。
1引言
洁净室空调系统经典的方案是采用中央空调和三级过滤器集中送风，通过大型风道将已经处理的空气送至过滤器的接联管道，然后经高效空气过滤器(HEPAFilter)或者超高效空气过滤器(ULPAFilter)送到洁净室。而另一种方案是采用室内循环风就地冷却，利用干冷却盘管解决新风不能提供全部冷负荷的问题，同时利用风机过滤器单元来进行空气循环。每种方式各有一定的适用范围，风机过滤器单元(FFU)因其灵活性大,即可通过置换盲板来提高局部区域的洁净度、占用空间较少等优点得到越来越多的应用，尤其适合于旧厂房的改造及技术更新较快的工程。虽然FFU系统成本较高，而从综合投资角度，分析认为采用FFU方式在末端过滤器铺设率为25%-30%时较为有利【1】。
ISO5级(百级)洁净室属于洁净室用暖通空调系统耗能大户，通常采用吊顶满布高效过滤器的送风方式，运行能耗较大。有关洁净室运行费用的文献指出，在某些欧洲国家，能源消耗的费用已占洁净室运行、维护年度总费用的65%～75%【2】，其主要影响因素是洁净室的空气流量和采暖通风空调系统如何有效地向洁净室分布经过净化和温湿度调节的空气，所以在保证洁净污染控制的条件下，合理选择送风速度，布置末端过滤器、回风口、减少送风量以便节能是人们关注的焦点。
另外国外对一些ISO5级洁净室实测数据表明，大部分换气次数远低于建议的下限值【2】，而在设计中存在系统风量过大的倾向，这可能与对气流缺乏了解，担心系统运行可靠性的保守思想有关，说明提高节省能源的机会确实存在。随着计算流体动力学(CFD)技术自身的发展，已广泛应用于暖通空调和洁净室等工程领域，通过计算机求解流体所遵循的控制方程，可以获得流动区域的流速、温度、浓度等物理量的详细分布情况，是一种较好的优化设计工具。其优势在于利用CFD技术对设计方案进行模拟可以在施工前发现失误并及时更正，避免经济损失;可以迅速发现提高系统运行效率的可能性;另外，通过模拟可以得到一系列运行的备选方案，以便在寻找最经济方案时有所依据。
本文利用CFD软件，对拟采用FFU净化空调系统的某微电子洁净厂房的ISO5级洁净室进行计算机模拟，通过几个设计方案相比较，利用所得到的速度场，分析评价其性能，利用理论计算验证其平衡态的洁净度，并提出一些应用中的注意事项，为实际工程应用提供参考。
2数值模拟及分析
2.1数学模型
从流动的雷诺数Re来考虑，洁净室的气流均为紊流【3】，空气的流动满足连续性方程，动量方程和能量方程。对于工程问题，我们不需要关心紊流的精细结构及其瞬时变化，而只关心紊流随机变量的有关平均值，因此，本文采用数值计算三类方法中雷诺时均方程中的紊流粘性系数法，流动模型采用暖通空调广泛采用的标准k-ε二方程模型，k-ε模型通过求解紊流动能与紊流动能耗散率的输运方程得到紊流粘性系数。
控制方程的通用形式为【4】：
式中：ρ为空气密度(kg/m3)，V为气流速度矢量(m/s)，Γφ，eff为有效扩散系数(kg/ms)，Sφ是源项，Φ代表1，u,v,w,k,ε中的一项，u,v,w为三个方向的速度分量(m/s)，k为紊流动能(m2/s2)，ε为紊流动能耗散率(m2/s3)，Φ=1时通用方程变为连续性方程。
边界条件：墙体边界设为无滑移边界条件。送风边界条件，送风速度取过滤器面风速平均值，速度方向竖直向下。回风边界条件，回风口满足充分发展段紊流出口模型。由于室内热负荷较小，不考虑温度浮升效应对气流的影响。采用混合迎风差分格式对偏微分方程进行离散，基于有限容积法的SIMPLEST算法进行求解。
2.2物理模型及计算结果分析方案一将风机过滤器单元(规格为1.2m×1.2m)成条型居中布置于天花板，满布比在25%，回风采用全地面均匀散布穿孔板作为回风口。物理模型平面图如图1。经模拟计算得到气流流场示于图3，由于送风口在Y方向呈对称布置，图中只给出一半流场。从图中可见，在送风口下方流线垂直向下，流线平行较好，而在送风口至墙体范围内有较大的涡流区，则主流区范围减少，不能使全室工作区达到较高级别。同时粒子也会被卷吸进入主流区，排除污染物的路径增长，增加污染的可能性。
方案二将FFU(规格为1.2m×1.2m)散布于天花板，满布比仍为25%，过滤器面风速在0.45m/s，回风采用全地面均匀散布高架格栅地板作为回风口。物理模型平面示意图如图2，气流流场分布如图4。模拟计算显示，对于均匀布置FFU方案，工作区1.2m及0.8m高度断面平均风速分别为0.1545m/s、0.1516m/s，可见散布末端过滤器送风口可以减小速度的衰减。虽然在送风口之间上部存在反向气流，形成小的涡流区，但在工作区0.8m-1.2m范围内已形成竖直向下的流线，时均流线平行较好，由于此洁净室产热量较小，热气流对流线影响可忽略，不会产生逆向污染，因此上部的涡流不会对主流区产生影响。空气中的微粒在重力、惯性和扩散三种作用力下运动速度和位移是微小的，直径在1μm时，微粒跟随气流运动的速度和气流速度相差不会大于10-3【3】。此设计中新风处理机组设三级过滤器，FFU中过滤器为U15≥99.9995%@MPPS，直径>1μm的微粒可视为零，因此，工作区产生的微粒能完全跟随气流一起运动，直接排出洁净室。
当进一步减小满布比时模拟计算可知，除送风口正下方—定区域外，其余部分已根本不能保证气流接近垂直向下，过滤器之间存在一个从天花板到地面贯通的巨大涡流区，污染物极易被卷吸进入涡流区内而不易排出。
经过模拟计算及分析，我们认为在送风口满布比为25%，均匀分布FFU，采用全地面均匀散布穿孔板回风，过滤器面风速在0.45m/s，相应换气次数为147次/小时，由于FFU可达到较大的送风面风速，以及均匀散布穿孔地板回风口的均流作用，因为如果采用侧墙下侧回风，就会在洁净室中间下部区域形成较大的涡流三角区【5】，因此，洁净室内能够形成比较合理的气流流形，在主流区内能形成基本垂直向下的流线，但在靠近四周墙壁处，由于形成受限射流，出现涡旋，因此在布置设备时，应避免将设备靠墙壁布置，而应留有一定距离，这是洁净室施工完毕，开始投入使用时应加以注意的。另外，此设计中虽然不能形成如传统满布高效过滤器送风口而形成的全室平行气流，但美国标准IES-RP-CC012.1【6】中已认为ISO5级洁净室也可采用非单向流流型或混合流型。
3理论计算洁净度
洁净室的洁净度级别由通风系统和室内污染源所决定。可以通过数学公式对其进行计算。根据粒子平衡理论，进入洁净室的粒子有室外新风带入、循环空气带入及室内污染源。对于电子厂房室内污染源主要是工作人员的产尘，而设备产尘很小可忽略不计。从洁净室排出的粒子有回风带出及由于室内正压而渗出的粒子。可得如下方程【7】：
以上式中：Q，送风量，m3/sq，渗出的空气量，m3/s;V，洁净室的容积，m3;x，循环风的比例，此处为1;c，洁净室的浓度，粒/m3;c0，洁净室的初始浓度，粒/m3;c∞，洁净室的平衡浓度，粒/m3;c1，渗出空气的浓度，粒/m3;cout，室外新风的浓度，粒/m3;t，时间;ηout，新风过滤器效率;ηrec，回风过滤器效率;S，室内污染源，粒/秒;ε,通风效率。
新风预过滤器为F5(η=55%)，中效过滤器为F9(η=95%)，高效过滤器为H12(η=99.5%)，FFU中过滤器为U15(η≥99.9995%@MPPS);新风含尘浓度天津地区取为3×107粒/m3(≥0.5μm);身着洁净服的工作人员走动时的产尘量为1×104粒/秒•人(≥0.5μm);设同时有3人在工作;通风效率取为90%;新风比为4.42%。计算得出此设计的洁净室稳定含尘浓度为2857粒/m3(即81粒/ft3)，达到ISO5级100粒/ft3的设计要求。
4结论
通过本文的研究可得到如下结论：
1)针对电子厂房洁净室发尘量较低，室内人员较少，热负荷较小的情况，通过选择级别较高的过滤器，合理布置末端高效过滤器的位置，回风方式后，即使设计的室内换气次数、断面平均风速低于规范建议的下限值，仍可有效地滤除粒子，满足空气洁净度要求。
2)CFD是一种较好的优化设计工具，结合工程实际情况，借助模拟工具进行辅助设计是必然趋势