NUMERICAL SIMULATION OF VACUUM ARC
UNDER DIFFERENT AXIAL MAGNETIC FIELDS
AND DIFFERENT EROSION RATES

Claire Simonnet, Pierre Freton, F Reichert, A Petchanka, Jean Jacques
Gonzalez

► To cite this version:

Claire Simonnet, Pierre Freton, F Reichert, A Petchanka, Jean Jacques Gonzalez. NUMERICAL SIMULATION OF VACUUM ARC UNDER DIFFERENT AXIAL MAGNETIC FIELDS AND DIFFERENT EROSION RATES. 2023. hal-04126550

HAL Id: hal-04126550

https://hal.science/hal-04126550

Preprint submitted on 13 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
NUMERICAL SIMULATION OF VACUUM ARC UNDER DIFFERENT AXIAL
MAGNETIC FIELDS AND DIFFERENT EROSION RATES

C. SIMONNET$^{a,*}$, P. FRETON$^a$, J.J. GONZALEZ$^a$, F. REICHERT$^b$,
A. PETCHANKA$^b$

<sup>a</sup> Laboratoire plasma et conversion d'énergie UMR 5213 - CNRS - Toulouse INP - UT3 Université Toulouse 3 - Paul Sabatier 118, route de Narbonne - bât 3R3 - 31062 Toulouse cedex 9, France

<sup>b</sup> Siemens Energy, Paulsternstr. 26, 13629, Berlin, Germany

* claire.simonnet@laplace.univ-tlse.fr

**Abstract.**

In this paper, we present a three-dimensional (3D) numerical model that describes a high current vacuum arc in the subsonic regime between copper contacts under the influence of an axial magnetic field (AMF). The model is based on a MHD approach to the inter-electrode plasma and is realized with the commercial software Ansys Fluent. Only the plasma region is modelled. The behaviour of the copper electrodes is taken into account by specific boundary conditions. We describe the model, the boundary conditions and present the simulation results. Temperature profiles and electron and ion densities will be presented in a reference case. Based on this arc modelling, a parametric study on arc characteristics will be performed, including varying the AMF strength, the ablation rate and the current intensity.

**Keywords:** plasma, numerical simulation, subsonic regime, vacuum circuit breaker.

# 1. Introduction

It is with a view to ecological transition that this study is being developed in order to reduce greenhouse gas emissions from circuit breakers and to comply with the environmental standards [1]. Siemens Energy is a central player in the field of electrical switching; the company develops high voltage circuit breakers (above 72.5 kV). In collaboration with the AEPPT team at Laplace, the company's objective is to study high-voltage vacuum arc breaking in order to develop its product range. Vacuum arc plasmas with an external axial magnetic field (AMF) have been successfully used in commercial vacuum switches for medium voltage (above 1 kV up to 52 kV included) distribution networks for several decades [2, 3]. Extending the technique of these switching arcs to high currents (several tens of kA) and high voltages (over 100 kV) is a technological challenge that requires the help of models. Improving the design and operation of vacuum switches depends largely on a good understanding of the physics of the vacuum arc. Over the last decade, the subject has received considerable attention, progress has been made in the physical modelling of vacuum arcs [4-8] and a large number of physical phenomena have been studied. Nevertheless, it is still difficult to conclude on the influence of each specific parameter and the basis of the model used is sometimes not well explained. Due to the hermetic configuration of the vacuum systems, it is difficult to account for the physico-chemical phenomena occurring within the plasma during arc extinction.

# 2. Numerical simulation

## 2.1. Configuration

We consider two flat, parallel copper electrodes of diameter D = 5.5 cm and spaced at a height h = 1 cm.

The electrodes are in contact in a vacuum at the initial moment, outside the switching off process. During the electrical overload, the moving contact is moved by the operation of drive, and an electric arc is established which burns in the vapours resulting from the erosion of the copper electrodes. The pressure is then of the order of some $10^{-3}$ Pa. This arc can adopt two operating states: the diffuse regime and the concentrated one (at high current intensities) [4]. To avoid constriction of the plasma column and heating of the metal contacts, it is possible to apply an AMF [9, 10]. This $B_{AMF}$ field makes it possible to widen the current range in which the arc is in diffuse mode, i.e. stable in the inter-electrode medium, and to reduce the ablation of the materials. During the switching-off process, the particle flows have a helical motion and are directed from the cathode to the anode.

This type of geometry is classically studied in the literature, the work of E.Schade [4] being one of the most
complete. It mentions that two regimes are generally considered depending on the average current density at
<figure><img src="image_1.png" /><figcaption>Figure 1. Schematic of the physical model.</figcaption></figure>

the cathode. If the current density is less than 3000 $kA/m^2$, the flow is assumed to be supersonic and subsonic if the current density is higher than the given limit.

## 2.2. Equations

The mathematical and physical study of the behaviour of such a medium has analogies with that of the representation of a two-temperature plasma. The set of magneto-hydrodynamic equations that describe the arc in this configuration will be treated for two components, on the one hand the electrons and on the other hand the heavy particles (neutrals and ions) [4, 8]. The implementation of a three-dimensional numerical model describing the plasma medium as a two-temperature fluid will make it possible to predict the profiles of the temperatures, densities, pressures, and velocities of the heavy particles and electrons [5]. Today, the challenge is to be able to make the transition from a supersonic to a subsonic regime (at current densities above 3000 $kA/m^2$). For the moment, the model in place only describes the subsonic in stationary regime. To understand the approach and the equations established to obtain the following results, it is necessary to be inspired by the work of Y.Langlois and E.Schade [4, 5].

## 2.3. Boundary conditions and hypotheses

Only the plasma part is described by the equations. It is not completely a vacuum medium because, due to the erosion of the cathode, the pressure increases rapidly. It can therefore be assumed that it is a fluid. The contribution of the AMF is imposed. According to the literature, a value of 5 mT/kA is generally used [4]. The electrodes are not really represented in the domain, but are taken into account by specific boundary conditions. As already mentioned, the plasma is a frozen mixture of ions and electrons with an average charge $Z = 1.9$ [11]. A two temperature model will describe this medium.

The particles in the plasma medium come from the erosion of the cathode. The cathode is therefore considered as a mass flow inlet. In the work of E.Schade [4], this ablated area (corresponding to the inlet) represents 75 % of the surface while the remaining 25 % on the cathode side is considered as unablated. In our case, as a first approach, we consider that the ablation is governed by a Gaussian profile imposed on the whole electrode. It is set so that the value of the profile on the cathode edge is close to 10 % of the value of the axis. As explained earlier, depending on the value of the current density, the flow will be considered subsonic or supersonic. In our case, as proposed by E.Schade [4], we assume that the electron and ion temperatures are the same $T_i = T_e = 3$ eV. A copper flux is imposed with a specific value of 50 µg/C. An ablation profile and the current density profile are directly related. The Gaussian profile is therefore also used for the current density. As the self-induced magnetic field is calculated from the potential vector equations, closed by the Biot and Savart boundary condition, the current density inside the cathode is required, and this same Gaussian profile is used. The anode collects the plasma flows. We consider the anode with an output pressure condition as proposed by Y.Langlois [5] and presented in figure 2. The anode collects the plasma flows and the previous equations are used to find the exit pressure between the pre-sheath (ps) and the plasma. The velocity of the ions in the sheath (s) is assumed to be equal to the Böhm velocity and this velocity is constant in the pre-sheath. A space charge zone with a potential drop tends to adjust in order to regulate the flow of electrons collected at the anode. The pressure of the ions at the exit of the pre-sheath is to be calculated.
<figure><img src="image_2.png" /><figcaption>Figure 2. Diagram of the physical model with annotations of the boundary conditions.</figcaption></figure>

This calculation strategy comes from Y.Langlois's thesis applied to our configuration [5]. First, we calculate the conservation equation of the energy of the ions in the pre-sheath. Then we calculate the ion density equation in the pre-sheath by the Boltzmann factor. Finally, we establish the law of conservation of mass through the pre-sheath. By linking these relationships we can deduce the exit pressure by assuming that the Böhm velocity is constant in the pre-sheath. From this combination, assuming that the velocity inside the duct is the Böhm velocity, we deduce the outlet pressure. It should be noted that it is not obvious how to re-demonstrate Y.Langlois' proposal, as it only works if the Böhm velocity is assumed to be constant in the pre-duct. Nevertheless, this assumption seems reasonable.

For the magnetic field boundary condition, it is difficult to report a given current density profile within the anode. We then decided to take the current density within the plasma in the anode.

## 2.4. Simulation results and discussion

### 2.4.1. Vacuum arc modelling

The electrodes and the inter-electrode space are represented by a fine mesh as shown in figure 3 a), of 280 000 cells, with a step of 0.1 mm in the vicinity of the electrodes and a step of 0.05 mm in the middle. The results will be presented on a frontal section plane to facilitate their understanding. In figure 3 b), we observe the mesh of the system with the anode (in red) at the top and the cathode (in blue) at the bottom. The data are presented without taking into account the radiation and in a stationary state.

<figure><img src="image_3.png" /><figcaption>Figure 3. (a) 3D mesh of the physical model, (b) 2D face plane of the mesh used to represent the simulation results.</figcaption></figure>

### 2.4.2. Reference case

First simulation results of the vacuum arc in subsonic regime and in stationary in the inter-electrode space are presented. The case used in the 2003 paper by E.Schade and D.L.Shmelev [4] is our reference for this study whose characteristics are given in the table 1.

From the numerical calculations we first obtain the temperature profiles of the electrons and ions in the inter-electrode space, and then the corresponding particle densities presented in figure 4.
<table><caption>Table 1. Characteristics of the reference case, Case - 0.</caption><thead><tr><th>Case</th><th>Intensity [kA]</th><th>Axial magnetic field [T]</th><th>Ablation rate [µg/C]</th><th>Limits conditions</th></tr></thead><tbody><tr><td>Case - 0</td><td>15</td><td>0.075</td><td>50</td><td>T<sub>e</sub> = T<sub>i</sub> = 3 eV</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>Figure 4. Fields of temperature and electron and ion densities in the reference case.</figcaption></figure>

The temperature simulation results indicate that the maximum reached by $T_e$ is $7.1 \times 10^4$ K which is a relatively high value compared to the literature and which we seek to understand. The important thing here is to realise the impact of some parameters on the physics of the arc. The electron and heavy densities are consistent with the literature results and are close to the cathode with values of $n_e = 1.53 \times 10^{16} \text{ cm}^{-3}$ and $n_i = 8.26 \times 10^{15} \text{ cm}^{-3}$.

### 2.4.3. Parametric study

It will also be interesting to look the influence of the variation of the ablation rate, the strength of the AMF and the current intensity on the metal vapour plasma. Figure 5 shows the ion density of the plasma for an ablation rate of 50 and 115 µg/C. This variation corresponds to an increase in the amount of metal vapour in the vacuum chamber. As a result, the amount of $Cu^+$ and $Cu^{++}$ ions is higher and the number of collisions between charged particles increases. We also expect the temperature to increase with ablation. However, quantifying the proportion of ablation in a vacuum arc is relatively complex, especially the method of estimating the electrode material losses.

<figure><img src="image_5.png" /><figcaption>Figure 5. Fields of the ion density for an ablation rate of a) 50 µg/C and b) 115 µg/C, $I = 15$ kA, $B_{AMF} = 0.075$ T.</figcaption></figure>

Figure 6 shows the temperature profiles of the electrons and heavy particles for several values of AMF, 0.075 T and 0.125 T. The higher the AMF, the more the arc constriction will be attenuated, resulting in a more homogeneous distribution of the ion density in the plasma.

The current intensity and the AMF intensity have a direct impact on $T_e$ and $T_h$ as can be shown in figure 7. There is an increase in electron temperature of about $\delta T = 2.4 \times 10^4$ K between 15 kA and 25 kA.
<figure><img src="image_6.png" /><figcaption>Figure 6. Fields of the electron and ion temperature for several values of the axial magnetic field, $B_{AMF} = 0.075$ T and $B_{AMF} = 0.125$ T, $I = 25$ kA, the ablation rate is 50 µg/C.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 7. Fields of the electron and ion temperature for several values of current intensity, $I = 15$ kA and $I = 25$ kA, $B_{AMF} = 0.075$ T, the ablation rate is 50 µg/C.</figcaption></figure>

The value of the current intensity will also have an impact on the physical components of the system, such as the temperature of the electrons and the heavies, see figure 8, with constant rate of 5 mT/kA.

<figure><img src="image_8.png" /><figcaption>Figure 8. Fields of electron and ion temperature for several values of current intensity $I = 10$ kA, $I = 15$ kA, $I = 25$ kA, $B_{AMF} = 0.075$ T, ablation rate is 50 µg/C, rate is 5 mT/kA.</figcaption></figure>

Whether electrons or ions, there is a significant increase in $T_e$ and $T_h$ with the value of the current intensity. This is due to the higher collisional energy transfer between the charged particles. The problems encountered in this simulation of the arc are still numerous, in particular the one related to the values reached for the temperature.

### 3. Conclusion

The simulation behaviors are consistent with the literature in a subsonic and stationary arc mode. The supersonic mode needs to be developed as well as the transition between the two regimes. Studies need to be carried out on the contribution of radiation and electrode ablation. A more complex geometry should be considered so that the results can be closer to reality. We will also try to vary the inter-electrode distance by a few centimetres to observe its impact on plasma parameters.

### References

[1] J. officiel de l'Union européenne. Règlement (ce) no842/2006 du parlement européen et du conseil du 17 mai 2006relatif à certains gaz à effet de serre fluorés, 2006. URL:
[2] R. Renz, D. Gentsch, H. Fink, P. Slade, and M. Schlaug. Vacuum interrupters-sealed for life. In *CIRED-19th Int. Conf. on Electricity Distribution*, number 0156, 2007.

[3] N. Wenzel, S. Kosse, A. Lawall, R. Renz, and W. Hartmann. Numerical simulation of multi-component arcs in high-current vacuum interrupters. In *2012 25th International Symposium on Discharges and Electrical Insulation in Vacuum (ISDEIV)*, pages 321–324. IEEE, 2012.

[4] E. Schade and D. L. Shmelev. Numerical simulation of high-current vacuum arcs with an external axial magnetic field. *IEEE Transactions on Plasma Science*, 31:890–901, 10 2003.

[5] Y. Langlois, P. Chapelle, A. Jardy, and F. Gentils. On the numerical simulation of the diffuse arc in a vacuum interrupter. *Journal of Applied Physics*, 109, 6 2011.

[6] W. Hartmann, A. Hauser, A. Lawall, R. Renz, and N. Wenzel. The 3d numerical simulation of a transient vacuum arc under realistic spatial amf profiles. In *24th ISDEIV 2010*, pages 285–288. IEEE, 2010.

[7] S. Jia, L. Zhang, L. Wang, B. Chen, Z. Shi, and W. Sun. Numerical simulation of high-current vacuum arcs under axial magnetic fields with consideration of current density distribution at cathode. *IEEE Transactions on Plasma Science*, 39(11):3233–3243, 2011.

[8] S. Braginskii. Transport processes in a plasma. *Reviews of plasma physics*, 1:205, 1965.

[9] L. Wang, S. Jia, L. Zhang, D. Yang, Z. Shi, F. Gentils, and B. Jusselin. Current constriction of high-current vacuum arc in vacuum interrupters. *Journal of Applied Physics*, 103(6):063301, 2008.

[10] L. Wang, D. W. Matson, E. Polikarpov, J. S. Swensen, C. C. Bonham, L. Cosimbescu, J. J. Berry, D. S. Ginley, D. J. Gaspar, and A. B. Padmaperuma. Highly efficient blue organic light emitting device using indium-free transparent anode ga: Zno with scalability for large area coating. *Journal of Applied Physics*, 107(4):043103, 2010.

[11] J. Kutzner and H. C. Miller. Integrated ion flux emitted from the cathode spot region of a diffuse vacuum arc. *Journal of Physics D: Applied Physics*, 25(4):686, 1992.