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

Boxue Du, Hucheng Liang and Jin Li

#### Abstract

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.

Keywords: 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–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 trappingdetrapping procedure within the bulk of the insulator [8–12]. Moreover, the plasma hydrodynamics models are adopted widely for the purpose of simulating the gas discharge procedure [13–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.

2.2 Corona discharge in air

process is going to be simulated in this section.

DOI: http://dx.doi.org/10.5772/intechopen.80635

and the electron density ne are presented below:

reactions are taken into consideration,

(1/m3

(m<sup>2</sup>

53

(V), respectively.

through working out the equations below,

is averaged diffusion coefficient of species k (m<sup>2</sup>

On the basis of the needle-plane model shown before, the corona discharge

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

The governing equations for the drift-diffusion of the electron energy density n

where the energy mobility με, the electron mobility μe, the energy diffusivity Dε, and the electron diffusivity De are computed through working out the two-term Boltzmann equation. R<sup>ε</sup> and Re are the energy loss and electron source because of inelastic collisions. In this simulation, M-N three-body reactions and N two-body

where xj2 xj1 and xj, are the species' mole fractions which are involved in reaction j; Δε<sup>j</sup> and Δnej are the energy loss (V) and electron increment of reaction j, respectively; p is the atmospheric pressure of air (1 atm); Nn is the total neutral density

where ε and Te is the mean electron energy (V) and the electron temperature

In terms of heavy species, every species' mass fraction is able to be gained

where E is E-field strength (V/m); μ<sup>k</sup> is averaged mobility of species k

/(V s)); zk is charge number of species k; Mn is s mean molar of air (kg/mol); Dk

/s or m<sup>6</sup>

); and kj is the rate coefficient of reaction j (m<sup>3</sup>

gained through working out the two-term Boltzmann equation as well.

ð1Þ

ð2Þ

ð3Þ

ð4Þ

ð5Þ

ð6Þ

ð7Þ

ð8Þ

ð9Þ

/s), which is able to be

/s); Rk is generation rate of species

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 trappingdetrapping 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.

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

#### 2.2 Corona discharge in air

combine the charge trapping-detrapping procedure and the plasma hydrodynamics

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 trappingdetrapping 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

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.

to simulate the surface charge accumulation procedure.

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

measurements.

Figure 1.

52

A schematic diagram of the needle-plane electrode system.

2. Simulation model

2.1 Geometric model

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 n and the electron density ne are presented below:

$$\frac{\partial}{\partial t}(\boldsymbol{n}\_{\boldsymbol{e}}) + \nabla \cdot [-\boldsymbol{n}\_{\boldsymbol{e}}(\boldsymbol{\mu}\_{\boldsymbol{e}} \cdot \vec{E}) - D\_{\boldsymbol{e}} \cdot \nabla \boldsymbol{n}\_{\boldsymbol{e}}] = \boldsymbol{R}\_{\boldsymbol{e}} \tag{1}$$

$$\frac{\partial}{\partial t}(\boldsymbol{n}\_{\boldsymbol{\varepsilon}}) + \nabla \cdot [-\boldsymbol{n}\_{\boldsymbol{\varepsilon}}(\boldsymbol{\mu}\_{\boldsymbol{\varepsilon}} \cdot \vec{E}) - D\_{\boldsymbol{\varepsilon}} \cdot \nabla \boldsymbol{n}\_{\boldsymbol{\varepsilon}}] + E \cdot \boldsymbol{\Gamma}\_{\boldsymbol{\varepsilon}} = \boldsymbol{R}\_{\boldsymbol{\varepsilon}} \tag{2}$$

where the energy mobility με, the electron mobility μe, the energy diffusivity Dε, and the electron diffusivity De are computed through working out the two-term Boltzmann equation. R<sup>ε</sup> and Re are the energy loss and electron source because of inelastic collisions. In this simulation, M-N three-body reactions and N two-body reactions are taken into consideration,

$$R\_c = \sum\_{j=1}^{N} \mathbf{x}\_j k\_j N\_n n\_a \Delta n\_{aj} + \sum\_{j=N+1}^{M} \mathbf{x}\_{j1} \mathbf{x}\_{j2} k\_j N\_n^{\ \ 2} n\_a \Delta n\_{aj} \tag{3}$$

$$\mathcal{R}\_e = \sum\_{j=1}^{P} \mathbf{x}\_j k\_j N\_n n\_e \Delta \varepsilon\_j \tag{4}$$

$$N\_n = pT \,/k\_B \tag{5}$$

where xj2 xj1 and xj, are the species' mole fractions which are involved in reaction j; Δε<sup>j</sup> and Δnej are the energy loss (V) and electron increment of reaction j, respectively; p is the atmospheric pressure of air (1 atm); Nn is the total neutral density (1/m3 ); and kj is the rate coefficient of reaction j (m<sup>3</sup> /s or m<sup>6</sup> /s), which is able to be gained through working out the two-term Boltzmann equation as well.

$$
\mathcal{E} = \frac{n\_c}{n\_s} = \frac{3}{2} T\_s \tag{6}
$$

where ε and Te is 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,

$$
\rho \frac{\partial}{\partial t}(\mathbf{w}\_k) - \nabla \cdot \vec{\mathbf{j}}\_k = \mathbf{R}\_k \tag{7}
$$

$$
\vec{j}\_k = \rho w\_k \vec{\mathcal{V}}\_k^\top \tag{8}
$$

$$\vec{V}\_k = D\_k(\nabla w\_k \, / \, w\_k + \nabla M\_n \, / M\_n) - z\_k \mu\_k \vec{E} \tag{9}$$

where E is E-field strength (V/m); μ<sup>k</sup> is averaged mobility of species k (m<sup>2</sup> /(V s)); zk is charge number of species k; Mn is s mean molar of air (kg/mol); Dk is averaged diffusion coefficient of species k (m<sup>2</sup> /s); Rk is generation rate of species k (kg/(m3 �s)); jk is flux of species k; wk is the mass fraction of species k; ρ is density of air (kg/m<sup>3</sup> ). In this simulation, M-N three-body reactions and N two-body reactions that change the species k's mass fraction are taken into consideration,

$$\begin{aligned} R\_k &= \frac{M\_k}{N\_A} \{ \sum\_{j=1}^N k\_j \mathbf{x}\_{jt} \mathbf{x}\_{j2} \mathcal{N}\_n^{-2} \Delta n\_{ky} \\ &+ \sum\_{j=N+1}^M k\_j \mathbf{x}\_{j1} \mathbf{x}\_{j2} \mathbf{x}\_{j3} \mathcal{N}\_n^{-3} \Delta n\_{ky} \} \end{aligned} \tag{10}$$
 
$$\mathbf{x}\_k = \frac{M\_n \mathbf{w}\_k}{M\_k} \tag{11}$$

ð19Þ

ð20Þ

ð21Þ

ð22Þ

0.5

�2.5

�0.8

). Some particular physicochemical

where γ<sup>p</sup> is the secondary emission coefficient; ε<sup>p</sup> is 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; T is regarded as the environment temperature (K) and vk,th is the species k 's thermal velocity (m/s). The definitions of the boundary conditions at the open boundary are as follows:

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

Ions may become neutral species because of the surface reactions. Just several

The plasma hydrodynamics model is made up of 19 reactions and 10 species

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

No. Formula Type Δε (eV) Δne Rate coefficient

<sup>+</sup> Ionization �15.6 1 —

<sup>+</sup> Ionization �12.06 1 —

+e ! 2O2 Reaction — �1 1.4 � <sup>10</sup>�<sup>12</sup> (300/Te)

+e ! 2O Reaction — �1 2.42 � <sup>10</sup>�<sup>13</sup> (300/Te)

� ! 3O2 Reaction — — <sup>1</sup> � <sup>10</sup>�<sup>13</sup>

� + M ! 3O2 + M Reaction — — <sup>2</sup> � <sup>10</sup>�<sup>37</sup>

� + M ! 2O2 + M Reaction — — <sup>2</sup> � <sup>10</sup>�<sup>37</sup>

<sup>+</sup> + N2 ! 2N2 Reaction — �1 6.07 � <sup>10</sup>�<sup>34</sup>Te

<sup>+</sup> ! N2 + e Reaction — �1 5.65 � <sup>10</sup>�<sup>27</sup>Te

<sup>16</sup>–17 O + O2 + M ! O3 + M Reaction — — 2.5 � <sup>10</sup>�<sup>46</sup>

<sup>+</sup> + M Reaction — — 2.04 � <sup>10</sup>�<sup>34</sup> <sup>T</sup>�3.2

� Reaction — �1 2 � <sup>10</sup>�<sup>41</sup> (300/Te)

�<sup>1</sup> for three-body reactions, K for Te and T; Notes: M=O2, N2.

+ , O2 +

typical surface reactions will be considered for simplification in the paper.

, N2, N2, O2, N2

1 N2 + e ! e+N2 Elastic 0 0 — 2 O2 + e ! e+O2 Elastic 0 0 —

5 N2 + e ! e+N2 Excitation �8 0 — 6 O2 + e ! e+O2 Excitation �5.5 0 —

(e.g.: O4

+

, O, O3, O2

3 N2 + e ! 2e + N2

4 O2 + e ! 2e + O2

+

+

+

10 2O2 + e ! O2 + O2

<sup>+</sup> + O2

<sup>+</sup> + O2

<sup>+</sup> + O2

�<sup>1</sup> for two-body reactions, m<sup>6</sup> s

Some typical physicochemical reactions in the corona discharge model.

+O2 + M ! O4

7 O2

8 O4

9 O2

11 O4

12–13 O4

14–15 O2

18 e + N2

19 2e + N2

Units: m<sup>3</sup> s

Table 1.

55

�, O2 +

DOI: http://dx.doi.org/10.5772/intechopen.80635

where xj3, xj2 and xj1 are the species' mole fractions which are involved in reaction j; Mk is species k's molar mass (kg/mol); Δnkj is species k increment of reaction j; NA is 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 Dk while Eq. (13) is the definition of the mixture averaged mobility μ<sup>k</sup> in accordance with the Relation of Einstein,

$$D\_k = \frac{\text{l} \cdot \text{w}\_k}{\sum\_{j \neq k}^{\mathcal{Q}} \text{x}\_j \mid D\_{k,j}} \tag{12}$$

$$
\mu\_k = \frac{eD\_k}{k\_B T} \tag{13}
$$

where T is gas temperature (K); kB is constant (J/K) of Boltzmann; e is unit charge (C); Dk,j is the binary diffusion coefficient between species k and j, which is able to be evaluated through the Fuller Formula,

$$D\_{k,j} = \frac{0.0101T^{1.75} \sqrt{\frac{1}{M\_k} + \frac{1}{M\_j}}}{P\|(\sum \nu\_k)^{1/3} + \sum \nu\_j)^{1/3}\|^2} \tag{14}$$

where ∑vj and ∑vk is the species k and j's diffusion volume (cm3 /mol); M<sup>j</sup> and Mk are species k and j's molar mass (kg/mol).

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

$$
\vec{m} \cdot \vec{\Gamma}\_e = \frac{1}{4} \nu\_{e,0} n\_e - \sum\_p \gamma\_p \left(\vec{\Gamma}\_p \cdot \vec{n}\right) \tag{15}
$$

$$
\vec{m} \cdot \vec{\Gamma}\_a = \frac{\mathbf{S}}{12} \nu\_{a,b} n\_a - \sum\_p \mathfrak{s}\_p \mathcal{V}\_p \left( \vec{\Gamma}\_p \cdot \vec{m} \right) \tag{16}
$$

$$\nu\_{e,th} = (\frac{8k\_B T\_e}{\pi m\_e})^{1/2} \tag{17}$$

$$
\vec{n} \cdot \vec{j}\_k = \frac{1}{4} \nu\_{k, \ell h} \rho \nu v\_k + \alpha \rho \nu v\_k \,\pi\_k \,\mu\_k \,(\vec{E} \cdot \vec{n}) \tag{18}
$$

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

$$\mathbf{v}\_{k,tb} = (\frac{8k\_{\theta}T}{\pi m\_{k}})^{1/2} \tag{19}$$

where γ<sup>p</sup> is the secondary emission coefficient; ε<sup>p</sup> is 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; T is regarded as the environment temperature (K) and vk,th is the species k 's thermal velocity (m/s).

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

$$
\vec{n} \cdot \nabla n\_a = 0 \tag{20}
$$

$$
\vec{n} \cdot \nabla n\_a = 0 \tag{21}
$$

$$
\vec{n} \cdot \nabla \mathcal{W}\_k = 0 \tag{22}
$$

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

$$\begin{array}{ccc} \mathrm{N}^{\perp}\rightharpoonup \mathrm{N}\_{2} & & \mathrm{O}z^{\perp}\rightharpoonup \mathrm{O}z \\\\ \mathrm{O}z^{\perp}\rightharpoonup \mathrm{O}z & & \mathrm{O}z^{\perp}\rightharpoonup \mathrm{O}z \end{array}$$

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


#### Table 1.

Some typical physicochemical reactions in the corona discharge model.

k (kg/(m3

of air (kg/m<sup>3</sup>

k's mass fraction and the mole fraction.

able to be evaluated through the Fuller Formula,

Mk are species k and j's molar mass (kg/mol).

Relation of Einstein,

54

�s)); jk is flux of species k; wk is the mass fraction of species k; ρ is density

ð10Þ

ð11Þ

ð12Þ

ð13Þ

ð14Þ

ð15Þ

ð16Þ

ð17Þ

ð18Þ

/mol); M<sup>j</sup> and

). In this simulation, M-N three-body reactions and N two-body

reactions that change the species k's mass fraction are taken into consideration,

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

where xj3, xj2 and xj1 are the species' mole fractions which are involved in reaction j; Mk is species k's molar mass (kg/mol); Δnkj is species k increment of reaction j; NA is Avogadro constant. Eq. (10) presents the relation between species

Eq. (12) is the definition of the mixture averaged diffusion coefficient Dk while Eq. (13) is the definition of the mixture averaged mobility μ<sup>k</sup> in accordance with the

where T is gas temperature (K); kB is constant (J/K) of Boltzmann; e is unit charge (C); Dk,j is the binary diffusion coefficient between species k and j, which is

where ∑vj and ∑vk is the species k and j's diffusion volume (cm3

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

papers [17, 18]. The initial electron density is set to 1 � <sup>10</sup><sup>9</sup> 1/m<sup>3</sup> and the O2 to N2 ratio is 1:4. Δn<sup>e</sup> means the electron increment of each reaction.

#### 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.

$$\frac{\partial n\_{mb}}{\partial t} + \nabla \cdot \vec{J}\_c^\* = n\_{tr}P\_{dc} - n\_{mb}P\_{tr} - R\_{l}n\_{mb}h\_{tr} - \frac{1}{\pi}\Delta n \tag{23}$$

$$\frac{\partial n\_{tr}}{\partial t} = n\_{mb}P\_{tr} - R\_2 n\_{tr}h\_{tr} - R\_3 n\_{tr}h\_{mb} \tag{24}$$

ð31Þ

ð32Þ

ð33Þ

ð34Þ

ð35Þ

ð36Þ

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

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

where ji and je are ion and electron flux through insulator's upper surface in

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

The p species (negative ions, positive ions and electrons) are made up by the corona discharge model and zk is the species k's charged number. Therefore, we can

gained from experimental measurements [12, 19, 20].

describe the Poisson equation in the air as follows:

/(V s))

Parameters Value

<sup>μ</sup><sup>e</sup> 1.0 � <sup>10</sup>�<sup>14</sup> <sup>μ</sup><sup>h</sup> 1.0 � <sup>10</sup>�<sup>14</sup>

/s) De 2.6 � <sup>10</sup>�<sup>16</sup> Dh 2.6 � <sup>10</sup>�<sup>16</sup>

Ptr 7.0 � <sup>10</sup>�<sup>3</sup> Pde 7.7 � <sup>10</sup>�<sup>5</sup>

<sup>R</sup><sup>1</sup> 8.0 � <sup>10</sup>�<sup>19</sup> <sup>R</sup><sup>2</sup> 8.0 � <sup>10</sup>�<sup>19</sup> <sup>R</sup><sup>3</sup> 8.0 � <sup>10</sup>�<sup>19</sup>

/s)

Some parameters of epoxy in the trapping-detrapping model.

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

corona discharge model (1/m2 s).

DOI: http://dx.doi.org/10.5772/intechopen.80635

2.4 Poisson equation

Charge carrier mobility (m<sup>2</sup>

Charge carrier diffusion coefficient (m2

Trapping and detrapping coefficients (1/s)

Recombination coefficients (m<sup>3</sup>

Table 2.

57

$$
\vec{J}\_c^\* = \mu\_c n\_{mb} \cdot \nabla V - D\_c \cdot \nabla n\_{mb} \tag{25}
$$

$$\frac{\partial h\_{mb}}{\partial t} + \nabla \cdot \vec{J}\_c^h = h\_\nu P\_{\text{dc}} - h\_{mb} P\_\nu \cdot R\_3 n\_\nu h\_{mb} - \frac{1}{\tau} \Delta h \tag{26}$$

$$\frac{\partial h\_{tr}}{\partial t} = h\_{m\flat} P\_{\iota r} \text{ - } R\_{\natural} n\_{m\flat} h\_{\iota r} \text{ - } R\_{\natural} n\_{\iota r} h\_{\iota r} \tag{27}$$

$$
\overline{J}^{h}\_{\cdot} = -\mu\_{h}\hbar\_{ab} \cdot \nabla V - D\_{h} \cdot \nabla h\_{mb} \tag{28}
$$

where τ is the non-equilibrium carriers' lifetime (s); R3, R2 and R1 are the recombination coefficients (m3 /s); Pde and Ptr are the detrapping and trapping probability of hole and electron (1/s); nmb is the mobile electron density (1/m3 ); μ<sup>h</sup> and μ<sup>e</sup> are the hole and electron's diffusion coefficient (m2 /s); V is potential (V); ntr is the trapped electron density (1/m3 ); Jc <sup>h</sup> and Jc <sup>e</sup> are the mobile electron and hole's flux (1/(m2 s)); htr is trapped hole density (1/m3 ); hmb is the mobile hole density (1/m3 ).

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

$$m\_0 h\_0 = N\_\upsilon N\_c e^{-E\_\kappa/k\_b T} \tag{29}$$

where Eg is the insulator's energy gap (eV); Nc and Nv are states' effective densities at conduction band bottom and the valence band top (1/m<sup>3</sup> ). Supposing the hole density and electron are hmb and nmb in an unbalanced condition, Δh and Δn is able to be obtained through working out Eq. (29). The procedure of the electron–hole recombination is represented by positive Δh and Δn while the electron-hole pairs' generation is represented by negative Δh and Δn.

$$n\_0 h\_0 = (n\_{mb} - \Delta n)(h\_{mb} - \Delta h) \tag{30}$$

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

$$
\Delta n = \Delta h = \frac{(n\_{mb} + h\_{mb}) - \sqrt{(n\_{mb} - h\_{mb})^2 + 4n\_0h\_o}}{2} \tag{31}
$$

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

$$-\vec{n} \cdot \vec{\Gamma}\_{n\_{\text{obs}}} = \vec{n} \cdot \vec{j}\_{\text{e}} \tag{32}$$

$$-\vec{n} \cdot \vec{\Gamma}\_{h\_{\text{obs}}} = \vec{n} \cdot \vec{j}\_{\text{t}} \tag{33}$$

where ji and je are 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:

$$
\vec{n} \cdot \vec{\Gamma}\_{n\_{ub}} = \vec{n} \cdot (-\mu\_e n\_{ub} \cdot \nabla V) \tag{34}
$$

$$
\vec{n} \cdot \vec{\Gamma}\_{h\_{nb}} = \vec{n} \cdot (-\mu\_h h\_{mb} \cdot \nabla V) \tag{35}
$$

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].

#### 2.4 Poisson equation

papers [17, 18]. The initial electron density is set to 1 � <sup>10</sup><sup>9</sup> 1/m<sup>3</sup> and the O2 to N2

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 distri-

where τ is the non-equilibrium carriers' lifetime (s); R3, R2 and R1 are the recom-

Inspired by the theories in the semiconductor physics on the non-equilibrium, the products of the hole density h<sup>0</sup> and mobile electron density are considered to be

constants in the insulator within the condition of the thermal equilibrium.

densities at conduction band bottom and the valence band top (1/m<sup>3</sup>

where Eg is the insulator's energy gap (eV); Nc and Nv are states' effective

the hole density and electron are hmb and nmb in an unbalanced condition, Δh and Δn is able to be obtained through working out Eq. (29). The procedure of the electron–hole recombination is represented by positive Δh and Δn while the electron-hole pairs' generation is represented by negative Δh and Δn.

hole and electron (1/s); nmb is the mobile electron density (1/m3

<sup>h</sup> and Jc

hole and electron's diffusion coefficient (m2

); Jc

/s); Pde and Ptr are the detrapping and trapping probability of

); hmb is the mobile hole density (1/m3

/s); V is potential (V); ntr is the trapped

<sup>e</sup> are the mobile electron and hole's flux (1/(m2 s));

ð23Þ

ð24Þ

ð25Þ

ð26Þ

ð27Þ

ð28Þ

ð29Þ

ð30Þ

). Supposing

); μ<sup>h</sup> and μ<sup>e</sup> are the

).

ratio is 1:4. Δn<sup>e</sup> means the electron increment of each reaction.

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

bution along the surface of the insulator.

bination coefficients (m3

electron density (1/m3

56

htr is trapped hole density (1/m3

2.3 Detrapping process and charge trapping in the epoxy insulator

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

$$-\nabla^2 V = e(\sum\_{1}^{p} n\_k z\_k) / \,\varepsilon\_0 \tag{36}$$


Table 2. Some parameters of epoxy in the trapping-detrapping model. where species k's density is nk (1/m<sup>3</sup> ). We can describe Poisson equation within insulator as follows:

$$-\nabla^2 V = e(h\_{m\flat} + h\_{\wr} - n\_{m\flat} - n\_{\wr}) / (\varepsilon\_0 \varepsilon\_{\r}) \tag{37}$$

where ε<sup>r</sup> is the insulator's relative dielectric constant.

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

$$\mathcal{V} = \mathcal{V}\_0 - I\_p \mathcal{R}\_b \tag{38}$$

$$I\_p = -\int (\vec{n} \cdot \vec{j}\_i + \vec{n} \cdot \vec{j}\_e)dS \tag{39}$$

where V<sup>0</sup> is the voltage of power supply (V); Ip is the discharge current (A), which can be obtained by integrating the current density on the boundary of the HV electrode; V is the potential on needle electrode (V).

The open boundary is as follows:

$$
\vec{n} \cdot \nabla V = 0\tag{40}
$$

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. Electron density's variation is shown in Figure 4. The electron density's variation

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

The simulated distribution of electric potential and force lines.

DOI: http://dx.doi.org/10.5772/intechopen.80635

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

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

. It is clear that the electron avalanche grows

accumulation.

59

Figure 3.

reached approximately 1 <sup>10</sup><sup>15</sup> 1/m<sup>3</sup>

#### 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.

#### 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

Figure 2. A schematic diagram of the surface potential measurement system.

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

Figure 3.

where species k's density is nk (1/m<sup>3</sup>

).

We add a resistor Rb for the purpose of limiting the discharge current Ip. The

where V<sup>0</sup> is the voltage of power supply (V); Ip is the discharge current (A), which can be obtained by integrating the current density on the boundary of the HV

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 measure-

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

ð37Þ

ð38Þ

ð39Þ

ð40Þ

We can describe Poisson equation within insulator as follows:

where ε<sup>r</sup> is the insulator's relative dielectric constant.

definition of the voltage on HV electrode is as follows:

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

electrode; V is the potential on needle electrode (V).

ment of 5 points' values from the center to the insulator's edge.

A schematic diagram of the surface potential measurement system.

The open boundary is as follows:

2.5. Surface potential measurement

3. Results and discussion

Figure 2.

58

The simulated distribution of electric potential and force lines.

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.

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.

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 <sup>10</sup><sup>15</sup> 1/m<sup>3</sup> . It is clear that the electron avalanche grows

#### Figure 4.

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

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 Efield strength at needle tip is presented in Figure 6b's margin for the purpose of

The distribution of electron and ion density along the symmetry axis: (a) 1 <sup>10</sup>–<sup>9</sup> s, (b) 1.8 <sup>10</sup>–<sup>8</sup> s,

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

DOI: http://dx.doi.org/10.5772/intechopen.80635

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

seeing the outcome obviously.

Figure 5.

61

(c) 4.6 <sup>10</sup>–<sup>7</sup> s, (d) 60 s.

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.

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

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

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 Efield 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

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

The variation of electron density with the development of corona discharge: (a) 1 <sup>10</sup><sup>9</sup> s, (b) 1.8 <sup>10</sup><sup>8</sup> s,

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.

reduces with time, which is presented in Figure 5d.

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

Figure 4.

60

(c) 4.6 <sup>10</sup><sup>7</sup> s, (d) 60 s.

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.

Figure 8.

63

electron and (b) trapped electron.

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

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

DOI: http://dx.doi.org/10.5772/intechopen.80635

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

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.

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

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

Figure 7.

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

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

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

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

escape at the time of transporting to the ground.

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

Figure 6.

62

center will spend more time to reach saturation.

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

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

Figure 9.

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

4. Conclusions

Figure 10.

65

major conclusion is shown below:

reaches saturation.

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

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.

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

DOI: http://dx.doi.org/10.5772/intechopen.80635

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

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

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.

surface potential takes a shorter time to reach saturation.

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.

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.

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

#### Figure 10.

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. 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.

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

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

Figure 9.

64

(b) t = 60 s.

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:


Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

References

1754-1763

2602-2610

2012;19(3):1084-1090

[1] Kumara S, Alam S, Hoque IR, et al. DC flashover characteristics of a polymeric insulator in presence of surface charges. IEEE Transactions on Dielectrics and Electrical Insulation.

DOI: http://dx.doi.org/10.5772/intechopen.80635

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

model with a constant mobility: Model prediction. Journal of Physics D: Applied Physics. 2003;37(2):298

[9] Min D, Li S. Simulation on the influence of bipolar charge injection and trapping on surface potential decay of polyethylene. IEEE Transactions on Dielectrics and Electrical Insulation.

[10] Chen G. A new model for surface potential decay of corona-charged polymers. Journal of Applied Physics.

2014;21(4):1627-1636

2010;43(5):55405-55411

2017;23(6):3742-3751

2012;51(51):6001

1995;28(4):738

47(8):6459-6463

[11] Zhou F, Li J, Yan Z, et al. Investigation of charge trapping and detrapping dynamics in LDPE, HDPE and XLPE. IEEE Transactions on Dielectrics and Electrical Insulation.

[12] Min D, Li S, Li G. The effect of charge recombination on surface potential decay crossover characteristics

of LDPE. IEEE International Symposium on Electrical Insulating Materials (ISEIM). 2014. pp. 104-107

[13] He W, Liu XH, Yang F, et al. Numerical simulation of direct current glow discharge in air with experimental validation. Journal of Applied Physics.

[14] Nahomy J, Ferreira CM, Gordiets B, et al. Experimental and theoretical investigation of a N2-O2 DC flowing glow discharge. Journal of Physics D.

[15] Choi J, Matsuo K, Yoshida H, et al.

atmospheric pressure air microplasma jet. Journal of Applied Physics. 2008;

[16] Tochikubo F, Arai H. Numerical simulation of streamer propagation and

Characteristics of a DC-driven

[2] Kumara S, Serdyuk YV, Gubanski SM. Simulation of surface charge effect on impulse flashover characteristics of outdoor polymeric insulators. IEEE Transactions on Dielectrics and Electrical Insulation. 2010;17(6):

[3] Du BX, Li ZL. Surface charge and DC flashover characteristics of directfluorinated SiR/SiO2 nanocomposites. IEEE Transactions on Dielectrics and Electrical Insulation. 2015;21(6):

[4] Nakanishi K, Yoshioka A, Arahata Y, et al. Surface charging on epoxy spacer at DC stress in compressed SF6 gas. IEEE Transactions on Power Apparatus and Systems. 1983;102(12):3919-3927

[5] Sato S, Zaengl WS, Knecht A. A numerical analysis of accumulated surface charge on dC epoxy resin spaces.

[6] Mohamad A, Chen G, Zhang Y, et al. Mechanisms for surface potential decay on fluorinated epoxy in high voltage DC applications. IEEE Conference on Electrical Insulation and Dielectric Phenomena (CEIDP). 2014. pp. 863-866

[7] Du BX, Zhang JW, Gao Y. Effect of nanosecond rise time of pulse voltage on the surface charge of epoxy/TiO2 nanocomposites. IEEE Transactions on Dielectrics and Electrical Insulation.

[8] Roy SL, Segur P, Teyssedre G, Laurent C. Description of bipolar charge transport in polyethylene using a fluid

2013;20(20):321-328

67

IEEE Transactions on Electrical Insulation. 1987;22(3):333-340

### Author details

Boxue Du, Hucheng Liang and Jin Li\* Key Laboratory of Smart Grid of Education Ministry, School of Electrical and Information Engineering, Tianjin University, Tianjin, China

\*Address all correspondence to: lijin@tju.edu.cn

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

Simulation on the Surface Charge Behaviors of Epoxy Insulator by Corona Discharge DOI: http://dx.doi.org/10.5772/intechopen.80635

#### References

[1] Kumara S, Alam S, Hoque IR, et al. DC flashover characteristics of a polymeric insulator in presence of surface charges. IEEE Transactions on Dielectrics and Electrical Insulation. 2012;19(3):1084-1090

[2] Kumara S, Serdyuk YV, Gubanski SM. Simulation of surface charge effect on impulse flashover characteristics of outdoor polymeric insulators. IEEE Transactions on Dielectrics and Electrical Insulation. 2010;17(6): 1754-1763

[3] Du BX, Li ZL. Surface charge and DC flashover characteristics of directfluorinated SiR/SiO2 nanocomposites. IEEE Transactions on Dielectrics and Electrical Insulation. 2015;21(6): 2602-2610

[4] Nakanishi K, Yoshioka A, Arahata Y, et al. Surface charging on epoxy spacer at DC stress in compressed SF6 gas. IEEE Transactions on Power Apparatus and Systems. 1983;102(12):3919-3927

[5] Sato S, Zaengl WS, Knecht A. A numerical analysis of accumulated surface charge on dC epoxy resin spaces. IEEE Transactions on Electrical Insulation. 1987;22(3):333-340

[6] Mohamad A, Chen G, Zhang Y, et al. Mechanisms for surface potential decay on fluorinated epoxy in high voltage DC applications. IEEE Conference on Electrical Insulation and Dielectric Phenomena (CEIDP). 2014. pp. 863-866

[7] Du BX, Zhang JW, Gao Y. Effect of nanosecond rise time of pulse voltage on the surface charge of epoxy/TiO2 nanocomposites. IEEE Transactions on Dielectrics and Electrical Insulation. 2013;20(20):321-328

[8] Roy SL, Segur P, Teyssedre G, Laurent C. Description of bipolar charge transport in polyethylene using a fluid

model with a constant mobility: Model prediction. Journal of Physics D: Applied Physics. 2003;37(2):298

[9] Min D, Li S. Simulation on the influence of bipolar charge injection and trapping on surface potential decay of polyethylene. IEEE Transactions on Dielectrics and Electrical Insulation. 2014;21(4):1627-1636

[10] Chen G. A new model for surface potential decay of corona-charged polymers. Journal of Applied Physics. 2010;43(5):55405-55411

[11] Zhou F, Li J, Yan Z, et al. Investigation of charge trapping and detrapping dynamics in LDPE, HDPE and XLPE. IEEE Transactions on Dielectrics and Electrical Insulation. 2017;23(6):3742-3751

[12] Min D, Li S, Li G. The effect of charge recombination on surface potential decay crossover characteristics of LDPE. IEEE International Symposium on Electrical Insulating Materials (ISEIM). 2014. pp. 104-107

[13] He W, Liu XH, Yang F, et al. Numerical simulation of direct current glow discharge in air with experimental validation. Journal of Applied Physics. 2012;51(51):6001

[14] Nahomy J, Ferreira CM, Gordiets B, et al. Experimental and theoretical investigation of a N2-O2 DC flowing glow discharge. Journal of Physics D. 1995;28(4):738

[15] Choi J, Matsuo K, Yoshida H, et al. Characteristics of a DC-driven atmospheric pressure air microplasma jet. Journal of Applied Physics. 2008; 47(8):6459-6463

[16] Tochikubo F, Arai H. Numerical simulation of streamer propagation and

Author details

66

Boxue Du, Hucheng Liang and Jin Li\*

Key Laboratory of Smart Grid of Education Ministry, School of Electrical and

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

Information Engineering, Tianjin University, Tianjin, China

Atmospheric Pressure Plasma ‐ From Diagnostics to Applications

\*Address all correspondence to: lijin@tju.edu.cn

provided the original work is properly cited.

radical reactions in positive corona discharge in N2/NO and N2/O2/NO. Journal of Applied Physics. 2002;41 (2A):844-852

[17] Itikawa Y, Ichimura A. Cross sections for collisions of electrons and photons with atomic oxygen. Journal of Physical and Chemical Reference Data. 1990;19(3):637-651

[18] Itikawa Y. Cross sections for electron collisions with nitrogen molecules. Journal of Physical and Chemical Reference Data. 2006;34(38): 1-20

[19] Zhou T, Chen G, Liao R, Xu Z. Charge trapping and detrapping in polymeric materials: Trapping parameters. Journal of Applied Physics. 2009;106(12):644-637

[20] Chen G, Zhao J, Zhuang Y. Numerical modeling of surface potential decay of corona charged polymeric material. IEEE International Conference Solid Dielectrics (ICSD). 2010. pp. 1-4

**69**

Section 2

Plasma Applications

## Section 2
