Estimation of photo-degradation of dielectrics surrounding the narrow channel due to PD activity

Partial discharge (PD) taking place within narrow channels produces a large number of charged particles, causing degradation of the polymer. Additionally, PD pulses also produce UV photons due to decay of excited states (radiative states) in air. These may have enough energy to break C–C bonds and thereby add to the degradation of the dielectrics surrounding the narrow channel. In this paper, a radiation transport (RT) model has been developed and integrated with a Particle-in-Cell/Monte Carlo collision (PIC-MCC) model to study the behavior of excited (radiative) states of air within discharge in the narrow channel. The radiative state atoms are described by a fluid model combined with the Holstein–Biberman equation. This model has the ablation to follow the spatial evolution of the radiative excited-state density throughout narrow channel. The effect of applied electric field, narrow channel dimensions, gas pressure on the extent of degradation of dielectrics surrounding the narrow channel is studied.


Introduction
Partial discharge (PD) pulses occurring within the narrow channels are responsible for causing degradation of the surrounding dielectric material. Degradation is defined as a process which leads to a deterioration of any physical property of a polymer [1,2]. Degradation of the narrow channel walls is a complex process, and some of principal mechanisms involved are mechanical erosion, thermal and chemical decomposition of the polymeric material.
Laurent et al. [3,4] investigated the growth of electrical trees through the simultaneous measurement of light intensity and PD. A photomultiplier tube (PMT) was used to detect light emission at the tip of a needle electrode embedded in LDPE. They studied radiation emission in insulating materials (LDPE) used for cables. Further, they [5] evaluated the gases appearing in the channels of bush and branch-like electrical trees in LDPE and cross-linked polyethylene (XLPE) when subjected to a strong alternating electric field. They report that the emission spectra obtained confirm the presence of hydrogen and carbon oxides.
A quantitative study was done by Champion et al. [6] on light emission during initiation and the early growth stages of electrical treeing in synthetic resins. The three types of light emission observed by them are as follows: (a) electroluminescence caused by charge injection based on Fowler-Nordheim or Schottky effect dependent on local electric field, (b) micro-discharge activity along with formation of narrow channels and (c) conventional PD activity within narrow channel. Champion et al. [7] found that under long-term continuous electrical stress, light emission is found to spread principally along the cone of the pin rather than into the resin. In another work [8,9], they observed the spatial distributions of light emission in various types of trees (bush, bush/branch, branch).
Kaneiwa et al. [10] attempted to understand the mechanism of degradation through electrical treeing and tree propagation using an artificial tree channel in LDPE, from which tree propagation occurred. It was seen that the region of light emission reaches the vicinity of the channel tip before tree propagation. In those samples where tree did not propagate, the region of light emission due to PD does not reach the channel tip. They concluded that tree propagation occurs only when large magnitude PDs are present; long luminous regions occur within the artificial channel and damage the channel tip before tree initiation. Kaneiwa et al. [11] further observed that for positive PD, the light emission is near the needle-tip and a positive correlation is to be found between the PD magnitude and the length of luminous region. They concluded that PD starts at the needle tip and the PD magnitude is determined by the PD length. For negative PD, there appeared to be no correlation between the PD magnitude and the length of luminous region.
Wu et al. [12][13][14][15] reported that in a long-narrow channel, the behavior of PD is similar to that in an actual tree and the range of light emission in the channel increases when the applied voltage increases.
The main source of light emission within air-filled narrow channels during PD activity is emission from the second positive system (2P) of N 2 (C 3 P u ? B 3 P g ) [16]. They are UV photons in the range of 300 (4.13 eV) to 406 nm (3.05 eV). As mentioned earlier, the UV light produced inside the channel is due to the decay of excited states and can cause photo-degradation of the surfaces surrounding the narrow channel. It would be expected that the predominant products of radiolysis of linear alkanes would be due to the rupture of C-C bonds [17][18][19]. This process would be particularly sensitive to electric fields if charged species such as radical ions are formed as intermediates, or if they bring about changes in the dielectric constant by altering the concentrations of polar species [17][18][19].
It is obvious that light emission plays an important role in PD pulses. As far as we are aware, there is no theoretical model to simulate generation and propagation of photons within narrow channels, and their subsequent effect on the surrounding dielectrics.
In previous works [20], a two-dimensional Particle in Cell-Monte Carlo Collision (PIC-MCC) model [21][22][23] to study the effect of similarities and dissimilarities of electrodes on the spark-type PD pulses in sub-millimeter sized cavities was developed. Moreover, the effect of applied electric field on the damage due to partial discharges activity into the surroundings dielectrics of a narrow channel encapsulated within the volume of a dielectric material was evaluated using this model [24]. Based on this model, the effects of different parameters of PD activity in the narrow channel on the damage into the surroundings dielectrics were studied [25]. In this paper, the above PIC-MCC model coupled with a Radiation Transport (RT) model [26][27][28][29][30][31][32][33][34][35][36][37][38][39] is used to describe the production of UV radiation inside narrow channel and its interaction with surrounding dielectrics. This RT-PIC model incorporates the resonance trapping effect of photons, to the kinetic model of the charged particles when a PD pulse develops inside a narrow channel. This model allows us to follow the spatial evolution of the density of radiative excited state (photons) within the narrow channel. The effect of applied electric field and micro-channel length on the extent of degradation of dielectric anode, cathode and side-walls, based on C-C bonds breaking due to UV photons impact, is studied.

Simulation model
Computational method for PD within narrow channel In the PIC scheme, particles are defined in a continuum position and velocity space. Field values are defined at discrete locations in space. Particle and field values are advanced sequentially in time, starting from initial conditions. The particle equations of motion are solved at every time step, using field values interpolated from the discrete grid to particle locations. The force on every particle is computed by interpolation of the field values from the grid position to the given particle positions. The position and velocities of each particle is next updated based on the solution of the classical equation of motion [40]. Next, particle boundary conditions are applied. For modeling collisions, the Monte Carlo collision (MCC) scheme is applied [41]. Source terms for the field equations are accumulated from the particle locations to the grid locations. The field values are then advanced by one time-step, and the time loop starts again (Fig. 1).
In the PIC-MCC scheme, each super-particle, representative of one or a much larger number of real particles, is designated in continuum space by its position and velocity. Initially, a Maxwellian distribution is assumed for electrons and ions. Further, the particles are uniformly distributed in the gap space. Particle boundary conditions such as absorption and emission are used to account for the relation between the discharge current in the gap and the current in the external circuit. When an electron or ion passes from the discharge into an end wall, it adds to the wall charge and is deleted from the list of active particles [23]. Additionally, secondary emissions occur when a charged particle impacts a surface (metallic or dielectric) causing ejection of electrons from the surface.
In the current work, the gas within the cavity is assumed to be air at atmospheric pressure. Therefore, photo-emission at the cathode is likely to be the dominant secondary emission process, rather than electron emission from cathode or anode due to ion or electron bombardment. However, all of the above processes have been accounted for in the model. After emission of a new electron, its energy is also assigned from an assumed Maxwell-Boltzmann distribution. The interaction of charged particles with neutral atoms and molecules, and other collisional processes are included by using a Monte Carlo Collision technique [42]. Three main charged particle species, viz. electrons, positive ions and negative ions have been considered. The major electron-neutral molecule impact reactions included are (i) elastic, (ii) excitation, (iii) ionization (including all important ionization reactions) and (iv) electron attachment. Related cross-sectional values are extracted from experimental data available in the literature [29]. Integral cross-section values for electron interaction in dry air have been determined from its molecular constituents by assuming a target composition of 78.09 % N 2 and 21.19 % O 2 [43]. Only elastic and charge transfer collisions between ions and neutrals are considered. Because ion-molecule cross-sections are not readily available, a constant cross-section (10 -20 m 2 ) has been assumed based on a hard-sphere collision model [44]. A significant contribution of this work is that the discharge within the cavity is simulated in conjunction with an external circuit comprising of a voltage source, and a blocking capacitor that stops the flow of any average current. The coupling between the charge flow within the micro-cavity and the external current is obtained via the parallel-plane electrodes, which may be metallic or dielectric, depending on the location of the discharge site.
The potentials and fields are obtained using Poisson's equation: r Á e 0 rUðx; tÞ ¼ qðx; tÞ ð 1Þ The potential (A) can be separated as follows [42,45,46]: where A P and A L represent the Poisson and Laplacian parts of electrical potential (A). Therefore, the field may be described by the combination of the following equations: The boundary condition for Eq. (3) is A = 0 on all boundaries as the Poisson field is solely due to charge in the medium. For each boundary with a Dirichlet condition, Eq. (4) is solved for A i = 0 on the equipotential surface, and A = 0 elsewhere which gives us A Li (the potential due to the ith Laplacian field). A L is obtained by the superposition of all the Laplacian fields. Neumann boundary conditions are included through A P = 0. This method neglects charges induced by a driven electrode on other boundaries which are connected to an external circuit. It is also possible to solve the field Poisson equation with boundaries and circuits [40,42,47].
In this work, identical cavities, as might occur at discharge sites, with two different electrode configurations have thus been considered. These are the two generic cases described below. The model is two-dimensional with Neumann boundary conditions on the side-walls. Dirichlet boundary conditions (potential drop specified) are assumed on the electrodes.
A narrow channel bounded by dielectrics on both sides is depicted in Fig. 2a. The equivalent circuit of this configuration is shown in Fig. 2b, where C g , C d1 , C d2 and C r represent the equivalent capacitances of the gap, the dielectric layers and dielectric material parallel to the micro-channel, respectively. The discharge current may be represented by I conv . Using KCL, it is possible to obtain the following relations [21]: where I ext and r T are the total external current and charge on the dielectric surface. I r is the leakage current flowing in the bulk of dielectric material parallel to the narrow channel. By combining the above equations, Equation (27) in ref. [21] may be rewritten as follows: where a = (1 ? 2(C ? C r )/C d ) -1 . r T , A, C, V(t), Q c , Q r , Q conv , q, u ij and N j are respectively the total surface charge density on the electrodes at a given time 't', cross-sectional area of electrode, capacitance of blocking capacitor, driving voltage, charge in the external circuit, charge in the bulk of dielectric material parallel to the narrow channel, charge convection due to discharge, space charge density, the electric potential at each grid point and the number of grid divisions along the y-axis. The derivation of Eq. (9) is presented in ''Appendix 1''.

Model for radiation transport
The PIC-MCC model for gas discharge within narrow channel is described in detail in ref. [20,24,25]. In this section, we describe the radiation transport model for photons dynamics. For both N 2 and air, the excitation coefficient is assumed to be that representing the second positive system of N 2 [48]. Hence, in this study, we consider the radiation trapping of the photons emitted from the second positive system (2P) of N 2 (C 3 P u ? B 3 P g ) [48,49]. The radiation transport equation is included with the PFM [34] and the neutral gas density profile is assumed to be uniform, which is valid for the pressure regime concerned. Because the cell is operated at atmospheric pressure, broadening with a Lorentz line shape is dominant. It has been assumed that the frequency of a re-emitted photon is completely independent of that of the absorbed photon.
The modified Holstein equation for the resonant (radiative) state density is as follows [39]: where n r (r,t) and P r (r,t) are the radiative excited state density and their production rate inside the micro-channel. s t is the vacuum radiative decay time, which is around 40 ns for air (N 2 (C 3 P u ? B 3 P g )) [40], G(r, ŕ) is the kernel function, which is the probability for a resonance photon emitted at the position ŕ to be absorbed at position r. The calculation of G(r, ŕ) is presented in ref. [39]. The total radiation flux at the dielectric anode and cathode surface is: and at dielectric side-walls: Derivation of above relations is presented in ''Appendix 2'', ''Appendix 3'' and ''Appendix 4'' [39].
Photo-degradation of surrounding dielectrics due to UV photons impact As mentioned above, the excitation coefficient of radiative state for air is assumed to be that representing the second positive system of N 2 (C 3 P u ? B 3 P g ) [50]. The experimental results reported in literature show that the photon yield from excited electrons of air is between 300 and 406 nm which is equivalent to photon energies between 3.05 and 4.13 eV (UV) [16].
Based on experimental results presented in Fig. 3, around 65.34 % of photons produced due to decay of radiative states of air (nitrogen) have energies more than 3.5 eV [51]. These photons have the ability to break C-C bonds and thus cause degradation in dielectric. Absorption of radiation is an essential first step to photo-degradation. Strongly absorbed radiation will be attenuated as it passes through the polymer and reaction will be concentrated in the surface layers. The first chemical step in photo-degradation is usually hemolytic bond scission to form free radicals.
Diglycidyl ether of bisphenol-A (DGEBA) is a typical, commercial epoxy resin and is chosen, in this work, as the base material to illustrate the mechanism of damage computation. Figure 4 shows the molecular structure of DGEBA.
Thus, each DGEBA monomer has '14' C-C bonds. The radius of a DGEBA monomer is approximately R (&13.28 Å ) [52], and its cross-sectional area is approximately A DGEBA = 5.546 9 10 -18 m 2 . The cross-sectional area of the dielectric anode is A anode = A cathode = L y 9 L z and the area of each side-wall is A side-walls = L x 9 L z . The number of DGEBA monomers on the surface of the dielectric anode and cathode is thus approximately (N DGEBA ) anode = (N DGEBA ) cathode = A anode /A DGEBA , and that on the surface of side-wall is (N DGEBA ) side-walls = A side-walls /A DGEBA . If N CC is number of C-C bonds which can be broken during one PD pulse, then the total pulses required to break all C-C bonds existing on the dielectric anode, cathode and side-walls (N pc ) are as follows: Here (N pc ) anode , (N pc ) cathode and (N pc ) side-walls are the total pulses required to break all C-C bonds existing on the surface of dielectric anode, cathode and side-walls.
The total number of cycles required to break all C-C bonds on surrounding dielectrics of narrow channel may be obtained as follows: Here (N cycle ) anode , (N cycle ) cathode and (N cycle ) side-walls are the total number of required cycles to break all C-C bonds on the surface of dielectric anode, cathode and side-walls, where N PDC indicates the number of PD pulses per cycle [53][54][55][56][57].
So, the time required to break all C-C bonds on the surface of surrounding dielectrics of narrow channel (t cc ) at a supply frequency of 50 Hz (T = 20 ms) can be calculated as follows: Here (t cc ) anode , (t cc ) cathode and (t cc ) side-walls are the times required to break all C-C bonds on surface of dielectric anode, cathode and side-walls.

Simulation results
The simulation results of the two-dimensional RT-PIC model are presented here. The results of the effect of several factors like the applied electric field, channel dimensions and gas pressure on the extent of degradation of dielectric anode, cathode and side-walls, due to C-C bond break by UV photon impact, are presented in the following sections. Results are presented here with simulation parameters presented in Table 1. It may be noted that a channel with a significantly large length/width ratio is chosen.
The scatter plots in phase space for velocities in xdirection i.e. (V x , x, y) for electrons, and positive ions inside the micro-channel are shown in Figs. 5 and 6 at particular instants of time. As shown here, drift velocity of electrons is much greater than that of ions (3 orders of magnitude higher) and quickly traverse the micro-channel length. All charged particles are distributed about the center-line y = 50 lm. Most of charged particles are created very close to the anode. The number of ions near the anode increases through ionization with passage of time. The positive ions slowly drift towards the cathode, till at about 600 ps the entire gap distance is almost evenly populated. From this time onwards, ions gradually leave the channel, thus depleting the ion density in the gap. This is similar to the evolution of the electrons within the gap, as seen in Fig. 5, except that electrons move faster and need travel only a small distance between the ionization region and the anode. As is clear from Fig. 6, ion velocities increase as they reach near the cathode, accelerated by the electric field. Generally, owing to high electric field strength along micro-channels, attachment process is weak and the production of negative ions is low.
The spatial variation of the x-component of electric field (E x (x, y)) is depicted in Fig. 7. The distribution of field within the narrow channel is originally uniform. The distribution of charged particles inside the narrow channel (in Figs. 5,6) shows that the effect of space charge on applied field is mostly due to positive ions. A positive space charge region forms ahead of the anode where they accumulate. This results in a reduced field close to the anode. Simultaneously the field near the cathode is enhanced in the direction of the applied field, due to the positive charge cloud now nearing the negatively charged cathode. Owing to the accumulation of negatively charged particles at the anode which causes the anode potential to reduce, the field at the anode never recovers original values. Finally, the field within the channel becomes uniform as the charge dissipates.

Applied electric field
In this subsection, the effect of applied electric field on the extent of degradation of dielectrics surrounding the narrow channel based on C-C bonds broken by PD pulses inside the narrow channel is studied. As electric field across narrow channel increases, ionization increases and consequently there is an increase in the number of excited states (radiative states) inside the narrow channel. Thus not only the total number of electrons, but additionally the number of UV photons reaching the surrounding dielectrics increases. As described before, every DGEBA monomer has '14' C-C bonds; photons and electrons with energies higher than 3.5 eV are capable of breaking these molecular bonds. The number of broken C-C bonds during a single PD pulse (N CC ) by electrons and photons are depicted in Fig. 8. Figure 8a shows the impact of electrons is higher on the dielectric anode than on the side-walls, while photons have a greater effect on sidewalls than on anode, as seen in Fig. 8b. This is because in the current simulation, the ratio of length to width (20) is high and this causes a larger number of photons to reach the dielectric side-walls rather than the anode or cathode. Generally, propagation of photons inside the micro-channel is not determined by the applied electrostatic field due to their mass-less and charge-less nature, while charged particle movement is governed by fields. Thus photons are likelier to reach the side-walls rather than the ends of the narrow channel, as that entails lesser distance of travel. Therefore, the side-walls are overwhelmingly likely to undergo photo-degradation. Interestingly, as shown in Fig. 8, the value of N CC due to UV photons is almost 1,000 higher than due to electrons. Hence, the resultant degradation based on C-C bond-break due to the effect of electrons is negligible.
As seen in Fig. 9a, the time required to break all the existing C-C bonds (t cc ) on the surface of dielectrics surrounding the narrow channel varies considerably with electric field. An increase of electric field from 8 to 50 kV/ mm causes t cc to vary from 347 to 0.2121 ms for the dielectric anode, 248 to 0.1102 ms for the dielectric sidewalls and 618 to 0.4537 ms for dielectric cathode. Figure 9b through d show the temporal variation in the number of UV photons reaching the dielectric anode, cathode and side-walls for different values of applied electric field. We see a higher number of UV photons at higher fields. Tanaka et al. [59] have shown that for a cavity with a diameter of 5 mm, the intensity of light emitted from a discharge within the cavity decreases with decrease in applied voltage.
It may be noted that the light emission from PD activity gives us information about spatial distribution of PD in Fig. 6 Scatter plot of positive ions in the phase-space (V x , x, y) inside the micro-channel at particular instants of time narrow channels. Figures 10, 11, 12 and 13 show that the radiative states (or photons) inside the narrow channel are mostly concentrated near the anode (ionization region). We see that with increase in electric field, the region of production and concentration of photons moves further towards the dielectric anode.
Laurent et al. [5] stated that because of the very small sizes of narrow channels, the results of the discharge physics for plane, parallel, metallic electrodes may not be easy to predict. It might not be possible to know definitely whether discharges would develop between or along the walls of the channel [60].
Our results corroborate the observations of Kaneiwa et al. [10]. They noted that light emission occurred due to PD in an artificial channel before a new narrow channel propagated. At lower applied voltage of 3 kV, light emission mainly occurs near the needle tip. At higher applied voltage of 6 kV, the light emission extended to the vicinity of channel tip.
Champion et al. [61] developed a partial discharge model (a pin-plane electrode arrangement) to simulate experimental partial discharge activity within electrical tree (narrow channel) structures grown in polymer based insulation. They also simulated partial discharges within an arbitrarily defined non-conducting electrical tree structure [62]. They report that the light emitted from partial discharge activity was initially observed from the whole tree structure. The photographs of light emission due to a PD pulse in an artificial tree channel of 5 mm length reported by Wu et al. [12] show that the region of light emission becomes larger at higher voltages. They report that the light is brighter near the needle tip and the region of light emission is much shorter than the channel length, i.e. light emission takes place only in a part of the channel closer to the needle electrode.

Micro-channel length
In this section, the effect of narrow channel length on the degradation of surrounding dielectrics due to PD pulse inside the narrow channel based on the C-C bond break is considered. Simulations are performed for the operating parameters presented in Table 1. As a single electron traverses larger distances between the dielectric electrodes in larger narrow channels, with increase in narrow channel length (applied electric field remaining the same), greater amounts of excitation, ionization and attachment occur due to larger number of collisions per electron. Increase in number of excitation collisions leads to increase in the total number of radiative states and consequently the number of UV photons reaching the dielectrics surrounding the narrow channel. Thus, as seen in Fig. 14a, with increase in length the time required to break all existing C-C bonds (t cc ) reduces. Noskov et al. [63] developed a self-consistent model for electrical tree propagation and PD activity within narrow tree channels and showed that the magnitude of individual PD event increases with the narrow channel length.
Temporal variations in the total number of UV photons reaching the dielectric anode, cathode and side-walls for specific narrow channel lengths are presented in Fig. 14b through d. As seen in these figures, the number of UV photons reaching the surrounding dielectrics is higher at greater channel lengths. The total number of UV photons arriving at the surface of dielectric side-walls is an order of magnitude higher than at the dielectric anode and two orders of magnitude greater than at dielectric cathode.
As seen in Figs. 15, 16 and 17, as the narrow channel length increases, the region where the UV photons are mostly produced moves towards the dielectric anode. Also the size of this region increases with channel length, and the photons are distributed over a wider area within the narrow channel. As may be seen, for the channel length of L x = 10 lm, the order of radiative states (photons) is 10 -18 m -3 and is mostly concentrated around x/L x = 0.3, while for greater narrow channel lengths, i.e. L x = 25 lm, the order of radiative states is 10 -20 m -3 and with a wider distribution which is concentrated around x/L x = 0.04.

Width
The effect of width is studied in this subsection. The simulation is performed for the operating parameters in Table 1.
With increase in width of narrow channel, the likelihood of random movement of charge carriers in a direction perpendicular to the applied electric field (y-direction) increases [64]. Hence, the rate of ionization, as the charges move under the influence of the electric field, is lower [20]. The likelihood of excitation reduces. Due to decrease in the number of exciting collisions, the density of radiative states and consequently the number of UV photons reaching the surrounding dielectrics decreases. Thus, as is depicted in  Figure 18a shows that for channel widths below 2 lm, the value of t cc for dielectric side-walls is lower than for dielectric cathode. Thus, the effect of width appears to be similar to the effect of electric field, as seen before. But above 2 lm, the value of t cc for dielectric side-walls is higher. As mentioned earlier, the movement of photons is not affected by applied electrostatic field, and they are likely to travel radially outwards from the site of production. In doing so, at greater narrow channel widths, the probability of such photons reaching the side-walls rather than the anode or even the cathode is lower. This would explain the above observation.
The effect of reduced excitation is lower production in radiative states (UV photons) in the narrow channel for wider narrow channels, from 10 -19 m -3 for L y = 2 lm to 10 -15 m -3 for L y = 80 lm, as seen in Figs. 19, 20 and 21. A comparison of Figs. 19 and 20 shows that photons are concentrated mostly around x/L x = 0.2 for L y = 2 lm and x/L x = 0.6 for L y = 80 lm inside the narrow channel. It may be concluded that the region of concentration of UV photons is seen to be located far away from the anode at greater narrow channel widths.

Pressure
Generally, gas pressure inside the narrow channel has considerable effect on the light emitted from gaseous discharges. In this subsection we investigate the effect of gas pressure inside the narrow channel. The simulation is performed for the operating parameters presented in Table 1.
The ionization process and consequent production of new electrons in partial discharges is affected by change in pressure of gas due to its direct effect on the probability of collision of charged particles with neutrals. Just as an increase in gas pressure causes greater ionization of neutrals, so does it cause a number of excited (radiative) states or UV photons to be produced inside the narrow channel and thus a larger number reaches the surrounding dielectrics. As a result, as shown in Fig. 22a, with increase in gas pressure inside the narrow channel, the time required to break all the existing C-C bonds (t cc ) by photon impact reduces. Figure 22b through d show the temporal variations of the total number of the UV photons reaching the dielectric anode, cathode and side-wall surfaces for three different values of gas pressure inside narrow channel. The total number of UV photons reaching the dielectrics surrounding the narrow channel is higher at higher pressures.
As seen in Figs. 23, 24 and 25, region of concentration of radiative states or UV photons inside the narrow channel occurs nearer the anode at higher pressures. Also, the UV photons inside the narrow channel are distributed over a more extensive region at higher pressures.

Discussion
From the results presented in this paper, we see that the impact of electrons on dielectric anode is significantly more than on side-walls. In contrast, the UV photons impact predominantly on the dielectric side-walls.
We have shown that the UV photons are mainly produced in ionization region and around y = 0.5 lm and production decreases nearer the dielectric side-walls as well as dielectric anode and cathode.
Increase in applied electric field, micro-channel length and gas pressure, and decrease in channel width cause an increase in the total number of UV photons reaching the dielectric anode. This is due to increased amount of excitation of neutrals inside the narrow channel. This leads to greater cumulative degradation of surrounding dielectrics over time.
Dissociation of C-C bonds can reduce progressively the mass of the fragments of epoxy that eventually may deteriorate the insulation. Moreover, cross-linking between different chains holding free radicals can occur, along with oxidation and graphitization or other chemical reactions. These processes are triggered by a production of such radicals and double bonds on the surrounding dielectrics of narrow tree channel under the action of energetic charged Fig. 10 a Surface plot and b contour map of radiative states (average respect to time) inside narrow channel; applied field across narrow channel E = 20 kV/mm, narrow channel length L x = 20 lm, narrow channel width L y = 1 lm, gas pressure inside narrow channel P g = 760 Torr Fig. 11 a Surface plot and b contour map of radiative states (average respect to time) inside narrow channel; applied field across narrow channel E = 30 kV/mm, narrow channel length L x = 20 lm, narrow channel width L y = 1 lm, gas pressure inside narrow channel P g = 760 Torr particles and UV photons produced during PD process. Unfortunately, experimental data for C-C dissociation cross-sections and the other possible chemical reactions are not readily available in literature. Thus, the studies undertaken in this work provide a comparative analysis of the general manner in which various parameters affect degradation due to photon impact, rather than an exact quantitative analysis.

Conclusion
The RT-PIC model presented in this paper successfully reproduces the salient features of light emission during PD pulse within a narrow channel, and yields detailed information about the extent of C-C bond breakage. This gives us information about the extent of degradation of the surrounding dielectrics due to UV photon impacts produced by PD pulses. C-C bond breaks due to impact of electrons is negligible and is mainly related to photodegradation.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

Appendix 1: Electrostatic equations for dielectric narrow channel
The total charge density r T on the anode metallic electrode may be expressed by [21]: where r t 0j is the surface charge density on anode electrode at time 't'. The simplified form of equation (26) of [21] for a 'simple metallic-metallic' configuration (ref. [20]) is: where, u t 0 ,u t 1j , N are the potential on the anode at 't', potential at x 1 (=1 9 Dx) and number of divisions along yaxis. Particle 'p' contributes to Q ij through the space charge density q ij (bilinear weighting: [21]).
The equivalent circuit for the narrow channel bounded by dielectrics is shown in Fig. 2b. C g , C d and C r represent the equivalent capacitances of the gap, the dielectric layers and dielectric material parallel to the micro-channel respectively, and J conv is the discharge current density. The time variation of r T , the total charge density on the driven Fig. 12 a Surface plot and b contour map of radiative states (average respect to time) inside narrow channel; applied field across narrow channel E = 40 kV/mm, narrow channel length L x = 20 lm, narrow channel width L y = 1 lm, gas pressure inside narrow channel P g = 760 Torr Fig. 13 a Surface plot and b contour map of radiative states (average respect to time) inside narrow channel; applied field across channel E = 50 kV/mm, channel length L x = 20 lm, channel width L y = 1lm, gas pressure inside channel P g = 760 Torr dielectric electrode, may be obtained from the Kirchhoff's Current Law: where I ext , r T and I r the total external current, charge on the dielectric surfaces and current in the dielectric material parallel to the narrow channel. The time variation of voltage drop (U d ) across the dielectric slab is obtained as: Figure 2b shows a simple external circuit consisting of a voltage source in series with a capacitor coupled to the left electrode. The voltage drop V c across the capacitor in the external circuit is obtained from voltage balance as, V(t) is the applied voltage source. The discrete finite differenced form of (23) can be expressed as: where R Idt ¼ Q c ¼ CV c and R I r dt ¼ Q r ¼ CU r and are the charge on the external capacitor and bulk of dielectric material parallel to the narrow channel, and Q t conv is the charge deposited on the electrode from the discharge in the cavity during the time interval (t -1, t). Combining Eqs. (24) and (25), the following relation can be obtained: Similarly, the discrete finite differenced form of (23) can be expressed as: If Eqs. (24) and (27) are combined, then we have: By inserting the value of U t d to the Eq. (26) and taking a = (1 ? 2(C ? C r )/C d ) -1 , the value of r t T can be written as follow: Equations (20), (21) for r t T and (29) can be combined and solved for u t 0 on dielectric anode to produce: Appendix 2: Derivation of radiation transport matrix The probability of a photon emitted at ŕ within the (j,k)th cell to be reabsorbed at the center of the (j,k)th cell, r jk , is defined as where R lm dr 0 is a volume integral over the (l,m)th cell, and the kernel function G(r,ŕ) is calculated from the transmission factor as With uniform ground-state density profile, the line shape, g(t, r), and the absorption coefficient, k(t, r), depend only on frequency of the photons, t, and thus the transmission factor is where R = |r -ŕ|. Pressure broadening which has the Lorentz line shape is the dominant broadening mechanism in a PDP cell, and for this case, Eq. (33) becomes [34]: where I 0 (x) is the zeroth-order modified Bessel function. The equation can be approximated at high opacity (k 0 R) ) 1 as From Eq. (32) with divergence theorem The surface integral in Eq. (36) has four segments in Cartesian coordinates as shown in Fig. 26 With the property of shift invariance and isotropy, it has been verified that A jk;lm is only a function of the distance between the positions of photon emission and re-absorption, and thus A jk;lm = A |j-l|,|k-m| . For the top boundary in Fig. 26: where R ? ¼ R cos h ¼ ½ðx j À x 0 Þ 2 þ ðy k À y m À Dy=2Þ 2 1=2 and z ¼ R ? tan h. With the high opacity approximation, Eq. (35) A ðaÞ jÀl j j; kÀm j j ¼ À which is a single integral over position x. Here, C(1/ 4) = 3.6256 is a specific value of the gamma function and comes from the relation For the bottom boundary In a similar way, for the right-hand side boundary ð kÀm j jÀ1=2ÞDy Appendix 3: Derivation of spectral radiation flux At the back plate, spectral photon flux is calculated as; where r = (x i , y 0 , 0) is a location on the phosphor surface, V t is the volume of the plasma region, and S r (ŕ) is the photon emission rate defined as: With the piecewise constant approximation: If we consider three-dimensional distribution of the excited state, the infinite range of z directional integral should be changed to a finite range. However, we assume infinitely long z length with only two-dimensional distribution. With R ? ¼ ½ðx i À x 0 Þ 2 þ y 02 1=2 , R ¼ R ? sec h and z 0 ¼ R ? tan h, With the relation that Fðx i ; mÞ ¼ gðmÞ 4p where Ki 2 (x) is the second repeated integral of the zerothorder modified Bessel function, K 0 (x) [65]. However, Eq. (49) needs huge computations, and thus we calculate it using average distance between the phosphor surface and Fig. 19 a Surface plot and b contour map of radiative states (average respect to time) inside narrow channel; applied field across narrow channel E = 40 kV/mm, narrow channel length L x = 20 lm, narrow channel width L y = 2 lm, gas pressure inside narrow channel P g = 760 Torr Fig. 20 a Surface plot and b contour map of radiative states (average respect to time) inside narrow channel; applied field across narrow channel E = 40 kV/mm, narrow channel length L x = 20 lm, narrow channel width L y = 10 lm, gas pressure inside narrow channel P g = 760 Torr Fig. 21 a Surface plot and b contour map of radiative states (average respect to time) inside narrow channel; applied field across narrow channel E = 40 kV/mm, narrow channel length L x = 20 lm, narrow channel width L y = 80 lm, gas pressure inside narrow channel P g = 760 Torr the center of the cell where a photon is emitted, which is defined as R 0 . If R 0 is enough larger than Dx and Dy, we can pull R Án out of the volume integral in Eq. (44) Appendix 4: Derivation of total radiation flux The total photon flux is: From the integral of Eq. (46) over frequency t: Fðr; mÞ ¼ 1 4p S r ðr 0 ÞTð r À r 0 j jÞ With a high opacity approximation, T(R) = (pk 0 R) -1/2 :   Simulation mesh diagram of radiation transport for a photon emitted at (l, m)th cell to be reabsorbed at the (j, k)th cell. With the piecewise-constant approximation, the density at each cell is assumed to be constant [39] Fðx i Þ¼