基于响应曲面法的纤维过滤性能预测模型.docxVIP

基于响应曲面法的纤维过滤性能预测模型 人类平均超过80%的时间都在室内。根据统计数据，室内污染通常超过室内2.5倍，有时超过100倍。室内污染种类繁多，来源不同，对人类的危害十分突出。如何有效控制室内污染，改善室内质量，是一个必须尽快解决的问题。因此，空气纤维过滤装置成立，并迅速发展。空气过滤理论的研究越来越深入。首次实验是在德国进行的。德国用纤维垫开始了气溶胶过滤机制的理论探讨。从那时起，许多国际科学家已经先后研究了纤维过滤的性能。许多科学家使用不同的影响因素来研究过滤器的性能指标。jagnathan分析了纤维直径（或多峰）分布对纤维过滤装置性能的影响，chen研究了填充密度的变化对纤维过滤装置性能的影响，sureh和smithhorsh分析了不同形状纤维对纤维过滤装置性能的影响。 区别于以上学者的研究,本文通过数值模拟的方法研究纤维过滤器内部结构和运行参数对纤维过滤器性能的影响. 改变其结构参数和运行条件往往重复计算很多次,如果不对计算次数进行优化设计,其工作量是很大的. 而响应曲面法(Response surface methodology)恰好能够解决正交试验所带来的麻烦,并得到了广泛应用,本文应用响应曲面法对计算次数进行设计,并通过Minitab软件拟合出纤维过滤器性能与纤维过滤器内部结构参数和运行条件之间的关系式,从而找出最佳的纤维过滤器内部结构排列,得到高效、低阻的纤维过滤器结构. 1 数值模型的构建 1.1 估计方程的建立 响应曲面设计是在试验设计中使用的一种数据处理方法,根据其工作原理,首先考虑选择响应曲面方程的形式. 本文选择二次响应曲面方程,并考虑所有的一次项,二次项和两两交叉项,响应曲面方程可表示为 Y=β0+∑i=1kβiXi+∑i=1kβiiX2i+∑i?jkβijXiXj+Y=β0+∑i=1kβiXi+∑i=1kβiiXi2+∑i?jkβijXiXj+ e(X1,X2,…,Xk) (1) 式中,Y为目标函数或响应;Xi为自变量;βi,βii,βij代表一次、二次交互作用项的回归系数;k为影响因素的数量;e为误差,有两个误差源,即试验误差和拟合不足误差,后者包括高阶项和交互项. 研究纤维过滤器性能时考虑的参数:a/dl,b/dl,df/dl,Re.a为纤维之间的横向距离,b为纤维之间的纵向距离,df为纤维直径,dl为单根纤维长度,雷诺数Re=νdl/η. 一种估计方程中各参数的方法是来研究在3个不同水平下所有因素的响应特性,这种全部析因设计将需43=64种不同的CFD模拟. 从时间和资源角度考虑,全部析因设计是不可取的. 本文所采用的4因素3水平面中心的中心复合设计仅需25次CFD模拟,和3水平的全部析因设计数量相比大大减少. 式(1)的回归系数由最小二乘法拟合得到. 为方便起见,将所有变量按下式规格化: xi=Xi?Xˉˉˉi1/2(XiH?XiL)(2)xi=Xi-Xˉi1/2(XiΗ-XiL)(2) 式中,XiH和xiL分别为第i个参数的最大值和最小值;XˉˉˉXˉi为Xi的平均值.表1给出了响应曲面中变量的范围. 1.2 压力出口边界 本文所研究的纤维过滤器的模型尺寸及其边界条件如图1所示. 假定空气以速度进口进入计算区域,以压力出口边界条件离开计算区域,对于风管边界,本文采用了对称边界条件,纤维表面处采用无滑移边界条件,过滤器厚度为h,单根纤维的长度dl取200 μm,出口距过滤器距离L′取3df,进口距过滤器距离为L,计算时为了保证气流在过滤器管道内充分发展,文中L取10df. 1.3 dp/df1模型 本模型采用稳态、层流,不可压缩模型,其控制方程如下: ?ρ?t+??(ρv)=0(3)?ρ?t+??(ρv)=0(3) ρDvDt=??p?η?×(?×v)+4/3η(??v)(4)ρDvDt=-?p-η?×(?×v)+4/3η(??v)(4) ρcpDTDt=?DpDt+Φ+k?2T(5)ρcpDΤDt=-DpDt+Φ+k?2Τ(5) 式中,η为流体的动力粘度;cp为流体比热;ρ为流体密度;v为速度;p为压力. 本模型为稀相,不考虑颗粒相对气相场的作用. 在常温常压条件下,当颗粒直径dp和纤维直径df量级相当时,才在过滤过程中考虑拦截作用. 当dp/df?1时,拦截作用可以忽略不计. 本文中,dp/df≤1,忽略拦截作用对过滤效率的影响. 颗粒运动方程为 dvipdt=Fd(vi?vip)+Fbi(6)dvipdt=Fd(vi-vip)+Fbi(6) 式中,vi和vip为x、y或z方向上气相和颗粒相速度;Fd为曳力;Fbi为布朗力. 在颗粒雷诺数Rep(Rep=ρvdp/η)小于1时,作用在颗粒上的曳力和布朗力分别为 Fd=18ηd2pρpCc(7)Fd=18ηdp2ρpC