UDC 621.3.011.5
doi: 10.20998/2074-272X.2020.2.09
O.O. Palchykov
DETERMINATION OF THE EFFECTIVE PERMITTIVITY OF A HETEROGENEOUS MATERIAL
Purpose. To develop a two-dimensional numerical-field model for determining the effective permittivity of a multicomponent material represented by a system of homogeneous volumes with known physical characteristics. Methodology. The model is based on the solution by the finite element method of an electrostatic problem with the subsequent determination of the energy contained in the volume under consideration. Than we have compared this result with the energy of a flat capacitor with a rectangular cross-section of the plates and determined the effective permittivity of test material. We also have used Rayleigh, Odelevsky and Lichtenecker models and the model with a perpendicular arrangement of layers relative to the main electric flux. Results. Based on the developed field model, the effective permittivities for dry, wet and transformer oil-soaked insulating papers of various grades, including taking into account ash, are determined. We have proved that a macroscopically homogeneous multicomponent material is well approximated by uniformly spaced cylindrical volumes with a substance of different nature in a matrix of another substance. We have showed a significant error of the layer model and the Rayleigh model relative to the proposed model. We have showed the equivalence of models with the location of inclusions in the nodes of a rectangular and parallelogram mesh. Originality. For the first time we have proposed wet paper models with an asymmetric arrangement of a cylindrical volume of water with a circular and segment cross-section in a cylindrical pore. For the first time we have proposed models of insulating paper with evenly spaced cylindrical inclusions of different volumes. Practical value. The proposed model allows to calculate the effective permittivity of an inhomogeneous material with a given accuracy without restricting the shape of the components. Based on the proposed field model, it is possible to determine the Lichtenecker index, which allows to calculate the effective permittivity for any ratio of the volumes of the components of a heterogeneous material. References 10, tables 3, figures 4.
Key words: effective permittivity, electrostatic field, energy, finite element method, cylindrical volumes.
ОбГрунтовано застосування методу, заснованого на чисельному розрахунку електростатичного поля для визначення ефективног дiелектричноí проникност гетерогенного матерiалу, який замшяеться системою однорiдних об'ежв з вiдомими фiзичними характеристиками. Для сухого, вологого i просоченого трансформаторним маслом iзоляцiйних паперiв рiзних марок визначеш дiелектричнi проникност^ Результати розрахунку за запропонованими моделями зктавлялися з результатами, отриманими на основi шаруватоí модет, моделей Релея i Оделевського. Запропоноват апроксимацiйнi залежностi на основi узагальненого виразу Л^хтенеккера для визначення дiелектричноí проникностi iзоляцiйних паперiв. Бiбл. 10, табл. 3, рис. 4.
Ключовi слова: ефективна дiелектрична проникшсть, електростатичне поле, енерпя, метод скшченних елеменив, цилвдричш об'еми.
Обосновано применение метода, основанного на численном расчете электростатического поля для определения эффективной диэлектрической проницаемости гетерогенного материала, который представляется системой однородных объемов с известными физическими характеристиками. Для сухой, влажной и пропитанной трансформаторным маслом изоляционных бумаг различных марок определены диэлектрические проницаемости. Результаты расчета по предложенным моделям сопоставлялись с результатами, полученными на основе слоистой модели, моделей Рэлея и Оделевского. Предложены аппроксимационные зависимости на основе обобщенного выражения Лихтенеккера для определения диэлектрической проницаемости изоляционных бумаг. Библ. 10, табл. 3, рис. 4.
Ключевые слова: эффективная диэлектрическая проницаемость, электростатическое поле, энергия, метод конечных элементов, цилиндрические объемы.
Introduction. In the manufacture of electrical machines, apparatus and power capacitors, various insulating materials are used, such as multicomponent (varnished fabrics and varnished paper, electrical insulating tapes, laminated plastics, electrical ceramics), and almost uniform ones, such as mica and cellulose papers. The latter, however, contain a large number of cavities that can be filled with atmospheric air and moisture or an impregnating dielectric. Therefore, the permittivity, conductivity and breakdown voltage of insulating materials is a complex function of the presence and placement of certain components in the volume of the dielectric. In this paper, we consider the effect of the composition of an insulating material on its effective permittivity.
Literature review. To calculate the effective permittivity (EP) £eq of two-component materials, the
formulas of Wiener, Rayleigh, Maxwell, generalized Lichtenecker formula for flat, cylindrical, spherical and granular inclusions [1, 2], as well as the Odelevsky formula for statistical mixtures can be used:
teq + k YV = (t + k )~1Vi + ( + к
(1)
(ffeq - t )/(ffeq +1 ) = (t - t )(t +1) 1V2 VL ^ (2)
tq = t[1 + (3V2(t -t)/V£)) + 2t -V(t)))] ; (3)
m.
П '
(4)
j=1
teq = A + J A2 +et2/2
(5)
where s1, s2 are the permittivities of the components of the material; V2, V1, V2 are the total volume and volumes of
© O.O. Palchykov
the components of the material, respectively; k is the Wiener coefficient, 0 < k < n is the number of components of the material; m is the Lichtenecker index, -1 < m < 1; A is the Odelevsky parameter
A = [(3 V / -1) + (3V2 / - 1)S2 ]/4 .
The use of formulas (1), (4) is difficult due to the unknown coefficients k and m. Formulas (2), (3) and (5) are derived on the basis of certain assumptions, the violation of which can lead to a significant error.
To determine the EP of impregnated paper, in [3] it is proposed to use a model of series-connected layers, which is a rather rough approximation.
In [4-6], using the double and triple-periodic Rayleigh model, effective material parameters are calculated based on the solution of the static field problem by the method of summing multipole interactions. The field models in these works for reducing infinite sums of interactions to their finite sum are limited only by the nearest neighbors with respect to the chosen inclusion.
Calculation of the EP can be based on the polarization mechanism [7]. This approach is mainly applied only for pure substances with consideration of processes at the ionic and molecular levels. In [8], based on field models, the electrophysical properties of a cement-based composite were explained.
Therefore, a sufficiently accurate determination of the averaged characteristics of the material is based on the solution of the corresponding field model of the system of solid substances located in the volume of the material in accordance with the technology of its manufacture.
The goal of the work is the development of a two-dimensional numerical-field model for determining the EP of a multicomponent material represented by a system of homogeneous volumes with known physical characteristics.
Object of study. A rectangular parallelepiped made of capacitor paper with the following dimensions: length l = 100 ^m; width b = 200 ^m; height h = 14 ^m for papers KOH 0,8; KOH 2; MKOH 1 or height h = 12 ^m for paper CKOH 3,5. Humidity, composition and density of the selected paper grades correspond to [9].
Mathematical model. The main physical properties of the materials needed to build mathematical models for determining the EP of a particular paper are given in Table 1. The temperature of the insulating paper is taken 60 °C.
Table 1
Physical properties of dielectric materials
Name Density, kg/m3 e20 a, i/°C
Cellulose i530 6.5 0.5-i0-3
Air i,06 i.00058 -i.5i0-6
Water 983 80.2 -4.02-i0-3
Transformer oil 856 2.25 -0.5-i0-3
The permittivity of the material at an arbitrary temperature t is determined as
S = S20 [1 + a(t - 20)]. (6)
The veracity of expression (6) is confirmed by the graphical dependencies shown in Fig. 1 for dry air and water, recommendations in [3] for cellulose and the linear
dependence of the permittivity of transformer oil on temperature within 20...90 °C. The assumption of the absence of thermal aging of the material also applies.
20 30 40 50 60t°C 20 40
a
80 100/,°C
b
Fig. 1. Dependence of the permittivity of dry air (a) and water (b) on temperature
Real insulating paper is presented as cellulose with cylindrical pores uniformly spaced in its volume. Cross sections of pores are in planes perpendicular to the length of the volume under consideration. Depending on the model formulation, air, air and water, transformer oil may be in the pores.
The volumes of cellulose and pores are found by solving the system of equations:
Vp Pp ^ - y) = Vc Pc + Va Pa;
■ Vw Pw = Vp Pp y;
Vp = Vc + Va + Vw,
where Vp, Vc, Va and Vw are the volumes of paper, cellulose, air and water, respectively; Pp, pc, Pa and pw are the densities of paper, cellulose, air and water, respectively; y is the relative mass content of water in the insulating paper [9].
The electric field in cross section relative to the length of the considered volume of insulating paper is described by the following differential equations [10]:
eV 2p = -p; E = -Vp; D = e E,
(7)
where q> is the scalar electrostatic potential; p is the volume charge density of the domain; E is the electric field vector; D is the electric induction vector.
At the interface, the following conditions are satisfied [10]:
in12 x[v^1 -V^2] = 0;
ni2 -{siV(f>i -eflVT.)); (Pi =f2;
(8)
ei
P 3(2
e2_
■ = <y,
3n12 3n12
where n12 is the normal from the first to the second medium; a is the surface charge at the interface between two media.
Equation (7) in a piecewise homogeneous dielectric medium is reduced in each homogeneous region to the Laplace equation, and the permittivities of the regions is included in the solution of the problem only by means of conditions (8).
The second equation of condition (8) at the interfaces between two dielectrics, as well as a dielectric and a conductor, respectively, takes the form:
an
cV2 n a^i — = 0; £1-— = CT-
12
an
12
an
12
On the upper and lower planes perpendicular to the height of the volume under consideration, Dirichlet conditions respectively ^d1 = 1.9 V; ^d2 = -1.9 V are set. The solution domain is limited to a cylindrical surface with a radius of 200 ^m, a length of 100 ^m, and zero potential on its surface.
After the numerical calculation of the electric field the energy enclosed in the volume of insulating paper is found
W = J-2E Ddv = -2 Je(v^)2dv,
V 2 2 V
where V is the volume over which the integration is performed; v is the elementary volume.
By comparing the calculated energy with the energy of a flat capacitor, we find the EP
2Wh
s =-
lb (1 -q>A2 )2
The numerical calculation of the electrostatic field by the finite element method is implemented in the FEMM code.
The main assumptions of the model: dielectric materials do not have conductivity; there are no charges at their interfaces; space charges in dielectric materials are
absent, and their volumes do not change under the influence of an electrostatic field; ideal conductors with infinitesimal height are sources of an electrostatic field.
The significant width of the investigated dielectric is explained by the desire to reduce the edge effect. The model was tested by comparing capacitor capacities with sizes 100x200x14 ^m and the same thickness of cellulose and air obtained on the basis of the equations of the electrostatic field by the finite element method and on the basis of the analytical expression without taking into account the edge effect. The error in this case was 0.23 %.
The results of the study. To calculate the EP of dry and transformer oil-soaked insulating papers, we used models of uniformly distributed cylindrical volumes filled with air and transformer oil, respectively, in cellulose matrices.
To calculate the EP of wet paper, a model of uniformly distributed cylindrical volumes filled with air and water was proposed, and in this total volume the volume of water was represented as cylindrical with a circular cross section in a cellulose matrix. When determining the volume content of a particular component, the author was based on the structure of wet paper (that is, paper with normal water content [9]), therefore, the pore volume (the sum of the volume of air and water for wet paper) in the models of dry and impregnated paper did not change. The calculated EPs for the indicated models, as well as for the layered model, Rayleigh and Odelevsky models are given in Table 2.
Table 2
Simulation results
e
Paper grade Proposed model Layer model Rayleigh model Odelevsky model Approximation dependence
Number of pores Number of model nodes ^eq Number of model nodes ^eq ^eq ^eq eeq
dry paper
Paper KOH 0,8 248 662632 3.052 4607 1.682 4.500 2.860 3.052
Paper KOH 2 248 702562 4.370 4749 2.543 5.300 4.456 4.507
Paper MKOH 1 248 693639 3.674 4920 2.048 4.895 3.667 3.783
Paper CKOH 3,5 248 750148 4.966 4872 3.150 5.665 5.114 5.123
wet paper
Paper KOH 0,8 248 4462060 3.403 5566 1.912 - - 2.740
Paper KOH 2 248 4877348 5.185 5652 3.503 - - 5.185
Paper MKOH 1 248 4744585 4.130 5509 2.448 - - 3.634
Paper CKOH 3,5 248 3782978 6.262 6711 5.095 - - 6.808
impregnated paper
Paper KOH 0,8 248 662632 3.949 4607 3.233 4.789 3.914 3.949
Paper KOH 2 248 702562 5.000 4749 4.212 5.525 5.068 5.06
Paper MKOH 1 248 693639 4.453 4920 3.686 5.155 4.503 4.516
Paper CKOH 3,5 248 750148 5.448 4872 4.755 5.835 5.540 5.525
The decrease in the EP values according to the the need for the entire electric flux to pass through the layered model relative to the proposed one is explained by region with low dielectric constant. The Rayleigh model
also cannot claim to be a good description of the change in the EP when the volumes and properties of the components of the material change and gives overestimated values. The Odelevsky model determines the EP with an accuracy of -2.98...6.3 % and -1.7...0.9 % for dry and impregnated paper, respectively. On the basis of field models and the generalized Lichtenecker formula, approximation dependences are proposed for calculating the EP of insulating papers. The approximation dependences indexes for dry, wet and transformer oil-soaked papers amounted to 0.495; -0.283 and 0.391, respectively. The EPs calculated by this method are given in Table 2.
The effect of the number of inclusions on the EP value of KOH 2 dry paper is shown in Fig. 2, from which it can be seen that with an increase in the number of uniformly distributed volumes of inclusions, the accuracy of the model increases until it begins to be limited by the accuracy of displaying the boundaries of these volumes.
Consequently, the proposed models for the specified in Table 2 numbers of pores accurately enough allow to determine the EP.
In addition to the above wet paper models, models with a cylindrical volume of water and a sector section, as well as with a cylindrical volume of water, in the cavity of which a cylindrical air volume is coaxially located are considered.
£eq 4.4
3.9 -I----
100 150 200 250 300 N
Fig. 2. Dependence of the permittivity on the number of pores of the model of dry paper grade KOH 2
The maximum discrepancy between the EPs calculated for models with a circular and sector section of the water volume was 2.9 %, therefore, these models are equivalent. Examples of the distribution of the electrostatic field in these models are shown in Fig. 3, 4.
Fig. 3. Distribution of the electrostatic field in the model of wet paper brand KOH 2 with a circular cross section of the volume
of water
Fig. 4. Distribution of the electrostatic field in the model of wet paper brand KOH 2 with a sector cross-section of the volume
of water
As can be seen from Fig. 3, 4, with almost identical EPs, the electric field strength modulus in the air in the first model is greater than in the second one. The EP values determined for models with a hollow cylindrical volume of water are 1.36...1.8 times greater than for a model with a cylindrical volume of circular cross section, which is explained by the formation of a path with an increased permittivity. This model was discarded due to
the presence of a local extremum of the EP value when changing paper grades.
For the KOH 2 dry paper model and pore number 100, the effect of the location of pores at the nodes of a rectangular and parallelogram mesh on the final result was studied. The maximum discrepancy was 0.8 %.
It is known that CKOH paper is characterized by a significant, up to 1 %, ash content. In the manufacture of
paper calcium sulfate can be considered as ash. The calculated values of the EP of CKOH insulating paper, taking into account ash, are given in Table 3. The error of not accounting of ash was 0.8...1.3 %.
Table 3
CKOH 3,5 paper simulation results with ash inclusions 1 %
Paper type Number of pores Number of inclusions of ash Number of model nodes ^eq
Dry 248 49 843735 4.925
Wet 248 49 3997067 6.182
Impregnated 248 49 843735 5.415
Conclusions.
1. The principal possibility of describing the permittivity of a mixture of components on the basis of field models with the known permittivity of each component, their volumetric content and the features of their location is shown. Field models have been developed for determining the EP of dry, wet, and transformer oil-soaked insulating papers.
2. The adequacy of the developed models is confirmed by the following: with an increase in the number of cylindrical volumes inside the matrix, the EP asymptotically tends to a certain finite value.
3. The layered model and the Rayleigh model give, respectively, underestimated and overestimated values of the EP relative to the proposed field model. These alternative models suitably describe only a mixture of components with close permittivities.
4. Approximate dependencies of the behavior of the EP of component mixtures with a change in their volume content, constructed on the basis of the generalized Lichtenecker model and the proposed field models, have comparable accuracy with the Odelevsky model for two-component materials and a 2-fold reduced error relative to the layered model for three-component materials.
5. Models with the arrangement of cylindrical volumes in the nodes of a rectangular and parallelogram mesh give the same result when calculating the EP of a macroscopically homogeneous heterogeneous material.
REFERENCES
1. Raju G. G. Dielectrics in electric fields. New York, Marcel Dekker, 2003. 578 p.
2. Moulson A.J., Herbert J.M. Electroceramics: materials, properties, applications. Chichester, John Wiley & Sons Ltd,
2006. 557 p.
3. Shcheglov N.V. Sovremennye vidy izoliatsii. Izoliatsiia silovykh kondensatorov [Modern types of isolation. Power capacitor isolation]. Novosibirsk, NSTU Publ., 2016. 116 p. (Rus).
4. Balagurov B.Ya., Kashin V.A. Electrophysical properties of the Rayleigh three-dimensional model. Technical Physics, 2007, vol. 52, no. 2, pp. 216-225. doi: 10.1134/s1063784207020119.
5. Emets Yu.P. Electrical characteristics of three-component dielectric composites with close-packed inclusions. Journal of Applied Mechanics and Technical Physics, 2001, vol. 42, iss. 4, pp. 704-713. doi: 10.1023/A:1019272201997.
6. Tolmachev S.T., Yukhymovych D.L. The problem of reducing the hollow circular cylinders with arbitrary grating of periods for regular system. Technical Electrodynamics, 2011, no. 2, pp. 11-17. (Rus).
7. He T. Effects of impurities on silicon dioxide dielectric properties. Advanced Materials Research, 2014, vol. 1022, pp. 56-59. doi: 10.4028/www.scientific.net/AMR.1022.56.
8. Chen R.J., Zhang Y., Wang B. Numerical simulation study on the cement-based absorbing material. Advanced Materials Research, 2013, vol. 853, pp. 169-173. doi: 10.4028/www.scientific.net/AMR.853.169.
9. DSTU 3467-96. Papir kondensatornyj. Zagal'ni tehnichni umovy. [DSTU 3467-96. State standard of Ukraine. Capacitor paper. General technical conditions]. Kyiv, Derzhstandart of Ukraine Publ., 1996. 34 p. (Ukr).
10. Stratton J.A. Electromagnetic theory. Hoboken, IEEE Press,
2007. 630 p.
Received 17.01.2020
O.O. Palchykov, Candidate of Technical Science, Associate Professor,
Admiral Makarov National University of Shipbuilding, 9, Heroyiv Ukraine Ave, Mykolaiv, Ukraine, 54025, Ten/phone +380 95 8934509; e-mail: ole2012hulk@gmail.com
How to cite this article:
Palchykov O.O. Determination of the effective permittivity of a heterogeneous material. Electrical engineering & electromechanics, 2020, no.2, pp. 59-63. doi: 10.20998/2074-272X.2020.2.09.