Three-dimensional reconstruction and filtration behavior analysis system for nanofiber membrane

纳米纤维膜三维重构与过滤行为分析系统 Three-dimensional reconstruction and filtration behavior analysis system for nanofiber membrane 1 2 first-author 郑高峰（1984－），男，福建泉州人，博士，副教授，2006年于武汉科技大学获得学士学位，2011年于厦门大学获得博士学位，主要从事静电纺丝技术和有机微纳米系统制造领域的研究。E-mail：zheng_gf@xmu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425979&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425955&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425981&type= 22.01333237 32.00399780 郑高峰 （1984－），男，福建泉州人，博士，副教授，2006年于武汉科技大学获得学士学位，2011年于厦门大学获得博士学位，主要从事静电纺丝技术和有机微纳米系统制造领域的研究。E-mail： zheng_gf@xmu.edu.cn zheng_gf@xmu.edu.cn 1 2 1 2 1 2 1 2 corresp E-mail： yfliu@xmu.edu.cn E-mail： yfliu@xmu.edu.cn 刘益芳（1976－），女，福建莆田人，博士，副教授，1999年、2002年于东南大学分别获得学士、硕士学位，2010年于厦门大学获得博士学位，主要从事MEMS传感器、微硅陀螺仪的控制与测试和微弱信号检测等方面的研究。E-mail：yfliu@xmu.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425957&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425990&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425989&type= 22.01332855 32.00399780 刘益芳 （1976－），女，福建莆田人，博士，副教授，1999年、2002年于东南大学分别获得学士、硕士学位，2010年于厦门大学获得博士学位，主要从事MEMS传感器、微硅陀螺仪的控制与测试和微弱信号检测等方面的研究。E-mail： yfliu@xmu.edu.cn yfliu@xmu.edu.cn 厦门大学 仪器与电气系，福建 厦门 361005 厦门大学 仪器与电气系，福建 厦门 361005 Department of Instrumental and Electrical Engineering， Xiamen University， Xiamen 361105， China Department of Instrumental and Electrical Engineering， Xiamen University， Xiamen 361105， China 福建省增材制造创新中心，福建 福州 350118 福建省增材制造创新中心，福建 福州 350118 Fujian Innovation Center of Additive Manufacturing， Fuzhou 350118， China Fujian Innovation Center of Additive Manufacturing， Fuzhou 350118， China 为了对纳米纤维空气过滤膜进行有效的理论数值计算和分析，建立了纳米纤维膜空气过滤性能计算的分析系统。对该系统所采用的三维重构、过滤效率计算及过滤压降计算等算法进行研究。通过去噪处理、二值化、细化算法和离散化进行纳米纤维膜的特征提取，并通过对纳米纤维膜的特征信息定义、存储和可视化渲染实现三维重构。通过对颗粒物与纳米纤维的碰撞检测和力学分析判断颗粒物被纳米纤维膜捕获的数量进行过滤效率的计算。最后，根据泊肃叶定律和压降原理实现了基于SEM图像的过滤压降计算模块。针对纳米纤维膜的拓扑结构和厚度进行了数值计算与实验比对，对比结果表明，过滤效率计算误差在10%以内，过滤压降计算误差在20%以内，计算结果趋势与实验结果相符，能够反映不同拓扑结构的过滤性能差异。 To realize effective theoretical calculation and analyses of a nanofiber air filtration membrane， an analysis system was developed to calculate the air filtration performance of the nanofiber membrane. The 3D topology reconstruction algorithm， the calculation module of the filtration efficiency， and the calculation method of pressure drops were investigated. First， feature extraction of the nanofiber membrane was conducted using denoising， binarization， thinning algorithm， and discretization， and the three-dimensional reconstruction was realized through the definition， storage， and visual rendering of the feature information of the nanofiber membrane. Subsequently， through collision detection and mechanical analyses of particles and nanofibers， the number of particles captured by the nanofiber membrane was determined to calculate the filtration efficiency. Finally， according to Poiseuille's law and the pressure drop principle， a filtration pressure drop calculation module based on SEM images was realized. The numerical results of both the topological structure and thickness of the nanofiber membrane were compared with experimentally determined values. The comparison results show that the calculation error of filtration efficiency is less than 10%， and the calculation error of filtration pressure drop is less than 20%. The trend of calculation results is consistent with the experimental results， which reflects the difference in the filtration performance of different topologies. 纳米纤维膜 三维重构 分析系统 过滤效率 过滤压降 nanofiber membrane 3D reconstruction analysis system filtration efficiency filtration pressure drop 1 引　言 1 随着工业经济的高速发展，大气污染日益严重，气候日益复杂，空气中的颗粒物特别是细粒子给人们的身体健康带来了严重的危害，对城市的能见度产生了巨大的影响 ［ 1 1 ］ ，空气过滤技术的需求和发展极其迫切。传统的无纺布、棉絮、玻璃纤维等过滤器难以有效过滤PM 2.5 ，纳米纤维具有比表面积大、吸附能力强等特点，在空气过滤领域中表现出了强大的应用潜力。 由直径小和比表面积大的纳米纤维 ［ 2 2 - 3 3 2-3 ］ 堆叠而成的滤膜具有孔径小、孔隙率高、密度低及通透性好等特点，科研人员对纳米纤维膜的制造方法、过滤机理和过滤性能等做了大量的研究。然而，当前研究主要通过实验的方式进行，纳米量级的不确定性增加了进一步探索的难度。为了深入理解纳米纤维膜的空气动力学行为，科研人员构建了相关数学模型对过滤机理进行分析，研究包括空气过滤过程粒子行为、纳米纤维膜空气过滤流场和纳米纤维膜结构优化的数值模拟。 在空气过滤过程粒子行为的数值研究过程中，Cheung等 ［ 4 4 ］ 采用半经验公式法模拟过滤过程中驻极纤维对颗粒物动力学行为的影响；Li等 ［ 5 5 ］ 研究颗粒物在湍流通道下的扩散与沉积行为，用数值表示了流体阻力、布朗效应和重力效应；Cleaver等 ［ 6 6 ］ 建立了湍流中颗粒物的惯性和扩散行为；Wang等 ［ 7 7 ］ 在格子玻尔兹曼方法的基础上进行优化，建立了空气颗粒物在单根纤维上的沉积行为。在纳米纤维膜空气过滤流场数值研究中，Wu等 ［ 8 8 ］ 利用计算机仿真技术模拟了纳米纤维膜的过滤过程，探究了不同纳米纤维直径对空气流场的影响；Dong等 ［ 9 9 ］ 使用CFD-DEM对颗粒物在纤维表面的沉积行为进行了数值研究；Pan等 ［ 10 10 ］ 利用纤维特性建立了逼真的纤维3D模型，通过仿真模拟发现当纤维直径小于5 μm时，纤维直径的变化对滤膜性能有很大影响，而当纤维直径大于12 μm时，纤维直径对于性能的影响会大幅减小。在纳米纤维滤膜结构优化方面，Zhao等 ［ 11 11 ］ 利用PVDF与负离子粉末进行混合制备出了能够释放负氧离子的纳米纤维，并通过数值计算对纳米纤维膜进行了优化；Sambaer等 ［ 12 12 ］ 为了探究过滤超细颗粒物（20~400 nm）的纳米纤维膜的结构特点，提出了一种数值计算模型，该计算模型能够有效预测出纳米纤维膜结构的过滤性能；Zhao等 ［ 13 13 ］ 研制了一种高透光率的纳米纤维膜，并借助计算机仿真技术探究纳米纤维膜结构中纤维直径对滑移效应的影响。 深层过滤是纳米纤维膜的特点与优势。纳米纤维膜的立体结构特征，如纳米纤维直径分布 ［ 14 14 ］ 、排列方式 ［ 15 15 ］ 、层叠结构 ［ 16 16 ］ ，是空气过滤品质的重要决定因素。但现有研究仍采用随机生成的方法进行重构，无法真实反映纳米纤维膜的结构特征，难以实现过滤性能的精确计算。另外，仅依靠单一纤维碰撞或捕获效应的计算方法，也限制了纳米纤维膜过滤仿真计算准确性的提升，不利于理论验证或为过滤膜结构优化提供有效的理论指导。因此，针对纳米纤维膜空气过滤构建专门的集成化仿真分析系统，成为纳米纤维过滤机制探究与高效滤膜结构优化生产的研究重点。 本文搭建了一个可用于纳米纤维膜空气过滤的仿真分析系统，该分析系统集成了纳米纤维膜重构、过滤效率计算和压降计算等功能，能够有效计算出不同纳米纤维膜的过滤性能，从而为科研人员提供便捷的纳米纤维膜空气过滤的分析工具。 2 纳米纤维膜空气过滤分析系统 1 2.1　分析系统的总体方案设计 2 在实现计算颗粒物过滤效率和计算过滤压降两大基本功能的基础上，系统要具有以下功能：能够构建纳米纤维膜的形貌，而非通过随机堆叠或手动创建的简易的等直径纳米纤维空气过滤膜模型；整个过滤过程可视化；可增加新功能模块。针对这些要求，本文搭建的仿真分析系统是一个能够通过SEM图像提取纳米纤维膜的拓扑结构并进行三维重构，通过三维重构的纳米纤维膜模型计算纳米纤维膜空气过滤的过滤效率和压降。 纳米纤维膜空气过滤仿真分析系统的整体设计方案如 图1 图1 所示。将整个系统分为可视化渲染和数值计算两个模块，每个模块采用独立线程进行处理。多线程处理能够有效解决可视化渲染与数值计算队列运算造成的时间浪费问题。在多线程的基础上，渲染模块和数值计算模块存在着处理速率不匹配的问题，通常来说渲染模块的效率要低于数值计算模块。为了解决速率不匹配问题，采用线程间的带锁缓存队列。在渲染模块与计算模块的交互中，数值计算模块将计算完成的颗粒物位置存入线程安全的缓存队列中，相应的渲染模块从缓存队列取出颗粒物更新数据，形成即存即用的机制。这种方式能够有效降低所需的计算存储资源，降低所搭建系统的使用要求。 可视化渲染模块分为滤膜渲染、颗粒物渲染和视角变换3个小模块。数值计算模块分为压降计算模块和过滤效率计算模块。 OpenGL作为最底层与GPU进行交互的图形标准，它独立于环境易于移植并且具备功能强大、效率高等特点，所以将它作为可视化渲染模块的主要工具。在图像处理部分，选取的工具是OpenCV，这也是当前广泛使用的开源图像处理函数库。 在整个程序中，关键模块主要有纳米纤维膜的三维重构、空气过滤效率计算和过滤压降计算。 2.2　纳米纤维膜三维重构模块 2 2 . 2 . 1　纳米纤维膜的特征提取 3 纳米纤维膜三维重构模块的功能是从SEM图中经过图像前处理提取出纳米纤维膜的拓扑信息和纳米纤维的表面拓扑形貌，并通过相应方法建立几何关系，最后通过计算机渲染进行可视化。为了能够构建逼真的纳米纤维膜模型，首先进行纳米纤维膜的特征提取。Sambaer等 ［ 17 17 ］ 提出了一种基于SEM图像的逼真纳米纤维膜构建方法，该方法主要有去噪处理、二值化、细化算法和离散化处理4个步骤。去噪处理用于去除原生SEM图像的噪点或模糊纳米纤维表面形貌的细节；二值化用于提取顶层的纳米纤维层膜，其原理在于SEM图像的灰度值与深度信息存在联系；细化算法用于提取出纳米纤维膜的中心线；离散化的作用是将连续纤维离散化成单位元素体的组合，通过半径提取算法赋予每个元素体对应的改变量，以还原纳米纤维的表面拓扑形貌。 去噪处理模块的流程为：首先，将原生SEM图像通过cvtColor函数转换为灰度图像；其次，将SEM灰度图像利用GaussianBlur函数进行高斯模糊，结果如 图2 图2 所示。 二值处理的功能在于将高斯模糊处理后的SEM图像进行“前景”与“后景”的分离。其中，前景是指处于观测方向最近的纤维图像，后景是指处于观测方向较远的纤维图像。首先，遍历图像的像素点计算出SEM图像的灰度值的均值（ M ）及其方差（ V ）。然后，计算顶层纤维所在深度的灰度阈值，即： <math display="block" id="M1"><mi>T</mi><mo>=</mo><mi>M</mi><mo>+</mo><mfrac><mrow><mi>V</mi></mrow><mrow><mi>n</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425999&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28425998&type= 19.38866615 7.95866632 . 最后，通过threshold函数对原始SEM图像进行二值化处理，得到的顶层纤维二值图像如 图3 图3 所示。其中，白色像素部分表示为处于顶层的纳米纤维，而黑色像素的部分为孔径部分。 中心线提取模块采用的是Zhang等 ［ 18 18 ］ 的快速并行细化算法，得到的纳米纤维膜中心线图像如 图4 图4 所示。 离散化处理是将顶层连续纤维离散化，并且获取纳米纤维的拓扑形貌信息以及保存三维重构的必要信息。离散化处理方法 ［ 19 19 ］ 如 图5 图5 所示，将纳米纤维离散成以中心线上的像素点为球心的不同大小的球体组合。通过半径检测算法检测所在像素点的纳米纤维半径，用以决定当前像素点的填充球体的半径。半径检测算法原理如 图6 图6 所示，当检索到中心线像素点时，通过截取感兴趣区域（Region of Interest， ROI）执行半径检测算法。 2 . 2 . 2　纳米纤维膜三维重构程序 3 纳米纤维膜的三维重构程序作为特征信息提取程序的下游处理程序，负责上游程序所提取数据的有效处理，主要分为纳米纤维膜的三维信息定义、三维信息的有效存储和纳米纤维膜的可视化三个方面。三维重构程序的逻辑框图如 图7 图7 所示。首先，将上游程序即特征信息提取程序的反馈信息进行存储，存储结构的名称定义为体素点。其次，将体素点重组，构建成八叉树结构 ［ 20 20 ］ 。最后，对重构的纳米纤维膜利用OpenGL开源图形库进行可视化渲染。 2.3　纳米纤维膜过滤效率计算 2 纳米纤维膜过滤效率计算程序包括颗粒物动力学行为计算及其可视化模块、颗粒物与纳米纤维碰撞检测以及力学处理模块、和过滤效率计算模块。 颗粒物动力学行为的计算及其可视化模块的功能是计算并动态显示空气中颗粒物的位置。该功能由两个线程协同完成，主线程用于实时更新渲染颗粒物并对不同的设备输入进行处理，子线程负责计算颗粒物在拖曳力和布朗力的作用下的空间位置。 颗粒物与纳米纤维碰撞检测和力学分析模块的功能是实时检测颗粒物是否与纳米纤维发生碰撞，并在碰撞发生时进行力学分析，判断颗粒物是否被纳米纤维捕获。根据当前颗粒物与所检测的填充球体的球心距离是否小于两者半径之和判断是否碰撞。若没有发生碰撞则将颗粒物的位置放入缓存队列中；当发生碰撞时，计算出摩擦力、升力等作用力的大小并判断颗粒物是否被捕获。若判定为捕获，则将该颗粒物标记，并结束渲染；当判定为非捕获状态时，计算其绕行轨迹，并将其放入缓存队列中。纳米纤维特征尺度小、电荷密度高，静电吸附对空气过滤起到更重要的作用 ［ 21 21 ］ 。但纳米纤维膜上的电荷分布暂时难以定量化描述，这限制了静电效应在过滤仿真中的准确表达，以及系统计算误差的减小。 过滤效率计算模块的功能是当设定条件满足时计算出颗粒物被标记的数量时计算出纳米纤维膜的过滤效率。该模块检测是否有颗粒物的位置超出所假定的通风管，若颗粒物超出管道，则把该颗粒物标记为已透射并在后续的渲染中不再出现。当标记为捕获和透射的颗粒物达到初始化颗粒物的数量时，则纳米纤维膜的过滤效率为： <math display="block" id="M2"><mi>E</mi><mo>=</mo><mfrac><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><mo>-</mo><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi></mrow></msub></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426016&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426021&type= 23.11400032 8.80533314 ， 其中： <math id="M3"><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426025&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426023&type= 4.48733330 3.80999994 表示进入模拟管道的颗粒物数量， <math id="M4"><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426027&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426033&type= 5.58799982 3.80999994 表示离开模拟管道的颗粒物数量。 图8 图8 为过滤效率计算模块的可视化界面。对颗粒物进行轨迹追踪，可以观察数值计算模拟的纳米纤维膜的空气过滤过程。由于测试的颗粒物繁多，所以 图8 图8 中仅对100个颗粒物进行追踪作为示意。图中箭头1所指区域是模拟测试管道的入口，在该区域中大量的颗粒物在跟随气流运动的同时做无规则布朗运动。箭头2所指区域为三维重构的纳米纤维膜，当气溶胶颗粒物与纳米纤维膜发生碰撞后会被捕获，使得随气流运动的颗粒物数量迅速减少。因此，只有少数颗粒物能够穿过纳米纤维膜达到箭头3所指的区域，即模拟测试管道的出口。 2.4　纳米纤维膜过滤压降计算 2 在空气过滤的过程中，气体流过纳米纤维膜时的流动状态属于层流状态，并且有序地依次通过纳米纤维所构建的孔径通道穿过纳米纤维滤膜，即气体实际上是通过变截面积的通道穿过纳米纤维滤膜的 ［ 13 13 ］ 。根据泊肃叶定律可以得到流体在圆形管道中层流时的压降计算公式为： <math display="block" id="M5"><mi mathvariant="normal">Δ</mi><mi>p</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">8</mn><mi>Q</mi><mi>η</mi><mi>L</mi></mrow><mrow><mi mathvariant="normal">π</mi><msup><mrow><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426041&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426039&type= 18.79599953 9.05933285 ， 其中： <math id="M6"><mi>Q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426046&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426044&type= 2.87866688 3.21733332 表示气体的体积流量， <math id="M7"><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426049&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426063&type= 3.72533321 3.72533321 表示圆心管道的半径。 Sambaer等通过对SEM图像的处理构建出欧式距离图 ［ 12 12 ］ ，如 图9 图9 所示。图中，黑色部分代表纳米纤维膜，亮度值不为零的部分为孔径，其值与像素点到最近纤维的距离有关，距离越远，其值越大，亮度越高。根据 式（4） 式（4） ： <math display="block" id="M8"><msub><mrow><mi>R</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo>=</mo><mroot><mrow><msup><mrow><mfrac><mrow><msup><mrow><mi>s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><mi>n</mi></mrow></munderover></mstyle><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>d</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><msub><mrow><mi>d</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow></mrow><mrow><mi mathvariant="normal">π</mi></mrow></mfrac></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426066&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426050&type= 31.15733528 16.17133331 ， 其中： <math id="M9"><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426069&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426068&type= 1.26999998 2.79399991 表示像素点的实际长度， <math id="M10"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426054&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426052&type= 1.86266661 2.79399991 表示灰度值不为零的像素点个数， <math id="M11"><msub><mrow><mi>d</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426058&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426056&type= 2.96333337 3.72533321 表示当前像素点的值， <math id="M12"><msub><mrow><mi>d</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426073&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28426071&type= 5.75733376 3.72533321 表示距离最远的像素点的值。计算出单个片层的等效半径，最后根据 式（3） 式（3） 得出单个片层的过滤压降。总压降为多个片层的叠加。 3 纳米纤维膜空气过滤分析系统计算结果与验证 1 3.1　不同特征纳米纤维膜的空气过滤效率计算 2 为了验证本系统的功能，本文依据实验制备的纳米纤维膜对样品进行SEM观测，结果如 图10 图10 所示。为了更好地提取纳米纤维膜的特征，需要动态选择合适倍数的SEM图像，在纤维拓扑特征提取时只需动态地调整阈值与像素点所表达的空间尺度即可。当确定SEM图像后，为了降低计算资源的需求，减少数值计算负担，并加快速度，将SEM图像进行切分，其范围只需能够反映纳米纤维拓扑结构即可。纳米纤维膜的重构规则如下：以一张SEM图作为基准，对比实验与数值计算的结果，确定其相应图像处理参数以及重构纳米纤维膜的层数，校准完毕后其余数值计算以同样条件进行。在本次计算中，所涉及的图像处理参数为3×3的图像切分规格和mean+variance×3/10的顶层纤维提取阈值。数值计算条件如下：温度设定为298 K；气压设定为101.325 kPa；空气动力黏度设定为1.79×10 - 5 Pa·s；空气密度为1.293 kg/m 3 ；根据自动过滤效率测试仪的参数可以确定空气流速为0.15 m/s。 PVDF/PAN纳米纤维不同比例纳米纤维膜的数值计算结果如 图11 图11 所示，最大误差为8%。计算结果表明，随着PVDF纳米纤维比例的减少，PAN纳米纤维的比例增加，纳米纤维膜的空气过滤效率逐渐上升，其趋势与实验数据相符。 对于纳米纤维膜的厚度计算，选取PVDF/PAN=2∶1的纳米纤维进行数值计算，通过增减重构纳米纤维膜的层数控制厚度的变化，计算结果如 图12 图12 所示。计算结果表明，随着纳米纤维膜厚度的增加，过滤效率升高，该趋势与实验结果相吻合。 3.2　不同特征纳米纤维膜的空气过滤压降计算 2 对纤维直径与过滤压降关系和纳米纤维膜厚度与过滤压降关系分别进行数值计算，用以验证过滤压降模块的功能，实验条件及数值计算条件均与3.1节相同。 对于过滤压降，PVDF/PAN纳米纤维不同比例纳米纤维膜的数值计算结果如 图13 图13 所示。计算结果显示，随着PAN纳米纤维比例的增加，即纤维直径的减小，纳米纤维膜的过滤压降增加，该现象与实验结果相匹配。 对纳米纤维膜厚度进行数值计算，其结果如 图14 图14 所示。计算结果表明，随着膜厚的增加，过滤压降增加，该结果也与实验相符。 4 结　论 1 本文搭建了可用于纳米纤维膜空气过滤的分析系统。该系统基于SEM图像对纳米纤维膜进行三维拓扑特征的重构，构建具备相同拓扑特征的3D纳米纤维膜，通过所构建的纳米纤维膜进行过滤效率以及过滤压降的计算，并对空气过滤的过程进行可视化显示。系统具有一定的交互性功能，能够追踪颗粒物的运行轨迹等。对所搭建的纳米纤维膜空气过滤分析系统进行功能验证，针对纤维直径以及膜厚两个纤维结构特征分别进行过滤效率和过滤压降的实际测试和数值计算，对比结果表明，过滤效率的计算误差在10%以内，过滤压降的计算误差在20%以内，由此验证了所搭建的纳米纤维膜空气过滤分析系统具有良好的实用性。 参考文献 李学彬 ， 宫纯文 ， 徐青山 ， 等 . 气溶胶细粒子与能见度的相关性 ［J］. 光学 精密工程 ， 2008 ， 16 （ 7 ）： 1177 - 1180 . doi: 10.3321/j.issn:1004-924X.2008.07.005 http://dx.doi.org/10.3321/j.issn:1004-924X.2008.07.005 气溶胶细粒子与能见度的相关性 LI X B ， GONG CH W ， XU Q SH ， et al . Relationship between aerosol fine particle and visibility ［J］. Opt. Precision Eng. ， 2008 ， 16 （ 7 ）： 1177 - 1180 . （in Chinese） . doi: 10.3321/j.issn:1004-924X.2008.07.005 http://dx.doi.org/10.3321/j.issn:1004-924X.2008.07.005 Relationship between aerosol fine particle and visibility 王凌云 ， 马思远 ， 吴德志 . 电纺压电聚偏二氟乙烯有序纳米纤维及其在压力传感器中的应用 ［J］. 光学 精密工程 ， 2016 ， 24 （ 10 ）： 2498 - 2504 . doi: 10.3788/OPE.20162410.2498 http://dx.doi.org/10.3788/OPE.20162410.2498 电纺压电聚偏二氟乙烯有序纳米纤维及其在压力传感器中的应用 WANG L Y ， MA S Y ， WU D ZH . Electrospinning of aligned PVDF nanofibers with piezoelectricity and its application in pressure sensors ［J］. Opt. Precision Eng. ， 2016 ， 24 （ 10 ）： 2498 - 2504 . （in Chinese） . doi: 10.3788/OPE.20162410.2498 http://dx.doi.org/10.3788/OPE.20162410.2498 Electrospinning of aligned PVDF nanofibers with piezoelectricity and its application in pressure sensors 赵扬 ， 姜佳昕 ， 张恺 ， 等 . 基于电纺直写的图案化微纳结构喷印技术 ［J］. 光学 精密工程 ， 2016 ， 24 （ 9 ）： 2224 - 2231 . doi: 10.3788/OPE.20162409.2224 http://dx.doi.org/10.3788/OPE.20162409.2224 基于电纺直写的图案化微纳结构喷印技术 ZHAO Y ， JIANG J X ， ZHANG K ， et al . Precision deposition of micro/nano pattern printed by electrohydrodynamic direct-write ［J］. Opt. Precision Eng. ， 2016 ， 24 （ 9 ）： 2224 - 2231 . （in Chinese） . doi: 10.3788/OPE.20162409.2224 http://dx.doi.org/10.3788/OPE.20162409.2224 Precision deposition of micro/nano pattern printed by electrohydrodynamic direct-write CHEUNG C S ， CAO Y H ， YAN Z D . Numerical Model for particle deposition and loading in electret filter with rectangular split-type fibers ［J］. Computational Mechanics ， 2005 ， 35 （ 6 ）： 449 - 458 . doi: 10.1007/s00466-004-0634-5 http://dx.doi.org/10.1007/s00466-004-0634-5 Numerical Model for particle deposition and loading in electret filter with rectangular split-type fibers LI A ， AHMADI G . Dispersion and deposition of spherical particles from point sources in a turbulent channel flow ［J］. Aerosol Science and Technology ， 1992 ， 16 （ 4 ）： 209 - 226 . doi: 10.1080/02786829208959550 http://dx.doi.org/10.1080/02786829208959550 Dispersion and deposition of spherical particles from point sources in a turbulent channel flow CLEAVER J W ， YATES B . A sub layer model for the deposition of particles from a turbulent flow ［J］. Chemical Engineering Science ， 1975 ， 30 （ 8 ）： 983 - 992 . doi: 10.1016/0009-2509(75)80065-0 http://dx.doi.org/10.1016/0009-2509(75)80065-0 A sub layer model for the deposition of particles from a turbulent flow WANG H M ， ZHAO H B ， GUO Z L ， et al . Numerical simulation of particle capture process of fibrous filters using Lattice Boltzmann two-phase flow model ［J］. Powder Technology ， 2012 ， 227 ： 111 - 122 . doi: 10.1016/j.powtec.2011.12.057 http://dx.doi.org/10.1016/j.powtec.2011.12.057 Numerical simulation of particle capture process of fibrous filters using Lattice Boltzmann two-phase flow model WU J J ， AKAMPUMUZA O ， LIU P H ， et al . 3D structure design and simulation for efficient particles capture： the influence of nanofiber diameter and distribution ［J］. Materials Today Communications ， 2020 ， 23 ： 100897 . doi: 10.1016/j.mtcomm.2020.100897 http://dx.doi.org/10.1016/j.mtcomm.2020.100897 3D structure design and simulation for efficient particles capture： the influence of nanofiber diameter and distribution DONG M ， LI J Y ， SHANG Y ， et al . Numerical investigation on deposition process of submicron particles in collision with a single cylindrical fiber ［J］. Journal of Aerosol Science ， 2019 ， 129 ： 1 - 15 . doi: 10.1016/j.jaerosci.2018.12.001 http://dx.doi.org/10.1016/j.jaerosci.2018.12.001 Numerical investigation on deposition process of submicron particles in collision with a single cylindrical fiber PAN Z Y ， LIANG Y ， TANG M ， et al . Simulation of performance of fibrous filter media composed of cellulose and synthetic fibers ［J］. Cellulose ， 2019 ， 26 （ 12 ）： 7051 - 7065 . doi: 10.1007/s10570-019-02605-8 http://dx.doi.org/10.1007/s10570-019-02605-8 Simulation of performance of fibrous filter media composed of cellulose and synthetic fibers ZHAO X L ， LI Y Y ， HUA T ， et al . Low-resistance dual-purpose air filter releasing negative ions and effectively capturing PM 2.5 ［J］. ACS Applied Materials & Interfaces ， 2017 ， 9 （ 13 ）： 12054 - 12063 . doi: 10.1021/acsami.7b00351 http://dx.doi.org/10.1021/acsami.7b00351 Low-resistance dual-purpose air filter releasing negative ions and effectively capturing PM 2.5 SAMBAER W ， ZATLOUKAL M ， KIMMER D . 3D air filtration modeling for nanofiber based filters in the ultrafine particle size range ［J］. Chemical Engineering Science ， 2012 ， 82 ： 299 - 311 . doi: 10.1016/j.ces.2012.07.031 http://dx.doi.org/10.1016/j.ces.2012.07.031 3D air filtration modeling for nanofiber based filters in the ultrafine particle size range ZHAO X ， WANG S ， YIN X ， et al . Slip-effect functional air filter for efficient purification of PM 2.5 ［J］. Scientific Reports ， 2016 ， 6 ： 35472 . doi: 10.1038/srep35472 http://dx.doi.org/10.1038/srep35472 Slip-effect functional air filter for efficient purification of PM 2.5 WANG L ， GAO Y F ， XIONG J P ， et al . Biodegradable and high-performance multiscale structured nanofiber membrane as mask filter media via poly（lactic acid） electrospinning ［J］. Journal of Colloid and Interface Science ， 2022 ， 606 ： 961 - 970 . doi: 10.1016/j.jcis.2021.08.079 http://dx.doi.org/10.1016/j.jcis.2021.08.079 Biodegradable and high-performance multiscale structured nanofiber membrane as mask filter media via poly（lactic acid） electrospinning LIU H ， ZHANG S C ， LIU L F ， et al . High-performance PM 0.3 air filters using self-polarized electret nanofiber/nets ［J］. Advanced Functional Materials ， 2020 ， 30 （ 13 ）： 1909554 . doi: 10.1002/adfm.201909554 http://dx.doi.org/10.1002/adfm.201909554 High-performance PM 0.3 air filters using self-polarized electret nanofiber/nets PAIS V ， MOTA C ， BESSA J ， et al . Study of the filtration performance of multilayer and multiscale fibrous structures ［J］. Materials ， 2021 ， 14 （ 23 ）： 7147 . doi: 10.3390/ma14237147 http://dx.doi.org/10.3390/ma14237147 Study of the filtration performance of multilayer and multiscale fibrous structures SAMBAER W ， ZATLOUKAL M ， KIMMER D . 3D modeling of filtration process via polyurethane nanofiber based nonwoven filters prepared by electrospinning process ［J］. Chemical Engineering Science ， 2011 ， 66 （ 4 ）： 613 - 623 . doi: 10.1016/j.ces.2010.10.035 http://dx.doi.org/10.1016/j.ces.2010.10.035 3D modeling of filtration process via polyurethane nanofiber based nonwoven filters prepared by electrospinning process ZHANG T Y ， SUEN C Y . A fast parallel algorithm for thinning digital patterns ［J］. Communications of the ACM ， 1984 ， 27 （ 3 ）： 236 - 239 . doi: 10.1145/357994.358023 http://dx.doi.org/10.1145/357994.358023 A fast parallel algorithm for thinning digital patterns GHASEMI-MOBARAKEH L ， SEMNANI D ， MORSHED M . A novel method for porosity measurement of various surface layers of nanofibers mat using image analysis for tissue engineering applications ［J］. Journal of Applied Polymer Science ， 2007 ， 106 （ 4 ）： 2536 - 2542 . doi: 10.1002/app.26949 http://dx.doi.org/10.1002/app.26949 A novel method for porosity measurement of various surface layers of nanofibers mat using image analysis for tissue engineering applications MEAGHER D . Geometric modeling using octree encoding ［J］. Computer Graphics and Image Processing ， 1982 ， 19 （ 2 ）： 129 - 147 . doi: 10.1016/0146-664x(82)90104-6 http://dx.doi.org/10.1016/0146-664x(82)90104-6 Geometric modeling using octree encoding SHAO Z G ， JIANG J X ， WANG X ， et al . Self-powered electrospun composite nanofiber membrane for highly efficient air filtration ［J］. Nanomaterials （Basel， Switzerland） ， 2020 ， 10 （ 9 ）： 1706 . doi: 10.3390/nano10091706 http://dx.doi.org/10.3390/nano10091706 Self-powered electrospun composite nanofiber membrane for highly efficient air filtration 10.37188/OPE.20223009.1071 TQ340.64 TP391.9 A 1004-924X（2022）09-1071-09 郑高峰,姜佳昕,邓世卿等.纳米纤维膜三维重构与过滤行为分析系统[J].光学精密工程,2022,30(09):1071-1079. 郑高峰,姜佳昕,邓世卿等.纳米纤维膜三维重构与过滤行为分析系统[J].光学精密工程,2022,30(09):1071-1079. ZHENG Gaofeng,JIANG Jiaxin,DENG Shiqing,et al.Three-dimensional reconstruction and filtration behavior analysis system for nanofiber membrane[J].Optics and Precision Engineering,2022,30(09):1071-1079. ZHENG Gaofeng,JIANG Jiaxin,DENG Shiqing,et al.Three-dimensional reconstruction and filtration behavior analysis system for nanofiber membrane[J].Optics and Precision Engineering,2022,30(09):1071-1079. 航空科学基金资助项目（No.20180868001）；福建省科技计划项目高校产学合作项目（No. 2020H6003）；福建省增材制造创新中心开放课题基金资助项目（No. ZCZZ202-31）；广东省自然科学基金资助项目（No. 2022A1515010949） 2021-12-08 2022-01-14 2022-05-10 2022-05-26 光学精密工程 9 30