Open access peer-reviewed chapter

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge

By Boxue Du, Hucheng Liang and Jin Li

Submitted: April 13th 2018Reviewed: July 31st 2018Published: November 5th 2018

DOI: 10.5772/intechopen.80635

Downloaded: 500


A majority of the high voltage (HV) electrical equipment which has solid-gas insulation has suffered greatly from the accumulation of the surface charges generated from the corona discharge. The local electric field may be distorted by the surface charge’s existence and in turn causes the surface flashover faults in excessive circumstances. Consequently, it’s significant to work out the mechanism of the procedure of the surface charge accumulation. A simulation model which combines both the charge trapping-detrapping procedure and the plasma hydrodynamics was created. The outcome of the simulation has agreed with the experimental results. The corona discharge intensity rises in the initial stage and then reduces as time goes by. There are various shapes of the surface potential distribution curves at various times. The central value increases quickly with time first and at last becomes saturated. Surface charges are observed in the epoxy insulator’s skin layer, some of them are mobile but some are captured by traps.


  • DC power transmission
  • epoxy insulator
  • corona discharge
  • FEM simulation
  • surface charge
  • trapping and detrapping

1. Introduction

A majority of the electrical equipment which has solid-gas insulation has suffered seriously from the existence of surface charges. Electrons and/or ions generated by the corona discharge migrate under electric force and accumulate on the insulator surface. In some cases, this may cause the local electric field distortion and even the surface flashover faults [1, 2, 3]. The surface charge distribution of GIS (gas insulated switchgear) spacer was presented in Ref. [4], which concluded that the distribution of surface charge always reaches its steady state after some time. Ref. [5] measured the surface potential decay (SPD) process of epoxy resin and found that this process takes several hours. Ref. [6] used the fluorination treatment to enhance the SPD rate of epoxy resin. Ref. [7] discussed the surface charge behaviors after various pulse application. In addition, several scholars have established the drift-diffusion equations for the purpose of describing the charge trapping-detrapping procedure within the bulk of the insulator [8, 9, 10, 11, 12]. Moreover, the plasma hydrodynamics models are adopted widely for the purpose of simulating the gas discharge procedure [13, 14, 15, 16]. Nevertheless, very few people have ever tried to combine the charge trapping-detrapping procedure and the plasma hydrodynamics to simulate the surface charge accumulation procedure.

This chapter presents a needle-plane model for the purpose of studying the surface charge accumulation procedure. It is the first time that the charge trapping-detrapping procedure combines with the plasma hydrodynamics. Compared with the existing simulation models for surface charge accumulation, our model has some advantages: (1) in the ionization region, the existing models used some simplified charge transport equations to simulated the generation and transport process of charged ions. Our model is closer to the reality with many physicochemical reactions (the collision ionization etc.) in consideration; (2) in the insulator bulk, the charge trapping-detrapping process is taken into consideration; (3) in our model, the charge transport parameters (carrier mobility, carrier diffusion coefficient etc.) are obtained from solving the Boltzmann equation, which is more reliable. This chapter aims at doing some fundamental researches on the procedure of the surface charge accumulation rather than guiding the engineering application. The needle-plane model has usually been adopted to do the SPD test in many published research papers, which is of great convenience for us to compare between the outcomes of experiment and simulation. This chapter may provide some help for readers to understand the surface charge accumulation process through some simulated details which is difficult to be gained from the experimental measurements.


2. Simulation model

2.1 Geometric model

According to Figure 1, the needle-plane electrode system’s schematic diagram is considered to be axisymmetric and therefore, it is simplified to a 2D issue. Within a lot of papers, the model often added with a mesh electrode was adopted for charging the insulators for the purpose of doing surface charge measurements. The needle electrode’s radius curvature is set to 50 μm in this chapter. The grounded (GND) electrode is placed which has a thickness of 0.5 mm and a radius of 20 mm. The axial distance between the insulator upper face and the needle tip is 3.5 mm.

Figure 1.

A schematic diagram of the needle-plane electrode system.

2.2 Corona discharge in air

On the basis of the needle-plane model shown before, the corona discharge process is going to be simulated in this section.

The governing equations for the drift-diffusion of the electron energy density nand the electron density neare presented below:


where the energy mobility με, the electron mobility μe, the energy diffusivity Dε, and the electron diffusivity Deare computed through working out the two-term Boltzmann equation. Rεand Reare the energy loss and electron source because of inelastic collisions. In this simulation, M-Nthree-body reactions and Ntwo-body reactions are taken into consideration,


where xj2 xj1and xj, are the species’ mole fractions which are involved in reaction j; Δεjand Δnejare the energy loss (V) and electron increment of reaction j, respectively; pis the atmospheric pressure of air (1 atm); Nnis the total neutral density (1/m3); and kjis the rate coefficient of reaction j(m3/s or m6/s), which is able to be gained through working out the two-term Boltzmann equation as well.


where εand Teis the mean electron energy (V) and the electron temperature (V), respectively.

In terms of heavy species, every species’ mass fraction is able to be gained through working out the equations below,


where Eis E-field strength (V/m); μkis averaged mobility of species k(m2/(V s)); zkis charge number of species k; Mnis s mean molar of air (kg/mol); Dkis averaged diffusion coefficient of species k(m2/s); Rk is generation rate of species k(kg/(m3·s)); jkis flux of species k; wkis the mass fraction of species k; ρis density of air (kg/m3). In this simulation, M-Nthree-body reactions and Ntwo-body reactions that change the species k’smass fraction are taken into consideration,


where xj3, xj2and xj1are the species’ mole fractions which are involved in reaction j; Mkis species k’s molar mass (kg/mol); Δnkjis species kincrement of reaction j; NAis Avogadro constant. Eq. (10) presents the relation between species k’s mass fraction and the mole fraction.

Eq. (12) is the definition of the mixture averaged diffusion coefficient Dkwhile Eq. (13) is the definition of the mixture averaged mobility μkin accordance with the Relation of Einstein,


where Tis gas temperature (K); kBis constant (J/K) of Boltzmann; eis unit charge (C); Dk,jis the binary diffusion coefficient between species kand j, which is able to be evaluated through the Fuller Formula,


where ∑vjand ∑vkis the species kand j‘s diffusion volume (cm3/mol); Mj and Mkare species kand j’smolar mass (kg/mol).

The boundary conditions on the surface of insulator and electrodes are defined as,


where γpis the secondary emission coefficient; εpis the secondary electrons’ mean energy (V); α·zk/|zk| = 1 if electric field can be directed towards the boundary; α·zk/|zk| = 0 if electric field can be directed away from the boundary; Tis regarded as the environment temperature (K) and vk,this the species k‘s thermal velocity (m/s).

The definitions of the boundary conditions at the open boundary are as follows:


Ions may become neutral species because of the surface reactions. Just several typical surface reactions will be considered for simplification in the paper.

The plasma hydrodynamics model is made up of 19 reactions and 10 species (e.g.: O4+, O, O3, O2, O2+, N2, N2, O2, N2+, O2+). Some particular physicochemical reactions are listed in Table 1 which are considered in the model after some reduction. These collision reactions’ energy losses and cross sections are extracted from papers [17, 18]. The initial electron density is set to 1 × 109 1/m3 and the O2 to N2 ratio is 1:4. ∆ne means the electron increment of each reaction.

No.FormulaTypeΔε(eV)ΔneRate coefficient
1N2 + e → e + N2Elastic00
2O2 + e → e + O2Elastic00
3N2 + e → 2e + N2+Ionization−15.61
4O2 + e → 2e + O2+Ionization−12.061
5N2 + e → e + N2Excitation−80
6O2 + e → e + O2Excitation−5.50
7O2++O2 + M → O4+ + MReaction2.04 × 10−34 T−3.2
8O4++e → 2O2Reaction−11.4 × 10−12 (300/Te)0.5
9O2++e → 2OReaction−12.42 × 10−13 (300/Te)
102O2 + e → O2 + O2Reaction−12 × 10−41 (300/Te)
11O4+ + O2 → 3O2Reaction1 × 10−13
12–13O4+ + O2 + M → 3O2 + MReaction2 × 10−37
14–15O2+ + O2 + M → 2O2 + MReaction2 × 10−37
16–17O + O2 + M → O3 + MReaction2.5 × 10−46
18e + N2+ + N2 → 2N2Reaction−16.07 × 10−34Te−2.5
192e + N2+ → N2 + eReaction−15.65 × 10−27Te−0.8

Table 1.

Some typical physicochemical reactions in the corona discharge model.

Units: m3 s−1 for two-body reactions, m6 s−1 for three-body reactions, K for Te and T; Notes: M = O2, N2.

2.3 Detrapping process and charge trapping in the epoxy insulator

Electrons are going to be injected into the discharge channel when it reaches the surface of the insulator. The positive ions are converted into the neutral particles through the reactions on the surface. It can be assumed that we can also inject some holes into the insulator. Both the holes and the electrons are possible to be released from traps through thermal excitation and captured by traps when transporting. It should be noticed that the conduction on the surface is not considered in this chapter because very little tangential component exists in the electric field distribution along the surface of the insulator.


where τis the non-equilibrium carriers’ lifetime (s); R3, R2and R1are the recombination coefficients (m3/s); Pdeand Ptrare the detrapping and trapping probability of hole and electron (1/s); nmbis the mobile electron density (1/m3); μhand μeare the hole and electron’s diffusion coefficient (m2/s); Vis potential (V); ntris the trapped electron density (1/m3); Jchand Jceare the mobile electron and hole’s flux (1/(m2 s)); htris trapped hole density (1/m3); hmbis the mobile hole density (1/m3).

Inspired by the theories in the semiconductor physics on the non-equilibrium, the products of the hole density h0 and mobile electron density are considered to be constants in the insulator within the condition of the thermal equilibrium.


where Egis the insulator’s energy gap (eV); Ncand Nvare states’ effective densities at conduction band bottom and the valence band top (1/m3). Supposing the hole density and electron are hmband nmbin an unbalanced condition, Δhand Δnis able to be obtained through working out Eq. (29). The procedure of the electron–hole recombination is represented by positive Δhand Δnwhile the electron-hole pairs’ generation is represented by negative Δhand Δn.


We can describe the boundary conditions at the insulator’s upper surface as follows:


where jiand jeare ion and electron flux through insulator’s upper surface in corona discharge model (1/m2 s).

We can describe the boundary conditions at the insulator’s lower surface as follows:


Some of the epoxy’s parameters adopted in the charge trapping-detrapping model are listed in Table 2. Among them, some referred to papers and some were gained from experimental measurements [12, 19, 20].

Charge carrier mobility (m2/(V s))
μe1.0 × 10−14
μh1.0 × 10−14
Charge carrier diffusion coefficient (m2/s)
De2.6 × 10−16
Dh2.6 × 10−16
Trapping and detrapping coefficients (1/s)
Ptr7.0 × 10−3
Pde7.7 × 10−5
Recombination coefficients (m3/s)
R18.0 × 10−19
R28.0 × 10−19
R38.0 × 10−19

Table 2.

Some parameters of epoxy in the trapping-detrapping model.

2.4 Poisson equation

The pspecies (negative ions, positive ions and electrons) are made up by the corona discharge model and zkis the species k’s charged number. Therefore, we can describe the Poisson equation in the air as follows:


where species k’s density is nk(1/m3).

We can describe Poisson equation within insulator as follows:


where εris the insulator’s relative dielectric constant.

We add a resistor Rbfor the purpose of limiting the discharge current Ip. The definition of the voltage on HV electrode is as follows:


where V0 is the voltage of power supply (V); Ipis the discharge current (A), which can be obtained by integrating the current density on the boundary of the HV electrode; Vis the potential on needle electrode (V).

The open boundary is as follows:


2.5. Surface potential measurement

Some surface potential measurement tests are implemented for the purpose of verifying this simulation’s validity. The measurement system of the needle-plane electrode surface potential is shown in Figure 2. The HV electrode was put 3.5 mm over the central insulator and supplied HVDC voltages on it through a power source with high voltage. Needle electrode’s radius curvature was ∼50 μm. The design of the insulator was like a disc of 0.5 mm thick. The insulator was slid to the position 3 mm on a rail below the probe after charging for a minute, which was adopted for measuring the distribution of the surface potential. It is not difficult to obtain the distribution of the surface potential along the radial distance through the measurement of 5 points’ values from the center to the insulator’s edge.

Figure 2.

A schematic diagram of the surface potential measurement system.

3. Results and discussion

With DC −5 kV applied on the needle, the electric potential and force line distribution are simulated and presented in Figure 3. The electric potential’s variation is represented by the gradient color which is presented in the color legend on the figure’s right side. The electric force line distribution can be obtained through the calculation of the potential gradient. It can be seen that all the force lines go from the needle to the ground and stronger field strength at the needle tip is indicated by denser force lines. Consequently, the corona discharge is produced like an avalanche. Electrons with negative charge transport from the needle tip to the insulator’s upper surface under electric force but ions with positive charge go oppositely, which causes the negative charges to accumulate on surface of insulator. In addition, very little tangential component of E-field exists on the surface of the insulator. Therefore, this chapter will not consider about the surface conduction.

Figure 3.

The simulated distribution of electric potential and force lines.

Electron density’s variation is shown in Figure 4. The electron density’s variation is represented by the gradient color which is presented in the color legend on figures’ right side. Here the E-field strength is high enough to produce electron-ion pairs and ionize the air due to the needle tip’s small curvature radius. According to Figure 4, it can be seen that the electron avalanche’s small crescent appears near the needle tip. Electrons gain higher speed from the needle tip under electric force, which produces more electrons through impact ionizations. We can see in Figure 5b that a clear increment exists in the electron avalanche’s density and size. The electron avalanche’s size becomes larger and larger at the time of moving to the insulator surface. The electron density’s distribution can be shown in Figure 4c. The electron avalanche’s head has reached the surface of the insulator at that time, which forms a clear discharge channel. Many electrons start accumulating on the surface of the insulator, which causes the E-field strength’s disadvantages and the surface potential’s rise. The discharge channel is barely able to be recognized from Figure 4d, which indicates the end of the corona discharge and surface charge accumulation.

Figure 4.

The variation of electron density with the development of corona discharge: (a) 1 × 10−9 s, (b) 1.8 × 10−8 s, (c) 4.6 × 10−7 s, (d) 60 s.

Figure 5.

The distribution of electron and ion density along the symmetry axis: (a) 1 × 10–9 s, (b) 1.8 × 10–8 s, (c) 4.6 × 10–7 s, (d) 60 s.

According to Figure 5, the ion and electron densities are shown along the symmetry axis, which echoes Figure 4 which is presented before. The axial coordinate z = −3.5 mm is the insulator surface and z = 0 mm is the needle tip. Figure 5a shows that the corona discharge is started from the needle tip. Positive charges concentrate at the tail and negative charges concentrate at the electron avalanche’s head. Electrons are accelerated away with time, which produces more ions and electrons. Figure 5b shows that the peaks of the ion and electron densities have reached approximately 1 × 1015 1/m3. It is clear that the electron avalanche grows sharply. The mobility of ions is much smaller than that of electrons. Therefore, they seldom move under electric force which results in the electron avalanche tail’s extension. Figure 5c shows that the electron avalanche has reached the surface of the insulator for the purpose of forming a discharge channel. Both the ion and electron densities are increased by more than magnitude’s four orders. Many negative charges have been accumulated as the discharge time goes by, which causes the field strength’s disadvantages and increases the surface potential. The procedure of the corona discharge becomes weak. Thus, both the ion and electron densities reduces with time, which is presented in Figure 5d.

According to Figure 6a, the surface potential distributions are shown with the corona discharge’s development. More electrons are able to be accumulated on the center of the surface because the needle electrode is above the center of the insulator directly. The peak of the surface potential appears in the central part and decreases to the edge from the central part. Additionally, the surface potential distributions’ curves have various kinds of shapes at various times. The surface potential is 2000 V at the center when time t = 1 s, which is much higher compared with that at the edge. One second is not very long that although many electrons have reached the surface of the insulator, hardly any electrons have enough time to diffuse. Consequently, the surface potential increases just by concentrating on the central surface’s small zone. When time t = 30 s, it is obvious to see not only a sharp reduction but also an increment in the center of the surface potential. It’s because an increasing number of surface charges have diffused with time from the central area. The surface potential gradient’s decrease is even clearer when time t = 60 s but the central surface potential’s growth becomes negligible.

Figure 6.

The surface potential and electric field distributions: (a) surface potential and (b) electric field.

According to Figure 6b, the E-field strength distributions are shown along the symmetry axis. It’s obvious that the field strength’s sharp decrease occurs at the interface of insulator and air because of the epoxy insulator’s higher relative dielectric constant compared with that of air. With the development of the corona discharge, an increasing number of electrons will exist on the insulator surface. The inside E-field strength’s increment in insulator is caused by higher surface potential. The E-field strength at needle tip decreases as time goes by due to the decrease of the potential difference between the insulator surface and the needle tip. The E-field strength at needle tip is presented in Figure 6b’s margin for the purpose of seeing the outcome obviously.

When the epoxy insulator is injected with electrons, they tend to transport to GND electrode across the bulk under the built-in electric area. Several mobile electrons are possible to be captures by traps at the time of transporting and later on gain de-trapped by thermal excitation. Owing to the exceeding existence and low mobility of trapped electrons, the injected electrons are able to be observed just in insulator’s skin layer. According to Figure 7, the distribution of trapped electrons and the mobile along symmetry axis is shown from GND electrode to insular upper surface. d = 0.5 mm is GND electrode while d = 0 mm is epoxy insulator’s upper surface. According to Figure 7a, it is quite obvious to see that mobile electron density decreases quickly from upper surface to ground. Mobile electron increases a bit deeper in the skin layer and decreases at the upper surface. Figure 6a concludes that the surface potential’s quick rise occurs at the beginning 1 s. In other words, most of the electrons will not be injected into epoxy insulator within the rest 59 s. Nevertheless, many electrons have gone deeper in skin layer under electric force. Consequently, electron density reduces at upper surface as time goes by. According to Figure 7b, the trapped electrons’ density increases obviously with time in the skin layer because an increasing number of mobile electrons get trapped but seldom escape at the time of transporting to the ground.

Figure 7.

The mobile and trapped electron density distribution along the symmetry axis of insulator bulk: (a) mobile electron and (b) trapped electron.

According to Figure 8, the central surface potential grows with time when various kinds of charging voltages are applied. It is clear central surface potential’s absolute value goes up quickly. The surface potential increases sharply at the beginning 1 s under −5 kV. Later on, the growth gradually slows down with time. At last, the surface potential may reach saturation with the dissipation procedure in dynamic equilibrium and the surface charge accumulation when there is enough discharge time. When the charging voltage reduces, the surface potential in the center will spend more time to reach saturation.

Figure 8.

The growth of central surface potential with time under different charging voltages.

According to Figure 9, the comparison between experimental and computational surface potential distributions is presented under −5 kV. The difference is very clear between the outcomes of the simulation and the experiment resources when t = 30 s which is presented in Figure 9a. The simulated value is higher than the experimental surface potential at the center. However, the experimental potential will become much higher within the area around ∼ 15 mm away from the center. Firstly, it is possible that the probe’s precision will not be high enough. Secondly, several seconds were taken to take off the power source and move the insulator before measuring its surface potential. Many accumulated electrons have spread in the surrounding areas from the center in the time delay, which increases the surface potential in the surrounding regions but decreases the surface potential at the center. According to Figure 9b, two similar curves of experimental and computational surface potential distributions are shown when t = 60 s. nevertheless, it is evident that the simulated values are higher than the experimental resources. The time delay is possible to be responsible before the surface potential measurement.

Figure 9.

The comparison between the computational and experimental surface potential distributions: (a) t = 30 s, (b) t = 60 s.

According to Figure 10, epoxy insulator’s central surface potentials are presented under various kinds of voltages after the charging for 60 s. It is obvious to see both the simulated and measured potential values grow up nearly in a linear way with the charging voltage at the insulator center. According to Figure 10b, the central potential grows with time under −5 kV. Experimental resources present a nice agreement with the outcomes of the simulation that central potential goes up quickly at the beginning seconds and at last reaches a fixed condition. Because of the diffusion of charges from center, the outcomes of the simulation are often higher than the measured surface potentials, particularly at the beginning 10 s. According to what has been discussed in Figure 5, the gradient of the surface potential distribution is much clearer along the radial distance at the beginning 10 s. Therefore, the central surface charges’ spread around is clearer to cause a greater effect on measured surface potential. Thus, the difference between the outcomes of the simulation and the experimental resources is wider at the beginning 10 s.

Figure 10.

The comparison between the computational and experimental central surface potentials: (a) the central surface potentials after 60s’ charging and (b) the growth of central surface potential with time.

4. Conclusions

The accumulation procedure of the surface charge under the needle-plane corona discharge has been explored in This chapter and some comparisons between experimental resources and the simulation outcomes have also been made. The major conclusion is shown below:

  1. In the accumulation procedure of the surface charge, the corona discharge intensity increases at the beginning and later reduces as time goes by. The epoxy insulator’s surface potential has increased sharply first and later slowly reaches saturation.

  2. The surface potential distributions’ curves have various kinds of shapes at various periods. As time goes by, the surface potential gradient along epoxy insulator’s radial distance reduces because of the surface charge’s diffusion.

  3. Meanwhile, a higher central surface potential exists in the epoxy insulator under higher charging voltage when it is being charged. Then, the central surface potential takes a shorter time to reach saturation.

  4. The epoxy insulator’s skin layer has a surface charge. These mobile electrons have a tendency to cross the insulator to the GND electrode. Among them, some get trapped and gradually get detrapped by thermal excitation.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Boxue Du, Hucheng Liang and Jin Li (November 5th 2018). Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge, Atmospheric Pressure Plasma - from Diagnostics to Applications, Anton Nikiforov and Zhiqiang Chen, IntechOpen, DOI: 10.5772/intechopen.80635. Available from:

chapter statistics

500total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Applications of Dielectric Barrier Discharge Microplasma

By Kazuo Shimizu, Jaroslav Kristof and Marius Gabriel Blajan

Related Book

First chapter

Introductory Chapter: Plasmonics

By Tatjana Gric

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us