Next Article in Journal
BPM Support for Patient-Centred Clinical Pathways in Chronic Diseases
Next Article in Special Issue
An Expert System for Rotating Machine Fault Detection Using Vibration Signal Analysis
Previous Article in Journal
Towards a Mobile Gait Analysis for Patients with a Spinal Cord Injury: A Robust Algorithm Validated for Slow Walking Speeds
Previous Article in Special Issue
Motion Magnification of Vibration Image in Estimation of Technical Object Condition-Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Dielectric Oil with a Low-Cost CMOS Imaging Sensor and a New Electric Permittivity Matrix Using the 3D Cell Method

by
José Miguel Monzón-Verona
1,2,
Pablo Ignacio González-Domínguez
1,2,*,
Santiago García-Alonso
3 and
Jenifer Vaswani Reboso
4
1
Electrical Engineering Department (DIE), University of Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
2
Institute for Applied Microelectronics, University of Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
3
Department of Electronic Engineering and Automatics (DIEA), University of Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
4
Process Engineering Department (DIP), University of Las Palmas de Gran Canaria, 35017 Las Palmas de Gran Canaria, Spain
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(21), 7380; https://doi.org/10.3390/s21217380
Submission received: 23 September 2021 / Revised: 3 November 2021 / Accepted: 4 November 2021 / Published: 6 November 2021

Abstract

:
In this paper, a new method for characterizing the dielectric breakdown voltage of dielectric oils is presented, based on the IEC 60156 international standard. In this standard, the effective value of the dielectric breakdown voltage is obtained, but information is not provided on the distribution of Kelvin forces an instant before the dynamic behavior of the arc begins or the state of the gases that are produced an instant after the moment of appearance of the electric arc in the oil. In this paper, the behavior of the oil before and after the appearance of the electric arc is characterized by combining a low-cost CMOS imaging sensor and a new matrix of electrical permittivity associated with the dielectric oil, using the 3D cell method. In this way, we also predict the electric field before and after the electric rupture. The error compared to the finite element method is less than 0.36%. In addition, a new method is proposed to measure the kinematic viscosity of dielectric oils. Using a low-cost imaging sensor, the distribution of bubbles is measured, together with their diameters and their rates of ascent after the electric arc occurs. This method is verified using ASTM standards and data provided by the oil manufacturer. The results of these tests can be used to prevent incipient failures and evaluate preventive maintenance processes such as transformer oil replacement or recovery.

1. Introduction

In electrical power systems, transformers are one of the elements with the highest economic cost, accounting for around 60% of the investment in high-voltage substations. This requires a set of monitoring and diagnostic techniques that affect the life cycle of these important elements. These techniques include the following: dissolved gas analysis, oil quality test, infrared thermography test, power factor test, dielectric dissipation factor test and dielectric oil breakdown test, among others [1,2].
In this paper, an analysis of the quality of the oil is made using a combination of electrical, physical and chemical tests [3]. The transformer oil corresponds to an oil sample without hours of use.
The most important and common tests are related to the dielectric breakdown voltage (BDV), water content, acidity and color. There are three main standards for the determination of the dielectric breakdown voltage of insulating liquids: ASTM D1816–12 (2019), ASTM D877/D877M–19 and IEC 60156:2018. In this paper, UNE-EN 60156, which is based on IEC 60156:2018, was followed.
The results of these tests are used to prevent incipient failures and evaluate preventive maintenance processes, such as transformer oil replacement or recovery [4].
On the one hand, mineral oils in transformers play an important role as an element of electrical insulation between the parts under voltage and, on the other hand, they help to evacuate the heat generated due to hysteresis losses and eddy currents in iron, as well as the losses due to the Joule effect in the transformer coils. This last condition requires the oil to have a high thermal conductivity and a low coefficient of dynamic viscosity.
The breakdown strength of dielectric oils for transformers will depend on the nature of the impurities present in their solid or gaseous state. Oil analysis is important for extending the life of the transformer.
The state of knowledge of dielectric breakdown voltages in insulating liquids is less developed than in the case of gas and solid dielectrics. The studies carried out may in some cases be contradictory [5].
Among these studies are those that explain the dielectric breakdown voltages of liquids based on an extension of the dielectric breakdown voltages in gases, in turn based on the avalanche ionization of the atoms caused by electron collisions in the applied field [6].
Dielectric breakdown voltages in different temperature ranges show little dependence on temperature. This suggests that the cathode emission process is one of field emission rather than thermionic emission [6].
Electronic theory predicts the relative magnitudes of the dielectric breakdown voltages well, but not the times at which these breakdowns occur in the insulating liquid. This last aspect—the temporal aspect—is partly explained by the presence of polluting particles inside the insulation. These give rise to local breakdowns, which in turn lead to the formation of small bubbles with a much lower dielectric strength and, hence, finally lead to breakdown.
Other phenomena that explain electrical breakdown include the electroconvection of dielectric breakdown, dielectric liquids subjected to high voltages and electrical conduction resulting mainly from charge carriers injected into the liquid from the electrode surface. The resulting space charge gives rise to Coulomb’s force which, under certain conditions, causes hydrodynamic instability, creating an eddy motion of the liquid which yields a convection current.
Thus, the charge transport will be largely via liquid motion and not ionic drift. The key condition for the instability onset is that the local low velocity exceeds the ionic drift velocity [5].
In most of the studies referenced in the bibliography dealing with analyses of the dielectric breakdown voltage in dielectric oils, approximate analytical equations are used, or numerical methods are employed based on a differential formulation such as the finite elements method (FEM).
In the present study, we propose the finite formulation (FF) method [7], together with the cell method (CM) [8,9] as an associated numerical method to analyze this type of device. Using this methodology, we consider the global magnitudes associated with space-oriented elements such as volumes, surfaces, lines and points of the discretized space, as well as temporal elements, instead of field magnitudes associated with independent variables, i.e., spatial and temporal coordinates [7].
In addition, equations of the constitutive type—equations of the medium—are clearly differentiated from those of the topological type—equations of balance. In FF, the physical laws that govern the electromagnetic equations are expressed in their integral form. In this way, the final system of equations is obtained directly, without the need to discretize the equivalent differential equations [8].
In the CM, the topological equations obtained directly from Maxwell’s laws are exact (balance equations), while the constitutive equations obtained from the discretization process are approximate. In the latter case, source-type quantities defined in the elements of the dual mesh must be related to the configuration quantities corresponding to the elements of the primal mesh [10]. Field magnitudes and the physical properties of the medium are assumed to be constant, at least in the primal mesh. This ensures that the discrete equations are consistent with the continuous constitutive equations, in the sense that the discrete constitutive equations approximate the continuous constitutive equations with an error that decreases with the mesh size [11].
Most of the research papers in the CM literature focus on the construction of discrete constitutive equations. Among those related to a quasi-electrostatic problem in 2D with plane symmetry, is [12]. In this paper, an isotropic and anisotropic electrostatic field is studied by means of the CM. In [13], the electrostatic problem is studied in 2D with plane symmetry. The constitutive equation is used with two approximations. The first approach assumes a uniform field with a triangular base inside each primal cell, and the second approach, which is more general, assumes the uniformity of the fields in subregions of each primal cell, with quadrilateral bases. In [14], a 2D analysis with axial symmetry (axisymmetric) is performed for a quasi-electrostatic problem regarding a gas-insulated line for an ITER neutral beam injector.
In [15], an electrostatic induction micro-motor is studied, using the CM in 2D.
In the literature, to the best of the authors’ knowledge, the 3D cell method has not hitherto been used to simulate dielectric breakdown voltage tests on transformer oils as a quasi-electrostatic problem. In this article, we propose the use of a geometric structure for the electrical permittivity constitutive matrix, analogous to the matrix that appears in [16] for an electromagnetic problem in 3D to calculate eddy currents.
The advantage of using the same geometric structure in the constitutive matrices of conductivity and permittivity in the quasi-electrostatic problem in 3D is that it reduces the complexity of the programmed source code and the execution times. This is because the constitutive matrix is calculated in the assembly of the system of equations, only once for each tetrahedron. The electrical conductivity and electrical permittivity properties of each tetrahedron are multiplied by the common matrix. This is done, element by element, until the final system of equations is complete.
In the present study, a new constitutive matrix is formulated that, using the CM, relates the differences in electrical potentials (magnitudes of configuration) due to primal mesh edges with the electrical flux (magnitudes of source) due to dual mesh planes. The magnitudes of configuration are associated with the edges of a primal mesh made up of tetrahedra, and the source-type magnitudes are associated with the surfaces of a dual mesh (the control volume) obtained in a barycentric division of the primal mesh.
This paper presents an experimental study of the dielectric strength of transformer oil based on the IEC 60,156 standard [17]. Our contribution consists of characterizing the behavior of the oil an instant before and after the electric arc rupture, combining a low-cost complementary metal oxide semiconductor (CMOS) imaging sensor and a new electrical permittivity matrix M ε , using the 3D CM [7,15]. In the standard test, only the effective value of the dielectric breakdown voltage is obtained. However, the information on the distribution of Kelvin forces [18] an instant before the dynamic behavior of the arc begins is lost, as is the information on the gases that are produced an instant after the moment of breakdown via the electric arc in the oil.
This last aspect was analyzed by recording images of the movement of the gas bubbles that are produced within the oil. This also allowed the diameter of these bubbles to be measured. The measurement of their magnitudes was used to indirectly obtain the viscosity of the oil. The physical property of viscosity could be obtained by analyzing the post-arc images using an equation to predict the terminal velocity of the rise of isolated bubbles in Newtonian liquids [19].
The data obtained with the sensors and the results of the simulations complement each other and offer information that would otherwise be lost when strictly following the standard test.
The use of low-cost camera systems in remote-sensing applications is not new. The use and study of low-cost cameras for engineering and scientific applications was addressed in detail in [20]. In [21] a study of the corona effect in aeronautical applications was performed, using low-cost Raspberry-Pi-type cameras. In [22], an array of single-board computers produced by Raspberry Pi, and their associated 8 MP cameras, was used at the University of Cambridge to capture the images required for particle image velocimetry analysis or analysis of the correlation of digital images.
There are two main objectives in this paper. One is to present a new matrix of electrical permittivity, M ϵ , that predicts the electric field before and after the electric rupture occurs. The other is to measure the kinematic viscosity of the dielectric oil using a low-cost CMOS imaging sensor to measure the distribution of bubbles, their diameters and their rates of ascent after the electric arc occurs. In addition, experiments were performed to estimate the dielectric breakdown voltage in order to obtain the boundary conditions for CM and FEM simulations. In this way, both objectives are related.
This paper has been divided into the following sections. Section 2 determines the distribution of E and E 2 in the dielectric strength test by applying the 3D CM using the new constitutive matrix M ϵ . Section 3 describes the low-cost 8 MP CMOS imaging sensor used in the experimental studies. Section 4 presents in detail the numerical results of the CM simulations with M ϵ . This matrix is verified by comparing the results obtained with those from the FEM analysis. Finally, Section 5 presents the experimental setup of the oil testing device and the results that were obtained. In addition, it lays out the test procedure for the dielectric strength and kinematic viscosity of the transformer oil and establishes the theoretical basis of the new procedure for the determination of the dynamic viscosity of the oil. Finally, the data obtained on the kinematic viscosity using the proposed methods are verified by comparing them with the manufacturer’s data.

2. Distribution of E and E 2 in the Dielectric Strength Test

The formation of the electric arc and the subsequent bubble formation are highly dependent on the estimation of the electric field distribution. Furthermore, its gradient determines the forces per unit volume that act on polluting particles and microbubbles within a dielectric [23,24].
It is important to determine the distribution of the field E and the gradient of its square E 2 as these are factors that determine the DBV.
The following section explains the proposed method for obtaining this field distribution using the 3D CM as an alternative method to the FEM. In this section, a new electric permittivity matrix M ϵ is proposed for use in the CM.

2.1. New Constitutive Matrix M ε . Discrete Constitutive Equations of Transformer Oil in the Finite Formulation

The electrical constitutive equation in transformer oils is a complex equation based on the Fowler–Nordheim theory [25]. In most dielectric materials, the conduction current of free carriers is relatively low, since their conductivity is usually several orders of magnitude lower than that of a metal or semiconductor. In new transformer oils at 50 °C it is usually of the order of 1 × 10−13 S/m, and in used oils of the order of 1 × 10−11 S/m [26].
In this paper, a conductive-type model is considered and the conductivity is assumed to be the same throughout the volume of the oil, where the volumetric current density J is directly proportional to the electric field E [25].
Taking into account the non-zero conductive properties of oil and the fact that it is subjected to an electric field, in the CM the current through the material can be described by its constitutive equation of current flow as a function of potential differences associated with the edges of the primal mesh ei i = 1:6, as shown in Figure 1. The constitutive equation shown in Equation (1) was developed in [16].
I ˜ f = M σ · U .
The electrical constitutive matrix M σ is given by the expression M σ = σ ν S ˜ i · S ˜ j   i ;   j = 1 : 6 . The matrix M σ is a function of the electrical conductivity of each tetrahedron; its volume and the dot product of the surface vectors correspond to the dual planes (green planes) of primal edges, edges ei and ej i j = 1:6, as shown in Figure 1.
In dielectrics, however, when the applied electric fields are variable with time, a new contribution to the free current appears, the so-called displacement current. This appears when there is variation in the electric flux with respect to time. The displacement current has two terms: I ˜ d = ψ ˜ t = ( ψ ˜ 0 + ψ ˜ P ) t . The first term within the parentheses only depends on the potential difference in a vacuum. It is independent of the characteristics of the material. The second term depends on the insulating material used. This depends exclusively on the polarization of the dielectric—in this case the transformer oil—and contains the response of the material, which will be different according to the polarization mechanisms that occur for each stimulus of the net applied potential difference (due both to free charges and to those of polarization). We can group both terms with the total empty and material medium permittivities in ε = ε r ε 0 , with the constitutive equation given in Equation (2).
ψ ˜ = M ε 0 U + M χ P U = M ε U .
Equation (2) relates the differences in electric potential U associated with the edges of the primal mesh (magnitudes of configuration) with the electric fluxes ψ ˜ of the dual planes S ˜ i   i = 1 : 6 , as shown in Figure 1. In Equation (2), the matrix M ε is the new electrical constitutive matrix proposed in this work. Given the analogy with Equation (1), the following expression is proposed:
M ε = ε ν S ˜ i · S ˜ j   i ;   j = 1 : 6 .
where M ε is a function of the electric permittivity of each tetrahedron, ε = ε 0 ε r , its volume ν and the dot product of the surface vectors that correspond to the dual planes of primal edges, edges ei and ej. The value of ε r depends on the type of oil. At 20 °C, this lies between the values of 2.1 and 3.5 [27]. Therefore, the displacement current is obtained from Equation (4).
  I ˜ d = ( M ε U ) t .
The total current will be given by the sum of both contributions, according to Equation (5). The total current I ˜ t is represented for the dual plane s ˜ 3 in Figure 1.
I ˜ t = I ˜ f + I ˜ d .

2.2. Maxwell’s Laws in Finite Formulation Applied to Transformer Oil

Maxwell’s laws, applied in their finite formulation in the dielectric strength test are, first of all, the laws corresponding to the configuration magnitudes. According to [8], these are as follows.
(a)
Gauss’s law for the magnetic field, Equation (6)
D ( Φ ) = 0 ,
where D is the volume–face incidence matrix of the primal mesh, which is equivalent to the standard divergence operator. The magnitude ϕ represents a vector with all the magnetic fluxes associated with the four faces of the primal mesh tetrahedron if i = 1:4, as shown in Figure 1.
(b)
Faraday’s law of induction, Equation (7)
C · U = ( ϕ ) t ,
where C is the face–edge incidence matrix of the primal mesh, which is equivalent to the standard rotational operator. U is a vector of potential differences extended to all edges of the primal mesh and t is time.
The next two laws correspond to the laws that operate with magnitudes of energy.
(c)
Generalized Ampere’s law, Equation (8)
C ˜ · F ˜ = I ˜ f + ( ψ ˜ ) t ,
where C ˜ is the face–edge incidence matrix in the dual mesh, the vector F ˜ is a vector of magnetomotive force associated with all the edges of the dual mesh, I ˜ f is a vector of electric currents extended to all planes of the dual mesh and, finally, ψ ˜ is the electric flux due to the polarization of the dielectric associated with the faces of the dual mesh.
(d)
Gauss’s law of the electric field, Equation (9)
D ˜ · ψ ˜ = Q f ,
where D ˜ is the incidence matrix of the volumes–faces of the dual mesh and Q f is the charge contained in each dual volume.
Finally, if the divergence is applied to Equation (8), taking into account Equation (9), Equation (10) is obtained.
D ˜ · I ˜ f + ( Q f ) t = 0 .

2.3. Maxwell’ Laws and Constitutive Equations

In this section, the constitutive equations and Maxwell’s laws given in Section 2.1 and Section 2.2 are combined to obtain the final equation in the time and frequency domain. The electrical scalar potential is used for this purpose, and it significantly reduces the number of unknowns.
As we mentioned in Section 2.1, we assume that the electrical conductivity of the oil is low, of the order of magnitude of σ = 1 × 10 13 S/m. The permittivity of oil is   ε = ε r 8.854187818 × 10 12 F/m. Three time constants can be considered to fit the type of problem to be solved. The first is defined as the charge-relaxation time, τ e = ε σ ε r 88   s . The second is the electromagnetic time constant, τ e m = l c 10 12   s , where l = 10 cm, the characteristic length of the domain. The constant c is the speed of light. The third constant is the magnetic time constant, τ m = τ e m 2 τ e 10 22   s . We can consider that if the frequency is 50 Hz, τ ≈ 20 ms and the conditions for the field to be quasi-electrostatic are β = ( τ e m τ ) 2 1 and τ m < τ e m < τ e . For further information on quasi-static laws and time rate expansions, see [28]. In this way, the field can be considered to be quasi-electrostatic, and Equation (7) can be written as
C · U 0
Since the field is, in this situation, almost electrostatic, it is possible to work with a single electric potential U = G φ , where φ is an electric scalar potential. This is imposed on the surface of the electrodes, as shown in Figure 2. Taking into account the constitutive Equations (2) and (9), this is as follows.
D · M ε ( G φ ) = Q f
Taking into account Equation (10), and substituting in this equation the free volumetric electric charge Qf from Equation (12) and the free current I f from Equation (1), using U = G φ , Equation (13) is obtained.
D ˜ · M σ G φ + ( D ˜ · M ε G φ ) t = 0
Equation (13) corresponds to the differential equation derived in [28]. If the electrodes work at a frequency of f = 50 Hz, with an angular frequency of ω = 2π50 s−1, the final equation in the frequency domain will be Equation (14). It can be observed that this equation involves electrical permittivity and electrical conductivity. This is the equation to be programmed, together with the global electrode current I ˜ t , which is calculated using Equation (15) where I c is an incidence vector of the relative cut between the edges of the oil volume mesh and the surface of one of the electrodes, as shown in Figure 3. The sum of all currents in that cut is equal to the total current entering or leaving at each of the electrodes.
D ˜ · M σ G φ + j ω D ˜ M ε G φ = 0
I ˜ t = I c · ( M σ G + j ω D ˜ M ε G ) φ
The matrix representation of both equations is shown in Equation (16).
[ G t M σ G + j ω G t M ε G 0 I c M σ G + j ω I c M ε G 1 ] [ φ I ˜ t ] = [ 0 0 ]
The unknowns are all the potentials φ of the primal mesh nodes and the global magnitude of the current I ˜ t associated with the surface of one of the electrodes.
Note that the total current has been defined twice by Equations (5) and (15) because Equation (5) shows the two main components of the total current, i.e., the displacement current and the conductive electric current, while Equation (15) details its explicit composition based on its physical and geometric properties.
In the matrix system (16), the first row is uncoupled from the second. However, we have preferred a compact equation system that includes the total intensity of the electrode. This avoids an additional post-processing calculus. It is also true that we have increased the dimensions of the system by one degree of freedom, corresponding to the total intensity through the electrode. In this way, by solving a single matrix system all the unknowns (degrees of freedom) are obtained at the same time, without post-processing.
Boundary conditions are stabilized on the electrode surfaces. All the nodes on one electrode surface have a value of the electric potential equal to zero and all the nodes on the surface of the other electrode have the dielectric breakdown voltage.
The electrodes have an axisymmetric geometry, but the oil container is cubic, making the global model non-axisymmetric. Hence, we must perform a 3D analysis. Furthermore, the new proposed electrical permittivity matrix is a 3D matrix, and these calculations give generality to the matrix M ε .

2.4. Kelvin Polarization Forces in Dielectric Materials

The force suffered by a supposed spherical particle or a micro-bubble of gas in suspension of radius r and relative permittivity ε , in a liquid with relative permittivity ε o i l , in the presence of an electric field E = ( E x , E y , E z ) , is calculated using the Kelvin polarization force formula, according to [23], as shown in Equation (17).
F e = 0 r ϵ 0 ( ϵ ϵ o i l ) E E   d v .
If the x component of F e is developed, Equation (18) is obtained.
F e x = ε 0 ( ε ε o i l ) ( E x E x x + E y E y y + E z E z z ) .
It is developed in a similar way for y coordinates as for z coordinates. As × E 0 , Equation (17), according to [24], is simplified and Equation (19) is obtained.
F e = 0 r ε 0 ( ε ε o i l ) 1 2 E 2   d v .
If ε > ε o i l , this force tends to move the particles to the area where the field is stronger, aligning them and forming a bridge that makes it easier for the current to cross the liquid dielectric via that path. The field in the area of the particles increases and the breakdown value is reached [5].
If the number of particles is not sufficient to bridge the gap, the particles give rise to a local field enhancement, and if the field exceeds the dielectric strength of the liquid, local breakdown will occur near the particles, thus resulting in the formation of gas bubbles which have a much lower dielectric strength and hence, finally, leading to breakdown.
It is important to highlight the dependence of the electric field strength and its gradient and hence the importance of the estimation of the distribution of this electric field and its gradient. This is achieved by applying the method proposed in this paper: the CM together with the proposed new electrical permittivity matrix M ε .

3. The Low-Cost 8 MP CMOS Imaging Sensor

3.1. Camera Features

Sony IMX219 cameras, version v2.1, with 8 megapixels and CMOS sensors, were used in all the experiments. They can be seen in Figure 4, and a detail of the central camera is shown in Figure 5. The plates of the camera weigh 3 g and measure 25 × 24 mm. They are much lighter than most other computer vision cameras and significantly cheaper, costing around USD 30. Each camera requires its own Raspberry Pi to control all its parameters. The small size of the camera and its various elements, make it easy to assemble at a low cost and to connect to its paired computer, compared with other more traditional computer vision cameras or digital cameras.
The camera, as a sensor, must be perfectly linear throughout the detection area. The ideal sensor produces the same response to a photon striking the surface in the center as it does to one at the edge of the detection area. In [29], the data plotted for the Raspberry Pi camera showed that, for the standard camera operating mode, the gain did not allow radiance values higher than 0.4 W m 2 s r , but the signal produced at these lower light intensities was very linear. Cameras without a near-infrared (NIR) filter have different gains with similar compensations [29].
The cameras can capture fixed photos at a resolution of 3280 × 2464 pixels or high-definition video at a resolution of 1920 × 1080 pixels and a rate of 30 frames per second (fps). Higher frame rates have been achieved via custom manipulation of the camera interface using the OpenCV software library [30], with 640 × 480 resolution and capture rates of up to 87 fps obtained in all the experiments carried out.

3.2. Camera Triggering with Raspberry Pi

A Raspberry Pi is a small board of a size similar to a credit card which has an ARM microprocessor with power of up to 1 GHz, integrated in a Broadcom BCM2835 chip. One of the features of the Raspberry Pi is the electronic general-purpose input/output (GPIO) interface [22].
This is a built-in 40-way connection that can be used to send and receive 3.3 VDC signals, as shown in Figure 6. It was used in the configuration of the cameras described in this paper to activate three cameras simultaneously. They were controlled through an interface programmed in Python, in the Raspberry Pi operating system. The operating system used by the Raspberry Pi is called Raspbian and it is a distribution of the GNU/Linux operating system based on Debian and therefore free to use. Synchronization of all the cameras was achieved by electronically and simultaneously activating the cameras through a switch that triggered them through the GPIO inputs on each Raspberry Pi. These were connected to each other and to a central computer to collect the data, as shown in Figure 4.
The camera–Raspberry Pi system is ideal for performing image analysis in near real time, using the numerical analysis libraries in Python, such as the fundamental package for scientific computing with Python, NumPy [31]. There is also extensive support for the OpenCV computer vision package.
A method of interacting with the Raspberry Pi camera board through a Python interface via the picamera software library was used [32]. This high-level software provides the user with a library of Python commands to manipulate the GPU on the Raspberry Pi and control the camera settings through software manipulation of the data captured by the CMOS image sensor. The Raspberry Pi computers are configured by a startup program that runs automatically on each of the three Raspberry Pi computers paired with their corresponding cameras, i.e., three Raspberry Pi computers and three cameras, as shown in Figure 4.

4. Numerical Results of the Simulations in CM with M ε vs. FEM

The system of equations shown in Equation (16), in a sinusoidal steady state, was programmed in C++. As the matrices are sparse and of large dimensions, the Krylov subspace was used. These algorithms were implemented in the PETSc software [33].
In particular, the linear solver employed was the generalized minimal residual algorithm (GMRES).
The main characteristics of the computer used to make the simulations were: computer model: X399 AORUS PRO; architecture: x86_64; total memory: 128 GB; processors: 24; cpu: 2185.498 MHz; thread(s) per core: 2; core(s) per socket: 12.
The numerical validation of the proposed constitutive equation as shown in Equation (2) and the system of equations shown in Equation (16) was carried out by comparing it with the standard FEM and its implementation within the GetDP software [34]. A very dense reference mesh was used in 3D with a number of tetrahedra equal to 2,790,589 volumes and 487,435 nodes. The Gmsh software [35] was used for the mesh and the data visualization. The nodes determine the number of unknowns in the system of equations, as shown in Equation (16).
A section of this mesh and the solution of the potential distribution corresponding to a potential difference between the electrodes of 19,488 V for the dielectric breakdown voltage, can be seen in Figure 7.
The results corresponding to the spatial distribution of the densities of the conductive and displacement currents are shown in Figure 8 and Figure 9, respectively.
The results of the distribution of the force density per unit volume scaled with the coefficient 1 2 ( ε ε o i l ) are shown in Figure 10 and Figure 11. These graphs were obtained by cutting the simulation result with the planes 1 X + 0 Y +0 Z − 0.0489 = 0 and 0 X + 1 Y + 0 Z − 0.035 = 0, respectively. Note that the maximum force occurs in a ring around the center.
Figure 12 shows the capture of the electric arc recorded by cameras A, B and C, from the left, center and right, respectively. The images show the breakdown point of the arc in an area that coincides with the one estimated in the model of Equation (19).
Furthermore, the usefulness of calculating Kelvin forces an instant before the dielectric breakdown voltage is reached, using CM, helps in understanding where the dielectric breakdown voltage occurs. Figure 10 and Figure 11 show that maximum forces are located in a ring around the electrode, but not at its center. This is coherent with the experimental images shown in Figure 12, obtained with the imaging sensors.

4.1. Validation of the Numerical Simulations

This section presents the numerical experiments to validate the results obtained by the proposed method, i.e., the cell method CM and the permittivity matrix M ε from Equation (3). This is compared with the FEM in 3D, with a high number of elements (tetrahedra). The latter adapt well to the surface of the electrodes. The geometry of the problem has planes of symmetry. In the end, this always becomes a 3D problem. The experiments solved the entire problem in 3D.
Table 1 summarizes the geometric and physical properties of the dielectric oil used in the numerical experiments. Each type of experiment E1, E2 and E3 was subdivided into two data analyses corresponding to cut A and cut B, which are two characteristic areas of the electrodes, as shown in Figure 7. Between these two zones, the potential gradients differ greatly and serve to compare the two numerical methods of CM and FEM–getdp using the same mesh.
The number of points in the cut used in all the experiments was 180. The analysis of the metrics was carried out on this basis.
Three numerical experiments were designed: E1, E2 and E3. These consisted of discretizing the electrodes and the volume around them, as shown in Figure 7. Each experiment had a different mesh density, as indicated in Table 1. The objective of these experiments was to check the convergence of the CM numerical method, as well as to determine the error produced when comparing it with another numerical reference method, FEM–getdp.

4.1.1. Results of Experiment E1

Figure 13 shows the results of the simulation using the CM and the FEM–getdp. These are the distributions of electric potential obtained in cuts A and B. The lines with a higher slope correspond to cut A and those with a lower slope to cut B.
Figure 14 shows the distribution of the electric field strength modules. A greater electric field is observed for cut A.
In experiment E1, a low-density mesh was used. In Figure 13, it can be seen that the voltage values in the FEM–getdp and CM are almost coincident for cut A. However, they differ more for cut B. In Figure 14 there is a notable difference between the values of the electric field given by the FEM–getdp and those given by the CM, in both cut A and cut B, because the mesh is not very dense.

4.1.2. Results of Experiment E2

Figure 15 and Figure 16 represent the results for the electric potential and the modules of the electric field strength, respectively, corresponding to experiment E2.
In experiment E2, a medium-density mesh was used. In Figure 15, it can be seen that the voltage values for the FEM–getdp and CM are practically the same. In Figure 16, there are still differences between the electric field values given by the FEM–getdp and CM, in both cut A and cut B.

4.1.3. Results of Experiment E3

The maximum convergence corresponds to Figure 17 and Figure 18, with a total number of tetrahedra of 2,790,589 for experiment E3, as shown in Table 1.
A high-density mesh was used in experiment E3. In Figure 17 and Figure 18, there are no differences in the values given by the FEM–getdp and CM for either the voltage or the electric field.
Therefore, it is observed that as the mesh density is increased, the two methods give coincident solutions. This confirms the validity of the new Mε matrix proposed in this paper for the CM.

4.2. Metrics of Numerical Experiments

To validate the proposed method ( M ε using the CM), a series of comparisons was established between the results obtained with the CM and with the FEM–getdp in the numerical experiments performed.
Metrics are applied to the validation of a model against a reference or pattern. In this case, the model to be validated was the results obtained with the CM and the reference or pattern was the results obtained with the FEM. The CM and FEM are approximate numerical methods, therefore not exact. In both methods a tolerable error is pre-set.
There are various sources of error. One is the truncation of the figures and the accumulation of errors due to the numerical operations performed. Another, which is the one that most affects our problem, is the layout of the cuts. The proximity of the cut to the nodes of the mesh makes the calculated value at the cut more accurate. In contrast, when the cut moves further away from the node, it will be necessary to obtain interpolated values, thus producing a greater error. This happens in both the FEM and CM.
In this study, an analysis was carried out in the FEM with a very dense mesh. These results were used as a reference. Different mesh densities were established in the CM. With the results obtained using the CM, different comparisons were made with the results of the high-mesh-density FEM model. Various statistical indicators (metrics) were used to check the validity of the CM versus the FEM.
The comparisons made are shown in Table 2.
Of the various statistics that can be used to measure the goodness of fit of a model, the following were chosen for the present study: the coefficient of determination (R2), the root mean square percentage error (RMSPE), the mean absolute percentage error (MAPE) and the percentage bias (PBIAS). Table 3 and Table 4 show the range of these statistics, their optimal values and the value obtained in each the comparisons carried out. When the value obtained is closer to the optimum of the statistic, the analyzed mathematical model has a better goodness of fit.
Figure 19 shows various error histograms for some of the comparisons carried out. Following error theory, an ideally random error distribution should have its mean around zero and a Gaussian or normal distribution.
If the comparisons are grouped following the order {C1, C5, C9}, {C2, C6, C10}, {C3, C7, C11} and {C4, C8, C12}, it is observed, in a generalized way, that increasing the mesh density improves the model for any of the numerical experiments analyzed.
Comparing Figure 13, Figure 15 and Figure 17 with Figure 14, Figure 16 and Figure 18, it can also be seen that the voltages in {C3, C7, C11} and {C4, C8, C12} were modeled, in relative terms, much better than the electric fields in {C1, C5, C9} and {C2, C6, C10}.
The values of the coefficient of determination (R2) indicate a good fit of the data for all the comparisons except those cases where the mesh was not very dense. The RMSPE, MAPE and PBIAS values are relatively high in the comparison. Nevertheless, all the indicators are in the optimal range. Even the largest error value of 0.36% for the RMPSE, as shown in Table 3, is a more than acceptable value.
The distribution of the errors, following error theory, conforms to a normal distribution centered on the zero value, as shown in Figure 19.
The detailed formulations of the metrics used appear in the Annex.

5. Experimental Setup of Oil Testing Device

5.1. Procedure for Dielectric Strength Test and Kinematic Viscosity of Oil

Dielectric oil strength test equipment [17] is known as a spark meter or dielectric oil tester. The equipment is in the form of a closed metal cabinet with a handle to facilitate its transport, as shown in Figure 20. Inside is the test cell, where the oil sample to be tested is placed. The tester is equipped with detachable electrodes for testing with different types of electrodes. The cell has a magnetic stirrer in the lower part of the oil deposit for perfect homogenization of the oil sample to be tested, and the rotating part is a small Teflon magnet, as shown in Figure 20.
The technical characteristics of the equipment are summarized in Table 5.

5.2. Performed Experiments

In order to obtain data prior to using the proposed method, a series of standard tests were performed. All the tests carried out on the transformer oil sample are shown in Table 6, together with the test methodology [39,40,41,42,43,44].
In view of the obtained results, it can be confirmed that the parameters analyzed in the oil used complied with the specifications recommended by the UNE 60296:2012 standard [44].
In addition to the previous tests, a set of additional experiments were carried out with the objective of carrying out a strength test and estimating the viscosity of the oil by measuring the diameters of the bubbles produced and their velocity. All tests used a sample with an oil volume of 450 mL. Tests were carried out at different rates of increase of the applied voltage: 0.5, 2, 3 and 5 kV/s. A total of 32 dielectric strength tests were conducted with the capture of the corresponding images. A waiting time of approximately 5 min was established between each cycle to guarantee the disappearance of possible bubbles. In each experiment, the rupture voltage was indicated as well as the number of frames per second (fps) necessary to calculate the velocity of the bubbles produced after the arc occurred. The average capture recording time was 16 s for each experiment, giving an average number of images per experiment of 1344, with a mean time between frames of 0.012 s. As an example, Table 7 shows the results corresponding to the voltage increase rate of 0.5 kV/s. Electric arc at the instant of the breakdown of the electrical strength is shown in Figure 21.
Here, fps is frames per second and R, C and L refer to the right, central and left cameras, respectively. In Figure 4, R, C and L correspond to A, B and C, respectively.
We used a particular experiment to establish the boundary conditions. More specifically, in these simulations, a dielectric breakdown voltage close to 20 kV was set as the boundary condition, as can be seen in Figure 13, Figure 15 and Figure 17. We considered this to be a number consistent with the experimental tests performed. Table 7 shows the dielectric breakdown voltage of a particular set of experiments for a voltage ramp of 0.5 kV/s.

5.3. Experimental Measurement of the Dynamic Viscosity of Oil

5.3.1. Theoretical Basis

In [19], an equation is proposed to predict the terminal velocity of the rise of isolated bubbles in Newtonian liquids.
The formulation combines a balance of forces obtained from the limit layer theory for spherical bubbles with an analytical equation from a mechanical energy balance. This is the model denoted model 2 (Mod2). The equation proposed in [19], which gives an approximation of the terminal velocity of the bubble inside the oil, is shown below as Equation (20).
V T = 1 ( 1 V T 1 2 + 1 V T 2 2 ) 0.5   ,
where V T 1 is calculated using Equation (21).
V T 1 = V T p o t [ 1 + 0.73667 ( g D i ) 0.5 V T p o t ] 0.5
V T p o t is calculated using Equation (22)
V T p o t = 1 36 Δ ρ D i 2 μ o i l
and the term V T 2 is calculated using Equation (23).
V T 2 = [ 3 σ b ρ o i l + g D i Δ ρ 2 ρ o i l ] .
The procedure proposed in this paper consists of solving for the dynamic viscosity μ o i l , knowing the terminal velocity of the bubble and its diameter. Using these equations, we start from the experimentally measured terminal velocity V T , the diameter of the bubble Di from the images of the sensors, the difference in gas and oil densities ∆ρ, the acceleration due to gravity g and the viscosity μ o i l . The surface tension of the bubble σ b is obtained from the Laplace law, as expressed in Equation (24).
σ b = P R 2 , P = P 0 + ρ o i l · g · h o i l ,
where R = D i 2 and P is the absolute pressure, which is the atmospheric pressure at a height of 326 m above sea level plus the mean pressure of the oil column.
From these equations, the viscosity is found as a function of the terminal velocity V T obtained from the second-degree algebraic equation given by Equation (25),
x 2 ( V T 2 k 4 k 2 2 k 2 2 ) + x ( V T 2 k 4 k 5 k 5 ) + V T 2 = 0 .
where x = 1 μ o i l is the unknown. The constants k i are as follows.
k 2 = Δ ρ D i 2 36 ,
k 3 = 0.73667 ( g D i ) 0.5 ,
k 4 = 1 V T 2 2 ,
k 5 = k 2 · k 3 .
A simpler approximation than Equation (25) is to assume that the acceleration of the bubble is equal to zero. It was found in all the experiments that the terminal velocities of the bubbles were substantially constant. This simplifies the equilibrium equation between the buoyancy force and the Stoke friction force, as indicated in Equation (30).
ρ o i l · g 4 3 π R 3 + 6 π R · μ o i l · V T 0 .
We also started from V T , solving for the dynamic viscosity coefficient from Equation (30). This model was denoted Mod1.

5.3.2. Experimental Results for a Rate of Increase of 0.5 kV/s

Figure 22 represents 14 positions in the movement of a bubble for the experiment corresponding to Table 8. The time between two consecutive positions is equivalent to four frames. The time between frames is 0.0118 s. The image was previously scaled using the ImageJ program [45], which is useful for the subsequent image treatment and better definition of the bubbles. As can be seen in Figure 23, the contour of the bubble is improved, enabling a better measurement of its diameter.
The experimental results corresponding to a dielectric breakdown voltage with a rate of increase of 0.5 kV/s, for a mean dielectric breakdown voltage of 25.725 kV with a standard deviation of ±2.8303 kV and a mean room temperature of 17.2 °C are shown in Table 8.

5.3.3. Comparison of the Kinematic Viscosity with the Proposed Method and Manufacturer’s Data

Figure 24 represents the variation in kinematic viscosity with temperature using the manufacturer’s data as starting data [46]. For this purpose, the following kinematic viscosities were used: 10.5 cSt, 2.7 cSt and 950 cSt, for temperatures of 40 °C, 100 °C and −30 °C, respectively. With these data and starting from the two- and three-parameter models [47], the continuous black and red lines, respectively, were obtained, giving the viscosity values for temperatures from 0 to 100 °C. The discrete data point represented by a circle corresponds to the experimental measurement carried out in our laboratory, as shown in Table 8, by the standard method [39]. The discrete data point represented by a triangle corresponds to the result of applying Equation (30) to the evolution over time of a bubble recorded with the image sensor, according to Table 8. The mean error was less than 0.50%.

6. Conclusions

This paper describes an experimental study of the dielectric strength of transformer oil based on the IEC 60156 standard. The contribution made in this paper consists of characterizing the behavior of the oil before and after the electric arc breakdown, by combining a low-cost CMOS imaging sensor and a new matrix of electrical permittivity Mε associated with the dielectric oil, using the 3D cell method. The root mean square percentage error compared to the finite element method was less than 0.36%.
The IEC 60156 standard test indicates the effective value of the breakdown voltage. However, information on the distribution of Kelvin forces an instant before the dynamic behavior of the arc begins is lost, as is information on the gases that occur an instant after the moment of the electric arc breakdown in the oil. In this paper, after analyzing the images after rupture with a low-cost CMOS imaging sensor, the dynamic viscosity of the oil was indirectly estimated by measuring the rate of rise of the bubbles. These results were compared with a standard method (ASTM D445), and an error of less than 0.5% was obtained.

7. Annex

R2: coefficient of determination.
R 2 = 1 i = 1 N ( X i Y i ) 2 i = 1 N ( X i Y ¯ ) 2
Value: 1 R 2 1 , the nearer to 1, the better.
Advantage: indicates the proximity to the regression line. The perfect regression line has slope of 1.
Disadvantages: it does not always indicate a linear correlation between the data. If the sample is small, the data may, when enlarged, indicate a nonlinear correlation.
Consult: [36].
RMSPE: root mean square percentage error.
R M S P E = ( 1 N i = 1 N ( X i Y i ) 2 ) × 1 Y ¯
Value: 1 R M S P E 1 , the nearer to 0, the better.
Advantage: it is dimensionless and can be used to compare models.
Disadvantages: may underestimate the true measurement as it tries to reproduce the actual data.
Consult: [37].
MAEP: mean absolute percentage error.
M A E P = 1 N i = 1 N | X i Y i | × 1 Y ¯
Value: 1 M A E P 1 , the nearer to 0, the better.
Advantage: it is a dimensionless and robust measure of error.
Disadvantages: may underestimate the true measurement as it tries to reproduce the actual data.
Consult: [37].
PBIAS: percent bias.
P B I A S = i = 1 N ( X i Y i ) i = 1 N Y i × 100
Value: 1 P B I A S 1 , the nearer to 0, the better.
Advantage: mean of the deviations, or difference in bias errors, or simply systemic error. It detects aberrant errors or extreme values. If the value is equal to or greater than 1 it is indicative of the existence of these extreme values. Measures the percentage tendency of the simulated data to be larger or smaller than the reference data.
Disadvantages: may underestimate the true measurement as it tries to reproduce the actual data.
Consult: [38].

Author Contributions

J.M.M.-V. supported the theory background, implemented the low-cost imaging sensors, developed experiments, analyzed the data and wrote the paper. P.I.G.-D., S.G.-A. and J.V.R. supported the theory background, analyzed the data and wrote the paper and J.V.R. analyzed the material properties of the oil. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We wish to acknowledge the Electrical Engineering Department, the Department of Electronic Engineering and the Automatics and Process Engineering Department of the University of Las Palmas de Gran Canaria.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

SymbolNameUnit
ω Angular frequencyrad/s
DiBubble diameterm
σ b Bubble surface tensionN/m
RBubble radiusm
τeCharge-relaxation times
NData number-
XiData obtained-
ΔρDensity difference between oil and gaskg/m3
VDBVDielectric breakdown voltageV
I ˜ d Displacement electric currentA
lDomain characteristic lengthm
S ˜ i , S ˜ j Dual facesm2
Q f Electric charge in dual volumesC
σElectric conductivityS/m
MσElectric conductivity constitutive matrixS
I ˜ f Conductive electric currentA
E Electric field strengthV/C
ψ ˜ Electric flowC
ψ ˜ 0 Electric flow in vacuumC
M ϵ Electrical permittivity constitutive matrixF
ε0Electrical permittivity in vacuumF/m
Mε0Electrical permittivity in vacuum constitutive matrixF
εElectrical permittivity of the mediumF/m
UElectric potential differenceV
f FrequencyHz
gGravity accelerationm/s2
I c Incidence vector of relative cohomology between oil volume and electrode surface.-
C, C ˜ = C t Incidence matrix faces–edges in primal and dual mesh-
D ,   D ˜ Incidence matrix faces–volumes in primal and dual mesh-
GIncidence matrix edges–nodes of primal mesh-
jImaginary unit-
F e Kelvin polarization forceN
ΦMagnetic flowWb
F ˜ Magnetomotive forceA
ρ o i l Oil densitykg/m3
μ o i l Oil dynamic viscosityNs/m2
h o i l Oil heightm
ϵ o i l Oil relative permittivity-
ψ ˜ P Polarization electric flowC
M χ P Polarization electrical permittivity constitutive matrixF
Y i Reference data-
Y ¯ Reference data mean
εrRelative electrical permittivity-
φScalar electric potentialV
rSpherical particle radiusm
τSinusoidal signal periods
cSpeed of lightm/s
TTemperature°C
VTTerminal velocity of the bubble in the oilm/s
νTetrahedron volumem3
τemTime electromagnetic constants
τmTime magnetic constants
I ˜ t Total currentA
J Volumetric density currentA/m2
J d Volumetric density of conductive currentA/m2

References

  1. Islam, M.; Lee, G.; Hettiwatte, S.N. A review of condition monitoring techniques and diagnostic tests for lifetime estimation of power transformers. Electr. Eng. 2017, 100, 581–605. [Google Scholar] [CrossRef]
  2. Mahanta, D.K.; Laskar, S. Electrical insulating liquid: A review. J. Adv. Dielectr. 2017, 7, 1730001. [Google Scholar] [CrossRef] [Green Version]
  3. IEEE. Guide for Acceptance and Maintenance of Insulating Mineral Oil in Electrical Equipment; IEEE: Piscataway, NJ, USA, 2015. [Google Scholar]
  4. CIGRE. Life Management Techniques for Power Transformers; CIGRE: Paris, France, 2003; Volume 49. [Google Scholar]
  5. Kuffel, E.; Zaengl, W.; Kuffel, J. High Voltage Engineering Fundamentals; Elsevier: Amsterdam, The Netherlands, 2000; ISBN 0750636343. [Google Scholar]
  6. Radjenovic, B.; Klas, M.; Bojarov, A.; Matejcik, S. The breakdown mechanisms in electrical discharges: The role of the field emission effect in direct current discharges in microgaps. Acta Phys. Slovaca 2014, 63, 105–205. [Google Scholar] [CrossRef]
  7. Tonti, E. Why starting from diferential equations for computational physics? J. Comput. Phys. 2014, 257, 12601290. [Google Scholar] [CrossRef]
  8. Tonti, E. The Mathematical Structure of Classical and Relativistic Physics; Birkhäuser: Basel, Switzerland, 2013; p. 17. ISBN 9781461474210. [Google Scholar]
  9. Tonti, E. A Direct Discrete Formulation of Field Laws: The Cell Method. CMES Comput. Model. Eng. Sci. 2001, 2, 237–258. [Google Scholar]
  10. Marrone, M. Computational Aspects of the Cell Method in Electrodynamics. Prog. Electromagn. Res. 2001, 32, 317–356. [Google Scholar] [CrossRef] [Green Version]
  11. Codecasa, L.; Trevisan, F. Constitutive equations for discrete electromagnetic problems over polyhedral grids. J. Comput. Phys. 2007, 225, 1894–1918. [Google Scholar] [CrossRef]
  12. Bullo, M.; Dughiero, F.; Guarnieri, M.; Tittonel, E. Isotropic and Anisotropic Electrostatic Field Computation by Means of the Cell Method. IEEE Trans. Magn. 2004, 40, 1013–1016. [Google Scholar] [CrossRef]
  13. Bettini, P.; Trevisan, F. Electrostatic analysis for plane problems with finite formulation. IEEE Trans. Magn. 2003, 39, 1127–1130. [Google Scholar] [CrossRef]
  14. Bettini, P.; Specogna, R.; Trevisan, F. Electroquasistatic Analysis of the Gas Insulated Line for the ITER Neutral Beam Injector. IEEE Trans. Magn. 2009, 45, 996–999. [Google Scholar] [CrossRef]
  15. Monzon-Verona, J.M.; Santana-Martin, F.J.; Garcia-Alonso, S.; Montiel-Nelson, J.A. Electro-quasistatic analysis of an electro-static induction micromotor using the cell method. Sensors 2010, 10, 9102–9117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Specogna, R.; Trevisan, F. Discrete constitutive equations in A-χ geometric eddy-current formulation. IEEE Trans. Magn. 2005, 41, 1259–1263. [Google Scholar] [CrossRef]
  17. IEC 60156 Insulating Liquids Determination of the Breakdown Voltage at Power Frequency. Test Method. Available online: https://webstore.iec.ch/publication/28297IEC (accessed on 23 September 2021).
  18. Haus, H.A.; Melcher, J.R. Electromagnetic Fields and Energy; (Massachusetts Institute of Technology: MIT OpenCourseWare). License: Creative Commons Attribution—NonCommercial—Share Alike; Prentice-Hall: Englewood Cliffs, NJ, USA, 1989; ISBN 9780132490207. Available online: https://ocw.mit.edu/resources/res-6-001-electromagnetic-fields-and-energy-spring-2008/chapter-3/03.pdf (accessed on 23 September 2021).
  19. Baz-Rodríguez, S.; Aguilar-Corona, A.; Soria, A. Rising velocity for single bubbles in Pure Liquids. Rev. Mex. Ing. Quím. 2012, 11, 269–278. [Google Scholar]
  20. Pagnutti, M.; Ryan, R.E.; Cazenavette, G.; Gold, M.; Harlan, R.; Leggett, E.; Pagnutti, J. Laying the foundation to use Raspberry Pi 3 V2 camera module imagery for scientific and engineering purposes. J. Electron. Imaging 2017, 26, 013014. [Google Scholar] [CrossRef] [Green Version]
  21. Riba, J.-R.; Gómez-Pau, Á.; Moreno-Eguilaz, M. Experimental Study of Visual Corona under Aeronautic Pressure Conditions Using Low-Cost Imaging Sensors. Sensors 2020, 20, 411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Eichhorn, G.N.; Bowman, A.; Haigh, S.K.; Stanier, S. Low-cost digital image correlation and strain measurement for geotechnical applications. Strain 2020, 56, e12348. [Google Scholar] [CrossRef] [Green Version]
  23. Book: Electromagnetics and Applications (Staelin). Available online: https://phys.libretexts.org/Bookshelves/Electricity_and_Magnetism/Book%3A_Electromagnetics_and_Applications_(Staelin) (accessed on 23 September 2021).
  24. Schmid, S.; Hierold, C.; Boisen, A. Modeling the Kelvin polarization force actuation of micro- and nanomechanical systems. J. Appl. Phys. 2010, 107, 54510. [Google Scholar] [CrossRef] [Green Version]
  25. Sha, Y.; Zhou, Y.; Nie, D.; Wu, Z.; Deng, J. A study on electric conduction of transformer oil. IEEE Trans. Dielectr. Electr. Insul. 2014, 21, 1061–1069. [Google Scholar] [CrossRef]
  26. Li, H.; Zhong, L.; Yu, Q.; Mori, S.; Yamada, S. The resistivity of oil and oil-impregnated pressboard varies with temperature and electric field strength. IEEE Trans. Dielectr. Electr. Insul. 2014, 21, 1851–1856. [Google Scholar] [CrossRef]
  27. Taslak, E.; Arikan, O.; Kumru, C.F.; Kalenderli, O. Analyses of the insulating characteristics of mineral oil at operating conditions. Electr. Eng. 2017, 100, 321–331. [Google Scholar] [CrossRef]
  28. Woodson, H.H.; Melcher, J.R. Electromechanical Dynamics; John Wiley and Sons, Inc.: New York, NY, USA, 1968; Reprint Krieger Publishing Company: Malabar, FL, USA, 1990, 1985, ISBN 9780894644597 (v.1), 9780898748475 (v. 2), 9780898748482 (v. 3).18. [Google Scholar]
  29. Craig, A.; Coburna, A.M.; Smitha, B.; Gordon, S.; Logieaand, P.K. Radiometric and spectral comparison of inexpensive camera systems used for remote sensing. Int. J. Remote Sens. 2018, 39, 4869–4890. [Google Scholar] [CrossRef]
  30. Home—OpenCV. Available online: https://opencv.org (accessed on 23 September 2021).
  31. NumPy. Available online: https://numpy.org/ (accessed on 23 September 2021).
  32. Picamara. Available online: https://picamera.readthedocs.io/en/release-1.13/ (accessed on 23 September 2021).
  33. Portable, Extensible Toolkit for Scientic Computation. Available online: http://www.mcs.anl.gov/petsc (accessed on 23 September 2021).
  34. Dular, P.; Geuzaine, C. GetDP Reference Manual: The Documentation for GetDP, a General Environment for the Treatment of Discrete Problems. Available online: http://getdp.info (accessed on 23 September 2021).
  35. Geuzaine, G.; Remacle, J.-F. Gmsh: A three-dimensional finite element mesh generator with built-in pre- and post-processing facilities. Int. J. Numer. Methods Eng. 2009, 79, 13091331. [Google Scholar] [CrossRef]
  36. Martínez, E.R. Errores frecuentes en la interpretación del coeficiente de determinación lineal. Anu. Jurid. Econ. Escur. 2005, 38, 317–331. [Google Scholar]
  37. Hyndman, R.; Koehler, A.B. Another look at measures of forecast accuracy. Int. J. Forecast. 2006, 22, 679–688. [Google Scholar] [CrossRef] [Green Version]
  38. Sanabria, J.; García, J.; Lhomme, J.P. Calibración y Validación de Modelos de Pronóstico de Heladas en el Valle del Mantaro. Master’s Thesis, Universidad Nacional Agraria La Molina, La Molina, Peru, 2006. [Google Scholar]
  39. ASTM D445. Standard Test Method for Kinematic Viscosity of Transparent and Opaque Liquids (and Calculation of Dynamic Viscosity); ASTM Volume 05.01 Petroleum Products, Liquid Fuels, and Lubricants (II): C1234—D4176; ASTM: West Conshohocken, PA, USA, 2019. [Google Scholar]
  40. ASTM. ASTM D4052. Standard Test Method for Density, Relative Density, and API Gravity of Liquids by Digital Density Meter; ASTM Volume 05.01 Petroleum Products, Liquid Fuels, and Lubricants (II): C1234—D4176; ASTM: West Conshohocken, PA, USA, 2018. [Google Scholar]
  41. ASTM. ASTM D92. Standard Test Method for Flash and Fire Points by Cleveland Open Cup Tester; ASTM Volume 05.01 Petroleum Products, Liquid Fuels, and Lubricants (II): C1234—D4176; ASTM: West Conshohocken, PA, USA, 2018. [Google Scholar]
  42. ASTM. ASTM D1533: Standard Test Method for Water in Insulating Liquids by Coulometric Karl Fischer Titration; ASTM Volume 10.03 Electrical Insulating Liquids and Gases; Electrical Protective Equipment; ASTM: West Conshohocken, PA, USA, 2020. [Google Scholar]
  43. ASTM. ASTM D1500—12 Standard Test Method for ASTM Color of Petroleum Products (ASTM Color Scale); ASTM Volume 05.01 Petroleum Products, Liquid Fuels, and Lubricants (II): C1234—D4176; ASTM: West Conshohocken, PA, USA, 2017. [Google Scholar]
  44. Norma UNE 60296:2012. Fluidos para Aplicaciones Electrotécnicas. Aceites Minerales Aislantes Nuevos para Transformadores y Aparamenta de Conexión; AENOR: Madrid, Spain, 2012. [Google Scholar]
  45. ImajeJ. Available online: https://imagej.net/Welcome (accessed on 23 September 2021).
  46. Aceites Dieléctricos. Available online: https://www.repsol.com/imagenes/global/es/af-catalogo-aceites-dielectricos_tcm13-48939.pdf (accessed on 23 September 2021).
  47. Stanciu, I. A new viscosity-temperature relationship for mineral oil SAE 10W. Analele Univ. Ovidius Constanta Ser. Chim. 2012, 23, 27–30. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Reference tetrahedron for equations programming. The elements of the primal mesh are represented in black and dual planes in green. A current associated with the dual plane S ˜ 3 is also shown.
Figure 1. Reference tetrahedron for equations programming. The elements of the primal mesh are represented in black and dual planes in green. A current associated with the dual plane S ˜ 3 is also shown.
Sensors 21 07380 g001
Figure 2. (a) Dimensions and coordinates of VDE-type electrodes. (b) Surface mesh, with 65,846 triangles, for the electrodes used in the tests.
Figure 2. (a) Dimensions and coordinates of VDE-type electrodes. (b) Surface mesh, with 65,846 triangles, for the electrodes used in the tests.
Sensors 21 07380 g002
Figure 3. The 51,679 edges corresponding to the relative cohomology between the oil volume and the electrode surface, to obtain the incidence vector Ic from Equation (15).
Figure 3. The 51,679 edges corresponding to the relative cohomology between the oil volume and the electrode surface, to obtain the incidence vector Ic from Equation (15).
Sensors 21 07380 g003
Figure 4. Cameras A, B and C connected with Raspberry Pi and router for the collection of images in a central computer.
Figure 4. Cameras A, B and C connected with Raspberry Pi and router for the collection of images in a central computer.
Sensors 21 07380 g004
Figure 5. Front view of central camera B.
Figure 5. Front view of central camera B.
Sensors 21 07380 g005
Figure 6. 40-pin GPIO connectors to synchronize the triggering of the three cameras.
Figure 6. 40-pin GPIO connectors to synchronize the triggering of the three cameras.
Sensors 21 07380 g006
Figure 7. Electrical scalar potential φ obtained in CM when solving system of equations shown in Equation (16), meshed with 2,790,589 tetrahedra and cuts A and B for comparison.
Figure 7. Electrical scalar potential φ obtained in CM when solving system of equations shown in Equation (16), meshed with 2,790,589 tetrahedra and cuts A and B for comparison.
Sensors 21 07380 g007
Figure 8. Module of the volumetric density of the conductive current J f corresponding to Equation (1).
Figure 8. Module of the volumetric density of the conductive current J f corresponding to Equation (1).
Sensors 21 07380 g008
Figure 9. Module of the volumetric density of the displacement current J d corresponding to Equation (4).
Figure 9. Module of the volumetric density of the displacement current J d corresponding to Equation (4).
Sensors 21 07380 g009
Figure 10. Kelvin force density divided by 1 2 ( ε ε o i l ) cut by the plane 1 X + 0 Y +0 Z − 0.0489 = 0.
Figure 10. Kelvin force density divided by 1 2 ( ε ε o i l ) cut by the plane 1 X + 0 Y +0 Z − 0.0489 = 0.
Sensors 21 07380 g010
Figure 11. Kelvin force density divided by 1 2 ( ε ε o i l ) cut by the plane 0 X + 1 Y +0 Z − 0.035 = 0.
Figure 11. Kelvin force density divided by 1 2 ( ε ε o i l ) cut by the plane 0 X + 1 Y +0 Z − 0.035 = 0.
Sensors 21 07380 g011
Figure 12. Capture of the images of the electric arc by cameras (AC).
Figure 12. Capture of the images of the electric arc by cameras (AC).
Sensors 21 07380 g012
Figure 13. Voltage in cut A (steeper curves) and cut B (curves with lower slope) for experiment E1. The dielectric breakdown voltage is 19,488 V.
Figure 13. Voltage in cut A (steeper curves) and cut B (curves with lower slope) for experiment E1. The dielectric breakdown voltage is 19,488 V.
Sensors 21 07380 g013
Figure 14. Electric field in cut A (upper curves) and cut B (lower curves) for experiment E1.
Figure 14. Electric field in cut A (upper curves) and cut B (lower curves) for experiment E1.
Sensors 21 07380 g014
Figure 15. Voltage in cut A (steeper curves) and cut B (curves with lower slope) for experiment E2. The dielectric breakdown voltage is 19,488 V.
Figure 15. Voltage in cut A (steeper curves) and cut B (curves with lower slope) for experiment E2. The dielectric breakdown voltage is 19,488 V.
Sensors 21 07380 g015
Figure 16. Electric field in cut A (upper curves) and cut B (lower curves) for experiment E2.
Figure 16. Electric field in cut A (upper curves) and cut B (lower curves) for experiment E2.
Sensors 21 07380 g016
Figure 17. Voltage in cut A (steeper curves) and cut B (curves with lower slope) for experiment E3. The dielectric breakdown voltage is 19,488 V.
Figure 17. Voltage in cut A (steeper curves) and cut B (curves with lower slope) for experiment E3. The dielectric breakdown voltage is 19,488 V.
Sensors 21 07380 g017
Figure 18. Electric field in cut A (upper curves) and cut B (lower curves) for experiment E3.
Figure 18. Electric field in cut A (upper curves) and cut B (lower curves) for experiment E3.
Sensors 21 07380 g018
Figure 19. Error histograms: (a) comparison C1; (b) comparison C4; (c) comparison C9; (d) comparison C12.
Figure 19. Error histograms: (a) comparison C1; (b) comparison C4; (c) comparison C9; (d) comparison C12.
Sensors 21 07380 g019aSensors 21 07380 g019b
Figure 20. Oil testing device with oil sample.
Figure 20. Oil testing device with oil sample.
Sensors 21 07380 g020
Figure 21. Electric arc at the instant of the breakdown of the electrical strength, with a value of the dielectric breakdown voltage of 33.9 kV.
Figure 21. Electric arc at the instant of the breakdown of the electrical strength, with a value of the dielectric breakdown voltage of 33.9 kV.
Sensors 21 07380 g021
Figure 22. Position of the bubble at 14 instants of time for the experiment in Table 8.
Figure 22. Position of the bubble at 14 instants of time for the experiment in Table 8.
Sensors 21 07380 g022
Figure 23. Improvement in the definition of the bubble using the ImageJ program for measuring the diameter.
Figure 23. Improvement in the definition of the bubble using the ImageJ program for measuring the diameter.
Sensors 21 07380 g023
Figure 24. Kinematic viscosity versus temperature in dielectric oil.
Figure 24. Kinematic viscosity versus temperature in dielectric oil.
Sensors 21 07380 g024
Table 1. Mesh characteristics and physical properties considered in the numerical experiments.
Table 1. Mesh characteristics and physical properties considered in the numerical experiments.
ExperimentNodesEdgesFacesVolumesConductivity
[S/m]
Relative
Permittivity
BDV
[V]
E1117564679687439410103.819,488
E2588934,73754,14225,29310103.819,488
E3487,4353,362,3855,665,5402,790,58910103.819,488
Table 2. Comparisons of developed numerical experiments.
Table 2. Comparisons of developed numerical experiments.
C1Numerical experiment E1-Cut A. Electric field.
C2Numerical experiment E1-Cut B. Electric field.
C3Numerical experiment E1-Cut A. Voltage.
C4Numerical experiment E1-Cut B. Voltage.
C5Numerical experiment E2-Cut A. Electric field.
C6Numerical experiment E2-Cut B. Electric field.
C7Numerical experiment E2-Cut A. Voltage.
C8Numerical experiment E2-Cut B. Voltage.
C9Numerical experiment E3-Cut A. Electric field.
C10Numerical experiment E3-Cut B. Electric field.
C11Numerical experiment E3-Cut A. Voltage.
C12Numerical experiment E3-Cut B. Voltage.
Table 3. Metrics of the proposed comparisons.
Table 3. Metrics of the proposed comparisons.
ComparisonC1C2C3C4C5C6References
R2 [0, +1] Optimum: +10.57720.00001.00000.99790.00000.1083[36]
RMSPE [−1, +1] Optimum: 00.36240.11270.00450.05120.01520.1302[37]
MAEP [−1, +1] Optimum: 00.24090.11190.00400.04050.01290.0959[37]
PBIAS [−1, +1] Optimum: 00.1447−0.12600.00000.02160.00010.0254[38]
Table 4. Metrics of the proposed comparisons.
Table 4. Metrics of the proposed comparisons.
ComparisonC7C8C9C10C11C12References
R2 [0, +1] Optimum: +11.00000.99980.99971.00001.00001.0000[36]
RMSPE [−1, +1] Optimum: 00.00500.01010.00030.00010.00000.0000[37]
MAEP [−1, +1] Optimum: 00.00440.00860.00020.00000.00000.0231[37]
PBIAS [−1, +1] Optimum: 0−0.00190.00540.00000.00000.00000.0000[38]
Table 5. Technical characteristics.
Table 5. Technical characteristics.
OT-40Value
Supply voltage230 V
Test voltage0–41 kV
Consumption100 VA/800 VA
Frequency50 Hz
Measurement tolerance±2%
Response time on disconnection<20 ms
Test standardsUNE EN 60156:1997
Work temperature+15/+25 °C
Exterior size385 × 300 × 400 mm
Weight38 kg
Table 6. Results of the oil sample tests.
Table 6. Results of the oil sample tests.
ParametersUnitMethodValue
Kinematic viscosity at 40 °CcStASTM D4458.75
Density at 15 °Cg/mLASTM D40520.858
Flashpoint°CASTM D92145.0
Water contentppmASTM D153320.40
Color-ASTM D15001.0
Total aciditymg KOH/gIEC 61125C0.06
Table 7. Breakdown voltage. Ramp of 0.5 kV/s.
Table 7. Breakdown voltage. Ramp of 0.5 kV/s.
Sample NumberfpsT [°C]VDBV [kV]Agitator [Y/N]
RCL
187.084.786.217.124.4Y
287.084.786.217.124.4Y
387.084.786.317.123.5Y
486.984.786.217.223.3Y
587.085.086.217.323.1Y
687.084.786.117.329.1Y
787.084.786.117.330.0Y
887.184.786.117.328.0Y
Table 8. Determination of the kinematic viscosity of the oil using gas bubbles.
Table 8. Determination of the kinematic viscosity of the oil using gas bubbles.
Experiment No.Bubble
Diameter [mm]
Bubble Mean Velocity [m/s]Voltage
[kV]
Viscosity
(cSt) Mod1
Viscosity
(cSt) Mod2
11.330.02824.433.654.7
21.340.03024.432.149.8
31.280.02623.534.660.4
41.330.02823.133.254.3
51.340.03029.131.949.3
61.260.03230.026.939.0
Mean 25.75Mean 32.05Mean 51.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Monzón-Verona, J.M.; González-Domínguez, P.I.; García-Alonso, S.; Vaswani Reboso, J. Characterization of Dielectric Oil with a Low-Cost CMOS Imaging Sensor and a New Electric Permittivity Matrix Using the 3D Cell Method. Sensors 2021, 21, 7380. https://doi.org/10.3390/s21217380

AMA Style

Monzón-Verona JM, González-Domínguez PI, García-Alonso S, Vaswani Reboso J. Characterization of Dielectric Oil with a Low-Cost CMOS Imaging Sensor and a New Electric Permittivity Matrix Using the 3D Cell Method. Sensors. 2021; 21(21):7380. https://doi.org/10.3390/s21217380

Chicago/Turabian Style

Monzón-Verona, José Miguel, Pablo Ignacio González-Domínguez, Santiago García-Alonso, and Jenifer Vaswani Reboso. 2021. "Characterization of Dielectric Oil with a Low-Cost CMOS Imaging Sensor and a New Electric Permittivity Matrix Using the 3D Cell Method" Sensors 21, no. 21: 7380. https://doi.org/10.3390/s21217380

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop