工程科学学报.第42卷.第2期：154-162.2020年2月 Chinese Journal of Engineering,Vol.42,No.2:154-162,February 2020 https://doi.org/10.13374/j.issn2095-9389.2019.02.24.004;http://cje.ustb.edu.cn 磁场形式及参数对单纤维捕集钢铁行业粉尘中 PM25性能影响 张俪安，刁永发四，庄加玮，周发山，沈恒根 东华大学环境科学与工程学院，上海201620 通信作者，E-mail:diaoyongfa@dhu.edu.cn 摘要目前钢铁行业已成为大气污染防治的重点，为解决现有钢铁行业对于PM2:细颗粒难以捕集的雄题，实现粉尘的超 低排放.基于CFD-DPM(computational fluid dynamics-.discrete phase model)方法对磁性纤维产生的磁场以及高梯度磁场等不同 磁场形式下单纤维对钢铁行业捕集PM2性能的影响进行研究，通过X射线衍射图谱分析可知钢铁行业生产过程产生的粉 尘因含有F©3O4以及单质Fe而具有磁特性，进而提出了利用磁场来增强单纤维捕集PM2s性能的方法.计算结果表明，在运 动轨迹方面，磁性纤维产生的磁场会在纤维周围形成引力区，高梯度磁场会在纤维周围形成2个引力区和2个斥力区：在捕 集性能方面，当粉尘粒径d,为0.5~1.0um,入口风速v≤0.2ms时，高梯度磁场下磁性纤维的捕集能力要强于单一磁性纤维 的捕集能力，若磁场强度H=0.5T,磁感应强度B=0.01T,1=0.1ms,高梯度磁场可以使单纤维的捕集效率提高为传统单纤维 捕集的28.32倍，若B=0.01T,1=0.1ms.磁性纤维产生的磁场可以使捕集效率提高为传统单纤维捕集的4.037倍：在磁性纤 维产生的磁场中，当磁感应强度B≥0.03T时，磁性单纤维对PM2、的捕集效率随着入口风速的增加而减小，后趋于稳定，当 B0.4ms时，捕集效率为0.B越大.捕集效率下降越快：捕集效率随着粉尘粒径增大 呈现先增加后减小的趋势. 关键词磁性纤维；高梯度磁场：PM25;铁磁性；捕集性能 分类号TF701.3 Performance of single fiber collection PM25 under different magnetic field forms in the iron and steel industry ZHANG Li-an,DIAO Yong-fa,ZHUANG Jia-wei.ZHOU Fa-shan,SHEN Heng-gen School of Environmental Science and Engineering College,Dong Hua University,Shanghai 201620,China Corresponding author,E-mail:diaoyongfa@dhu.edu.cn ABSTRACT At present,the steel industry has become the focus of air pollution prevention and control.To solve the difficulty in collecting PM2s fine particles and achieving ultra-low emission of dust,based on the method of computational fluid dynamics-discrete phase model (CFD-DPM),the influence of different magnetic field forms,such as magnetic field generated by magnetic fiber and high- gradient magnetic field,on the performance of PM2s collection in the iron and steel industry was studied.Through X-ray diffraction (XRD)analysis,it was found out that the dust produced in the iron and steel industry production process has magnetic characteristics due to the presence of FeO and elemental Fe,furthermore,the method of using magnetic field to enhance the PM2s collection performance of single fiber was proposed.The results show that the magnetic field generated by the magnetic fiber will form a gravitational region 收稿日期：2019-02-24 基金项目：国家重点研发计划资助项目(2018YFC0705300):中央高校基本科研业务费重点资助项目(2232017A-09)
磁场形式及参数对单纤维捕集钢铁行业粉尘中 PM2.5 性能影响 张俪安，刁永发苣，庄加玮，周发山，沈恒根 东华大学环境科学与工程学院，上海 201620 苣通信作者，E-mail: diaoyongfa@dhu.edu.cn 摘 要 目前钢铁行业已成为大气污染防治的重点，为解决现有钢铁行业对于 PM2.5 细颗粒难以捕集的难题，实现粉尘的超 低排放. 基于 CFD-DPM（computational fluid dynamics-discrete phase model）方法对磁性纤维产生的磁场以及高梯度磁场等不同 磁场形式下单纤维对钢铁行业捕集 PM2.5 性能的影响进行研究，通过 X 射线衍射图谱分析可知钢铁行业生产过程产生的粉 尘因含有 Fe3O4 以及单质 Fe 而具有磁特性，进而提出了利用磁场来增强单纤维捕集 PM2.5 性能的方法. 计算结果表明，在运 动轨迹方面，磁性纤维产生的磁场会在纤维周围形成引力区，高梯度磁场会在纤维周围形成 2 个引力区和 2 个斥力区；在捕 集性能方面，当粉尘粒径 dp 为 0.5～1.0 μm，入口风速 v≤0.2 m·s−1 时，高梯度磁场下磁性纤维的捕集能力要强于单一磁性纤维 的捕集能力，若磁场强度 H=0.5 T，磁感应强度 B=0.01 T，v=0.1 m·s−1，高梯度磁场可以使单纤维的捕集效率提高为传统单纤维 捕集的 28.32 倍，若 B=0.01 T，v=0.1 m·s−1，磁性纤维产生的磁场可以使捕集效率提高为传统单纤维捕集的 4.037 倍；在磁性纤 维产生的磁场中，当磁感应强度 B≥0.03 T 时，磁性单纤维对 PM2.5 的捕集效率随着入口风速的增加而减小，后趋于稳定，当 B0.4 m·s−1 时，捕集效率为 0，B 越大，捕集效率下降越快；捕集效率随着粉尘粒径增大 呈现先增加后减小的趋势. 关键词 磁性纤维；高梯度磁场；PM2.5；铁磁性；捕集性能 分类号 TF701.3 Performance of single fiber collection PM2.5 under different magnetic field forms in the iron and steel industry ZHANG Li-an，DIAO Yong-fa苣 ，ZHUANG Jia-wei，ZHOU Fa-shan，SHEN Heng-gen School of Environmental Science and Engineering College, Dong Hua University, Shanghai 201620, China 苣 Corresponding author, E-mail: diaoyongfa@dhu.edu.cn ABSTRACT At present, the steel industry has become the focus of air pollution prevention and control. To solve the difficulty in collecting PM2.5 fine particles and achieving ultra-low emission of dust, based on the method of computational fluid dynamics-discrete phase model (CFD-DPM), the influence of different magnetic field forms, such as magnetic field generated by magnetic fiber and highgradient magnetic field, on the performance of PM2.5 collection in the iron and steel industry was studied. Through X-ray diffraction (XRD) analysis, it was found out that the dust produced in the iron and steel industry production process has magnetic characteristics due to the presence of Fe3O4 and elemental Fe, furthermore, the method of using magnetic field to enhance the PM2.5 collection performance of single fiber was proposed. The results show that the magnetic field generated by the magnetic fiber will form a gravitational region 收稿日期: 2019−02−24 基金项目: 国家重点研发计划资助项目（2018YFC0705300）；中央高校基本科研业务费重点资助项目（2232017A-09） 工程科学学报，第 42 卷，第 2 期：154−162，2020 年 2 月 Chinese Journal of Engineering, Vol. 42, No. 2: 154−162, February 2020 https://doi.org/10.13374/j.issn2095-9389.2019.02.24.004; http://cje.ustb.edu.cn
张俪安等：磁场形式及参数对单纤维捕集钢铁行业粉尘中PM25性能影响 ·155· around the fiber,and the high-gradient magnetic field will form two gravitational regions and two repulsive regions around the fiber.In terms of the collection ability,when particle diameter d between 0.5 and 1.0 um,inlet velocity v2 m-s,the collection ability of magnetic fiber under the high-gradient magnetic field is stronger than that of the single magnetic fiber.If magnetic field intensity /=0.5 T,magnetic induction intensity B=0.01 T,and v=0.1 m-s,the high-gradient magnetic field can improve the single fiber collection efficiency by 28.32 times as much as the original;if B=0.01 T,v=0.1 ms,the magnetic field generated by the magnetic fiber can improve the single fiber collection efficiency by 4.037 times as much as the original.In terms of the collection law,in the magnetic field generated by the magnetic fiber,when the magnetic flux density B20.03 T,the collection efficiency of magnetic single fiber on PM2s decreases with the increase of inlet velocity speed and then tends to be stable.When B<0.03 T,the collection efficiency decreases with the inlet velocity speed.The collection efficiency increases with the increase of dust particle size.For the high-gradient magnetic field, the single fiber collection efficiency of PM2s particles also decreases with the increase of inlet velocity speed.When v.4 ms,the collection efficiency is 0.The larger B is,the faster the collection efficiency decreases.The collection efficiency increases first and then decreases with a increase in dust particle size. KEY WORDS magnetic fiber;high-gradient magnetic field;PM2s;ferromagnetism;collection performance 当前我国正处于工业发展的快速阶段，在发 方向会形成一些碎片链，这些链状团簇不断拉长 展的同时存在着产能过剩、重复建设严重、产业 对齐，最终接近稳定状态.Zhao等uo研究表明当 环保压力加大等深层次的矛盾及问题.尤其是钢 颗粒受磁场影响在流体中运动时，大颗粒易受磁 铁等有色金属冶炼行业，属于典型的重污染行业， 场影响，小颗粒易受流场影响.与此同时，高梯度 在烧结、炼铁、炼钢、轧钢等各个过程中都会产生 磁分离技术同样是一种回收弱磁性微粒的有效方 大量的粉尘四.现今，解决钢铁行业粉尘排放主要 法，而高梯度磁场的形成则是在均匀的背景磁场 以传统的袋式除尘器为主，其对粗颗粒收集效率 中填充饱和聚磁介质以产生高磁场梯度磁场 高达99.9%以上，但在PM25范围内存在穿透窗 Bak等]对影响高梯度磁场性能的参数进行了分 口，纤维难以捕集)，而PM25作为空气质量中最重 析研究.结果表明，在流体中作用于颗粒的磁场力 要的指标之一，高PM25指数的环境对人类健康的 与磁通密度和磁场梯度成正比.Zheng等研究 危害极为严重) 了高梯度磁场中特殊形状横截面矩阵的磁特性以 国内外有关纤维捕集细颗粒性能的研究，顾 及颗粒的被捕集情况.结果表明，高梯度磁场中椭 丛汇等经实验研究对比发现，同一时刻“Y”型截 圆形横截面矩阵表现出比圆形横截面矩阵更好的 面单纤维比圆形截面单纤维对颗粒截留率要高 磁特性，且捕集颗粒能力提高.通过上述研究发 Bao等研究表明，较高的纤维线密度，较低的弹 现，外磁场的形式以及参数对颗粒的行为有着重 性模量都可以增加纤维的捕集能力.Yang等通 要的影响 过考虑第一次碰撞的能量损失、线性阻尼和滚动 然而，利用磁场来增强纤维捕集钢铁行业生 摩擦阻力，较好的预测了单纤维捕集效率随斯托 产过程产生的铁磁性细粉尘的研究，其相关研究 克斯数的变化规律.Hosseini与Tafreshil通过数 文献鲜见报道.且上述研究中忽略了颗粒在磁场 值模拟确定单纤维捕集效率和颗粒沉积位点以及 力作用下被捕集过程的运动规律.不同的磁场形 三维树枝状颗粒沉积物的形成过程.通过上述研 式或者同种磁场形式在不同入口风速、粉尘粒 究发现，无论是纤维结构或者参数的变化，方法的 径、磁场强度等工况下所带来的增强效果有待研 改善只能提高单纤维对粉尘颗粒的捕集能力或增 究.为了在选择最优的磁场形式使得在同等条件 加对纤维沉积捕集的认知程度，但是不能从根本 下对PM25的捕集效果最好，因此本文基于CFD- 上解决PM25等微细颗粒物捕集效率低的难题.而 DPM方法.通过加入两种磁场形式的UDF添加磁 在传统捕集颗粒基础上利用磁场提高颗粒捕集的 场力，对其进行数值模拟.通过PM2.5被单纤维捕 研究，Huang等研究表明，在高磁场强度下低速 集过程的运动轨迹以及被捕集规律两个方向与传 运动时，纤维结构越复杂，沉降现象越明显.而在 统形式下单纤维捕集PM2s过程进行对比，探究出 低磁场强度下低速运动时，沉降现象与之相反 两种磁场形式下捕集性能的动力学行为和特征， Ke等对垂直和横向的均匀磁场中的流体-粒子 为用磁性纤维滤料控制钢铁等有色金属冶炼行业 相互作用进行了数值模拟.结果表明，颗粒沿流动 排放提供理论基础
around the fiber, and the high-gradient magnetic field will form two gravitational regions and two repulsive regions around the fiber. In terms of the collection ability, when particle diameter dp between 0.5 and 1.0 μm, inlet velocity v≤0.2 m·s−1, the collection ability of magnetic fiber under the high-gradient magnetic field is stronger than that of the single magnetic fiber. If magnetic field intensity H=0.5 T, magnetic induction intensity B=0.01 T, and v=0.1 m·s−1, the high-gradient magnetic field can improve the single fiber collection efficiency by 28.32 times as much as the original; if B=0.01 T, v=0.1 ms−1, the magnetic field generated by the magnetic fiber can improve the single fiber collection efficiency by 4.037 times as much as the original. In terms of the collection law, in the magnetic field generated by the magnetic fiber, when the magnetic flux density B≥0.03 T, the collection efficiency of magnetic single fiber on PM2.5 decreases with the increase of inlet velocity speed and then tends to be stable. When B0.4 ms−1, the collection efficiency is 0. The larger B is, the faster the collection efficiency decreases. The collection efficiency increases first and then decreases with a increase in dust particle size. KEY WORDS magnetic fiber；high-gradient magnetic field；PM2.5；ferromagnetism；collection performance 当前我国正处于工业发展的快速阶段，在发 展的同时存在着产能过剩、重复建设严重、产业 环保压力加大等深层次的矛盾及问题. 尤其是钢 铁等有色金属冶炼行业，属于典型的重污染行业， 在烧结、炼铁、炼钢、轧钢等各个过程中都会产生 大量的粉尘[1] . 现今，解决钢铁行业粉尘排放主要 以传统的袋式除尘器为主，其对粗颗粒收集效率 高达 99.9% 以上，但在 PM2.5 范围内存在穿透窗 口，纤维难以捕集[2] ，而 PM2.5 作为空气质量中最重 要的指标之一，高 PM2.5 指数的环境对人类健康的 危害极为严重[3] . 国内外有关纤维捕集细颗粒性能的研究，顾 丛汇等[4] 经实验研究对比发现，同一时刻“Y”型截 面单纤维比圆形截面单纤维对颗粒截留率要高. Bao 等[5] 研究表明，较高的纤维线密度，较低的弹 性模量都可以增加纤维的捕集能力. Yang 等[6] 通 过考虑第一次碰撞的能量损失、线性阻尼和滚动 摩擦阻力，较好的预测了单纤维捕集效率随斯托 克斯数的变化规律. Hosseini 与 Tafreshi[7] 通过数 值模拟确定单纤维捕集效率和颗粒沉积位点以及 三维树枝状颗粒沉积物的形成过程. 通过上述研 究发现，无论是纤维结构或者参数的变化，方法的 改善只能提高单纤维对粉尘颗粒的捕集能力或增 加对纤维沉积捕集的认知程度，但是不能从根本 上解决 PM2.5 等微细颗粒物捕集效率低的难题. 而 在传统捕集颗粒基础上利用磁场提高颗粒捕集的 研究，Huang 等[8] 研究表明，在高磁场强度下低速 运动时，纤维结构越复杂，沉降现象越明显. 而在 低磁场强度下低速运动时，沉降现象与之相反. Ke 等[9] 对垂直和横向的均匀磁场中的流体–粒子 相互作用进行了数值模拟，结果表明，颗粒沿流动 方向会形成一些碎片链，这些链状团簇不断拉长， 对齐，最终接近稳定状态. Zhao 等[10] 研究表明当 颗粒受磁场影响在流体中运动时，大颗粒易受磁 场影响，小颗粒易受流场影响. 与此同时，高梯度 磁分离技术同样是一种回收弱磁性微粒的有效方 法，而高梯度磁场的形成则是在均匀的背景磁场 中填充饱和聚磁介质以产生高磁场梯度磁场. Baik 等[11] 对影响高梯度磁场性能的参数进行了分 析研究. 结果表明，在流体中作用于颗粒的磁场力 与磁通密度和磁场梯度成正比. Zheng 等[12] 研究 了高梯度磁场中特殊形状横截面矩阵的磁特性以 及颗粒的被捕集情况. 结果表明，高梯度磁场中椭 圆形横截面矩阵表现出比圆形横截面矩阵更好的 磁特性，且捕集颗粒能力提高. 通过上述研究发 现，外磁场的形式以及参数对颗粒的行为有着重 要的影响. 然而，利用磁场来增强纤维捕集钢铁行业生 产过程产生的铁磁性细粉尘的研究，其相关研究 文献鲜见报道. 且上述研究中忽略了颗粒在磁场 力作用下被捕集过程的运动规律. 不同的磁场形 式或者同种磁场形式在不同入口风速、粉尘粒 径、磁场强度等工况下所带来的增强效果有待研 究. 为了在选择最优的磁场形式使得在同等条件 下对 PM2.5 的捕集效果最好，因此本文基于 CFDDPM 方法，通过加入两种磁场形式的 UDF 添加磁 场力，对其进行数值模拟. 通过 PM2.5 被单纤维捕 集过程的运动轨迹以及被捕集规律两个方向与传 统形式下单纤维捕集 PM2.5 过程进行对比，探究出 两种磁场形式下捕集性能的动力学行为和特征， 为用磁性纤维滤料控制钢铁等有色金属冶炼行业 排放提供理论基础. 张俪安等： 磁场形式及参数对单纤维捕集钢铁行业粉尘中 PM2.5 性能影响 · 155 ·
·156 工程科学学报，第42卷，第2期 1 数值计算模型 样品第1次测试 样品第2次测试 11气固两相流模型 样品第3次测试 平均值 对于内部气-固两相流动进行数值模拟需先 4 计算气相场.可采用标准k-ε标准模型、稳态、及 不可压缩模型进行数值模拟.控制方程1如下： 7.(pey)=0 (1) V.(spvv)=-sVp+V.(sT)+psg-F (2) 0 0.01 0.1 10 1001000 式中，p为流体的密度，kgm3;ε为体积分数项；v是 d/um 流体的速度，ms;p为计算单元的压力，N;t为流 图1钢铁厂排放粉尘的动态光散射技术测试结果 体黏附性应力张量；g为重力加速度，ms子；F为网 Fig.1 DLS test results of dust emission from steel works 格单元内受到的综合作用力，N 7000 颗粒的运动平衡方程表达式如下： 6000 8(Pp-P) CaC03(29.12) dVp=Fp(v-vp)+ +FM+FB+Fot (3) 5000 Pp 照4000 式中，。为颗粒运动速度，ms:FD为流体的曳 00 Fe04(35.42) Fez03(39.32) 力，N;pp为颗粒的密度，kgm3;FB为颗粒所受到 2000 Fe(44.14) 的布朗力，N;FM为通过UDF编程在高梯度磁场 1000 中以及磁性纤维产生中所受的磁场力，N.F指的 10 20 30 405060 708090 是在气固两相流中受到压力梯度力、虚拟质量 28M°) 力、Basset力、Magnus力、Saffman升力等一系列 因2钢铁厂排放粉尘的X射线衍射图谱 的作用力但是在单纤维捕集颗粒模型中，由于 Fig.2 XRD pattern of dust emission from steel works 各力在适用条件下对颗粒的影响很小，为了简化 图3(d)为P84纤维，从图中可以清晰的看出纤维 计算，均可忽略不计 呈现出圆柱形结构，为了达到更好的磁场叠加效 1.2粉尘特性 果，因此选用P84纤维为依据建立高梯度磁场下 钢铁厂排放粉尘的动态光散射技术(dynamic 单纤维捕集颗粒的模型.同时，磁性纤维在水基环 light scattering,.DLS)测试结果如图1所示，粉尘中 境制备负载时，铁酸钴纳米粒子等铁氧体很容易 含有直径d,小于等于2.5um的颗粒，因此钢铁等 将水中的羟基带到样品中，表面的活性基团 有色金属行业利用现有的除尘技术无法满足日趋 一OH上的氧原子和P84纤维上的C=O构成可p 严格的环保要求.图2给出了钢铁厂排放粉尘的 π共轭，导致O=C一O的体系稳定，在共沉反应时 X射线衍射图谱，CaCO3含量高是由于转炉炼钢 更容易在纤维表面负载 中用作造渣剂的Ca0在高温下会和CO2生成 1.4边界条件 CaC03,当衍射角20-35.42°和20=44.14时分别出 图4(a)为P84纤维的电镜扫描图，根据电镜 现了FeO4和Fe的特征峰，说明钢铁厂排放的粉 扫描结果，并经过简化创建了图4(b)的单纤维结 尘含有FeO4和Fe等成分，具有一定的铁磁性，容 构模型.具体边界条件设置如下：计算域入口边界 易被磁化.因此可以根据钢铁行业粉尘含有铁磁性 设为速度进口，出口边界设为压力出口：在模型计 的特性，利用铁磁性颗粒在磁场中磁化后会被磁 算中，粉尘颗粒简化为球形颗粒，颗粒的磁化率 性纤维吸引的优势来增强纤维捕集PM25的性能. Kp=0.025,a为填充率，是纤维横截面面积与假想控 1.3纤维的扫描形貌 制面的比值，根据Davies结合的实验结果，a=d,? 图3分别为PPS、FMS、PSA、P84等4种纤维 在0.6%~30%的范围内都是正确的7，当入口高 滤料的电镜扫描图像.由图可知，4种滤料中纤维 度h=5.3d时，则填充率a=0.035,因此模型的长、 分布以及结构是不同的.高梯度磁场的形成则是 宽、高分别定为240、120和80um是可行的 在均匀的背景磁场中填充饱和聚磁介质以产生 磁性纤维可直接通过纺丝或基体纤维的物 高磁场梯度磁场，而聚磁介质以圆柱形为主2. 理、化学改性制备，且纤维直径d=l5m.计算区
1 数值计算模型 1.1 气固两相流模型 对于内部气−固两相流动进行数值模拟需先 计算气相场. 可采用标准 k-ε 标准模型、稳态、及 不可压缩模型进行数值模拟. 控制方程[13] 如下： ∇ ·(ρεν) = 0 （1） ∇ ·(ερvv) = −ε∇p+∇ ·(ετ)+ρεg− F （2） ρ ε τ 式中， 为流体的密度，kg·m−3 ； 为体积分数项；v 是 流体的速度，m·s−1 ；p 为计算单元的压力，N； 为流 体黏附性应力张量；g 为重力加速度，m·s−2 ；F 为网 格单元内受到的综合作用力，N. 颗粒的运动平衡方程表达式[14] 如下： dvp dt = FD ( v−vp ) + g ( ρp −ρ ) ρp + FM + FB + Fot （3） 式中， vp 为颗粒运动速度，m·s−1 ；FD 为流体的曳 力，N；ρp 为颗粒的密度，kg·m−3 ；FB 为颗粒所受到 的布朗力，N；FM 为通过 UDF 编程在高梯度磁场 中以及磁性纤维产生中所受的磁场力，N. Fot 指的 是在气固两相流中受到压力梯度力、虚拟质量 力、Basset 力、Magnus 力、Saffman 升力等一系列 的作用力[15] . 但是在单纤维捕集颗粒模型中，由于 各力在适用条件下对颗粒的影响很小，为了简化 计算，均可忽略不计. 1.2 粉尘特性 钢铁厂排放粉尘的动态光散射技术（dynamic light scattering， DLS）测试结果如图 1 所示，粉尘中 含有直径 dp 小于等于 2.5 μm 的颗粒，因此钢铁等 有色金属行业利用现有的除尘技术无法满足日趋 严格的环保要求. 图 2 给出了钢铁厂排放粉尘的 X 射线衍射图谱，CaCO3 含量高是由于转炉炼钢 中用作造渣剂 的 CaO 在高温下会 和 CO2 生 成 CaCO3，当衍射角 2θ=35.42°和 2θ=44.14°时分别出 现了 Fe3O4 和 Fe 的特征峰，说明钢铁厂排放的粉 尘含有 Fe3O4 和 Fe 等成分，具有一定的铁磁性，容 易被磁化. 因此可以根据钢铁行业粉尘含有铁磁性 的特性，利用铁磁性颗粒在磁场中磁化后会被磁 性纤维吸引的优势来增强纤维捕集 PM2.5 的性能. 1.3 纤维的扫描形貌 图 3 分别为 PPS、FMS、PSA、P84 等 4 种纤维 滤料的电镜扫描图像. 由图可知，4 种滤料中纤维 分布以及结构是不同的，高梯度磁场的形成则是 在均匀的背景磁场中填充饱和聚磁介质以产生 高磁场梯度磁场，而聚磁介质以圆柱形为主[12， 16] ， 图 3（d）为 P84 纤维，从图中可以清晰的看出纤维 呈现出圆柱形结构，为了达到更好的磁场叠加效 果，因此选用 P84 纤维为依据建立高梯度磁场下 单纤维捕集颗粒的模型. 同时，磁性纤维在水基环 境制备负载时，铁酸钴纳米粒子等铁氧体很容易 将水中的羟基带到样品中 ，表面的活性基团 ―OH 上的氧原子和 P84 纤维上的 C＝O 构成可 p- π 共轭，导致 O＝C―O 的体系稳定，在共沉反应时 更容易在纤维表面负载. 1.4 边界条件 χp 图 4（a）为 P84 纤维的电镜扫描图，根据电镜 扫描结果，并经过简化创建了图 4（b）的单纤维结 构模型. 具体边界条件设置如下：计算域入口边界 设为速度进口，出口边界设为压力出口；在模型计 算中，粉尘颗粒简化为球形颗粒，颗粒的磁化率 =0.025，α 为填充率，是纤维横截面面积与假想控 制面的比值，根据 Davies 结合的实验结果，α=dp 2 /h 2 在 0.6%～30% 的范围内都是正确的[17] ，当入口高 度 h=5.3df 时，则填充率 α=0.035，因此模型的长、 宽、高分别定为 240、120 和 80 μm 是可行的. 磁性纤维可直接通过纺丝或基体纤维的物 理、化学改性制备，且纤维直径 df=15 μm. 计算区 0.01 0.1 1 10 100 1000 0 1 2 3 4 5 6 体积加权/% dp /µm 样品第1次测试 样品第2次测试 样品第3次测试 平均值 图 1 钢铁厂排放粉尘的动态光散射技术测试结果 Fig.1 DLS test results of dust emission from steel works 0 10 20 30 40 50 60 70 80 90 0 1000 2000 3000 4000 5000 6000 7000 Fe2O3 (39.32°) Fe (44.14°) CaCO3 (29.12°) 相对强度 2θ/(°) Fe3O4 (35.42°) 图 2 钢铁厂排放粉尘的 X 射线衍射图谱 Fig.2 XRD pattern of dust emission from steel works · 156 · 工程科学学报，第 42 卷，第 2 期
张俪安等：磁场形式及参数对单纤维捕集钢铁行业粉尘中PM25性能影响 ·157 (a (b) 100μm 100um (c) (d) 100μm 100um 图3不同纤维滤料的电镜扫描形貌.(a)PPS;(b)FMS:(c)PSA:(d)P84 Fig.3 SEM images of different fiber filter materials:(a)PPS;(b)FMS;(c)PSA;(d)P84 (b) 磁性纤维 磁场强度H 压力出口 速度进口 15 um 80μm 120μm 240μm 对称边界 图4计算区域及边界条件设置图.(a)P84纤维电镜扫描：(b)单纤维模型计算边界条件 Fig.4 Calculation area and boundary condition setting diagram:(a)P84 fiber SEM image;(b)single fiber model used to calculate the boundary conditions 域四周的边界则根据单纤维结构特点设为对称边 的，且与压力经验公式(4)对比误差都在5%的范 界条件：纤维表面边界设为无滑移边界条件，且数 围内.选取其中一种工况进行效率计算，当速度 值计算区域的网格采用的是六面体结构化网格 =0.2ms,d,=2.5um时，单纤维捕集效率与Davies 传统单纤维捕集、磁性单纤维捕集以及高梯度磁 效率修正经验公式(6)的误差分别为8.51%、 场下磁性单纤维的捕集示意图如图5所示，图中 3.50%、1.50%,后两者的误差都在5%范围内，根 H为磁场强度，B为磁感应强度 据网格数量和误差综合考虑选取55万的网格用 1.5模型验证 于数值模型的计算 为了去除网格数量对数值模拟计算准确性的 (4) 影响，对模型进行网格独立性验证，计算不同网格 Ap=f(a)und 密度下的压力损失，数值模拟计算如图6所示，当 f(a)=64a2(1+56a） (5) 网格数为14万、55万、90万左右时，单纤维模型 结构的压力损失随入口风速的变化规律是一致 7=70(1+4.5a) (6)
域四周的边界则根据单纤维结构特点设为对称边 界条件；纤维表面边界设为无滑移边界条件，且数 值计算区域的网格采用的是六面体结构化网格. 传统单纤维捕集、磁性单纤维捕集以及高梯度磁 场下磁性单纤维的捕集示意图如图 5 所示，图中 H 为磁场强度，B 为磁感应强度. 1.5 模型验证 为了去除网格数量对数值模拟计算准确性的 影响，对模型进行网格独立性验证，计算不同网格 密度下的压力损失，数值模拟计算如图 6 所示，当 网格数为 14 万、55 万、90 万左右时，单纤维模型 结构的压力损失随入口风速的变化规律是一致 的，且与压力经验公式（4）对比误差都在 5% 的范 围内. 选取其中一种工况进行效率计算，当速度 v=0.2 m·s−1 ，dp=2.5 μm 时，单纤维捕集效率与 Davies 效 率 修 正 经 验 公 式 （ 6） 的 误 差 分 别 为 8.51%、 3.50%、1.50%，后两者的误差都在 5% 范围内，根 据网格数量和误差综合考虑选取 55 万的网格用 于数值模型的计算. ∆p = f (α) µvd d 2 f （4） f (α) = 64α 3 2 ( 1+56α 3 ) （5） η = η0 (1+4.5α) （6） 100 µm (b) 100 µm (a) (c) 100 µm (d) 100 µm 图 3 不同纤维滤料的电镜扫描形貌. （a） PPS；（b） FMS；（c） PSA；（d） P84 Fig.3 SEM images of different fiber filter materials: (a) PPS; (b) FMS; (c) PSA; (d) P84 对称边界 磁场强度H 磁性纤维 压力出口 速度进口 15 µm df (a) (b) 120 µm 240 µm 80 µm h 100 µm 100 µm 6.0 mm Scios WD 6.0 mm HFW 207 µm det ETD mag 1000× HV 5.00 KV 12/5/2018 2:27:01 PM 图 4 计算区域及边界条件设置图. （a） P84 纤维电镜扫描；（b） 单纤维模型计算边界条件 Fig.4 Calculation area and boundary condition setting diagram: (a) P84 fiber SEM image; (b) single fiber model used to calculate the boundary conditions 张俪安等： 磁场形式及参数对单纤维捕集钢铁行业粉尘中 PM2.5 性能影响 · 157 ·
·158 工程科学学报，第42卷，第2期 (a) 与流场垂直区域 (b) 与流场水平区域 普通纤维 (c) (d) 磁性纤维 磁性纤维 图5不同形式磁场示意图.(a)粉尘颗粒被纤维捕集区域示意图：(b)无磁场：(c)磁性纤维产生的磁场：(d)高梯度磁场 Fig.5 Schematics of different magnetic fields:(a)area where PMs dust particles are collected by fibers.(b)no magnetic field;(c)magnetic field generated by magnetic fiber,(d)high-gradient magnetic fields 2.5 是可以结合磁选领域中的宏观模型来间接验证模 ■一Davies经验公式 ●一940754网格数 拟的正确性.数值模拟结果与河南理工大学的王 2.0F ▲一551300网格数 发辉与铁占续和东南大学的杨荣清例在磁选 T-148272网格数 1.5 领域中模拟金属磁介质捕集颗粒的运动轨迹进行 ed/d 1.0 了对比，运动轨迹相同；此外与熊大和2研究高梯 度磁选中垂直磁场和水平磁场对捕集影响的实验 0.5F 结果进行对比，通过颗粒的堆积分布发现数值模 拟的颗粒物分布和实际情况一致，因此说明采用 00 0.10.20.30.40.50.60.70.80.9 v/m-s) 的数学模型和模拟方法是正确的 14r 图6网格独立性检验 R 一Davies修正公式 Fig.6 Mesh independence test ●一Davies经验公式 10 ▲一500颗粒数 0=0.16R,+(0.25+0.4Rp)S1-0.0263R,·S2(7) -200颗粒数 式中：△p为压力降，Pa;o为计算捕集效率，%； Rp为直接碰撞系数；S为斯托克斯数；n为计算修 正后的捕集效率，%：d为单纤维厚度，m. 为了验证单纤维结构捕集粉尘颗粒数值模拟 的准确性，图7为不同S下，数值模拟值与经验公 0.0 0.4 0.8 1.21.62.0 St 式对比图.结果表明，数值模拟计算结果与 困7数值模拟值与经验公式对比 Davies修正经验公式(6)吻合度较高，适合颗粒物 Fig.7 Comparison between numerical simulation values and empirical 捕集效率的计算.过滤效率的误差主要来源于 formula Davies效率经验公式考虑了颗粒的拦截和惯性效 2试验结果 应，而当粒径在0.5~1.0um时，忽略了颗粒碰撞后 受范德华力的团聚情况，因此与经验公式有一定 2.1不同磁场形式下单纤维对PM25运动轨迹的 差距，而对于而对于Fluent用户自定义函数UDF 影响 正确性的验证，虽然在纤维微观领域很少研究，但 图8为=0.1ms,d。=1.0μm时，在图5所示
η0 = 0.16[ Rp + ( 0.25+0.4Rp ) St] −0.0263Rp · St2 （7） η0 η 式中 ： Δp 为压力降 ， Pa； 为计算捕集效率 ， %； RP 为直接碰撞系数；St 为斯托克斯数； 为计算修 正后的捕集效率，%；d 为单纤维厚度，μm. 为了验证单纤维结构捕集粉尘颗粒数值模拟 的准确性，图 7 为不同 St 下，数值模拟值与经验公 式 对 比 图 . 结 果 表 明 ， 数 值 模 拟 计 算 结 果 与 Davies 修正经验公式（6）吻合度较高，适合颗粒物 捕集效率的计算. 过滤效率的误差主要来源于 Davies 效率经验公式考虑了颗粒的拦截和惯性效 应，而当粒径在 0.5～1.0 μm 时，忽略了颗粒碰撞后 受范德华力的团聚情况，因此与经验公式有一定 差距，而对于而对于 Fluent 用户自定义函数 UDF 正确性的验证，虽然在纤维微观领域很少研究，但 是可以结合磁选领域中的宏观模型来间接验证模 拟的正确性. 数值模拟结果与河南理工大学的王 发辉与铁占续[18] 和东南大学的杨荣清[19] 在磁选 领域中模拟金属磁介质捕集颗粒的运动轨迹进行 了对比，运动轨迹相同；此外与熊大和[20] 研究高梯 度磁选中垂直磁场和水平磁场对捕集影响的实验 结果进行对比，通过颗粒的堆积分布发现数值模 拟的颗粒物分布和实际情况一致，因此说明采用 的数学模型和模拟方法是正确的. 2 试验结果 2.1 不同磁场形式下单纤维对 PM2.5 运动轨迹的 影响 图 8 为 v=0.1 m·s−1 ，dp=1.0 μm 时，在图 5 所示 (a) (b) (c) (d) 与流场垂直区域 与流场水平区域 v X Z v X Z 普通纤维 X Z B 磁性纤维 X Z H B 磁性纤维 图 5 不同形式磁场示意图. （a） 粉尘颗粒被纤维捕集区域示意图；（b） 无磁场；（c） 磁性纤维产生的磁场；（d） 高梯度磁场 Fig.5 Schematics of different magnetic fields: (a) area where PM2.5 dust particles are collected by fibers; (b) no magnetic field; (c) magnetic field generated by magnetic fiber; (d) high-gradient magnetic fields 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 0.5 1.0 1.5 2.0 2.5 ∆p/Pa v/(m⋅s −1) Davies 经验公式 940754 网格数 551300 网格数 148272 网格数 图 6 网格独立性检验 Fig.6 Mesh independence test 0.0 0.4 0.8 1.2 1.6 2.0 0 2 4 6 8 10 12 14 η/% St Davies 修正公式 Davies 经验公式 500 颗粒数 200 颗粒数 图 7 数值模拟值与经验公式对比 Fig.7 Comparison between numerical simulation values and empirical formula · 158 · 工程科学学报，第 42 卷，第 2 期
张俪安等：磁场形式及参数对单纤维捕集钢铁行业粉尘中PM25性能影响 ·159· 的3种工况下的数值模拟计算结果，PM25在3种 合图10可知，v≤0.2ms时，高梯度磁场的捕集能 不同工况下运动轨迹以及被捕集区域截然不同 力要强于磁性纤维的捕集能力，当>0.3ms单 在两种磁场形式下，单纤维对PM25的捕集强度要 一磁性纤维的捕集能力要强于高梯度磁场中磁性 远大于传统的单纤维捕集，这是因为传统的单纤 纤维的捕集能力，由此说明高梯度磁场更适用于 维捕集机理仅依靠布朗、惯性、拦截三种机制进 低过滤风速的环境 行捕集，捕集效率低，但是当PM25经磁场磁化后 20 会受磁场力的作用向磁性纤维表面运动并被纤维 18 16 捕集，捕集效率高，为了对比两种磁场效果采用 14 B和H+B的对比方式，当场强分别为0.05T和0.5T+ 0.05T时，结合示意图5(a)可知，在高梯度磁场 10 只 .80.0s7 B-00LI 0- (H+B)中，在磁性纤维周围存在两个引力区和两个 6 斥力区叫，引力区位于与流场垂直区域，斥力区位 2 于与流场水平区域，其捕集区域集中在引力区，这 与磁选理论中的结论是一致的；而在磁性纤维产 0.5 0.6 0.70.80.9 1.0 dWμm 生的磁场中，引力区环绕纤维周围 图9不同磁场形式下PM5捕集效果对比 c Fig.9 Comparison of PM2s trapping effects under different magnetic field forms 2.3入口风速对PM25捕集性能的影响 图10(a)为不同入口风速下磁性单纤维捕集 PM25效率关系图，如图所示，当B≥0.03T时，磁性 单纤维对PM2s的捕集效率随着入口风速的增加 而减小，后趋于稳定.当B1.0μm时，单纤维捕集以 影响很小；随着磁性单纤维的磁感应强度逐渐减 惯性和拦截作用为主)，捕集效率都较高.因此对 小，当磁场力对颗粒的影响很小时，若流动占主导 于0.5~1.0m粒径段颗粒的捕集研究显得尤为重 地位，被捕集的颗粒数会持续下降 要.当=0.1ms,d,=0.5~1.0m时，两种磁场形 图10(b)为PM25在高梯度磁场分离过程中， 式下，单纤维对粉尘颗粒捕集效率对比如图9所 不同入口风速下磁性纤维对PM25效率关系图，如 示，高梯度磁场中磁性纤维的捕集能力要强于单 图所示，当H=0.5T,B=0.01~0.05T时，在高梯度 一磁性纤维的捕集能力.取其中一组数据，当H= 磁场中磁性纤维对PM25的捕集效率随入口风速 0.5T,B=0.01T时，磁性单纤维捕集效率的加权平 逐渐减小，磁性纤维磁感应强度越大，捕集效率下 均值为5.943%，而数值模拟计算得到传统未加磁 降越快.当v>0.4ms时.捕集效率为0.根据图11(b) 场单纤维捕集效率的加权平均值为0.2099%，高梯 运动轨迹图可知，由于在高梯度磁场中捕集区域 度磁场可以使单纤维的捕集效率提高为传统单纤 位于与流场垂直的磁性纤维表面，风速对于颗粒 维捕集的28.32倍；当B=0.01T时，磁性单纤维捕 的捕集影响很大，同时，气流速度越大，相应的颗 集效率的加权平均值为0.8473%，磁性纤维产生的 粒在磁场中的作用时间越短，颗粒所受的磁场作 磁场可以使单纤维的捕集效率提高为传统单纤维 用的效果就会在一定程度上减弱，使得落在引力 捕集的4.037倍.随着粉尘粒径的增加，单纤维对 区的粉尘颗粒逐渐减小，因此捕集效率减少直至 粉尘颗粒的捕集效率都随粒径呈递增趋势.当结 为0.且随着入口风速的增加，斥力区的“空腔”缩
的 3 种工况下的数值模拟计算结果，PM2.5 在 3 种 不同工况下运动轨迹以及被捕集区域截然不同. 在两种磁场形式下，单纤维对 PM2.5 的捕集强度要 远大于传统的单纤维捕集，这是因为传统的单纤 维捕集机理仅依靠布朗、惯性、拦截三种机制进 行捕集，捕集效率低，但是当 PM2.5 经磁场磁化后 会受磁场力的作用向磁性纤维表面运动并被纤维 捕集，捕集效率高，为了对比两种磁场效果采用 B 和 H+B 的对比方式，当场强分别为 0.05 T 和 0.5 T+ 0.05 T 时，结合示意图 5（ a）可知，在高梯度磁场 （H+B）中，在磁性纤维周围存在两个引力区和两个 斥力区[21] ，引力区位于与流场垂直区域，斥力区位 于与流场水平区域，其捕集区域集中在引力区，这 与磁选理论中的结论是一致的；而在磁性纤维产 生的磁场中，引力区环绕纤维周围. 2.2 不同磁场形式下单纤维对 PM2.5 捕集效率对比 从机理分析，当 dp1.0 μm 时，单纤维捕集以 惯性和拦截作用为主[23] ，捕集效率都较高. 因此对 于 0.5～1.0 μm 粒径段颗粒的捕集研究显得尤为重 要. 当 v=0.1 m·s−1 ，dp=0.5～1.0 μm 时，两种磁场形 式下，单纤维对粉尘颗粒捕集效率对比如图 9 所 示，高梯度磁场中磁性纤维的捕集能力要强于单 一磁性纤维的捕集能力. 取其中一组数据，当 H= 0.5 T，B=0.01 T 时，磁性单纤维捕集效率的加权平 均值为 5.943%，而数值模拟计算得到传统未加磁 场单纤维捕集效率的加权平均值为 0.2099%，高梯 度磁场可以使单纤维的捕集效率提高为传统单纤 维捕集的 28.32 倍；当 B=0.01 T 时，磁性单纤维捕 集效率的加权平均值为 0.8473%，磁性纤维产生的 磁场可以使单纤维的捕集效率提高为传统单纤维 捕集的 4.037 倍. 随着粉尘粒径的增加，单纤维对 粉尘颗粒的捕集效率都随粒径呈递增趋势. 当结 合图 10 可知，v≤0.2 m·s−1 时，高梯度磁场的捕集能 力要强于磁性纤维的捕集能力，当 v>0.3 m·s−1 单 一磁性纤维的捕集能力要强于高梯度磁场中磁性 纤维的捕集能力，由此说明高梯度磁场更适用于 低过滤风速的环境. 2.3 入口风速对 PM2.5 捕集性能的影响 图 10（a）为不同入口风速下磁性单纤维捕集 PM2.5 效率关系图，如图所示，当 B≥0.03 T 时，磁性 单纤维对 PM2.5 的捕集效率随着入口风速的增加 而减小，后趋于稳定. 当 B0.4 m·s−1 时，捕集效率为0. 根据图11（b） 运动轨迹图可知，由于在高梯度磁场中捕集区域 位于与流场垂直的磁性纤维表面，风速对于颗粒 的捕集影响很大，同时，气流速度越大，相应的颗 粒在磁场中的作用时间越短，颗粒所受的磁场作 用的效果就会在一定程度上减弱，使得落在引力 区的粉尘颗粒逐渐减小，因此捕集效率减少直至 为 0. 且随着入口风速的增加，斥力区的“空腔”缩 v v v (a) (b) (c) 图 8 3 种工况下 PM2.5 的运动轨迹（v=0.1 m·s−1 ， dp=1.0 μm）. （a） 无 磁场；（b） 磁性纤维产生的磁场；（c） 高梯度磁场 Fig.8 PM2.5 movement trajectory under three working conditions: (a) no magnetic field; (b) magnetic field generated by magnetic fiber; (c) high-gradient magnetic fields 0.5 0.6 0.7 0.8 0.9 1.0 0 2 4 6 8 10 12 14 16 18 20 η/% dp /µm H=0.5 T, B=0.05 T H=0.5 T, B=0.03 T H=0.5 T, B=0.01 T B=0.05 T B=0.01 T 图 9 不同磁场形式下 PM2.5 捕集效果对比 Fig.9 Comparison of PM2.5 trapping effects under different magnetic field forms 张俪安等： 磁场形式及参数对单纤维捕集钢铁行业粉尘中 PM2.5 性能影响 · 159 ·
.160 工程科学学报，第42卷，第2期 10 (a) 30 (b) 8 B=0.05T ◆-B=0.04T -■-=0.5T,B=0.05T B=0.03T ●-H=0.5T.B=0.04T 6 -B=0.02T ▲-=0.5T,B=0.03T B=0.011 10 H=0.5T,B=0.02T -=0.5T.B=0.01T 2 0 0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.5 v/(m-s-) v/(m-s-) 图10不同人口风速与捕集效率关系图(d,=1.0m.(a)磁性纤维产生的磁场：(b)高梯度磁场 Fig.10 Relation between different inlet velocity speeds and collection efficiencies (d=1.0 um):(a)magnetic field generated by magnetic fiber;(b)high- gradient magnetic field (a) 0.1ms 0.2m-s 0.3ms 0.4m-s 0.5ms 一一 (b) 0.1m-s1 0.2m-s1 0.3ms 0.4m-s 0.5ms 图11不同入口风速下粉尘在磁场中的运动轨迹(d。=1.0um).(a)磁性纤维产生的磁场-0.05T:(b)高梯度磁场=0.5T,B=0.05T Fig.11 Movement trajectory of dust in magnetic field at different inlet velocity speeds(d,=1.0 um):(a)magnetic field generated by magnetic fiber (B=0.05 T);(b)high-gradient magnetic field (H=0.5 T,B=0.05 T) 小，这是因为当颗粒接近磁性纤维时，虽然在斥力 因为，当磁感应强度较小时，此时磁场虽然很弱 区所受的斥力相同，但是颗粒速度越大，运动状态 但磁场力的出现会显著增加纤维对颗粒的捕集， 越不易改变，在斥力区运动的时间越长，导致斥力 因此增长率不会随粒径增加而减小.随着磁感应 区的“空腔”缩小 强度的继续增加，磁场增强，但磁场力的变化相比 2.4粉尘粒径对PM25捕集性能的影响 之前减弱，且纤维捕集面的面积一定，因此捕集效 图12(a)为磁性纤维对不同粉尘粒径的捕集 率的增长率会减小 效率关系图，如图所示，当B=0.01~0.05T时，磁性 图12(b)为高梯度磁场分离过程中，磁性纤维 单纤维对PM25的捕集效率随粒径的增加而增大 对不同粉尘粒径捕集效率关系图，由图可知，当 这是因为磁性纤维周围都是引力区，捕集面大，增 H=0.5T,B=0.01~0.05T时，随着粉尘粒径的增加， 加了粉尘颗粒受磁场作用被磁性单纤维捕集的概 磁性纤维对PM25的捕集效率呈现先增大后减小 率，且粒径越大，粉尘颗粒所受的磁场力越大，因 的趋势.根据图13(b)运动轨迹图可知.当粉尘粒 此捕集效率越高.根据图13(a)运动轨迹图可知， 径较小时，气流携带颗粒能力强，被气流携带的颗 粉尘粒径越大，落在纤维上的颗粒越多，且分布面 粒在经过引力区时很少被磁性纤维捕获，因此捕 积越大.同时，捕集效率随粒径的增长率由于磁性 集效率较小.随着粒径的增大，在引力区由于受到 纤维磁感应强度增加会有一定程度的减小.这是 磁场力增大使得捕集效率提高，但是随着粒径继
小，这是因为当颗粒接近磁性纤维时，虽然在斥力 区所受的斥力相同，但是颗粒速度越大，运动状态 越不易改变，在斥力区运动的时间越长，导致斥力 区的“空腔”缩小. 2.4 粉尘粒径对 PM2.5 捕集性能的影响 图 12（a）为磁性纤维对不同粉尘粒径的捕集 效率关系图，如图所示，当 B=0.01～0.05 T 时，磁性 单纤维对 PM2.5 的捕集效率随粒径的增加而增大. 这是因为磁性纤维周围都是引力区，捕集面大，增 加了粉尘颗粒受磁场作用被磁性单纤维捕集的概 率，且粒径越大，粉尘颗粒所受的磁场力越大，因 此捕集效率越高. 根据图 13（a）运动轨迹图可知， 粉尘粒径越大，落在纤维上的颗粒越多，且分布面 积越大. 同时，捕集效率随粒径的增长率由于磁性 纤维磁感应强度增加会有一定程度的减小. 这是 因为，当磁感应强度较小时，此时磁场虽然很弱， 但磁场力的出现会显著增加纤维对颗粒的捕集， 因此增长率不会随粒径增加而减小. 随着磁感应 强度的继续增加，磁场增强，但磁场力的变化相比 之前减弱，且纤维捕集面的面积一定，因此捕集效 率的增长率会减小. 图 12（b）为高梯度磁场分离过程中，磁性纤维 对不同粉尘粒径捕集效率关系图，由图可知，当 H=0.5 T，B=0.01～0.05 T 时，随着粉尘粒径的增加， 磁性纤维对 PM2.5 的捕集效率呈现先增大后减小 的趋势. 根据图 13（b）运动轨迹图可知，当粉尘粒 径较小时，气流携带颗粒能力强，被气流携带的颗 粒在经过引力区时很少被磁性纤维捕获，因此捕 集效率较小. 随着粒径的增大，在引力区由于受到 磁场力增大使得捕集效率提高，但是随着粒径继 H=0.5 T, B=0.05 T H=0.5 T, B=0.04 T H=0.5 T, B=0.03 T H=0.5 T, B=0.02 T H=0.5 T, B=0.01 T 0.1 0.2 0.3 0.4 0.5 0 2 4 6 8 10 η/% η/% (a) B=0.05 T B=0.04 T B=0.03 T B=0.02 T B=0.01 T 0.1 0.2 0.3 0.4 0.5 0 5 10 15 20 (b) v/(m⋅s −1) v/(m⋅s −1) 图 10 不同入口风速与捕集效率关系图（dp=1.0 μm）. （a） 磁性纤维产生的磁场；（b） 高梯度磁场 Fig.10 Relation between different inlet velocity speeds and collection efficiencies (dp=1.0 μm): (a) magnetic field generated by magnetic fiber; (b) highgradient magnetic field v v (a) (b) 0.1 m⋅s −1 0.2 m⋅s −1 0.3 m⋅s −1 0.4 m⋅s −1 0.5 m⋅s −1 0.1 m⋅s −1 0.2 m⋅s −1 0.3 m⋅s −1 0.4 m⋅s −1 0.5 m⋅s −1 图 11 不同入口风速下粉尘在磁场中的运动轨迹（dp=1.0 μm）. （a） 磁性纤维产生的磁场 B=0.05 T； （b）高梯度磁场 H=0.5 T，B=0.05 T Fig.11 Movement trajectory of dust in magnetic field at different inlet velocity speeds (dp=1.0 μm): (a) magnetic field generated by magnetic fiber (B=0.05 T); (b) high-gradient magnetic field (H=0.5 T, B=0.05 T) · 160 · 工程科学学报，第 42 卷，第 2 期
张俪安等：磁场形式及参数对单纤维捕集钢铁行业粉尘中PM2.5性能影响 ·161· 25(a 20F(b) ■-H=0.5TB=0.05T ●-H=0.5T.B-=0.04T 30 量-B=0.05T ▲-H=0.5TB=0.031 ●-B=0.04T 15 -H=0.5T.B=0.02T ▲一B=0.03T ◆-H=0.5TB=0.01T 15 -B=0.02T -B=0.01T 10 10 05 1.0 1.5 2.0 2.5 05 1.0 1.5 2.0 )5 d/um d/um 图12不同粒径与捕集效率关系图(1-0.1ms).(a)磁性纤维产生的磁场：(b)高梯度磁场 Fig.12 Relation between particle size and collection efficiencies (1=0.1 m's):(a)magnetic field generated by magnetic fiber;(b)high-gradient magnetic field (a) 0.5um 1.0μm 1.5um 2.0um 2.5μm 6 0.5um 1.0um 1.5m 2.0μm 2.5um 图13不同粒径粉尘在磁场中的运动轨迹(-0.1ms).(a)磁性纤维产生的磁场(B-0.05T):(b)高梯度磁场(H=0.5T,B=0.05T) Fig.13 Movement trajectory of dust with different particle sizes in magnetic field(v=0.1 m's):(a)magnetic field generated by magnetic fiber(B=0.05 T):(b)high-gradient magnetic field (H=0.5 T,B=0.05 T) 续增大，在接近纤维的过程中由于受斥力作用增 在纤维周围形成引力区，高梯度磁场会在纤维周 强使颗粒向远离纤维的方向运动，造成“空腔”增 围形成2个引力区和2个斥力区.当H=0.5T,B= 大.当颗粒再经过引力区时，由于磁场强度会随着 0.01T,=0.1ms时，高梯度磁场可以使单纤维的 颗粒距离纤维距离的增大而减弱.再加上在引力 捕集效率提高为原来的28.32倍；当B=0.01T时， 区颗粒受引力运动方向与流场几乎垂直，受流场 磁性纤维产生的磁场可以使单纤维的捕集效率提 影响很大，此时很难再被纤维捕集，捕集效率减小 高为原来的4.037倍 当d,=2.5m时，H=0.5T,B≤0.03T时，=0： (2)在0.5um≤d,≤1.0um的粒径范围内，当v≤ H=0.5T,B>0.03T时，>0.这是因为当d,≥1.5um时， 0.2ms时，PM25在高梯度磁场中(H=0.5T, 粉尘粒径的增加导致的“空腔”大小增加缓慢，即 B=0.01~0.05T,Xp=0.025)捕集效率比单一磁性纤 粉尘颗粒距离单纤维的距离变化较小，当外磁场 维(B=0.01~0.05T,Xp=0.025)的捕集效率大，> 一定，磁感应强度增加时，叠加后的磁场强度增 03ms时，单一磁性纤维的捕集能力比高梯度磁 强，此时磁力增加，而此时颗粒又会被捕集 场中磁性纤维的捕集能力强，由此说明高梯度磁 场更适用于低风速的环境 3结论 (3)对于磁性纤维产生的磁场(B=0.01~0.05T, (1)在运动轨迹方面，磁性纤维产生的磁场会 Kp=0.025),当B≥0.03T时，磁性单纤维对PM2s颗
续增大，在接近纤维的过程中由于受斥力作用增 强使颗粒向远离纤维的方向运动，造成“空腔”增 大. 当颗粒再经过引力区时，由于磁场强度会随着 颗粒距离纤维距离的增大而减弱. 再加上在引力 区颗粒受引力运动方向与流场几乎垂直，受流场 影响很大，此时很难再被纤维捕集，捕集效率减小. 当 dp=2.5 μm 时 ， H=0.5 T， B≤0.03 T 时 ， η=0； H=0.5 T，B>0.03 T 时，η>0. 这是因为当 dp≥1.5 μm 时， 粉尘粒径的增加导致的“空腔”大小增加缓慢，即 粉尘颗粒距离单纤维的距离变化较小，当外磁场 一定，磁感应强度增加时，叠加后的磁场强度增 强，此时磁力增加，而此时颗粒又会被捕集. 3 结论 （1）在运动轨迹方面，磁性纤维产生的磁场会 在纤维周围形成引力区，高梯度磁场会在纤维周 围形成 2 个引力区和 2 个斥力区. 当 H=0.5 T，B= 0.01 T，v=0.1 m·s−1 时，高梯度磁场可以使单纤维的 捕集效率提高为原来的 28.32 倍 ；当 B=0.01 T 时 ， 磁性纤维产生的磁场可以使单纤维的捕集效率提 高为原来的 4.037 倍. χp χp （2）在 0.5 μm≤dp≤1.0 μm 的粒径范围内，当 v≤ 0.2 m·s−1 时 ， PM2.5 在 高 梯 度 磁 场 中 （ H=0.5 T， B=0.01～0.05 T， =0.025）捕集效率比单一磁性纤 维 （ B=0.01～ 0.05 T， =0.025）的捕集效率大 ， v> 0.3 m·s−1 时，单一磁性纤维的捕集能力比高梯度磁 场中磁性纤维的捕集能力强，由此说明高梯度磁 场更适用于低风速的环境. χp （3）对于磁性纤维产生的磁场（B=0.01～0.05 T， =0.025），当 B≥0.03 T 时，磁性单纤维对 PM2.5 颗 0.5 1.0 1.5 2.0 2.5 0 5 10 15 20 25 dp /μm dp /μm 0.5 1.0 1.5 2.0 2.5 0 5 10 15 20 H=0.5 T, B=0.05 T H=0.5 T, B=0.04 T H=0.5 T, B=0.03 T H=0.5 T, B=0.02 T H=0.5 T, B=0.01 T η/% η/% (a) B=0.05 T B=0.04 T B=0.03 T B=0.02 T B=0.01 T (b) 图 12 不同粒径与捕集效率关系图（v=0.1 m·s−1）. （a） 磁性纤维产生的磁场；（b） 高梯度磁场 Fig.12 Relation between particle size and collection efficiencies (v=0.1 m·s−1): (a) magnetic field generated by magnetic fiber; (b) high-gradient magnetic field v v (a) (b) 0.5 µm 1.0 µm 1.5 µm 2.0 µm 2.5 µm 0.5 µm 1.0 µm 1.5 µm 2.0 µm 2.5 µm 图 13 不同粒径粉尘在磁场中的运动轨迹（v=0.1 m·s−1）.（a） 磁性纤维产生的磁场（B=0.05 T）； （b） 高梯度磁场（H=0.5 T，B=0.05 T） Fig.13 Movement trajectory of dust with different particle sizes in magnetic field (v=0.1 m·s−1): (a) magnetic field generated by magnetic fiber (B=0.05 T); (b) high-gradient magnetic field (H=0.5 T, B=0.05 T) 张俪安等： 磁场形式及参数对单纤维捕集钢铁行业粉尘中 PM2.5 性能影响 · 161 ·
.162 工程科学学报，第42卷，第2期 粒的捕集效率随着人口风速的增加而减小，后趋 particle within a high gradient magnetic separator.Physica C, 于稳定；当B0.4ms时，捕集效率 [14]Zhu H J.Fluid Analysis and Engineering Simulation of Fluent 12. 为0：单纤维对PM2.5的捕集效率随着粉尘粒径呈 Beijing:Tsinghua University Press,2011 现先增加后减小的规律. (朱红钧.Fluent12流体分析及工程仿其.北京：清华大学出版社， 2011) 参考文献 [15]Tripathy S K,Bhoja S K,Kumar C R,et al.A short review on hydraulic classification and its development in mineral industry. [1]Zhang Y.Eco-design of steel industry and policy options in China Powder Technol,2015.270:205 China's Popul Resour Environ,2012,22(7):162 [16]Eisentrager A,Vella D,Griffiths I M.Particle capture efficiency in (张雅.钢铁产业生态化设计与政策选择.中国人口·资源与环 a multi-wire model for high gradient magnetic separation.Appl 境，2012,22(7)：162) Ph3Lei,2014,105(3):033508 [2]Xiong G L,Li S Q,Chen S,et al.Development of advanced [17]Zhao HL,Fu H M,Lei C L,et al.Effect of fiber's cross-sectional electrostatic precipitation technologies for reducing PM2s shape on fiber collection efficiency and pressure drop.J Donghua emissions from coal-fired power plants.Proc CSEE,2015,35(9): Univ Nat Sci,2016,42(1):86 2217 (赵洪亮，付海明，雷陈磊，等，纤维藏面形状对纤维捕集效率及 (熊桂龙，李水清，陈晟，等.增强PM,脱除的新型电除尘技术的 压力损失的影响.东华大学学报：自然科学版，2016,42(1)：86) 发展.中国电机工程学报，2015,35(9)：2217) [18]Wang F H,Tie Z X.Numerical simulation for high gradient [3]Qu Y,Qian X,Song HQ,et al.Machine-learning-based model and magnetic field located single magnetic medium in entrapping simulation analysis of PM2.5 concentration prediction in Beijing magnetism particles.Coal Prepar Technol,2012(2):20 Chin J Eng,2019,41(3):401 (王发辉，铁占续.高梯度磁场中单根磁介质捕集磁性微粒的数 (曲悦，钱旭，宋洪庆，等.基于机器学习的北京市PM2.5浓度预 值模拟.选煤技术，2012(2)：20) 测模型及模拟分析.工程科学学报，2019,41(3)：401) [19]Yang R Q.Investigation on Kinetic Characteristic of Magnetic [4]Gu C H.La S W.Li R.et al.Influence of fiber on filtration Fine Particles in High Gradient Magnetic Field [Dissertation] performance for PM2 s.C/ESC J,2014,65(6):2137 Nanjing:Southeast University,2006 (顾从汇，吕士武，李瑞，等.纤维对PM,过滤性能的影响.化工 (杨荣清.高梯度磁场中磁性可吸人颗粒物动力学特性研究学 学报，2014,65(6)：2137) 位论文].南京：东南大学，2006) [5]Bao L,Musadiq M,Kijima T,et al.Influence of fibers on the dust [20]Xiong D H.Study on comparison between vertical and horizontal dislodgement efficiency of bag filters.Text Res J,2014,84(7) magnetic fields in pulsating high gradient magnetic separation. 764 Metal Mine,2004(10):24 [6]Yang MM,Li S Q,Yao Q.Mechanistic studies of initial (熊大和，脉动高梯度磁选垂直磁场与水平磁场对比研究.金属 deposition of fine adhesive particles on a fiber using discrete. 矿山，2004(10)：24) element methods.Powder Technol,2013,248:44 [21]Sun Z Y.Magnetic Separation Theory.Changsha:Central South [7] Hosseini S A.Tafreshi H V.Modeling particle-loaded single fiber University Press,2007 efficiency and fiber drag using ANSYS-Fluent CFD code.Compu (孙仲元.磁选理论.长沙：中南大学出版社，2007) Fluids,2012,66:157 [22]Qian F P.Wang H G.Numerical analysis on particle capture [8] Huang S,Zhang X M,Tafu M,et al.Study on subway particle characteristics of fibrous filters with random structure.J Civil capture by ferromagnetic mesh filter in nonuniform magnetic field. Architect Environ Eng,2010,32(6):120 Sep Purif Technol,2015,156:642 (钱付平，王海刚.随机排列纤维过滤器颗粒捕集特性的数值研 [9]Ke C H,Shu S,Zhang H,et al.LBM-IBM-DEM modelling of 究.土木建筑与环境工程，2010,32(6)：120) magnetic particles in a fluid.Powder Techmnol,2016,314:264 [23]Zhu H,Fu H M,Kang Y M.Numerical analysis of pressure drop [10]Zhao L,Li X L,Sun W Q,et al.Experimental study on bag and inertial collection efficiency of a single fiber.China Environ filtration enhanced by magnetic aggregation of fine particles from Sc,2017,37(4):1298 hot metal casting process.Powder Technol,2018,327:255 (朱辉，付海明，亢燕铭.单纤维过滤阻力与惯性捕集效率数值 [11]Baik S K.Ha D W,Kwon J M,et al.Magnetic force on a magnetic 分析.中国环境科学，2017,37(4)：1298)
粒的捕集效率随着入口风速的增加而减小，后趋 于稳定；当 B0.4 m·s−1 时，捕集效率 为 0；单纤维对 PM2.5 的捕集效率随着粉尘粒径呈 现先增加后减小的规律. 参 考 文 献 Zhang Y. Eco-design of steel industry and policy options in China. China's Popul Resour Environ, 2012, 22（7）: 162 （张雅. 钢铁产业生态化设计与政策选择. 中国人口•资源与环 境, 2012, 22（7）：162 ） [1] Xiong G L, Li S Q, Chen S, et al. Development of advanced electrostatic precipitation technologies for reducing PM2.5 emissions from coal-fired power plants. Proc CSEE, 2015, 35（9）: 2217 （熊桂龙, 李水清, 陈晟, 等. 增强PM2.5脱除的新型电除尘技术的 发展. 中国电机工程学报, 2015, 35（9）：2217 ） [2] Qu Y, Qian X, Song H Q, et al. Machine-learning-based model and simulation analysis of PM2.5 concentration prediction in Beijing. Chin J Eng, 2019, 41（3）: 401 （曲悦, 钱旭, 宋洪庆, 等. 基于机器学习的北京市PM2.5浓度预 测模型及模拟分析. 工程科学学报, 2019, 41（3）：401 ） [3] Gu C H, Lü S W, Li R, et al. Influence of fiber on filtration performance for PM2.5. CIESC J, 2014, 65（6）: 2137 （顾从汇, 吕士武, 李瑞, 等. 纤维对PM2.5过滤性能的影响. 化工 学报, 2014, 65（6）：2137 ） [4] Bao L, Musadiq M, Kijima T, et al. Influence of fibers on the dust dislodgement efficiency of bag filters. Text Res J, 2014, 84（7）: 764 [5] Yang M M, Li S Q, Yao Q. Mechanistic studies of initial deposition of fine adhesive particles on a fiber using discreteelement methods. Powder Technol, 2013, 248: 44 [6] Hosseini S A, Tafreshi H V. Modeling particle-loaded single fiber efficiency and fiber drag using ANSYS-Fluent CFD code. Comput Fluids, 2012, 66: 157 [7] Huang S, Zhang X M, Tafu M, et al. Study on subway particle capture by ferromagnetic mesh filter in nonuniform magnetic field. Sep Purif Technol, 2015, 156: 642 [8] Ke C H, Shu S, Zhang H, et al. LBM-IBM-DEM modelling of magnetic particles in a fluid. Powder Technol, 2016, 314: 264 [9] Zhao L, Li X L, Sun W Q, et al. Experimental study on bag filtration enhanced by magnetic aggregation of fine particles from hot metal casting process. Powder Technol, 2018, 327: 255 [10] [11] Baik S K, Ha D W, Kwon J M, et al. Magnetic force on a magnetic particle within a high gradient magnetic separator. Physica C, 2013, 484: 333 Zheng X Y, Wang Y H, Lu D F. Investigation of the particle capture of elliptic cross-sectional matrix for high gradient magnetic separation. Powder Technol, 2016, 297: 303 [12] Qian F P, Huang N J, Zhu X J, et al. Numerical study of the gassolid flow characteristic of fibrous media based on SEM using CFD-DEM. Powder Technol, 2013, 249: 63 [13] Zhu H J. Fluid Analysis and Engineering Simulation of Fluent 12. Beijing: Tsinghua University Press, 2011 （朱红钧. Fluent12流体分析及工程仿真. 北京: 清华大学出版社, 2011） [14] Tripathy S K, Bhoja S K, Kumar C R, et al. A short review on hydraulic classification and its development in mineral industry. Powder Technol, 2015, 270: 205 [15] Eisenträger A, Vella D, Griffiths I M. Particle capture efficiency in a multi-wire model for high gradient magnetic separation. Appl Phys Lett, 2014, 105（3）: 033508 [16] Zhao H L, Fu H M, Lei C L, et al. Effect of fiber’s cross-sectional shape on fiber collection efficiency and pressure drop. J Donghua Univ Nat Sci, 2016, 42（1）: 86 （赵洪亮, 付海明, 雷陈磊, 等. 纤维截面形状对纤维捕集效率及 压力损失的影响. 东华大学学报: 自然科学版, 2016, 42（1）：86 ） [17] Wang F H, Tie Z X. Numerical simulation for high gradient magnetic field located single magnetic medium in entrapping magnetism particles. Coal Prepar Technol, 2012（2）: 20 （王发辉, 铁占续. 高梯度磁场中单根磁介质捕集磁性微粒的数 值模拟. 选煤技术, 2012（2）：20 ） [18] Yang R Q. Investigation on Kinetic Characteristic of Magnetic Fine Particles in High Gradient Magnetic Field [Dissertation]. Nanjing: Southeast University, 2006 （杨荣清. 高梯度磁场中磁性可吸入颗粒物动力学特性研究[学 位论文]. 南京: 东南大学, 2006） [19] Xiong D H. Study on comparison between vertical and horizontal magnetic fields in pulsating high gradient magnetic separation. Metal Mine, 2004（10）: 24 （熊大和. 脉动高梯度磁选垂直磁场与水平磁场对比研究. 金属 矿山, 2004（10）：24 ） [20] Sun Z Y. Magnetic Separation Theory. Changsha: Central South University Press, 2007 （孙仲元. 磁选理论. 长沙: 中南大学出版社, 2007） [21] Qian F P, Wang H G. Numerical analysis on particle capture characteristics of fibrous filters with random structure. J Civil Architect Environ Eng, 2010, 32（6）: 120 （钱付平, 王海刚. 随机排列纤维过滤器颗粒捕集特性的数值研 究. 土木建筑与环境工程, 2010, 32（6）：120 ） [22] Zhu H, Fu H M, Kang Y M. Numerical analysis of pressure drop and inertial collection efficiency of a single fiber. China Environ Sci, 2017, 37（4）: 1298 （朱辉, 付海明, 亢燕铭. 单纤维过滤阻力与惯性捕集效率数值 分析. 中国环境科学, 2017, 37（4）：1298 ） [23] · 162 · 工程科学学报，第 42 卷，第 2 期