FEA-Based Simulation of Accelerated Ageing in a Power Cable Due to Sustained Partial Discharge Activities in a Spherical Cavity

The reliability of electrical assets is greatly influenced by the quality of their insulations. Key power installations such as power cables are manufactured with polymer-based materials as part of their insulation system. However, accelerated ageing of equipment insulations due to manifestation of defect(s), and partial discharges (PDs) can offset the operation of these systems or even lead to breakdowns. In this study, a non-deterministic model to simulate the phenomenon of repetitive discharges in a spherical air-filled cavity within a practical power cable has been investigated. In addition, the work contributes to the understanding of PD behaviour and field distribution under different ageing conditions considering changes in cavity surface conductivity. First, a section of the practical XLPE cable containing the cavity is developed in 2D using COMSOL software, and a finite element analysis (FEA) of the electric field distribution within the cable insulation is performed. The magnitude of the cavity local field, that is enough to ignite a PD, is investigated. Alongside the COMSOL model, the activity of sustained internal PD is simulated in MATLAB by introducing a random sample generating factor and adjusting the model’s parameters to obtain something close to the practical results. Furthermore, the impact of continuous PD in the power cable under different cavity dimensions and surface conductivity is likewise investigated, and a phase resolved PD (PRPD) pattern is established. The result shows that the magnitude and number of PDs per cycle increase as the cavity size and its surface conductivity increase. Finally, when the cavity surface conductivity rises, the amplitude of the electric field generated by the surface charge distribution and the number of PDs per cycle approach their maximum values.


Introduction
Cross-linked polyethylene (XLPE) cables are widely acceptable variants of polymer-based power cables used in electric power grids.These cables are largely used because of their good mechanical, electrical and thermomechanical properties.In addition, the ease of installation and maintenance accord the XLPE cables added advantage, particularly for underground power transmission and distribution purposes [1,2].According to earlier research, high-voltage (HV) and medium-voltage (MV) power cables account for over 11% of the entire power distribution network.For such reasons, PD studies, as it relates to power cables, has become imperative for effective and reliable operation of the overall power system [3].
Partial discharges (PDs) have been explained by researchers from varying perspectives.According to IEC 60270, PDs are stochastic events that describe a localized breakdown in insulations between electrodes [4][5][6].The knowledge of PD is critical for assessing the level of degradation in insulation systems of HV plants as its emergence signifies ageing [7].In electric power cables, PD often occurs in the insulation-defect(s) by the availability of a seed charge and sufficiently high field [8][9][10].Sustained discharge activities are likely to be maintained in the defect(s) with the availability of residual charges either from past PD events or defect(s) in the wall.
The geometrical and physical properties of the discharging domain significantly influence PD behaviour [11][12][13].The greatest physical property that impacts PD is the local electric stress as its magnitude affects the ignition and extinction of the PD phenomenon [14].In a cavity-defect, the size and entrapped gaseous composition, mostly from manufacturing imperfections, influence the discharge parameters.For example, the rate of discharge among other PD associated parameters is highly dependent upon the cavity size [14,15].When in service, power cables are usually subjected to high and continuous electrical excitation [16,17].This excitation, in turn, ensures constant external fields in the insulation bulk, capable of sustaining an ongoing discharge phenomenon within the cavity.The result of such a process is a rapid erosion and degradation of the insulation system, finally leading to the cable's complete failure [17,18].Replacing a failed cable in power grid can last up to a few hours or days, and during this period the affected areas are subjected to blackout.Preventive actions such as condition monitoring of the facility will go a long way in improving the overall system resiliency.A critical diagnostic tool used for assessing HV installations is PD [19][20][21][22].PDs are both a cause and a consequence of dielectric ageing, activating one of the quicker degradation processes in solid polymeric dielectrics like XLPE [23][24][25].The changes on the behaviour of PDs characteristics, such as the φ-q-n patterns, q max , q mean , q min , among others during ageing, can be correlated with the variations on the cavity surface [26,27].The pressure and gas composition, and the inception of trees [28], can be used for implementing diagnosis and prognosis analyses.The φ-q-n patterns are also known as PRPD patterns, which corresponds to a graphical representation, in 2d or 3d, of the distribution in phase of the PD pulses, magnitudes and rate, taking as reference a period of the sinusoidal wave of the applied voltage.They present particular distributions depending on the PD source.Some authors have proposed life models and health indexes based on the characteristics of PD behaviour in simulations and measurements [23,24,29].Thus, an in-depth understanding of sustained discharge activities is critical for the analysis and interpretation of PD severity in the equipment.
Prognosis of sustained PD activities in HV apparatus, based on tests, are at an advanced stage [30][31][32][33].However, the case is different through simulation as only a few PD parameters can only be obtained by measurements.This scenario has made it difficult to investigate the physical processes of PD in a cavity from modelling perspective.To this moment, handful of models aimed at simulating PD in an insulation-cavity have been proposed in the literature [18,[34][35][36][37].Although all models, to some extent, have provided appropriate insights into the behaviour of PD, the FEA-based PD models have proven to be more accurate [38].This model utilizes a number of experimentally determined physical parameters to explain the behaviour of sustained PD in systems.Section 2.2 summarizes the current state of the art for models for simulating PDs in cavities within solid dielectrics.
In this paper, a computer simulation model of PD activity in a single-core, 11 kV, XLPE power cable containing a cavity is presented.In addition, simulation results of sustained PD activity and electric field distribution, considering different cavity sizes and changes in the surface conductivity, due to the PD induced ageing, are presented.This paper is organized as follows: after the introductory part, already presented at Sect. 1, the theoretical background of the PD modelling is presented in Sect.2; then, the results and discussion with the case study are presented in Sects.3 and 4, respectively, and lastly, some conclusions and future work are depicted in Sect. 5.

FEA PD Modelling
PDs are a phenomenon that has been modelled using several approaches, the key difference being how the electric field created by prior PDs is modelled [5].In FEA models, the electric fields are calculated by solving the partial differential equations in the domains of interest using the conductance and electrostatic models [8,9].Below, the modelling process in FEA of PDs in gaseous cavities within the solid dielectric is described in detail.

Multiphysics Model for Field Calculation
For a non-dispersive, bounded, homogenous and isotropic material, the basic field equations that describe the behaviour of electric field variables (electric scalar potential and field strength) are derived from the well-known Maxwell's equations as follows [3,11]: where → E is the electric field strength, → D is the electric field displacement, σ is electric conductivity, → J represents the free current density, ε 0 ε r represent the relative and vacuum permittivity of the insulation system, and U denotes the general location's scalar potential.

Stochastic Model for PD Calculation
For the sake of completeness, the evolution of PD simulation models, in cavities within solid insulation, is considered in the first part of this section, after which the theoretical framework of the PD model used for simulation is presented.
The analytical model is based on the electromagnetic theory and some considerations that, under practical conditions, are valid [39].A generalization of the analytical model was presented in [8] to simulate PDs in cavities and protrusions and has been utilized to obtain the theoretical magnitude of PD charge within insulated cables [40].In [41], the analytical model was used to study the characteristics of PDs under direct current (DC) voltage.Deterministic variations of analytical models for PD have also been proposed that can approximately represent the PD charge magnitudes.However, they are unable to reproduce the phase distribution of PD pulses that have been found experimentally [42,43].Although the analytical stochastic models show good results when compared to experimental measurements, they do not consider all the physical processes and phenomena involved during PD, such as the memory effect, charge decay, temperature and pressure variations during and after PD activity [2].A critical discussion of the analytical model considerations, as well as, an improved version that increases the modelling accuracy for cavities of different sizes, are presented in [10].In [21], a multiphysics version of the analytical model was presented that allows calculating the pressure and temperature variations within the cavity subjected to sustained PD activity.
Gemant and Philippoff [44] first introduced the threecapacitance model to evaluate power losses in cables owing to PD in cavities.Later, Whitehead [34] proposed a modified equivalent circuit that is known as the "abc" model.In this model, unlike what happens in the analytical and FEA models, analysis are made as a function of currents and voltages in the equivalent capacitances' terminals, allowing the simulations to be easily implemented in a dedicated circuit and network analysis software, making this model widely used [45][46][47][48][49].A variable resistance technique may also be used to represent the surface charge dynamics [50].Nevertheless, the high non-linearity introduced into the circuit by the streamer equivalent resistance causes numerical challenges, necessitating extensive computing effort to solve the resultant circuit.In [2], a modification to the equivalent "abc" circuit presented in [51] was proposed to consider the stochasticity of the PD phenomenon and a method was proposed to calculate the equivalent capacitance of the cavity that allows to reconcile the models of capacitors and electromagnetic fields.
FEA models numerically solve the field equations in partial derivatives and has practically no limitations regarding the geometry of the test object, the waveform of the applied voltage and even the linearity or isotropy of the materials.In addition, multiple-physics can be considered; electrical, thermal and mechanical, interrelated during PD [15], the charges and currents, induced and real, are calculated numerically in the electrodes using boundary conditions and field solutions and in the cavity, without the need to use analytical expressions and electrostatic approximations.These models have been used to investigate the behaviour of PD in cavities under various frequency conditions [12,52,53], amplitude and waveform of the applied voltage [41,[54][55][56], cavity size and location [57,58], temperature [15,59] and number of cavities [60].Electrostatic FEA models have been used to analyse changes in PD behaviour during the ageing process [53], and in [14] a 3D electrostatic model was used to study the PD phenomenon in a geometry similar to that found in a power cable of practical dimensions.In [5], a generalized electrostatic model for PD simulation using a surface continuous charge distribution was presented.In [61], a new hybrid model was proposed that combines the advantages of the electric current and electrostatic FEA models and allows to overcome the difficulties about the surface charge distribution of the existing electrostatic FEA models [62].
There are other advanced models in which the physical processes of impact ionization, recombination, addition, drift of charge carriers and diffusion are quantitatively defined by fluid equations to describe cavity plasma dynamics during PD [63].In these models, the chemical and physical phenomena that occur during PD in cavities are considered in a precise and detailed manner [64].However, the solution of these models requires high computational consumption which makes them impractical for multiple PD analysis and ageing where simulation for hundreds of cycles of the AC voltage signal is required.
Experimental studies have revealed the PD phenomenon to be non-deterministic.For a cavity containing gaseous substance, such as air, two conditions must be satisfied for PD to begin.First, the electric field strength at the cavity centre must be strong enough to ignite and maintain an avalanche (streamer inception); and then, there must be a seed electron within the cavity to kickstart the ionization process.While the magnitude of the local electric field strength can be computed from FEA, the condition about the availability of the seed electron is evaluated through a critical avalanche criterion and the electron generation rates from different sources [8].

The Streamer Inception and Extinction Fields
The minimum electric field strength required for initiating a streamer-type PD within a cavity, denoted as E inc , depends on the cavity diameter as well as the pressure and the ionization parameters of the gas inside the cavity.Once the radius and the pressure of the gas filling the cavity are known, the inception magnitude of the electric field is obtained using the expression [8,9]: where r cav and p represent the radius and cavity pressure, respectively, and B, n and (E/p) cr are parameters related to the gas ionization process.For air-filled cavities: n = 0.5, B = 8.6 m 1/2 .Pa 1/2 and (E/p) cr = 24.2V.Pa −1 .m−1 [9].After a PD event, the electric field strength produced by the chargers left by the streamer at the cavity surface, Poisson's type field, opposes the electric field strength established by the HV source, Laplacian Field, and when its magnitude is lower than an extinction value, E ext , the streamer propagation inside a discharging cavity ceases.The electric field extinction magnitude is estimated using the expression [9]: where γ is a dimensionless factor that changes accordingly with variation in voltage polarity.The value of γ is usually obtained from experimental measurements.The PD occurrence in a cavity entails increasing the value of the cavity gas conductivity.At a predefined value of conductivity, σ cav_max , the intensity of the field inside the cavity reduces drastically.By maintaining the conductivity at this higher value, the cavity field recedes further, and the PD signals cease to exist at certain charge magnitude.

Electron Generation
Seed electrons required to initiate an avalanche and are produced by two basic processes: volume production and surface emission.In the volume electron generation process, the seed electron is assumed to emerge from ionization activities caused by background radiation in a virgin cavity.In the case of surface emission, first electrons are produced through detrapping of electrons, previously left by the PD activity, from the cavity wall [8].Therefore, the overall electron generation rate for a specific instant of time, N e,tot (t), is defined as the sum of the charge generating processes caused by volume ionization, N e,vol (t) and surface emission, N e,surf (t): Therefore, the overall electron generation rate, N e,tot (t), is the summation of the charge generation processes due to volume ionization, N e,vol , and surface emission, N e,surf (t) [8,9].
Starting with the volume ionization process, the initial electron generation rate is computed using [9]: where ψ rad represents the radiation quantum flux density, C rad indicates the interaction between radiation and gas molecules, β 0 is an exponent in approximation for effective ionization coefficient, v is the quotient of the applied voltage, U 0 , to the inception voltage, U inc , and (ρ/p) 0 denotes the pressure reduced gas density [8].
In the surface emission process, the electron de-trapping mechanism is highly dependent on two quantities: the temperature within the cavity and the local electric field strength magnitude at the cavity centre.Niemeyer argued that the said process of charge generation obeys the Richardson-Schottky scaling and is governed by [8,9]: where T represents the absolute temperature, k b is the Boltzmann constant, v 0 is the fundamental photon frequency of insulation material, φ is the effective de-trapping work function, and N e,det is the total available charges to be de-trapped for an immediate PD event and is obtained using [9]: where τ dec indicates the effective time decay constant, e is the elementary charge, q is the real charge on the cavity surface after the previous PD event, t PD is the time passed since the last discharge, and ξ is a proportionality constant that describes the fraction of q in a de-trappable condition.
Owing to the variation in permittivity of the surrounding material and the gas trapped inside the cavity, the applied field, E 0 , is multiplied by a dimensionless field enhancement factor, f c , to step up the local field inside the cavity.The value of f c is obtained using the expression [65]: where K is a parameter that depends on the cavity geometry and ε r,surf is the relative permittivity of the dielectric material.K is given by [65]: where b and a are the geometrical axes of an ellipsoidal cavity perpendicular and parallel, respectively, to E 0 .Since the insulation-cavity considered in this work is of spherical cross-sections (i.e. a = b), the value of K used is 3.After a PD event, aside the Laplacian field, f c E 0 , a second component of the electric field within the cavity is the field, E s , and is due to charges on the cavity wall from previous PD activity.Both fields (f c E 0 and E s ) combine to form the net cavity electric field, E cav (t). Figure 1 demonstrates the movement of trapped charges within the cavity due the 123 Fig. 3 Case study: a 2D geometry of the cable section with a cavity parallel.b Meshed pattern of the PD system applied field.When the direction of E cav (which is the same as f c E 0 ) is opposed by Poison's type field (E s ), the liberated charges travel along the orientation of f c E 0 and become more compact at the bottom of the cavity surface (Fig. 1a).As the direction of E 0 changes, the directions of E cav and E s become the same as shown in Fig. 1b.At this instant, charges of opposite polarity may drift through or propagate along the surface of the cavity, resulting in charge recombination [9,65].This process reduces the surface charge density and is thought to be the principal source of charge decay in the air-filled cavity.

The PD Process
The possibility of a PD owing to the presence of a seed charge in the cavity is believed to be a function of the time step and the total electron generation rate, N e,tot (t), when no PD activity occurs.At first, the likelihood of a seed electron, given that the magnitude of the electric field strength in the cavity is higher than the inception value, P(t), at each time step, is calculated from [9]: To ascertain the availability of electron at a particular time step, a randomly generated number R (between 0 and 1) is compared with P(t).If P(t) > R, the presence of an electron is ascertained at that time step; else, the time is increased by a time step before a discharge occurs, and field magnitudes and parameters of media are calculated and updated.

Cavity Conductivity
The discharge inside the cavity is dynamically modelled by switching the cavity's state from non-conducting to conducting.The initial value of the cavity conductivity, σ cav0 , is increased to its highest value.Once the peak cavity conductivity, σ cav_max , is reached, the cavity field, E cav (t), at any instantaneous time keeps dropping, and the discharge finally quenches at a certain stage when E cav (t) < E ext [9].σ cav_max is calculated using formula (14-15): N e = q max 4 3π e r 3 cav (15) where α denotes the electron energy distribution, λ e is the electron's mean free path, m e is the mass of the electron, c e represents the thermal velocity of the electron, r cav is the cavity radius, q max represents the maximum charge magnitude, and N e corresponds to the density of the electrons.

Charge Decay
At the end of each discharge, free charges deposited on the cavity wall decay before the next PD event.The charges produced might decay through charge neutralization, ion drift or surface conduction.The surface conductivity, σ surf , is dependent upon the temperature and the field strength polarity of the cavity surface.When the direction of the applied field changes, σ surf (see Figure 1b) is enhanced from its original value, σ surf0 , to model the movement of charges due to conduction within the surface of the cavity [8,9].Assuming a constant cavity temperature, the surface charge, σ surf (t),  decays at each time step and is simulated using: Since the cavity field, E cav , is understood to be the summation of the electric field in the cavity produced by the surface charge distribution from the previous PD event and the enhanced applied field, E cav0 , and, E s , then the instantaneous value of E s is obtained from [8,9]:

Charge Magnitudes
FEA is utilized to access the current density needed for the computation of the PD charge magnitudes.Once this quantity is calculated, it is integrated over specific regions over time in order to obtain the transient current I(t) flowing in the discharging system.For each PD event, the apparent and real charge magnitudes are calculated by integrating I(t) through the ground electrode and cavity centre, respectively, over the discharge time interval [9]: where t inc and t ext are the PD inception and extinction time, respectively, S G Electr and S cav are ground electrode and surface of the cavity, respectively, dt is the time integration differential, and d → S explains an infinitesimal area over which the field flows.Figure 2 depicts the flowchart for FEA modelling of PD.

Case of Study
In this work, the case study's geometry alongside a representation of its mesh is shown in Fig. 3. Its configuration is the same as that in [3], with the cavity situated right at the middle (x_axis = 3.34 mm, y_axis = 5.10 mm) of the insulation bulk.The PD model and free parameters are, respectively, listed in Tables 1, and 2. Simulations are executed considering a stressing period between 0 and 0.5 h.Because the authors were unable to measure the free parameters of the case study in this research, they were taken from the test arrangement in [66,67] but with slight adjustments.It should be noted that the adjustments were achieved through extrapolations involving the geometrical and material parameters of the test and case models, respectively.The results of those adjustments are presented in Table 1.Such adjustments are critical in PD simulations and are carried to ensure a close agreement between simulation and measurements.PD is simulated under different conditions of cavity radius and cavity surface conductivity listed in Table 2.At the beginning of each condition, it is considered a virgin cavity without previous PD activity.A summary of the PD simulation conditions performed in this paper is therefore presented in Table 2.The magnitude of the electric field strength inside the cavity is considered homogeneous; however, it changes with time owing to sinusoidal excitation during the simulation process; the magnitudes of the inception and extinction fields are determined for each PD event.
The PD model, which was developed in COMSOL, is summarized as follows: i.The magnitude of the electric field strength across the cavity is lower than the magnitude of the ignition field before the occurrence of a PD.The scalar potential and the electric field magnitudes can be determined using the field Eqs. ( 1) and (3) in COMSOL as described in Sect.2.1; ii.If the streamer inception condition is satisfied, the second inception pre-requisite is modelled using Eqs.( 4)-( 10), to check for the occurrence of PD during the subsequent time step using Eq. ( 13); iii.Once the PD inception conditions are fulfilled, the time step is reduced, and the initial surface conductivity of the cavity is enhanced to a higher magnitude defined by Eq. ( 16); iv.Using Eqs. ( 1) and (3), the electric field strength and the electric scalar potential field across the entire discharging sub-domain are recalculated until the local electric field magnitude at the cavity centre recedes below the quenching field value defined by Eq. ( 5); and v.When the discharge process reaches the streamer quenching condition, the time step is reset to a higher value, and the apparent and real charges are estimated using Eqs.( 18) and ( 19) throughout the time interval of the PD activity.

Results and Discussion
The activity of sustained PD in a cavity within the insulation system is modelled and analysed under different cavity sizes and surface conductivity.The associated PD parameters such as the number of PDs per cycle, real charge and inception field were investigated.Alongside the field distribution in the case study, the PD outputs were also seen to match both the simulated and experimental data in [10,14,67].In those papers, the implemented model, presented in Fig. 2, was validated

Electric Field
Figure 4a and b depicts the variation of electric field across the geometry of the case study before and just after a discharge event, respectively.From Fig. 4a, the electric field strength reaches its maximum inside the cavity, specifically, at the dielectric-cavity boundary nearest to the HV conductor.PD charges propagate from this boundary along the straight path, crossing the cavity, towards the cavity wall closest to the ground conductor.When a PD occurs, the cavity field reduces drastically below the inception field (E cav < E inc ) as illustrated in Fig. 4b.At this stage, there is a delay between the Laplacian field established by the resultant field inside the cavity and the HV source.The variation in field strengths in both stages is due to the changes in material properties.

Effect of Cavity Radius on PD
To investigate the impact of cavity size on sustained PD activity, simulations were implemented considering different cavity radius as it is shown in Table 2. Figure 5 depicts the obtained PRPD distributions and the typical electric field strength distributions per cycle for this case study.Statistical data from simulation results are provided in Table 3.
In Figs.5a, c, e and 5(g), it can be observed that all of the PRPD patterns exhibit "turtle-like" distributions.Despite the appearance of large PD pulses (positive and negative) in each half cycle, their phase distributions remain symmetrical.
The number of PDs per cycle rises as cavity size increases, which can be ascribed to a faster electron production rate inside the cavity.Furthermore, the real PD charge magnitude increases with greater cavity radius, which is explained by a long streamer propagation path.Referring to values in Table 3, the peak values of number of PDs per cycle and the PD magnitude were found to be approximately 15.50 and 422.17 pC respectively in the cavity with radius of 1.6 mm.When these results are compared to that of a cavity with radius 0.4 mm, it can be seen that the number of PDs per cycle and the maximum PD charge magnitude have risen by 43.74% and 88.57%, respectively, when the radius increases to 1.6 mm.In the same way, it was observed that while the streamer inception field magnitude increases as the cavity size becomes greater, the PD extinction field magnitude remains constant at any cavity radius.The change in inception field is mainly due to the variations in the properties of trapped cavity gas.However, it can be observed that the field magnitude of the cavity reduces as the cavity increases in size.The maximum cavity field magnitude obtained was 4.02 kV m −1 within the smallest cavity, 0.4 mm radius.With this value, the magnitude of the cavity field strength increased by 37.81% as compared to the 1.6 mm radius cavity.As the cavity size becomes larger, the very first phase of PD event is brought forward (see Fig. 5g).Although the average magnitude of PD charge and the number of PDs per cycle have risen due to increase in the cavity size, and the PRPD structure remains almost uniform, similar to a "turtle-like" pattern.As it is shown in Fig. 5, the PD pulses are distributed in the phases in such a way that rounded hemispheres, like the shell of a turtle, are obtained.This is because the dispersion between the values of the magnitudes of the pulses is low, and the detrapping work function and pressure remain constant along all the simulations.6a also exhibits a "turtle-like" structure very similar to those discussed in Fig. 5a, c, e, g.This is due to the cavity surface conductivity having a negligible impact on the cavity's sustained  4, when the cavity surface conductivity reaches 1 × 10 -13 S m −1 , the maximum PD magnitude and number of PDs per cycle are approximately 58.91 pC and 9.73, respectively.Although the stated peak PD charge amplitude is the same as that of the cavity with surface conductivity of 9 × 10 -13 S m −1 , and its magnitude is reduced by 6.36% when compared to the corresponding case with cavity surface conductivity of 2 × 10 -11 S m −1 .Similarly, number of PDs per cycle is reduced by 7.25% considering the earlier comparison.

PD as a Function of Cavity Surface Conductivity
As the conductivity of the surface wall is increased to 9 × 10 -13 S m −1 and 2 × 10 -11 S m −1 , the PRPD pattern distribution changes from "turtle-like" to "wing-like" as brought in Fig. 6c and e because in both, the dispersion of PD charge 123 magnitudes and the maximum value of PD charge, increase with the conductivity of the surface wall.Among Fig. 6a, c and e, the highest number of PDs per cycle and the maximum PD charge magnitude were found to appear in Fig. 6f and e, respectively.The electric field strength magnitude produced by the surface charge distribution and the number of PDs per cycle reach their maximum values when the cavity surface conductivity reaches 2 × 10 -11 S m −1 .Such outcomes are anticipated, owing to the ability of the cavity surface conductivity, 2 × 10 -11 S m −1 , to contribute in the generation of higher charge (electron) in the discharging cavity.Similarly, the PRPD patterns exhibit a sort of symmetrical distribution in the various half-cycles.Despite the changes in structures of PD patterns seen in Fig. 6a, c, e, the local, streamer inception and extinction fields, respectively, remains constant within the cavity across all values of cavity surface conductivity.

Conclusion
A FEA-based non-deterministic model was employed to study the sustained PD activity in an air-filled spherical cavity within the insulation of an MV XLPE cable.The experimental results as reported in previous studies are in agreement with the simulation results.The discharge analysis revealed that the PRPD patterns and PD statistical quantities vary for each simulation condition.The number of PDs per cycle and their magnitudes were observed to increase as the conductivity of the cavity surface and the cavity size increased.In addition, it was observed that the PD activity increases for both, due to a larger cavity and cavities with greater surface conductivity.The greatest PD magnitude, 477.17 pC, and highest PD number per cycle, 15.50, were obtained in the cavity with 1.6 mm radius.As compared to the same values in a smaller cavity (0.4 mm in radius), these numbers represent increases of 88.57% and 43.74%, respectively.Furthermore, the influence of cavity size on sustained PD activity was investigated and simulations with varying cavity dimensions were run.The results show that all PRPD patterns have "turtle-like" distributions.Despite the appearance of significant positive and negative PD pulses in each half cycle, their phase distributions remain symmetrical.The number of PDs per cycle rises as cavity size increases, which can be ascribed to a larger electron production rate inside the cavity.Furthermore, the size of the apparent PD charge grows by increasing the cavity radius, which is accounted by a lengthy streamer propagation path.Similarly, the approximate values for the PD number per cycle and the highest PD magnitude are 10.49 and 62.91 pC, respectively, as obtained for the highest cavity surface conductivity value of 2 × 10 -11 S m −1 .The stated cavity conductivity value translates to a 6.36% and 7.25% increase in the number of PD per cycle and the peak PD charge magnitude, and when compared to similar values, the cavity surface conductivity corresponds to 1 × 10 -13 S m −1 .These findings are explained by the "winglike" distribution of the PD patterns.The results of electric field distributions and their PRPD pattern as presented in this study can be utilized in the power industry to determine the status of vital HV apparatus, and also for implementing diagnosis and prognosis studies in order to minimize equipment failure and outages.

Fig. 1 Fig. 2
Fig.1Charge migration inside cavity during two PD events with a same b opposite directions of fields[31]

Fig. 4
Fig.4 The variation of electric field across the discharging PD model a before initial PD activity and b immediately after the PD activity for 1.6 mm cavity

Figure 6 depicts
Figure 6 depicts PRPD plots from prolonged PD simulations under different conditions of cavity wall conductivity.With a streamer inception field magnitude of 871.37 kV mm −1 , the distribution of the PRPD pattern in Fig.6aalso exhibits a "turtle-like" structure very similar to those discussed in Fig.5a, c, e, g.This is due to the cavity surface conductivity having a negligible impact on the cavity's sustained

Table 4
Summary of PD quantities for different cavity surface conductivity