1. 引言
往复式压缩机是一种常见的气体增压设备，广泛应用于动力、化工、制冷等领域。而气阀作为压缩机的重要组件之一，起着控制气体进出压缩机的作用，对气体压缩程度和压缩机性能起着至关重要的作用。然而，在实际应用中，压缩机还存在阀片启闭不及时、颤振和撞击等问题。上述问题的存在促使人们对压缩机气阀展开研究。
2. CFD数值模拟进展
自上世纪50年代开始，往复式压缩机的热力计算都采用气阀理论进行相关结构的设计计算。该理论以理想气体的一元绝热流动假设为前提，建立描述阀片运动和汽缸压力比的方程组。实践表明，该理论适用于气阀结构的初步设计，但也存在许多局限性和不足，主要表现为：1) 一维假设忽略了汽缸内压力、温度的差异性，无法得到温度随曲柄转角的变化；2) 由于忽略了阀隙的流动细节，该理论对于改进气阀结构的指导作用不大。
随着计算流体动力学(CFD)学科在上世纪80年代后的快速发展，应用CFD方法模拟热力设备的工作过程成为可能。分析热力设备可以发现，他们具有共同的显著特点，即气体温度场、流场和运动部件的运动相互耦合。对于往复式压缩机的气阀，各物理场互相影响，紧密联系，其分布特点反映了气阀结构与运行参数的匹配程度，分析多场的耦合特点成为深入分析其热力过程的关键。在1994年，Cyklis [1] 对往复式压缩机阀隙的二维稳态流场进行了CFD数值模拟，针对阀片的不同开度，建立了不同的二维轴对称计算模型。在此基础上，验证CFD数值模拟在压缩机气阀中的可行性及提高CFD数值计算的精度是研究的重点，主要通过引进实际气体的状态方程和高级紊流模型来实现。
A. P. Peskin [2] 于1999年研究了不同气体状态方程对往复式压缩机数值模拟结果的影响，引进了四种气体状态方程，亥姆霍兹能量状态方程式最复杂，计算的密度、焓值最准确，在压力不太高时采用第二维里系数的维里方程产生的计算误差仅百分之几，但压力较高时，引起的误差与理想气体方程相当，P-R方程计算误差更大，密度值高达6%，焓值在5%和15%之间，理想气体状态方程与亥姆霍兹方程相比，各计算数据中排气温度的差别最大。自1999年开始，针对阀隙的低雷诺数紊流流场进行了相关研究，Pérez-Segarra等人[3] 分别应用RANS k-C和k-w双方程紊流模型进行了数值模拟，Rovaris等人[4] 对雷诺数为2.5 × 104的紊流进行了LES大涡模拟，而Colaciti等人[5] 应用RANS紊流模型模拟了马赫数接近0.2的高雷诺数(1.6 × 106)流场。随着计算机的进一步发展，J. Rigola等人[6] 在2012年开始应用LES大涡模型对低雷诺数流动(3 × 105)进行数值模拟，分析了气体压缩性对模拟结果的影响，认为入口马赫数为0.1时，流体的压缩性可以忽略不计，但马赫数达到0.3时，必须考虑压缩性的影响。
提高往复式压缩机CFD模拟准确性的另外一项重要工作是改进模拟方法，即从二维发展到三维、建立复杂流道模型并改进网格质量，应用动网格技术。在2005年，Roland Aigner [7] 用数值模拟方法研究了阀片的动态性能和压缩机内的压力脉动，比较了二维、三维和根据压缩机厂家经验开发的准一维模型的计算结果，认为应用商业化软件FLUENT 6.2进行数值模拟得到的二维计算结果与实测值最接近，但计算精度与网格大小、离散式精度紧密相关，三维模拟时由于网格太多，计算两周也未收敛。准一维模型考虑了气缸内的压力脉动，描绘的气体流动特点、压力分布规律和阀片运动基本可满足工程精度要求。Yogesh V. Birari等人[8] 利用CFD方法设计和改进了一台往复式制冷压缩机，确定了由定子、转子、排气管、吸液管、吸气消声器和上盖组成的计算区域，在忽略热辐射，假设制冷剂不可压缩、稳定流动的条件下，进行结构化、块状结构化、非结构化网格划分，然后设置边界条件，求解，获得热力场，显示两种工质的温度分布，与实验数据相比，最大误差仅5%；计算得到进排气通道压降，与实验比较吻合较好，优化了进排气通道结构；用动网格技术模拟了压缩过程，计算区域包括余隙容积在内的整个气缸容积，得到活塞运动参数。CFD技术的应用，大大缩短了实验时间，与实验结果的比较，验证了模拟的正确性。Akira Nakano等人[9] 研究了CFD方法在改变往复式压缩机吸气消音器压力、提高制冷量和制冷效率中的应用，模拟过程分为两个阶段，首先采用转移矩阵法计算出存在进气阀片振动时的缸内压力变化，然后将工质视为三元非稳定可压缩粘性流体，以前面计算的缸内压力为边界条件，应用FLUENT分析吸气消音器中的压力变化，从而优化制冷剂的质量流量，降低吸气消音器中的热量损失和压力损失。2010年，Raimund Almbauer等人[10] 运用计算流体动力学方法对往复式压缩机的热力过程进行二维数值模拟，通过比较不同结构的气缸压力分布，改进了压缩机结构，认为CFD和有限元方法(FEM)是进行热力设备优化设计的有效快速工。同一年，Kenji Kinjo等人[11] 运用CFD方法对一台直径20.6 mm、气缸行程18 mm、曲轴转速3000 rpm、压缩比为10的往复式压缩机进行了数值计算，建立笛卡尔网格模型，运用虚流量法，模拟了往复式压缩机的工作过程，获得活塞运动引起的压缩腔内气体流动规律以及由气流推力引起的簧片阀片开启和关闭规律，比较计算结果与实验数据，验证了流–固耦合模拟的有效性。2012年，Guido Pratelli [12] 等人对往复式压缩机的整体模型进行了稳态的CFD模拟数值，该模型包含16个气阀、气缸头，气阀减荷器，计算的缸内P-V分布与测量数据吻合较好，Hiteshkumar MISTRY等人[13] 建立了簧片阀的三元瞬态计算模型，弹性簧片分别按悬臂梁和圆盘进行处理，运用CFD对含有压力脉冲装置的压缩机工作过程进行数值模拟，与实测的压力曲线相比，按悬臂梁计算的结果具有更好的吻合性。
21世纪以来，数值模拟得到广泛应用。辽宁石油化工大学赵斌采用CFD方法对压缩机气阀的内流场进行仿真计算[14] ，分析了气阀内的速度场和压力场分布特点。广西大学卢朝霞等人运用有限元理论分析了往复式压缩机环状阀阀片与升程限制器及阀座正碰撞与倾侧碰撞产生的冲击应力[15] ，北京化工大学诊断与自愈过程研究中心武炳等人运用CFD方法对往复式压缩机进气阀的受力进行了数值计算[16] ，通过建立不同开度的进气阀和气缸模型进行流场模拟，分析了阀片受气体力与阀片开启高度的关系。
3. 结语
当前，往复式压缩机的相关研究主要集中在气阀阀片运动规律的理论计算、气阀弹簧参数的合理设计方面，对于压缩机流场、温度场、阀片运动的耦合数值计算研究不多。虽然在压缩机的CFD数值计算方面有些研究，但大多限于二维计算，对计算模型也作了较多简化，不能准确地模拟压缩机的实际工作过程。压缩机的工作过程具有显著的流–热–固多场耦合特点，进行多场耦合计算是目前研究的重点问题之一。CFD广泛应用于压缩机的数值模拟，成为提高压缩机设计效率的有效工具，随者CFD模拟技术的发展，对压缩机内流场进行三维动态模拟是当前和今后的发展趋势，将为深入分析气阀工作规律提供很好的基础。