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

This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/ 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

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

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.
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.
_{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:
where
ρ
is the density,
is the velocity,
S_{m}
is the mass source,
p
is the pressure,
is the stress tensor,
μ
is the molecular viscosity,
I
is the unit tensor,
h
is the enthalpy,
k
is the thermal conductivity,
S_{e}
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]
.
Thermal conductivity of SF_{6} gas according to temperature.
ṁ_{T}
, in accordance with the arc energy, is calculated by
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:
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]
.
E_{crit}
to the electric-field intensity for the TRV
E_{trv}
.
E_{crit}
is the dielectric strength of the SF
_{6}
gas in accordance with the temperature and pressure. Breakdown can occur when
E_{trv}
is greater than
E_{crit}
.
Fig. 3
shows
E_{crit}
as a function of both temperature and pressure.
E_{crit}
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
.
Critical electric-field intensity for SF_{6} gas as function of both temperature and pressure.
It is noted that
E_{crit}
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
E_{crit}
is practically zero at 3,000 K.
For
E_{crit}
used in this work, PTFE ablated from the nozzle is not considered. If the PTFE properties are considered,
E_{crit}
might change. However,
E_{crit}
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
E_{crit}
due to erosion of the arcing contacts is small and negligible.
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.
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.
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}
:
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.
Output variables and inlet boundary for the 3-D analysis.
_{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.
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
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
E_{trv}/E_{crit}
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.
Contours of hot gas greater than 2000K before and after CZ for a 145kV self-blast gas CB.
Contours of E_{trv}/E_{crit} 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
E_{crit}
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.
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.

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

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
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 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
PPT Slide

Lager Image

PPT Slide

Lager Image

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
PPT Slide

Lager Image

- 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

PPT Slide

Lager Image

- 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
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

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
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

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.

BIO

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

Citing 'Prediction of Insulation Capability for Ground Fault to Consider Asymmetry in SF6 Circuit Breaker
'

@article{ E1EEFQ_2015_v10n5_2046}
,title={Prediction of Insulation Capability for Ground Fault to Consider Asymmetry in SF6 Circuit Breaker}
,volume={5}
, url={http://dx.doi.org/10.5370/JEET.2015.10.5.2046}, DOI={10.5370/JEET.2015.10.5.2046}
, number= {5}
, journal={Journal of Electrical Engineering and Technology}
, publisher={The Korean Institute of Electrical Engineers}
, author={Oh, Yeon-Ho
and
Song, Ki-Dong
and
Kim, Hong-Kyu
and
Lee, Hae June
and
Hahn, Sung-Chin}
, year={2015}
, month={Sep}