Prediction of Insulation Capability for Ground Fault to Consider Asymmetry in SF6 Circuit Breaker
Prediction of Insulation Capability for Ground Fault to Consider Asymmetry in SF6 Circuit Breaker
Journal of Electrical Engineering and Technology. 2015. Sep, 10(5): 2046-2051
Copyright © 2015, The Korean Institute of Electrical Engineers
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : June 10, 2014
  • Accepted : May 13, 2015
  • Published : September 01, 2015
Export by style
Cited by
About the Authors
Yeon-Ho, Oh
Corresponding Author: Korea Electrotechnology Research Institute, Korea. (
Ki-Dong, Song
Korea Electrotechnology Research Institute, Korea. ({kdsong, kimhk}
Hong-Kyu, Kim
Korea Electrotechnology Research Institute, Korea. ({kdsong, kimhk}
Hae June, Lee
Dept. of Electrical Engineering, Pusan University, Korea. (
Sung-Chin, Hahn
Dept. of Electrical Engineering, Dong-A University, Korea. (

Currently, most high-voltage gas circuit breakers (CBs) include asymmetrical geometries in the shield, the tank, the hot-gas exhaust, and the connection parts for bushings. For this reason, a 3-dimensional (3-D) analysis of the insulation capability is necessary, rather than a 2-D analysis. However, a 3-D analysis has difficulties due to the computational time and complex modeling. This paper presents a 3-D analysis considering the asymmetry in high-voltage gas CBs and a technique to reduce the calculation time. In the proposed technique, the arc plasma requiring the most computational time is first calculated by a 2-D analysis. Then, the results such as pressure, temperature, and velocity are input as a source for the 3-D analysis. This technique is applied to a 145kV self-blast-type CB and the analysis result exhibits good agreement with the experimental result.
1. Introduction
When an arc plasma is extinguished just after a fault current interruption at current zero (CZ) of alternating current, hot gas, which is several thousand degrees, fills the inside of the gas CB. Because of the rapid reduction in the insulation capability due to the hot gas, the probability of a ground fault is increased during the transient recovery voltage (TRV) after the interruption.
Thus, the insulation capability for a ground fault must be considered in the dielectric design of a high-voltage CB [1] , [2] . A 3-D analysis should be performed to estimate the dielectric characteristic after current interruption if the CB has an asymmetric region such as a shield, a tank, a hot-gas exhaust, or the connection parts of a bushing, as shown in Fig. 1 .
PPT Slide
Lager Image
Schematic diagram for a gas CB. The shield, exhaust, and connection parts are asymmetrical parts in general.
In general, the calculation time of the 3-D analysis of an arc plasma depends on the accuracy of the physical arc models. The physical phenomena in a gas CB, which are the radiative heat transfer, turbulence effects, electromagnetic field, metal vapor, nozzle ablation, and so on, are mainly focused on the arc physics [1] . However, the 3-D analysis has difficulty simulating these arc phenomena because of the significant computational time compared to a 2-D analysis.
This paper presents a novel technique to reduce the computational time of a 3-D analysis.
In this study, the calculation of the arc plasma requires the most computational time and is first conducted by a 2-D axisymmetrical analysis. These 2-D analysis results such as pressure, temperature, and velocity are used as an input source for the 3-D analysis. Finally, the 3-D analysis of the insulation capability is calculated for the entire domain, except the arc plasma region. For the effects of turbulent flow, the standard κ–ε model is considered. The commercial computational fluid dynamic (CFD) code ANSYS-Fluent is used for simulation of these analyses.
In order to verify the proposed technique, experimental data for a 145kV self-blast-type CB is used. The analysis result exhibits good agreement with the experimental result.
2. Transient Behavior of the Thermal Fluid Dynamics
- 2.1 Governing equations and arc model
In a gas CB, a switching arc is generated when the arcing contacts are separated, heating the surrounding insulation gas such as SF 6 . To extinguish the arc, hot gas between the arcing contacts should be cooled down to recover the insulation capability. For this purpose, the gas inside the puffer cylinder is compressed mechanically and blasted into the arc to exhaust the hot gas. In order to simulate the transient thermal behavior, the commercial CFD code ANSYS-Fluent is used with a combination of arc physics.
The thermal plasma can be treated as a one-temperature fluid model because of the assumption of local thermodynamic equilibrium [3] .
The governing equations for analyzing the transient thermal fluid dynamics in the CB are represented by the mass, momentum, and energy conservation:
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where ρ is the density,
PPT Slide
Lager Image
is the velocity, Sm is the mass source, p is the pressure,
PPT Slide
Lager Image
is the stress tensor, μ is the molecular viscosity, I is the unit tensor, h is the enthalpy, k is the thermal conductivity, Se is the energy source, σ is the electrical conductivity, E is the electric field strength, R is the radiation loss, and T and r are the gas temperature and arc radius, respectively.
To consider the compressible effect and the turbulent model included in the real configuration of the nozzle and exhaust, the standard κ–ε turbulence model is adopted in this work [4] .
The arc energy is considered to be ohmic heating equal to σE 2 . With a finite element method (FEM), which was developed by using the user-defined-function of ANSYS-Fluent, E was obtained in this study. In this work, the radiative heat transfer is treated as a loss of arc energy. The radiation loss can be calculated by using empirical arc model as function of both arc temperature and arc radius. In this model, the emitted energy from arc is 90% of the total arc energy and 80% of the emitted energy is reabsorbed in the arc edge [1] . This simplified model can save the computational time significantly.
The deviation of the insulation gas from ideal behavior was also considered by using the thermodynamic and transport properties of the gas. Fig. 2 shows the thermodynamic conductivity of a insulation medium of SF 6 ; further details are presented in [5] .
PPT Slide
Lager Image
Thermal conductivity of SF6 gas according to temperature.
- 2.2 Self-Blast-type circuit breaker and nozzle ablation
Fig. 1 shows a schematic diagram of the geometry of the 145kV self-blast-type CB, which consists of two almost symmetrical parts separated by the throat of the back-flow channel. In a self-blast-type CB, the nozzle ablation plays an important role in the pressure increase in the heating volume, which strongly affects the interruption capability. During the high-current phase, the inten0se radiation of the high-temperature arc leads to ablation of the nozzle surfaces. This ablated material flows back into the heating volume during the arc, clogging the nozzle throat. The backflow of the ablation gas creates a very high pressure in the heating volume and blows hot gas between the arcing contacts.
In this study, a simplified model based on the two-zone model is suggested by the following procedure [6] .
1) The total flux of ablated mass T , in accordance with the arc energy, is calculated by
PPT Slide
Lager Image
where δ is the total ablation factor and has a value of approximately 17 mg/kJ, U is the arc voltage, and i is the arc current.
2) The ablated materials, which are composed mainly of PTFE (C 2 F 4 ), are mixed with cold SF 6 gas. The PTFE mass fraction can be calculated by an additional mass-fraction equation:
PPT Slide
Lager Image
where αPTFE is the PTFE mass fraction, and Γ is the diffusion coefficient of PTFE, which is approximately 10 −4 kg m −1 s −1 .
The thermodynamic and transport properties of a mixture of these gases are given in [5] .
3. Calculation of the Dielectric Characteristics
- 3.1 Critical electric field,Ecrit
The insulation capability of a gas CB under hot-gas conditions can be quantified by the ratio of the critical electric-field intensity Ecrit to the electric-field intensity for the TRV Etrv . Ecrit is the dielectric strength of the SF 6 gas in accordance with the temperature and pressure. Breakdown can occur when Etrv is greater than Ecrit .
Fig. 3 shows Ecrit as a function of both temperature and pressure. Ecrit has a non-linear characteristic above approximately 800 K because the thermal dissociation and ionization of SF 6 begin to be generated at this temperature. These thermal properties, which are characterized by drastic dissociation and ionization starting at 800 K, especially at 1,500-2,500 K, affect the thermal conductivity, as shown in Fig. 2 .
PPT Slide
Lager Image
Critical electric-field intensity for SF6 gas as function of both temperature and pressure.
It is noted that Ecrit strongly depends on temperature but weakly depends on pressure above 2,000 K. This indicates that the probability of dielectric breakdown is very high above 2,000 K at the peak of the TRV with almost no dependence on pressure. Uchii [7] assumed that Ecrit is practically zero at 3,000 K.
For Ecrit used in this work, PTFE ablated from the nozzle is not considered. If the PTFE properties are considered, Ecrit might change. However, Ecrit of pure SF 6 gas is practically used to evaluate the insulation capability for a ground fault because there is small amount of PTFE vapor between the exhaust and the tank. For the same reason, the amount of variation in Ecrit due to erosion of the arcing contacts is small and negligible.
- 3.2 Algorithm for the evaluation of the dielectric characteristics
The entire computational domain can be divided into 2-D axisymmetrical and 3-D asymmetrical analyses, as shown in Fig. 4 . The arc plasma region, puffer cylinder, and heating volume are included in the 2-D analysis, which are completely axisymmetrical geometries, and the 3-D analysis contains the other asymmetrical regions. The 2-D analyzed results of pressure, temperature, velocity, and density are used as the 3-D analysis input data.
PPT Slide
Lager Image
Entire computational domain for the 2-D axisymmetrical and 3-D asymmetrical analysis geometries.
Fig. 5 shows the flowchart for the prediction algorithms of the dielectric strength under hot-gas conditions.
PPT Slide
Lager Image
Flowchart for the prediction algorithms of the dielectric strength.
The physical phenomena of the arc plasma can be calculated via a 2-D axisymmetrical analysis including nozzle ablation. The hot gas generated from the arc plasma flows toward the nozzle exhaust, including pressure, velocity, density and temperature components.
These 2-D output variables then are converted to 3-D input variables, and the 3-D thermal behavior and insulation capability are calculated, except the arc plasma.
- 3.3 Inlet boundary and conditions for the 3-D Analysis
Fig. 6 shows the 2-D output and 3-D inlet boundaries in the real CB geometry. The pressure P ( t ) 2d from the 2-D analysis is used as the inlet source of the 3-D analysis with the temperature T ( t ) 2d , velocity v ( t ) 2d , and density ρ ( t ) 2d :
PPT Slide
Lager Image
PPT Slide
Lager Image
where P ( t ) s2d is the static pressure, and P ( t ) d2d is the dynamic pressure from the 2-D analysis. In this study, the entire computational domain except for the nozzle area is assumed to be a subsonic flow with a flow-speed Mach number less than one.
For conversion from the 2-D analysis results to the source of the 3-D analysis, there is an absence in the input dimensions due to the expansion of dimensions, namely, the z-axis components do not exist. For the 2-D analysis, the arc plasma between the arcing contacts can be fully preformed in an axisymmetrical analysis. In the nozzle exhaust, the hot gas from the arc collides with the arcing contacts and shields. Then, there could be z-axis components outside of the nozzle.
The axial and radial components in the 2-D axisymmetrical analysis are used as the sources of x- and y-axis components, respectively, in the 3-D analysis. However, zeroes are used as the inputs for the z components in the 3-D analysis owing to the absence of z-axis components from the 2-D analysis.
Thus, the inlet boundary in 3-D should be the surface on which the z-components of the output variables, pressure and velocity, are minimized. In this study, the inlet boundary in 3-D is the surface near the nozzle exhaust in Fig. 6 . These z-components in the inlet surface are expected to be negligible because of the significantly strong axial velocity component. This assumption could be validated by a comparison with experimental results.
PPT Slide
Lager Image
Output variables and inlet boundary for the 3-D analysis.
4. Comparison with Experimental Results
- 4.1 145-kV Self-blast-type gas circuit breaker
Fig. 7 shows the 145kV self-blast CB installed in the testing room. There are the asymmetric parts inside for the 145kV self-blast CB as shown in Fig. 8 . The amplitude of the testing current is 40kA rms with a frequency of 60Hz. The arcing time, which is the time after the separation of the arcing contacts to CZ, is 20.3ms. As a result of the test a dielectric breakdown occur between a shield and tank at 73.6μs after CZ.
PPT Slide
Lager Image
Installation of a 145kV self-blast gas CB in the testing room.
PPT Slide
Lager Image
Asymmetric parts for a 145kV self-blast CB.
In general, the pressure increase inside a self-blast gas CB is greater than a puffer-type CB because the hot gas from the arc plasma flows backward, heating the volume by nozzle ablation. This heats the gas and increases the pressure inside the heating volume. Thus, there would be many weak points in the dielectric capability of a self-blast gas CB owing to the large amount of hot gas.
In a self-blast CB, nozzle ablation considerably affects the pressure distribution during the high-current phase and also the movement of hot gas. Considering nozzle ablation and the arc energy the movement of the contacts are simulated.
As a result of simulation contours of the hot gas before and after CZ are shown in Fig. 9 . After CZ, the temperature of the hot gas is greater than 2000 K, and there is a high probability of breakdown in the exhaust. Fig. 10 (a) and (b) show the contours of Etrv/Ecrit and the location of the breakdown after the test respectively. The test results show that there is a dielectric breakdown at t = 73.6μs after CZ, and the location of the breakdown is shown in Fig. 10 (b) . The location having a higher probability of breakdown, as predicted by the analysis, is similar to that of the test results. The experimental results demonstrate the reliability of the proposed technique.
PPT Slide
Lager Image
Contours of hot gas greater than 2000K before and after CZ for a 145kV self-blast gas CB.
PPT Slide
Lager Image
Contours of Etrv/Ecrit and the location of the breakdown.
The computational time including the modeling time is reduced to 48 h in comparison with 220 h for the full 3-D analysis
It should be emphasized that Ecrit is affected by temperature rather than pressure. Thus, the design to adequately pass the hot gas or to extend the passage of hot gas from the exhaust to the outlet is necessary to decrease the temperature of the hot gas from the exhaust.
5. Conclusion
This paper presents a technique to characterize the insulation capability for a ground fault considering the asymmetry in a gas CB. The results of this work are summarized as follows.
  • 1) In order to consider the asymmetric configuration of a CB, a 3-D analysis based on the 2-D analysis results is performed.
  • 2) To reduce the calculation time, a 2-D axisymmetrical analysis is conducted for the physical phenomena of the arc plasma, which requires considerable computational time.
  • 3) The output variables from the 2-D analysis are used as the source of the 3-D analysis as function of time. In the 3-D analysis, Etrv/Ecritis obtained as a function of pressure and temperature, and the characteristics of the dielectric strength are evaluated.
  • 4) The simulation results are compared with experimental results and demonstrate the reliability of the proposed technique.
It is noted that the total computational time is drastically reduced in comparison with a full 3-D analysis. However, this technique does not consider the z-axis components of the inlet boundary in the 3-D analysis. Thus, the proposed technique is only applied to a gas CB that has a strong axial flow.
Yeon-Ho Oh He received the B.S. and M.S. degrees in electrical engineering from Dong-A University, Busan, Korea, in 1991 and 1993, respectively. He is currently a Research Engineer with the Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute (KERI).
Ki-Dong Song He received the B.S. and M.S. degrees in electrical engin- eering from InHa University, Incheon, Korea, in 1988 and 1990, respectively. He got his Ph.D. degree in electrical engineering from Busan National Uni- versity in 2003. He worked for Korea Electrotechnology Research Institute (KERI) from 1990 and is now the director of the Power Apparatus Research Center, HVDC Research Division. His research interests focus on analysis and design for AC Power Devices and HVDC circuit breakers.
Hong-Kyu Kim He was born in Yeongchon, Korea, in 1969. He received his B. Sc., M. Sc. and Ph.D. degrees in 1995, 1997 and 2001, respectively, from the School of Electrical Engineering, Seoul National University, Korea. He has been working for the Korea Electrotechnology Research Institute (KERI) since 2001. His research field covers the numerical analysis of electromagnetic fields, and the analysis of cold and hot gas flow in high voltage gas circuit breakers. He is a regular member of the CIGRE working group A3.24.
Hae June Lee Professor Hae June Lee received his B.S. degree with from the Department of Nuclear Engineering at the Seoul National University in 1994, and his M.S. and Ph.D. degrees in physics in 1996 and 1998 from POSTECH. He worked as a post doc. in the Department of Electrical Engi- neering at the UC Berkeley in 2000 and 2001 and as a research scientist in Korea Electrotechnology Research Institute (KERI) after then. Since 2004, Prof. Lee has been a faculty member in the Department of Electrical Engineering at Pusan National University (PNU), South Korea.
Sung-Chin Hahn He received the B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National Uni- versity, Seoul, Korea, in 1979, 1981, and 1992, respectively. He is currently a Professor with the Department of Electrical Engineering, Dong-A Uni- versity, Busan, Korea. His current research covers multi-physics analysis and optimal design of power apparatus and electric machines.
Kim Hong-Kyu , Chong Jin-Kyo , Song Ki-Dong 2010 “Analysis of Dielecric Breakdown of Hot SF6 Gas in a Gas Circuit Breaker,” JEET 5 (2) 264 - 269
Takana H. , Uchii T. , Kawano H. , Nishiyama H. 2007 “Real-Time Numerical Analysis on Insulation Capabiity Improvement of Compact Gas Circuit Breaker,” IEEE Trans. Power Delivery 22 (3) 1541 - 6    DOI : 10.1109/TPWRD.2007.899796
Yousfi M. , Robin-Jouan Ph. , Kanzari Z. 2005 “Breakdown Electric Field Calculations of Hot SF6 for High Voltage Circuit Breaker Applications,” IEEE Trans. on Dielect. Elect. Insulation 12 (6) 1192 - 1200
Song Ki-Dong , Lee Byeong-Yoon , Park Kyong-Yop 2004 “Calculation of the Post-Arc Current in a Supersonic Nozzle by Using the κ-ε Turbulence Model,” Journal of the Korean Physical Society 45 (6) 1537 - 43
Zhang Jin Ling , Yan Jiu Dun , Murphy Anthony B. , Hall W. , Fang M.T.C. 2002 “Computational Investigation of Arc Behavior in an Auto-Expansion Circuit Breaker Contaminated by Ablated Nozzle Vapor,” IEEE Trans. on Plasma Science 30 (2) 706 - 19    DOI : 10.1109/TPS.2002.1024273
Claessens M , Möller K , Thiel H G 1997 “A computational fluid dynamics simulation of high- and low-current arcs in self-blast circuit breakers,” J. Phys. D: Appl. Phys. 30 1899 - 1907    DOI : 10.1088/0022-3727/30/13/011
Uchii T. , Nakamoto T. , Nishiwaki S. , Toyoda M. , Boggs S. A. 2004 “Optimization of Dead Tank Gas Circuit Breaker Design Based on Quantification of Hot Gas Flow and Dielecgtric Properties,” IEEE Trans. Power Delivery 19 (1) 181 - 85    DOI : 10.1109/TPWRD.2003.820217
Müller L 1993 “Modelling of an ablation controlled arc,” J. Phys. D:Appl. Phys. 25 1253 - 59