作者：张毅
为了简化计算，对洁净室的实际条件进行了如下假设：
①室内气流流动为稳态流动，室内气流为不可压缩流体，物性为常数，忽略质量力；
②室内无内热源，围护结构绝热，对于洁净室来说，可假设为无温差送风，而且将室内温度场视为均匀温度场；
⑧室内污染源的发尘速率恒定。
1.1物理模型
洁净手术室布置图1所示。
为保证工作面流速在0.2～0.3m/s[4]，中心区出口流速选为0.38m/s，对于带气幕的手术室具体边界条件如表1所示。
根据Woods等人估计[5]：手术台前手术人员散发典型碎屑平均直径约7.5 um，因此，模拟中散发的携带微生物的粒子粒径按7.5u m考虑，散发率采用2.5 x104粒／s的面污染源，这与一般洁净室设计规定在人穿着洁净服剧烈活动时，全身散尘率为5.6x105粒/min相比还是比较安全的。同时，污染源的位置对局部洁净室浓度场有一定关系。为了更好地说明问题，考虑污染场最不利情况，选择将面污染源放在送风口两侧的吊顶上，这也可减少面污染源对流型的影响。
1.2数学模型
洁净手术室的换气次数较大，气流组织基本可视为强制对流流型且不考虑有温差送风，所以本文的计算模型采用标准的高雷诺数k--
二方程模型其控制方程如下：
式中：妒是普适标量，可以表示三个速度分量
，湍流动能k，耗散率
是空气密度，kg/m3；
为耗散项。
颗粒物使用拉格朗日法计算，由于分散相的总体积分数较小，CFD采用DPM模型进行计算。颗粒物浓度使用质量密度，kg/m3。靠近壁面处采用压力壁面函数。离散方法采用有限差分法。在划分网格时，使用交错网格，使用非结构网格对控制体进行离散。在方组求解时，对耦合方程组使用SIMPLE算法，单个方程组使用ADI逐行迭代法。
1.3计算模型准确性的比较
为验证模型的准确性，将无气幕时工作区断面的模拟计算的速度与实验数据[2]进行比较，实验一共采集30个测点的数据，断面平均速度计算值为0.251m/s，实测值为0.234m/s，偏差率仅为7.2%，由图2可以看出两者数据变化趋势基本相同（四周高，中心低）。经分析，误差的主要原因是在计算中送风口出风假定完全均匀，但是在实际情况下送风口出风并非全均匀（高效过滤器自身滤纸厚度和密度不均匀，同时均流层均流效果不足）。
2 数值模拟的结果与分析
2.1气幕速度的影响
空气幕到达地面时才能完全发挥其隔离作用，此时由气幕把污染区和洁净区完全遮断；进一步提高空气幕出口风速时，隔离效果也不会再提高，在此条件下的空气幕的出口风速被称为遮断风速[1]。但是对于手术洁净室，气幕的主要作用不仅体现在隔离作用上，同时也体现在抗外界干扰能力上。此时需对气幕速度选择做进一步的分析。
现对左侧气幕施加一横向干扰，由图3(a)~图3(d)可以看出，在相同的横向扰动下，没有气幕时横向扰动入侵主流区十分严重，在相同气幕宽度下随着气幕射流力的提高，气幕射流与竖直方向的夹角减小，隔离作用在增强，从图3(d)和图3(e)可以看出在相同射流力时，宽口低速气幕与窄口高速气幕的射流与竖直方向的夹角差不多（射流轴线末端基本都在Z方向2.2m处），但窄口高速气幕在末端已形成遮断点（Z方向2.15m处），基本隔离了横向扰动，而宽口低速气幕的末端却未能形成很好的遮断点。
根据流体射流理论，射流在运动过程中各断面的动量保持平衡，即：
式中：M为气幕射流的总动量，N- s;L0、Lx、Le为气幕射流的出口流量、某断面的流量、末端流量，m3/s；v0、vx、ve为气幕射流的出口流速、某断面的流速、末端流速，m/s；F为气幕射流的射流力，N。
取一微元控制体，在某一断面上对于一微元体将有
2.2气幕宽度的影响
手术洁净室气幕作用不仅体现在其抗横向扰动能力上，而且体现在其隔离污染物能力上。先以射流力作为单独考量参数。
从图4和图5可以看出在对于不同的气幕宽度，随着射流力的加大气幕内外区的浓度比在逐渐降低。但同时在加大到一定值后，再加大对降低内外区的浓度比的效果将不再明显。
从图6可以看出在射流力相同时内外区的浓度比随着气幕宽度的增加而减少，结合图4与图5可以发现尽管射流力表达式从形式上综合考虑了流量与流速，但这不能使射流力作为气幕设计的一个独立参数。因为空气幕的隔离作用并不是像固体壁一样阻挡微粒的穿透，而是由它不断卷吸两侧空气，不断去稀释和带走卷吸进来的脏空气（因为空气幕喷口送出的也是洁净空气），使得微粒不能穿透气幕，只有少数微粒可能由气流的横向脉动进入中心区[1]。增加气幕速度的作用之一就是加大稀释的风量使卷吸进来的脏空气更快地被带走，但是通过增加气幕出口速度的方式来加大风量的同时势必会使气幕射流的卷吸作用加剧，卷吸进气幕的脏空气的量会大大增加，这就解释了为什么在气幕宽度一定时增加射流力到一定程度气幕的隔离效果便不会显著提高，而在射流力相同时增加气幕宽度却可以显著提高气幕的隔离效果，虽然两者都是增加风量，但后一种方法由于采用了较低的射流速度，故引起的卷吸作用较小。
气幕的隔离效果可由下式计算[1]：
式中：Na为主流区含尘浓度，粒／L；Nb为涡流区含尘浓度，粒／L；NS为主流区与气幕区送风浓度，粒／L;f为气幕洁净棚隔离效果固有特性项。
2.4气幕送风装置的有效性比较
为验证带气幕的送风装置的有效性，在相同风量，相同散尘量时将带气幕与不带气幕的送风装置进行比较，二者参数如表3。
经过模拟，带气幕的送风装置下，在工作面（距地面高度0.8m）的颗粒物浓度为2.156x10-9kg/m3，不带气幕的工作面颗粒物浓度为1.253x10-8kg/m3，同时由图7可以看出，不带气幕的送风装置下的颗粒物在工作面上的分布要比带气幕的多很多，不利于颗粒物的排出，增加了颗粒物入侵主流区的可能，同时虽然两者的洁净区域是差不多的（基本在X方向2到4m处），但带气幕的送风装置送风而积只是不带气幕的64.4%。
3 结论
本文通过理论分析与数值模拟的方法分析了气幕各参数对气幕抗干扰能力与颗粒物隔离能力的影响，可以得到以下结论：
1)射流力对于气幕的抗横向干扰能力有着重要影响，但并非是隔离能力的决定因素，气幕的出口流速对于射流力的影响更为重要。
2)气幕颗粒物隔离能力的决定因素为气幕射流的引带比，采用宽口低速的气幕由于可以更有效地降低引带比，故隔离效果比窄口高速气幕更好。
3)对于洁净手术室流场采用一定角度的射流可以有效降低工作区流场的不均匀度从而可以降低引带比增加气幕的隔离效果。
4)对于I级洁净手术室，若采用气幕隔离方式扩大主流区建议的气幕参数宽度为150mm，出口流速为1.10～1.24m/s．射流角度为5。。
4摘 要：
本文使用标准k-E紊流模型，基于I级标准手术室，对带气幕的洁净手术室流和污染物浓度场进行了数值模拟，并做了相应的理论分析。得出了射流力是气幕抗干扰能力的主要影响因素，引带系数是气幕隔离能力的主要影响因素的结论，对气幕风速、气幕风口尺寸、气幕倾斜角度等工作参数的影响特性做了深入探讨，并对比了相同风量下带气幕与不带气幕送风装置的有效性。最后针对I级洁净手术室内的流场给出了建议气幕设计参数。
