Nonequilibrium Thermodynamics and Entropy Production in Simulation of Electrical Tree Growth

In the present work we applied the nonequilibrium thermodynamic theory in the analysis of the dielectric breakdown (DB) process. As the tree channel front moves, the intense field near the front moves electrons and ions irreversibly in the region beyond the tree channel tips where electromechanical, thermal and chemical effects cause irreversible damage and, from the nonequilibrium thermodynamic viewpoint, entropy production. From the nonequilibrium thermodynamics analysis, the entropy production is due to the product of fluxes Ji and conjugated forces Xi: σ = ∑iJiXi ≥ 0. We consider that the coupling between fluxes can describe the dielectric breakdown in solids as a phenomenon of transport of heat, mass and electric charge.


Introduction
All materials conduct electricity to a greater or lesser extent and all suffer some form of breakdown in a sufficiently strong electric field.
For low field strengths the conduction process is ohmic, but as the field strength is increased, the conductivity usually becomes field-dependent; if the field strength is increased further, some form of destructive irreversible conduction takes place [1]: the dielectric breakdown.
Gases and liquids due to the mobility of their molecules do not present permanent damage (channels) when the application of the field ceases, something that does happen in solids. The permanence of these channels in the dielectric breakdown of solids is attributed to thermal effects: local fusion of the material that originates the formation of the channel and gives rise to the transport of ion or electron carriers within it.
Dielectric solids usually exhibit a permanently damaged discharge track. In a recently published work [2], we have made an analysis of the entropy associated with the structure of electric trees through the behavior of four entropies in the field of statistical mechanics in the framework of the DBM Dielectric Breakdown Modell [3].
Historically, the treatment of the simulation of the dielectric breakdown has placed emphasis only on the electric field, particularly in the local field, as the sole cause of the dielectric breakdown, omitting perhaps other factors such as the thermal effect or the presence of moisture. The latter usually affects the structure of the insulating material in high-voltage transformers.
In this work, we will try to correlate the entropy production associated with the simulation process of the structure formed by means of the DBM model with observables of real systems. In other words, we will try to give physical meaning to the simulated discharge within the framework of an "abstract" structure, where the discharge structure, the electric tree, is characterized by the parameters of fractal geometry with fractal dimension df: 1 <df < 2 for branched trees.
Given that the dielectric breakdown is a clear example of a nonequilibrium process, we believe that the originality of this work lies in the combination of a stochastic model based on the DBM for the description of the local field and the purely electrical contribution and others effects associated with the dielectric breakdown process, such as the ionic or electronic conduction, or the presence of moisture on the insulator material, or thermal instability. For this purpose we will use the methods of linear nonequilibrium thermodynamics.

Electrical Conductivity of Dielectrics
The electrical conductivity of the dielectrics can be ionic, electronic or both. It is very difficult, from an experimental point of view, to be able to separate these contributions, particularly to high fields.
The ionic conductivity (σ) is simply due to the migration of positive and negative ions. The basic expression of all electrical conduction is given by: Where ni is the density of carriers i and ei and μi are the corresponding charge and mobility, respectively.
In the usual analysis of low-field ionic conductivity one can write: where σ0 and Φ are determined experimentally (constants in a certain range of T).
Leaving aside the theoretical models that can be adopted to explain conduction in insulating solids, experimental data can be adequately described by Equation (2).
The basic equation of steady-state driving is: Where J is the current density, E is the field strength and σ is the conductivity defined by Equation (1).The current and the potential difference V are experimentally measured. A simple procedure establishes the link between theory and experiment: Where A is the cross section and d the thickness of the sample. Since the rupture path usually extends over a small area of the cross section, the current density in the rupture channel is much greater than that found with the set of Equations (4a) and (4b).
It is well known that the relaxation time associated with the current distribution in the steady state is given by: where ε is the dielectric constant and σ is the conductivity. If the current is due to charge carriers of different species, the conductivity is given by Equation (1). In this case, a relaxation time is obtained: The relaxation time is associated with each species of cargo carriers where the eventual approach to the steady state will be determined by the value of the corresponding time constants. It can be easily verified that small concentrations of relatively immobile species can lead to relaxation times of hours. This is not important if the spatial charge density due to the carriers is small enough to produce a weak distortion of the electric field.
To estimate this effect, let us see, for example: which gives a value of 10 7 V cm −2 for the speed of field variation.

The Model
We think of a dielectric breakdown model due not only to the effect of pure dielectric breakdown but also to the contribution of other effects, with the effect of moisture on the disturbance of the electric field during the pre-breaking process coming to mind.
In a rigorous analysis of the rupture process, we should take into account all possible causes: electric field, thermal effect and transport of matter.
In the formalism of linear nonequilibrium thermodynamics [4,5], the phenomenological equations that would describe the flows and coupled forces that characterize this complex phenomenon would have this form: In the set of Equation (9) we will consider T = Tm. The simulation of the electric field effect was done using the classic DBM model in a 2D grid [2]. The used model is a two-dimensional square lattice (Figure1) in which the opposite sides represent the two electrodes. Breakdown starts at a point of a high local field, and this is usually attributed to electrical conducting inclusions [6]. The conducting inclusions are represented by electrode pins.
For the model, see [2]. When we apply the model according to Equation (10) we obtain several results related to the fractal structure of the electric tree. In this opportunity we will try to relate some physical results of the model, such as the applied field strength and the number of branches of the structure formed with the characteristics of a real insulating system. Then we will relate these observables of applied field and conductivity with the characteristics of a given material (it can be cellulose with oil).
In the same way, we will relate the effect of the transport of water through the insulating material, which reduces the breaking stress with the number of branches (channels) formed.
In the framework of nonequilibrium thermodynamics, the entropy production is a bilinear form of conjugated fluxes and conjugated forces.
For isothermal systems and low values of the electric field E, the current density J is given by Ohm's law (Equation (3) For the system to obey the second law: If we know the conductivity of the insulating material and the diffusion coefficient of the species that diffuses within it for certain field values, in the model given by the value of η we can estimate the entropy production by Equation (11a).
In the next figures (Figures 2 and 3) we see the behavior of the field strength and the number of channels formed according to the η parameter.

Conclusions
The presented model constitutes an attempt to associate the pure dielectric rupture with other effects that can also affect the phenomenon of rupture in insulators, such as the presence of small amounts of water that accumulate or perhaps are formed by degradation of products in the insulators of high voltage. We see through a model based on nonequilibrium thermodynamics presented in Equation (14) as it influences the presence of the contaminant in the production of entropy of the system by the use of the Equation (11a).