基于脉管制冷机预冷的1 K大冷量低温系统

上一篇

下一篇

刘旭明, 查奎帆, 马帅, 韩丽明, 谢晓霖, 郭伟杰, 潘长钊. 基于脉管制冷机预冷的1 K大冷量低温系统[J]. 物理学报, 2025, 74(11): 110701-1. doi: 10.7498/aps.74.20250181
引用本文: 刘旭明, 查奎帆, 马帅, 韩丽明, 谢晓霖, 郭伟杰, 潘长钊. 基于脉管制冷机预冷的1 K大冷量低温系统[J]. 物理学报, 2025, 74(11): 110701-1. doi: 10.7498/aps.74.20250181
Xuming LIU, Kuifan ZHA, Shuai MA, Liming HAN, Xiaolin XIE, Weijie GUO, Changzhao PAN. A high-capacity 1-K cryogenic system pre-cooled by pulse tube cryocooler[J]. Acta Physica Sinica, 2025, 74(11): 110701-1. doi: 10.7498/aps.74.20250181
Citation: Xuming LIU, Kuifan ZHA, Shuai MA, Liming HAN, Xiaolin XIE, Weijie GUO, Changzhao PAN. A high-capacity 1-K cryogenic system pre-cooled by pulse tube cryocooler[J]. Acta Physica Sinica, 2025, 74(11): 110701-1. doi: 10.7498/aps.74.20250181

基于脉管制冷机预冷的1 K大冷量低温系统

    作者简介: 刘旭明.E-mail: liuxm3@sustech.edu.cn .
    通讯作者: E-mail: pancz@iqasz.cn.
  • 中图分类号: 07.20.Mc

A high-capacity 1-K cryogenic system pre-cooled by pulse tube cryocooler

    Corresponding author: E-mail: pancz@iqasz.cn
  • MSC: 07.20.Mc

  • 摘要: 1 K低温系统是进一步实现mK温区及更低温度的基础, 广泛应用于量子计算、凝聚态物理研究、低温科学仪器等领域. 目前国内的1 K低温系统大多使用GM (Gifford-McMahon)制冷机进行预冷, 系统在实现更低振动控制、更低电噪声干扰、更低预冷温度和更高液化效率等方面存在一定难题, 而基于脉管制冷机预冷的1 K系统在解决这些问题方面具有先天优势. 本文发展了一台全国产化的4 K GM脉管制冷机, 获得了2.14 K的最低制冷温度, 并可同时提供1.5 W@4.2 K和45 W@45 K的制冷量. 将其作为预冷制冷机, 设计并搭建了1 K低温系统, 最终获得了1.1 K的最低制冷温度, 并可在1.6 K提供100 mW的制冷量. 本研究为后续开展更大冷量稀释制冷技术奠定了重要基础.
  • 加载中
  • 图 1  工质氦4压焓图

    Figure 1.  The pressure-enthalpy diagram of the helium-4.

    图 2  基于脉管制冷机预冷的1 K低温系统结构示意图

    Figure 2.  Schematic diagram of the 1 K cryogenic system pre-cooled by a pulse tube refrigerator.

    图 3  实验系统实物照片

    Figure 3.  Photograph of the experimental system.

    图 4  1 K气路系统流程示意图

    Figure 4.  Schematic diagram of the 1 K gas circuit system.

    图 5  4 K脉管制冷机典型制冷性能 (a)降温曲线, (b)制冷量

    Figure 5.  Typical cooling performance of the developed 4 K pulse tube refrigerator: (a) Cooling curve; (b) cooling capacity.

    图 6  1 K低温系统典型制冷性能 (a)降温曲线; (b)制冷量

    Figure 6.  Typical cooling performance of the 1 K cryogenic system: (a) Cooling curve; (b) cooling capacity.

    图 7  制冷量和工质流量与制冷温度关系变化曲线

    Figure 7.  The relationship curve between cooling capacity, working flow and cooling temperature.

    表 1  自研脉管制冷机与国外产品比较

    Table 1.  Comparison between the developed prototype and foreign 4 K GM-type PTRs.

    时间 最低温度 一级制冷量 二级制冷量 功耗 备注
    Cryomech PT415-RM <60 min <2.8 K 40 W @ 45 K 1.35 W @ 4.2 K 9.2 kW 阀分离
    住友RP-182B2S <60 min <2.8 K 36 W @ 48 K 1.5 W @ 4.2 K 11.8 kW 阀分离
    本文 <40 min 2.14 K 45 W @ 45 K 1.5 W @ 4.2 K 14 kW 阀分离
    下载: 导出CSV
  • [1] Zhao Z Y, Wang C 2019 Cryogenic Engineering and Technologies: Principles and Applications of Cryogen-Free Systems (CRC Press) p233
    [2] 李珂, 王亚男, 刘萍, 禹芳秋, 戴巍, 沈俊 2023 物理学报 72 190702 doi: 10.7498/aps.72.20231102 Li K, Wang Y N, Liu P, Yu F Q, Dai W, Shen J 2023 Acta Phys. Sin. 72 190702 doi: 10.7498/aps.72.20231102
    [3] 俎红叶, 程维军, 王亚男, 王晓涛, 李珂, 戴巍 2023 物理学报 72 080701 doi: 10.7498/aps.72.20222257 Zu H Y, Cheng W J, Wang Y N, Wang X T, Li K, Dai W 2023 Acta Phys. Sin. 72 080701 doi: 10.7498/aps.72.20222257
    [4] 王昌, 李珂, 沈俊, 戴巍, 王亚男, 罗二仓, 沈保根, 周远 2021 物理学报 70 090702 doi: 10.7498/aps.70.20202237 Wang C, Li K, Shen J, Dai W, Wang Y N, Luo E C, Shen B G, Zhou Y 2021 Acta Phys. Sin. 70 090702 doi: 10.7498/aps.70.20202237
    [5] Zheng M W, Guo H W, Wei L J, Pan Z J, Zou J R, Li R X, Zhao M G, Chen H L, Liang J T 2024 Acta Phys. Sin. 73 230701 (In Chinese) [郑茂文, 郭浩文, 卫铃佼, 潘子杰, 邹佳润, 李瑞鑫, 赵密广, 陈厚磊, 梁惊涛 2024 物理学报 73 230701] doi: 10.7498/aps.73.20241211 Zheng M W, Guo H W, Wei L J, Pan Z J, Zou J R, Li R X, Zhao M G, Chen H L, Liang J T 2024 Acta Phys. Sin. 73 230701 (In Chinese) doi: 10.7498/aps.73.20241211
    [6] Guan X, Fan J, Bian Y B, Cheng Z G, Ji Z Q 2024 Chin. Phys. B 33 070701 doi: 10.1088/1674-1056/ad3b87
    [7] Jahromi A E, Miller F K 2014 Cryogenics 61 15 doi: 10.1016/j.cryogenics.2014.02.001
    [8] DeMann A, Mueller S, Field S B 2016 Cryogenics 73 60 doi: 10.1016/j.cryogenics.2015.11.008
    [9] Cao H 2021 J. Low Temp. Phys. 204 175 doi: 10.1007/s10909-021-02606-7
    [10] Bluefors and Cryomech 1 K systems products https://bluefors.com/products/1k-systems/
    [11] Oxford Instruments 1 K cryostats products https://nanoscience.oxinst.com/dry-systems/products/teslatronpt
    [12] Quantum Design 1 K measurement systems products https://www.qd-china.com/zh/pro/detail3/1/1912091422155/1909260926498
    [13] Wang L G, Qu Q X, Chen H, Dai N N, Zhao W Y, Jia P, Xu D, Li L F 2025 Cryogenics 145 103992 doi: 10.1016/j.cryogenics.2024.103992
    [14] Pengli 1 K cryostats products https://isite.baidu.com/site/wjzru1zo/96a2066c-8800-4095-9ec4-a0489538571f?ch=48&wid=6dfbf96df3554e288101d75dc1ec8912_0_0&uniqId=c4db65aa10714a8697008f3c034cf058
    [15] ZL Cryogenic 1 K cryostats products http://zlcryogenic.com/display/141459.html
    [16] Radebaugh R. 2009 J. Phys. Condens. Matter 21 164219 doi: 10.1088/0953-8984/21/16/164219
    [17] Liu X M, Chen L B, Wu X L, Yang B, Wang J, Zhu W X, Wang J J, Zhou Y 2020 Sci. China Technol. Sci. 63 434 doi: 10.1007/s11431-019-1471-7
    [18] Wang C 2016 Cryocoolers 19 299
    [19] Wang C, Hanrahan T, Johnson M 2018 Cryogenics 95 64 doi: 10.1016/j.cryogenics.2018.07.010
    [20] Wang C, Lichtenwalter B, Friebel A, Tang H X 2014 Cryogenics 64 5 doi: 10.1016/j.cryogenics.2014.07.013
    [21] Qu Q X, Wang L G, Chen H, Dai N N, Jia P, Xu D, Li L F 2024 Cryogenics 138 103797 doi: 10.1016/j.cryogenics.2024.103797
    [22] Wu S G, Zhao B J, Tan J, Zhao Y J, Zhai Y J, Xue R J, Tan H, Ma D, Wu D R, Dang H Z 2023 Energy 277 127691 doi: 10.1016/j.energy.2023.127691
    [23] Shen Y W, Liu D L, Chen S F, Zhao Q Y, Liu L, Gan Z H, Qiu M 2020 Appl. Therm. Eng. 166 114667 doi: 10.1016/j.applthermaleng.2019.114667
    [24] Li X, Xu D, Wang W, Lin P, Liu H M, Nishimura A, Shen F Z, Li L F 2019 Cryogenics 102 50 doi: 10.1016/j.cryogenics.2019.06.007
    [25] Uhlig K 2002 Cryogenics 42 73 doi: 10.1016/S0011-2275(02)00002-4
    [26] Wang C 2001 Cryogenics 41 491 doi: 10.1016/S0011-2275(01)00116-3
    [27] Liu X M, Pan C Z, Zhang Y, Liao Y, Guo W J, Yu D P 2023 Acta Phys. Sin. 72 190701 (In Chinese) [刘旭明, 潘长钊, 张宇, 廖奕, 郭伟杰, 俞大鹏 2023 物理学报 72 190701] doi: 10.7498/aps.72.20230910 Liu X M, Pan C Z, Zhang Y, Liao Y, Guo W J, Yu D P 2023 Acta Phys. Sin. 72 190701 (In Chinese) doi: 10.7498/aps.72.20230910
  • 加载中
图( 7) 表( 1)
计量
  • 文章访问数:  30
  • HTML全文浏览数:  30
  • PDF下载数:  1
  • 施引文献:  0
出版历程
  • 收稿日期:  2025-02-15
  • 刊出日期:  2025-06-05

基于脉管制冷机预冷的1 K大冷量低温系统

    通讯作者: E-mail: pancz@iqasz.cn.
    作者简介: 刘旭明.E-mail: liuxm3@sustech.edu.cn
  • 1. 深圳国际量子研究院, 深圳 518048
  • 2. 深圳鲲腾酷领科技有限公司, 深圳 518048

摘要: 1 K低温系统是进一步实现mK温区及更低温度的基础, 广泛应用于量子计算、凝聚态物理研究、低温科学仪器等领域. 目前国内的1 K低温系统大多使用GM (Gifford-McMahon)制冷机进行预冷, 系统在实现更低振动控制、更低电噪声干扰、更低预冷温度和更高液化效率等方面存在一定难题, 而基于脉管制冷机预冷的1 K系统在解决这些问题方面具有先天优势. 本文发展了一台全国产化的4 K GM脉管制冷机, 获得了2.14 K的最低制冷温度, 并可同时提供1.5 W@4.2 K和45 W@45 K的制冷量. 将其作为预冷制冷机, 设计并搭建了1 K低温系统, 最终获得了1.1 K的最低制冷温度, 并可在1.6 K提供100 mW的制冷量. 本研究为后续开展更大冷量稀释制冷技术奠定了重要基础.

English Abstract

    • 极低温制冷技术通常是指能够获取1 K及以下温度的机械式制冷技术, 目前广泛应用于量子计算、凝聚态物理研究、低温科学仪器等领域[16]. 1 K温区作为关键温度区域, 能够为多个领域提供稳定且必要的低温环境[79]: 在基础研究领域, 1 K温区是研究低温物理中量子现象(如量子霍尔效应、拓扑相变等)的理想条件; 在技术应用领域, 1 K温区是一些量子器件(如超导量子干涉仪、单光子探测器等)实现高灵敏度工作的必要条件; 在极低温技术领域, 1 K温区是稀释制冷机等制冷技术的前级预冷阶段, 是进一步实现mK温区及更低温度的基础. 因此, 开发高效、稳定的1 K低温系统对于推动低温科学和先进技术发展至关重要.

      近年来, 国外在1 K低温系统领域取得了显著进展, 特别是在系统效率、低振动设计、系统智能化与模块化方面. 例如, Cryomech, Bluefors, Oxford Instruments, Quantum Design等公司已开发出集成式1 K系统, 并广泛应用于量子计算、精密测量等领域[1012]. 在这些系统中, 4 K GM (Gifford-McMahon)脉管制冷机被用来作为预冷制冷机来对整个系统进行预冷. 在国内, 1 K低温系统的研究近年来也取得了长足进步, 多家科研机构围绕1 K制冷机的国产化、低噪声、高效热交换等关键技术展开深入研究[1315]. 但由于近年来国外4 K GM脉管制冷机产品被全面禁运, 这些系统大多使用国产4 K GM制冷机进行预冷. 但由于GM制冷机自身回热器结构和排出器振动的原因, 使得1 K系统在实现更低振动控制、更低电噪声干扰、更低预冷温度和更高液化效率等方面存在一定困难, 相应的解决方案包括采用主动振动补偿、可调流阻设计等, 但系统工程复杂, 成本较高. 脉管制冷机在解决这些问题方面具有先天优势, 这也是国外1 K系统发展基于脉管制冷机的主要原因. 因此, 开展基于国产脉管制冷机预冷的1 K低温系统, 具有重要的战略价值和现实意义[1626].

      为打破国外技术封锁, 保障我国在量子科技等领域的发展, 本课题组于2021年开展液氦温区大冷量GM脉管制冷机研究, 并于2023年报道了一台自研的脉管制冷机冷头, 使用美国压缩机和旋转阀, 获得了3.1 K的最低温度, 并可在4.2 K提供0.8 W制冷量[27]. 基于此, 我们发展了全国产化的脉管制冷机系统, 并优化提升了其制冷性能. 更进一步地, 我们将其作为预冷制冷机, 设计并搭建了全国产1 K低温系统. 本文将从基本原理、系统设计、关键部件及实验结果等多个方面对该系统进行详细介绍.

    • 1 K低温系统的基本原理是通过液氦蒸发循环实现制冷, J-T(Joule-Thomson)流阻为其关键部件, 如图1所示, 假设工质通过J-T流阻发生绝热膨胀, 因此工质进出流阻前后的焓值相等. 假设经过脉管制冷机预冷, 进入J-T流阻前的工质为过冷液相, 流出J-T流阻后的工质, 其中一部分在蒸发室中以液态形式存在, 另一部分则转变为气态被干泵抽离. 基于能量守恒, 以上过程可表示为

      其中, hinlet为进入J-T流阻工质氦的比焓, 单位为J/kg; h2, l为膨胀后液态氦的比焓, 单位为J/kg; h2, g为膨胀后气态氦的比焓, 单位为J/kg; f 为膨胀后保持液态的氦的质量分数, 1 – f为膨胀后转变为气态的氦的质量分数.

      由上述方程可解得液态氦的质量分数f:

      假设膨胀后的气态氦被立即抽离, 其对制冷量的贡献可忽略不计. 同时, 外界热负荷仅引起蒸发室内液态氦的相变, 而不会导致液体温度的升高. 在稳态条件下, 系统的制冷量可由以下关系式计算:

      其中, $ \dot Q $为制冷量, W, $ \dot m $为工质的质量流量, 单位为kg/s. 由此表明, 1 K低温系统的制冷能力取决于工质氦的质量流量及其在相变过程中释放的汽化潜热. 使用脉管制冷机预冷可以最大程度释放工质氦气的液化显热和潜热, 将其液化至更低的过冷温度. 增大流阻有利于获得更大压差, 实现更低的制冷温度, 但也限制了工质流量的增大, 不利于系统获取更大的制冷量, 因此设计时需要综合考虑. 此外, 当制冷温度足够低, 还会出现超流爬膜现象(superfluid creeping film), 表现为超流氦会沿着容器壁面蠕动, 直至铺满容器表面. 对于1 K系统而言, 超流爬膜会导致: 部分液体不参与蒸发室的蒸发制冷, 降低制冷效率; 爬膜在蒸发室与较高温度的部件间产生热连接, 增大系统热负载, 因此在设计时也要避免该现象.

    • 基于脉管制冷机预冷的1 K低温系统基本结构如图2所示. 其中, 1为干泵, 用于实现气路循环, 关键因素为其抽速; 2为储气罐, 内充工质氦4气体, 根据理论制冷量计算其容积为50 L; 3为自研脉管制冷机, 二级回热器上通过真空钎焊绕制不锈钢盘管, 再通过银焊焊接至铜柱表面, 并螺栓固定在4 K冷头上. 预冷制冷机关键因素为其在4 K温区和40 K温区的制冷性能. 相比先前的报道, 我们进一步使用国产大排量氦气压缩机和自研旋转阀, 并重新设计了冷端换热器结构和室温端调相机构, 得以大幅提升其制冷性能; 4为减振法兰, 以进一步降低系统振动; 5为热桥, 使用铜辫子实现软连接; 6为冷阱, 内填活性炭, 用于过滤工质气体中的杂质; 9为逆流换热器, 使用管套管结构; 10为自研组合式热开关, 该热开关在传统气隙式热开关的基础上加上了铟片和氧化铝垫片. 铟的延展性好, 能够填补固体表面缝隙, 使得在连通时增强导热, 而在极低温区, 铟和氧化铝会逐渐进入深度超导状态, 其热导率大大降低, 能够有效减小热开关断开状态下的残余漏热; 11为J-T流阻, 其主体是用直径0.15 mm的不锈钢丝插入内径0.2 mm的不锈钢管, 通过改变长度使其达到合适阻抗, 实验测量流阻阻抗约为3.8×1010 cm–3, 符合经验值; 12为蒸发室, 根据储气罐容积换算液体体积设计其容积约80 cc. 使用CF16波纹管来对蒸发室进行抽气, 内设刀口防超流结构, 设计时需要综合考虑刀口法兰内径及其距液面的高度. 为了下一步能够适配9 T低温超导磁体, 该样机设计了更大的冷盘尺寸和冷却空间, 其中, 7为一级冷盘, 直径超过450 mm; 8为二级冷盘, 直径超过360 mm; 13为1 K冷盘, 直径超过220 mm. 更大冷盘尺寸和空间, 也导致了系统需要设置更大的辐射屏和承载更多的信号控制及输出线缆, 这对自研脉管制冷机的预冷性能提出了更高要求. 图中吊杆、防辐射屏、多层绝热材料等均未示出.

    • 图3为整机系统实物照片, 主要包含以下4个关键系统.

      1) 4 K预冷系统. 包括一台自研的4 K GM型脉管制冷机, 主要用于对泵送回来的工质氦气进行预冷和液化, 并负责冷却一级和二级防辐射屏, 阻挡外界环境的辐射漏热.

      2) 1 K气路系统. 具体如图4所示, 主要包括干泵、压力计、流量计、电磁阀、手阀、针阀、储气罐等. 整个实验过程分为预冷阶段、液化阶段、常规循环阶段和收气回温阶段. 预冷阶段: 1 K系统主要通过自研的脉管制冷机进行预冷, 脉管利用氦气的西蒙膨胀原理(绝热放气)使得冷头降到液氦温区, 为避免冷头振动影响系统温度和测量, 冷头通过铜辫子的软连接将冷量传递给一级和二级冷盘. 而1 K冷盘通过自研的组合式热开关实现与二级冷盘的连接. 预冷过程伊始, 打开电磁阀V2, V3, V5和V6, 开启干泵, 待储气罐中少量工质气体进入主循环回路后, 关闭V5, 通过循环少量氦气实现加速降温. 启动冷头压缩机, 打开热开关加热器, 当二级冷盘温度低于4 K, 1 K冷盘温度低于5 K时预冷结束, 开始液化. 液化阶段: 通过反复开关V1从气罐取气, 室温氦气先经过40 K冷阱, 少量杂质气体被吸附并进行初步预冷, 接着经过缠绕在脉管制冷机上的盘管, 利用二级回热器的冷量释放显热, 然后经过布置在二级冷头的换热器, 释放潜热实现液化. 最后经过流阻实现进一步降温, 液体累积在蒸发室里, 气体被干泵抽离. 当气罐的压力低于500 Pa时, 液化结束进入常规循环. 常规循环阶段: 蒸发出来的氦气不断被干泵抽离, 并通过逆流换热器对泵送回来经冷头再次液化的液体进一步预冷, 形成闭式循环. 当液氦累积到一定程度, 液面上方的饱和蒸气压维持平衡, 液面也维持稳定, 1 K冷盘达到最低温度, 此时关闭V1. 收气回温阶段: 当系统停机时, 关闭冷头压缩机, 关闭V1, V6, 打开V2, V3和V5, 先将氦气收回气罐, 打开手阀V10可以加速收气, 打开加热器可以加速升温. 当温度高于100 K时, 可以往真空腔里通入少量的氦气, 通过气体对流换热可迅速将室温热量传递进系统内部, 有助于加速升温. 当所有温度高于100 K时, 关闭所有阀门, 关闭干泵.

      3)测控系统. 自主开发了一套集数据采集和控制于一体的测控板卡. 该板卡涉及压力传感器、流量计、真空计等多种设备的数据采集, 并具备对所有电磁阀门开关的精确控制. 此外, 结合Lakeshore温控仪, 实现了对1 K系统的温度监测与控制. 最终通过开发上位机测控软件, 实现了包括一键降温等智能功能. 系统中各级冷盘温度通过Lakeshore标定过的CX-1030-CU-HT-1.4 L温度计进行测量, 不确定度为5 mK.

      4)真空绝热系统. 主要包括真空罩、高真空分子泵、真空管路、防辐射屏和多层绝热材料等, 保证系统内的真空度, 并阻挡外界环境的辐射漏热.

    • 图5给出了发展的脉管制冷机的典型制冷性能, 该样机使用国产压缩机和自研旋转阀, 整机具备阀分离和电隔离特性. 室温开机40 min内二级冷头温度即可降至3 K以下, 最终二级冷头无负荷温度为2.14 K, 一级冷头无负荷温度为28.8 K, 并可同时在4.2 K提供1.5 W的制冷量和在45 K提供45 W的制冷量, 性能超越美国Cryomech PT415 RM和日本住友RP-182B2S产品, 具体见表1.

    • 基于上述脉管制冷机作为预冷级搭建的1 K低温系统典型制冷性能如图6所示. 该系统4 K冷盘直径超过360 mm, 1 K冷盘直径超过220 mm, 具有足够大的冷却空间. 室温开机12 h内开始液化, 最终1 K冷盘的最低制冷温度为1.1 K, 此时一级冷盘的温度为41 K, 二级冷盘的温度为2.7 K. 该样机可在1.5 K提供67.5 mW的制冷量或者在1.6 K提供100 mW的制冷量, 单冲程最低温度可达1.01 K, 并可维持大于0.5 h.

      图7给出了系统制冷量与循环回路工质氦气的摩尔流量比值随1 K冷盘温度的变化曲线. 可以看出, 随着1 K冷盘温度由1.1 K升高至1.7 K, 系统制冷量与工质摩尔流量的比值从0逐渐增大至0.088 mW/(μmol·s–1), 基本符合(3)式的变化规律. 这主要是因为1 K冷盘温度越高, 对应的饱和蒸气压也越高, 更多的气体被干泵抽离, 致使泵后压力增大. 在系统流阻保持不变下, 泵前后压差增大, 泵送的流量也会增大.

    • 本文设计并研制了一台基于脉管制冷机预冷的1 K大冷量低温系统. 全国产化的脉管制冷机最低获取了2.14 K的制冷温度, 并可同时在4.2 K提供1.5 W的制冷量和在45 K提供45 W的制冷量, 性能超越美国Cryomech PT415 RM和日本住友RP-182B2S产品. 将其作为预冷制冷机, 搭建的1 K低温系统4 K冷盘直径超过360 mm, 1 K冷盘直径超过220 mm, 最低获取了1.1 K的制冷温度, 并可在1.6 K提供100 mW的制冷量. 单冲程最低温度可达1.01 K, 并可维持大于0.5 h. 后续将进一步提高4 K脉管制冷机的性能, 将该1 K系统搭载上9 T超导磁体, 并将进一步开展更大冷量mK稀释制冷技术的研究.

    参考文献 (27)

目录

/

返回文章
返回