（气象学专业论文）生产车间内纳米二氧化钛粉尘扩散的数值模拟研究.pdf

生产下问内纳米二氧化钛粉尘扩散的数值模拟研究生产车间内纳米二氧化钛粉尘扩散的数值模拟研究摘要本文主要研究车问内纳米二氧化钛粉尘的扩散与堆积问题，通过分析车间内纳米粉尘的受力情况，结合车间内环境参数，选择最适合粉尘扩散的物理模型和求解方法。根据现场观测数据建立起生产车间的物理模型，并运用计算流体力学相应的软件g a m b i t 、f l u e n t 、t e c p l o t 等对模型进行网格划分、求解和结果处理等，分析了不同通风条件下车间内纳米二氧化钛的扩散问题，发现通风虽能显著降低室内粉尘的浓度，但是风速较低时，在排放源的下风向仍有大范围污染区存在。运用宽范围颗粒粒径谱仪对模拟结果进行了对比验证，对比显示，模型值与采样值具有一致性，粉尘的浓度量级与扩散趋势较为一致。关键词：粉尘扩散，通风条件，模式验证生产车问内纳米二氧化铁粉尘扩散的数值模拟研究n u m e r i c a ls i m u l a t i o no ft i t a n i u md i o x i d ed i s p e r s i o ni nt h ep r o d u c t i o nw o r k s h o pa b s t r a c td i s p e r s i o na n da c c u m u l a t i o no fn a n o m e t e rt i t a n i u md i o x i d ep a i r t i c l e sw e r ei n v e s t i g a t e di n t h i sp a p e r t h en a n o m e t e rd u s ts t r e s ss i t u a t i o na n dt h ew o r k s h o pe n v i r o n m e n tp a r a m e t e r sw e r ea l s oa n a l y z e di no r d e rt oc h o o s et h eb e s tp h y s i c a lm o d e ia n ds o l u t i o nm e t h o d b a s e du p o nt h em e a s u r i n gd a t a ，t h ep h y s i c a lm o d e lo ft h ep r o d u c t i o nw o r k s h o pw a se s t a b l i s h e d c o n e s p o n d i n 9 1 y ；c o m p u t a t i o n a lf l u i dd y n 锄i c s( c f d )s o r w a r es u c ha sg a m b i t ，f l u e n t ，t e c p l o tw e r ee m p l o y e dt oc o m p l e t et h eg r i dd i v i s i o n ，s o l u t i o n ，a n dt h ep o s t p l o c e s s i n go f1 e s u l t s 1 1 1 ed i s p e r s i o no fn a n o m e t e rt i t a n i u md i o x i d ep a r t i c l e su n d e rd i 行e r e mv e n t i 】a t i o ns i t l j a t i o nw a si n v e s t i g a t e d t h er e s u l t ss h o wt 1 1 a tv e n t i l a t i o ni se ! 矗吝c t i v er e s 0 1 u t i o nf o rp a r t i c l er e m o v a l ，b u tt 1 1 e r es t i l lh a dw i d er a n g ep 0 1 1 u t i o na r e ad o w n w i n dt h ee m i s s i o n ss o u r c e b yu t i l i z i n gt 1 1 ew i d e r a n g ep 叭i c l es p e c 咖m e t e r ( w p s ) ，t h ec o n t r a s tb e 铆e e nt h es i m u l a t i o nr e s u l t sa n dm em e a s u r i n gd a t aw a sp e r f o m e d i t i ss h o w nt h a tt h em o d e lr e s u l t sa n dt h es a m p l i n gv a l u ea g r e e dw e l lw i t he a c ho t h e r t h em a g n i t u d eo f d u s td e n s i t yh a dt h es i m i l a rt e n d e n c yw i t ht h ep r o l i f e r a t i o nd i s p e r s i o nk e yw o r d s ：p a r t i c l ed i s p e r s i o n ，v e n t i l a t i o ns i t u a t i o n ，m o d e lc o n f i r m a t i o n2独创声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。据我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含未获得( 洼；麴遗查基丝盂要挂型壹堕的! 奎拦亘窒!或其他教育机构的学位或证书使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。兰兰兰竺兰兰篁丑塑一兰兰三兰鲨兰塑一学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，有权保留并向国家有关部门或机构送交论文的复印件和磁盘，允许论文被查阅和借阅。本人授权学校可以将学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存、汇编学位论文。同时授权中国科学技术信息研究所将本学位论文收录到中国学位论文全文数据库，并通过网络向社会公众提供信息服务。( 保密的学位论文在解密后适用本授权书)学位论文作者签名：0 聊砍约修签字日期：6 弱7 目导师签字：彳签字开期：汐7 年月) 日生产车间内纳米二氧化钛粉尘扩散的数值模拟研究1 1 研究的背景与意义第一章引言悬浮颗粒物是指漂浮在空气中的固态和液态颗粒物的总称，j 其粒径范困约为0 1 1 0 0 岬。有些颗粒物因粒径大或颜色黑可以为肉眼所见，比如烟尘。有些则小到使用电子显微镜才可观察到。一般根据空气动力学直径( a e r o d y m m i cd i 锄e t e r ，简写为a d ) 将空气中的颗粒物分为3 类：粗颗粒物( 2 5 “m e a n图4 - 3 风速为l0 f n 如时车问内的风场和颗粒物浓度分椰f 培4 - 3w i n dn c l da n dt h ed s c r | b u e i o no f p a n i c u l a i cc o n c c n t r a t j o nl n s i d et h ep r o d u c l j o nw o r k s h o pw h c nt h ew i n ds p e e di s 。q u a l t o1om 止一憾：黜恐古警窖。一，赢。躐麓f 澎图4 _ 4 风速为2om 瓜时车i i l j 内的风场承嘲粒物浓度分巾f i g4 _ 4w i n dn e i da n dm ed 随r 洳n l l o no f p a n ；c u l a t ec o n c e n i r a t i o ni n s i d et h e$ 产f n = 氧化仕粉尘扩肚的& m 拟w 究p r o d u c i l o nw o r k s h o p w h e n t h c w i n d3 p e e d i s e q u a l l o20 m bw 州o ir c t o nr 啪嗯羽嘲牌：凳，。戴图4 5 风速为4om 旭时车川内的风场和颗粒物浓度分巾f 1 94 5w l n dn e l da n dl l l ed l s t r i b u t l o no f p a n i cu i a l ec o n c e n i r a “o nl n s i d el h ep r o d u c t i o nw o r k s h o pw h e nt h ew i n ds p c c di se q u a l t o40n s图4 - 6 n l 速为8on 以时牟f i i 】内的风场和颗 物浓度分靠f l g4 6w j n dn e l da n dt h ed l s t r i b u i j o no f p 盯 l c u l a t cc o n c en i 哺t i o nm s l d ct h cp f o d u c t l o nv 帕r k s h o pw h e l li h ew l n ds p e e di se m l a l l o80 毗甏fjll、i戮生产车间内纳米二氧化钛粉尘扩散的数值模拟研究风速m ，s毫j篓图4 7 最大浓度和平均浓度随风速变化曲线( z _ 1 6m )f 培4 7t h ec h a l l g ec u r v eo fm a x i m u mc o n c e n t r a t i o na n da v e r a g ec o n c e n t r a t i o na l o n gw i t hw i n ds p e e d ( z = 1 6m ) 比较最大浓度和平均浓度与风速的关系( 图4 7 ) 后发现，风速越大，平均浓度越小，即通风能明显降低车间内纳米二氧化钛的污染水平，但是，最大浓度并不总随着风速的加大而降低，当风速增大到4m s 和8m s 时，最大浓度反而有所增加，可能因为较大的风速加速了粉尘向排放源下风向区域堆积的趋势。从粉尘颗粒的运动轨迹( 图4 8 ) 来看，风速较低时( v = o 2 5 州s 、v = o 5 州s ) ，粉尘仍然可以扩散到车间的其他区域，随着风速的增大，粒子扩散受到流场的影响十分明显，在排放源附近随气流做一定周期的螺旋运动后被清除出车间外。2 9生产车间内纳米二氧化钛粉尘扩散的数值模拟研究图4 8 不同风速下粒子的运动轨迹f i g 4 8t h et r a j e c t o r i e so fp a r t i c l e sw i t ht h ed i f l f e r e n tw i n ds p e e di se q u a lt oz e r o 3 0，问 纳米= 机化蚀耢尘扩散的数值筷拟研究42 不同风向时车间内纳米t i 0 2 浓度分布情况物理模型、网格划分同3l 节，模拟为空气温度2 85 ，使用r n gh 两方程模型采用二阶迎风格式对控制方程进行离散，利用s i m p l e 算法对流场进行求解。风速为05 毗、4o m ，s 。05n 以属十一级风，也是人体感觉比较舒适的风速，在这里代表低风速40 m a 属于= 级风，在陆面卜可掀起大量扬半而加重空气污染，在这里代表高风速。考虑东、南、凸、北四种风向。421 低风速时粉尘的分布情况x e a s t图4 - 9 东风风速为o5m 居时车间内的风场和颗粒物浓度分布f i g4 9w i n dn e l da n dt h ed i s t r i b 。no f p a i c u l a 【ec o n c e n t r 矾i o ni n s l d e h cp r o d u c t l o nw o r k s h o pw h e nt l l cc a s tw i n ds p c c dl se q u a l i o05m s086420十p 1 j m n 米氘* 仕耪m r 数值撞m 究唔想。，。；。瓣黉j ，爹x e as l幽4 i o 南风m 述为05 m 居时车1 1 j 内的风场和颗粒物浓度分mf 增4 - l ow i n dn e l d l n dc h cd l s t r l b u t o ip n n l c u l ac cc o n c e n t r 柚n s i d cc h ep r o d u c t j o r k s h o p “h e nt h es o l i i hw m ds p e e dl s 。ql i a i 吣05 州sx 一- e a s t吲4 1 l 两j 扎肛l 述为05 m 州1 ：f i 】内的肛l 场和颗粒物浓艘分却f 咕4 一l lw l n dn e l da 川l d l s t r m u t l o n0 r 坤m c u l a t ec o n c c 胁a “0 1 1m s i d ci h ep r o m l c t i o nw o r k s h o pw h e nt h ew c “w i n ds p e e di se q 协0sm 瓜086420生p j 目自蚺米一氧化钛特尘”触教值徭拟目f 究黑尊想烹烹。，。；。辫j i 鬻：j 麓繁。黎i ：，懑j i j i 10 誉jx - ) e a s l图4 1 2 北肛l 风速为o5r r i ，s 时年问内的风场和颗粒物浓度分布i 。l g4 - 1 2w i n d6 e l da di h ed l s l b u l l o no f p a n l c u l a kc o n c e m f 砒l o ni n s i d et h ep r o d u c t i o nw o r k s h o pu m e n 山e n o n hw i n d3 p e e di s 。q u a l 如o5 u s幽4 9 4 - 1 2 是低风速时z _ l6 m 高度上纳米二氧化钛的风场和浓度场分布图巾的箭头代表水平方向的风矢量，阴影等值线区代表纳米二氧化钛的堆积浓度。总的来看，低风速不能很好地清除车间内堆积的纳米氧化钛粉尘，不州风向的情况又柏所不同，西风时( 图4 1 1 ) 的最大浓度最高，为9 35 g ，m 3 ，堆积范田也最大，北风时( 图4 1 2 ) 的堆积范围屉小，浓度较西风低，这与车删内的布局有关，南北向窗户鞍多，并且位鼹对称，排放源不位下南北窗户的连线上，容易形成穿堂风，排放源的下风向难以形成较大范围的堆积区。牟问的东西两侧通风口较少，井东西距离较长，加上排放源对西侧通风口通风的阻碍作用，东西风向时，难以形成强劲的穿堂风，敛使纳米二氧化钛粉尘扯排放源下风向大范围堆积。从平均浓度来看( 图4 1 3 ) ，o5m 愚的情况f ，西风时平均浓度最高，北风时的平均浓度屉低说明在低口l 速是，【e 风对车问内堆积的纳米二氧化钛粉乍的消洁作丌】最立f 西风的清洁作用最井。生产车间内纳米二氧化钛粉尘扩散的数值模拟研究东南风向西北图4 1 3 最大浓度和平均浓度随风向变化曲线( z = 1 6m ，v = 0 5 州s )f i g 4 - i3t h ec h a n g ec u r v eo fm a x i m u mc o n c e n t r a t i o na n da v e r a g ec o n c e n t r a t i o nw i t hd i f f e r e n tw i n dd i r e c t i o n ( z = 1 6m ，v = 0 5n 1 s ) 3 4舌；约的on乓t1雠塔产 相内蚺米一氰化钍粉尘f 瞅的敷横拟q 究图4 1 4 不同风向下粒子的运动轨迹( v = 05m 居)f 1 94 - 1 4t h e t r 面e c t o n e so f p 甜】c l e sw i t h m ed i 虢r e n t w l n dd l r e c t l o n ( v = 05 州s )4 22 高风速时粉尘的分巾情况哐圈一t 0 3 。4 0 “m e 4 m，：毫；。芍：誓量笔量譬釜= 量善罨；襄暑爹蟊摹爹x 一 e a s t| 墨| 4 - 1 5 乐风风速为40 m ，s 时车问内的皿畅和颗粒物浓度分布f i g4 1 5w i n dn e l da n dt h cd i s t r i b u t l o no f p a n l c u l a kc o n c e n t r a l l o ni n s i d et h ep r o d u c t j o nw o r k s h o pw h e nt h ee a s tw i n ds p e c di se q u a it o4or n ，sx 一 e a s t幽4 一1 6 南风肛l 述为40 m s 时乍m 内的风场和颗 t 物浓度分布f 碴4 1 6w i n dn e l da n dt h ed l s m b u l l o no f p a n l c u l a t cc o n c e n i r a t i o ni n s l d e 【h cp r o d u c 【i u nw o r k s h o p “h e ni h es o u t l 】w m ds p e e di s 。q u a l t o40m 肛p f 目m 纳= 氧化镕尘* 艘的鼓自摸m q 究唧黑烹耐，。蒜耄耋萋差囊警0 | 篙董蔫芋戛彳，j i _ 一- i 耋囊鍪嬖| ；| | 一i0：爹x 一 e a s t罔4 - 1 7 西风风速为4om 如时车| i l 内的眦场和颗粒物浓度分布f 吨4 - 1 7w j n dn e l da n dt kd l s t r l b u l i o no f p a n l c u l a c cc o n c e n f 哺l i o l li n s i d et h ep r o d u c t i o nw o r k s h o pw 、e nt h ew e scw ln ( is p e e dl so q u a lt o40m ，sx 一 e a s t罔4 - 1 8 北风风速为40 m ，s 时车问内的肛l 场和顺粒物浓度分撕f i g4 - 1 8w l n dn e l da n di h ed i s t 曲u i l o no r p a n 川l a l cc o n c e 呲r a l i o nm s l d ct h ep 巾d u d l o nw o r k s h o pw h c nm cn o n hw i n ds p c c dl se q u a l 吣40m s斟4 1 5 4 1 8 是高风速时z - l6m 高度上纳米二氧化钛的m 场和浓度场分布，棚比于05m 氐的风速，40m 居的风速使颗粒物的浓度和扩散范| 】；| 有显著减小，最大浓度减少到5 5 ”g ，m 3 ( 图4 1 9 ) 以下，污染物也只在污染源的下风向随风场有小范围的分布，东西风向时污染物帕堆捌范周比南北风向的堆积范围丈，这主要与车间的布局有关。从截面的平均浓度( 图4 1 9 ) 来看，北风的清洁效果弘08642o086420生产车间内纳米二氧化钛粉尘扩散的数值模拟研究仍然是最好的，西风的情洁效果还是最差的。从粒子的运动轨迹( 图4 2 0 ) 来看，用于风场的作用，粒子在车间的项部被气流吹起，做螺旋运动，直至被风带出车间。c oj越螺东南西北风向图4 - 1 9 最大浓度和平均浓度随风向变化曲线( z = 1 6m ，v = 4 0 州s )f i g 4 - l9t h ec h a n g ec u r v eo fm a x i m u mc o n c e n t r a t i o na n da v e r a g ec o n c e n t r a t i o nw i t hd i f f e r e n tw i n dd i r e c t i o n ( z = 1 6m ，v = 4 om s ) 3 7印如 寻加。生产车间内纳米二氧化钛粉尘扩散的数值模拟研究图4 2 0 不同风向下粒子的运动轨迹( v = 4 0 耐s )f i g 4 - 2 0t h et r a j e c t o r i e so fp a r t i c l e sw i t ht h ed i f i f e r e n tw i n dd i r e c t i o n ( v = 4 。om s )4 3 车间内纳米t i 0 2 浓度的问题中华人民共和国国家标准室内空气质量标准中规定，室内可吸入颗粒物日均值不能超过o 1 5m m 3 ( 1 5 0p m 3 ) ，与世界卫生组织健康住宅标准的标准相同。从模拟和观测结果来看，通风时，车间内纳米t i 0 2 浓度值都低于1 5 0 州m 3这一数值，但是由于纳米二氧化钛粉尘的粒径很小，长期在纳米二氧化钛生产车间工作的工人和其他办公人员是否安全，还需要进一步论证。在不通风的情况下，车间内纳米二氧化钛粉尘的浓度很容易超过1 5 0 肛咖3 这一限值，因此应该经常开门开窗通风，使车间内纳米二氧化钛粉尘得到及时有效的清除。3 8生产车间内纳米二氧化钛粉尘扩散的数值模拟研究4 4 小结通风对纳米材料生产车间内颗粒物的浓度和扩散范围有显著的影响。较大的风速能显著降低粉尘的浓度，缩小颗粒物的扩散范围，对纳米材料生产车间有很好的清洁作用。车间内粉尘泄漏点( 排放源) 的位置位于两侧通风口的中间时，污染物容易大范围堆积。但较大的风速也使从业人员感到不适甚至使他们的健康受到损害，尤其是在秋冬季节，同时较大的风速也加速了超细颗粒物向外界的排放。因此，改进设备和生产工艺、尽可能减少车间内颗粒物的泄露才是根本。3 9生产车间内纳米二氧化钛粉尘扩散的数值模拟研究第五章结论与展望车间内纳米二氧化钛粉尘的扩散与堆积问题应当引起足够的重视，因为粉体不仅容易引起爆炸，还会对身处其中的工作人员造成健康的危害，本文用计算流体力学f l u e n t 软件研究了车间内粉尘扩散的问题，主要完成了以下工作：( 1 ) 初步建立起室内纳米粉尘的扩散模型，并与现场采样的结果进行对比，对比显示：模型值与采样值比较吻合，粉尘的浓度量级与扩散趋势较为一致。( 2 ) 系统分析了不同通风条件下车间内纳米二氧化钛的扩散问题，发现通风虽能显著降低室内粉尘的浓度，但是风速较低时，在排放源的下风向仍有大范围污染区存在。但在研究中还是存在一些问题：a ，由于时间的限制，论文未就不同湿度、温度、气压等空气条件对粉尘扩散的影响进一步展开分析研究。论文是以相关的纳米课题丽展开的，课题的计划是分析车间内纳米粉尘在不同的环境空气条件的浓度分布，进而得到车间内纳米粉尘扩散的日变化特征和年变化特征，但是这个项目还在开展中，所以这部分内容还不能表现在论文中。b ，车间内纳米粉尘可能对从业人员的人身安全和健康造成威胁，但是粉尘的浓度值多少才是安全，却需要医学科研人员对粉尘从业人员长期跟踪对比分析给出，也需要防爆评估行业的专家进行评定。c ，模型的调试过程中，我们利用车间的现场采样结果对模型进行了不断的修正，但是修正后的模型是否同样适用于其他类型的生产车间，还需要更进一步的验证。生产车间内纳米二氧化钛粉尘扩散的数值模拟研究参考文献【l 】c h e o n gk ，e d 口u n 扯d y t p o h ，e ta 1 m e a s u r e m e n t sa n dc o m p u t a t i o n so fc o n t a m i n a m sd i s 仃i b u t i o ni n 锄o m c ee n v i r o n m e n t b u i l d i n g 锄de n v i r o 啪e n t ，2 0 0 3 ，3 8 ：1 3 5 1 4 5【2 】h a n o gj ，gh o e ka p e t e r s ，e ta 1 e 髓c t so ff i n ea n du l 仃a f i n ep a n i c l e so nc a r d i o r e s p i r a t o d rs y m p t o m si ne l d e r l ys u b j e c t sw i t hc o r o n a 叫h e a nd i s e a l s e ：t h eu l t r as t u d y a m e r i c a nj o u m a lo fe p i d e m i o l o g 弘2 0 0 3 ，15 7 ( 7 ) ：613 6 2 3【3 】i b a l d - m u l i ia ，h w i c h m a n n ，w k r e y l i n g e p i d e m i o l o g i c a le v i d e n c eo nh e a l t he 毹c t so fu l t r a f i n ep a n c l e s j o u m a lo f a e r o s o lm e d i c i n e ，2 0 0 2 ，l5 ( 2 ) ：18 9 2 01【4 】k h a nj ，c f e i g l e y ，e l e e ，e ta 1 e f f e c t so fi n l e ta n de x h a u s tl o c a t i o n sa n de m i t t e dg 够d e n s i t ) ，o ni n d o o ra i rc o n t a m i n a n tc o n c e n t r a t i o n s b u i l d i n ga n de n v i r o n m e n t ，2 0 0 6 ，41 ：8 5 l 8 6 3【5 】l iy ，z c h e n ab a l a n c e - p o i n tm e 也o df o ra s s e s s i n gt h e e f r e c to fn a t u r a lv e n t i l a t i o no ni n d o o rp a r t i c l ec o n c e n t r a t i o n s a t m o s p h e r i ce n v i r o n m e n t ，2 0 0 3 ，3 7 ：4 2 7 7 4 2 8 5【6 】o b e r d 6 r s t e rgp u l m o n a r ye f r e c t so fi n h a l e du l t r a f i n e i n ta r c ho c c u pe n v i r o nh e a l t h ，2 0 0l ，7 4 ( 1 ) ：1 8 7 】o d a rf v r e r i n c a t i o no fp r o p o s e de q u a t i o nf o rc a l c u i a t i o no ff o r c e so nas p h e r ea c c e l e r a t i n gi nav i c o u sn u i d j f l u i dm e c h 19 6 6 ，2 5 ( 3 ) ：5 9 1 5 9 2 8 】r i l e yw ，a r o b i n s o n ，a g a d g i l ，e ta 1 e 疏c t so fv a r i a b l ew i n ds p e e da n dd i r e c t i o no nr a d o nt r a n s p o r tf r o ms o i li n t ob u i l d i n g s ：m o d e ld e v e l o p m e n ta n de x p l o r a t o r yr e s u l t s a t m o s p h e r i ce n v i r o n m e n t ，1 9 9 9 ，3 3 ：2 1 5 7 2 1 6 9【9 】s p a l d i n gd ag e n e r a lp u r p o s ec o m p u t c rp r o g r a n lf o rm u l t i d i m e n s i o n a lo n e - a n dt w o p h a s ef l o w m a t h e m a t i c sa n dc o m p u t e r si ns i m u l a t i o n ，1 9 8 1 ，2 3 ( 3 ) ：2 6 7 2 7 6【10 】t h o m a sk ，k y u e f f e c t so fa rc o n d i t i o n i n g ，d e h u m i d i f i c a t i o na n dn a t u r a lv e n t i l a t i o no ni n d o o rc o n c e n t r a t i o n so f2 2 2 r na n d2 2 0 r n j o u r n a lo fe n v i r o n m e n t a ir a d i o a c t i v i t y ，2 0 0 0 ，4 7 ：1 8 9 1 9 9【1l 】t i a nz ，j t u ，gy e o h o nt h en u m e r i c a ls t u d yo fc o n t a m i n a n tp a r t i c l ec o n c e n t r a t i o ni ni n d o o ra i r f l o 、b u i l d i n ga n de n v i r o n m e n t ，2 0 0 6 ，4l ：15 0 4 l5l4【12 】v e r s t e e gh ，w m a l a l a s e k e r a a ni n t r o d u e t i o nt oc o m p u t a t i o n a lf l u i dd y n a m i c s t h ef i n i i ev o j u m em e t h o d 。u s a ：a d d j s o n w e s l e yp u bc o ，j9 9 6 4 7 4 8【l3 】w a n gl ，t b u r t o n ，d s t o e k c h a o t i cd y n a m i c so fh e a v yp a n i c l ed i s p e r s i o n ：f r a c t a ld i m e n s i o nv e r s u sd i s p e r s i o nc o e 衔e i e n t s p h y s ，f l i u d sa ，2 ，1 9 9 0 ：1 3 0 5 1 3 0 8【l4 】z h a ob ，c y a n g ，x y a n g ，e ta 1 p a n i c l ed i s p e r s i o na n dd e p o s i t i o ni nv e n t i i a t e dr o o m s ：t e s t i n ga n de v a l u a t i o no fd i 仃e r e n te u l e r i a na n dl a g r a n g i a nm o d e l s b u i l d i n ga n de n v i r o n m e n t ，2 0 0 8 ，4 3 ：3 8 8 3 9 7【15 】z h a ob ，j w u p a n i c l ed e p o s i t i o ni ni n d o o re n v i r o n m e n t s ：a n a l y s i so fi n n u e n c i n gf 2 比t o r s j o u m a lo f h a z a r d o u sm a t e r i a l s ，2 0 0 7 ，1 4 7 ：4 3 9 4 4 8 1 6 】程易气固两相流动数值模拟及其非线性动力学分析： 尊士学位论文】北京：清华大学4 1生产车问内纳米二氧化钛粉尘扩散的数值模拟研究化工系2 0 0 0【1 7 】丛晓春，张旭等开放性尘源粉尘运动轨迹的数值计算同济大学学报( 自然科学版) ，2 0 0 5 ，3 3 ( 3 ) ：3 3 0 3 3 3【1 8 】樊建人，罗坤，金晗辉等宣接数值模拟三维气固两相混合层中颗粒与流体的双向耦合中国电机工程学报，2 0 0 3 ，2 3 ( 4 ) ：15 3 l5 7【1 9 】方丁酉两相流动力学长沙：国防科技大学出版社，1 9 8 8 4 2 4 4 2 6 2 0 】方坤计算流体力学的几种常用软件煤炭技术，2 0 0 6 ，2 5 ( 12 ) ：l1 1 2【2 1 】韩占忠，王敬，兰小平f l u e n t 流体工程仿真计算实例与应用北京：北京理工大学出版社，2 0 0 4 2 2 2 5【2 2 】黄社华、李炜、程良骏任意流场中稀疏颗粒运动方程及其性质应用数学和力学，2 0 0 0 ，2 l ( 3 ) ：2 6 5 2 7 6 2 3 】黄优，桂祥友，梦絮屹等筛分车间内粉尘颗粒在气流中的受力分析及其运动轨迹的研究煤炭工程，2 0 0 8 ，8 ：8 8 9 0 【2 4 】李爱民，姜秀明颗粒对边界层及对流换热影响机理的探讨华北电力大学学报1 9 9 9 ，2 6 ( 1 ) ：4 6 5l【2 5 】李洪钟，郭