流体动力学在往复式压缩机管道系统气流脉动计算中应用.docVIP

流体动力学在往复式压缩机管道系统气流脉动计算中应用 　　摘 要：参照实验室所搭建的管道系统实验平台，根据计算流体动力学（CFD）方法建立管道内气体的二维非定常流动模型。利用流体仿真软件FLUENT计算了缓冲器及孔板前后的气流脉动，通过分析气流脉动曲线及流场的分布情形验证了缓冲器及孔板对气流脉动的消减作用。通过对比实验数据验证了利用CFD技术研究管道系统气流脉动是准确可靠的。 　　关键词：管道系统 CFD技术 气流脉动 FLUENT 孔板 　　往复式压缩机是石油、天然气、化工及电力等工业生产中的重要机械设备，其管道系统又是实现物质运输的主要途径，然而管道系统的振动会对安全生产造成很大的威胁，众多生产实践表明压缩机管路的绝大多数振动问题都是由气流脉动引起的，而压缩机吸排气的间歇性、周期性特点是产生气流脉动的主要原因。因此研究气流脉动的产生机理，建立合理的流体动力学模型进行管道中气流脉动的预测具有重要的理论意义和工程实用价值。 　　现有研究气流脉动较为成熟的方法大多基于平面波动理论[1]或一维非定常流动理论[2]，它们均未考虑流体流动时湍流的影响，同时对缓冲器、孔板、冷却器、分离器等管路元件的气流脉动计算精度也较差。随着计算机速度的提高和近年来CFD技术的发展，选用有限元方法[3，4]及有限容积法[5]计算管系的气流脉动取得了一定的成效。CFD方法[6]应用于稳态的工业流场模拟已有较多的报道，但对非稳态的脉动流场研究较少。 　　本文基于CFD方法建立管道系统流体动力学模型。在考虑湍流的情况下[7]，模拟了含空冷器及孔板管道等管路原件的管道系统非定常流动时气流脉动及流场特性。通过和实验数据对比验证了CFD方法计算管道系统气流脉动的合理性及准确性。 　　一、CFD模拟计算理论 　　目前广泛用于计算流体力学的数值方法有有限差分法、有限元法、有限体积法等，其目的都是将控制方程离散化，本文用到的CFD软件FLUENT[8-9]采用有限体积法将非线性偏微分方程转变为网格单元上的线性代数方程，然后通过求解线性方程组得出流场的解。因此对于所有流动，FLUENT都求解质量和动量守恒方程；对于包含传热或可压性流动，还需要增加能量守恒方程；如果是湍流问题，还有相应的输运方程需要求解；我们称以上各方程为控制方程。根据模型特点，本文所用FLUENT中的标准 湍流模型对低速可压管流问题有良好的表现。 　　1.统一控制方程形式 　　式中， 为通用变量，可以代表 ， ， ， 等求解变量； 为广义扩散系数； 为广义源项 ，式中各项依次为瞬态项（transient term）、对流项（convective term）、扩散项（diffusive term）和源项（source term）。 　　2.标准 方程 　　标准 模型需要求解湍动能 和耗散率 方程。该模型假设流体流动为完全湍流，忽略分子的粘性影响。该模型的湍动能 和耗散率 方程为： 　　-由于平均速度梯度引起的湍动能； 　　-由于浮力影响引起的湍动能； 　　-可压缩湍流脉动膨胀对总的耗散率的影响； 　　-湍流粘性系数， 。 　　此模型中的各个常数在FLUENT中的默认值为 ， 　　，湍动能 和 耗散率 的湍流普朗特数分别为： ， 　　该模型适合对完全湍流流动的模拟。 　　二、气流压力脉动实验及管道系统计算模型 　　1.孔板消减气流脉动实验 　　为了对管道系统气流脉动的模拟结果进行验证，依照图1搭建了往复式压缩机管道系统压力脉动测试实验平台。测试系统主要由一台往复式压缩机、变频器、与其相连的管道系统、直流稳压电源、压力传感器、示波器、数据采集及分析系统组成。 　　2.管道系统的计算模型 　　根据搭建的管道系统实验平台，在不影响模拟结果的前提下，为提高计算效率采用如图2所示（尺寸单位为mm）的二维计算模型。将模型导入Gambit中划分单元网格，网格类型采用非结构化（Pave） 　　的四边形网格（Quad），共得到17644个单元。最后在Gambit中对管道系统添加边界条件，入口边界条件采用非定常压力入口边界条件，出口边界条件为压力出口边界。 　　三、气流脉动模拟及实验分析 　　1.模拟分析 　　将在Gambit中化分好的网格模型导入Fluent中，在Fluent中根据管道内部介质的性质采用基于密度（耦合式）的隐式2ddp（二维双精度）求解器，并指定其计算模式为非稳态；计算模型选择标准 　　双方程湍流模型；流体类材料为理想空气。根据现场实验条件，出口边界条件采用非定常压力边界条件： （根据其公式用C语言编译非定常速度UDF函数）即施加一个频率为7.3Hz（模拟转速为438r·min-1压缩机）、压力不均匀度为24%的脉动压力条件，入口边界条件采用定常压力边界条