[1]王炳棋,畅云峰.压缩机工作过程模拟中数值振荡问题的成因分析及改善方法[J].西安交通大学学报,2020,54(06):161-167.[doi:10.7652/xjtuxb202006021]
 WANG Bingqi,CHANG Yunfeng.Analysis and Improvement of Numerical Oscillation in Simulation of Compressor Working Process[J].Journal of Xi'an Jiaotong University,2020,54(06):161-167.[doi:10.7652/xjtuxb202006021]
点击复制

压缩机工作过程模拟中数值振荡问题的成因分析及改善方法
分享到:

《西安交通大学学报》[ISSN:0253-987X/CN:61-1069/T]

卷:
54
期数:
2020年第06期
页码:
161-167
栏目:
出版日期:
2020-06-10

文章信息/Info

Title:
Analysis and Improvement of Numerical Oscillation in Simulation of Compressor Working Process
文章编号:
0253-987X(2020)06-0161-07
作者:
王炳棋 畅云峰
西安交通大学能源与动力工程学院, 710049, 西安
Author(s):
WANG Bingqi CHANG Yunfeng
School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China
关键词:
压缩机 数值刚性 后向差分算法
Keywords:
compressor stiffness backward difference formula
分类号:
TB652
DOI:
10.7652/xjtuxb202006021
文献标志码:
A
摘要:
为解决数值模拟压缩机工作过程中存在的数值振荡和发散问题,采用四阶后向差分(BDF-4)算法和经典四阶龙格库塔(RK-4)算法分析压缩机基元容积控制方程的稳定性,并对2种算法的计算结果进行对比。结果表明:在吸排气过程中,压力计算的微小偏差会导致质量流量产生较大的误差,从而使基元容积控制方程产生较大的数值刚性,限制了RK-4算法的最大稳定步长; 在稳定性上,BDF-4算法的绝对稳定域包含整个负半x轴,可以避免数值刚性引起的振荡和发散,其计算步长不再受数值刚性的限制; 对计算效率而言,每步BDF-4算法所消耗时间为同步长RK-4算法的1.7~1.8倍; 计算精度上,BDF-4算法与RK-4算法代数精度相同,实际计算时两者的表现也十分接近。BDF-4算法在保证工程精度的前提下提高了压缩机工作过程模拟的稳定性,解决了数值振荡问题,对容积式压缩机的工作过程模拟、新型泄漏通道的研究等有实际意义。
Abstract:
To solve the problems of numerical oscillations and divergence in the simulation of compressor process, 4th-order backward difference formula(BDF-4)and classical 4th-order Runge-Kutta(RK-4)method are adopted to analysis the stability of control equations of element volume, and the results of two methods are compared. The analysis results show that a small fluctuation of pressure in the process of suction and exhaust would induce a large fluctuation in the calculation of the mass flow, which would result in a large stiffness of the governing equation and is the key factor to cause numerical oscillations. The results indicate that there is no stability constraint in the calculation using the BDF-4 method, so that the oscillations and divergence due to the numerical stiffness can be avoided. However, the calculation time of the BDF-4 method is 1.7~1.8 times of the RK-4 method at one step, and the simulation accuracy of the BDF-4 method is nearly the same as that of the RK-4 method. The BDF-4 method can greatly improve the simulation stability of the compressor working process without loss of accuracy, and give a reference for simulation of compressor thermal process and the leakage simulation through clearance.

参考文献/References:

[1] 邢子文, 吴华根, 束鹏程. 螺杆压缩机设计理论与关键技术的研究和开发 [J]. 西安交通大学学报, 2007, 47(7): 755-763.
XING Ziwen, WU Huagen, SHU Pengcheng. Research and development on design theory and key technology of screw compressors [J]. Journal of Xi’an Jiaotong University, 2007, 47(7): 755-763.
[2] 吴华根, 邢子文, 束鹏程. R134a螺杆制冷压缩机工作过程数值模拟及实验研究 [J]. 西安交通大学学报, 2004, 38(5): 487-490.
WU Huagen, XING Ziwen, SHU Pengcheng. Theoretical and experimental investigation on working process of screw refrigeration compressor with R134a [J]. Journal of Xi’an Jiaotong University, 2004, 38(5): 487-490.
[3] 邢子文, 邓定国, 束鹏程. 双螺杆压缩机工作过程的计算机模拟研究 [J]. 流体工程, 1991(12): 11-16.
XING Ziwen, DENG Dingguo, SHU Pengcheng. Research on computer simulation of twin screw compressor working process [J]. Fluid Machinery, 1991(12): 11-16.
[4] 卢阳, 郭蓓, 周瑞鑫. 双螺杆真空泵工作过程数值模拟 [J]. 西安交通大学学报, 2015, 49(7): 67-71.
LU Yang, GUO Bei, ZHOU Ruixin. Numerical simulation for working process in twin screw vacuum pump [J]. Journal of Xi’an Jiaotong University, 2015, 49(7): 67-71.
[5] WU Huagen, XING Ziwen, SHU Pengcheng. Theoretical and experimental study on indicator diagram of twin screw refrigeration compressor [J]. International Journal of Refrigeration, 2004, 27(4): 331-338.
[6] 邢子文. 螺杆压缩机: 理论、设计及应用 [M]. 北京: 机械工业出版社, 2000: 117-132.
[7] 齐元渠, 余岳峰. 双螺杆膨胀机泄漏过程的数值模拟与分析 [J]. 上海交通大学学报, 2016, 50(4): 496-501.
QI Yuanqu, YU Yuefeng. Numerical simulation and analysis of leakage process in twin screw expander [J]. Journal of Shanghai Jiao Tong University, 2016, 50(4): 496-501.
[8] SALIKEEV S, BURMISTROV A, BRONSHTEIN M, et al. Conductance calculation of channels in laminar gas flow regime at an arbitrary pressure difference [J]. Vacuum, 2014, 107: 178-183.
[9] SALIKEEV S, BURMISTROV A, BRONSHTEIN M, et al. Conductance calculation of slot channels in laminar gas flow regime at small pressure difference [J]. Vacuum, 2014, 99: 303-306.
[10] BURMISTROV A, SALIKEEV S, BRONSHTEIN M, et al. Conductance calculation of slot channels with variable cross section [J]. Vakuum in Forschung und Praxis, 2015, 27(1): 36-40.
[11] LAMBERT J D. Numerical methods for ordinary differential systems: the initial value problem [M]. New York, USA: John Wiley & Sons Ltd., 1991: 98-103, 149-157, 213-254.
[12] HAIRER E, N?RSETT S P, WANNER G. Solving ordinary differential equations: I Nonstiff problems [J]. SIAM Review, 1990, 32(3): 485.
[13] ENRIGHT W H, HULL T E, LINDBERG B. Comparing numerical methods for stiff systems of O.D.E:s [J]. BIT Numerical Mathematics, 1975, 15(1): 10-48.
[14] ASHINO R, NAGASE M, VAILLANCOURT R. Behind and beyond the MATLAB ODE suite [J]. Computers & Mathematics with Applications, 2000, 40(4/5): 491-512.
[15] SHAMPINE L F. Evaluation of a test set for stiff ODE solvers [J]. ACM Transactions on Mathematical Software, 1981, 7(4): 409-420.
[16] HAIRER E, WANNER G. On the instability of the BDF formulas [J]. SIAM Journal on Numerical Analysis, 1983, 20(6): 1206-1209.
[17] BUI T D, BUI T R. Numerical methods for extremely stiff systems of ordinary differential equations [J]. Applied Mathematical Modelling, 1979, 3(5): 355-358.
[18] SHAMPINE L F. Measuring stiffness [J]. Applied Numerical Mathematics, 1985, 1(2): 107-119.
[19] SHAMPINE L F, REICHELT M W. The MATLAB ODE suite [J]. SIAM Journal on Scientific Computing, 1997, 18(1): 1-22.
[20] HAIRER E, WANNER G. Solving ordinary differential equations II [M]. Berlin, Germany: Springer-Verlag, 1996: 16-18, 240-247.
[本刊相关文献链接]
李元,全恺,吴睿,等.往复压缩机管路大气流脉动的数值模拟和实验研究.2019,53(11):20-26.[doi:10.7652/xjtuxb 201911003]
赵家毅,王志恒,席光,等.大预旋下离心压缩机进口结构与导叶掠角的改进研究.2017,51(11):1-6.[doi:10.7652/xjtuxb201711001]
吴小华,杨俊玲,越云凯,等.一种新型热压缩机运行特性实验研究.2017,51(11):144-149.[doi:10.7652/xjtuxb201711 020]
王晟旻,琚亚平,张楚华.离心压缩机级通流与盘盖侧泄漏流的耦合分析.2017,51(7):73-77.[doi:10.7652/xjtuxb201707 011]
史云昊,陈雪飞,宁远钊,等.采用关联模型的离心压缩机喘振工况数值研究.2017,51(5):128-133.[doi:10.7652/xjtuxb201705018]
武晓昆,陈文卿,周明龙,等.双螺杆制冷压缩机气流脉动衰减器的研究与开发.2017,51(4):23-29.[doi:10.7652/xjtuxb201704004]
王枫,陈征,李花,等.采用多体动力学的压缩机曲轴结构优化.2017,51(3):7-13.[doi:10.7652/xjtuxb201703002]
王增丽,申迎峰,王振波,等.单螺杆压缩机工作腔喷液过程几何特性及工作机理研究.2017,51(2):128-134.[doi:10.7652/xjtuxb201702020]
慕光宇,王枫,米小珍.压缩机舌簧排气阀运动模型对比及动态特性影响因素分析.2017,51(9):69-76.[doi:10.7652/xjtuxb201709010]
王亚洲,秦国良,和文强,等.时空耦合谱元方法求解一维Burgers方程.2017,51(1):45-50.[doi:10.7652/xjtuxb2017 01008]

相似文献/References:

[1]张金权,畅云峰.动磁直线电机驱动压缩机的建模与仿真[J].西安交通大学学报,2007,41(07):820.[doi:10.7652/xjtuxb200707015]
 Zhang Jinquan,Chang Yunfeng.Modeling and Simulation of Compressor Driven by MovingMagnet Linear Motor[J].Journal of Xi'an Jiaotong University,2007,41(06):820.[doi:10.7652/xjtuxb200707015]
[2]冯诗愚,高秀峰,刘卫华,等.内啮合转子压缩机齿间啮合效率研究[J].西安交通大学学报,2007,41(09):1040.[doi:10.7652/xjtuxb200709008]
 Feng Shiyu,Gao Xiufeng,Liu Weihua,et al.Investigation to Meshing Efficiency of InternalMeshing Rotary Compressor[J].Journal of Xi'an Jiaotong University,2007,41(06):1040.[doi:10.7652/xjtuxb200709008]
[3]费继友,曹锋,邢子文,等.吸气喷液对空气源热泵热水器性能的影响[J].西安交通大学学报,2008,42(07):818.[doi:10.7652/xjtuxb200807008]
 FEI Jiyou,CAO Feng,XING Ziwen,et al.Influence of Suction Stream Liquid Refrigeration Injection on the Performance of an AirSource Heat Pump Water Heater[J].Journal of Xi'an Jiaotong University,2008,42(06):818.[doi:10.7652/xjtuxb200807008]
[4]冯诗愚,高秀峰,刘卫华,等.内啮合摆线转子压缩机型线连续性条件研究[J].西安交通大学学报,2009,43(01):034.[doi:10.7652/xjtuxb200901008]
 FENG Shiyu,GAO Xiufeng,LIU Weihua,et al.Continuity Condition of Profiles for InternalMeshing Cycloid Rotary Compressor[J].Journal of Xi'an Jiaotong University,2009,43(06):034.[doi:10.7652/xjtuxb200901008]
[5]武晓昆,陈文卿,周明龙,等.双螺杆制冷压缩机气流脉动衰减器的研究与开发[J].西安交通大学学报,2017,51(04):023.[doi:10.7652/xjtuxb201704004]
 WU Xiaokun,CHEN Wenqing,ZHOU Minglong,et al.Development of Gas Pulsation Damping Device for Twin Screw Refrigeration Compressor[J].Journal of Xi'an Jiaotong University,2017,51(06):023.[doi:10.7652/xjtuxb201704004]

备注/Memo

备注/Memo:
收稿日期: 2020-01-02。作者简介: 王炳棋(1996—),男,硕士生; 畅云峰(通信作者),男,副教授,硕士生导师。
更新日期/Last Update: 2020-06-10