搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

环形磁场金属等离子体源冷却流场的数值模拟与优化

陈国华 石科军 储进科 吴昊 周池楼 肖舒

引用本文:
Citation:

环形磁场金属等离子体源冷却流场的数值模拟与优化

陈国华, 石科军, 储进科, 吴昊, 周池楼, 肖舒

Numerical simulation and optimization of cooling flow field of cylindrical cathode with annular magnetic field

Chen Guo-Hua, Shi Ke-Jun, Chu Jin-Ke, Wu Hao, Zhou Chi-Lou, Xiao Shu
PDF
HTML
导出引用
  • 环形磁场金属等离子体源作为一种全新的等离子体源结构, 可用于产生高度离化、无大颗粒、高密度的离子束流, 但传统流道结构不能保证其高效、均匀散热, 大功率工作时可能引起密封胶圈的烧蚀失效, 需对其冷却流场进行优化设计. 利用Solidworks Flow Simulation软件对等离子体源冷却流道进行模拟, 分析出入水孔分布角度、孔数、孔径以及入水孔高度对冷却效果的影响规律, 并对流道结构参数进行优化. 结果表明, 增大水孔的周向分布范围, 有利于提高散热的均匀性; 入水孔设置在结构上层有利于减少冷却水的温度分层现象, 使铜套和密封胶圈都处于较好的冷却状态; 适当减小孔径有利于增大冷却水射流速度, 增大湍流程度强化传热, 提高换热效率. 优化后的流场结构可以提高冷却水的利用率, 在相同流量条件下获得更好的冷却效果, 改善等离子体源的放电稳定性, 为环形磁场金属等离子体源的冷却结构设计提供理论依据.
    In order to solve the problems of unstable discharge, low deposition rate and large difference in ionization rate between different targets in high power impulse magnetron sputtering, a novel cylindrical cathode with annular magnetic field based on hollow cathode effect is proposed, which can be used to produce ion beam with high ionization rate, high plasma density and no large particles. However, the traditional channel structure could not guarantee its high efficiency and uniform heat dissipation. The sealing ring may be damaged by ablation due to high power density, which restricts the further improvement of power density. Therefore, it is necessary to optimize the design of the channel structure. SolidWorks flow simulation software is used to simulate the cooling channel of plasma source. The influence of water hole structure parameters on cooling effect is analyzed, including distribution angle, hole number, diameter and inlet hole height. And the channel structure parameters are optimized. The results show that the increasing of the circumferential distribution range of the water hole is beneficial to the uniformity of heat dissipation, ensuring a large temperature difference between cooling water and copper sleeve, and strengthening heat exchange. The water inlet hole set in the upper layer of the structure is conducive to alleviating the temperature stratification phenomenon of the cooling water, so that the copper sleeve and sealing ring are in good cooling condition. Appropriately reducing the aperture is beneficial to increasing the cooling water jet velocity, enhancing the jet impact effect, and then increasing the turbulence degree, strengthening the heat transfer and improving the heat transfer efficiency. By systematically studying the influencing factors, the optimized cooling flow field structure of cylindrical cathode with an annular magnetic field is obtained. The distribution angle is 30°, the number of holes is 6, the aperture is 4 mm, and the height of water inlet hole is 36 mm. The optimized channel structure can improve the utilization rate of cooling water, obtaining better cooling effect at the same flow rate, and improving the discharge stability of the plasma source, which provides a basis for designing the cooling structure of the cylindrical cathode with annular magnetic field.
      通信作者: 吴昊, wuhao@scut.edu.cn ; 肖舒, xiaos@scut.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 52005187, 51705157, 51905177)、广东省基础与应用基础研究基金(批准号: 2019A1515110065)、中国博士后科学基金(批准号: 2019M662909)和中央高校基本科研业务费专项资金(批准号: 2019MS063)资助的课题
      Corresponding author: Wu Hao, wuhao@scut.edu.cn ; Xiao Shu, xiaos@scut.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 52005187, 51705157, 51905177), the Basic and Applied Basic Research Foundation of Guangdong Province, China (Grant No. 2019A1515110065), the China Postdoctoral Science Foundation (Grant No. 2019M662909), and the Fundamental Research Fund for the Central Universities, China (Grant No. 2019MS063)
    [1]

    D'Avico L, Beltrami R, Lecis N, Trasatti S P 2019 Coatings 9 7

    [2]

    Maksakova O, Simoẽs S, Pogrebnjak A, Bondar O, Kravchenko Y O, Koltunowicz T, Shaimardanov Z K 2019 J. Alloys Compd. 776 679Google Scholar

    [3]

    Vu T D, Chen Z, Zeng X, Jiang M, Liu S, Gao Y, Long Y 2019 J. Alloys Compd. C 7 2121

    [4]

    Kouznetsov V, Macak K, Schneider J M, Helmersson U, Petrov I 1999 Surf. Coat. Technol. 122 290Google Scholar

    [5]

    吴忠振, 田修波, 李春伟, Ricky K Y Fu, 潘锋, 朱剑豪 2014 物理学报 63 175201Google Scholar

    Wu Z Z, Tian X B, Li C W, Fu R K Y, Pan F, Chu P K 2014 Acta Phys. Sin. 63 175201Google Scholar

    [6]

    吴忠振, 田修波, 潘锋, Ricky K Y Fu, 朱剑豪 2014 物理学报 63 185207Google Scholar

    Wu Z Z, Tian X B, Pan F, Fu R K Y, Chu P K 2014 Acta Phys. Sin. 63 185207Google Scholar

    [7]

    Biswas B, Purandare Y, Sugumaran A, Khan I, Hovsepian P E 2018 Surf. Coat. Technol. 336 84Google Scholar

    [8]

    Biswas B, Purandare Y, Khan I, Hovsepian P E 2018 Surf. Coat. Technol. 344 383Google Scholar

    [9]

    Vitelaru C, Aijaz A, Parau A C, Kiss A E, Sobetkii A, Kubart T 2018 J. Phys. D: Appl. Phys. 51 165201Google Scholar

    [10]

    Meier S M, Hecimovic A, Tsankov T V, Luggenhölscher D, Czarnetzki U 2018 Plasma Sources Sci. Technol. 27 035006Google Scholar

    [11]

    Diyatmika W, Liang F K, Lou B S, Lu J H, Sun D E, Lee J W 2018 Surf. Coat. Technol. 352 680Google Scholar

    [12]

    Ganesan R, Akhavan B, Dong X, McKenzie D, Bilek M 2018 Surf. Coat. Technol. 352 671Google Scholar

    [13]

    Rudolph M, Brenning N, Raadu M A, Hajihoseini H, Gudmundsson J T, Anders A, Lundin D 2020 Plasma Sources Sci. Technol. 29 05LT01Google Scholar

    [14]

    吴忠振, 田修波, 潘锋, Ricky K Y Fu, 朱剑豪 2014 金属学报 10 1279Google Scholar

    Wu Z Z, Tian X B, Pan F, Fu R K Y, Chu P K 2014 Acta Meta. Sin. 10 1279Google Scholar

    [15]

    肖舒, 吴忠振, 崔岁寒, 刘亮亮, 郑博聪, 林海, 傅劲裕, 田修波, 潘锋, 朱剑豪 2016 物理学报 65 185202Google Scholar

    Xiao S, Wu Z Z, Cui S H, Liu L L, Zheng B C, Lin H, Fu J Y, Tian X B, Pan F, Chu P K 2016 Acta Phys. Sin. 65 185202Google Scholar

    [16]

    崔岁寒, 吴忠振, , 肖舒, 刘亮亮, 郑博聪, 林海, 傅劲裕, 田修波, 朱剑豪, 谭文长, 潘锋 2017 物理学报 66 095203Google Scholar

    Cui S H, Wu Z Z, Xiao S, Liu L L, Zheng B C, Lin H, Fu J Y, Tian X B, Chu P K, Tan W C, Pan F 2017 Acta Phys. Sin. 66 095203Google Scholar

    [17]

    崔岁寒, 吴忠振, , 肖舒, 陈磊, 李体军, 刘亮亮, 傅劲裕, 田修波, 朱剑豪, 谭文长 2019 物理学报 68 195204Google Scholar

    Cui S H, Wu Z Z, Xiao S, Chen L, Li T J, Liu L L, Fu R K Y, Tian X B, Chu P K, Tan W C 2019 Acta Phys. Sin. 68 195204Google Scholar

    [18]

    Cui S H, Wu Z Z, Xiao S, Zheng B C, Chen L, Li T J, Fu R K Y, Chu P K, Tian X B, Tan W C Fang D N, Pan F 2020 J. Appl. Phys. 127 023301Google Scholar

    [19]

    Marin F, de Miranda J R, de Souza A F 2018 Polym. Eng. Sci. 58 552Google Scholar

    [20]

    Liu X, Xu X, Liu C, Bai W, Dang C 2018 Energy 147 1Google Scholar

    [21]

    曹珍恩, 巩春志, 田修波, 杨士勤 2008 真空 2 70Google Scholar

    Cao Z E, Gong C Z, Tian X B, Yang S Q 2008 Vacuum 2 70Google Scholar

    [22]

    Bräuer G, Szyszka B, Vergöhl M, Bandorf R 2010 Vacuum 84 1354Google Scholar

  • 图 1  环形磁场金属等离子体源简化模型

    Fig. 1.  Simplified model of cylindrical cathode with annular magnetic field.

    图 2  计算模型的网格划分

    Fig. 2.  Grid of computing mode.

    图 3  不同功率下等离子体源的放电情况

    Fig. 3.  Discharge of plasma source at different power.

    图 4  出口冷却水温升的实验值与模拟值对比

    Fig. 4.  Comparison of experimental and simulated water temperature rise.

    图 5  不同孔分布角度的出水孔高度(40 mm)冷却水温度切面

    Fig. 5.  Water temperature section with outlet height (40 mm) of different hole distribution angles.

    图 6  不同出入水孔数量的铜套表面温度

    Fig. 6.  Surface temperature of copper sleeve with different number of holes.

    图 7  不同水压下出入水孔孔径对胶圈温度的影响

    Fig. 7.  Influence of hole diameter on apron temperature under different water pressure.

    图 8  出入水孔孔径对入水孔射流速度的影响

    Fig. 8.  Influence of hole diameter on water jet velocity at water inlet.

    图 9  入水孔高度对冷却效果的影响 (a)胶圈温度; (b)出水孔侧铜套温度

    Fig. 9.  Influence of hole height on cooling effect: (a) Apron temperature; (b) copper sleeve temperature of outlet side.

    图 10  入水孔高度分别为6与36 mm时的冷却水温度云图

    Fig. 10.  Temperature nephogram of cooling water when h is 6 and 36 mm.

    表 1  模型中各零部件材料物性参数

    Table 1.  Material property parameters of each component in the model.

    材料密度ρ/kg·m–3比热J/kg·K热导率W/m·K部件
    纯铁787444780磁短路环
    紫铜8960384401铜套、铜靶
    软钢787047251.9上板、上盖板、
    周向套、底板
    钕铁硼740050232.2环形磁铁
    氟橡胶125017900.2胶圈
    100041820.58冷却水
    下载: 导出CSV

    表 2  不同出入水孔分布角度、孔数时的胶圈最高温度 (℃)

    Table 2.  Maximum temperature of aprons with different distribution angles and number of holes (℃).

    角度/(°)51015202530
    数量/个
    160.4860.4860.4860.4860.4860.48
    259.8359.5159.9359.2058.2759.63
    357.4258.7957.7557.0457.6255.84
    457.4057.1855.5253.9553.4053.24
    557.1556.1353.4049.7551.2650.47
    655.9354.8253.0752.7549.8746.50
    下载: 导出CSV
  • [1]

    D'Avico L, Beltrami R, Lecis N, Trasatti S P 2019 Coatings 9 7

    [2]

    Maksakova O, Simoẽs S, Pogrebnjak A, Bondar O, Kravchenko Y O, Koltunowicz T, Shaimardanov Z K 2019 J. Alloys Compd. 776 679Google Scholar

    [3]

    Vu T D, Chen Z, Zeng X, Jiang M, Liu S, Gao Y, Long Y 2019 J. Alloys Compd. C 7 2121

    [4]

    Kouznetsov V, Macak K, Schneider J M, Helmersson U, Petrov I 1999 Surf. Coat. Technol. 122 290Google Scholar

    [5]

    吴忠振, 田修波, 李春伟, Ricky K Y Fu, 潘锋, 朱剑豪 2014 物理学报 63 175201Google Scholar

    Wu Z Z, Tian X B, Li C W, Fu R K Y, Pan F, Chu P K 2014 Acta Phys. Sin. 63 175201Google Scholar

    [6]

    吴忠振, 田修波, 潘锋, Ricky K Y Fu, 朱剑豪 2014 物理学报 63 185207Google Scholar

    Wu Z Z, Tian X B, Pan F, Fu R K Y, Chu P K 2014 Acta Phys. Sin. 63 185207Google Scholar

    [7]

    Biswas B, Purandare Y, Sugumaran A, Khan I, Hovsepian P E 2018 Surf. Coat. Technol. 336 84Google Scholar

    [8]

    Biswas B, Purandare Y, Khan I, Hovsepian P E 2018 Surf. Coat. Technol. 344 383Google Scholar

    [9]

    Vitelaru C, Aijaz A, Parau A C, Kiss A E, Sobetkii A, Kubart T 2018 J. Phys. D: Appl. Phys. 51 165201Google Scholar

    [10]

    Meier S M, Hecimovic A, Tsankov T V, Luggenhölscher D, Czarnetzki U 2018 Plasma Sources Sci. Technol. 27 035006Google Scholar

    [11]

    Diyatmika W, Liang F K, Lou B S, Lu J H, Sun D E, Lee J W 2018 Surf. Coat. Technol. 352 680Google Scholar

    [12]

    Ganesan R, Akhavan B, Dong X, McKenzie D, Bilek M 2018 Surf. Coat. Technol. 352 671Google Scholar

    [13]

    Rudolph M, Brenning N, Raadu M A, Hajihoseini H, Gudmundsson J T, Anders A, Lundin D 2020 Plasma Sources Sci. Technol. 29 05LT01Google Scholar

    [14]

    吴忠振, 田修波, 潘锋, Ricky K Y Fu, 朱剑豪 2014 金属学报 10 1279Google Scholar

    Wu Z Z, Tian X B, Pan F, Fu R K Y, Chu P K 2014 Acta Meta. Sin. 10 1279Google Scholar

    [15]

    肖舒, 吴忠振, 崔岁寒, 刘亮亮, 郑博聪, 林海, 傅劲裕, 田修波, 潘锋, 朱剑豪 2016 物理学报 65 185202Google Scholar

    Xiao S, Wu Z Z, Cui S H, Liu L L, Zheng B C, Lin H, Fu J Y, Tian X B, Pan F, Chu P K 2016 Acta Phys. Sin. 65 185202Google Scholar

    [16]

    崔岁寒, 吴忠振, , 肖舒, 刘亮亮, 郑博聪, 林海, 傅劲裕, 田修波, 朱剑豪, 谭文长, 潘锋 2017 物理学报 66 095203Google Scholar

    Cui S H, Wu Z Z, Xiao S, Liu L L, Zheng B C, Lin H, Fu J Y, Tian X B, Chu P K, Tan W C, Pan F 2017 Acta Phys. Sin. 66 095203Google Scholar

    [17]

    崔岁寒, 吴忠振, , 肖舒, 陈磊, 李体军, 刘亮亮, 傅劲裕, 田修波, 朱剑豪, 谭文长 2019 物理学报 68 195204Google Scholar

    Cui S H, Wu Z Z, Xiao S, Chen L, Li T J, Liu L L, Fu R K Y, Tian X B, Chu P K, Tan W C 2019 Acta Phys. Sin. 68 195204Google Scholar

    [18]

    Cui S H, Wu Z Z, Xiao S, Zheng B C, Chen L, Li T J, Fu R K Y, Chu P K, Tian X B, Tan W C Fang D N, Pan F 2020 J. Appl. Phys. 127 023301Google Scholar

    [19]

    Marin F, de Miranda J R, de Souza A F 2018 Polym. Eng. Sci. 58 552Google Scholar

    [20]

    Liu X, Xu X, Liu C, Bai W, Dang C 2018 Energy 147 1Google Scholar

    [21]

    曹珍恩, 巩春志, 田修波, 杨士勤 2008 真空 2 70Google Scholar

    Cao Z E, Gong C Z, Tian X B, Yang S Q 2008 Vacuum 2 70Google Scholar

    [22]

    Bräuer G, Szyszka B, Vergöhl M, Bandorf R 2010 Vacuum 84 1354Google Scholar

  • [1] 宋利伟, 石颖, 陈树民, 柯璇, 侯晓慧, 刘志奇. 地下黏弹性介质波动方程及波场数值模拟. 物理学报, 2021, 70(14): 149102. doi: 10.7498/aps.70.20210005
    [2] 邹雄, 漆小波, 张涛先, 高章帆, 黄卫星. 惯性约束聚变靶丸内杂质气体抽空流洗过程的数值模拟. 物理学报, 2021, 70(7): 075207. doi: 10.7498/aps.70.20201491
    [3] 牛越, 包为民, 李小平, 刘彦明, 刘东林. 大功率热平衡感应耦合等离子体数值模拟及实验研究. 物理学报, 2021, 70(9): 095204. doi: 10.7498/aps.70.20201610
    [4] 叶欣, 单彦广. 疏水表面振动液滴模态演化与流场结构的数值模拟. 物理学报, 2021, 70(14): 144701. doi: 10.7498/aps.70.20210161
    [5] 左娟莉, 杨泓, 魏炳乾, 侯精明, 张凯. 气力提升系统气液两相流数值模拟分析. 物理学报, 2020, 69(6): 064705. doi: 10.7498/aps.69.20191755
    [6] 姜春华, 赵正予. 化学复合率对激发赤道等离子体泡影响的数值模拟. 物理学报, 2019, 68(19): 199401. doi: 10.7498/aps.68.20190173
    [7] 喻明浩. 非平衡感应耦合等离子体流场与电磁场作用机理的数值模拟. 物理学报, 2019, 68(18): 185202. doi: 10.7498/aps.68.20190865
    [8] 危卫, 张力元, 顾兆林. 工业中粉体颗粒的荷电机理及数值模拟方法. 物理学报, 2015, 64(16): 168301. doi: 10.7498/aps.64.168301
    [9] 徐肖肖, 吴杨杨, 刘朝, 王开正, 叶建. 水平螺旋管内超临界CO2冷却换热的数值模拟. 物理学报, 2015, 64(5): 054401. doi: 10.7498/aps.64.054401
    [10] 高启, 张传飞, 周林, 李正宏, 吴泽清, 雷雨, 章春来, 祖小涛. Z箍缩Al等离子体X特征辐射谱线数值模拟及考虑叠加效应后的修正. 物理学报, 2014, 63(12): 125202. doi: 10.7498/aps.63.125202
    [11] 靳冬欢, 刘文广, 陈星, 陆启生, 赵伊君. 三股互击式喷注器及燃烧室流场的数值模拟. 物理学报, 2012, 61(6): 064206. doi: 10.7498/aps.61.064206
    [12] 蔡利兵, 王建国, 朱湘琴. 强直流场介质表面次级电子倍增效应的数值模拟研究. 物理学报, 2011, 60(8): 085101. doi: 10.7498/aps.60.085101
    [13] 庞学霞, 邓泽超, 贾鹏英, 梁伟华. 大气等离子体中氮氧化物粒子行为的数值模拟. 物理学报, 2011, 60(12): 125201. doi: 10.7498/aps.60.125201
    [14] 陈华俊, 米贤武. 强耦合光机械腔中的简正模式分裂和冷却. 物理学报, 2011, 60(12): 124206. doi: 10.7498/aps.60.124206
    [15] 邓峰, 赵正予, 石润, 张援农. 中低纬电离层加热大尺度场向不均匀体的二维数值模拟. 物理学报, 2009, 58(10): 7382-7391. doi: 10.7498/aps.58.7382
    [16] 朱昌盛, 王智平, 荆 涛, 肖荣振. 二元合金微观偏析的相场法数值模拟. 物理学报, 2006, 55(3): 1502-1507. doi: 10.7498/aps.55.1502
    [17] 郭文琼, 周晓军, 张雄军, 隋 展, 吴登生. 等离子体电极普克尔盒电光开关单脉冲过程数值模拟. 物理学报, 2006, 55(7): 3519-3523. doi: 10.7498/aps.55.3519
    [18] 袁行球, 李 辉, 赵太泽, 王 飞, 郭文康, 须 平. 超音速等离子体炬的数值模拟. 物理学报, 2004, 53(3): 788-792. doi: 10.7498/aps.53.788
    [19] 訾炳涛, 姚可夫, 许光明, 崔建忠. 脉冲磁场下金属熔体凝固流场的数值模拟. 物理学报, 2003, 52(1): 115-119. doi: 10.7498/aps.52.115
    [20] 于艳梅, 杨根仓, 赵达文, 吕衣礼, A. KARMA, C. BECKERMANN. 过冷熔体中枝晶生长的相场法数值模拟. 物理学报, 2001, 50(12): 2423-2428. doi: 10.7498/aps.50.2423
计量
  • 文章访问数:  1521
  • PDF下载量:  44
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-08-20
  • 修回日期:  2020-10-26
  • 上网日期:  2021-03-19
  • 刊出日期:  2021-04-05

环形磁场金属等离子体源冷却流场的数值模拟与优化

  • 1. 华南理工大学机械与汽车工程学院, 广州 510641
  • 2. 香港城市大学物理与材料科学系, 香港 999077
  • 3. 寰球工程项目管理公司(北京)有限公司, 北京 100029
  • 通信作者: 吴昊, wuhao@scut.edu.cn ; 肖舒, xiaos@scut.edu.cn
    基金项目: 国家自然科学基金(批准号: 52005187, 51705157, 51905177)、广东省基础与应用基础研究基金(批准号: 2019A1515110065)、中国博士后科学基金(批准号: 2019M662909)和中央高校基本科研业务费专项资金(批准号: 2019MS063)资助的课题

摘要: 环形磁场金属等离子体源作为一种全新的等离子体源结构, 可用于产生高度离化、无大颗粒、高密度的离子束流, 但传统流道结构不能保证其高效、均匀散热, 大功率工作时可能引起密封胶圈的烧蚀失效, 需对其冷却流场进行优化设计. 利用Solidworks Flow Simulation软件对等离子体源冷却流道进行模拟, 分析出入水孔分布角度、孔数、孔径以及入水孔高度对冷却效果的影响规律, 并对流道结构参数进行优化. 结果表明, 增大水孔的周向分布范围, 有利于提高散热的均匀性; 入水孔设置在结构上层有利于减少冷却水的温度分层现象, 使铜套和密封胶圈都处于较好的冷却状态; 适当减小孔径有利于增大冷却水射流速度, 增大湍流程度强化传热, 提高换热效率. 优化后的流场结构可以提高冷却水的利用率, 在相同流量条件下获得更好的冷却效果, 改善等离子体源的放电稳定性, 为环形磁场金属等离子体源的冷却结构设计提供理论依据.

English Abstract

    • 目前, 等离子体镀膜技术在表面处理领域得到广泛应用[1-3]. 针对磁控溅射粒子离化率不高、可控性差、阴极弧存在“大颗粒”等问题, Kouznetsov等[4]在1999年提出了高功率脉冲磁控溅射(high power impulse magnetron sputtering, HiPIMS)技术, 离化率可达到60%—90%[5,6], 且束流中不含“液滴”[7,8], 但其放电不稳定[9,10]、沉积速率低[11-13], 并且对于不同靶材料, 其离化率差异显著[14]. 针对HiPIMS存在的不足, 吴忠振课题组[15]提出了一种基于空心阴极效应的环形磁场金属等离子体源结构, 可将溅射约束在筒形靶材内部, 溅射粒子在腔内反复碰撞、离化, 有效提高离化率. 通过引出栅将腔内离子引出获得完全离化的束流实现工件表面的沉积, 而将“打弧”产生的“金属液滴”留在腔内. 同时引出的束流离开了靶电压鞘层, 不容易被回吸到靶面, 可提高薄膜沉积速率.

      前期研究已对这种环形磁场金属等离子体源结构进行设计优化. 针对磁场分布不均导致的靶面刻蚀不均和电子快速逃逸问题, 通过在等离子体源两端添加电子阻挡屏极和磁场补偿, 提高等离子体放电强度、范围和靶材刻蚀均匀性[16]. 然而, 这种装置仅通过浓度差扩散实现等离子体向基体表面的输送沉积, 输送过程中离子损失较大, 薄膜的沉积速率没有明显改善. 为此, 利用外加扩展电磁场对离子的运动方向加以约束, 实现等离子体的聚焦和可控输运, 从而减少离子损失, 提高沉积效率[17]. 最后, 对环形磁场金属等离子体源中HiPIMS的放电过程和粒子演化进行了系统研究分析[18].

      以上研究表明环形磁场金属等离子体源能够有效改善HiPIMS技术存在的缺陷, 实现薄膜的高质量高效率制备. 然而, 等离子体源在放电过程中会产生大量的热量, 在高功率密度条件下工作若不能及时散热冷却将会对等离子体源内部的密封胶圈等结构件造成烧蚀损坏, 制约了等离子体源功率密度的进一步提高, 因此提高其冷却效率尤为重要. 本文对其冷却系统流场结构进行设计优化, 借鉴随形冷却技术和模具型腔冷却结构的设计思路[19,20], 将环形磁铁与铜套等腔体结构都浸泡在冷却水中进行直接换热冷却, 通过Solidworks Flow Simulation软件建立等离子体源冷却流场的三维数值模型, 对其短路环上入水孔分布角度、数量、孔径以及高度进行系统研究, 分析不同结构参数对等离子体源内部密封胶圈温度与铜套温度的影响, 为环形磁铁金属等离子体源的冷却结构设计提供理论依据, 以实现在相同冷却水条件下获得更高的冷却效果, 改善等离子体源的放电稳定性, 为进一步提高其功率密度提供硬件支撑.

    • 环形磁场金属等离子体源简化模型如图1所示, 在短路环上设置出入水孔, 入水孔与出水孔数量一致且沿周向对称分布, 冷却水由进水管进入底板凹槽并沿着周向套竖直矩形槽分成数股经入水孔进入等离子体源腔内, 最后从另一侧的出水孔流出. 其中, h为入水孔高度, d为出入水孔直径, α为出入水孔分布角度. 为了获得较高的溅射效率, 提高束流密度, 溅射靶面平均功率达到50 W/cm2, 其中约60%转化成热量并通过铜套传导至等离子体源内部, 通过冷却系统进行散热. 为了保证等离子体源的正常工作, 需要保证: 结构件不被熔化(< 1000 ℃), 磁体不能退磁(< 300 ℃), 密封胶圈不被烧蚀($ \leqslant $ 200 ℃), 筒形铜套受热较均匀. 在输入热量、冷却水条件相同的情况下, 通过改变孔数量、孔分布角度提高等离子体源的散热性能, 在此基础上再对孔径、入水口高度进行优化, 进一步降低胶圈温度和提高温度均匀性.

      图  1  环形磁场金属等离子体源简化模型

      Figure 1.  Simplified model of cylindrical cathode with annular magnetic field.

    • 在不影响等离子体源功能结构以及冷却效果的前提下, 作以下假设:

      1)因流道需要在短路环上开孔, 认为孔数量、孔径、孔高度、孔分布角不会对等离子体源磁场产生影响.

      2)采用水作为冷却介质, 为不可压缩的牛顿流体, 定常流动.

      3)溅射发生在靶材内表面, 溅射区域呈带状分布, 占靶面64%左右, 其导热部件是紫铜, 导热性良好, 因此假设热源为均匀面热源.

      流体的流动和换热满足连续性、质量守恒定律、能量守恒定律和动量定律, 采用的控制方程为

      $ \frac{{\partial \left( {\rho \phi } \right)}}{{\partial t}} + \nabla \cdot \left( {\rho {{\mathit{\boldsymbol{U}}}}\phi } \right) = \nabla \cdot \left( {{\varGamma _\phi }\nabla \phi } \right) + {S_\phi }, $

      式中, $ \phi $为通用变量; $ {S}_{\phi} $为广义源项; $ \varGamma_{\phi} $为广义扩散系数, 所述广义是数值计算模型方程中的一种定义, 不同求解变量间的区别除了边界条件与初始条件以外, 就在于$ \varGamma_{\phi} $$ {S}_{\phi} $的表达式不同; 其中$ \rho $为密度; t为时间; $\nabla $表示梯度; $\nabla \cdot$表示散度; U为速度矢量.

      热功耗均匀分布在铜套内表面, 从铜套内表面至外表面的热传导符合单层圆筒壁传热模型:

      $ Q = 2\pi L\lambda \frac{{{t_1} - {t_2}}}{{{\rm{\ln}}({{{r_2}}}/{{{r_1}}}})}, $

      式中, Q为热量; L为铜套高度; λ为铜套的导热系数; t1t2分别为铜套内、外表面的温度; r1r2分别为铜套内、外表面的半径.

      热量传导到铜套外表面后与腔内的冷却水发生对流换热, 冷却水吸收的热量为

      $ Q = \alpha A\Delta {T_{\rm{m}}}, $

      式中, α为平均对流传热系数; A为有效传热面积; ΔTm为传热平均温差.

    • 对模型各部件材料物性参数进行定义, 如表1所列. 设置边界条件, 进水端采用压力入口条件, 如无特别说明, 入口压力都设置为5.0 × 105 Pa水压, 初始温度为20 ℃, 出水端采用压力出口条件, 设置为环境压力(1.013 × 105 Pa). 热源选取铜套内表面, 在不同电源功率条件下设置相应的不同热源热功耗. 本文选取其中一种工况进行模拟研究, 热功耗加载为2 kW.

      材料密度ρ/kg·m–3比热J/kg·K热导率W/m·K部件
      纯铁787444780磁短路环
      紫铜8960384401铜套、铜靶
      软钢787047251.9上板、上盖板、
      周向套、底板
      钕铁硼740050232.2环形磁铁
      氟橡胶125017900.2胶圈
      100041820.58冷却水

      表 1  模型中各零部件材料物性参数

      Table 1.  Material property parameters of each component in the model.

    • 采用SolidWorks软件建立模型, 并用SolidWorks Flow Simulation对模型进行自动网格划分, 初始网格级别选择为3级, 手动定义最小缝隙尺寸为0.4 mm. 以入水孔高度32 mm、孔径4 mm、分布角度30°的六孔模型网格划分为例, 总网格数为1481963, 其中流体网格数564680, 固体网格数917283, 如图2所示.

      图  2  计算模型的网格划分

      Figure 2.  Grid of computing mode.

    • 为了验证数值模拟的正确性, 加工得到入水孔高度32 mm、孔径4 mm、分布角度30°的六孔环形磁场等离子体源, 采用直流磁控溅射在自主设计的真空系统中对其进行放电测试[15,16]. 真空腔长和宽各600 mm, 高500 mm, 腔体材料为不锈钢. 采用直流溅射电源, 溅射靶采用内径120 mm、外径132 mm、高50 mm的环形铜靶(纯度为99.9%). 将背景压力抽真空至1 × 10–3 Pa, 通入纯度为99.99%的高纯氩气, 测试不同功率下等离子体源的放电情况, 并通过红外测温仪监测3, 6, 9和12 kW这四种功率加载下出水端的冷却水温度. 在实际放电过程中, 由于开关管的饱和压降、变压器、电缆、电极及其他阴极结构件上的内阻损耗, 电源功率中只有约75%的能量供应到靶材表面作为实际的溅射功率[21], 而这部分能量中约80%会转变成为靶材的热量[22], 此时靶材热量Q约为电源功率的60%. 因此, 在实验对应的模拟中将热源热功耗设置为1.8, 3.6, 5.4 和7.2 kW.

      不同功率下等离子体源的放电辉光如图3所示, 辉光呈绿色, 采用等离子体发射光谱仪在等离子体源中央轴上探测到406.3和404.3 nm的光谱波长信号, 分别为Cu原子与Cu+的典型特征波长. 随着功率增大辉光更加强烈, 即等离子体放电强度更大、密度更高. 将不同功率下出水端冷却水温度的实验测试值与模拟值进行对比, 结果如图4所示, 图中纵坐标值表示冷却水出口温度与入口初始温度的温升差值. 结果表明模拟值与实验结果误差较小, 符合度高, 证明了数值模拟的正确性. 存在的偏差主要是由于实验在25 ℃室温条件下进行, 初始温度20 ℃的冷却水在输入过程中存在一定的温升, 且出口高于室温的冷却水存在降温的现象. 此外, 模拟结果还表明当功率超过18 kW时胶圈最高温度超过200 ℃, 发生烧蚀失效, 等离子体源将无法正常工作.

      图  3  不同功率下等离子体源的放电情况

      Figure 3.  Discharge of plasma source at different power.

      图  4  出口冷却水温升的实验值与模拟值对比

      Figure 4.  Comparison of experimental and simulated water temperature rise.

    • 根据环形磁场金属等离子体源结构, 冷却水的引入必须经过短路环, 在短路环上设置适当的孔以保证良好的冷却效果. 同时为了保证等离子体源上部的冷却, 出水孔高度设置为40 mm, 其数量、分布角度、孔径与入水孔一致. 当采用一个入水孔, 孔径为4 mm时, 其最高温度出现在出水孔附近, 密封胶圈内侧, 达到稳态时温度为60.48 ℃. 为了研究开孔分布角度、数量与冷却效果间的关系, 在短路环上按分布角度5°, 10°, 15°, 20°, 25°和30°分别开1—6个直径为4 mm的孔, 监测胶圈的最高温度, 结果如表2所列. 当分布角度为30°, 孔数量为6时对应的胶圈最高温度最低, 冷却效果最好, 此时胶圈上的稳态最高温度为46.50 ℃. 胶圈上的最高温度随着孔数量的增加呈逐渐下降的趋势. 而随着分布角度的增加, 1—3个水孔时的胶圈最高温度没有表现出明显的规律性, 但4—6个水孔时呈现下降的趋势.

      角度/(°)51015202530
      数量/个
      160.4860.4860.4860.4860.4860.48
      259.8359.5159.9359.2058.2759.63
      357.4258.7957.7557.0457.6255.84
      457.4057.1855.5253.9553.4053.24
      557.1556.1353.4049.7551.2650.47
      655.9354.8253.0752.7549.8746.50

      表 2  不同出入水孔分布角度、孔数时的胶圈最高温度 (℃)

      Table 2.  Maximum temperature of aprons with different distribution angles and number of holes (℃).

      针对表2结果, 结合不同孔分布角度时的冷却水温度切面进行分析, 如图5所示, 随着分布角度的增大, 整体水温趋向均匀, 局部的高温区域温度降低且面积逐渐减少直至消失. 局部高温会减弱该区域的换热效果, 而初始冷却水与铜套表面的较大温差有利于强化传热, 因此当入水孔分布跨度更广时, 铜套表面接触初始冷却水的范围更广, 换热效果更好, 且冷却更加均匀. 部分冷却水从最两端的出水孔流出, 在腔体中的停留时间短, 水温相对较低, 因此换热充分、温度较高的冷却水不会在中央出水孔处大量富集, 避免了局部高温的产生.

      图  5  不同孔分布角度的出水孔高度(40 mm)冷却水温度切面

      Figure 5.  Water temperature section with outlet height (40 mm) of different hole distribution angles.

      铜套表面温度更加直观地反映了冷却水的冷却效果, 与胶圈最高温度有着相同的变化趋势, 铜套温度降低有利于减少传递给胶圈的热量, 使胶圈温度降低. 图6为分布角为30°时不同出入水孔数量时的出水侧铜套表面温度, 随着孔数的增加铜套表面温度明显降低, 且上端接近胶圈的位置降温更明显. 孔数量的增加同样增大了初始冷却水的分布范围, 使冷却更加均匀, 减少局部高温.

      图  6  不同出入水孔数量的铜套表面温度

      Figure 6.  Surface temperature of copper sleeve with different number of holes.

    • 根据前一节模拟分析得到的最优冷却效果, 采用分布角度为30°的6孔冷却结构, 保持入水孔高度32 mm、出水孔高度40 mm不变, 对孔径大小为3, 4, 5, 6和7 mm五种结构的模型进行模拟计算.

      不同入水端压力下的胶圈最高温度如图7所示, 在实际工况下(入口水压5.0 × 105 Pa), 胶圈最高温度随着入水孔孔径增大呈现先减小后增大再减小的趋势, 在孔径4 mm时实现最优的冷却效果, 胶圈稳态最高温度为46.71 ℃. 当入水端压力恒定时, 随着孔径改变冷却水流速与流量都会发生变化, 此时系统冷却效果受到流量与流速的综合影响. 由图8可见, 冷却水经过入水孔发生射流, 冲击铜套壁面, 射流速度随孔径减小而增大, 冲击面积也随之增大. 射流冲击作用是强化局部换热的有效方法, 提高射流速度有利于增强冷却水的湍流度, 提高对流传热系数, 同时加剧对壁面的冲击, 冲击面积的增大有利于扩大换热强化区域范围, 进而提升换热效果. 因此, 胶圈最高温度呈现随着孔径减小而降低的趋势.

      图  7  不同水压下出入水孔孔径对胶圈温度的影响

      Figure 7.  Influence of hole diameter on apron temperature under different water pressure.

      图  8  出入水孔孔径对入水孔射流速度的影响

      Figure 8.  Influence of hole diameter on water jet velocity at water inlet.

      随着孔径不断增大, 射流冷却水变得相对发散, 冷却射流越难直接冲击到壁面上, 局部换热减弱, 但是由于流道内的总流量增大, 沿壁面的周向流动更加剧烈, 因此整体壁面的换热效果有所提升, 胶圈最高温度开始降低. 受到结构尺寸的限制, 在一定孔径范围内, 通过减小孔径可以充分利用射流强化作用实现最优的冷却效果. 孔径的选取还需要根据实际的水压条件, 由图7所示, 在水压高于或低于5 × 105 Pa的条件下, 孔径3 mm时冷却水的射流效应仍占据主导作用, 因此可以适当进一步减小孔径.

    • 本节保持孔径4 mm、孔分布角度30°、出水孔高度40 mm不变, 对6, 10, 14, 18, 32和36 mm六种不同入水孔高度对冷却效果的影响进行模拟分析, 其中6—18 mm高度的入水孔位于下端磁铁和中间磁铁之间, 32—36 mm高度的入水孔位于中间磁铁和上端磁铁之间, 中间磁铁高度22 mm、厚6 mm, 结果如图9所示.

      图  9  入水孔高度对冷却效果的影响 (a)胶圈温度; (b)出水孔侧铜套温度

      Figure 9.  Influence of hole height on cooling effect: (a) Apron temperature; (b) copper sleeve temperature of outlet side.

      图9可见, 当入水孔高度为36 mm时冷却效果最好, 胶圈上的最高温度最低, 此时胶圈上的最高温度达到稳态时为46.2 ℃, 入水孔位于上层时的入口侧胶圈低温段温度相比于位于下层时的温度差异为ΔT > 5 ℃, 铜套表面温度分布随入水孔高度的变化趋势与胶圈温度变化一致. 结合图10冷却水温度云图, 选择最低高度6 mm与最高高度36 mm的结构进行对比分析, 高度为6 mm时冷却水上层与下层温差明显, 温度较高的水往上层流动, 使温度高的水流聚集在上层, 而温度较低的水则处于底层, 导致明显的分层现象, 冷却效果被削弱; 当高度为36 mm时, 冷却水从上部流入, 位置更接近顶部胶圈, 较低的初始温度能够很好地冷却入水口侧的铜套和胶圈, 且冷却水带动上层水向下流动并不断带走上层热水的热量, 有效避免冷却水的温度分层, 使铜套和胶圈都处于较好的冷却状态.

      图  10  入水孔高度分别为6与36 mm时的冷却水温度云图

      Figure 10.  Temperature nephogram of cooling water when h is 6 and 36 mm.

      胶圈上的最高温度随入水口高度变化可分成两部分: 高度6—18 mm的入水孔介于底端环形磁铁与中间磁铁间, 在这个区间内胶圈的最高温度先减小后增大, 高度为10 mm时胶圈最高温度为50.61 ℃, 为该区间最低温度; 高度32—40 mm的入水孔介于中间磁铁与上端磁铁间, 相比于前一部分温度显著降低. 当入水孔高度靠近中间磁铁时, 冷却水进入流道后冲刷中间磁铁下表面, 形成涡流, 对下表面边界层增加了扰动, 使边界层变薄, 换热增强, 冷却水经过中间磁铁后温度相对更高, 导致对铜套和胶圈的冷却效果下降, 因此入水孔高度高于10 mm后胶圈最高温度又开始升高.

    • 本文基于计算流体力学方法, 采用Solidworks Flow Simulation软件对环形磁场金属等离子体源冷却流场进行模拟分析, 根据不同结构下胶圈、铜套以及冷却水的温度分布, 分析了短路环上出入水孔分布角度、数量、孔径以及入水孔高度对冷却效果的影响规律, 得出以下结论:

      1)增大水孔的周向分布范围, 有利于提高散热的均匀性, 使流体均匀维持在较低温度, 保证与铜套温度间的更大温差, 强化换热冷却;

      2)入水孔应设置在结构上层, 有利于减少冷却水的温度分层现象, 使铜套和胶圈都处于良好的冷却状态;

      3)适当减小出入水孔孔径, 有利于增大冷却水射流速度, 增强冷却水的射流冲击作用, 进而增大湍流程度强化传热, 提高换热效率;

      4)通过系统地研究各影响因素, 在0.5 MPa水压条件下获得优化的环形磁场金属等离子体源冷却流场结构, 即分布角度为30°, 孔数为6孔, 孔径为4 mm, 入水孔高度为36 mm.

参考文献 (22)

目录

    /

    返回文章
    返回