Af25175180

Narapareddy Rama Rao, J.Amarnath / International Journal of Engineering Research and
Applications (IJERA) ISSN: 2248-9622 www.ijera.com
Vol. 2, Issue 5, September- October 2012, pp.175-180
Electric Field Effect On Particle Trajectory In Three Phase
Dielectric Coated Gas Insulated Busduct Using Numerical
Methods
Narapareddy Rama Rao*, J.Amarnath**
*
Department of Electrical and Electronics Engineering, Nigama Engineering College, Karim Nagar, AP, INDIA
**
Department of Electrical and Electronics Engineering, JNTUH College of Engineering, Hyderabad, AP,
INDIA

ABSTRACT
Free metallic conducting particles in Gas manufacturing of GIS equipment, care should be
Insulated Busduct(GIB) may cause loss of as much taken to ensure that all components are free from free
as 90% of the gas dielectric strength. Research metallic particles. However, metallic contaminants
studies revealed that 20% of failures in Gas are inevitable in installed systems. The most common
Insulated Substations (GIS) are due to presence of causes are mechanical vibrations during shipment and
free metallic contaminations in the form of loose service, thermal expansion/contraction at expansion
particles. Therefore, the presence of free metallic joints. Several methods of conducting particle
particles in Gas Insulated Busduct(GIB) operating control and deactivation have been proposed[4,5] and
at high electric fields can be a problem. If these one of them is dielectric coating on the inner surface
particles are eliminated from GIB, then the of the outer enclosure.
reliability of compressed Gas Insulated Substation The dielectric coatings can improve the
can be improved. The purpose of this paper is to insulation performance by reducing the charge on the
develop techniques for formulating the basic particle and thereby increase the lift off field[2,6,7].
equations that will control the movement of The work reported in this paper deals with the
metallic particles in a dielectric coated Gas techniques, which will formulate the basic equations
Insulated Busduct. Analytical and numerical that will govern the movement of metallic particles
techniques have been used to describe the effect of like aluminum and copper and also the maximum
the ambient electric field on the movement of movement is to be reduced for increasing the
metallic particle in GIB. The simulation has been dielectric strength of gaseous Insulation. The results
carried out to obtain the particle trajectories at are analyzed and presented.
various voltages for particles of Aluminum and
Copper. The results have been presented and II. MODELLING TECHNIQUE
analyzed in this paper. For this study a typical three phase common
enclosure horizontal busduct comprising of three
Keywords - Dielectric coat, Gas Insulated inner conductors A, B and C with dielectric coated
Busduct, Contaminated metallic particle, Finite outer enclosure filled with SF6 gas as shown in fig.1
Element Method, Metallic particle trajectory. is considered.
A wire like particle is assumed to be at rest
I. INTRODUCTION on the dielectric coated inner surface of enclosure.
Gas Insulated Substation (GIS) suffer from When a three phase voltage is applied to three phase
certain drawbacks and one of them is the outage due GIB, the particle acquires charge in the presence of
to seemingly innocuous conducting particles, which high electric fields and low gas pressures mainly due
accounts for nearly 20% of the GIS failures. These to two different particle charging mechanisms. They
particles may have any shape or size, may be are 1. Conduction through a dielectric coating 2.
spherical or filamentary (wire like) or in the form of Micro discharges between the particle and coating
fine dust. Flash over in a GIS is, in general, [7].
associated with longer outage times and greater costs An appropriate particle charge and electric
than in a conventional air insulated substation. The field causes the particle to lift and begins to move in
conducting particles can either be free to move in the the direction of the electric field after overcoming the
Gas Insulated Busduct(GIB) or they may stick to an forces due to its own weight and drag[5,6,8].
energized electrode or to an enclosure surface. Free
conducting particles in GIB could reduce the
insulation strength drastically as they can short-circuit
a part of the insulation distance, and thereby initiates a
breakdown and if these particles could be eliminated,
then this would improve the reliability of compressed
gas insulated substations[1,2,3]. At the time of

175 | P a g e

Cga ,Cgb and Cgc represents capacitances between the
three phase conductors and the particle, „Cc‟
represents capacitance between the particle and the
enclosure. The conductance „G‟ represents the part of
the dielectric coating where the Charging current is
flowing.
By using particle motion equation, the lift-
off field „Elo‟ of the metallic particle can be obtained
as,

(3)
Where „K‟ is a constant. „Cg‟ is effective
capacitance between three phase conductors and
metallic particle. „ω‟ is angular velocity, „T‟ is
thickness of dielectric coating, „S‟ is contact area
Figure 1. Typical three phase common enclosure Gas between particle and dielectric coating, „ρc‟ is
Insulated Busduct resistivity of dielectric material, „R‟ is resistance
between particle and GIB enclosure and „Cc‟ is
The simulation considers several parameters capacitance between particle and GIB enclosure. It
like the macroscopic field at the location of the
particle, its length, radius and weight, viscosity of the
gas, Reynold‟s number, drag coefficient and
coefficient of restitution[6] on its impact to the
enclosure. During the return flight, a new charge on
the particle is assigned, based on the instantaneous
electric field.
A conducting metallic particle moving under
the external electric field will be subjected to
Electrostatic force (Fe), Gravitational force (Fg) and
Drag force (Fd).
The equation of the motion for can be expressed as, can be noted that „Elo‟ is approximately proportional
to square root of the thickness and resistivity of the
dielectric.
(1)
Where m = mass of the particle, y = displacement in Fig.2 Circuit model of particle charging through the
vertical direction, g = gravitational constant. dielectric coating.
The motion equation using all forces can therefore The particle motion equation is a second
be expressed as[4-8]: order non-linear differential equation and it is solved
by using Runge- Kutta 4th Order Method.

III.SIMULATION OF METALLIC
PARTICLE
For the study of the motion of moving
metallic particle in Gas Insulated Substation requires
magnitude of the charge acquired by the particle and
0.25 electrostatic field present at the metallic particle
 C 
2

 c T 
0.5
1 location. The electric field at the metallic particle
Elo  K 1  c      location is calculated analytical method, Finite
 C g  R 2 2 C g 
2
 S 
   Difference Method[9 ] and also by using Finite
Element Method[9-11].

A. FINITE ELEMENT METHOD:
The following Figure 3 depicts the basic
concept for application of Fnite Element Method in
the required space for finding potentials at the
(2) element nodes.
Where E(t) is electric field intensity at time „t‟ at the
particle location.
The circuit model of particle charging
through the dielectric coating is as shown in Fig. 2.

176 | P a g e

region into grid points and finite difference five node
molecule.
1
𝑉𝑖,𝑗 = 𝑉𝑖+1,𝑗 + 𝑉𝑖−1,𝑗 + 𝑉𝑖,𝑗 +1 + 𝑉𝑖,𝑗 −1 (9)
4

Fig. 3: Finite element mesh for calculating potentials
at finite element nodes.
The Total Energy(W) associated with the
assemblage of all elements in Gas Insulated Busduct
is,

The finite difference approximation to poission‟s
equation applied to Gas Insulated Busduct solution
region for five-node is:
1
(4) 𝑉𝑜 = 𝑉 + 𝑉2 + 𝑉3 + 𝑉4 (10)
4 1
The above equation shows the average value property
Where „N‟ is number of elements, „V‟ is of Laplace equation. Using Band Matrix
node voltage matrix of „n‟ nodes and „C‟ is overall or Method[105] to all free nodes in solution region
global coefficient matrix. results in a set of simultaneous equations of the form
In Finite Element Method, the solution 𝐴 𝑉 = 𝐵 (11)
region has minimum total energy satisfying the Where 𝐴 is sparce matrix, 𝑉 is free node
laplace‟s or poission‟s equation. So, partial
potentials and 𝐵 is column matrix of potentials at
derivatives of „W‟ with respect to each nodal value of
fixed nodes of solution region. Therefore unknown
potential must be zero.
potentials at free nodes is,
𝑉 = 𝐴 −1 𝐵 (12)
Electric Field intensity at any point in Gas Insulated
Busduct is calculated by using following equation,
(5) 𝐸 = −∇ 𝑉 (13)
In general, simplifying the finite element mesh,
C. ANALYTICAL METHOD:
Analytically ambient electric field „E‟ at any time
in three phase GIB can be calculated by using
(6) following equation,
Where i is number of nodes and k=1,2,3,………n. 1   1   Cos 2  
E  Va  (Vb  Vc )
So, a set of „n‟ simultaneous equations are obtained ln( h / Rc )   h  x   Rbx
    
 

and solving the above simultaneous equations using
band matrix method for unknown node voltages(Vf), (14)

Where Va, Vb, and Vc are phase voltages of
A, B, and C conductors respectively, Rc is the high
(7) voltage conductor radius, Rbx is distance between B
Where Vf is free node voltage matrix, V p phase conductor and particle location, „Ө2” is the
prescribed or fixed node voltage matrix, Cff free node angle between Rbx and vertical axis at B or C phase
global coefficient matrix and Cfp is free to prescribed conductor and „x‟ is the distance from enclosure inner
nodes global coefficient matrix. surface to the position of the particle which is moving
Electric Field intensity at any point in Gas Insulated upwards.
Busduct is calculated by using following equation,
IV. RESULTS AND DISCUSSION
The motion equation of metallic particle is
solved by using RK 4th Order method and it gives
(8) movement in the radial direction only. The Axial
B. FINITE DIFFERENCE METHOD: movement of the metal particle is calculated by using
Figure 3.6 Finite difference solution pattern Monte-Carlo Technique based on the works of
The figure 3.6 shows the finite difference J.Amarnath et al[5]. The Electric fields are
solution pattern representing division of solution determined by using equations(12) and (13) for Finite
Difference Method, Finite Element Method as given

177 | P a g e

by equations (7) and (8) and with Analytical Method FEM calculated fields respectively for 600kV. Table
using equation (14). I show the maximum radial movements for
Computer simulations of motion for the aluminium and copper particles for different voltages
metallic wire particles were carried out using with analytical, Finite Difference and Finite Element
Advanced C Language Program for GIB with high calculated electric fields.
voltage conductors diameter 64mm and enclosure
diameter of 500mm for 220kV, 300kV, 400kV, Table II Maximum Axial Movements of Aluminium
500kV and 600kV applied voltages. Aluminum and and Copper particles.
copper wire particles were considered for Max.Axia
Type Max.Axi
simulations. l Max.Axia
Table I and Table II are showing the Vol of al
Movemen l
movement patterns of aluminium and copper particles tag Particl Moveme
t (mm) Movemen
e in e nt (mm)
for various power frequency voltages. 289 nodes are with t (mm)
kV Materi with
considered for finite difference and finite element analytical with FEM
al FDM
methods in dielectric coated three phase Gas method
Insulated Busduct space for calculating node Al 52.98 192.26 195.40
potentials and thereby electric field at respective 220
Cu 8.48 49.69 59.04
nodes. The radius of aluminium and copper particles
Al 73.55 275.43 273.32
in all cases are considered as 0.01mm, length of the 300
Cu 15.32 84.04 105.18
particle as 12mm, restitution coefficient is 0.9 and
Al 95.90 418.75 476.38
SF6 gas pressure is 0.4MPa. 400
During application of power frequency Cu 22.59 146.32 197.61
voltage, the moving metallic particle makes several Al 127.09 503.80 555.92
500
impacts with the enclosure and the maximum radial Cu 36.24 218.78 237.59
movement increases with increase of applied voltage. Al 149.10 423.38 541.78
600
For Aluminium metallic particles the maximum radial Cu 50.11 248.72 385.48
movement is 2.39mm, 8.69mm and 10.75mm with
analytical, FDM and FEM calculated field for 220kV Similarly for Aluminium particles the
and the radial movement is increasing with increase maximum axial movement is 52.98mm, 192.26mm
of applied voltage and reaching maximum movement and 195.40mm with analytical, FDM and FEM
of 8.38mm, 34.96 and 35.66mm with analytical, calculated fields respectively for 220kV. The
FDM and FEM calculated field at 600kV maximum axial movements of Aluminium particles
respectively. are increasing with increase of voltages upto 500kV
Table I Maximum Radial Movements of aluminum and after that it is observed that slight decrease in
and copper particles. movements. For 800kV these maximum movements
Type Max. are reaching 149.10mm, 423.38mm and 541.78mm
Max. Max. for electric fields calculated with analytical, FDM and
of Radial
Volta Radial Radial FEM methods respectively. For Copper particles the
Partic Movement
ge in Movement Movemen maximum axial movements are 8.48mm, 49.69mm
le (mm) with
kV (mm) with t (mm) and 59.04mm with analytical, FDM and FEM
Mater Analytical
FDM with FEM calculated fields respectively. The maximum axial
ial Method
Al 2.39 8.69 10.75 movements of Copper particles are increasing with
220 increase of voltage and for 600kV these maximum
Cu 0.57 3.11 4.14
Al 3.19 13.91 16.29 movements are reaching 50.11mm, 248.72mm and
300 385.48mm for analytical, FDM and FEM calculated
Cu 1.01 5.76 7.63
fields respectively. Table II represents Aluminium
Al 4.29 21.06 25.63
400 and Copper particle maximum axial movement for
Cu 1.78 9.58 12.35
different voltages with analytical, FDM and FEM
Al 6.41 27.57 32.93
500 calculated electric fields.
Cu 2.42 14.31 17.66 Fig. 3 to Fig. 8 show the radial movement
Al 8.38 34.96 35.66 patterns of aluminum and copper particles using
600
Cu 3.19 18.61 24.36 Analytical, FDM and FEM calculated fields at power
frequency voltage of 400kV.
For copper metallic particles the radial
movement for electric fields at 220kV is 0.57mm,
3.11mm and 4.14mm with analytical, FDM and FEM
calculated fields. For copper particles also, the radial
movement is increasing with increase of applied
voltage and reaching maximum value of 3.19mm,
18.61mm and 24.36mm with analytical, FDM and

178 | P a g e

Radial Movement in mm

5 15
4
3 10
2
5
1
0 0
0 0.2 0.4 0.6 0.8 1 1.2
0 0.2 0.4 0.6 0.8 1 1.2
Time in Seconds
T ime in Seconds

Fig.3 Al particle radial movement for 400kV with
Fig.8 Cu particle radial movement for 400kV with
analytically calculated field
FEM calculated field

25
IV. CONCLUSION
20
A mathematical model has been formulated
15 to simulate the movement of wire like particle under
10 the influence of electric field in dielectric coated three
5 phase Gas Insulated Busduct. When an electrostatic
0 force on metallic particle exceeds the gravitational
0 0.2 0.4 0.6 0.8 1 1.2 and drag forces, the particle lifts from its position. A
T ime in Seconds
further increase in the applied voltage makes the
particle to move into the inter electrode gap in the
direction of applied field. The influence of increased
Fig.4 Al particle radial movement for 400kV with
voltage level on the motion of the metal particles is
FDM calculated field
also investigated. For instance, it can be noted that

30
aluminum particles are more influenced by the
voltage than copper particles due to their lighter
20 mass. This results in the aluminum particle acquiring
10 greater charge-to-mass ratio. Monte-Carlo simulation
is also adopted to determine axial as well as radial
0
0 0.2 0.4 0.6 0.8 1 1.2
movements of particle in the busduct. Also it is
observed that particle maximum movement in electric
Time in Seconds
field calculated by using Finite Element Method is
more than the maximum movement in electric field
Fig.5 Al particle radial movement for 400kV with calculated by using analytically calculated electric
FEM calculated field field. The dielectric performance is improved by
coating the inner surface of the enclosure with a light

shade with micron thickness epoxy resin. This
2
dielectric coating improves the insulation
1.5
performance by reducing the charge on the particle
1
and the maximum movement is reduced. All the
0.5
0
above investigations have been carried out for various
0 0.2 0.4 0.6 0.8 1 1.2
voltages under power frequency. The results
obtained are analyzed and presented.
T ime in Seconds

ACKNOWLEDGEMENT
Fig. 6 Cu particle radial movement for 400kV with
analytically calculated field The authors are thankful to the managements
of Nigama Engineering College, Karim Nagar, and

12
JNTUH University, Hyderabad, for providing
10 facilities and to publish this work.
8
6 REFERENCES
4
2
[1] Anis H.and Srivastava K.D.,1989, “
0 Breakdown characteristic of dielectric
0 0.2 0.4 0.6 0.8 1 1.2 coated electrodes in SF6 as with particle
contamination,”, sixth Intl symposium on
T ime in Seconds
HVE, New Orleans, LA,USA, paper No.32-
06
Fig.7 Cu particle radial movement for 400kV with [2] K.D. Srivastava and R.G. Van Heeswi jk
FDM calculated field “Dielectric coatings Effect of break down

179 | P a g e

and particle movement in GITL systems”,
IEEE transactions on Power Apparatus and
Systems, vol. PAS- 104, No.1, pp. 22-31,
January 1985.
[3] J. R. Laghari and A.H.Qureshi “A Review
of Particle Contaminated Gas Breakdown”
IEEE Transactions on Electrical Insulation,
Vol. EI-16 No.5, pp. 388-398, October
1981.
[4] H.Anis and K.D.Srivastava “Free
conducting particles in compressed gas
insulation”, IEEE Tranctions on electrical
insulation, Vol. EI-16, pp.327-338, August
1995.
[5] M.M. Morcos, S.Zhang, K.D. Srivastava,
and S.M. Gubanski,“Dynamics of Metallic
particle contamination in GIS with
dielectric coating electrodes”, IEEE Trans.
Power Delivery Vol.15, pp. 455-460, 2000.
[6] J. Amarnath, B.P. Singh, C. Radhakrishna
and S. Kamakshiah, “Determination of
Particle trajectory in a Gas Insulated
Busduct predicted byMonte-Carlo
technique”, IEEE Conference on Electrical
Insulation and Dielectric Phenomena
(CEIDP), Texas, Austin, USA, 1991 Vol.1,
pp. 399-402,1991.
[7] G.V.Nagesh Kumar, J.Amarnath and
B.P.Singh, “Behavior of Metallic Particles in
a Single Phase Gas Insulated System with
Dielectric Coated Electrodes”, International
Conference on Condition Monitoring and
Diagnosis, Beijing, China, April 21-24,
2008.
[8] K.B.Madhu Sahu and J.Amarnath, “Effect of
Various Parameters on the Movement of
Metallic Particles in a Single Phase Gas
Insulated Busduct with Image Charges and
Dielectric Coated Electrodes”, ARPN
Journal of Engineering and Applied
Sciences, Vol.5, No.6, June 2010, pp.52-60.
[9] Sadiku Matthew, “Elements of
Electromagnetics”, 4th Edition, (New
Delhi,Oxford Publishers,2007), pp 694-711.
[10] J.Weiss and Z.J.Csendes “A One step
Finite Element Method for Multiconductor
Skin Effect Problems”, IEEE Transactions
on Power Apparatus and Systems, Vol.
PAS-101, No.10 October 1982.
[11] S.Cristina and M.Feliziani, “A Finite
Element Technique for Multiconductor
Cable Parameters Calculation”, IEEE
Transactions on Magnetics, Vol.25, No.4,
July 1989.

180 | P a g e

Af25175180

More Related Content

What's hot

Viewers also liked

Similar to Af25175180

Af25175180