机械密封是一种用于流体机械或者动力机械中保证设备运行可靠的零部件。其工作原理是通过动环与静环之间相对高速的旋转产生一层薄液膜，这层极薄的液膜将密封介质与缓冲液，或者将缓冲液与大气隔离开，从而达到密封的效果。机械密封工作可靠、适用范围广，炼油装置中８５％以上的机泵 使用了机械密封［１］。
随着泵内介质变得复杂，对于泵机械密封的要求也越来越高。文献［２－４］通过对机械密封失效分析发现，机械密封产生泄漏的主要原因中包括：①机械密封工作期间有固体颗粒物质进入到密封摩擦端面，导致摩擦端面过早磨损。②波纹管和小弹簧因为介质颗粒的聚合失去弹性，无法补偿轴向的位移。
对于介质中含有固体颗粒的机械密封，通常采用适当的冲洗方案来减少固体颗粒或污染物［５］。目前国内 常 采 用 ＡＰＩ６８２—２０１２《Ｐｕｍｐｓ－ｓｈａｆｔ ＳｅａｌｉｎｇＳｙｓｔｅｍｆｏｒＣｅｎｔｒｉｆｕｇａｌａｎｄＲｏｔａｒｙＰｕｍｐｓ（４ｔｈＥｄｉ－ｔｉｏｎ）》［６］中的冲洗方案［７］。ＡＰＩ６８２—２０１２中的冲洗方案需要在机械密封外部设置额外的循环管路，在某些不允许改变外部循环管路的工况下，可选择在机械密封内部增加流动改变装置来改善密封端面的工作状况。
目前关于机械密封腔内的数值研究主要集中于温度场、压力场和速度场［８］。任立朝等［９］利用计算流体动力学对复杂工况下机械密封腔内流场进行数值模拟，研究了冲洗口位置对密封腔内流场的影响，并得到了最佳冷却位置。张明明等［１０］的模拟研究结果显示，适当增加冲洗液的流量可以降低密封端面的温度，而且密封端面温度随介质压力或旋转速度的增大而升高。目前对于密封腔内含有固体颗粒的固－液两相流的研究相对较少。Ａｚｉｂｅｒｔ等［１１］运用ＦＬＵＥＮＴ 软件模拟了不同工况条件下机械密封腔内固体粒的分布情况，设计了一种锥形的旋转结构，减少了固体颗粒在密封端面的分布。
文中运用 ＦＬＵＥＮＴ 软件对机械密封介质侧密封腔内固－液两相流场进行数值模拟，分析弹簧和泵送环对密封端面处颗粒流动速度的影响以及颗粒体积分数的分布状态。
１ 控制方程
流 体 动 力 学 的 基 本 控 制 方 程 包 括 质 量 守 恒 方程、动量守恒方程和能量守恒方程［１２］。连续性方程即质量守恒方程，按照质量守恒定律，单位时间内流出控制体的质量应该等于同样时间内控制体内减少的质量，由此可推导出流体流动连续性方程的微分形式为：
动量方程可描述为，对于一给定的流体微元体，其动量对时间的变化率等于外界作用在该微元体上的各种力之和。根据这一定律，可导出ｘ、ｙ、ｚ 这３个方向上的动量方程为：
２ 机械密封腔内流体三维计算模型
２．１ 物理模型建立与网格划分
机械密封介质侧密封腔 结 构 示 意 图 见 图１，图中动环与静环之间为密封摩擦端面，弹簧可以补偿泵轴轴向窜动的位移量。
由于研究对象为流体，因此将机械密封腔内流体区域的模型提取出来。为便于观察，将模型沿轴线半剖，得到的几何模型见图２ａ。流体区域为封闭结构，通过弹簧和泵送环的搅动作用改变颗粒的速 度与体积分数。
进行网格划分时，需要对流域中旋转的部件和壁面进行旋转交界面设置。交 界 面 越 靠 近 旋 转 部件，旋转部件对周围流体区域的影响越小。交界面越靠近壁面，旋转部件对周围流体区域的影响越大。
因此将交界面位置选择在壁面与旋转部件之间，靠近旋转部件的１／４处［１４］。由于模型较为复杂，故采用非结构网格进行划分，划分网格数１４４４７７８，节点数２６７１６９（图２ｂ）。
２．２ 数值计算模型和求解方法选择
机械密封腔内液体与固体颗粒的流动属于多相流，欧拉模型 计 算 结 果 精 度 相 对 较 高［１５］，故 选 择 欧拉模型为数值计算模型。旋转作用下机械密封腔内的流动较为复杂，在计算时作如下假设：①将液相介视为不可压缩流体，且每一相的物理特性为常数，不发生改变。②视固相颗粒为球形，不考虑相变。③液相和固相皆为物性不变的连续相。
模拟时，泵轴的转速不发生改变。机械密封腔内的流动属于定常流动，选择标准ｋ－ε 湍流模型和ＳＩＭＰＬＥ算法进行定常数值计算，使用 ＰｈａｓｅＣｏｕ－ ｌｅｄＳＩＭＰＬＥ 算法求解压力－速度耦合方程组。
２．３ 边界条件设置
机械密封腔内流场数值模拟的边界条件包括泵送环壁面、弹簧壁面、轴套和叶轮的壁面以及密封腔体的内壁面。其中泵送环壁面和弹簧壁面与交界面形成１个圆柱形的转动区域，而其他区域则为静止区域。这样在动静区域之间存在一组动静耦合的交界面，２个区域通过交界面进行数据交换，可以通过多重参考系法（ＭＲＦ）来描述。
将轴套壁面和叶轮壁面设置为运动边界，将密封腔体的内壁面设置为静止边界。对于所有的壁面，均采用无滑移壁面边界条件［１６－１８］。
２．４ 介质物性
液体介质为精丙烯酸，其密度为１０５０ｋｇ／ｍ３，黏度为０．００１１４９Ｐａ·ｓ。
固体颗粒选用的是丙烯酸树脂，其颗粒密度为２１７０ｋｇ／ｍ３，颗粒直径为０．０５ｍｍ。
３ 机械密封腔内三维流场模拟结果与分析
将模型导入 ＦＬＵＥＮＴ 软件，设置泵轴转速为１５００ｒ／ｍｉｎ，重力加速度为９．８ｍ／ｓ２，方向为ｙ 轴的负方向。机械密封腔内的介质为液相丙烯酸和固相丙烯酸树脂颗粒的混合物，设置丙烯酸树脂颗粒的体积分数为１０％。分别在无弹簧和泵送环、有弹簧无泵送环及有弹簧和泵送环这３种条件下对机械密封腔内三维流场进行数值模拟。由于模型较为复杂，为了便于观察，分别取密封端面所在平面和与轴线平行的平面位置进行结果分析。
３．１ 固体颗粒流场
不同条件下机械密封腔内固体颗粒流线分布见图３。从图３ａ可以看出，在 叶 轮 旋 转 作 用 下，靠近叶轮侧的机械密封腔体内形成１个较大的环流，而其他位置的颗粒都围绕在密封端面附近流动，颗粒聚集在密封端面附近。从图３ｂ可以看出，机械密封腔内颗粒的流线形成３个环流，最左侧靠近密封端面处形成１个较小的环流，在该环流的作用下，密封端面处的颗粒可以沿着该环流循环流动，从而减少颗粒在此处的聚集。中间位置处的环流主要是弹簧的旋转产生的，并且该环流能够带走部分左侧小环流内的颗粒。右侧最大的环流则是由叶轮的旋转产生的，由于叶轮的径向尺寸较大，所以在叶轮 边缘位置颗粒的速度最大，达到了１９．６ｍ／ｓ。从图３ｃ可以看出，左侧密封 端 面 处没 有 形 成 环 形 流 域，且在泵送环的影响下，中间的环流范围较图３ｂ中的环流范围大。在增加了泵送环后，颗粒不容易在密封端面处堆 积 且 更 容 易 流 向 叶轮侧的机械密封腔内部。
３．２ 固体颗粒速度
３．２．１ 分布云图
在３种不同条件下机械密封腔内固体颗粒轴向速度场分布 云 图 见 图４。在 速 度 场 中，正 值 代 表 流速方向为ｘ 轴的正方向，负值代表流速方向为ｘ 轴的负方向。
对比图４ａ和图４ｂ可以发现，弹簧的旋转对流体有一定的推动作用，且最大速度出现在弹簧座与壁面的间隙处，此处径向尺寸减小所以流速增加。
图４ｃ中最大流速出现在泵送环的齿形间隙内，且该位置处颗 粒 的 轴 向 速 度 高 于 图 ４ｂ 中 颗 粒 的 轴 向速度。
３．２．２ 变化曲线
不同条件下机械密封密封端面所处平面流场在ｘ＝０．００１５ｍ、ｚ＝０ ｍ 且 ｙ 从 ０．０４１５ ｍ 变化到０．０６２５ｍ直线上的颗粒速度变化曲线见图５，机械密封轴向平面流场在ｙ＝０．０５５ｍ、ｚ＝０ｍ 且ｘ 从０ｍ变化到０．１４６ｍ 直线上的颗粒轴向速度变化曲线见图６。
从图５和图６可以看出，在无弹簧和泵送环的条件下，颗粒的速度基本在０．０１ｍ／ｓ左右，由于没有弹簧和其他旋转零件的影响，颗粒的速度很低且变化很小。增加了弹簧之后，颗粒的速度明显提升。
在密封端面处颗粒的轴向速度最高达０．０３ｍ／ｓ、径向速度最高为０．０４ｍ／ｓ。在有泵送环的流场中，其密封端面上颗粒轴向速度最大值达到０．０８ｍ／ｓ，相较于没有泵送环的最大流速增加了０．０５ｍ／ｓ。在其他位置处，有泵送环时颗粒的速度基本都大于没有泵送环时颗粒的速度。有泵送环时颗粒的径向速度最大值为０．１１ｍ／ｓ，比无泵送环时同一点处的速度增加０．０７ｍ／ｓ。轴向平面上泵送环产生的影响更明显， 泵送环使颗粒的最高速度达到０．５５ｍ／ｓ， 而没有泵送环时颗粒的最高速度为０．２ｍ／ｓ。在有泵送环和无泵送环条件下，颗粒速度的最高值出现在０．０５～０．０７５ｍ这段距离，该段处于泵送环与密封腔体内壁面之间，由于此处间隙相对较小，且泵送环对流体具有推动作用，故颗粒速度在此处达到最大值。当轴向距离远离泵送环时，泵送环对颗粒的作用逐渐减小，颗粒的轴向速度也随之降低并最终趋于０。
３．３ 固体颗粒体积分数
３．３．１ 分布云图
不同条件下机械密封轴向平面颗粒体积分数分布云图见图７。
从图７ａ可以看出，由于没有弹簧和泵送环的旋转搅动作用且受自身重力的影响，颗粒最终聚集在密封腔的底部。从图７ｂ可以看出，密封腔内的颗粒体积分数在０．１左右。弹簧在自身的旋转作用下产生一定的搅动作用，弹簧中间段颗粒的体积分数减小到约０．０３。在下方的锥形壁面上，由于弹簧对流体的推动作用不够大，同时在重力作用下，颗粒体积分数会有所升高。结合前述对颗粒速度的分析可知，在泵送环的作用下，密封腔内的液体被泵送环推送至密封端面附近，对密封端面形成一定的冲洗效果，且在冲洗过程中可以带走弹簧两端的颗粒，使弹簧所处位置的颗粒体积分数进一步减小（图７ｃ）。
不同条件下机械密封密封端面所在平面颗粒体积分数分布云图见图８。泵轴及轴上的旋转件在高速旋转的过程中会产生一定的离心力，减少颗粒物质在密封端面上的附着。比较图８ｂ和图８ｃ可以发现，在离心力作用下，颗粒体积分数随密封腔径向尺寸的增大而增大，在泵送环影响下，该现象更明显。
３．３．２ 变化曲线
不同条件下机械密封密封端面所在平面在ｘ＝０．００１５ｍ、ｚ＝０ ｍ，并 且 ｙ 从 ０．０４１５ ｍ 变 化 到０．０６２５ｍ直线上的颗粒体积分数变化曲线见图９。在无弹簧和泵送环条件下，颗粒的体积分数变化趋势明显不同于其他２种情况，其颗粒体积分数呈逐渐减小的趋势，在最靠近密封端面处，体积分数达到０．０９，为３种条件中的最高值。在有弹簧的条件下，颗粒体积分数随着径向尺寸的增大而逐渐增大。对于没有泵送环的情况，靠近密封端面处颗粒的体积分数为０．０５，壁面处颗粒的体积分数约０．１３。结合前述对颗粒速度分析的结果可知，增加泵送环后，泵送环提高了颗粒的径向速度，使颗粒在密封端面处的体 积 分 数 趋 于 ０，且 在 ｙ＝０．０４１５ ｍ 至 ｙ＝０．０５５ｍ这段距离内，颗粒体积 分 数 一 直 小 于 没 有泵送环的情况。随着径向尺寸的增加，颗粒体积分数逐渐升高到０．２左右。
不同条件下机械密封轴向平面在ｙ＝０．０５５ｍ、ｚ＝０，ｘ从０变化到０．１４６ｍ 直线上的颗粒体积分数变化曲线见图１０。从图１０可以看出，在该直线上，有泵送环情况下的颗粒体积分数基本高于没有泵送环的情况。在泵送环的推动和旋转作用下，固体颗粒更多分布在靠近壁面的位置以及右侧的腔体内。在无弹簧和泵送环条件下，颗粒的体积分数最 小。在有弹簧和泵送环条件下，颗粒体积分数介于其他２种情况之间。这种现象说明弹簧及泵送环旋转产生的离心力使颗粒围绕密封腔的内壁面均匀分布，且密封端面处颗粒的体积分数在离心力和搅动作用下明显减小。
４ 结语
运用多重参考系法和欧拉模型对机械密封腔内流场进行了三维数值模拟，模拟结果显示，在无弹簧和泵送环的条件下，密封腔内的颗粒主要受重力的影响，容易 附 着 在 密 封 端 面 处 以 及 密 封 腔 的 底 部。
增加弹簧后，弹簧旋转对流体产生一定的推动作用，密封端面的颗粒明显减少，且颗粒的轴向和径向速度增加，有利于阻止颗粒进入密封端面，从而减少密封端面的磨损。增加泵送环后，密封端面处的颗粒体积分数进一步降低，固体颗粒在弹簧和泵送环旋转及推动作用下向壁面处和右侧腔体靠近。
在机械密封的设计和改进中，需要充分考虑密封摩擦副端面附近颗粒的流动状态和分布情况，可以通过增加额外的流动改变装置来减小密封端面附近颗粒的体积分数。
参考文献：
［１］ 顾永泉．机械密 封 实 用 技 术［Ｍ］．北 京：机 械 工 业 出 版社，２００１．ＧＵ Ｙ Ｑ．Ｐｒａｃｔｉｃａｌｔｅｃｈｎｏｌｏｇｙ ｏｆ ｍｅｃｈａｎｉｃａｌｓｅａｌｓ［Ｍ］．Ｂｅｉｊｉｎｇ：ＣｈｉｎａＭａｃｈｉｎｅＰｒｅｓｓ，２００１．
［２］ 李国富，徐龙甫，高 洪 志．机 械 密 封 失 效 原 因 分 析［Ｊ］．设备管理与维修，２０１６，４７（Ｓ１）：９１－９２．ＬＩＧＦ，ＸＵ ＬＦ，ＧＡＯ Ｈ Ｚ．Ｆａｉｌｕｒｅａｎａｌｙｓｉｓｏｆｍｅ－ｃｈａｎｉｃａｌｓｅａｌｓ［Ｊ］．Ｐｌａｎｔ ｍｉｎｔｅｎａｎｃｅ ｅｎｇｉｎｅｅｒｉｎｇ，２０１６，４７（Ｓ１）：９１－９２．
［３］ 来永斌，陈秀．机械密封失效的原因分析及装配［Ｊ］．润滑与密封，２００５（２）：１６３－１６４．
ＬＡＩＹ Ｂ，ＣＨＥＮ Ｘ．Ｍｅｃｈａｎｉｃａｌｓｅａｌｆａｉｌｕｒｅａｎａｌｙｓｉｓａｎｄａｓｓｅｍｂｌｙ［Ｊ］．Ｌｕｂｒｉｃａｔｉｏｎｅｎｇｉｎｅｅｒｉｎｇ，２００５（２）：１６３－１６４．
［４］ 李强，高海容．机械密封的密封失效原因分析［Ｊ］．煤矿机械，２００６（４）：６８５－６８６．
ＬＩＱ，ＧＡＯ ＨＲ．Ｍｅｃｈａｎｉｃａｌｓｅａｌｓｅａｌｆａｉｌｕｒｅｒｅａｓｏｎａ－ｎａｌｙｓｉｓ［Ｊ］．Ｃｏａｌｍｉｎｅｍａｃｈｉｎｅｒｙ，２００６（４）：６８５－６８６．
［５］ 张营营．机械 密 封 冲 洗 方 案［Ｊ］．化 工 设 计 通 讯，２０１８，４４（７）：１３０－１３１．
ＺＨＡＮＧＹ Ｙ．Ｍｅｃｈａｎｉｃａｌｓｅａｌｆｌｕｓｈｉｎｇｓｏｌｕｔｉｏｎ［Ｊ］．Ｃｈｅｍｉｃａｌｅｎｇｉｎｅｅｒｉｎｇｄｅｓｉｇｎｃｏｍｍｕｎｉｃａｔｉｏｎｓ，２０１８，４４（７）：１３０－１３１．
［６］ Ｐｕｍｐｓ－ｓｈａｆｔｓｅａｌｉｎｇｓｙｓｔｅｍｆｏｒｃｅｎｔｒｉｆｕｇａｌａｎｄｒｏｔａｒｙｐｕｍｐｓ（４ｔｈｅｄｉｔｉｏｎ）：ＡＰＩ６８２—２０１２［Ｓ］．
［７］ 蔡仁良，陈伟．ＡＰＩ６８２标准的新进展———第４版的主要变化［Ｊ］．化工设备与管道，２０１４，５１（４）：１－１０．
ＣＡＩＲＬ，ＣＨＥＮ Ｗ．ＮｅｗｐｒｏｇｒｅｓｓｏｆＡＰＩ６８２—Ｐｒｉ－ｍａｒｙｃｈａｎｇｅｓｉｎ４ｔｈｅｄｉｔｉｏｎ［Ｊ］．Ｐｒｏｃｅｓｓｅｑｕｉｐｍｅｎｔ＆ｐｉｐｉｎｇ，２０１４，５１（４）：１－１０．
［８］ 赵龙，宋鹏云，王 红 利．机 械 密 封 腔 流 场 实 验 和 数 值 模拟研究进展［Ｊ］．机械，２００８，３５（１０）：８－１１．
ＺＨＡＯＬ，ＳＯＮＧ Ｐ Ｙ，ＷＡＮＧ Ｈ Ｌ．Ｐｒｏｇｒｅｓｓｉｎｔｈｅｅｘｐｅｒｉｍｅｎｔａｌａｎｄｎｕｍｅｒｉｃａｌｓｔｕｄｙｏｆｔｈｅｆｌｏｗｉｎｓｅａｌｃｈａｍｂｅｒ［Ｊ］．Ｍａｃｈｉｎｅｒｙ，２００８，３５（１０）：８－１１．
［９］ 任立朝，穆塔里夫·阿赫迈德，高强，等．复杂工况机械密封腔内两相流 的 流 场 特 性 分 析 ［Ｊ］．机 械 设 计 与 制造，２０１３（１）：６９－７１．
ＲＥＮＬＣ，ＡｈｍａｔＭｕｔｅｌｌｉｐ，ＧＡＯ Ｑ，ｅｔａｌ．Ｔｗｏ－ｐｈａｓｅｆｌｏｗｆｉｅｌｄｃｈａｒａｃｔｅｒｉｓｔｉｃｓａｎａｌｙｓｉｓｏｎｓｅａｌｃａｖｉｔｙｆｏｒｍｅ－ｃｈａｎｉｃａｌｓｅａｌｕｎｄｅｒｃｏｍｐｌｅｘｃｏｎｄｉｔｉｏｎｓ［Ｊ］．Ｍａｃｈｉｎｅｒｙｄｅｓｉｇｎ＆ ｍａｎｕｆａｃｔｕｒｅ，２０１３（１）：６９－７１．
［１０］张明明，郝木明．机械密封腔内流场及摩擦副温度场的数值计算［Ｊ］．石油化工设备技术，２００８，２９（１）：６３－６６．
ＺＨＡＮＧ Ｍ Ｍ，ＨＡＯ Ｍ Ｍ．Ｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｏｆｆｌｕｉｄｆｉｅｌｄａｎｄｔｅｍｐｅｒａｔｕｒｅｆｉｅｌｄｏｆｆｒｉｃｔｉｏｎｉｎｔｅｒｆａｃｅｉｎａｍｅｃｈａｎｉｃａｌｓｅａｌ［Ｊ］．Ｐｅｔｒｏ－ｃｈｅｍｉｃａｌｅｑｕｉｐｍｅｎｔｔｅｃｈｎｏｌ－ｏｇｙ，２００８，２９（１）：６３－６６．
［１１］ＡｚｉｂｅｒｔＨ，ＣｌａｒｋＲ．ＵｓｉｎｇＣＦＤｔｏｉｍｐｒｏｖｅｐｅｒｆｏｒｍ－ａｎｃｅａｎｄｅｘｔｅｎｄｌｉｆｅｏｆｍｅｃｈａｎｉｃａｌｓｅａｌｓｉｎｓｌｕｒｒｙａｐｐｌｉ－ｃａｔｉｏｎｓ［Ｃ］／／１７ｔｈＩｎｔｅｒｎａｔｉｏｎａｌＣｏｎｆｅｒｅｎｃｅｏｎ ＦｌｕｉｄＳｅａｌｉｎｇ．Ｙｏｒｋ，ＵＫ：ＢＨＲＧｒｏｕｐ，２００３．
［１２］王福军．计算 流 体 动 力 学 分 析———ＣＦＤ 软 件 原 理 与 应用［Ｍ］．北京：清华大学出版社，２００４：４６８．ＷＡＮＧ Ｆ Ｊ．Ｃｏｍｐｕｔａｔｉｏｎａｌｆｌｕｉｄ ｄｙｎａｍｉｃｓ ａｎａｌｙ－ｓｉｓ———ＰｒｉｎｃｉｐｌｅａｎｄａｐｐｌｉｃａｔｉｏｎｏｆＣＦＤｓｏｆｔｗａｒｅ［Ｍ］．Ｂｅｉｊｉｎｇ：ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙＰｒｅｓｓ，２００４：４６８．
［１３］张明明，黄琦．基 于 ＦＬＵＥＮＴ 的 机 械 密 封 腔 内 流 场 数值计算［Ｊ］．通用机械，２００８（３）：４０－４２，４７．
ＺＨＡＮＧ Ｍ Ｍ，ＨＵＡＮＧ Ｑ．Ｎｕｍｅｒｉｃａｌｃａｌｃｕｌａｔｉｏｎｏｆｆｌｏｗｆｉｅｌｄｉｎ ｍｅｃｈａｎｉｃａｌｓｅａｌｃｈａｍｂｅｒｂａｓｅｄｏｎＦＬＵ－ＥＮＴ［Ｊ］．Ｇｅｎｅｒａｌｍａｃｈｉｎｅｒｙ，２００８（３）：４０－４２，４７．
［１４］王文全，郝栋伟，张 立 翔．滑 移 界 面 位 置 对 叶 轮 机 械 内部旋转流场 的 影 响［Ｊ］．北 京 理 工 大 学 学 报，２０１２（９）：８９０－８９４．
ＷＡＮＧ Ｗ Ｑ，ＨＡＯ Ｄ Ｗ，ＺＨＡＮＧ Ｌ Ｘ．Ｅｆｆｅｃｔｓｏｆｓｌｉｄｉｎｇｍｅｓｈｉｎｔｅｒｆａｃｅｐｏｓｉｔｉｏｎｏｎｓｗｉｒｌｉｎｇｆｌｏｗｆｉｅｌｄｓｉｎｉｍｐｅｌｌｅｒｍａｃｈｉｎｅｒｙ［Ｊ］．ＴｒａｎｓａｃｔｉｏｎｓｏｆＢｅｉｊｉｎｇｉｎｓｔｉ－ｔｕｔｅｏｆｔｅｃｈｎｏｌｏｇｙ，２０１２（９）：８９０－８９４．
［１５］张锦洲，高陈，易先中，等．基于十字搅拌器的湿煤灰搅拌数值模拟［Ｊ］．长江大学学报（自然科学版），２０１８，１５（９）：６－９．
ＺＨＡＮＧＪＺ，ＧＡＯＣ，ＹＩＸＺ，ｅｔａｌ．Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａ－ｔｉｏｎｏｆｗｅｔｃｏａｌａｓｈｓｔｉｒｒｉｎｇｂａｓｅｄｏｎｃｒｏｓｓｓｔｉｒｒｅｒ［Ｊ］．ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅｕｎｉｖｅｒｓｉｔｙ（Ｎａｔｕｒａｌｓｃｉｅｎｃｅｅｄｉｔｉｏｎ），２０１８，１５（９）：６－９．
［１６］牟童，郑源，郑建坤，等．Ｔ 型搅拌器数值模拟研究［Ｊ］．南水北调与水利科技，２０１４（３）：８９－９４．
ＭＯＵ Ｔ，ＺＨＥＮＧ Ｙ，ＺＨＥＮＧＪＫ，ｅｔａｌ．ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆＴ－ｓｈａｐｅｄｓｔｉｒｒｅｒ［Ｊ］．Ｓｏｕｔｈ－ｔｏ－ｎｏｒｔｈｗａ－ｔｅｒｔｒａｎｓｆｅｒｓａｎｄ ｗａｔｅｒｓｃｉｅｎｃｅ ＆ ｔｅｃｈｎｏｌｏｇｙ，２０１４（３）：８９－９４．
［１７］莫惠珺，范海涛，王明玥，等．型式结构对双曲面搅拌器流场影响的数值模拟［Ｊ］．环境工程学报，２０１９，１３（５）：１１１９－１１２７．
ＭＯ ＨＪ，ＦＡＮ ＨＴ，ＷＡＮＧ Ｍ Ｙ，ｅｔａｌ．Ｆｌｏｗｆｉｅｌｄｎｕ－ｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｈｙｐｅｒｂｏｌｉｃ ｍｉｘｅｒｗｉｔｈｄｉｆｆｅｒｅｎｔｓｔｒｕｃｔｕｒｅｓ［Ｊ］．Ｃｈｉｎｅｓｅｊｏｕｒｎａｌｏｆｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉ－ｎｅｅｒｉｎｇ，２０１９，１３（５）：１１１９－１１２７．
［１８］栾晓坤，郑深晓，张盛峰，等．错位六弯叶桨搅拌假塑性流体流场特 性 数 值 模 拟 ［Ｊ］．石 油 化 工 设 备，２０１６，４５（６）：１３－１９．
ＬＵＡＮＸＫ，ＺＨＥＮＧＳＸ，ＺＨＡＮＧＳＦ，ｅｔａｌ．Ｎｕｍｅｒ－ｉｃａｌｓｉｍｕｌａｔｉｏｎｏｎｆｌｏｗｆｉｅｌｄｆｅａｔｕｒｅｓｏｆｐｓｅｕｄｏｐｌａｓｔｉｃｆｌｕｉｄｓｓｔｉｒｒｅｄ ｂｙｔｈｅｉｍｐｅｌｌｅｒｏｆｐｅｒｔｕｒｂｅｄｓｉｘ－ｂｅｎｔ－ｂｌａｄｅｄｔｕｒｂｉｎｅ［Ｊ］．Ｐｅｔｒｏ－ｃｈｅｍｉｃａｌｅｑｕｉｐｍｅｎｔ，２０１６，４５（６）：１３－１９．
文章属扬州市白浦江气动件有限公司原创，版权归其所有，转载请注明出处：http://www.zhichenggan.com，欢迎电询本司产品气弹簧│不锈钢气弹簧│压缩气弹簧│可控气弹簧。