Next Article in Journal
Improvement of the Photocatalytic Activity of Au/TiO2 Nanocomposites by Prior Treatment of TiO2 with Microplasma in an NH3 and H2O2 Solution
Previous Article in Journal
Quantum Matter Overview
Previous Article in Special Issue
On the Tree Gauge in Magnetostatics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling of Positive Streamers in SF6 Gas under Non-Uniform Electric Field Conditions: Effect of Electronegativity on Streamer Discharges

1
G2Elab, CNRS/Grenoble INP/Université Grenoble Alpes, 38000 Grenoble, France
2
Department of Electronic and Electrical Engineering, University of Strathclyde, Glasgow G1 1XW, UK
*
Author to whom correspondence should be addressed.
Submission received: 2 March 2022 / Revised: 22 April 2022 / Accepted: 28 April 2022 / Published: 9 May 2022
(This article belongs to the Special Issue Computation of Electromagnetic Fields)

Abstract

:
The use of SF6 in electrical insulation and fast-switching applications cannot be overemphasized. This is due to its excellent dielectric properties and high breakdown voltage, which are especially important for practical applications such as gas-insulated switchgears and pulsed power switches where pressurized SF6 is used. Breakdown in the gas occurs via streamer–leader transition; however, this transition is difficult to quantify numerically at atmospheric pressure because of the electronegativity of the gas. In the present work, streamer discharges in SF6 gas at pressures of 10 and 100 kPa were studied using a plasma fluid model implementation. Analysis of the electric field in the streamer body, streamer velocity, diameter, and the effect of the high electronegativity of the gas on streamer parameters are presented for positive polarity in a point-to-plane geometry. The streamers in SF6 for non-uniform background fields are compared to those in air, which have already been studied extensively in the literature.

1. Introduction

SF6 gas is used for different HV applications, including gas-insulated switchgear (GIS, circuit breakers, switches) and gas-insulated transmission lines (GIL), because of its excellent dielectric and arc quenching properties [1,2]. Reference data relating to the breakdown properties of the gas in these applications have been extensively investigated and presented in the literature. Notably, it has also been widely adopted in fast-breaking plasma-closing switches for pulsed power applications because the switching element exerts a strong influence on the rising time and amplitude of the pulse [3,4]. These switches generally work under high voltage and current conditions and are prone to breakdown and discharge processes both thermal and non-thermal because of the extreme operational conditions and thus require a careful consideration of all discharge phases [5] (and references there-in after).
On finding possible replacements that satisfy all necessary technical and ecological requirements for the greenhouse gas, a better understanding of the pre-breakdown and breakdown processes in SF6 is required. For the present contribution, the main focus is placed on the numerical modelling of fast-transient plasma discharge processes—streamers, their evolution in SF6 gas, and how that compares with a reference gas. A comparative numerical study has been carried out on the plasma discharge behaviour in air and SF6 in conditions relating to pulsed power applications: short separation distances, high instantaneous voltage application and non-uniform configuration, to minimize the time delay for discharge formation. The upper bounds of pressure, especially for SF6 and its influence on discharge behaviour, has also been analyzed, highlighting major physical concepts and some challenges in modelling of streamer discharges at 1 bar in the gas.
Limited literature on numerical modelling of streamers in SF6 under non-uniform fields is available, highlighting the complexity of such simulations and justifying the aims and objectives of the current study. The insights obtained from the comprehensive behaviour of an electronegative gas under high voltage stress, the streamer initiation, and propagation physics will help in identifying new gas mixtures for use in electro-technical applications.
A review of numerical modelling of streamers in SF6 is given in Section 2. Section 3 provides a brief description of the mathematical model with the swarm parameters used in the present study. The obtained results for streamers in air and SF6 at 10 and 100 kPa are reported and discussed in Section 4, and the work is concluded in Section 5.

2. Review of Numerical Simulations of Streamers in SF6

From an experimental point of view, previous studies have established that dielectric breakdown in SF6 develops through formation and propagation of a stepped leader (streamer-leader transition) [6,7,8,9]. This is further highlighted by Chalmers et al. in [10], where it was found that the development of leaders in a point to plane SF6 gap at 500 kPa produced a streamer during the leader step development. In the recent work by Bujotzek et al. in [11], positive and negative streamer radii and their propagation lengths have been obtained at gas pressures between 50 and 100 kPa using strong and weak non-uniform background electric fields.
From a numerical modelling point of view, limited literature exists on streamers in pure SF6 and SF6 gas mixtures.
In [12], Morrow reviewed the dominant physical processes that affect streamer formation in SF6 and highlighted the very high value of the electron attachment coefficient leading to a rapid formation of negative ions in the gas at atmospheric conditions. Since corona discharge formation is a time-dependent process, the effect of the high values of the attachment coefficient leading to the creation of negative ions was taken into account by evaluating the characteristic attachment time and subsequently the rate of change of electron density due to the attachment process.
One-dimensional continuity equations for electrons, cations, and anions including ionization, photoionization, attachment, recombination, and electron diffusion terms were used with the continuity equation for electrons being of the 2nd order and that of the ions, 1st order. The considered cases were for a uniform electric field with both negative and positive streamers, propagating in a 5 mm gap, and a non-uniform configuration with a 5 mm diameter needle electrode, 65 mm away from the grounded electrode. In the latter, a constant voltage of 50 kV was applied, and the discharge evolution studied for 3 ns before it stalled. The influence of a positive high voltage impulse was also studied as an applied voltage of 200 kV with a linear rise time of 15 ns enabled a longer propagation in the gap. To illustrate the re-illumination of the streamer channel and possibly stepped leaders, the current in the external circuit due to the movement of electrons and ions was computed using a modified version (with inclusion of negative ions and electron diffusion) of the Sato equation [13]. The computed current fell as the streamer propagated and, after a certain time period, started to pulse regularly. This was, however, inconclusive as the model did not take into account background heating of the gas.
Two-dimensional representation of streamers in SF6 at atmospheric pressure in a uniform field was studied in [14]. Similarly, the continuity equations were used in [14], but no photoionization was included. Instead, a background ionization was used in this model, which was implemented by introducing background charged particles (electrons and positive ions) uniformly distributed throughout the gap with the number density of 1010–1014 m−3. This provided the benefit of understanding the dependence of the streamer propagation on the ionization density in front of it. A Flux Corrected Transport (FCT) technique coupled with a finite difference method for discretization was used in the solving of the higher order equations. The Poisson’s equation, which was used for obtaining the electric field, was solved using a fast Fourier transform in the z-direction and a cubic spline interpolation in the r-direction.
A neutral ionization density was introduced in the gap and an applied voltage of 50 kV was used. Both cathode and anode streamers were observed, and the maximum field enhancement was observed at the front of the streamers. More qualitative occurrences, such as the increase in the velocity as the streamer approached an electrode and the increase of negative ions as electrons attached rapidly due to a lower electric field in the bulk of SF6, were also observed.
The constant streamer radius model developed in [15] accounts for the heating and radial expansion of the gas in the channel and describes the stepped pattern in propagation of positive streamer at 200 kPa in a sphere to plane configuration. In this work also, the continuity equations were used with an extra reaction term to describe the generation of precursor electrons ahead of the streamer front due to ionization of gas molecules by photons produced in the streamer head. An examination of the discharge current showed a series of short pulses, with each pulse resulting in an increase in gas temperature in the channel and with the most intense heating taking place at the anode. The gas heating resulted in a radial expansion of the streamer channel and decrease in the gas density inside the streamer, but this was limited by the initial streamer radius that was specified, resulting in the constant streamer radius model. With the reduction in the gas density, an increase in the reduced electric field and subsequently in the ionization integral is attained. When the ionization integral reached a specific value, the second streamer formed and propagated through the gas with lower density. The second streamer thus propagated further than the preceding one.
The studies conducted in [16,17,18] report on modelling of 1D and 2D streamers in mixtures of SF6 with air at 0.1 kPa, with N2 at 100 kPa and 200 kPa. For the 1D implementation, the swarm parameters of SF6 and the gas mixtures were computed by the Monte Carlo Simulation (MCS) method where only the collisions between electrons and neutral particles were considered. The implementation of the continuity equations, however, duplicates the work done in [15]. The 2D admixture approach developed in [17] for negative streamers was implemented in COMSOL™ Multiphysics for a non-uniform gap. The convection-diffusion equations were coupled with the Poisson’s equation, and the influence of photoionization was omitted in this work in order to simplify the model presented [17]. Further information on the evolution of the negative streamer, the effect of the mixed gas ratio, and effect of electrode shape has been provided. In [18], the fluid model was coupled with 47 chemical reactions, and the reason for streamer decay in SF6/N2 mixtures with high and low SF6 content was assessed.

3. Simulation Model

In the present work, the continuity equations for electrons, positive ions (cations), and negative ions (anions) in space and time are solved [19,20], taking into account drift velocity, diffusion, ionization, attachment, recombination, and background ionization (used in place of photoionization) terms.
n e t + · D e   n e µ e n e E = S b + α η n e   µ e   E β e p   n e   n p ,
n p t + · µ p n p E = S b + α   n e   µ e   E β e p   n e   n p β p n   n p   n n ,
n n t + · µ n n n E = η   n e   µ e E β p n   n p   n n .
The subscripts e , p , and n denote electrons, cations, and anions respectively; n is the particle number density (cm−3); µ is the mobility of the charged species (cm2/V·s); D is the diffusion coefficient (cm2/s); E is the electric field (kV/cm); t is the time (s); α   is the Townsend’s ionization coefficient (1/cm); η   is the attachment coefficient (1/cm); and β   is the respective recombination coefficient (cm3/s). S b is the background ionization source used to model the photoionization processes. The logarithmic values of the charged species were used. This is useful when steep gradients are involved in the transported scalar quantity and reduces the need for artificial diffusion terms while ensuring that density of species always remain positive. Validation of the logarithmic implementation of the drift diffusion equations can be found in [19]. The convention to use a background ionization term as opposed to a full photoionization model is already established in literature notably in [21] where the importance of photoionization and background ionization in air and N2-O2 mixtures for pulsed repetitive discharges were investigated; in [22] for N2-O2 mixtures comparing different streamer codes and also used as a convention for gases where the photoionization process is not well established [14,23].
At every time step, the electric field is also computed using Poisson’s equation for electric potential ϕ .
2 ϕ = e ϵ 0 n e + n n n p ,
E = ϕ ,
where e is the elementary charge of an electron and ε 0   represents the permittivity of vacuum.
Figure 1, Figure 2 and Figure 3 show the reduced effective ionization α e f f = α η N , reduced electron mobility µ e N , and reduced electron diffusivity D e N as functions of the reduced electric field E / N for SF6 and air. Air has been included in this analysis as a reference gas since streamer behaviour in air is well documented and understood. As shown in Figure 1, the critical reduced electric field of SF6 is ~362 Td as opposed to ~120 Td for air. This is the first evidence of the high attachment rate in the highly electronegative gas. The full list of parameters for the gases can be found in Appendix C and Appendix D.
To mimic the non-uniformity in the electric field distribution that may occur in electrical equipment due to defects, a needle—plane geometry was simulated in COMSOL™ Multiphysics [28] with the needle (radius of curvature of 0.1 mm) serving as the high voltage electrode and the plane electrode grounded. The plasma module of COMSOL™ was used for all physics implementation and calculations. This module couples the drift-diffusion equations, heavy species transport, and electrostatic interface into an integrated multiphysics model, which allows for comprehensive analysis of plasma discharges. The discretization formulation used was finite element, log formulation (linear shape function). This uses the Galerkin method to discretize the equations.
An axisymmetric implementation of a point-plane electrode geometry in 2D space was chosen for the model. This resolves the model in cylindrical coordinates, and the initial conditions are mostly chosen in a way that the streamer develops along the symmetry axis, and therefore only the positive half axis of the model has been represented. This approach reduces the demand for computational resources in comparison to a full 3D model. However, for the post processing and output presentation, a mirror approach was adopted, and both the positive and negative half axes are shown.
As the boundary conditions used for the simulations, a voltage of positive polarity was applied to the needle electrode with the planar electrode grounded. The needle served as an outflow for the electrons and negative ions (Neumann condition) and the plane functioned as the outflow for positive ions (Neumann condition). If an electrode does not serve as an outflow for a charge specie, a Dirichlet zero condition was imposed. Open (zero flux) boundary conditions were applied to the outside planes of the computational domain.
In the present simulations, the computational domain was 5 × 7 mm (r × z) with the gap distance between the needle electrode and the plane fixed at 5 mm. The needle radius of curvature was 0.1 mm. An adaptive mesh was used to follow the propagation of the streamer in this gap. Being a finite element method simulation, triangular meshes were used for the discretization. To be able to run simulations for relatively larger gaps and computational domains, the computational domain was sub-divided. This provided the advantage of utilizing different mesh sizes depending on activity levels. In a 1.5 mm domain in the r-direction, a cell grid ranging from 2 μm to 8 μm was utilized; beyond this region, the grid expanded according to geometric progression. An adaptive mesh was used to efficiently follow the propagation of the streamer and to accurately resolve the charged layers in the vicinity of the electrodes and in areas with high gradients of charge density. The error estimation for the refinement was done with the electron reaction rate as the indicator. A minimum mesh size of 2 μm was used for this purpose. (See Figure 4).
For the time stepping method, a BDF formula with adaptive time stepping was used. The adaptive time stepping allows the solver to take larger or smaller time steps as required to satisfy a specific tolerance. The maximum and minimum time steps were 10−9 and 10−15, respectively. The steps taken by the solver were automatically adapted within these bounds at each iteration to find a suitable solution. An output step time of 0.1 ns was used for post processing of the time dependent simulation.

4. Results and Discussion

Streamer discharges have been modelled in both air and pure SF6 gas; the needle electrode was stressed with a constant positive voltage. Development of streamers have been studied at the ambient gas temperature of 293 K and pressure values of 10 kPa and 100 kPa corresponding to the number density values of 0.2472 × 1025 m−3 and 2.472 × 1025 m−3. A background ionization in the present model is taken into account by the charge density term in Equation (1), S b = 1023 m−3 s−1. No Gaussian concentration of charge species was introduced for the initiation of the streamer. The gap between the point and the plane electrode was 5 mm.
For the results presented, the streamers are characterized in terms of the streamer channel length, l , which is the separation between the needle electrode and the streamer head, streamer velocity v = d l d t , maximal field enhancement in the streamer head E m a x , streamer radius r , which was obtained by measuring the radial extension of the electric field in the head of the streamer, diameter d = 2 r , and the electron density. Streamer velocities presented were computed when the streamer crossed the gap, and the radius/diameters were measured at axial position 2 mm.

4.1. Atmospheric Pressure (100 kPa)

Streamers in electropositive and weak electronegative gases such as N2 and air at atmospheric pressure have been extensively studied both numerically and experimentally, however very little information on the development of streamers in highly electronegative gases such as SF6 is available. To adequately understand the effect of the electronegativity on the streamer initiation and propagation process, a comparison has been made between the inception and propagation phases in air and SF6 for an applied voltage of 20 kV with a maximum electric field of 610 kV/cm at the tip of the needle electrode. With a 5 mm gap, the nominal average electric field in the computational domain was obtained by dividing the applied voltage, V, by the gap distance, d, E a v = V d , is E a v = 40 kV/cm. The high applied voltage and the non-uniform electrode configuration ensures that the electric field in the domain, specifically near the tip of the needle, is greater than the dielectric breakdown field for air. The average electric field E a v ,   is 40 kV/cm, which is equal to or higher than the stable electric field required for streamer or leader propagation in SF6 which is between 30–40 kV/cm [29,30].
Figure 5, Figure 6 and Figure 7 show the 2D surface plots for the log of the electron density, the electric field distribution, and the ion densities in air at 0.5, 1 and 2 ns. The logarithmic point 1 represents 10 m−3 while the point 20 represents 1020 m−3. The equivalent surface plots for SF6 are shown in Figure 8, Figure 9 and Figure 10.
For air, the initiation and propagation of the full plasma channel is attained. In Figure 5, the inception of the plasma front and subsequent propagation in the 5 mm air gap can be observed. The electron density in the ionized gas channel behind the propagating plasma front is ~1020 m−3. The different discharge phases are highlighted by the electric field diagram in Figure 6 where the electric field achieves its peak values at the earlier times steps, then reduces when the electrostatic coupling between the streamer head and anode reduces and increases again as the streamer approaches the cathode. The head of the streamer in the stable propagation mode has an electric field of ~150 kV/cm. A comparison between the plots for the number densities of cations and anions at 2 ns (Figure 7) show the higher concentration of positive ions to negative ions. The rate of production of positive ions is equivalent to the rate of production of electrons as both of these parameters are governed by the ionization coefficient. The rate of production of negative ions, on the other hand, is defined by the attachment coefficient which is lower in a weak electronegative gas such as air. This accounts for the difference in the ion densities.
For SF6, at the early time steps when the discharge initiates, the electric field is strengthened by the applied voltage and, consequently, the rate of ionization far exceeds the attachment rate. As the streamer propagates, the electrostatic coupling between the streamer head and the needle electrode reduces. This causes a strong affinity for attachment resulting in a high concentration of negative ions in the bulk of the streamer (See Figure 10). In the electron density profile diagram shown in Figure 8 (in logarithmic scale), a discontinuity can be observed between the head and the tail of the streamer resulting from the high level of attachment in the electronegative gas. In the streamer head however, due to high electric field and high ionization rate, the electron density is high and has almost the same order of magnitude as the positive ion density. The corresponding electric field plots in Figure 9 show the enhancement of the electric field in the streamer head. A thinner streamer head with the electric field concentrated in the tip is observed as the attachment process dominates in the bulk of the streamer.
Figure 10 highlights the ion distributions in the streamer and, owing to the low drift velocities of ions and the high electron attachment rate, similar profiles are observed for both the positive and negative ions except at the tip.
Unlike in air where the plasma front crosses the gap, in SF6, the streamer stops propagating at about 1 mm away from the needle electrode. This can be due to the limitations of the computational model as it was shown experimentally that in SF6 at 100 kPa and 1 cm, a leader discharge is formed, (See Figure 11); thus, the continuity equations used in the present work cannot accurately describe the complete breakdown process. In the case of leaders, the heating of the gas should also be taken into account.

4.2. Sub-Atmospheric Pressure (10 kPa)

The simulations were repeated in air and SF6 at the lower pressure of 10 kPa and for an applied step voltage of 5 kV with a Laplacian electric field at the tip of the needle electrode of 155 kV/cm. Similar to streamers in air at atmospheric pressure (100 kPa), the plasma front in air at 10 kPa has a radius of 1.6 mm and was observed to cross the complete inter-electrode space (5 mm gap). The streamer radius was obtained by measuring the radial extension of the electric field in the head of the streamer. Figure 12 shows the line diagram of the electron density profile in its logarithmic form at various time steps and the corresponding electric field distribution. At steady state propagation of the streamer with an average velocity of 1.3 mm/ns, the electric field in the head of the streamer is 23 kV/cm. The electron density values attained are an order of magnitude lower than that obtained for the streamers in air at 100 kPa. Nevertheless, the general behaviour of the streamers at both pressures is the same. The flat profile of the electron density in the streamer body signifies a steady conduction path between the point electrode and the streamer head. Comparing the positive and negative ion densities, the higher ionization rate as compared to the attachment rate is highlighted in Figure 13 where the line diagrams of the ion densities are presented in a logarithmic form for identical time steps. This correlates well with results obtained in air at 100 kPa, where a similar order of magnitude difference (2) between the positive ion density and the negative ion density was realized (see Figure 7).
In SF6 at pressure of 10 kPa, streamers that behave similarly to those in air were attained in the 5 mm gap when the applied voltage was 5 kV. The 1D line diagrams for the log of the electron density profiles and the electric field are shown in Figure 14. In Figure 15, the corresponding ion densities are presented. It was also confirmed experimentally that at gas pressure of 10 kPa, streamer discharges are obtained (see Figure 16) as opposed to leaders observed at atmospheric pressure.
Focusing on the profile of the electron density of the streamer, the peaks represent the streamer head and the troughs represent the reduced electron density in the streamer body. As the streamer head approaches the cathode, there is a gradual increase in the electron density, a consequence of the field difference between the head and the cathode causing an increase in electron ionization. This is different from the electron density profile obtained for a streamer in air, which has a pseudo flat profile apart from when it interacts with the cathode (refer to Figure 12).
The electric field in the streamer head in the stable propagation mode is ~27 kV/cm (Figure 14). Thus, the reduced electric field at 10 kPa is ~1092 Td, which is well above the critical reduced electric field of 362 Td required for possible streamer propagation in SF6 at atmospheric pressure. The high reduced electric field is required to maintain the higher ionization effect in the electronegative gas. The computed streamer velocity and radius are 0.96 mm/ns and 0.8 mm, respectively.
The ion densities in Figure 15 show a flat profile for the positive ion density, representing the constant nature and low drift of positive ions in the streamer body as the streamer propagates towards the cathode. Similar to the electron density, there is an increase in the positive ion density close to the cathode due to increased ionization rate. The number density of negative ions, on the other hand, increases in the streamer body with each passing time step due to the high electron attachment rate. Both ion densities are of the same order of magnitude unlike in air, where the positive ion density was higher than the negative ion density (refer to Figure 13).

4.3. Influence of Voltage

The effect of the applied voltage on the development of streamers in SF6 at a pressure of 10 kPa has been studied. The needle electrode was stressed with 5 kV and 10 kV, and the gap distance was 5 mm. As expected, the streamer forms early and travels faster with an increased applied field demonstrating a streamer velocity of 4.17 mm/ns for the applied voltage of 10 kV as opposed to 0.96 mm/ns seen in the preceding section for 5 kV. Additionally, it is interesting to discuss the electron density in the streamer body. As the voltage increases, resulting in an increase in the electric field, the rate of ionization far exceeds the rate of attachment, thus higher electron densities are realized in the streamer channel. The shape of the electron density profile resembles the flat profile characteristic of streamers in air (see Figure 12). This is shown in Figure 17 where the line diagrams of the positive ion density (blue dash-point lines) and electron density (solid gray lines) are plotted on the same figure for 5 kV (left) and 10 kV (right). The streamer head in this case is thicker at lower voltages than at higher ones.
Comparing the electric field in the streamer head in the steady state propagation mode, for a voltage of 5 kV, the electric field in the streamer head at steady state propagation was 27 kV/cm as opposed to 54 kV/cm for the applied voltage of 10 kV. Figure 18 highlights this effect at 10 kPa.
The streamer radius and diameter are directly proportional to the applied voltage, and these parameters increase with an increase in the applied voltage, which is a typical dependency for streamers in other gases. This is because most of the swarm parameters used in the computational model are functions of the reduced electric field and increase as the voltage and electric field increase while the particle number density remains the same.
When the applied voltage was scaled accordingly to have the same reduced electric field at 10 kPa as at 100 kPa, similar numerical results were obtained for both gas pressures. Recalling, the propagation of the streamer was halted at 20 kV and 100 kPa as the electron density and subsequently conduction in the streamer channel reduced (see Figure 8). The 1D plot of the log of the electron density at 10 kPa and for an applied voltage of 2 kV in Figure 19 highlights this effect as the electron density in the streamer channel reduces substantially with time.

4.4. Comparison with Experimental Results

For streamer characterization, one of the most commonly used parameters is the streamer radius, especially where there is the need to compare numerical results with experimental ones. The radius of the streamer is controlled by the diffusivity of the electrons. For constant temperature, as electron diffusivity depends on the gas pressure and subsequently on the gas density, the streamer radius r is inversely proportional to the gas density as in Equation (6), [8,11],
r = C + P ,
where P is the pressure and C + is a constant. The proportionality factor C + is determined under low pressure since the radial extension reduces with increasing pressure. Different values determined using the Schlieren technique are reported in previous studies; however, the uncertainty for these experimentally obtained values of C + are not provided. In [8], C + was estimated as 5 m Pa and in more recent work [32], a proportionality factor of 2 m Pa was used in the modelling of leader propagation in uniform background field. In [11], a confirmation for the value of C+ is provided for both strong and weak non-uniform fields at 50 kPa and 100 kPa. For strong non-uniform electric fields, the radius of the streamer is in the range between 41 ± 6 µm and 53 ± 6 µm at 50 kPa, and 17 ± 3 µm and 22 ± 3 µm at 100 kPa. In weak non-uniform electric fields, the radius of the streamer ranged between 48 ± 7 µm and 59 ± 5 µm at 50 kPa, and 21 ± 6 µm and 28 ± 5 µm at 100 kPa. Using Equation (6) and C+ = 2 (m·Pa) the expected radius of streamer at 50 kPa was 40 µm and 20 µm at 100 kPa. Taking only the average maximum values, the percentage error in strong non-uniform fields is 32.5% at 50 kPa and 10% at 100 kPa. For weak non-uniform fields, the percentage error is 47% at 50 kPa and 40% at 100 kPa for the average maximum radii. This is summarized in Appendix E.
Based on the numerical simulation of streamers in SF6 at 10 kPa, a similar error analysis has been conducted. For proportionality factor C+ = 2 (m·Pa), the expected streamer radius is 0.2 mm and for C+ = 5 (m·Pa), a radius of 0.5 mm is estimated. Both values are less than the 0.8 mm (for 5 kV stress) obtained in the current simulation yielding a 30% or 60% deviation depending on the proportionality factor. Even the experimental results demonstrate an increase in the deviation as the pressure decreases, which could be a possible explanation for the difference. Another possible reason for this discrepancy could be the background ionization term used in the model instead of photoionization. However, this investigation is beyond the scope of the present work as very limited information is available on the photoionization mechanisms in SF6. For more comparable values, the applied voltage can be varied. This is, however, a trade-off between the streamer velocity and its radius.

5. Conclusions

A comprehensive behaviour of streamer discharges in a highly electronegative gas under high voltage stress has been presented and compared to a possible ecological replacement especially in electro-technical applications.
Streamer discharges in SF6 at atmospheric pressure are difficult to model numerically because of the stepped leader transition leading to the complete breakdown. In the present work, the streamer initiation is observed however the high attachment rate in the streamer channel exceeding the ionization rate leads to the reduction and eventual interruption of the conduction channel due to decreasing electron density. When this occurred, the streamer radius decreases, it halts, and an increase in the electric field in the streamer head is observed. This could be a result of the limitations of the current computational model. At low pressures, however, the streamer discharge can be characterized using the drift-diffusion partial differential equations. This modelling has been conducted for a pressure of 10 kPa with varying voltages, and the quantitative parameters have been compared to that of air at the same pressure. The results show that the electronegativity of SF6 is particularly important for the cathode directed streamer discharges as it affects the electron density in the streamer channel with both increasing and decreasing gas densities. The number densities of ions in the streamer channel are approximately equivalent with a deviation occurring in the streamer head as shown in the spatial diagrams. It was found that streamers in SF6 are characterized by smaller radii in comparison to air. The streamer radius in air obtained in the present work shows a twofold increase as compared with that in SF6. However, the velocities of streamer propagating in air and SF6 at 10 kPa are comparable with each other; ~1 mm/ns for a 5 kV applied voltage. The results presented in this paper highlight the influence of increased electronegativity on streamer discharges in gases and provide a foundation for gas mixtures with SF6 reducing streamer initiation and propagation probabilities. In future studies, the possibility of using the fluid approach to quantify the streamer to leader transition and propagation of leaders in SF6 gas will be explored. Also, the influence of field utilization factor, η on streamer dynamics in SF6 will be studied.

Author Contributions

Conceptualization: F.B.-M., N.B. and R.H.; Methodology: F.B.-M., R.H. and I.N.; Writing—original draft preparation: F.B.-M., N.B., R.H. and I.T.; Writing—Reviewing and editing: F.B.-M., N.B. and I.T.; Supervision and project administration: N.B., R.H., I.N. and I.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the French National Research Agency in the framework of the “Investissements d’avenir” program (ANR-15-IDEX-02).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express appreciation to the eit InnoEnergy PhD school for supporting a mobility to the University of Strathclyde, Glasgow—UK during the preparation of this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Influence of Background Ionization Seed in Air

In most of the results obtained in this work, a background ionization of S b = 10 23   m 3 s 1 has been used. In the case of air specifically, S b replaces the photoionization mechanism that allows streamers to effectively propagate in the medium against the electron drift direction. Here, a comparison has been made between background ionization levels of 1015 m−3 s−1 and 1023 m−3 s−1 at atmospheric pressure. Figure A1 shows the 2D surface plot of the electric field at axial position 2 mm for both S b values.
For the streamer to reach an axial position of 2 mm from the point electrode, the lower background ionization level required a time of 3.6 ns as opposed to 3.2 ns for the high background ionization level. Also of note is the electric field in the head of the streamer with the background ionization level of 1015 m−3 s−1 having a maximum electric field of 154 kV/cm and background ionization level of 1023 m−3 s−1 having a maximum electric field of 129 kV/cm in the head of the streamer. The ease of propagation in the higher background ionization is further highlighted by the average velocity of the streamer, 1.2 mm/ns as compared to 1.11 mm/ns for the lower background ionization. The radius, however, is ~0.02 mm smaller with the lower background ionization.
Figure A1. 2D surface plot of electric field for S b = 10 15 m 3 s 1 (left) and S b = 10 23 m 3 s 1 (right).
Figure A1. 2D surface plot of electric field for S b = 10 15 m 3 s 1 (left) and S b = 10 23 m 3 s 1 (right).
J 05 00018 g0a1

Appendix B. Model Verification

To verify the model developed in COMSOL™ Multiphysics for streamer simulations, the propagation of a streamer in a 1 mm gap has been simulated and compared to simulation of a similar condition by the PRHE team in Laplace Laboratory. For these simulations, the conditions were an applied voltage of 3.5 kV, needle electrode radius of curvature of 0.1 mm, background ionization S b = 0   m 3 s 1 and an initial electron density of 1015 m−3.
Two test conditions are presented in the results below. The first being without adaptive mesh with the results for the electron density and electric field presented in Figure A2 and the second being with the inclusion of the adaptive mesh. Generally, the two streamer codes are consistent with each other however the COMSOL™ implementation shows a relatively lower adherence to the conservation of charge which is even more accentuated when the adaptive mesh refinement is used.
Figure A2. Electron density variation (left) and electric field distribution (right) on streamer axis at times t = 1, 2, and 3 ns for COMSOL ™ model (solid line) and reference work by PRHE team, Laplace lab [33] (dotted line). No background ionization and no adaptive mesh used.
Figure A2. Electron density variation (left) and electric field distribution (right) on streamer axis at times t = 1, 2, and 3 ns for COMSOL ™ model (solid line) and reference work by PRHE team, Laplace lab [33] (dotted line). No background ionization and no adaptive mesh used.
J 05 00018 g0a2
Figure A3. Electron density variation (left) and electric field distribution (right) on streamer axis at times t = 1, 2, and 3 ns for COMSOL™ model (solid line) and reference work by PRHE team, Laplace lab [33] (dotted line). Adaptive mesh used in COMSOL™ model.
Figure A3. Electron density variation (left) and electric field distribution (right) on streamer axis at times t = 1, 2, and 3 ns for COMSOL™ model (solid line) and reference work by PRHE team, Laplace lab [33] (dotted line). Adaptive mesh used in COMSOL™ model.
J 05 00018 g0a3

Appendix C. Cross Sectional Species and Swarm Parameters of SF6 Used in Simulations

Table A1. Cross Sectional Reactions.
Table A1. Cross Sectional Reactions.
Reaction TypeCross Sectional Reactions
Electron Impact Ionization e + SF 6 SF 6 + + 2 e
Electron Attachment e + SF 6 SF 6
Electron—Ion Recombination e + SF 6 + SF 6
Ion—Ion Recombination SF 6 + SF 6 + SF 6 + SF 6
Table A2. Swarm parameters.
Table A2. Swarm parameters.
ParameterEquationReference
Reduced   Electron   Diffusivity ,   D e N (1/m s) 3.553 · 10 2 ( | E ] / N ) 0.2424 × N
for   E N < 650   Td
[24,27]
Reduced   Electron   Mobility ,   μ e N (1/(m V s)) 1.027 · 10 19 ( | E ] / N ) 0.7424 × N E
for   10 < E N < 2000   Td
[24]
Positive   Ion   Mobility ,   μ p (m2/(V.s))) 6.0 · 10 5
For   E N < 120   Td
1.216 · 10 5 · ln ( | E ] / N ) + 5.89 · 10 4
For   120     E N < 350   Td
1.897 · 10 5 · ln ( | E ] / N ) 7.346 · 10 4
For   E N 350   Td
[24]
Negative   Ion   Mobility ,   μ n (m2/(V.s))) 1.69 · 10 32 ( | E ] / N ) 2 + 5.3 · 10 5
For   E N 500   Td
[24]
Reduced   Electron   Impact   Ionization ,   α N (m2) 3.4473 · 10 34 ( | E ] / N ) 2.985
for   E N 460   Td
11.269 · ( | E ] / N ) 1.159
for   E N > 460   Td
[24]
Reduced   Electron   Attachment ,   η N (m2) 2.0463 · 10 20 0.25379 ( | E ] / N ) + 1.4705 · 10 18 ( | E ] / N ) 2 3.0078 · 10 36 ( | E ] / N ) 3
For   50   <   E N     200   Td
7 · 10 21 exp ( 2.25 · 10 18 ( | E ] / N ) )
for   E N > 200   Td
[24]
Reduced Effective Ionization (m2) α η N
Electron Ion   Recombination ,   β e p (m3/s)1e13[34]
Ion Ion   Recombination ,   β p n (m3/s)1e12[35]

Appendix D. Cross Sectional Species and Swarm Parameters of Air Used in Simulations

Table A3. Cross Sectional Reactions.
Table A3. Cross Sectional Reactions.
Reaction TypeCross Sectional Reactions
Electron Impact Ionization e + A A + + 2 e
Electron Attachment e + A A
Electron—Ion Recombination e + A + A
Ion—Ion Recombination A + A + A + A
Table A4. Swarm Parameters.
Table A4. Swarm Parameters.
ParameterEquationReference
Reduced   Electron   Diffusivity ,   D e N (1/cm s) 0.3341 · 10 9 ( | E ] / N ) 0.54069 · μ e N [25]
Reduced   Electron   Mobility ,   μ e N (1/(cm V s)) E / E · ( 10 5.5236702 + 0.7822439 · l o g 10 E / N ) × N E
for   9.8   Td E N 1000   Td
E / E ( 10 5.8692884 + 0.4375671 · l o g 10 E / N ) × N E
for   E N < 9.8   Td
[26]
Positive   Ion   Mobility ,   μ p (m2/(V.s))) 2 e 4 [36]
Negative   Ion   Mobility ,   μ n (m2/(V.s))) 2.2 e 4 [36]
Reduced   Electron   Impact   Ionization ,   α N (cm2) 2 · 10 16 · exp 7.248 10 15 E / N
for   150   Td < E N
6.619 · 10 17 · exp 5.593 10 15 E / N
for   E N 150   Td
[25]
Reduced   Electron   Attachment ,   η N (cm2) 6.56041 · 10 19 1.45181 · 10 21 ( E / N ) + 1.45951 · 10 24 E / N 2 5.69565 · 10 28 E / N 3
for   600   Td E N < 1000   Td
6.23261 · 10 19 1.17646 · 10 21 ( E / N ) + 7.51103 · 10 25 E / N 2
for   170   Td E N < 600   Td
3.611 · 10 19 + 1.01192 · 10 20 ( E / N ) 3.17875 · 10 23 E / N 2
for   69   Td E N < 170   Td
3.10976 · 10 19 9.41213 · 10 21 ( E / N ) + 1.09693 · 10 22 E / N 2
for   23   Td E N < 69   Td
1.2409 · 10 19 + 8.9497 · 10 18 e x p ( | E | / N / 1.0931 ) + 1.3216 · 10 18 E / N / 6.05148 2
for   1   Td E N < 23   Td
[26]
Electron Ion   Recombination ,   β e p (m3/s) 2 e 13 [25]
Ion Ion   Recombination ,   β p n (m3/s) 2 e 13 [25]

Appendix E

Table A5. Summary of Experimental Radii Corresponding to Different Proportionality Factors in SF6 at 50 and 100 kPa.
Table A5. Summary of Experimental Radii Corresponding to Different Proportionality Factors in SF6 at 50 and 100 kPa.
C + Electric Field StrengthPressureRadius ExperimentalRadius ComputedReference
5 m Pa 100 kPa 50 µm[8]
2 m PaStrong non-uniform electric field50 kPa41 ± 6 µm → 53 ± 6 µm40 µm[11,32]
100 kPa17 ± 3 µm → 22 ± 320 µm
Weak non-uniform electric field50 kPa48 ± 7 µm → 59 ± 5 µm40 µm
100 kPa21 ± 6 µm → 28 ± 5 µm20 µm

References

  1. Christophorou, L.G.; Olthoff, J.K.; Green, D.S. Gases for Electrical Insulation and Arc Interruption: Possible Present and Future Alternatives to Pure SF6, Technical Note (NIST TN)-1425; NIST: Gaithersburg, MD, USA, 1997. [Google Scholar]
  2. Beroual, A.; Haddad, A.M. Recent Advances in the Quest for a New Insulation Gas with a Low Impact on the Environment to Replace Sulfur Hexafluoride (SF6) Gas in High-Voltage Power Network Applications. Energies 2017, 10, 1216. [Google Scholar] [CrossRef] [Green Version]
  3. Stewart, G.; Wilson, M.P.; Timoshkin, I.V.; Given, M.J.; MacGregor, S.J.; Sinclair, M.A.; Thomas, K.J. The suitability of N2 to replace SF6 in a triggered spark-gap switch for pulsed power applications. In Proceedings of the 44th International Universities Power Engineering Conference (UPEC), Glasgow, UK, 1–4 September 2009; pp. 1–5. Available online: https://0-ieeexplore-ieee-org.brum.beds.ac.uk/document/5429548 (accessed on 1 March 2022).
  4. McKee, R.G. Pulsed Power and SF6; NNSA: Washington, DC, USA, 2016. Available online: https://www.osti.gov/servlets/purl/1334943 (accessed on 22 December 2021). [CrossRef]
  5. Larsson, A. Gas-discharge closing switches and their time jitter. IEEE Trans. Plasma Sci. 2012, 40, 2431–2442. [Google Scholar] [CrossRef]
  6. Gallimberti, I.; Wiegart, N. Streamer and leader formation in SF6 and SF6 mixtures under positive impulse conditions. II. Streamer to leader transition. J. Phys. D Appl. Phys. 1986, 19, 2363. [Google Scholar] [CrossRef]
  7. Gallimberti, I. Breakdown mechanisms in electronegative gases. In Gaseous Dielectrics V, Proceedings of the Fifth International Symposium on Gaseous Dielectrics, Knoxville, TN, USA, 3–7 May 1987; Christophorou, L., Bouldin, D., Eds.; Pergamon Press: Oxford, UK, 1987; pp. 61–80. [Google Scholar]
  8. Niemeyer, L.; Ullrich, L.; Wiegart, N. The mechanism of leader breakdown in electronegative gases. IEEE Trans. Electr. Insul. 1989, 24, 309–324. [Google Scholar] [CrossRef]
  9. Niemeyer, L. Leader breakdown in compressed SF6: Recent concepts and understanding. In Gaseous Dielectrics VI, Proceedings of the Sixs International Symposium on Gaseous Dielectrics, Knoxville, TN, USA, 23–27 September 1990; Christophorou, L., Sauers, I., Eds.; Plenum Press: New York, NY, USA, 1991; pp. 49–60. [Google Scholar]
  10. Chalmers, I.D.; Gallimberti, I.; Gibert, A.; Farish, O. The development of electrical leader discharges in a point-plane gap in SF6. Proc. R. Soc. Lond. A Math. Phys. Sci. 1987, 412, 285–308. [Google Scholar]
  11. Bujotzek, M.; Seeger, M.; Schmidt, F.; Koch, M.; Franck, C. Experimental investigation of streamer radius and length in SF6. J. Phys. D Appl. Phys. 2015, 48, 245201. [Google Scholar] [CrossRef] [Green Version]
  12. Morrow, R. Theory of electrical corona in SF6. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrom. Detect. Assoc. Equip. 1996, 382, 57–65. [Google Scholar] [CrossRef]
  13. Sato, N. Discharge current induced by the motion of charged particles. J. Phys. D Appl. Phys. 1980, 13, L3. [Google Scholar] [CrossRef]
  14. Dhali, S.K.; Pal, A.K. Numerical simulation of streamers in SF6. J. Appl. Phys. 1998, 63, 1355–1362. [Google Scholar] [CrossRef]
  15. Babaeva, N.Y.; Naidis, G.V. Simulation of stepped propagation of positive streamers in SF6. J. Phys. D Appl. Phys. 2001, 35, 132. [Google Scholar] [CrossRef]
  16. Fuliang, L.I.; Feng, W.; Guoli, W.; Pfeiffer, W.; Rongtao, H. Study of Formation and Propagation of Streamers in SF6 and Its Gas Mixtures with Low Content of SF6 Using a One-Dimensional Fluid Model. Plasma Sci. Technol. 2012, 14, 187. [Google Scholar]
  17. Wang, L.; Ou, X.; Zheng, Y.; Liu, J.; Lin, X.; Zhang, T. Numerical simulations of the SF6-N2 mixed gas streamer discharge development process. AIP Adv. 2019, 9, 055320. [Google Scholar] [CrossRef] [Green Version]
  18. Luo, B.; He, H.; Cheng, C.; Wu, Y.; Huang, Y.; Chen, W. Simulation of the positive streamer propagation in SF6/N2 mixtures with different content of SF6. In Proceedings of the IEEE International Conference on High Voltage Engineering and Application (ICHVE), Beijing, China, 6–10 September 2020. [Google Scholar]
  19. Singh, S. Computational Framework for Studying Charge Transport in High-Voltage Gas-Insulated Systems; Chalmers University of Technology: Gothenburg, Sweden, 2017. [Google Scholar]
  20. Quast, M.; Lalic, N.R. Streamer Propagation in a Point-to-Plane Geometry. Excerpt from the Proceedings of the COMSOL Conference. 2009. Available online: https://www.comsol.fr/paper/download/101137/Quast.pdf (accessed on 17 March 2018).
  21. Wormeester, G.; Pancheshnyi, S.; Luque, A.; Nijdam, S.; Ebert, U. Probing photo-ionization: Simulations of positive streamers in varying N2:O2 mixtures. J. Phys. D Appl. Phys. 2010, 43, 505201. [Google Scholar] [CrossRef] [Green Version]
  22. Bagheri, B.; Teunissen, J.; Ebert, U.; Becker, M.M.; Chen, S.; Ducasse, O.; Eichwald, O.; Loffhagen, D.; Luque, A.; Mihailova, D.; et al. Comparison of six simulation codes for positive streamers in air. Plasma Sources Sci. Technol. 2018, 27, 095002. [Google Scholar] [CrossRef] [Green Version]
  23. Bagheri, B.; Teunissen, J.; Ebert, U. Simulation of positive streamers in CO2 and in air: The role of photoionization or other electron sources. Plasma Sources Sci. Technol. 2020, 29, 125021. [Google Scholar] [CrossRef]
  24. Morrow, R. A survey of the electron and ion transport properties of SF6. IEEE Trans. Plasma Sci. 1986, 14, 234–239. [Google Scholar] [CrossRef]
  25. Morrow, R.; Lowke, J.J. Streamer propagation in air. J. Phys. D Appl. Phys. 1997, 30, 614. [Google Scholar] [CrossRef] [Green Version]
  26. Ducasse, O.; Papageorghiou, L.; Eichwald, O.; Spyrou, N.; Yousfi, M. Critical analysis on two-dimensional point-to-plane streamer simulations using the finite element and finite volume methods. IEEE Trans. Plasma Sci. 2007, 35, 1287–1300. [Google Scholar] [CrossRef]
  27. Liu, T.; Timoshkin, I.V.; MacGregor, S.J.; Wilson, M.P.; Given, M.J.; Bonifaci, N.; Hanna, R. Field-Time Breakdown Characteristics of Air, N2, CO2, and SF6. IEEE Trans. Plasma Sci. 2020, 48, 3321–3331. [Google Scholar] [CrossRef]
  28. COMSOL 5.4, Plasma Module User Guide. 2018. Available online: https://doc.comsol.com/5.4/doc/com.comsol.help.plasma/PlasmaModuleUsersGuide.pdf (accessed on 17 March 2020).
  29. Farish, O. Breakdown in SF6 and its Mixtures in Uniform and Nonuniform Fields. In Proceedings of the International Symposium on Gaseous Dielectrics, Knoxville, TN, USA, 6–8 March 1978; pp. 60–83. [Google Scholar]
  30. Dupuy, J.; Gilbert, A. Comparison of point-to-plane discharges in air and SF6. J. Phys. D Appl. Phys. 1982, 15, 655. [Google Scholar] [CrossRef]
  31. Soulié, S. Study of the Dielectric Properties of HFO Gas, and Its Application to Reduce the Environmental Impact of Medium-Voltage Systems; Université de Grenoble: Grenoble, France, 2021. [Google Scholar]
  32. Seeger, M.; Niemeyer, L.; Bujotzek, M. Leader propagation in uniform background fields in SF6. J. Phys. D Appl. Phys. 2009, 42, 185205. [Google Scholar] [CrossRef]
  33. Eichwald, O.; Bensaad, H.; Ducasse, O.; Yousfi, M. Effects of numerical and physical anisotropic diffusion on branching phenomena of negative-streamer dynamics. J. Phys. D Appl. Phys. 2012, 45, 385203. [Google Scholar] [CrossRef]
  34. Gleizes, A.; Habib, A.A.M.; Vacquié, S. Calculation of electron-ion recombination in an SF6 arc plasma. J. Phys. D Appl. Phys. 1989, 22, 1464. [Google Scholar] [CrossRef]
  35. Christophorou, L.G.; Olthoff, J.K. Electron interactions with SF6. J. Phys. Chem. Ref. Data 2000, 29, 267–330. [Google Scholar] [CrossRef] [Green Version]
  36. Serdyuk, Y. Propagation of Cathode-Directed Streamer Discharges in Air. 2013. Available online: https://www.comsol.com/paper/download/181801/serdyuk_paper.pdf (accessed on 17 March 2018).
Figure 1. Reduced Effective Ionization as a function of Reduced Electric Field (E/N) for SF6 [24] and air [25,26].
Figure 1. Reduced Effective Ionization as a function of Reduced Electric Field (E/N) for SF6 [24] and air [25,26].
J 05 00018 g001
Figure 2. Reduced Electron Mobility as a function of Reduced Electric Field (E/N) for SF6 [24] and air [26].
Figure 2. Reduced Electron Mobility as a function of Reduced Electric Field (E/N) for SF6 [24] and air [26].
J 05 00018 g002
Figure 3. Reduced Electron Diffusivity as a function of Reduced Electric Field (E/N) for SF6 [24,27] and air [25].
Figure 3. Reduced Electron Diffusivity as a function of Reduced Electric Field (E/N) for SF6 [24,27] and air [25].
J 05 00018 g003
Figure 4. Mesh Sizes and Refinement.
Figure 4. Mesh Sizes and Refinement.
J 05 00018 g004
Figure 5. 2D surface plot of electron density profile at times 0.5, 1, and 2 ns in air at atmospheric pressure. Gap distance 5 mm, Sb = 1023 m−3 s−1. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 5. 2D surface plot of electron density profile at times 0.5, 1, and 2 ns in air at atmospheric pressure. Gap distance 5 mm, Sb = 1023 m−3 s−1. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g005
Figure 6. Surface plot of electric field distribution (kV/cm) corresponding to the electron density profile in Figure 5.
Figure 6. Surface plot of electric field distribution (kV/cm) corresponding to the electron density profile in Figure 5.
J 05 00018 g006
Figure 7. Surface plot of positive ion density profile (left) and negative ion density profile (right) at time 2 ns in air at atmospheric pressure. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 7. Surface plot of positive ion density profile (left) and negative ion density profile (right) at time 2 ns in air at atmospheric pressure. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g007
Figure 8. 2D surface plot of electron density profile at times 0.5, 1 and 2 ns in SF6 at atmospheric pressure. Gap distance 5 mm, Sb = 1023 m−3 s−1. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 8. 2D surface plot of electron density profile at times 0.5, 1 and 2 ns in SF6 at atmospheric pressure. Gap distance 5 mm, Sb = 1023 m−3 s−1. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g008
Figure 9. Surface plot of electric field distribution (kV/cm) corresponding to the electron density profile in Figure 8.
Figure 9. Surface plot of electric field distribution (kV/cm) corresponding to the electron density profile in Figure 8.
J 05 00018 g009
Figure 10. Surface plot of positive ion density profile (left) and negative ion density profile (right) at time 2 ns in SF6 at atmospheric pressure. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 10. Surface plot of positive ion density profile (left) and negative ion density profile (right) at time 2 ns in SF6 at atmospheric pressure. Scaling parameter: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g010
Figure 11. Photo (experiments) [31] of leader discharges in SF6 at 1 bar in point to plane configuration. Needle radius of curvature = 0.2 mm, gap distance = 50 mm. Applied voltage = 107 kV.
Figure 11. Photo (experiments) [31] of leader discharges in SF6 at 1 bar in point to plane configuration. Needle radius of curvature = 0.2 mm, gap distance = 50 mm. Applied voltage = 107 kV.
J 05 00018 g011
Figure 12. Line diagrams along the symmetry axis of the log of the Electron Density profile (left) and its corresponding Electric Field distribution (right) for a streamer propagating in a 5 mm air gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 12. Line diagrams along the symmetry axis of the log of the Electron Density profile (left) and its corresponding Electric Field distribution (right) for a streamer propagating in a 5 mm air gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g012
Figure 13. Line diagrams along the symmetry axis of the log of the Positive Ion Density profile (left) and the Negative Ion Density profile (right) for a streamer propagating in a 5 mm air gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 13. Line diagrams along the symmetry axis of the log of the Positive Ion Density profile (left) and the Negative Ion Density profile (right) for a streamer propagating in a 5 mm air gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g013
Figure 14. Line diagrams along the symmetry axis of the log of the Electron Density profile (left) and its corresponding Electric Field distribution (right) for a streamer propagating in a 5 mm SF6 gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 14. Line diagrams along the symmetry axis of the log of the Electron Density profile (left) and its corresponding Electric Field distribution (right) for a streamer propagating in a 5 mm SF6 gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g014
Figure 15. Line diagrams along the symmetry axis of the log of the Positive Ion Density profile (left) and the Negative Ion Density profile (right) for a streamer propagating in a 5 mm SF6 gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 15. Line diagrams along the symmetry axis of the log of the Positive Ion Density profile (left) and the Negative Ion Density profile (right) for a streamer propagating in a 5 mm SF6 gap at a pressure of 10 kPa. Applied voltage of 5 kV. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g015
Figure 16. Photo (experiments) [31] of streamer discharges in SF6 at 0.1 bar in point to plane configuration. Needle radius of curvature = 0.2 mm, gap distance = 50 mm. Applied voltage = 37 kV.
Figure 16. Photo (experiments) [31] of streamer discharges in SF6 at 0.1 bar in point to plane configuration. Needle radius of curvature = 0.2 mm, gap distance = 50 mm. Applied voltage = 37 kV.
J 05 00018 g016
Figure 17. Line diagram of the positive ion density (blue dash-point lines) and electron density (solid gray lines) in SF6 at 10 kPa when a voltage of 5 kV (left) and 10 kV (right) is applied. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 17. Line diagram of the positive ion density (blue dash-point lines) and electron density (solid gray lines) in SF6 at 10 kPa when a voltage of 5 kV (left) and 10 kV (right) is applied. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g017
Figure 18. 1D line graph of the electric field for a streamer in SF6 at 10 kPa when a voltage of 10 kV is applied.
Figure 18. 1D line graph of the electric field for a streamer in SF6 at 10 kPa when a voltage of 10 kV is applied.
J 05 00018 g018
Figure 19. Line diagram along the symmetry axis of the log of the Electron Density profile for an applied voltage of 2 kV at pressure 10 kPa. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
Figure 19. Line diagram along the symmetry axis of the log of the Electron Density profile for an applied voltage of 2 kV at pressure 10 kPa. Scaling parameter for density profiles: 20 on logarithmic plot represents a density value of 1020 m−3.
J 05 00018 g019
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boakye-Mensah, F.; Bonifaci, N.; Hanna, R.; Niyonzima, I.; Timoshkin, I. Modelling of Positive Streamers in SF6 Gas under Non-Uniform Electric Field Conditions: Effect of Electronegativity on Streamer Discharges. J 2022, 5, 255-276. https://0-doi-org.brum.beds.ac.uk/10.3390/j5020018

AMA Style

Boakye-Mensah F, Bonifaci N, Hanna R, Niyonzima I, Timoshkin I. Modelling of Positive Streamers in SF6 Gas under Non-Uniform Electric Field Conditions: Effect of Electronegativity on Streamer Discharges. J. 2022; 5(2):255-276. https://0-doi-org.brum.beds.ac.uk/10.3390/j5020018

Chicago/Turabian Style

Boakye-Mensah, Francis, Nelly Bonifaci, Rachelle Hanna, Innocent Niyonzima, and Igor Timoshkin. 2022. "Modelling of Positive Streamers in SF6 Gas under Non-Uniform Electric Field Conditions: Effect of Electronegativity on Streamer Discharges" J 5, no. 2: 255-276. https://0-doi-org.brum.beds.ac.uk/10.3390/j5020018

Article Metrics

Back to TopTop