1、第 卷 第 期 年 月西安科技大学学报 .吴燕妮，姜华，张浩睿，等 矿用对旋风机气固两相流磨损特性 西安科技大学学报，（）：，（）：收稿日期：基金项目：国家重点基础研究发展计划（）第一作者：吴燕妮，女，陕西汉中人，硕士研究生，：通信作者：姜华，女，河北良香人，博士，副教授，：矿用对旋风机气固两相流磨损特性吴燕妮，姜 华，张浩睿，宫武旗（西安科技大学 能源学院，陕西 西安；西安交通大学 能源与动力工程学院，陕西 西安）摘 要：为探索对旋风机在矿井含尘环境中叶片的磨损情况，以及颗粒长久性磨损造成叶片表面粗糙度增加的后果，基于 湍流模型和 磨损模型，利用多相流求解技术对 矿用对旋式局部通风机进行气固
2、两相流数值模拟，分析不同颗粒浓度和粒径条件下叶片磨损分布特征，探究叶片表面粗糙度对风机内部流场和整机性能的影响。结果表明：气固两相流场内颗粒浓度对叶片磨损率的影响较大，磨损率随浓度的增加而增大；叶片表面粗糙度的均匀和非均匀分布形式导致前缘和尾缘不同程度的能量耗散，在进行对旋风机气固两相非定常流动特性的分析时，应采用粗糙度非均匀分布的方法真实反映叶片粗糙度变化的影响；风机性能随叶片粗糙度的增大而降低，大流量工况下的影响更为显著，但随工作时间的增长，叶片表面粗糙度均匀性提高，流场恶化情况减缓。因此，适时修复或更换磨损叶片，有利于风机安全高效运行。关键词：对旋风机；气固两相流；叶片磨损；叶片粗糙度；
3、风机性能中图分类号：文献标志码：文章编号：（）：开放科学（资源服务）标识码（）：，（，；，）：，：；引 言作为矿井通风系统的重要设备，对旋风机工作环境恶劣，长期输送含有煤粉的气体。当含尘气流作为工作介质通过风机时，被气流挟带的固体颗粒将对所流经的壁面产生磨损作用，颗粒的物理磨损将使叶片表面粗糙度增大。经过除尘设备净化后，巷道内悬浮的粉尘粒径集中在 ，粉尘浓度不超过 ，对旋风机输送含小粒径低浓度颗粒气流，颗粒相流动特征是研究叶片磨损的出发点。而对旋风机以其风量大、风压高和性能好等特点，已被广泛应用于工矿企业等气流含尘量较大领域，含尘气流中固体颗粒对风机性能及流动特性影响的研究具有实际意义和理论价
4、值。近年来，众多学者通过试验测定的方法对叶轮机械两相流进行研究，等用 硅砂进行压气机磨损试验，得出颗粒碰撞角与叶片磨损率呈正相关。李德顺等搭建了风洞试验台，探究了风沙冲击速度、冲击角度对风力机叶片磨损特性的影响。等基于射流试验研究了不同叶片材料的抗磨损特性，得出降低叶片表面粗糙度有利于提高叶轮机械性能。随着数值模拟技术的发展，国内外针对叶轮机械两相流的研究主要集中在离心泵、离心风机、离心压缩机、轴流风机和轴流压气机等方面。在压气机、燃气轮机等叶片粗糙度与气动性能的研究中，虽然能够体现粗糙度的不利影响，但已有文献指出，叶片前缘、吸力面较流场其他区域对粗糙度的变化更敏感，说明叶轮机械流场与粗糙度的
5、分布形式之间有一定关联。目前，有关对旋风机的研究主要集中在性能改善、振动噪声、旋转失速等方面，关于矿井通风机故障分析的文章中，虽有提及叶片磨损的成因，但几乎没有气固两相流磨损特征的探讨。关于叶轮机械颗粒磨损表征方法的研究中，叶片表面粗糙度的处理方法还未形成比较一致的意见。为此，文中以安设在掘进工作面进风侧的 矿用对旋式局部通风机为研究对象，分析气固两相流场内，不同颗粒浓度和颗粒粒径条件下叶片磨损分布特征，以及叶片粗糙度分布形式和粗糙度大小对风机性能的影响。相关结论在对旋风机防磨修复和安全运行方面具有工程意义，可为风机设计提供理论指导。计算模型及方法 计算模型图 为 矿用对旋式局部通风机模型，主
6、要参数见表。该风机由集流器、前级叶轮、后级叶轮和出口流道组成，前后两级叶轮的旋转方向相反。计算域包括集流器、前级叶轮、后级叶轮和出口流道 个部分，其中叶轮区域为运动域，其他区域为静止域。图 矿用对旋式局部通风机结构模型 网格划分与无关性验证图 为对旋风机模型的网格划分，考虑到叶顶间隙流动的复杂性，对叶片周围网格进行局部加密。为排除网格数量对计算结果的影响，分别以 种数量不同的网格进行无关性验证。由图 可知，网格数由 万增加到 万后，对旋风机的静 西 安科技大学学报 年第 卷第 期吴燕妮，等：矿用对旋风机气固两相流磨损特性压、功率和效率分别稳定在 ，和 附近。结合计算精度和时间成本，最终确定网格
7、总数为 万，进口区网格数为 万，前级叶轮区网格数为 万，后级叶轮区网格数为 万，出口区网格数为 万。表 矿用对旋式局部通风机主要参数 参数设计流量（）设计转速（）前级叶片数前级叶片叶顶间隙 前级叶片安装角（）后级叶片数后级叶片叶顶间隙 后级叶片安装角（）轮毂比数值 图 计算模型网格（局部）示意 （）图 网格无关性验证 气固两相流的计算方法采用气相稳态流场作为初场，计算两相流非稳态流场，对计算条件作出如下假设。）气固两相不可压缩，具有相同的温度场。）固体颗粒为粒径相同、密度均匀且无旋转效应的球体。）颗粒受力主要考虑曳力和重力，忽略热泳力、虚拟质量力以及压力梯度力等。）经过除尘设备净化后的气体含尘
8、浓度较低，故流动为稀相两相流，忽略颗粒对气相的作用力和颗粒之间的接触力。颗粒轨迹的求解采用颗粒轨道模型在拉格朗日坐标系下描述离散相颗粒的运动轨迹。根据作用在颗粒上各种力的平衡关系获得的单个颗粒运动方程如下。（）（）（）（）（）（）式中 为颗粒质量，；为颗粒速度，；为曳力，；为浮力，；为旋转作用力，；为无量纲曳力系数；为空气密度，；为颗粒有效面积，；为气流速度，；为颗粒密度，；为重力加速度，；为旋转域转速，；为颗粒所在位置，。磨损模型采用 等提出，并进行过试验验证的磨损关系式，以及 等改进的 磨损模型计算磨损量 （）（）其中（）（）|（）式中 为磨损率，（）；为碰撞材料之间的作用力之比；为壁面材
9、料的特征指数；为颗粒运动方向与壁面的夹角，（）。叶片表面粗糙度粗糙度是描述磨损表面形貌特征的参数，对于粗糙表面，仍可借助壁面函数描述近壁面的流动情况。处于粗糙表面对数律层的流速分布如下。（）（）其中 （）式中 为无量纲流速；为 常数；为无量纲距离；与表面粗糙度有关，一般取 ；为无量纲粗糙度偏移量；为无量纲砂粒粗糙度。将叶片的粗糙表面看作一层等粒径砂粒紧密排列形成的峰谷，以砂粒圆心所在平面为基准面，砂粒高出基准面的高度即为叶片表面粗糙度，用无量纲砂粒粗糙度表示。（）式中 为叶片表面粗糙度，；为壁面摩擦速度，；为空气动力黏性系数，（）。模拟设置连续相为 ，标准大气压下的空气，湍流模型采用 模型，近
10、壁面采用标准壁面函数。离散相为 固体颗粒，摩尔质量为 ，密度为 ，颗粒随空气从风机入口进入，位置随机分布。对旋风机计算模型的验证为验证计算模型的准确性，根据 工业通风机用标准化风道性能试验规定，采用 型管道进口和自由出口法搭建的 矿用对旋式局部通风机试验台如图 所示。在设计转速和设计流量下，测得风机静压、功率和效率分别为 ，和.。图 型对旋风机试验台 通过调节集流器多孔栅板改变风机工作状态，对设计转速下 个不同工况的风机性能进行测试，模拟值与测试值的偏差如图 所示。可以看出，模拟结果与试验数据相近，静压、功率和效率的最大偏差分别为 ，和 ，均低于，说明计算模型和网格数量能够真实反映该风机的气动
11、性能。模拟结果与分析 叶片磨损分布特征为探究对旋风机在含尘矿井中的磨损分布特征，对不同颗粒条件下的风机气固两相流场进行模拟。两级叶片在不同颗粒浓度下的磨损分布如图 所示。磨损率随着颗粒浓度的增加而增大，颗粒浓度为，时，叶片最大磨损率分别为 ，（）。前级叶片吸力面前缘的磨损范围比后级叶片小，压力面的磨损范围比后级叶片大。前级叶片压力面前缘和中部靠近尾缘的磨损较为严重，后级叶片压力面前缘和 叶高处的磨损较为严重。这是因为颗粒在叶轮旋转域内受离心力作用，运动方向沿径向向叶顶发生偏移，在前级叶轮区域内，颗粒还未加速运动至叶顶就已经与叶片发生碰撞，所以磨损分布较为均匀。而受前级叶轮的加速和级间非定常流动
12、的影响，颗粒在后级叶轮区域内的运动速度提高，与叶片的碰撞位置靠近叶顶，说明离心力对后级的作用更显著。图 对旋风机性能偏差 风机叶片在不同颗粒粒径下的磨损分布如图 所示。随着粒径由 增加至 ，叶片最大磨损率基本不变，分别为 ，（），但叶片的磨损形态发生了较大改变。当粒径大于 后，前级叶片压力面中部和后级叶片压力面 叶高及以上部位开始出现磨损，后级叶片吸力面前缘的磨损开始发散。比较粒径 和 时压力面的磨损形态，可以发现随着粒径的增大，前级叶片的磨损开始向尾缘移动，后级叶片的磨损开始向叶顶移动。比较图 和图，发现随着颗粒浓度的增大，磨损形态基本不变，最大磨损率变化较大。而随着颗粒粒径的增大，磨损形态
13、变化较大，最大磨损率变化很小。这说明颗粒浓度对叶片磨损率的影响较大，对碰撞位置的影响较小，颗粒粒径对碰撞位置的影响较大，对叶片磨损率的影响较小。这是因为粒径相同的颗粒，浓度越大，与叶片碰撞的 西 安科技大学学报 年第 卷第 期吴燕妮，等：矿用对旋风机气固两相流磨损特性次数越多，造成的叶片磨损率越大。粒径不同的颗粒，受到的作用力大小不同，在相同浓度下，颗粒运动轨迹和碰撞位置不同，造成了不同的叶片磨损形态。因此，为了减缓叶片磨损，应尽量确保风机气固两相流场内的颗粒浓度不超过 ，颗粒粒径不超过 。图 不同颗粒浓度下叶片磨损率分布 图 不同颗粒粒径下叶片磨损率分布 叶片表面粗糙度分布形式对流场的影响磨
14、损对风机的直接影响是改变了叶片的粗糙度。为反映磨损风机的实际情况，参考对旋风机磨损分布特征对叶片进行分区划分。前级叶片划分为 个区域，后级叶片划分为 个区域，如图 所示。为保证模拟计算结果的收敛性，对区域连接处的网格进行局部加密处理。两级叶片工作单位时间后的非均匀粗糙度及均匀粗糙度见表。非均匀粗糙度按照对旋风机在颗粒粒径 ，浓度 条件下各区域的平均磨损率确定。均匀粗糙度用各区域的非均匀粗糙度按面积占比相加得到的平均粗糙度表示。对旋风机的设计转速较高，会在叶片表面形图 叶片粗糙度的分区 成旋涡，涡量的大小可以反映旋涡的强度。了解边界涡量在叶片表面的分布，可分析对旋风机的流场特性。涡量在直角坐标系
15、下的表达式为|式中 为涡量，；，分别为，坐标轴方向上的速度，。图 为 种粗糙度分布形式下，前级叶片表面叶高处涡量沿流向的分布。图中对旋风机涡量最大的位置在叶片尾缘，其次是吸力面前缘附近。这是因为尾缘涡与粗糙叶片和机壳相互作用，使气流在叶片附面层内的旋涡不断增强。种粗糙度分布形 式 下，叶 片 前 缘 涡 量 最 大 相 差.，这是由于非均匀分布的叶片前缘粗糙度更大，导致的附面层增厚，涡量明显增加。而粗糙度均匀分布时，叶片附面层不存在梯度变化，涡量的变化相对稳定。表 不同分布形式的叶片粗糙度 前级叶片第 区第 区第 区第 区后级叶片第 区第 区第 区第 区面积 占比 非均匀粗糙度 均匀粗糙度 图
16、 叶高处叶片表面涡量分布 由图 前级叶片 叶高截面上湍流动能的分布情况可知，气流经过叶片时的湍流动能会大于其他部位，在文献中得到了证实。，这 处的流场恶化严重，这与前级叶轮处存在一定尺度的漩涡，以及叶片前缘和尾缘涡量的扩大有关。与均匀粗糙度分布相比，采用非均匀粗糙度分布后，叶片前缘和尾缘出现流动分离，负向涡量在吸力面尾缘和压力面前缘聚集，使，这 处的湍流动能增大，种粗糙度分布形式的湍流动能最大相差 。由此可见，粗糙度的不同分布形式导致叶片周围涡量强度不同，叶片表面粗糙度的突变使湍流动能加剧，气流在相同位置处的能量损耗增加。图 叶高截面上的湍流动能分布 叶片表面粗糙度分布形式对性能的影响图 为不同粗糙度分布形式下的对旋风机性能曲线。粗糙度分布形式对整机性能的影响并不显著，静压和效率的最大偏差分别为 和.。与叶片表面粗糙度均匀分布形式相比，非均匀分布形式在小流量工况下的静压和功率略高，但设计工况的性能参数基本一致。与前文研究结果比较，发现不同粗糙度分布形式下叶片附近流场的差异显著，极有可能影响非定常计算。因此，在以对旋风机整机性能为研究重点时，可以采用表面粗糙度均匀分布的方法，加快计算速度
