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.
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
. 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
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
. 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
. 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
The governing equations for analyzing the transient thermal fluid dynamics in the CB are represented by the mass, momentum, and energy conservation:
is the density,
is the velocity,
is the mass source,
is the pressure,
is the stress tensor,
is the molecular viscosity,
is the unit tensor,
is the enthalpy,
is the thermal conductivity,
is the energy source,
is the electrical conductivity,
is the electric field strength,
is the radiation loss, and
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
The arc energy is considered to be ohmic heating equal to
. With a finite element method (FEM), which was developed by using the user-defined-function of ANSYS-Fluent,
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
. 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.
shows the thermodynamic conductivity of a insulation medium of SF
; further details are presented in
Thermal conductivity of SF6 gas according to temperature.
- 2.2 Self-Blast-type circuit breaker and nozzle ablation
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
1) The total flux of ablated mass
, in accordance with the arc energy, is calculated by
is the total ablation factor and has a value of approximately 17 mg/kJ,
is the arc voltage, and
is the arc current.
2) The ablated materials, which are composed mainly of PTFE (C
), are mixed with cold SF
gas. The PTFE mass fraction can be calculated by an additional mass-fraction equation:
is the PTFE mass fraction, and
is the diffusion coefficient of PTFE, which is approximately 10
The thermodynamic and transport properties of a mixture of these gases are given in
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
to the electric-field intensity for the TRV
is the dielectric strength of the SF
gas in accordance with the temperature and pressure. Breakdown can occur when
is greater than
as a function of both temperature and pressure.
has a non-linear characteristic above approximately 800 K because the thermal dissociation and ionization of SF
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
Critical electric-field intensity for SF6 gas as function of both temperature and pressure.
It is noted that
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
is practically zero at 3,000 K.
used in this work, PTFE ablated from the nozzle is not considered. If the PTFE properties are considered,
might change. However,
of pure SF
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
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
. 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.
Entire computational domain for the 2-D axisymmetrical and 3-D asymmetrical analysis geometries.
shows the flowchart for the prediction algorithms of the dielectric strength under hot-gas conditions.
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
shows the 2-D output and 3-D inlet boundaries in the real CB geometry. The pressure
from the 2-D analysis is used as the inlet source of the 3-D analysis with the temperature
, and density
is the static pressure, and
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
. 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.
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
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
. The amplitude of the testing current is 40kA
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.
Installation of a 145kV self-blast gas CB in the testing room.
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
. 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)
show the contours of
and the location of the breakdown after the test respectively. The test results show that there is a dielectric breakdown at
= 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.
Contours of hot gas greater than 2000K before and after CZ for a 145kV self-blast gas CB.
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
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.
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.
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.
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.
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.
“Analysis of Dielecric Breakdown of Hot SF6 Gas in a Gas Circuit Breaker,”
“Real-Time Numerical Analysis on Insulation Capabiity Improvement of Compact Gas Circuit Breaker,”
IEEE Trans. Power Delivery
DOI : 10.1109/TPWRD.2007.899796
“Breakdown Electric Field Calculations of Hot SF6 for High Voltage Circuit Breaker Applications,”
IEEE Trans. on Dielect. Elect. Insulation
“Calculation of the Post-Arc Current in a Supersonic Nozzle by Using the κ-ε Turbulence Model,”
Journal of the Korean Physical Society
Zhang Jin Ling
Yan Jiu Dun
Murphy Anthony B.
“Computational Investigation of Arc Behavior in an Auto-Expansion Circuit Breaker Contaminated by Ablated Nozzle Vapor,”
IEEE Trans. on Plasma Science
DOI : 10.1109/TPS.2002.1024273
Thiel H G
“A computational fluid dynamics simulation of high- and low-current arcs in self-blast circuit breakers,”
J. Phys. D: Appl. Phys.
DOI : 10.1088/0022-3727/30/13/011
Boggs S. A.
“Optimization of Dead Tank Gas Circuit Breaker Design Based on Quantification of Hot Gas Flow and Dielecgtric Properties,”
IEEE Trans. Power Delivery
DOI : 10.1109/TPWRD.2003.820217
“Modelling of an ablation controlled arc,”
J. Phys. D:Appl. Phys.