双向进气型脉冲管膨胀机能流分析研究_贾启明.pdf

2023 年第 1 期总第 251 期低温工程CRYOGENICSNo.1 2023SumNo.251双向进气型脉冲管膨胀机能流分析研究贾启明1,2龚领会1,2崔文慧1,2,3朱伟平1,2谢秀娟1,2李正宇1,2张宇1,2张梅梅1,2(1航天低温推进剂技术国家重点实验室 100190北京)(2中国科学院低温工程学重点实验室(理化技术研究所)100190北京)(3中国科学院大学 100049北京)摘要:通过建立双向进气型的脉冲管膨胀机系统焓流模型,分析识别了能流在系统内部流向及变化规律,构建了热端换热器换热物理模型,分析出双向进气结构在模型中对热端换热器的焓流熵流影响。建立系统内小孔阀计算模型,并开展实验测试,获得焓流随时间变化关系。该研究有助于其进行性能预测。关键词:脉冲管双向进气能流分析小孔阀中图分类号:TB653文献标识码:A文章编号:1000-6516(2023)01-0006-07收稿日期:2022-11-26;修订日期:2023-02-08基金项目:国家重点研发计划(项目号:2020YFB1506201)资助。作者简介:贾启明,男,31 岁,博士研究生,助理研究员。通讯作者:龚领会,男,60 岁,研究员。Research on energy flow analyses of double-inletpulse tube expander systemJia Qiming1,2Gong Linghui1,2Cui Wenhui1,2,3Zhu Weiping1,2Xie Xiujuan1,2Li Zhengyu1,2Zhang Yu1,2Zhang Meimei1,2(1State Key Laboratory of Technologies in Space Cryogenics Propellants,Beijing 100190,China)(2Chinese Academy of Sciences Key Laboratory of Cryogenics,Technical Institute of Physics and Chemistry,Beijing 100190,China)(3University of Chinese Academy of Sciences,Beijing 100049,China)Abstract:A double-inlet pulse tube expander system model is established.The energy flowdirection and variation law inside the system are analyzed.A physical heat-exchange model at thehot end of the heat exchanger is built.The influence of double inlet structure on the enthalpy flowand entropy flow at the hot end exchanger is analyzed.The calculate model of the orifice valve inthe system is establised,and the experimental tests is carried out to obtain the relationship of en-thalpy flow with time.This research contributes to the performance prediction of a double-inletpulse tube expander system.Key words:pulse tube;double-inlet;energy flow;orifice valve1引言低温制冷系统在新能源利用、稀有气体提取等诸多领域都具有广泛应用,然而关于液氦温区提供十到百瓦级制冷量的低温制冷系统的研究鲜有报道。本研究开展适用于此冷量范围的脉冲管膨胀机研究。第 1 期双向进气型脉冲管膨胀机能流分析研究现有 2 级 GM 脉管制冷机可达液氦温区,文献列举了相关回热式制冷机研究进展1-6。以上回热式循环在约 40 K 以上温区效率高于间壁式循环,间壁式循环在 10 K 以下 温区显示出 独 特 优 势7。周 远等8提出了一种低温旋转阀脉冲管制冷机,采用间壁式换热器替换脉冲管制冷机中的回热器,选用冷端旋转阀实现高低压流路切换。研究发现当高低压压比为 11 时,绝热效率达最大值 42%,在使用液氮预冷后可实现 79 K 的最低制冷温度;分析发现影响膨胀效率的关键因素在于内部换热效率、低温旋转阀流动效率以及调相结构设计。De Waele 等9提出一种反向流脉管制冷机,其内部由两个子系统组成,各子系统中脉管共用旋转阀实现内部压力调节,实验中该系统获得 215 K 的最低制冷温度;其数值研究发现周期运动过程中脉冲管内部存在气体活塞,并模拟获得了 105 K 的最低制冷温度。朱绍伟10提出了一种脉冲管膨胀机的发明专利,其包括热端调相器、脉管、高压气体口等结构,采用多孔结构保证冷端待膨胀气体均匀流入脉冲管,避免湍流对内部气体活塞的破坏影响,通过脉冲管中部留存的气体活塞实现冷端流体膨胀降温。回热式循环经过多年发展在较小制冷量的工况下具备效率高的优势,间壁式低温循环在十到百瓦级制冷量范围中具有优势。对系统开展能量流动分析是探究其内部机理并实现优化的重要方式,以往研究已建立了较为完善的焓流理论11。基于脉冲管膨胀机系统内无回热器的特点,将焓流理论应用于双向进气型脉冲管膨胀机系统研究,可定性分析出不同部件对系统的性能影响,定量计算出理想条件下不同部件内的焓流及熵流值,该方法有助于分析系统内部不同组件对膨胀性能的影响方式,因此该方法具有明显优势。本研究基于脉冲管膨胀机已有研究12,针对双向进气结构对带有脉冲管膨胀机的间壁式循环影响机理展开研究,开展双向进气型脉冲管膨胀机焓流分析,研究了冷端流体流入脉冲管释放焓流的过程,表明冷端阀门流通性能对膨胀效率存在影响。能流分析及实验计算有助于对双向进气型脉冲管膨胀机进行性能分析及预测。2焓流理论建模根据双向进气型脉冲管膨胀机系统内部各部分进行分析研究,并在其中建立质量守恒、动量守恒以及能量守恒方程。本研究将双向进气型脉冲管膨胀机系统划分为7 个控制体单元,结构如图 1 所示。控制体包含压缩机系统 1,控制体包含高压入流 2、低压回流 3、冷端换热器 4,控制体包含冷端入口阀 5、冷端回气阀 6、热端进气阀 14、热端回气阀 15,控制体包含脉管 8,控制体包含层流化元件 7,控制体包括热端换热器 9,控制体包括小孔阀 10。低压氦气经压缩机系统后变为高压氦气,随后分为两股流体分别进入热端进气阀和冷端入口阀,冷端流体流入脉冲管膨胀,随后经冷端回气阀流出脉冲管,冷端低温流体在冷端换热器处吸热随后经换热器流出。热端流体经热端换热器及热端回气阀流回,并与冷端低温流体会和流回压缩机。对各控制体单元开展时均值计算,计算公式中的下角标 1-15 分别表示图 1 中数字对应的模型部件。图 1双向进气型脉冲管膨胀机能流分析计算模型1.压缩机系统;2.高压入流;3.低压回流;4.冷端换热器;5.冷端入口阀;6.冷端回气阀;7.层流化元件;8.脉管;9.热端换热器;10.小孔阀;11.气库;12.高压热端进气;13.低压热端回气;14.热端进气阀;15.热端回气阀。Fig.1Model for calculating enthalpy,entropyand exergy flow of double inlet valve PTE2.1控制体针对控制体进行分析,可获得式(1)(4),其中焓流流向与下标数字顺序相符,并假设内部放出热量设为正值。H1011=10?m1011h(T11)dt(1)S1011=10?m1011s(p11,T11)dt(2)H910=H1011+Q10(3)S10+S910=Q10T10+S1011(4)式中:?m 为工质的质量流率,kg/s;H 为时均焓流,W;S 为时均熵流,W/K;Q 为时均热流,W;T 为温度,K;h 为工质焓值,J/kg;s 为工质熵值,7低温工程2023 年J/(kgK);p 为压力,Pa;为周期时间,s;t 为时间,s。将式(4)化简可得式(5)。S10=?m910cplnT11T9+R?m910lnp910p0-R?m910lnp11p0(5)式中:cp为定压比热,J/(kgK);R 为气体常数,J/(kgK)。计算过程基于等熵流动等假设7,稳态运行时气库部分无环路振荡,周期内进出气库质量流率为 014。在控制体左边界,压力变化为一阶变化,温度为定值,可获得式(6)和(7)。H1011=0(6)S1011=0(7)2.2控制体假设脉冲管内部为等熵流动12,可获得式(8)(10)。H89+H129=H913+H910+Q9(8)S129+S89=S913+S910+Q9T9(9)Q9=T9(?m129-?m913+?m89-?m910)s(p0,T0)+cplnT9T0()+T9R?m913lnp913p0(-?m129lnp129p0+?m910lnp910p0-?m89lnp89p0)(10)当热端进回气质量流率相等时,可获得式(11)和式(12),表明冷端焓流经脉管从热端换热器处流出系统。Q9=T9R?m913lnp913p0-?m912lnp912p0+?m910lnp910p0()(11)H89=Q9(12)2.3控制体在稳定流动周期内,控制体中的冷端层流化元件无吸热效果,可建立方程获得式(13)和(14)。H76=H57-H89(13)S57+S129=S76+S910+S913+Q9T9(14)当热端进回气质量相等时,可得式(15),表明此时内部实现等熵膨胀。s(p5,T5)=s(p6,T6)(15)2.4控制体在控制体内,有式(16)和式(17)。S57=S76+S78(16)S78=(?m57-?m76)sp0,T0()+cplnT9T0()-(?m57-?m76)Rlnp89p0(17)化简可得式(18)。?m57 cplnT5T9-R?m57lnp57p0=?m76 cplnT6T9-R?m76lnp76p0+?m57lnp89p0-?m76lnp89p0()(18)当冷端流入流出脉管质量相等时,可获得下式(19),表明脉冲管冷端无熵流流入流出。S78=0(19)2.5控制体由于阀门开闭过程极其短暂,工作流体的熵流和焓流在阀门前后相等,可获得式(20)和式(21)。H57=10?m57h(T5)dt(20)H76=10?m76h(T6)dt(21)2.6控制体假设控制体的换热器中无换热损失,可计算获得式(22),在冷端换热器中一部分熵流伴随热量流入系统,可计算得出冷端换热量式(23)。H12+H64=H31+H25+Q4(22)Q4=T4(?m12-?m31)s(p0,T0)+cplnT1T0()+T4(?m64-?m25)s(p0,T0)+cplnT5T0()+T4R?m31lnp31p0-?m12lnp12p0+(?m25lnp25p0-?m64lnp64p0)(23)当流入流出脉冲管质量相等时,可获得式(24)。Q4=-TaR?m913lnp913p0-?m912lnp912p0+?m910lnp910p0()(24)8第 1 期双向进气型脉冲管膨胀机能流分析研究2.7控制体压缩机系统包括压缩机和水冷器,流体在压缩机系统中实现等温压缩。可推导出式(25)和式(26)。S31+S151=S12+S114+Q1T1(25)Q1T1=-Q4T4-Q9T9-S910(26)当流入流出脉冲管质量相等时,获得水冷器放热量如式(27)所示和压缩机压缩功如式(28)所示。其中 W 为时均功流,W。Q1=T1TaRT4?m913lnp913p0-?m912lnp912p0+(?m910lnp910p0)-T1R?m913lnp913p0-?m912lnp912p0()(27)W=R Ta0?m1lnp12p13dt(28)因此理想情况下焓流和熵流在双向进气型脉冲管膨胀机系统中变化如图 2 所示。首先流体经压缩机压缩后流出并分为两路,分别为流向换热器的制冷流路和流向脉冲管热端的调相流路。两股流体各自携带焓流,其中调相流路内焓流大小取决于热端调相阀开度和流路流量,制冷流路的焓流变化取决于换热器换热效率和流路流量,如文献12 中所述。热量在恒温的热端换热器处向外界释放,因此在热端换热图 2双向进气型脉冲管膨胀机焓流和熵流变化示意图Fig.2Change of enthalpy flow and entropy flowfor double inlet valve PTE器中流出与流入的焓流量相等。随后调相流路焓流和制冷流路焓流合并后流回压缩机,完成循环。对熵流分析发现,双向进气型脉冲管膨胀机的熵产主要发生在小孔阀处,其产生后一部分在热端换热器处以熵流的形式流出系统,另一部分在压缩机水冷器处以热量形式流出系统。热端的高压入流与低压回流减小了热端换热器的换热负担,因此双向进气型通常较基本型有更好的性能。而双向进气带回压缩机的熵流同样以热量的形式在水冷器处向外界释放,此过程增加了水冷器的换热负担。小孔阀处熵产是系统的主要熵产来源,其内部焓流变化也直接影响热端焓流分布情况,因此对小孔阀开展进一步实验分析研究,验证计算小孔阀内部焓流及熵流变化,有助于提高系统性能。3小孔阀分析建模针对小孔阀开展理论计算,将阀门开启闭合过程中的流道变化视为喉部直径可变的缩放喷管模型。阀门内部流通面积可视为缩放喷管的喉部截面积并可基于阀芯结构计算获得。其中当阀门开启时间为 0.236 s 时流通面积最大,此时为 41.8 mm2。将喉部截面积随时间变化关系写成式(29)。Athroat(t)=at2+bt1+c(29)式中:a 为二 次系数,mm2/s2;b 为一次系数,mm2/s;c 为常数值,mm2;Athroat为喉部截面积,mm2。对小孔阀及阀前腔体进行进一步研究,可以获得图 3,对其中不同部分进行离散,获得 ai 的离散单元,其中 ac 为阀前单元,d,e 和 f 为阀门单元,gi为阀后容积单元。每部分控制体中心点视为模块内物性参数点,模块间交接面为质量流、焓流与熵流等交换界面,内部中心点参数随时间变化。在计算过程中假设以下条件:控制体内部与外部之间为绝热状态;阀门模块流动截面随时间发生变化;阀后容积模块出口为常压;阀前容积模块入口封闭。图 3小孔阀数值简化模型Fig.3Simplified energy flow model ofvalve experiment part9低温工程2023 年基于守恒方程对小孔阀进行建模。基于能流计算方程,可获得任意模块内质量流式(30)和焓流式(31)。?mij=1tijuijAijdt(30)Hij=10?mijhTij()dt(31)式中:为工质密度,kg/m3;u 为流速,m/s;A 为流通面积,m2。将小孔阀内初始物性参数带入模型,可计算获得内部不同控制体单元压力随时间变化关系。4小孔阀实验测试理论分析发现小孔阀对系统内部焓流变化会产生重要影响,为研究其实际影响效果,本研究开展阀门内部流动及焓流变化实验研究。为实现对内部压力及焓流的有效分析,基于小孔阀测试平台,对阀门开闭时序进行 调控,采用传 感器 对 阀 前 压 力 进 行检测。小孔阀测试流程示意图如图 4 所示,通过钢瓶对测试阀前容积进行充气,其中阀前容积为 0.194 L,利用控制面板实现电机运转调控阀门开闭时序,在气库和测试阀门间布置有压力传感器,用于测量阀前后压力随时间的变化关系。假设腔体内部压力为均匀变化,采用气体均一性假设及压力脉冲传播模型15获得腔体内部压力数值。开展阀门压力测试实验。针对阀门开展入口压力 pin分别为 1 MPa、1.5 MPa 和2 MPa,出口背压为 0.1 MPa 的性能实验,阀门开启频率为 1 Hz。图 4小孔阀测试流程示意图Fig.4Schematic diagram of small orifice valve testing阀门开启后记录腔内压力波随时间的变化关系,并通过泛函分析法数值拟合,获得阀前压力变化曲线。阀前压力变化曲线如图 5 所示。其中压降百分比 按照式(32)定义为瞬时压力 pv与 pin的比值。结果表明阀前压力在阀门开启后压降速率呈先增后减的趋势,当阀前压力为 1 MPa 时,阀门开启 0.4 s 后,压力降低 85%至 0.15 MPa,并最终稳定于 0.123 MPa。=pvpin(32)图 5不同压力情况下的阀前压力变化曲线Fig.5Graph of pressure variation beforevalve with different inlet pressure5结果分析通过带入阀芯流通面积可对小孔阀数值模型进01第 1 期双向进气型脉冲管膨胀机能流分析研究行完善,进而实现对阀内压力的计算,获得不同初始压力情况下阀前压力随时间的变化关系,如图 6 所示。图中压力计算数据与实验数据拟合较好,阀前压力都随着小孔阀开启而逐渐降低,且降低速率先升高后降低,表明计算模型可靠。图 6不同实验中实验拟合值与计算值对比Fig.6Graph of calculated pressure and experimentalpressure before valve with different inlet pressure基于实验数据和式(30)和式(31),获得流经小孔阀的焓流为 346.31 W。而计算结果表明在2 MPa阀前压力的情 况下,开启时 流经 小 孔 阀 的 焓 流 为532.8 W,表明实际过程中,阀处损失较大或实验过程中存在部分误差。进而分析误差来源,主要是压力传感器数据采集、数值拟合过程中的数据误差和采用均一性假设过程导致的误差。压力传感器的测量精度为 0.2 级,测量误差在 0.2%。而数值拟合过程中的数据误差来源是多组数据平均以及拟合方法选用导致的拟合偏差,分析发现最大偏差为 5.71%,拟合优度大于 0.9 994,认为拟合过程误差在允许范围内。因此判断小孔阀内存在流动损失,进而采用该模型对阀内焓流进行分析计算,对比不同压力下阀内焓流变化。如图 7 所示,通过计算分别获得不同初始压力下小孔阀阀前单元和阀后单元焓流随时间变化关系。发现在阀门开启 0.1 s 后内部焓流基本保持稳定。当出现第一种情况,即阀前初始压力 2 MPa 时,阀前后压差较大,内部流速较快,且阀前气体焓值较高,此时阀前焓流较大;阀后流体初始流速较小但流速逐渐增加,因此阀后焓流逐渐增大。当出现第 2 种情况,即阀前初始压力为 1 MPa 时,由于阀前后压差减小,内部流速较第一种情况低,且阀前气体焓值低于前者,因此阀前焓流低于前者;阀后焓流与第一种情况相似,随着速度升高逐渐增加。计算结果表明初始情况下阀前阀后压差越大则单方向流动期间流过阀内的焓流越多。图 7不同压力下阀内焓流变化Fig.7Graph of enthalpy flow inner valve withdifferent inlet pressure6结论本研究基于焓流理论建立了带有小孔及双向进气的脉冲管膨胀机系统模型,分析识别了焓流与能流在系统内部的流向及变化规律。研究表明在理想情况下,流体在小孔阀处产生不可逆熵产,其以熵流的形式从脉冲管热端换热器伴随热量流出系统。当开启热端双向进气阀后,热端进回气将部分不可逆熵产从脉管热端带回压缩机,减小热端放热量,降低换热负担。冷端高压流体流入脉冲管膨胀降压随后流出脉冲管,膨胀过程产生部分焓流流向脉冲管热端,以热量形式释放。进一步基于测试平台对小孔阀焓流数值研究开展分析研究,建立小孔阀压降及焓流分析模型,并结合实验数据分析获得开启过程中小孔阀内焓流变化。参考文献1 Gifford W E,Longsworth R C.Pulse-tube refrigerationJ.Journal ofEngineering for Industry,1964,86(3):264-268.2 Snodgrass R,Ullom J.Direct cooling from the regenerators of Gifford-McMahon cryocoolers with comparison to pulse tube refrigeratorsJ.Cryogenics,2022,124:103473.3 Feng T S,Chen H L,Liu X T,et al.Effect of nitrogen contaminationon the performance of pulse tube cryocoolers J.Applied ThermalEngineering,2022,201:117810.4 Pang X M,Wang H,Wang X T,et al.Numerical investigation on theinfluence of radial thermal conduction in a co-axial pulse tube coolerJ.International Journal of Refrigeration,2022,139:128-135.5 Guo Z M,Pfotenhauer J M,Miller F K,et al.A research of micro-pulse tube cryocooler with displacer phase shifter J.Applied Ther-11低温工程2023 年mal Engineering,2022,205:117995.6 Cao Q,Qiu L M,Gan Z H.Real gas effects on the temperature profileof regenerators J.Cryogenics,2014,61:31-37.7 党海政,张涛,赵帮健,等.以氦-4 为唯一工质的 1.8 K 复合制冷机及其应用验证J.科学通报,2022,67(9):896-905.Dang Hanzheng,Zhang Tao,Zhao Bangjian,et al.A hybrid cryocoolerachieving 1.8 K with He-4 as the only working medium and its appli-cation verification J.Chinese Science Bulletin,2022,67(9):896-905.8 卞诗瑶,梁惊涛,杨鲁伟,等.低温旋转阀脉冲管制冷机实验研究J.工程热物理学报,2000(1):5-8.Bian Shiyao,Liang Jingtao,Yang Luwei,et al.Study of pulse tuberefrigerators with low temperature switching valve J.Journal of En-gineering Thermophysics,2000(1):5-8.9 Will M E,de Waele A T A M.Analytical treatment of counterflowpulse-tube refrigerators J.Cryogenics,2006,46(6):421-431.10朱绍伟.一种脉管膨胀机 P.中国专利,CN201310058359.X.2013-07-17.Zhu Shaowei.Pulse Tube ExpanderP.China,CN201310058359.X.2013-07-17.11Kuriyama F,Radebaugh R.Analysis of mass and energy flow rates inan orifice pulse-tube refrigerator J.Cryogenics,1999,39(1):85-92.12Jia Q M,Gong L H,Feng G C,et al.Thermodynamic analyses and theexperimental validation of the pulse tube expander system J.Cryo-genics,2018,91:118-124.13Jia Q M,Li Z Y,Gong L H,et al.Experimental investigation and nu-merical calculation of the cryogenic ejector in a liquid nitrogen systemJ.Applied Thermal Engineering,2021,184:116322.14罗二仓,刘 浩,吴剑锋.回热式热机中计算时均焓流的新公式兼论 Radebaugh 关于 脉 冲 管 中 时 均 焓 流 计 算 的 重 要 错 误J.低温与超导,2001,29(3):8-11.Luo Ercang,Liu Hao,Wu Jianfeng.A new formula for calculating av-erage enthalpy flow in regenerative engineJ.Cryogenics and Super-conductivity,2001,29(3):8-11.15李翠,武定航,杜飞平,等.液氧管路燃气射流冷凝压力脉动传播特性研究J.西安交通大学学报,2022,56(7):127-135.Li Cui,Wu Dinghang,Du Feiping,et al.Propagation of pressure fluc-tuation in jet condensing flow of liquid oxygen feed system J.Jour-nal of Xian Jiaotong University,2022,56(7):127-135.(上接第 6 页)参考文献1 陈国邦,汤珂.小型低温制冷机原理M.北京:科学出版社,2009.Chen Guobang,Tang Ke.The Principle of Miniature RefrigeratorM.Beijing:Science Press,2009.2 Wang N,Zhao M,Ou Y,et al.A high efficiency coaxial pulse tubecryocooler operating at 60KJ.Cryogenics,2018:48-50.3 连新昊,高慧莲.制冷技术及 其在航天遥感领 域的应用综述C.第二十三届全国空间探测学术交流会论文集,2010:1-9.Lian Xinhao,Gao Huilian.Summarize of refrigeration technology andapplication in space remote sensingC.Proceedings of the 23rd Na-tional Symposium on Space Exploration,2010:1-9.4 Liu S S,Chen X,Zhang A K,et al.Investigation on phase shifter of a10 W/70K inertance pulse tube refrigeratorJ.International Journalof Refrigeration,2017,74:448-455.5 Potratz S A,Nellis G F,Maddocks J R,et al.Design and test of 300W,65K pulse tube C.Asme International Mechanical EngineeringCongress&Exposition,2004.6 Potratz S A,Abbott T D,Johnson M C,et al.Stirling-type pulse tubecryocooler with 1 kW of refrigeration at 77KJ.American Institute ofPhysics,2008,985(1):42-48.7 Hu J Y,Zhang L M,Zhu J,et al.A high-efficiency coaxial pulsetube cryocooler with 500 W cooling capacity at 80 KJ.Cryogenics,2014,62:7-10.8 Dang H,Tan J,Zha R,et al.Review of recent advances in Stirling-type pulse tube cryocoolersJ.IOP Conference Series Materials Sci-ence and Engineering,2019,502(1):012034.9 Chao G U,Tang J B,Wang J J,et al.Advance in research of severaltypes of streaming of pulse tube refrigeratorsJ.Science China Tech-nological Sciences,2013(11):2690-2701.10顾超.脉冲管制冷机第三种直流效应研究D.北京:中国科学院大学,2013.Gu Chao.Studies on the third type of streaming in pulse tube refriger-ator D.Beijing:The University of Chinese Academy of Sciences,2013.11Dietrich M,Yang L W,Thummes G.High-power Stirling-type pulsetube cryocooler:Observation and reduction of regenerator temperature-inhomogeneitiesJ.Cryogenics,2007,47:306-314.12Andersen S K,Dietrich M,Carlsen H.Numerical study on transverseasymmetry in the temperature profile of a regenerator in a pulse tubecoolerJ.Int J Heat Mass Transfer,2007,50:2795-2804.13刘冬辉,邱利民,甘智华,等.基于并联回热器模型的大功率斯特林脉管制冷机回热器温度不均匀性研究 J.低温 工程,2011(6):1-5.Liu Donghui,Qiu Limin,Gan Zhihua,et al.Study on temperatureinhomogeneity in regenerator of high-power Stirling-type pulse tubecryocoolers based on two parallel regenerator modelJ.Cryogenics,2011(6):1-5.14王乃亮,赵密广,陈厚磊,等.72 Hz 大冷量空间脉冲管制冷机J.科学通报,2018,63(36):3968-3974.Wang Nailiang,Zhao Miguang,Chen Houlei,et al.Study of a 72 Hzhigh-capacity pulse tube cryocooler for space applications J.Chi-nese Science Bulletin,2018,63(36):3968-3974.21