A parametric study on the PD pulses activity within micro-cavities

A two-dimensional kinetic model has been used to parametric investigation of the spark-type partial discharge pulses inside the micro-cavities. The model is based on particle-in-cell methods with Monte Carlo Collision techniques for modeling of collisions. Secondary processes like photo-emission and cathode-emission are considered. The micro-cavity may be sandwiched between two metallic conductors or two dielectrics. The discharge within the micro-cavity is studied in conjunction with the external circuit. The model is used to successfully simulate the evolution of the discharge and yield useful information about the build-up of space charge within the micro-cavity and the consequent modification of the applied electric field. The phase-space scatter plots for electrons, positive, and negative ions are obtained in order to understand the manner in which discharge progresses over time. The rise-time and the magnitude of the discharge current pulse are obtained and are seen to be affected by micro-cavity dimensions, gas pressure within the micro-cavity, and the permittivity of surrounding dielectrics. The results have been compared with existing experimental, theoretical, and computational results, wherever possible. An attempt has been made to understand the nature of the variations in terms of the physical processes involved.


Introduction
Some of the principal factors that can influence partial discharges (PD's) and other such discharges inside micro-cavities are the dimension of micro-cavity, the pressure of the gaseous environment and gas composition, the magnitude and frequency of the applied voltage, the electrode arrangement and geometry, and the properties, condition, and age of the insulation. The different character of the discharges affects in turn the response of conventional PD pulse detectors, which is dependent to the pulse shape and its rise-time [1][2][3].
The rise-time of PD pulses is one of the most important subjects. The rise-time of the discharge pulses determines the type of the PD detection circuit that must be employed to the detection of the discharge pulses and record their pulse shape attributes with specified precision and accuracy [4]. While it is evident that the rise-time of the detected PD pulses is depend upon the bandwidth of the discharge detection equipment as well as location of PD source, the rise-time is also very much dependent to the discharge mechanism as such and the availability of free electrons to initiate that [4]. In contrast to the long rise-times characterizing the discharge pulse by glow discharges, spark-type discharges exhibit typically short rise-times and large magnitude pulses [2,4]. Bartnikas et al. [2] studied the effect of over-voltage on PD rise-time. They reported that with greater over-voltages, the PD pulses assume shorter rise-times and pulse widths as well as increasing magnitudes, thereby inferring larger charge transfers.
It has been demonstrated earlier that the discharge mechanisms in a micro-cavity are essentially of the Townsend type and the form of the ensuing discharge current pulse at the time of a spark-type breakdown is controlled by the electric field created by the ion space charge in the micro-cavity [1]. Ganjovi alirezaganjovi@yahoo.com; ganjovi@kgut.ac.ir A few studies deal with cavities in the submillimeter range; these include modeling of avalanches inside cylindrical micro-cavities by Burgener et al. [5] and experimental studies by Gu et al. [6]. As already reported in the literature, in large disk-shaped voids, transition takes place from streamer mechanism to Townsend mechanism due to a conducting layer formed on the dielectric walls by discharge by-products [7][8][9][10][11]. However, an important point to note is that streamer discharges do not occur in submillimeter cavities, in the absence of the critical length required for formation of a streamer head [1].
Kinetic PIC-MCC simulation methods are known to be better suited for problems in which the particle distribution may deviate significantly from Maxwellian, such as when trapping or particle stochastic heating occurs [12][13][14]. Additionally, the PIC-MCC method guarantees high accuracy for both transient and stationary problems [15], and yields more detailed information about the evolution of the discharge. Radjenovic et al. [16,17] used a one-dimensional PIC-MCC model to study Paschen's law for gas breakdown at small electrode separations.
A two-dimensional kinetic model based on PIC-MCC technique to study the behavior of spark-type PDs inside the micro-cavities has been developed [18]. The same model was used to investigate the filamentary distribution of spark-type PD's inside a micro-cavity [19]. In another work, a radiation transport (RT) model integrated with the PIC-MCC model was developed to study the behavior of excited (radiative) states of air within discharge in the narrow channel [20]. The radiative state atoms were described by a fluid model combined with the Holstein-Biberman equation.
Using this model, in this paper, the influences of different parameters of micro-cavity on the different features of a spark-type PD pulse such as discharge current, current rise-time, etc., with two different electrode configurations, i.e., metallic-metallic (MM) and dielectric-dielectric (DD) are examined [18]. The model has been employed for performing an extensive parametric study relevant to effect of quantities such that micro-cavity dimensions, gas pressure inside the micro-cavity, and dielectric permittivity on the PD development within MM and DD micro-cavity. The intent is to discuss the results obtained systemically using the model to establish the effect of above quantities on the nature and shape of PD pulses in micro-cavities undergoing a spark-type discharge in air.

Computational techniques
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 [21]. Next, particle boundary conditions are applied. For modeling collisions, the Monte Carlo collision (MCC) scheme is applied [22]. 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 [13]. 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 Fig. 1 Schematic representation of the computation sequence for particle in cell-Monte Carlo (PIC-MCC) method 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 using a Monte Carlo Collision technique [12]. Three main charged particle species, viz. electrons, positive ions, and negative ions, have been considered. The major electron-neutral molecule impact reactions included are (1) elastic, (2) excitation, (3) ionization (including all important ionization reactions), and (4) electron attachment. Related cross-sectional values are extracted from experimental data available in the literature [23]. 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 [24]. 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 ) has been assumed based on a hard-sphere collision model [25].
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Þ: The potential (U) can be separated as follows [12,26,27]: where U P and U L represent the Poisson and Laplacian parts of electrical potential (U). Therefore, the field may be described by the combination of the following equations: r Á e 0 rU Li ¼ 0: The boundary condition for Eq. (3) is U = 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 U i = 0 on the equi-potential surface, and U = 0 elsewhere which gives us U Li (the potential due to the ith Laplacian field). U L is obtained by the superposition of all the Laplacian fields. Neumann boundary conditions are included through U 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 [12,21,28].
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.

Metallic-metallic configuration
The first case considers a cavity between MM electrodes. For this configuration, the potential at the anode can be obtained as follows [29]: where 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', electrode cross-sectional area, capacitance of the blocking capacitor, driving voltage, charge in the external circuit, charge convection due to discharge, space charge density within the discharge, the electric potential at each grid point, and the number of grid divisions along the y-axis.

Dielectric-dielectric configuration
The second case we implement is that of DD where the micro-cavity is bounded in both sides by dielectrics. This necessitates the incorporation of a second dielectric slab between the metallic cathode and the void wall at the cathode end (Fig. 2a). The equivalent circuit is shown in Fig. 2b, where C g , C d1 , and C d2 represent the equivalent capacitances of the micro-cavity and the dielectric layers, respectively. Assuming C d1 = C d2 = C d , using KCL, it is possible to obtain the following relations: where U d , I ext , r s , J Conv , and V c are, respectively, the potential across each dielectric slab, total external current, charge density on the dielectric surface, discharge current density (I conv = I d = AJ Conv ), and voltage across blocking capacitor. By combining the above equations, Eq. (5) may be rewritten for DD case as follows

Simulation results
This section is devoted to the presentation of the results obtained by numerical modeling. A grid-size of Dx = Dy = 2 9 10 -6 m, and a time-step of Dt = 10 -14 s are used; typical simulation times of 24 h are required with a dual processor SUN workstation. As a first step, in this section, we trace the evolution of the distribution of charged particles within the submillimeter void. The scatter plots in phase space for velocities in x direction i.e., (V x , x, y) for electrons and positive ions in Figs. 3 and 4 at particular instants of time. As shown here, drift velocity of electrons formed during the discharge is much higher than ions (10 3 times) and quickly traverse the micro-cavity. All charged particles are distributed about the center-line y = 50 lm. Most of the charged particles are created very close to the anode. In the early stages, positive ions are distributed evenly over the inter-electrode distance, except for a small region near the cathode. It is also uniformly distributed along y-direction. The number of ions near the anode increases by ionization with passage of time. The positive ions slowly drift toward the cathode, till at about 13 ns, the entire gap distance is almost evenly populated. From this time onward, ions gradually leave the micro-cavity, thus depleting the ions resident in the gap. This is similar to the evolution of the electrons within the gap, as seen in Fig. 3. Electrons move faster and travel only the distance between the ionization region and the anode. As is clear from Fig. 4, the ion velocities increase as they reach near the cathode, accelerated by the electric field.
In following, the effects of micro-cavity dimensions, gas pressure, and dielectric permittivity on the PD pulse features such as discharge current magnitude and its rise-time have been investigated. Obviously, not all PD pulses are characterized by very a rapid rise-time; with a high probability of occurrence of a free electron within the gas of occluded cavities, only small or moderate over-voltages may be required to precipitate a breakdown event. Unless otherwise mentioned, all simulations reported in this work are performed for operating conditions specified in Table 1.

Micro-cavity length
For a relatively short anode-cathode gap distance, the effect of the micro-cavity dimensions on the discharge characteristics, from experimental as well as theoretical findings, is evident [7,[30][31][32][33][34][35][36][37][38][39][40][41]. In this paper, the length of the micro-cavity is varied from 50 to 600 lm, while the applied electric field across micro-cavity is held constant at E = 10 6 V/m. As the length of the micro-cavity increases, for a given pressure, the residence time of ionizing particles inside the micro-cavity is higher. Each electron suffers a larger number of ionizing collisions as it travels from the cathode to the anode; as a result, the production of charged particles inside the micro-cavity is much higher at higher lengths. This is evident as peak value of charged particles    [41]. Aging and life models by Montanari and coworkers [34,35] have similarly shown that the number of electrons generated inside the cavity increases with cavity length. The electrons generated inside the cavity are the main carriers of PD current [40,41]. Increase in the number of charged particles is reflected in an increase in the magnitude of discharge current shown in Fig. 5d. This finding is in agreement with experimental results reported by Devins [33], which shows that the discharge current is sensitive to length of micro-cavity, increases with microcavity length. Bartnikas [32] varied gap length from 0.1 mm to 1 mm and reports that a reduction in the gap length results in the disappearance of the pulse discharge (spark-type) and the onset of a diffused glow (pseudo-glow and glow discharge). Figure 5e shows the variation of PD pulse current risetime as a function of micro-cavity length. The rise-time increases as length increases. This is understandable as the ionization process in longer micro-cavities require a longer time to reach peak values. These finding corroborates the experimental results by Morshuis [36]. He studied the pulse shape of PD in disk-shaped cavities in polyethylene and other polymers. He varied the length of the cavities between 30 lm and 3.7 mm. He reports that for cavities with a diameter/length ratio )1, PD rise-time is lower [7,40]. Ramachandra [37].
The variation of charge accumulation at the dielectric anode and cathode at the end of discharge pulse is shown as a function of micro-cavity length in Fig. 5f. The amount of accumulated charge density is calculated as the difference between the charges at the beginning and end of discharge, (Dr s = r si -r sf ). As larger numbers of charged particles are deployed in the discharge with increased length, the amount of charge accumulated on the dielectric surfaces also increases. This value is always negative for the anode, as the absorbed negative charges lower the applied potential on the anode. For the MM case, the accumulated charge is soon dissipated into the external circuit, and there is no charge accumulation at discharge extinction. Figure 6 depicts the number of charged particles for both MM and DD configurations and for two different microcavity lengths. This is in agreement with simulation results obtained by Novak et al. [43] using a fluid model, for a 0.5mm gap and DD configuration. Negative ions are depleted fast. As is seen here, the number of negative ions is much lower than that of electrons or positive ions. For the DD configuration, a smaller number of charged particles are involved. Interestingly, for the smaller micro-cavity length, the rise-time for the DD is smaller, and so is the fall-time.
The time evolution of the anode potential for different micro-cavity lengths, u 0 , is presented in Fig. 7a, b. A similar decrease in anode potential with time was observed by Bartnikas et al. for 0.5 mm gap and DD configuration. It may be noted that the effect of external circuit is ignored by Bartnikas et al. [1]. The values remain constant with time for MM configuration as expected, because the Ion-induced emission coefficient k i 1. conductors are maintained at a constant applied potential by the external circuit. However, for DD configuration, due to the role of dielectric on the electric field next to dielectric anode in the discharge medium and also initially presence of positive-charged particles in this region, the dielectric anode electrode potential first increases. Then, because of accumulation of negative-charged particles (i.e., electrons and negative ions) on anode dielectric, the anode potential decreases, settling at considerably decreased values at the end of the discharge-period. The anode dielectric potential, for DD configuration, at the end of the discharge (u 0f ) for micro-cavity length of 200 lm is reduced to 99.9 % of its initial value (u 0i ) and for length of 400 lm, it has reached to 97 % of its initial value. Thus, while the applied electric field across micro-cavity is kept constant at different lengths, the drop (Du 0 = u 0i -u 0f ) at higher micro-cavity lengths is higher. The change in potential is brought about by charge accumulation at the dielectric electrodes, and as it is clear in Fig. 7c, d, the charge accumulated on the dielectrics for 200-lm micro-cavity length (10 -4 mC/m 2 ) is much \400 lm micro-cavity length (0.01 mC/m 2 ). Nikonov presented the variation of the field at the anode with time [44], and this shows a trend similar to change in r. It may be noted that r = e 0 E anode , where E anode is the field at the anode, and e 0 the permittivity of free space. The variation of E anode in Nikonov's work may therefore be taken to be a corroboration of the general trend in variation of surface charge density accumulated on the anode surface.
It is interesting to study the effect of length on the PD pulse by looking at the evolution of the electric field inside the micro-cavity. Figure 8 shows the variation in the axial field (component in the direction of the field) along the micro-cavity length, at different moments of time. As may be seen, the field within the micro-cavity is originally uniform for both lengths i.e., L X = 200 lm and L X = 400 lm. Due to less strength of ionization at lower micro-cavity lengths, the effect of space charge on the applied electric field in case of L X = 200 lm is much less than L X = 400 lm (Fig. 8a, b). Wu et al. [45] showed that as length of micro-cavity increases, the effect of the field due to charged particles inside the micro-cavity becomes more evident. Generally, for the MM case, the field at the dielectric anode is reduced considerably within a very short time. Positive charges take longer to travel the gap distance, and a positive space charge region forms ahead of the anode where they accumulate. This results in a negative field ahead of the anode. The field in the rest of the microcavity is low and positive, being directed from the positive space charge region toward the cathode. As positive ions gradually travel toward the cathode, the field in the near vicinity of the anode is relieved and gradually reaches initial values. Simultaneously the field near the cathode is enhanced in the direction of the applied electric field, due  Fig. 8c, d, the evolution of the field values is similar to the previous case till about 6 ns. However, the field at the anode never recovers the original values owing to accumulation of negative charges at the anode, lowering the potential thereof. Here too, the field within the micro-cavity becomes uniform; the ultimate field value is, however, lower than in the MM case. Again due to higher effect of space charges at higher micro-cavity lengths, change in the applied electric field is much higher at higher micro-cavity lengths. This is because of the space charge field (produced due to negative charge accumulation at the dielectric anode) opposes the applied field. Figure 9 presents the electric field (E x ) as a function of time at different axial distances from the anode (along the length of the micro-cavity) for two lengths of micro-cavity i.e., L X = 200 lm and L X = 400 lm. In the MM case, the field variation is similar in both cases; the magnitudes though are different, the field level near the anode falls initially, followed by a rise which is due to lowering of field values owing to an initial positive charge accumulation in front of the anode electrode. While the electric field near the cathode is almost close to its initial value, the electric field near the anode is at a minimum (at t = 6 ns) when most of the negative charges have left the microcavity. After this moment, the breakdown mechanism becomes much more space charge controlled, and the electric field near the cathode will be under influence of the positive charges arriving in its vicinity and consequently increases. At distances near the cathode, there is a considerable enhancement in the electric field between the positive ions and cathode. Similar observations were made by Serdyuk et al. in modeling of electric discharges in a 3-mm air-gap between parallel plate solid dielectric barriers [46]. They showed that the interaction of the avalanche with the solid dielectric led to charging of the barriers. Optical studies by Morshuis et al. revealed a large number of separate discharging sites dependent on the void length [7,40]. They have reported that the current pulses of these separate discharges were integrated by the relatively slow detector, resulting in extremely high discharge values. However, in this work, it has been assumed that only one discharging site on dielectric anode (DD configuration).
Thus, it could be concluded that greater PD pulse current rise-time and charge accumulation on anode and cathode dielectric surfaces occurs at higher micro-cavity lengths.

Micro-cavity width
It has been shown that on increasing the micro-cavity width, there is a greater likelihood of movement of charge carriers in the direction perpendicular to applied electric field (y-direction) [18,19]. This results in increased dissipation of energy of the PD pulse. Ionization and production of new charged particles therefore, reduces, as was seen in the peak values of charged species. Figure 10 shows the PD pulses for MM and DD configurations for two different micro-cavity widths. For a width of L y = 300 lm, the PD pulse current peak occurs at 350 and 300 lA for MM and DD configurations, respectively. For a width of L y = 600 lm, the values for PD pulse current occur at 273 and 256 lA for MM and DD configurations which are lower than the corresponding values for the case of L y = 300 lm. Thus in the wider micro-cavity, reduced peak values are observed. The current is lower for a DD configuration. The variation in micro-cavity width does not affect the resultant PD pulse current shape. Gu et al. varied the micro-cavity diameter in the range of 50-600 lm [6] and found that the ionization coefficient increases rapidly as micro-cavity diameter decreases. This is in agreement with the presented results in this paper. Figure 11a, b shows the time evolution of the anode potential, u 0 , for different micro-cavity widths. The anode dielectric potential, for DD configuration, at the end of the discharge for micro-cavity width of 300 lm is reduced to 92.2 % of its initial value (u 0i ) and for width of 600 lm, it has reached to 95.2 % of it is initial value. Thus, a lower potential drop (Du 0 = u 0i -u 0f ) occurs at greater widths, and the probability of re-strike is therefore higher. Consequently, in Fig. 11c, d, the charge 0.0694 mC/m 2 accumulated in the 300-lm wide micro-cavity is higher than that in the 600-lm wide micro-cavity (0.0432 mC/m 2 ). These findings are in agreement with experimental results reported by Auckland et al. [47][48][49]. He studied the influence of diameter on the quantity of charge accumulating on the side-walls of a micro-cavity in specimens with (c) (d) Fig. 11 Time history of the electric potential on the anode dielectric for a width of L y = 300 lm, b width of L y = 600 lm, and the charge density on the anode dielectric for c width of L y = 300 lm, and d width of L y = 600 lm diameters in the range 50-126 lm and length of 10 mm. While the model of discharge propagation in their experiments is different (streamer as opposed to pulses), he observed a reduction in the amount of accumulated charges with increasing in diameter of micro-cavity. These results corroborate the experimental results obtained by Auckland et al. [48], indicating that the electrical strength of gas confined in narrow micro-cavities increases with the length of the micro-cavity and decreases with micro-cavity diameter.
Morshuis et al. [7] reports that tests on voids with small diameters (l-4 mm) show discharge mechanisms similar to larger widths. These findings also show that the effect of width on discharge is negligible.

Gas pressure inside the micro-cavity
The effect of gas pressure on the nature of PD current pulse is studied in this section. The gas pressure varies from 76 to 1260 Torr, while other parameters are kept constant ( Table 1).
The discharge process occurring in cavities, depending upon gas pressure, may assume either a spark, glow or pseudo-glow character. As Morshuis [36] pointed out, in small cavities, the gas plays an important role in the ionization process. The effective ionization coefficient is a property of the gas and is therefore affected by the changing gas chemistry and by a change of pressure.
It is seen in Fig. 12a-c that the maximum number of charged particles residing inside the micro-cavity increases with gas pressure. Therefore in Fig. 12d, the magnitude of PD pulse increases. The rise-time of the PD pulse versus gas pressure is presented in Fig. 12e. The effect of gas pressure on the rise-time is significantly high. At lower gas pressures, e.g., 76 Torr, the resultant PD pulse current amplitudes are too low but the rise-time is too rapid and it might be possible to detect them by conventional PD detectors. Gaxiola et al. has reported that with increasing pressure, the time-to-breakdown rapidly increases over two orders of magnitude [30].
However, the experimental findings by Kasten et al. for a cavity with length of 20 mm show that the time to peak of discharge current decreases as pressure increases [50]. His findings relate to cavities in the millimeter range, where streamers may be taken to control discharge propagation and transfer the results correspond to a different physical mechanism.
The variation of charge accumulated on the anode and cathode dielectrics (Dr s = r si -r sf ) versus gas pressure inside the micro-cavity is depicted in Fig. 12f. Owing to higher ionization at higher pressures, larger number of charges (mostly electrons) will reach the dielectric at higher pressures.
As is clear in Fig. 13, for a gas pressure of 260 Torr, the maximum values of PD pulse current is around 416 lA for MM configuration and 256 lA for DD configuration, respectively, while for a gas pressure of 860 Torr, the corresponding values are 520 and 288 lA, respectively. It may be noted that lower pressure values lead not only to reduced peak values of the PD pulse current, but also cause shorter delays in the formation as well as reduction in the width of PD pulse.
In summary, at higher gas pressures, greater production of charged particles within micro-cavity over time and consequently increase in discharge current as well as charge accumulated on the anode and cathode dielectrics occurs.

Dielectric permittivity
In homogenous insulating materials, the electric field is distributed uniformly within the entire volume. In nonhomogenous insulating materials, due to manufacturing defects, aging effects, etc., the distribution of the electric field is no longer uniform.
The electric field distributes across two layers of dielectric connected in series in the ratio of their permittivities and thickness (Fig. 2a). In this section, the effects of dielectric permittivity on the behavior of PD pulse are studied. The applied electric potential between the metallic electrodes has been kept constant (e.g., E a = 20 kV/mm) for the DD configuration, i.e., the electric potential across the dielectric-gap combination. The permittivity of surrounding material is varied, while keeping the dimensions and geometry constant. The applied electric field inside the micro-cavity (E a ) increases and ionization inside the microcavity becomes stronger. Owing to this, as is clear from Fig. 14a-c, the maximum number of charged particles produced inside the micro-cavity increases. The magnitude of the PD pulse as a function of dielectric permittivity is presented in Fig. 14d. Figure 14e shows that the rise-time of the PD pulse falls with e r , probably due to increased space charge effects at the dielectric cathode, and decreases with dielectric permittivity as effective electric field along micro-cavity increases.
The variation of charge accumulated on the dielectric anode and cathode (Dr s = r si -r sf ) with dielectric permittivity is depicted in Fig. 14f. Owing to increase in the strength of ionization with permittivity of dielectric, the amount of charge accumulated on the dielectrics surface increases.
A point to be noted is that the effect of dielectric permittivity on charge accumulation on the dielectric anode and cathode, PD pulse magnitude, and its rise-time is significant. All simulations reported in the previous sections have been performed with a dielectric permittivity of e r = 5.

Discussions
An important contribution of this work is its ability to investigate of the effects of different parameters on PDs inside the micro-cavities for two different dielectric configurations (MM and DD). The calculated PD pulse forms are very much affected by values of micro-cavity length, gas pressure, and dielectric permittivity, while the effect of variations in micro-cavity width is relatively less. These results show that for a kept fix each one of the above parameters, the detected PD pulse form, amplitude and rise-time can exhibit considerable variation due to possible changes in each every above-mentioned parameter. Since the response of a conventional PD detector is a function of the rise-time and pulse amplitude of the PD pulse, variations in the detected pulse would result.
The accuracy of peak pulse measurement would thus be affected in that PD pulse detector would respond more favorably to pulses having rapid rise-times. Consequently, it is quite conceivable that a lower amplitude discharge pulse with a more rapid rise-time may produce a larger indication on that detector than a higher amplitude pulse with a slower rise-time. It would be difficult to resolve this measurement inaccuracy in practice.
By increasing some parameters such as micro-cavity length, gas pressure the space charge density accumulated at the cathode increases and leads to less rapid breakdowns while by increase micro-cavity width and dielectric permittivity charge accumulation on the cathode leads to more rapid breakdowns and results in pulses with correspondingly reduced rise-time and augmented magnitudes.

Conclusions
A two-dimensional kinetic model based on PIC-MCC has been developed to a parametric analysis of nature of PD pulse within micro-cavities. Indeed, this parametric analysis focused especially on micro-cavity dimensions, gas pressure inside the micro-cavity, and dielectric permittivity. It proves the ability of the model in reproducing what is known about the dependence of the nature of PD pulse on the various physical parameters. All the parameters that characterize the proposed model have physical meaning and they can be measured resorting analytical to analytical techniques. So, determining the dependency of the nature of highly localized intense spark-type PDs to different parameters governing them gives an idea of the physical processes involved in their operations within the micro-cavities. It is expected that the employed model will be a useful tool in designing the solid insulating materials.
In summary, while the peak magnitude of PD pulse current increases with micro-cavity length, gas pressure, and dielectric permittivity, it decreases with increase in width. The rise-time of the PD pulse increases with microcavity length and pressure, but it decreases with microcavity width and dielectric permittivity.
Note that this is Eq. (5) in the paper and is used for the MM case.
The equivalent circuit for DD configuration is shown in Fig. 3b. C g and C d represent the equivalent capacitances of the gap and the dielectric layers, respectively, and J Conv is the discharge current density. The time variation of r T , the total charge density on the driven dielectric electrode, may be obtained from the Kirchhoff's Current Law: where I ext and r T are the total external current and charge on the dielectric surfaces. The time variation of voltage drop (U d ) across the dielectric slab is obtained as: Figure 3b 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 Eq. (A4) can be expressed as: where R Idt ¼ Q c ¼ CV c is the charge on the capacitor 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. (A6) and (A7), we obtain: Similarly, the discrete finite differenced form of Eq. (A5) can be expressed as: If we combine Eqs. (A6) and (A9), then we have: By inserting the value of U t d to the Eq. (A8), the value of r t T can be written as follow: Equations (A1), (A2) for r t T and (A11) can be combined and solved for u t 0 on dielectric anode to produce: