1、1流体动力学（CFD）在往复式压缩机管道系统气流脉动计算中的应用摘 要：参照实验室所搭建的管道系统实验平台，根据计算流体动力学（CFD）方法建立管道内气体的二维非定常流动模型。利用流体仿真软件 FLUENT 计算了缓冲器及孔板前后的气流脉动，通过分析气流脉动曲线及流场的分布情形验证了缓冲器及孔板对气流脉动的消减作用。通过对比实验数据验证了利用 CFD 技术研究管道系统气流脉动是准确可靠的。关键词：管道系统 CFD 技术 气流脉动 FLUENT 孔板 往复式压缩机是石油、天然气、化工及电力等工业生产中的重要机械设备，其管道系统又是实现物质运输的主要途径，然而管道系统的振动会对安全生产造成很大的威
2、胁，众多生产实践表明压缩机管路的绝大多数振动问题都是由气流脉动引起的，而压缩机吸排气的间歇性、周期性特点是产生气流脉动的主要原因。因此研究气流脉动的产生机理，建立合理的流体动力学模型进行管道中气流脉动的预测具有重要的理论意义和工程实用价值。 现有研究气流脉动较为成熟的方法大多基于平面波动理论1或一维非定常流动理论2，它们均未考虑流体流动时湍流的影响，同时对缓冲器、孔板、冷却器、分离器等管路元件的气流脉动计算精度也较差。随着计算机速度的提高和近年来 CFD 技术的发展，选用有限元方法3，42及有限容积法5计算管系的气流脉动取得了一定的成效。CFD 方法6应用于稳态的工业流场模拟已有较多的报道，但
3、对非稳态的脉动流场研究较少。 本文基于 CFD 方法建立管道系统流体动力学模型。在考虑湍流的情况下7，模拟了含空冷器及孔板管道等管路原件的管道系统非定常流动时气流脉动及流场特性。通过和实验数据对比验证了 CFD 方法计算管道系统气流脉动的合理性及准确性。 一、CFD 模拟计算理论 目前广泛用于计算流体力学的数值方法有有限差分法、有限元法、有限体积法等，其目的都是将控制方程离散化，本文用到的 CFD 软件FLUENT8-9采用有限体积法将非线性偏微分方程转变为网格单元上的线性代数方程，然后通过求解线性方程组得出流场的解。因此对于所有流动，FLUENT 都求解质量和动量守恒方程；对于包含传热或可压
4、性流动，还需要增加能量守恒方程；如果是湍流问题，还有相应的输运方程需要求解；我们称以上各方程为控制方程。根据模型特点，本文所用 FLUENT中的标准 湍流模型对低速可压管流问题有良好的表现。 1.统一控制方程形式 式中， 为通用变量，可以代表 ， ， ， 等求解变量； 为广义扩散系数； 为广义源项 ，式中各项依次为瞬态项（transient term） 、对流项（convective term） 、扩散项（diffusive term）和源项（source term） 。 2.标准 方程 3标准 模型需要求解湍动能 和耗散率 方程。该模型假设流体流动为完全湍流，忽略分子的粘性影响。该模型的湍动
5、能 和耗散率 方程为： -由于平均速度梯度引起的湍动能； -由于浮力影响引起的湍动能； -可压缩湍流脉动膨胀对总的耗散率的影响； -湍流粘性系数， 。 此模型中的各个常数在 FLUENT 中的默认值为 ， ，湍动能 和 耗散率 的湍流普朗特数分别为： ， 该模型适合对完全湍流流动的模拟。 二、气流压力脉动实验及管道系统计算模型 1.孔板消减气流脉动实验 为了对管道系统气流脉动的模拟结果进行验证，依照图 1 搭建了往复式压缩机管道系统压力脉动测试实验平台。测试系统主要由一台往复式压缩机、变频器、与其相连的管道系统、直流稳压电源、压力传感器、示波器、数据采集及分析系统组成。 2.管道系统的计算模型
6、 根据搭建的管道系统实验平台，在不影响模拟结果的前提下，为提高计算效率采用如图 2 所示（尺寸单位为 mm）的二维计算模型。将模型导入 Gambit 中划分单元网格，网格类型采用非结构化（Pave） 的四边形网格（Quad） ，共得到 17644 个单元。最后在 Gambit 中对管道系统添加边界条件，入口边界条件采用非定常压力入口边界条件，出口边界条件为压力出口边界。 4三、气流脉动模拟及实验分析 1.模拟分析 将在 Gambit 中化分好的网格模型导入 Fluent 中，在 Fluent 中根据管道内部介质的性质采用基于密度（耦合式）的隐式 2ddp（二维双精度）求解器，并指定其计算模式为
7、非稳态；计算模型选择标准 双方程湍流模型；流体类材料为理想空气。根据现场实验条件，出口边界条件采用非定常压力边界条件： （根据其公式用 C 语言编译非定常速度 UDF 函数）即施加一个频率为 7.3Hz（模拟转速为 438rmin-1 压缩机） 、压力不均匀度为 24%的脉动压力条件，入口边界条件采用定常压力边界条件（由储气罐上的背压阀控制） ，文中设定为 101325Pa。图 3 为未加孔板时缓冲器前后压力脉动曲线图，由图可以看出缓冲器对消减气流脉动有一定的作用；图 4 为加入孔径比为 0.5 的孔板后缓冲器前后压力脉动曲线图，对比图 4 发现孔板对消减气流脉动有明显的作用效果，但在加入孔板
8、时需要选择适当的孔径比，由图可知 0.5 的孔径比衰减效果很理想。 2.实验测试结果 图 5 为当压缩机以 438 r？min-1 转速运行时，添加孔径比为 0.5，厚度为 8mm 的孔板后，孔板前后主管线上压力脉动的对比图；与模拟结果 （图 4 所示曲线）十分吻合。从二者对比结果可看出利用 CFD 的方法计算管道系统的脉动压力具有直观、形象、准确的特点。 四、结论 1.基于 CFD 技术，通过利用 FLUENT 软件对管道系统内部流场的计算，5直观形象的验证了缓冲器、孔板等管路原件对气流脉动的消减作用。 2.通过与实验数据的对照，发现利用 CFD 技术计算气流脉动更加可靠、准确。 参考文献
9、1党锡淇，陈守五.活塞式压缩机气流脉动与管道振动M.陕西西安：西安交通大学出版社，1984：7.143. 2陈守五，黄幼玲.往复式压缩机一维不稳定气流方程的数值解法J.西安交通大学学报，1982，1：5566. 3Enzo Giacomelli，Marco Passeri.Control of pulsation and vibrations in hyper compressors for LDPE plants. ASME PVP2004 Pressure vessel & piping conference， San Diego California 2004.7 19-22. 4Ing
10、.Attilio Brighenti，Ing.Andrea Pavan.ACUSCOMP and ACUSYS-A powerful hybrid linear/nonlinear simulation suite to analyses pressure pulsations in pipingM.Italy：SATE Systems and Advanced Technology Engineering，Santa Croce 664/A，301 35 Venezia，2006，7 23-27. 5 苏永生，王恒杰. 应用 CFD 消除气流脉动J.华东理工大学学报（自然科学版） ， 2006， 32 （4） ： 480 - 483. 6 任玉新，陈海昕.计算流体力学基础M.北京：清华大学出版社，2006：1-11. 7 王福军.计算流体动力学分析-CFD 软件原理与应用M.北京：清6华大学出版社，2004. 8 韩占忠，王敬，兰小平.流体工程仿真计算实例与应用.北京理工大学出版社，2004. 9 温正，石良臣，任毅如.FLUENT 流体计算应用教程.清华大学出版社.2009. 10王济，胡晓.MATLAB 在振动信号处理中的应用M.北京：中国水利水电出版社，2006：112-117.
