高温环境下料浆输送管流速度场特性研究.docxVIP

高温环境下料浆输送管流速度场特性研究 自充采法具有控制地压、减少废物排放和环境保护等优点以来，已在近年来得到了广泛的研究和应用。 为保证安全、高效地进行料浆输送，众多专家学者进行了深入研究 本研究在获取高温条件下料浆流变参数的基础上，利用COMSOL软件模拟高温条件下L管大流量输送过程中管内速度场的分布特征，得出高温环境下料浆输送最佳管径、初始流速等参数取值，为矿山充填管道输送设计提供借鉴，同时也为进一步分析温度变化及大流量输送料浆对管流特性的影响提供理论依据。 1 试验条件 1.1 料浆流场模拟 本研究采用COMSOL软件进行模拟分析，在料浆流动过程中进行以下假设 (1）本次数值模拟只考虑料浆自身温度，选取40℃、50℃、60℃作为模拟温度的变量，并且外界环境温度设置成室温25℃，流动过程中不与外界发生热交换。 (2）模拟过程将料浆视为不可压缩流体。 (3）料浆在管道中处于静止或者流动状态时无间隙，即料浆是连续的。 (4）流动过程满足常态稳定流，且各向力学特征一致。 1.2 黏度及屈服切应力 充填料浆流变模型符合宾汉模型 式中，τ为管道壁面的剪切应力，Pa；τ 计算了不同温度下充填料浆的塑性黏度以及屈服切应力，结果如表1所示。由表1可知：随着温度增加，充填料浆的塑性黏度随之减小，屈服切应力随之降低，但并非呈线性变化，在40℃时达到最大塑性黏度和屈服切应力（0.57 Pa·s和202.63 Pa），在60℃时达到最小塑性黏度和屈服切应力（0.54Pa·s和162.93 Pa）。 1.3 纳维-斯托克斯方程 在充填料浆输送过程中，一般将其认定为不可压缩流体，流动过程中应遵循能量方程、连续性方程以及动量方程，且料浆密度保持不变 (1）连续方程。方程表达式为 式中，ρ为密度，kg/m3；t为时间，s；??Vue04c为速度散度。 (2）动量方程。纳维-斯托克斯方程（NavierStokes Equation）是描述黏性不可压缩流体动量守恒的运动方程，简称N-S方程。其矢量形式为 式中，p为流体压力，Pa；Δ为速度梯度；F为流体上的质量力，N；?p为压力散度；μ为动力黏度，Pa·s。 (3）能量方程。充填料浆流动过程中的能量遵循伯努利方程，其表达式为 式中，z、z 1.4 管道直径及初始速度设计 本研究利用COMSOL软件建立输送L管模型，竖直管道长3 m，水平管道长9 m，直径选取200 mm、240 mm，料浆进口为竖直管道上口，初始速度选取2.5 m/s、3 m/s，出口为水平管道末端，具体参数如表2所示。 2 结果和分析 2.1 料浆流场分布调整 为分析管道输送过程中料浆流态变化规律，选取L管模型进行模拟，研究受重力影响下料浆流态分布。图1表示温度为40℃、50℃、60℃时L管模型的料浆速度变化分布情况。 由图1可知：不同温度下料浆流动速度云图变化基本相似，在竖直管道处，呈现结构流向下流动，在弯管处表现为塞流向前流动；自弯管转向水平管后，料浆速度分布发生改变，出现管道内侧速度变小、外侧速度变大的趋势，料浆内部的颗粒在重力影响下，粒径较小的颗粒会在外侧管道处形成一层能够促进料浆向前流动的润滑层，从而使其流动速度加快；在内侧流动的料浆颗粒之间，由于要改变运动轨道自弯管而进入水平管处，在流动过程中会发生比较剧烈的碰撞，并且颗粒与管道也会存在较大的摩擦，从而会加大颗粒与颗粒、颗粒与管道的作用力，阻碍料浆向前流动，导致流速呈现减小趋势。 由于受到重力影响，料浆在竖直管道处以初始速度和重力作为动力向下流动，各层料浆之间产生相互作用摩擦，但受到重力作用均匀，流态稳定，料浆经弯管流入水平管道后，重力不再是料浆向前流动的动力，颗粒间挤压作用逐渐变小，料浆整体向前流动。因此在流动一段时间后，料浆速度分布又会慢慢恢复成塞流形式，即中心存在一个速度最大的流核区，然而温度的改变不会对流态变化带来严重的影响，流动比较稳定，不会出现湍流区域。 2.2 温度对管速度剖面的影响 为研究稳定状态下料浆某一断面速度受温度的影响，选择7.5 m处断面进行研究分析。不同温度下料浆7.5 m处速度剖面如图2所示，不同温度下料浆速度径向变化如图3所示。 随着温度提高，料浆中心最大流动速度增大，会加剧对管道的磨损，缩短管道使用年限，本研究建议料浆输送过程控制温度为40～50℃，既可保证输送流速，又能减小料浆与管道的摩擦阻力损失。 2.3 管道直径、孔径 为分析稳定状态下管径变化对料浆流速的影响，选取管道直径分别为200 mm和240 mm进行模拟研究，结果如图4和图5所示。 由图5可知：管道直径为200 mm时，料浆输送在管道径向0.05 m处开始出现最大流核区域，在0.15 m处达到区域最大化，中心最大流核区径向长度为0.10 m；管道直径为240 mm时，料浆输送在管