通风管道内温湿度对颗粒沉积的影响

第 3 2 卷 第 4期 2 0 1 0年 0 8月 土 木 建 筑 与 环 境 工 程 J o u r n a l o f Ci v i l Ar c h i t e c t u r a l& En v i r o n me n t a l En g i n e e r i n g V o 1 3 2 N O 4 Aug 2 01 0 通风管道 内温湿度对颗粒沉积 的影 响 韩 云龙 ， 胡永梅 ， 钱付 平 ( 安徽工业大学 建筑工程 学院， 安徽 马鞍 山 2 4 3 0 0 2 ) 摘要 ： 于矩 形断 面 0 3 m 0 2 m、 长度 3 m 的通 风 管道 中 ， 对 完全 发展 湍 流 中的 颗粒 沉 积进 行 了模拟计算。由于空调通风管道 内空气温、 湿度有所变化， 因此综合考虑 了热泳力及湿度对颗粒沉 积 的影 响 ， 采 用拉 格 朗 日法随机轨道 模 型进行 了 3 0 0 0个颗粒沉积 的模 拟研 究 。在 粒径 范 围0 0 1 5 0 m、 管道风速 5 m s的条件 下 ， 得 到 了无 因次 沉积速 率 随无 因次松 弛 时 间变化 曲线并 与 前人 的 研 究结 果进行 了比较 ， 结果表 明 曲线呈现“ V” 形 分布 ， 模 拟 结 果与前 人 研 究结 果基 本 一致 。针 对 粒 径 1“ m的颗粒 ， 调查 了温度及湿度对颗粒沉积的影响， 结果表明由气流与管壁 的温差产生的热泳 力加速 了颗粒沉积， 并随温差的增加 而呈提 高趋势； 随着空气相对湿度的提 高， 颗粒沉积速 率也相 应增加 。 关 键词 ： 空 气品质 ； 沉积 ； 热 泳力 ； 通 风管道 中图分类 号 ： TU8 3 4 6 文献标 志码 ： A 文 章编号 ： 1 6 7 4 4 7 6 4 ( 2 0 1 0 ) 0 4 0 0 6 6 0 5 Ef f e c t s o f Ai r Te mp e r a t u r e a nd Hu m i d i t y o n Pa r t i c l e De p o s i t i o n i n a Ve n t i l a t i o n Du c t HAN Yu n l o n g，HU Yo n g - me i ， QI AN Fu pi n g ( S c h o o l o f Ci v i l En g i n e e r i n g a n d Ar c h i t e c t u r e ，An h u i Un i v e r s i t y o f Te c h n o l o g y，M a a n s h a n，2 4 3 0 0 2，P R Ch i n a ) Ab s t r a c t ：Pa r t i c l e de po s i t i o n i n a f ul l y de v e l o pe d t u r bu l e nt f l o w o f v e n t i l a t i o n du c t wi t h a r e c t a n gul a r c r o s s s e c t i o n o f 0 3 m 0 2 m a n d 3 m i n l e n gt h wa s i n v e s t i ga t e d The e f f e c t s o f t he r mo ph or e t i c f o r c e a n d a i r h u mi d i t y we r e a l s o c o n s i d e r e d Tr a j e c t o r i e s o f 3 0 0 0 p a r t i c l e s we r e t r a c k e d o n t h e b a s i s o f t h e r a n d o m wa l k mod e l f r om La gr a n gi a n a pp r oa c hThe p r o f i l e of di me n s i o nl e s s d e p o s i t i o n v e l o c i t y wi t h r e l a xa t i o n t i me wa s a c hi e v e d a nd t he c o mpa r i s o n wi t h t h e pr e v i o us s t u di e s wa s m a de wi t h p a r t i c l e s i z e r a n ge o f 0 0 1 5 0 um a n d a i r f l o w v e l o c i t y o f 5 m s I t i s f o u n d t h a t t h e t r e n d o f d i me n s i o n l e s s d e p o s i t i o n v e l o c i t y i s i n a g r e e me n t wi t h t h e p r e v i o us s t ud i e s a nd t he pr o f i l e o f di me ns i o nl e s s de po s i t i o n v e l o c i t y s h owe d a V s h a pe d c u r ve Th e e f f e c t s o f t e mp e r a t u r e a nd a i r hu mi di t y we r e a l s o i nv e s t i ga t e d wi t h p a r t i c l e s i z e o f 1 , m The r m o p ho r e t i c f or c e c a u s e d by t e m p e r a t ur e gr a d i e nt b e t we e n a i r f l o w a n d wa l l s u r f a c e a c c e l e r a t e d p a r t i c l e d e p o s i t i o n M e a n whi l e，i t i s a l s o f ou nd t ha t t he d e p o s i t i o n v e l oc i t y i nc r e a s e s wi t h t he i nc r e a s e o f a i r r e l a t i v e h umi d i t y Ke y wo r d s ：a i r q u a l i t y；d e p o s i t i o n；t he r m o p ho r e t i c f o r c e；v e n t i l a t i o n du c t 室内空气品质 问题 ， 如 C O。浓度高或新风 不 足、 挥发性有机气体的排放 、 细( 病) 菌等更易获得大 家的关注， 尤其对于那些 大部分时间在封闭的室内 工作和生活的人来说 ， 更是如此 ， 可是大家往往忽略 了空气 中可 吸入颗粒 物浓度 的影 响 。对 于 中央空调 系统而言， 如果通风系统只是装备了低效或 中效过 滤 器 ， 那 么 室 内可 吸 入 颗 粒 物 浓 度 将 会 受 到 影 响。 室外空气 中的灰尘会通过过滤器进入送风系统 ， 其 收稿 日期 ： 2 0 1 0 0 1 0 7 基金项 目： 国家 自然科学基金资助项 目( 5 0 7 0 8 0 0 1 ) 作者简介 ： 韩云龙 ( 1 9 7 4 一 ) ， 男 ， 副教授 ， 博 士， 主要从事空气净化及污染物处理研究 ， ( E - ma i l ) h a n l o n g y 1 6 3 t o m。 学兔兔 w w w .x u e t u t u .c o m 第 4 期 韩 云龙 ， 等 ： 通风 管道 内温 湿度对 颗粒沉 积的影 响 中一部分 会沉 积在 管 道 内， 而送 风 管 道 内的 温湿 度 非常适合病菌、 微生物的生存繁殖 ， 在空调系统启动 过程中或管道振动则可能会促使气流将这些附有病 菌的积尘重新扬起并送至室内。若人们长期暴露于 此环境下 ， 呼吸道疾病或病态建筑综合症将不可避 免 。因此 ， 诸 多研究 者l_ 1 对 管道 内颗 粒沉 积 进行 了 相关的研究工作， 他们的研究结果表明风速、 颗粒粒 径、 管壁表面粗糙度等对颗粒物沉积具有重要的影 响 ， 一般 而言 ， 随风速 、 粒径 的增加 ， 颗 粒沉 积 速率 提 高。对于水平的通风管道 ， 气流中的颗粒物在重力 沉积作 用下部 分 沉积 于 管 道底 部 ， 同时 惯 性 和 湍 流 扩 散也会 加速颗 粒 向管 道周 壁的沉 积 。 对 于空调房 间 ， 为 了维持 室 内一定 的温 、 湿 度环 境 ， 要 进行 一个热 湿交换 过程 ， 经过 热湿 处理 后 的空 气经 由通风 管 道送 往 空 调房 间 。因此 ， 若 管道 保 温 措施不当， 则气流与管壁间的温差产生的热泳力会 加速颗 粒 的沉积 ， 一些 研究 者I 6 8 _ 认 为 即使气 流 与管 壁之间存在一个小的温差也会加速颗粒沉积速率。 然而 ， 前 人 。 2 的研 究 多基 于小 断 面通 风 管道 内 的 颗粒物 沉积 ， 同 时空 气 湿度 变 化 对通 风 管 道 内颗粒 物沉积 的影 响还 未有 报道 。该文采 用 雷诺 应力 模 型 ( R S M) 模 拟管道 内湍 流 ， 在 接 近实 际 空调 通 风 管道 尺寸断面的管道内应用拉格朗 日随机轨道模型对完 全发 展湍 流 中颗粒 物 沉 积进 行 模 拟研 究 ， 对 于管 道 内气 流与管 壁温 差所产 生 的热 泳力 对颗 粒 物沉 积 的 影响进行探讨 ， 同时研究空气相对湿度对颗粒物沉 积 的影响 。 1 模型 1 1流体 与颗粒 运动方 程 通风管道 内空气与颗粒之间的运动属于典型的 气 固两 相流动 问题 ， 目前 ， 对 此 问题 的解决 方法 主 要 就是欧拉法与拉格朗日法。拉格朗 日法可追踪单个 颗 粒 的运 动轨 迹 ， 但 是 需 要 花 费 较 多 的 计 算 时 间 。 一 般而言， 空调通风管道内的灰尘浓度较低 ， 由于过 滤器 的装备 ， 粒 径也 较 小 ， 因此 ， 空 气 与 灰 尘 可分 别 处理 为连续 相 与分 散 相 。通 常 来 说 ， 由于 灰尘 粒 径 较小及浓度较低 ， 连续相对分散相有重要 的影响， 而 灰尘颗粒对连续相的影 响可 以被忽 略， 所 以采用拉 格朗 日法 的单 相耦合 计算 方法 。 空调通风管道内湍流流动的模拟计算应用商业 软件 F L UE NT 6 1 ， 一些 研究 者 1 采用 雷诺 平 均 的 Na v i e r S t o k e s( R ANS ) 方程 ， 即 k 一双方程模 拟不可压缩湍流流动 。Ti a n和 Ah ma d i E 1 6 于湍流管 道 内纳米 和微 米级 的颗 粒 物沉 积 进 行 了模 拟 研 究 ， 认 为 R S M 湍 流模 型 及 “ t wo l a y e r ” 边 界 条 件 的使 用， 能合理 的预 测颗粒 物 的沉积 。模 拟结 果表 明 0 0 1 5 0 m 范 围的无 因次颗 粒 沉积 速率 呈 现 “ V” 形 分布 ， 因此 ， 采 用 R S M 模 型 和 “ t w o l a y e r ” 模 型 的边 界处 理条件 。 1 2沉 积速率 及热泳 力 为了与前人的研究结果进行 比较 ， 采用 了无因 次颗 粒沉积 速率 ， 通常 表示为 ： 一 ( 1 ) 其 中 C 为 时均 颗 粒 浓 度 ； 为 磨擦 速度 ， 可 写为 下式 ： “ 一 U ( 2) 其 中 U 为 平 均 风 速 ；f 为 范 宁 系 数 ， 采 用 wh i t e l 1 得 出 的关 联式计算 ，。 1 _ 厂=一3 6 1 g 6 9 +( k 3 7 D) ( 3 ) 其中 k为壁面粗糙度， 光滑壁面为 0 i R e为雷诺 数 ； D为水 力 直径 。 无 因次 松弛 时间 r ， 如下示 ： 一 ( 4) 其 中 和v 分别为流体的动力粘度与运动粘度 ； 和 d 分 别为 颗粒密 度 和直 径 ； C 为 坎 宁汉 系 数 ， 如 下示 ： C 一 1 + K 1 2 5 7+0 4 e x p ( - ) ( 5 ) 其 中 K 为克努 森数 ： K 一 _ ,A ( 6) “。 其 中 为气体 分子 自由程 ， 在 1大气压 ， 2 5 时 为 0 0 6 5 m。 如果气 流与 管壁 之 间存 在 温 差 ， 颗 粒物 会 沿 着 温度降低的方向运动 ， 即为热泳现象。颗粒所受热 泳力 可表示 为如下 ： F 一一 D ( 7 ) H l n l d 工 t 其 中 为热泳力 ； D 为 热泳 系数 ， 其确 定采用 Ta l b o t 等 提出的关联式； 。 为颗粒质量 ； 丁为流 体 温度 。 2 计算域 于矩形断面 0 3 m0 2 m， 长度为 3 m接近实 际空 调通风 管道尺 寸 的管道 内模拟 完全 发展 湍 流颗 粒沉积。采用网格生成工具 Ga mb i t 生成六面体结 构网格 ； 对于管道内颗粒沉积而言， 边壁处高精度 网 格 的生 成 十分 重要 ， 因此采 用 T i a n和 Ah ma d i E “ 的 学兔兔 w w w .x u e t u t u .c o m 6 8 土 木 建 筑 与 环 境 工 程 第3 2 卷 边界层划分 方法， 即第一个 网格 距离边壁 处 0 0 5 mm， 其后 以 1 2的倍 率增 加 直 至 1 1 2个单 位 长 度 为止 ， 近壁处 网格 较 为 密 集 ， 见 图 1 ， 整 个 计 算 域 共划分 1 5 0 0 0 0个网格单元。 图 1矩形 断面 网格 划 分 应 用拉格 朗 日随机轨道模 型追 踪离 散 相颗 粒 的 运动轨迹 ， 颗粒 密度为 9 0 0 k g m。 ， 个数 为 3 0 0 0个 。 为 了和以前研 究者 的实 验 或计 算 结 果进 行 比较 ， 对 于粒径 范围 0 0 1 5 O u m 的球形 颗粒沉 积进行 了模 拟计算 ， 入口风速为 5 m s ， 温度 1 5 ， 管道出 口边 界条件为“ o u t f l o w” ， 管壁热绝缘并认为水力光滑 ， 为“ t r a p ” ； 忽 略颗 粒间 的合 并、 碰撞 及与壁 面的反 弹。 3 结果与讨论 3 1流场分 析 为 了和 前人 的研 究 结 果进 行 比较 ， 对 小 断 面 尺 寸 ( O 0 7 5 m0 0 5 m) 的通 风管 道 内低 雷诺 数 流 动 及颗粒 沉积 也进行 了模拟计 算 ， 流场 见 图 2 。2种 断 面尺 寸的管 道进 口风 速 均 为 5 m s ， 由于 断 面 尺 寸 的 减 小 ， 相 应 的 雷 诺 数 也 减 小 ， 分 别 为 8 2 1 3 5 、 5 1 3 6 。 由图 2可 以看 出 ， 虽 然 进 风风 速相 同 ， 但 断 面 尺寸 的减小导 致雷诺 数大 幅降低 。可 以看 出大 断 面 管道 内流速梯度 较小 ， 断面流 速较 为均 匀 ； 而小 断 面 管道流速梯 度较 大 ， 断 面越小 ， 管道 流 动核心 区流速 越 大 ， 而 近壁 区域 流速 明显 低 于大 断 面 管 道处 的流 速 。而前 人如 T i a n等 1 。 主要 基 于 小 断面 管道 内颗 粒沉积的研究 ， 由此可知这和实际通风管道内流体 流 动情况稍 有差距 。 3 2模型 验证 图 3 为 无 因次 沉积 速率 随无 因次 松 弛时 间 的变 化曲线， 尘粒粒径范围为 0 0 1 5 O m， 并将模拟结 果 与前人数 据进 行 了 比较 。由 于粒 径 不 同， 所 受 扩 散力 、 惯性力大小也不同， 致使颗粒沉积运动呈现 3 个 区 ： 即扩散 区 、 过 渡 区和惯 性 区 。对 于次 微米 及 纳 米级颗粒而言， 颗粒所受重力很小 ， 惯性力可 以忽 略 ， 而扩散 力则 是 颗粒 沉 积 的主 要 机 理 。对 于 大 粒 径颗粒， 重力 及惯性力是颗粒沉积 的主要作 用力。 由图 3可 以看 出 ， 模 拟 结果 基 本 符 合 前人 的实 验 数 据及计算结果， 无 因次沉积速率随无因次松弛时间 呈现 “ V” 形变化 ， 满 足颗 粒 沉积 的扩散 区、 过渡 区及 惯性区分布规律。但是计算域为矩形断面( O 3 i n 0 2 m) 的三 维计算 ， 更 接 近 于 实 际空 调 管道 的断 面 尺寸 ， 雷诺数 ( 8 2 1 3 5 ) 也要远高 于前人的实验或计 算 中采用的雷诺数口 “ ( T i a n等 采用的雷诺数为 6 6 6 7 ) 。但雷诺数 的提 高不是 由风速 的提高 而引 起 ， 而仅 由断 面尺寸 的增加 导致 的 。一般 而言 ， 随风 速 的提 高 ， 颗粒 的沉积 速率 会相应 提 高 ， 但 由断 面尺 寸增加而导致雷诺数提高对颗粒沉积的影响并不显 著 。 i 围 -瞄 ( b) 0 0 7 5m 0 0 5m 图 2不 同通风管道断面流速场 r 图 3 模型计算 结果 的比较验证 3 3温度对 颗粒 沉积 的影响 对于实际空调系统通风管道而言， 送风温度 、 湿 度因季节的变化而不同， 因主气流温度与管壁的温 差而导致的热泳力将会促进次微米级小微粒 的沉 积。有关热泳力对颗粒沉积的影 响， 一些研究者曾 加 以研 究 ， 如 Ya n g等 l_ 1 l _ 认 为 P M2 5的沉 积效 率 主 矾 拍 卵 勰 6 5 5 5444 3 3 3 32 2 2 ll OOO 学兔兔 w w w .x u e t u t u .c o m 第4 期 韩云 龙 ， 等 ： 通风 管道 内温 湿度 对颗 粒沉 积的影 响 要受主气流 与冷壁 面的温差 影 响， 热泳 力是 导致 P M2 5沉积 的关键 因素。He和 Ah ma d i _ 2 。 。 认为 只 有颗粒 进入边 界 层 的粘性 底 层 时 ， 布 朗 扩 散才 是 颗 粒沉积 的控制 机理 ， 热泳力 对小 颗粒 ( 如几 个 微米 级 或以下) 的传输 、 沉积有重要影响。并且沉积速率及 粒径范 围皆随着 主 气 流与 壁 面 温差 的提 高 而 增 加 。 其 它研究 者_ 2 船 也都 认 为温差导 致 的热泳 力 能促 进 颗粒的沉积。该研究针对 1 g m 的颗粒探讨 主气流 与管壁 之间温 差对颗 粒沉 积 的影 响 。图 4为 入 口空 气 温度 2 4 时 不 同壁 面 温 度 所 形 成 主 气 流 与 冷壁 面的温差 对颗 粒沉 积 的影 响 ， 温度 梯 度 均 为 沿 着壁 面的法线 方 向 ， 温 度 向管 壁 处 降 低 。 由 图 4可 以 看 出， 当固定入口温度 时， 随着管壁温度的提高， 颗粒 沉积速率有下降的趋势 ， 无论底面 、 侧壁还是顶壁面 都是 如此趋 势 ， 表 明 了主 气 流 与壁 面 温 差 所 产 生 的 热泳 力促进 了颗粒 的沉 积 。 同时 可 以看 出 ， 侧 壁 颗 粒沉 积速率 同顶 、 底 面 相 比 ， 其 沉 积速 率 略小 ， 底 面 的沉积速率最大。然而对于实际空调系统通风管道 而言 ， 都会有 保 温措 施 ， 管 道 壁 面 可视 为 热 绝 缘 ， 主 气流与管壁之 间的温差很小 ， 那么热泳力对颗粒沉 积所起 的作用 可 以被忽 略 。 管 壁 溢 度， 图 4热泳力对颗粒沉积的影晌 然而当入 口空气温度提高， 而管壁保持绝缘时 ， 发现 不 同之 现 象 ， 如 图 5示 。可 以发 现 随 着入 口温 度的提高， 底面、 侧壁的颗粒沉积速率均略呈下降趋 势， 而顶面的颗粒沉积速率则呈现升高趋势。虽然 设立 了管壁 绝缘 的边 界 条件 ， 但 管 道 顶 部 的气 流 温 度略 高于管道 底部 ， 会产生 一个 升力 。因此 ， 温差 产 生的升 力 作 用 使 较 多 的颗 粒 分 散 于 管 道 空 间 的 上 部， 所以， 管道顶部 的颗粒沉积速率有升高趋势， 底 部 及侧壁 的颗 粒 沉 积速 率 则 略 呈 下 降 趋 势 。另 外 ， 随着温度 的升 高 ， 空 气 的粘 度 增加 而导 致 管 道 内气 体流动的雷诺数有所减小 ， 湍动能有所 降低致使颗 粒 的湍 流扩散作 用减 弱也有 关 。 图 5 管道气流温度对颗粒沉积的影响 3 4湿度对 颗粒 沉积 的影 响 空调系统为了满足 室内一定的湿度要求 ， 在新 风处理 单元 一般 需 要加 湿 过 程 ， 经 过 加湿 处 理 后 的 空气 ， 其 中的颗 粒沉 积 是 否 有 影 响 尚需 研 究 。在 该 研究 中， 于人 口空气 1 6 时， 对通 风管道 内空气相 对 湿度 4 0 8 O 的环 境 条件 颗 粒 沉 积 进行 了模 拟 研 究 ， 见 图 6 。空气 相 对 湿 度 的 提 高会 改 变 空气 的粘度 、 密度 、 热导率等物性， 将空气相对湿度改变 所 引起 的物性 变化计 算 出来 ， 在模 拟计 算 过程 中 ， 认 为每个物性都对应着一个相对湿度。由图 6可以看 出， 随着气流相对湿度的提高， 颗粒沉积速率均相应 增加 。这是 由于空气相对 湿度 的提高 ， 湿空气的粘 度也会相应提高， 增加了湍流边界层的厚度 ， 使得更 多 的颗粒进 入 这个 低 速 的粘 性 边 界 层 ， 因此提 高 了 颗 粒沉积 速率 。该研 究只是 对于绝 缘 管壁 的情 况进 行的计算 ， 在实际运行 中， 倘若管壁保温措施不 当， 则较高湿度空气 中的水会凝结于管壁， 这就会相应 提 高管壁 的粗糙 度 ， 按 照 L a i 等 的研究 结果 可 知 ， 粗糙度的提高会促进颗粒 的沉积， 所以沉积速率会 有相应 的提 高 ， 而潮 湿 的管道 内环境 ， 适 宜 的温度 则 为真 菌 、 细 菌等提供 了繁 殖 的温床 ， 对 室 内空 气 品质 会 有严 重 的影响 。 图 6管道 内气流湿度对颗粒 沉积 的影响 4 结 论 应用拉格朗日随机轨道模型模拟分析了通风管 5 5 4 4 3 3 2 2 l 一 ； g一 一 瓣 器 学兔兔 w w w .x u e t u t u .c o m 土 木 建 筑 与 环 境 工 程 第3 2 卷 道 内温 湿度变 化对 颗 粒 物 沉积 的 影 响 ， 得 出 了如 下 主要结论 ： 1 ) 集 中式 空调 系 统通 风 管 道 内高温 送 、 回风 气 流与低 温管壁 形成 较 大 温差 时 ， 热 泳 力 会加 速 小 颗 粒 ( 亚微 米 、 纳米 级颗粒 ) 的沉 积速率 ， 对 于实 际 工程 应用 ， 若 管道具 有 良好 的保 温措施 ， 热泳 力对 颗粒 沉 积 的作 用则可 以被忽 略 。 2 ) 通风管 道 内送 、 回风 湿度 的提 高会 加 速 颗 粒 的沉积速 率 ， 对 于缺乏 良好保 温措施 的通风 管道 ， 若 湿度较 大 ， 管壁 上凝结水 的 出现会增 加 粗糙 度 ， 更 会 加速颗粒 物 的沉 积 。 3 ) 集 中式空 调 系统 通风 管道 内颗粒 物 的沉 积及 其卫生状 况 对 室 内 空 气 品 质 的 影 响 已逐 渐 引 起 重 视 ， 若通 风管道维 护不 当 ， 送 回风温 湿度 的变 化会 影 响颗粒 物的沉 积 ， 并 恶化 室 内空 气 品质 。通 过 文 章 的分析研 究 ， 有 助 于 空调 运行 管 理 人 员 提高 通 风 管 道 内颗粒 物沉 积 的认 识 ， 促进 其 对 空 调 系统 通 风 管 道 的维护 、 清洁工 作 ， 同时可 为评价 室 内人员 对颗 粒 物 的接触 、 暴露情 况提供 参考 。 参考文献 ： 1C HE O NG K W D e p o s i t i o n o f a e r o s o l p a r t i c l e s i n d u e t w o r k J Ap p l i e d E n e r g y ，1 9 9 7 ， 5 7 ( 4 ) ： 2 5 3 2 6 1 2 wu J ， Z HAO B E f f e c t o f v e n t i l a t i o n d u c t a s a p a r t i c l e f i l t e r J B u i l d i n g a n d E n v i r o n me n t ，2 0 0 7 ，4 2 ： 2 5 2 3 2 5 29 r 3I AI A C K， BYRNE M A， GODDARD A J H Enh an c e d pa r t i c l e l os s i n v e nt i l a t i o n d uc t wi t h r i bbe d s u r f a c e J B u i l d i n g a n d E n v i r o n me n t ， 2 0 0 0 ， 3 5 ： 4 2 5 43 2 4 Z HAO B ，wU J Mo d e l i n g p a r t i c l e d e p o s i t i o n f r o m f u l l y d e v e l o p e d t u r b u l e n t f l o w i n v e n t i l a t i o n d u c t J At mo s ph e r i c Env i r o nme n t ，2 0 06，4 0：4 57 46 6 5 Z HAN G H，AHMAD I G Ae r o s o l p a r t i c l e t r a n s p o r t a nd de p o s i t i o n i n v e r t i c a l a n d ho r i z o nt a l t ur bul e nt d uc t f l o ws J J o u r n a l o f F l u i d Me c h a n i c s ， 2 0 0 0 ，4 0 6 ：5 5 8 0 r 6B AE G N，L EE C S，P ARK S O Me a s u r e me n t s a n d c o nt r ol o f p a r t i c l e de p os i t i o n ve l o c i t y o n a ho r i z on t a l wa f e r w i t h t h e r r n o p h o r e t i c e f f e c t E J Ae r o s o l S c i e n c e a n d Te c hn ol o gy，1 9 95，2 3：3 21 - 3 30 r 7W I RZ BE RGER H，L EKHLANAKHER S，SHAP I RO M ，e t a 1 Pr e ve nt i on o f pa r t i c l e de p os i t i o n b y me a ns of h e a t i n g t h e d e p o s i t i o n s u r f a c e J J o u r n a l o f A e r o s o l Sc i e n c e，1 99 7，2 8( S1)：$83 - $8 4 8 TS A I R，I I ANG L J C o r r e l a t io n f o r t h e r mo p h o r e t i c d e p o s i t i o n o f a e r o s o l p a r t i c l e s o n t o c o l d p l a t e s J J o u r n a 1 o f Ae r o s o l S c i e n c e ，2 0 0 1 ，3 2：4 7 3 4 8 7 r 9CHANG Y C， RANADE M B， GE NTRY J W The r mo ph or e t i c d e p os i t i on i n f l o w a l on g a n a n nul a r 1 0 c r o s s s e c t i o n ：e x p e r i me n t a n d s i mu l a t i o n J J o u r n a l of Ae r os o I Sc i e nc e，1 99 5，2 6( 3)：4 07 - 42 8 LUO X W 。 YU S Y De p os i t i on of p a r t i c l e s i n t u r b u l e n t p i p e f l o w J C h i n a P a r t i c u o l o g y 2 0 0 6 ，4 ( 1)： 3 1 - 3 4 r 1 1 YANG R C，LI U R L，Z HOU T，e t a 1 Ki n e ma t i c c ha r a c t e r i s t i c s a nd t he r moph or e t i c de po s i t i on of i n h a l a b l e p a r t i c l e s i n t u r b u l e n t d u c t f l o w J C h i n e s e J o u r n a l o f C h e mi c a l E n g i n e e r i n g，2 0 0 8 ，1 6( 2) ： 1 9 2 一 l 9 7 1 2 D E HB I AA s t o c h a s t i c L a n g e v i n mo d e l o f t u r b u l e n t pa r t i c l e d i s pe r s i o n i n t he p r e s e nc e o f t he r mo ph or e s i s J I n t e r n a t i o n a l J o u r n a l o f Mu h i p h a s e F l o w，2 0 0 9 ， 35： 21 9 - 2 2 6 1 3 Z HANG Z ， C HE N QE x p e r i me n t a l me a s u r e me n t s a n d nu me r i c a l s i mul a t i on s of pa r t i c l e t r a ns p or t a n d d i s t r i b u t i o n i n v e n t i l a t e d r o o ms J A t mo s p h e r i c En vi r on me nt ，2 006，4 O：33 96 - 3 40 8 1 4 Z HA NG J P ，u A GS t u d y o n p a r t i c l e d e p o s i t i o n i n v e r t i c a l s qu a r e v e nt i l a t i o n d uc t f l o ws b y di f f e r e n t mo d e l s J E n e r g y C o n v e r s i o n a n d Ma n a g e me n t ， 2 00 8，4 9：1 0 0 8 1 O1 8 1 5 L AI A C K，C HE N F Z Mo d e l i n g p a r t i c l e d e p o s i t i o n a nd di s t r i b ut i on i n a c ha mb e r wi t h a t wo e qu a t i on R e y n o l d s a v e r a g e d Na v i e rS t o k e s mo d e l J J o u r n a l o f Ae r os o1 Sc i e n c e，20 06，37：1 7 7 0 1 78 0 1 6 T I AN L， AHM