输气管道弯头颗粒侵蚀数值模拟方法比较与选型

文献概况与研究意义

本文由中国海洋大学工程学院张日、朱东、董胜三位学者合著，发表于2020年《船舶力学》第24卷第10期，受科技部国家重点研发计划课题和中央高校基本科研业务费专项资助。研究针对输气管道弯头处固体颗粒侵蚀问题，系统比较了不同湍流模型、颗粒追踪方法和冲蚀模型的计算精度，为工程实践中侵蚀预测的数值方法选型提供了科学依据。

输气管道中携带的固体颗粒（如砂粒、粉尘、水合物碎屑等）在弯头处的侵蚀是管道完整性管理的重要问题。弯头处由于流场急剧变化，颗粒撞击管壁的频率和能量显著增加，导致管壁减薄甚至穿孔。数值模拟方法可以预测侵蚀分布和速率，为管道设计、维护和更换提供依据。

湍流模型选择与流场模拟

文献使用五种不同的湍流模型计算输气弯头处的速度场，通过与实验数据对比评估各模型的精度：

研究结果表明，剪切应力传输（SST）k-ω模型是计算输气弯头内流场的最优湍流模型。这是因为SST模型在近壁区采用k-ω模型的形式，在远场采用k-ε模型的形式，能够更好地捕捉弯头处的流动分离和二次流特征。弯头处的流场具有明显的逆压梯度和流动分离，SST模型对这类流动的描述能力优于其他模型。

颗粒追踪与冲蚀模型比较

颗粒追踪方法

文献比较了离散相与连续相的单向耦合和双向耦合方法：

• 单向耦合：颗粒运动受流场影响，但颗粒对流场无反馈作用。适用于颗粒浓度低、质量分数小的情况。
• 双向耦合：颗粒运动受流场影响，同时颗粒对流场产生反馈作用。适用于颗粒浓度较高、对气流产生显著影响的情况。

研究结果表明，双向耦合模型对于颗粒运动状态的预测更为精确。这是因为在输气管道弯头处，颗粒浓度可能较高，颗粒与气流之间的动量交换不可忽略，单向耦合会低估颗粒对流场的扰动，导致颗粒轨迹预测偏差。

冲蚀模型比较

文献使用四种颗粒冲蚀模型计算弯头侵蚀，并与实验数据进行对比：

研究结果表明，Oka颗粒冲蚀模型能够更为准确地计算各流场条件下的固体颗粒侵蚀结果。Oka模型综合考虑了颗粒撞击速度、撞击角度、靶材硬度和颗粒硬度等因素，其经验公式拟合精度较高，适用于金属管道的侵蚀预测。

相似性准则与工程应用

文献还基于相似性准则，计算了不同模型尺度下侵蚀速率的无量纲数，验证了数值方法的尺度不变性。研究表明，固体颗粒侵蚀的数值方法符合相似性准则，在低输沙率条件下计算精度更高。

这一发现对工程实践具有重要指导意义：

1. 模型验证策略：在进行全尺寸管道的侵蚀预测之前，可以先在缩尺模型上进行实验验证，通过无量纲数进行尺度转换，降低实验成本。
2. 低输沙率适用性：在低输沙率条件下，颗粒之间的相互作用可忽略，数值模拟的精度较高，预测结果更可靠。高输沙率条件下，颗粒-颗粒碰撞和颗粒群效应需要更复杂的模型来描述。
3. 工程应用建议：对于新建管道的侵蚀预测，建议采用SST k-ω湍流模型加双向耦合颗粒追踪加Oka冲蚀模型的计算组合，以获得最精确的预测结果。

学习心得

该文献给我最大的启示是：数值模拟方法的选择不是随意的，而是需要通过与实验数据的系统对比来确定最优组合。在实际工程中，很多技术人员在选择湍流模型和冲蚀模型时缺乏依据，往往凭经验或默认设置，导致计算结果不可靠。文献提供的模型比较方法，为我们建立了一套科学的选型流程。

另一个重要启示是相似性准则在工程中的应用价值。在管道完整性管理中，我们经常需要在不同管径、不同工况下进行侵蚀预测，相似性准则为这种尺度转换提供了理论依据。在实际工作中，我们可以利用这一原理，在实验室小尺寸弯头上进行侵蚀实验，然后通过无量纲数转换到全尺寸管道，大幅降低实验成本和时间。

此外，该文献的研究方法也提醒我们，数值模拟的每一步都需要验证：流场模拟需要与实验速度场对比，颗粒轨迹需要与实验轨迹对比，侵蚀结果需要与实验侵蚀深度对比。只有每一步都经过验证，最终的侵蚀预测才具有工程可信度。