Mathematical Modeling on Dynamic Characteristics of the Breakdown Process in Narrow-Gap of SF 6 Based on the FCT Algorithm

The two-dimensional and self-consistent fluid model of the SF 6 discharge was established based on the electron and ions continuity transfer equations coupled to Poisson’s equation, and also simultaneously considered the photoionization event, and then the flux corrected transport technique (FCT) was employed to numerically solve the particle flux-continuity equations, and some significant microphenomena were achieved that the dynamic behaviors of the charged particles, the spatio-temporal evolution of the discharge channel and the transformation law of the avalanche-streamer for the SF 6 nar-row-gap were revealed in this paper.


Introduction
Usually, the approaches to explore gas discharge phenomenon can be classified into experimental and theoretical methods. Experimental method can directly observe the discharge events, but the discharge mechanisms have not been clearly discovered in detail. However, the theoretical method not only gains key data for the dielectric breakdown of gases, but also some microcosmic parameters on controlling dynamic behaviors of charged particles produced in the time of the gas discharge process are founded [1] [2] [3]. Recently, the computational fluid dynamics has been widely used in the gas discharge field due to its advantages in explicit physical concept, legible image exhibition and highly calculating efficiency [4] [5] [6].
Despite the fact that the electron and ion densities can attain very steep gradients and make the shock fronts in particule fluids along the discharge channel when the gas discharge is triggered and developed in the overall process under the stimulation of electric field, the FCT technique (flux corrected transport) shows high efficiency and accuracy dealing with the intricacy situation [7] [8] [9]. By virtue of the excellent thermodynamic, dielectric and transport propeties, SF 6 is widely used in the eletrical industry, especially in the high voltage circuit breaker technology, although it might be replaced by the environment friendly insulating gases. Nevertheless, it is necessary to accurately find microcosmic mechanisms of the particle dynamical behaviors, the spatio-temporal characteristics of the electric field and the track evolution of the discharge channel along the direction of the specific discharge development for broader application domain.
This paper is organized as following: the mathematical model, the FCT algorithm and constraint conditions are carefully presented in Section 2, and the results of the mathematical simulation are analyzed in Section 3 including the following contents: 3.1, avalanche phase; 3.2, streamer formation phase; 3.3, the discharge channel evolution and the photoionization effect. Some valuable conclusions are given in Section 4.

Model of the Gas Discharge
The model for the narrow-gap with parallel plate electrodes which are filled with SF 6 gas has been presented in detail [10] [11], and the spatio-temporal evolution of the SF 6 discharge overall process is mathematically modeled by a set of equations governing the transport of particles, moment and energy for ions and electrons together with the electric field equation as follows [12]   reaches the ionization critical value, then the photoionization process occurs and produces a certain amount of the photo electron. Generally speaking, the number of photon electron is much smaller than that of the electrons generated by impact ionization. However, these photon electrons would result in the formation of the secondary avalanche, which moves towards the head of the first avalanche, composes a big avalanche and accelerates the development of the discharge process. Hence the photoionization effect plays a significant role in the gas discharge process. The term ph S , being a source term due to photoionization, is represented as shown below: here p γ , * α and µ are second ionization, excitation and absorption coefficients, and Ω is the solid angle subtended at z′ by the disk charge at the point z. A detail solution about ph S can be referred to the literature [15].
Taking into account the distortion of space charge effects on the electric field, the Poisson equation is given by: where ϕ is the electric potential; q is the electronic charge; 0 ε is the permittivity of the free space. The current I in the external circuit due to the motion of electrons and ions between the electrodes is calculated by the Sato formula [16]: here, r is the radius of the discharge channel, q is the electronic charge.
The convection term of particle's continuity Equations (1)

FCT Algorithm
The convective terms of Equations (1), (2) and (3) are written as where symbol N shows the density of the particle species and v is their velocities.

Constraint Conditions
In this paper, schematic diagram and the constraint conditions for parallel-plate electrodes discharge under atmospheric pressure are set as shown in Figure 1, and the calculation model is converted to a two-dimensional structure by rotating the axis of symmetry. In the present study, the gap distance between the electrodes filled with SF 6 gas is 5 mm under the pressure 0.1 MPa and temperature × is taken as time step, which is significantly smaller than that required for stability of the used numerical scheme [3]. At the initial moment of time for gas discharge, the quasi-neutral plasma spot of Gaussian shape in the radial and axial directions is placed at the front of cathode: where , r z is the radial and axial coordinates respectively; the origin of coor- ) is positioned at the center of cathode surface, the peak value density of particles (seed electrons and positive ions) is The solution of Poisson's equation is subject to the following boundary conditions: where 0 V is the applied voltage, r is the radius of the computational domain.

Modeling Results
The space between parallel-plane electrodes 0.5 cm apart is filled SF 6 gas with standard atmospheric pressure and commercial purity, and the 46 kV DC is applied on the anode plane, then the cathode grounding. The incepting discharge of the model is trigged by the seed electrons having a Gaussian distribution near the cathode at 0 ns t = (located at 0.1 cm from the cathode as mentioned earlier). Under the electric field, the seed electrons obtain energy to migrate, impact the neutral molecules and produce much more charged particles, so that the current through parallel-plane electrodes is shown in Figure 2. The transition time of the electron for the model is about 10 ns, according to the principle of gas discharge, the whole discharge process undergoes two phases from the avalanche to the streamer phase while the discharge incepted.

Avalanche Phase
The electric field stress 92 kV/cm applied to the SF 6 gap, which slightly larger than the threshold value 89.6 kV/cm, easily renders the discharge happen smooth-ly. With the electrons migrating to the anode and impacting with neutral molecules, an electron swarm is quickly made up as shown in Figure 3. When about 7.7 ns moment, the peak value of the electron swarm is up to the 18 3 1.01 10 cm − × at the point 0.25 cm z = apart from the cathode, and corresponding positive and negative ions densities are 18 3 3.97 10 cm − × and 18 3 3.22 10 cm − × respectively as shown in Figure 7 and Figure 8, then the distortion of the space electric field is not remarkable as shown in Figure 4 at 7.5 and 7.7 ns moments, therefore the interval from 0 ns t = to the 7.7 ns t = usually is known as the avalanche phase according to the gas discharge theory [18].

Streamer Phase
When the avalanche volume reaches to the critical value, it instantly changes the avalanche phase into the streamer phase, namely, the streamer phase formation [19] [20]. It is demonstrated in Figure 3 and Figure 4 that the electron densities are visibly larger than that of the avalanche phase in the interval 7.7 -8.0 ns, although the distorting of space electric field is still slight enhancement, it is considered an initial period of the streamer phase. From the 8.0 ns moment, the streamer discharge is rapidly development, its volume is promptly expanded and the length is also stretched quickly, and simultaneously the electric field distortion is also exacerbated as shown in Figure  5, then the mainly changes of the density distribution of charged particles for the electron, positive and negative ions on the axial direction are shown in Figures  6-8, the procedure described is rather intricacy including the dynamical behaviors of the charged particles in the electric field [21] [22].

1) The anode-directed streamer
When the transformation avalanche into the streamer phase, the peak value of the electron densities from location 0.25 cm z = at 7.7 ns to the point 0.44 cm z = at 9.2 ns t = is shown in Figure 6, and the streamer head moves to the anode at an average speed 6 1.33 10 cm s × . The positive and negative ion densities in the anode-directed streamer head increase steadily from 18 3 4.43 10 cm − × to 19 3 1.09 10 cm − × and from 18 3 3.62 10 cm − × to 19 3 9.41 10 cm − × respectively shown in Figure 7 and Figure 8. The head of the anode-directed streamer shows the electro-negativity due to gathering most electrons and greatly enhanced field intensity of the short distance between the head with the anode plate, so that SF 6 molecule ionization in this small region is dramatically aggravated, and many more electrons are reproduced and the head radius continuously becomes bigger and bigger, the iterative process finally ends 2) The cathode-directed streamer As shown in Figure 6, the curve of electron distribution at 0.18 cm z = and 8.3 ns t = presents the clearly escalating trend, that is to say, the event of the cathode-directed streamer obviously happens, then the appeared moment of the cathode-directed streamer is later than that of the anode-directed one. In contrast with anode-directed streamer, the head of the cathode-directed streamer presents the electro-positivity, because the electric field is strengthened in the short space between the cathode plate and the head, then the ionization rate of the SF 6 molecule is greatly accelerated and much stronger, the velocity of the cathode-directed streamer at 9.2 ns is 6 0.92 10 cm s which is only about 70% of the velocity of the anode-directed streamer in the same moment.

Discharge Channel and Photoionization
As mentioned above, the discharge process in the narrow-gap of the SF 6 undergoes transformation from the avalanche to the streamer phase, once the streamer is triggered it soon develops respectively towards the anode and cathode plate, namely, the anode-directed streamer and the cathode-directed streamer, at the same time accompanied by the photoionization appeared, the streamer volume not only grows quickly presenting a near cylindrical shape and but also the streamer length becomes much longer till through the anode and cathode plate, then the discharge path of the both electrode is formed and known as a breakdown channel. The fact has been proved by the experiment results and theoretical demonstration, the photoionization plays an important role and not been overlook when the whole discharge process in the narrow-gap of the SF 6 undergoes transformation from the avalanche to the streamer phase. The results of the mathematical modeling in this paper shows in Figure 9, considering and ignoring photoionization could achieve different results, if you ignore the photoionization effect, you might get the false result [23].

Conclusions
Based on the fluid model of gas discharge and the FCT algorithm, dynamic characteristics of the SF 6 breakdown process in narrow-gap has mathematically modeled and some important facts has demonstrated that the FCT algorithm is an efficient theoretical way to deal with the troublesome problems having shock fronts in the discharge channel. The results are both shown via easy visualization for the complicated course of the SF 6 discharge and revealed dynamic characteristics of the charged particles during the SF 6 discharge process in narrow-gap.
According to the mathematical model and numerical analysis results in this paper, the breakdown process in narrow-gap of SF 6 still presents two phases at the standard atmospheric pressure, and parallel-plane electrodes 0.5 cm apart and the 46 kV DC applied. Moreover, some facts are indicated that when electrons increase to a certain amount and the avalanche phase changes into the streamer phase. On the one hand, in the electron avalanche phase, the collision ionization is the key rule for producing electrons and making them grow fast until the electron swarms up to the 18 3 1.01 10 cm − × . The microcosmic mechanisms of the particle dynamical behaviors, the spatio-temporal characteristics of the electric field and the track evolution of the discharge channel are shown in Figure 3 and Figure 4 in the interval from 0 ns t = to the 7.7 ns t = . In this interval, the collision between SF 6 molecule and seed electrons is only triggered by the external electric field. When time ranges from 7.7 ns t = to the 8.0 ns, a dramatic change happens in the electron number and electric field of the gas-gap; that is to say, the streamer discharge would be coming. On the other hand, in the streamer phase, from the 8.0 ns moment, the streamer discharge is rapidly developing; its volume is promptly expanded and the length is also stretched quickly, and simultaneously the electric field distortion is also exacerbated as shown in Figure 5, then the main changes of the density distribution of charged particles for the electron, and positive and negative ions on the axial direction are shown in Figures 6-8.
In a nutshell, the photoionization effect plays a pivotal role in the process of