Post-Arc Current Simulation in Vacuum Circuit Breakers with CuCr Contacts Based on the Particle-in-Cell Monte Carlo Collision Method
Li Xubin1,2, Jiang Yuan1,3, Wang Wen4, Liu Xiaoming5, Zou Jiyan2
1. Department of Electrical Engineering and Electronics University of Liverpool Liverpool L69 3GJ U.K.; 2. School of Electrical Engineering Dalian University of Technology Dalian 116024 China; 3. School of Automation and Electrical Engineering University of Science and Technology Beijing Beijing 100083 China; 4. China Electric Power Research Institute Beijing 100192 China; 5. State Key Laboratory of Smart Power Distribution Equipment and System Hebei University of Technology Tianjin 300401 China
摘要 深入研究电流过零后触头间隙内的剩余等离子体动态演变过程对于提升真空开关的开断性能至关重要。以往基于粒子模拟-蒙特卡洛(PIC-MCC)方法的弧后仿真多采用铜触头数据,而对实用的铜铬材料的弧后介质恢复研究较少,且缺少实验验证。该文基于文献中的弧后电流实验数据,对铜铬触头材料的真空开关弧后电流进行仿真研究。首先,建立简化的二维PIC-MCC模型,设置初始条件进行粒子模拟,将PIC-MCC仿真结果代入连续过渡模型电流方程,得到弧后电流仿真曲线;其次,分析了设置的经验原子密度引起的击穿过程,并给出改进初始原子密度的方法;然后,研究了过零时刻等离子体漂移速度对弧后电流的影响,发现弧后电流峰值随初始漂移速度的增加而增大,当初始漂移速度为4 000 m/s 时,仿真结果与实验数据吻合良好,验证了仿真模型和参数设置的合理性;最后,对比分析了纯铜及不同铬含量的铜铬材料对弧后电流的影响规律,为后续真空开关弧后介质恢复特性的研究提供理论参考。
Abstract:Understanding the post-arc dielectric recovery process in the contact gap after the current zero is essential for improving the interruption capability of vacuum circuit breakers (VCBs). While previous studies have extensively employed the particle-in-cell Monte Carlo collision (PIC-MCC) method to simulate post-arc behavior in copper contacts, few have addressed the widely used CuCr contact materials, and experimental validation remains limited. This study aims to simulate the post-arc current (PAC) in VCBs with CuCr contacts using a simplified two-dimensional PIC-MCC model, validated against experimental data. The research focuses on improving the accuracy of initial parameter settings and analyzing key factors affecting PAC, such as initial atom density and plasma drift velocity. A simplified 2D PIC-MCC model was developed to simulate the dynamic behavior of residual plasma in the contact gap after current zero. The model incorporates electrons, monovalent Cu+ and Cr+ ions, and Cu/Cr atoms uniformly distributed in the gap, with quasi-neutrality assumed initially. The composition of Cu and Cr ions/atoms was set according to the mass ratio of the contact material. Key collision processes, including ionization, excitation, elastic and charge exchange, were considered using cross-section data from established databases. The PAC was computed by integrating PIC-MCC simulation results into the continuous transition model (CTM) current equation, which accounts for both conduction and displacement currents contributed by Cu and Cr ions. Initial parameters, including plasma density, were derived from experimental conditions and post-arc charge measurements, while atom density was estimated using an analytical model based on erosion rate and thermal velocity. Other initial settings, including the recovery voltage and gap length, are the same as the experimental conditions. Simulation results revealed that empirically set high atom densities (e.g., 1×1021 m-3) led to transient breakdown due to enhanced ionization collisions, inconsistent with experimental observations. By recalibrating the initial atom density to 3×1018 m-3 using an analytical model, the simulation accurately reproduced the experimental trend of successful interruption without breakdown. Furthermore, incorporating an initial plasma drift velocity significantly improved PAC waveform accuracy. The peak value of the PAC increased and its duration shortened with higher drift velocities, and a value of 4 000 m/s yielded the best agreement with measured data. Additional simulations comparing pure Cu, CuCr25, and CuCr50 contacts showed minimal influence of Cr content on PAC characteristics under identical conditions, suggesting that material composition has a limited impact on axial plasma dynamics during the post-arc phase and simplified assumptions using pure copper can be adequate when the detailed alloy composition is unavailable. The study demonstrates that the combination of PIC-MCC simulation with CTM-based current calculation offers an effective method for PAC prediction in VCBs with CuCr contacts. Accurate estimation of initial atom density and inclusion of plasma drift velocity are critical for replicating experimental behavior. The findings provide practical insights for optimizing post-arc recovery performance and support the use of simplified copper-based models when material-specific data are unavailable. Future work should address limitations such as radial diffusion effects, multivalent ions, and surface interactions to further enhance model fidelity.
李旭彬, 蒋原, 王雯, 刘晓明, 邹积岩. 基于粒子模拟-蒙特卡洛模型的真空铜铬触头弧后电流仿真研究[J]. 电工技术学报, 2026, 41(7): 2446-2456.
Li Xubin, Jiang Yuan, Wang Wen, Liu Xiaoming, Zou Jiyan. Post-Arc Current Simulation in Vacuum Circuit Breakers with CuCr Contacts Based on the Particle-in-Cell Monte Carlo Collision Method. Transactions of China Electrotechnical Society, 2026, 41(7): 2446-2456.
[1] 程显, 白青林, 葛国伟, 等. 真空与气体一体化串联机械开关及其直流快速转移应用[J]. 电工技术学报, 2024, 39(5): 1464-1473. Cheng Xian, Bai Qinglin, Ge Guowei, et al.High speed switch with series-connected vacuum and gas interrupters and its DC fast transfer application[J]. Transactions of China Electrotechnical Society, 2024, 39(5): 1464-1473. [2] 李元钊, 丁健刚, 刘志远, 等. 并联电容对真空灭弧室内部电位分布的改善作用[J]. 电工技术学报, 2025, 40(2): 610-624. Li Yuanzhao, Ding Jiangang, Liu Zhiyuan, et al.The improving effect of parallel capacitors on the internal potential distribution of vacuum interrupters[J]. Transactions of China Electrotechnical Society, 2025, 40(2): 610-624. [3] 葛国伟, 王文博, 程显, 等. 基于两间隙异步联动的一体化高压真空灭弧室电场设计[J]. 电工技术学报, 2024, 39(17): 5555-5564. Ge Guowei, Wang Wenbo, Cheng Xian, et al.Electric field design of integrated high-voltage vacuum interrupter based on two-gap asynchronous linkage[J]. Transactions of China Electrotechnical Society, 2024, 39(17): 5555-5564. [4] 马飞越, 姚晓飞, 刘志远, 等. 2/3匝线圈式纵磁触头大开距真空断路器分闸速度设计[J]. 电工技术学报, 2024, 39(13): 4139-4152. Ma Feiyue, Yao Xiaofei, Liu Zhiyuan, et al.Design of opening velocities for large-gap vacuum circuit breakers with 2/3 coil-type axial magnetic field contacts[J]. Transactions of China Electrotechnical Society, 2024, 39(13): 4139-4152. [5] 张登奎, 张立岩, 李志兵, 等. 基于经验小波变换的真空断路器关合预击穿燃弧时间分析[J]. 电工技术学报, 2024, 39(增刊1): 106-116. Zhang Dengkui, Zhang Liyan, Li Zhibing, et al.Analysis of pre-strike arcing time for vacuum circuit breaker closure based on empirical wavelet transform[J]. Transactions of China Electrotechnical Society, 2024, 39(S1): 106-116. [6] Andrews J G, Varey R H.Sheath growth in a low pressure plasma[J]. The Physics of Fluids, 1971, 14(2): 339-343. [7] 刘晓明, 黄翀阳, 邹积岩. 多断口直流真空断路器介质恢复强度研究[J]. 中国电机工程学报, 2016, 36(19): 5357-5364, 5418. Liu Xiaoming, Huang Chongyang, Zou Jiyan.Investigations on dielectric recovery strength for the DC-VCB with multi-breaks[J]. Proceedings of the CSEE, 2016, 36(19): 5357-5364, 5418. [8] Birdsall C K.Particle-in-cell charged-particle simulations, plus Monte Carlo collisions with neutral atoms, PIC-MCC[J]. IEEE Transactions on Plasma Science, 1991, 19(2): 65-85. [9] Wang Zhenxing, Wang Haoran, Zhou Zhipeng, et al.Fully kinetic model of breakdown during sheath expansion after interruption of vacuum arcs[J]. Journal of Applied Physics, 2016, 120(8): 083301. [10] Shi Qilin, Wu Hao, Yuan Zhao, et al.The influence of weak transverse magnetic field on plasma dissipation process in the post-arc phase in a vacuum interrupter[J]. Plasma Science and Technology, 2022, 24(5): 055501. [11] Cheng Xian, Chen Hui, Ge Guowei, et al.Influence of the vapor shield potential control on the post-arc sheath development in vacuum interrupters[J]. Vacuum, 2024, 222: 112957. [12] 王丹, 王立军, 王智炜, 等. 金属蒸气作用下弧后等离子体扩散过程的二维粒子模拟[J]. 中国电机工程学报, 2021, 41(15): 5391-5401. Wang Dan, Wang Lijun, Wang Zhiwei, et al.Two-dimensional particle simulation of the post-arc plasma decay process in metal vapor[J]. Proceedings of the CSEE, 2021, 41(15): 5391-5401. [13] 李旭彬, 冷爽, 黄智慧, 等. 基于PIC-MCC模型的真空开关弧后鞘层仿真研究[J]. 真空科学与技术学报, 2019, 39(11): 1033-1038. Li Xubin, Leng Shuang, Huang Zhihui, et al.Post-arc sheath of vacuum circuit breaker: a simulation study in PIC-MCC method[J]. Chinese Journal of Vacuum Science and Technology, 2019, 39(11): 1033-1038. [14] 李旭彬. 真空开关弧后鞘层发展阶段的仿真研究[D]. 大连: 大连理工大学, 2019. Li Xubin.Simlation of post-arc sheath development phase in vacuum circuit breakers[D]. Dalian: Dalian University of Technology, 2019. [15] 冷爽. 基于PIC-MCC模型的真空电弧零区等离子体仿真[D]. 大连: 大连理工大学, 2021. Leng Shuang.Simulation of plasma in current zero zone of vacuum arc based on PIC-MCC model[D]. Dalian: Dalian University of Technology, 2021. [16] Kaumanns J.Measurements and modelling in the current zero region of vacuum circuit breakers for high current interruption[C]//Proceedings of 17th International Symposium on Discharges and Electrical Insulation in Vacuum, Berkeley, CA, USA, 1996: 276-280. [17] Mo Yongpeng, Shi Zongqian, Jia Shenli.Study of post-arc residual plasma dissipation process of vacuum circuit breakers based on a 2D particle-in-cell model[J]. Plasma Science and Technology, 2022, 24(4): 045401. [18] Brown I G, Godechot X.Vacuum arc ion charge-state distributions[J]. IEEE Transactions on Plasma Science, 1991, 19(5): 713-717. [19] Nikolaev A G, Oks E M, Savkin K P, et al.Charge state, angular distribution, and kinetic energy of ions from multicomponent-cathodes in vacuum arc devices[J]. Journal of Applied Physics, 2014, 116(21): 213303. [20] Wang Zhenxing, Tian Yunbo, Ma Hui, et al.Decay modes of anode surface temperature after current zero in vacuum arcs: part II: theoretical study of dielectric recovery strength[J]. IEEE Transactions on Plasma Science, 2015, 43(10): 3734-3743. [21] Tkachev A N, Fedenev A A, Yakovlenko S I.Townsend coefficient, escape curve, and fraction of runaway electrons in copper vapor[J]. Laser Physics, 2007, 17(6): 775-781. [22] Pitchford L C, Alves L L, Bartschat K, et al.LXCat: an open-access, web-based platform for data needed for modeling low temperature plasmas[J]. Plasma Processes and Polymers, 2017, 14(1/2): 1600098. [23] Cullen D E, Perkins S T, Seltzer S M.Tables and graphs of electron-interaction cross-sections from 10 eV to 100 GeV derived from the LLNL Evaluated Electron Data Library (EEDL), Z=1-100[R]. Berkeley: University of California Radiation Laboratory, 1991. [24] Mason E A, McDaniel E W. Transport Properties of Ions in Gases[M]. New York: John Wiley & Sons, 1988. [25] Jia Shenli, Mo Yongpeng, Shi Zongqian, et al.Post-arc current simulation based on measurement in vacuum circuit breaker with a one-dimensional particle-in-cell model[J]. Physics of Plasmas, 2017, 24(10): 103511. [26] Sarrailh P, Garrigues L, Hagelaar G J M, et al. Sheath expansion and plasma dynamics in the presence of electrode evaporation: Application to a vacuum circuit breaker[J]. Journal of Applied Physics, 2009, 106(5): 053305. [27] Takahashi S, Arai K, Morimiya O, et al.A PIC-MCC simulation of the high-voltage interruption ability of a vacuum interrupter[J]. IEEE Transactions on Plasma Science, 2007, 35(4): 912-919. [28] Glinkowski M T, Stoving P.Numerical modeling of vacuum arc interruption based on the simplified plasma equations[J]. IEEE Transactions on Magnetics, 1995, 31(3): 1924-1927. [29] Lins G.Measurement of the neutral copper vapor density around current zero of a 500-a vacuum arc using laser-induced fluorescence[J]. IEEE Transactions on Plasma Science, 1985, 13(6): 577-581. [30] Rich J A, Farrall G A.Vacuum arc recovery phenomena[J]. Proceedings of the IEEE, 1964, 52(11): 1293-1301. [31] Wang Zhenxing, Geng Yingsan, Liu Zhiyuan.Stepwise behavior of free recovery processes after diffused vacuum arc extinction[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2012, 19(2): 582-590. [32] Lins G.Evolution of copper vapor from the cathode of a diffuse vacuum arc[J]. IEEE Transactions on Plasma Science, 1987, 15(5): 552-556. [33] Mo Yongpeng, Shi Zongqian, Bai Zhibin, et al.Influence of residual plasma drift velocity on the post-arc sheath expansion of vacuum circuit breakers[J]. Physics of Plasmas, 2016, 23(5): 053506. [34] Anders A, Yushkov G Y.Ion flux from vacuum arc cathode spots in the absence and presence of a magnetic field[J]. Journal of Applied Physics, 2002, 91(8): 4824-4832. [35] Glinkowski M T, Greenwood A, Stoving P.Numerical simulation of the high-current vacuum arc[C]//XVI International Symposium on Discharges and Electrical Insulation in Vacuum, Moscow-St. Petersburg, Russia, 1994: 153. [36] Budde M, Horn A, Körner F, et al.Dielectric behaviour of vacuum circuit-breakers[C]//14th Inter-national Symposium on High Voltage Engineering, Beijing, China, 2005: 1-6.