无油双涡圈空气涡旋压缩机的数值模拟及试验研究
Numerical Simulation and Experimental Study of Oil-Free Double-Warp Air Scroll Compressor
通讯作者: 彭斌,男,教授,博士生导师,电话(Tel): 0931-5172154; E-mail:pengb2000@163.com.
责任编辑: 石易文
收稿日期: 2021-04-6
基金资助: |
|
Received: 2021-04-6
作者简介 About authors
孙健(1991-),男,甘肃省张掖市人,博士生,主要研究方向为涡旋机械及低温余热发电.
为了分析无油双涡旋齿压缩机工作腔内部流场的流动特性,基于计算流体动力学方法对压缩机工作腔内部流场进行了三维非稳态数值模拟;为了验证数值模拟与实际工作过程的偏差,搭建了以空气为工质的涡旋压缩机试验测试平台,得到了无油双涡圈压缩机的实际输出性能参数.研究结果表明:在相邻工作腔之间的质量交换,对工作腔内以及沿着齿高方向的温度和速度的分布影响较大,对压力的影响较小;进出口质量流量和流速呈周期性变化,压缩机转速越高,理论容积流量越大.
关键词:
In order to analyze the flow characteristics of the internal flow field in the working chamber of the oil-free double-warp scroll compressor, a three-dimensional unsteady numerical simulation of the internal flow field in the working chamber of the compressor is conducted based on the computational fluid dynamics (CFD) method. In order to verify the deviation between the numerical simulation and the actual working process, a test platform for scroll compressors using air as the working medium is built, and the actual output performance parameters of the oil-free double-warp scroll compressor are obtained. The results show that the mass exchange between adjacent working chambers has a greater impact on the working chamber and the distribution of temperature and velocity along the tooth height direction, but has a little impact on pressure. The inlet and outlet mass flow and velocity rate change periodically, and as the compressor speed increases, the theoretical volume flow increases.
Keywords:
本文引用格式
孙健, 彭斌, 朱兵国.
SUN Jian, PENG Bin, ZHU Bingguo.
从已有的文献可以看出,目前对于双涡旋齿压缩机的研究大多只是集中在对其型线和几何模型的研究.文献[3-4]发现相较于单涡圈压缩机,双涡圈压缩机工作腔内部的传热特性更加复杂,并提出了多齿头型线啮合的基本准则. 文献[5]建立了双涡旋齿压缩机的热力学模型,并定量分析了基圆半径等基本几何参数对压缩机性能的影响. 文献[6]基于计算流体动力学(CFD)方法,建立了涡旋压缩机二维非稳态流动和传热的热力学模型,分析发现在工作腔内除压力外其他的性能参数都有着较强的空间分布特性. 文献[7⇓-9]基于CFD方法发现由于动涡旋齿对进排气口的影响,气体在进出口均存在气流脉动现象. 文献[10]发现在涡旋压缩机的两个中心腔之间存在着压力差,由于动涡旋齿会周期性的遮蔽排气口,所以出口质量流量存在着明显的波动趋势. 文献[11-12]提出了一种新的结构化动网格(SDM)生成方法,通过SDM可以生成高质量的流体域网格,并通过等距曲线法建立了3种不同齿头修正下涡旋压缩机齿头完全啮合的几何模型,并通过CFD方法数值模拟分析了这3种几何模型在相同吸气流量下排气压力的变化规律. 文献[13]对涡旋压缩机进行了二维非定常流动的数值模拟,并建立了一个新的包括Nusselt数泄漏模型的热力学模型. 文献[14]通过对涡旋压缩机的三维非稳态数值模拟,发现由于排气口的设置,在中心腔的两个腔室之间存在着压力差. 文献[15]采用基于粒子群算法的程序结合CFD方法对涡旋压缩机的基本几何参数进行了优化和分析,发现工作腔的形状、尺寸和理论压缩比都对压缩机的性能有着一定的影响. 文献[16]利用中线法设计了一种新的不对称双涡圈涡旋压缩机,并对其进行了三维非稳态数值模拟,发现优化后压缩机能够有效的提高排气量和内容积比. 然而在已有的研究中大都是通过单一的数学模型或数值模拟方法,对双涡圈涡旋压缩机进行理论研究和分析,并未将数学模型、数值模拟和试验研究相结合来综合评价涡旋压缩机的输出性能.
本文以一种无油双涡旋齿涡旋压缩机为研究对象,分析和比较了单双涡旋齿压缩机工作腔容积以及容积效率的变化规律;基于CFD方法对压缩机工作腔内部流场进行了三维非稳态数值模拟,分析和研究了变工况下工作腔内流体温度、压力和速度的变化规律;通过所搭建的涡旋压缩机试验平台验证了数值模拟的准确性,对无油双涡圈压缩机的改进和优化具有一定的指导意义.
1 数学模型
1.1 几何模型
研究和试验所用样机的涡旋齿型线由圆渐开线构成,其基本几何参数如表1所示.其中:rb为基圆半径;α为渐开线发生角;t为涡旋齿齿厚;h为涡旋齿齿高;φb为渐开线起始角;φe为渐开线结束角.
表1 涡旋齿基本几何参数
Tab.1
参数 | 取值 | 参数 | 取值 |
---|---|---|---|
rb/ m | 9.231×10-3 | α/(°) | 13.98 |
t/m | 4.450×10-3 | φb/(°) | 45 |
h/m | 65.00×10-3 | φe/(°) | 765 |
1.1.1 双涡旋齿压缩机的运行原理
图1
图1
双涡旋齿压缩机运行原理图
Fig.1
Operation principle diagram of double-warp scroll compressor
1.1.2 双涡旋齿几何结构
试验样机和研究样机的几何结构示意图如图2所示.其中:pt为涡旋齿节距.
图2
双涡旋齿压缩机的动涡旋齿可以看做是由一条单涡旋齿型线旋转180°所形成的,将动涡旋齿旋转90°可以形成静涡旋齿;设OP1和OP2分别为第1和第2条动涡旋齿,则双涡旋齿压缩机的动涡旋齿型线分别如下式所示[5]:
式中:下标d1和d3为涡旋齿内侧型线,d2和d4为涡旋齿外侧型线.
1.2 工作腔容积
无油双涡旋齿压缩机的工作腔容积计算公式如下[5].
压缩腔容积:
吸气腔容积:
rb[2(1-cosθ)-
2(φe-π)sinθ- πsin2θ]
排气腔容积:
式中: i为压缩机第i个压缩腔.
1.3 容积流量
涡旋压缩机的转速和基本几何参数对压缩机的理论容积流量有着极为重要的影响,转速越高理论容积流量也就越大,理论容积流量计算公式如下[1]:
式中: n为压缩机转速.
涡旋压缩机的容积效率ηv不但反映了容积的利用效率,而且还可以用来衡量压缩机的内泄漏程度,其计算公式如下:
式中: qv为实测容积流量.
1.4 讨论与分析
为了分析单双涡旋齿压缩机的性能,比较了在相同基本几何参数下,压缩机的基本输出性能.压缩机的基本几何参数如表2所示,其中: ror为回转半径.
表2 单双涡旋齿压缩机基本几何参数
Tab.2
参数 | 单涡旋齿 | 双涡旋齿 |
---|---|---|
rb/m | 9.231×10-3 | 9.231×10-3 |
α/(°) | 13.980 | 13.980 |
t/m | 4.45×10-3 | 4.45×10-3 |
ror/m | 24.500×10-3 | 10.000×10-3 |
1.4.1 压缩机工作腔容积
单双涡旋齿压缩机工作腔容积随主轴转角的变化规律如图3所示.其中: V为工作腔容积;ΔV为单双涡旋齿压缩机容积之差.在基本几何参数一致时,压缩机工作腔容积的变化趋势基本一致,单双涡旋齿压缩机的最大吸气容积之差为4.3×10-6 m3.
图3
1.4.2 容积效率
单双涡旋齿压缩机容积效率随转速的变化规律如图4所示,其中:Δηv为单双涡旋齿压缩机容积效率之差.在转速低于 2000 r/min时,双涡旋齿压缩机的容积效率一直高于单涡旋齿压缩机,随着转速的继续升高,单双涡旋齿压缩机的容积效率都不再变化.压缩机在低转速下工作时,由于流体工质在工作腔内停留的时间相对较长,使得相邻工作腔之间的质量交换量增加,所以降低了压缩机的容积效率.
图4
2 数值模拟
2.1 控制方程
在涡旋压缩机的实际工作过程中,工作腔内流体运动较为复杂,因此在数值模拟时应该符合质量、动量和能量守恒方程,以上方程表示如下所示[17].
质量守恒方程:
动量守恒方程:
能量守恒方程:
+ +
+Sr
式中:
2.2 网格划分
图5
表3 各计算域网格数量
Tab.3
名称 | 体网格单元数 | 网格节点数 |
---|---|---|
进气管 | 23,337 | 10,189 |
工作腔 | 5,613,089 | 2,608,404 |
排气管 | 80,100 | 13,762 |
总计 | 5,716,526 | 2,632,355 |
2.3 计算方法及边界条件
2.3.1 计算方法
在进行涡旋压缩机三维非稳态流场数值模拟时,采取动网格中自定义宏命令(UDF)驱动的方式,选择DEFINE_CG_MOTION来描述动涡旋齿的运动,动涡盘的速度控制方程如下所示:
式中: ω为动涡旋齿旋转角速度.
2.3.2 计算模型和边界条件
涡旋压缩机在实际工作过程中处于高速状态和高温状态,选择重正化群(RNG)理论κ-ε湍流模型来描述工作腔内的运动,工质传热采用高阶二阶迎风模式,选择Piso算法进行仿真计算,结合试验结果设置如下所示的初始条件: 初始进气压力ps=0.1 MPa,初始进气温度Ts=300 K,转速n=3000 r/min,流体工质为理想空气.
3 分析与结果
以涡旋压缩机进出口质量流量的变化趋势来衡量数值模拟是否达到了稳定状态,在计算了2个循环后,工作腔内流体的运动已经达到了稳定状态,因此以第3个循环来分析工作腔内温度、压力和速度变化云图,以第4和第5个循环来分析和研究进出口质量流量和速度的变化规律.由于试验样机在进行试验时,测得其排气压力pd=0.8 MPa,所以为了研究和分析变工况下双涡旋齿压缩机的进出口质量流量和流速的变化规律,设置如下工况条件: ① 相同转速n=3000 r/min,不同排气压力pd=0.7, 0.8,0.9 MPa; ② 相同排气压力pd=0.8 MPa,不同转速n=2000,3000,5000 r/min.
3.1 温度云图
双涡旋齿压缩机工作腔内工质流体温度云图的变化规律如图6所示,其中:T为工作腔内温度.随着涡旋压缩机的不断转动,工质流体被不断地压缩,工作腔内的温度会逐渐升高;由于内泄漏的存在,在相邻两个工作腔之间存在着质量的交换,会使得温度在工作腔内分布不均匀,相邻工作腔之间最大温差为127 K;沿着齿高方向,温度场分布也并不均匀在, 涡旋齿上最大温差为161 K.
图6
3.2 压力云图
涡旋压缩机工作腔内工质流体压力云图的变化规律如图7所示,其中:p为工作腔内压力.可知,越靠近中心腔压力越高,泄漏对工作腔内工质压力几乎没有太大的影响;在相邻工作腔之间,最大压差近0.44 MPa;沿着齿高的方向,压力分布均匀.
图7
3.3 速度云图
图8
3.4 进出口质量流量
双涡旋齿压缩机进出口质量流量随主轴转角的变化规律如图9所示,其中:qm,in和qm,out分别为进出口质量流量;qm为质量流量;入口负值表示压缩机吸气,出口正值表示压缩机排气.进口最大质量流量出现在270°附近,出口最大质量流量出现在180°附近.图9(a)所示为在额定转速下,涡旋压缩机进出口质量流量在两个工作周期内随主轴转角的变化规律;由于在相邻工作腔之间存在着质量交换,所以出口质量流量会大于进口质量流量;由于动涡旋齿会周期性的遮蔽排气口,所以出口质量流量的波动变化较大.图9(b)所示为在额定转速不同排气压力时,进出口质量流量随主轴转角的变化规律;当排气压力过大时,压缩机需要克服气体力逆向做功,因此排气压力越大进出口处质量流量越小.图9(c)所示为涡旋压缩机在相同排气压力不同转速时,进出口质量流量随主轴转角的变化规律;转速越高,进出口质量流量越大;当转速较低时,由于工质流体在压缩机工作腔内停留的时间较长,所以会增加工作腔之间的内泄漏量.
图9
图9
双涡旋齿压缩机进出口质量流量
Fig.9
Inlet and outlet mass flow of double-warp scroll compressor
3.5 进出口流速
无油双涡旋齿压缩机进出口流速随主轴转角的变化规律如图10所示,其中: vin和vout分别为压缩机进出口流速,进口最大流速出现在270°附近,出口最大流速出现在180°附近.图10(a)所示为涡旋压缩机进出口流速在两个工作周期内的变化规律,由于单涡旋齿会周期性的遮蔽排气口,所以排气口速度变化波动较大.图10(b)所示为涡旋压缩机在额定转速不同排气压力时进出口流速随主轴转角在一个周期内的变化规律,如图所示排气压力越低进出口流速越高;过大的排气压力,会使得压缩机出口出现回流现象,压缩机在较高的排气压力下工作时,进出口流速都会有所降低.图10(c)所示为涡旋压缩机在相同排气压力不同转速时,进出口流速随主轴转角的变化规律.由图10可知,转速越高进出口流速越高.
图10
图10
双涡旋齿压缩机进出口流速
Fig.10
Inlet and outlet velocity of double-warp scroll compressor
4 试验研究
为了验证数值模拟的准确性,搭建了以空气为工质的双涡旋齿压缩机试验测试平台.试验测试系统如图11所示.整个试验测试系统包括: 涡旋压缩机组、冷却装置和数据采集系统等.试验样机的额定转速为 3000 r/min,进气压力为0.1 MPa,环境温度为22 ℃.
图11
图11
试验测试图
1—电机冷却管道; 2—驱动电机; 3—涡旋压缩机; 4—压缩机冷却管道; 5—进气管道; 6—排气管; 7—单向截止阀; 8—温度传感器
Fig.11
Test chart
4.1 排气流量
图12
4.2 排气温度和噪声
无油双涡旋齿压缩机排气温度和机体噪声随转速的变化规律如图13所示,其中:Tout为压缩机排气温度;Nnoi为压缩机产生的噪声.所测排气温度为经过液体强制冷却之后的温度,因此排气温差较小只有12 ℃;在转速低于 1800 r/min时,由于液体强制对流换热的影响,排气温度基本没有太大的变化;当转速高于 1800 r/min时,由于冷却系统液体箱内液体温度会逐渐升高,所以排气温度也会开始逐步缓慢的升高;当转速达到 3200 r/min时,最高排气温度为152 ℃.随着压缩机转速的逐渐升高,机体产生的噪声值也会随之而增大,过高的噪声会对周围环境产生噪声污染,因此在优化时应该采取降噪措施.
图13
图13
双涡旋齿压缩机的排气温度和机体产生的噪声
Fig.13
Discharge temperature and body noise of double-warp scroll compressor
4.3 驱动电机电流和功率
无油双涡旋齿压缩机驱动电机功率和电流随转速的变化规律如图14所示,其中:I为驱动电机电流;P为驱动电机功率.随着压缩机转速的不断增加,驱动电机功率和电流也随之而增大;在额定转速3000 r/min时,驱动电机电流和功率分别为8.1 A和13.1 kW.
图14
5 结论
(1) 双涡旋齿压缩机是一种低压比、大排量的涡旋式机械,其容积效率高于单涡旋齿压缩机.
(2) 工作腔之间的内泄漏会扰乱相邻工作腔之间的流场分布;双涡旋齿压缩机相邻工作腔之间的传热更加复杂,内泄漏会使得出口质量大于进口质量;动涡旋齿对排气口的影响会使得出口质量和流速的波动较大;转速与流量和流速呈正相关,排气压力越高,流量和流速越低.
(3) 在额定转速下,无油双涡旋齿压缩机排气温度为148 ℃,机体产生的噪声为83.9 dB,容积流量的理论值和试验值的最大差值为0.158 m3/min,偏差在4%以内,从而验证了所构建的热力学模型的精确性;通过试验测试发现,液体冷却方法对大排气量的无油涡旋压缩机有着极佳的冷却效果.
参考文献
新型无油涡旋压缩机性能
[J]. ,
Performances of new oil-free scroll compressors
[J]. ,
双涡圈及多涡圈涡旋机械几何特性的理论研究
[J]. ,
Improved characteristics of scroll machines with multiple spiral
[J]. ,
多齿涡旋机械啮合型线的构建理论
[J]. ,
Construction theories of meshing profile for multi-wraps scroll machinery
[J]. ,
涡旋齿数对多齿涡旋压缩机性能的影响
[J]. ,
Influence of wrap number on performance of multi-wrap scroll compressor
[J]. ,
型线参数对双涡圈涡旋压缩机性能的影响
[J]. ,
Effects of various profile parameters on the performance of a double-scroll compressor
[J]. ,
Convective heat transfer in a scroll compressor chamber: A 2-D simulation
[J]. ,DOI:10.1016/j.ijthermalsci.2003.11.005 URL [本文引用: 1]
Effects of suction port arrangements on a scroll expander for a small scale ORC system based on CFD approach
[J]. ,DOI:10.1016/j.apenergy.2015.04.046 URL [本文引用: 1]
Unsteady flow in the suction process of a scroll expander for an ORC waste heat recovery system
[J]. ,DOI:10.1016/j.applthermaleng.2015.01.010 URL [本文引用: 1]
The impact of a bilateral symmetric discharge structure on the performance of a scroll expander for ORC power generation system
[J]. ,DOI:10.1016/j.energy.2018.06.053 URL [本文引用: 1]
Analysis of the three-dimensional transient flow in a scroll refrigeration compressor
[J]. ,DOI:10.1016/j.applthermaleng.2017.08.115 URL [本文引用: 1]
Novel structured dynamic mesh generation for CFD analysis of scroll compressors
[J]. ,DOI:10.1177/0957650915601927 URL [本文引用: 1]
Design methodology and geometric modeling of complete meshing profiles for scroll compressors
[J]. ,DOI:10.1016/j.ijrefrig.2018.05.011 URL [本文引用: 1]
Internal flow field and heat transfer investigation inside the working chamber of a scroll compressor
[J]. ,DOI:10.1016/j.energy.2020.117700 URL [本文引用: 1]
Unsteady characteristic and flow mechanism of a scroll compressor with novel discharge port for electric vehicle air conditioning
[J]. ,DOI:10.1016/j.ijrefrig.2020.06.022 URL [本文引用: 1]
CFD-based optimization of scroll compressor design and uncertainty quantification of the performance under geometrical variations
[J]. ,DOI:10.1016/j.energy.2020.118382 URL [本文引用: 1]
涡旋压缩机不对称双涡旋齿型线的构建与模拟
[J]. ,
Construction and simulation of an asymmetric twin-wrap profile for scroll compressors
[J]. ,
Scroll compressor with an arcuate and a logarithmic spiral sections: US 9605675 B2
[P].
/
〈 | 〉 |