Mathematical Modeling on AC Pollution Flashover Performance of Glass and Composite Insulator
Mathematical Modeling on AC Pollution Flashover Performance of Glass and Composite Insulator
Journal of Electrical Engineering and Technology. 2015. Jul, 10(4): 1796-1803 This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
• Received : July 10, 2014
• Accepted : March 31, 2015
• Published : July 01, 2015 PDF e-PUB PubReader PPT Export by style
Article
Author
Metrics
Cited by
TagCloud
N.B., Prakash
Dept. of Electrical and Electronic Engineering, National Engineering College, Kovilpatti, India. (nbprakash@gmail.com)
M., Parvathavarthini
Dept. of Electrical and Electronic Engineering, National Engineering College, Kovilpatti, India. (nbprakash@gmail.com)
Corresponding Author: Dept. of Electrical and Electronic Engineering, National Engineering College, Kovilpatti, India. (srmadavan@gmail.com)

Abstract
While considering the current scenario, in this world power demand goes on increases day by day. Frequent power outages occur in high voltage transmission line due to the deprived performance of polluted insulators; this affects overall operation of power system and may indirectly impinge on the growth of production sector. Many researchers are keenly taking efforts to provide highly reliable and stable power to neediest. In this paper, A.C pollution flashover performance of disc type glass insulator and composite long rod insulators investigation under various artificial pollutions by varying Equivalent Salt Density Deposition (ESDD) levels. Here, we use different types of pollution methods like binding method, dipping method and spraying methods with different types of pollutants concentration. Based on dimensional analysis, four different Mathematical models have been developed to predict the A.C pollution Flashover Voltage (FOV) of insulators. Both the experimental and mathematically modeled results are compared; it’s observed that mathematical model 3 yields better results.
Keywords
1. Introduction
The two broad categories of insulators are ceramic and polymeric. The ceramic insulators are constructed from porcelain and glass. Polymeric insulators referred to as non-ceramic insulators. These are typically composed of a central rod made of fibre reinforced plastic and an outer weather shed made of silicone rubber or Ethylene Propylene Diene Monomer rubber (EPDM). The composite insulators replace porcelain insulators because of its excellent anti-pollution performance, economical maintenance, light weight and suitable for highly polluted areas [1 - 3] . To ensure the safe and reliable operation of ac transmission lines, a great number of research on ac flashover performance of polluted insulators have been conducted [1 - 15] .
The environment in which an insulator is installed can have a significant impact on the unit’s performance. When insulators are situated in areas where they are exposed to contaminate, their performance can deteriorate significantly. Natural dust, industrial waste, coastal fog etc... forms dry pollution layer on the surface of the insulators. It does not affect the operation of the insulators. While rain drops, moisture, fog etc... forms surface of insulators as wet, resistivity of the pollution layer reduced and leakage current increases exponentially and it leads to flashover of the insulator, which become an important problem for the safe operation of transmission lines and the design of the external insulation [4 - 6] .
In order to effectively assess the degree of contamination present on an insulator surface the dimensions of the insulator must taken into account. The relevant dimensions include creepage distance, diameter and height. The troublesome pollution comes in two forms. The soluble and nonsoluble components – e.g. salts from the sea and industrial gases that result in weak acids being formed, it produce an electrolyte when dissolved in water, this can be expressed as equivalent salt density deposition. And the nonsoluble materials, the nonsoluble part of the pollutant, can be expressed as the nonsoluble deposit density  .
Artificial pollution testing on insulators classified into Salt Fog and Solid Layer methods. In this solid layer method is widely used because the contaminant conditions are significant for the selection of external insulation of transmission lines and substations in pollution regions. Test procedures of Solid Layer Method are described in IEC 60507  . Table 1 shows different types of pollution test methods for polymeric insulators.
Overview of different pollution test methods applied for polymeric insulators (adopted from). PPT Slide
Lager Image
Overview of different pollution test methods applied for polymeric insulators (adopted from ).
Solid layer pollution methods are classified as Brushing Method (BM), Spraying Method (SM) and Dipping Method (DM)  . Researchers adopted different polluting methods in the pollution tests.
One of the key steps in the process of mathematical modeling is to determine the relationship between the variables. Dimensional analysis is a method to determine how variables are related for simplifying a mathematical model. Dimensional analysis alone does not give the exact form of an equation, but it can lead to a significant reduction of the number of variables. Units must be taken into consideration when collecting the data as well as when making the list of factors impacting the model and when testing the model. There is need to check that all the equations in a model are dimensionally consistent. [11 - 15] . This paper has been preceded with the proposed method which includes the experimental method of testing the insulators, developing the mathematical model with the different dimensionality functions and matching the accuracy of both the results obtained.
2. Experimental Method
- 2.1 Test equipment
Experimental investigations are carried out in the artificial pollution chamber. Test power is supplied by 100 kV test transformer. Both the AC and DC voltage can be applied for the pollution tests. In this paper, the pollution tests are carried by using AC supply voltage and the artificial pollution chamber with a height and length of 60cm.
- 2.2 Test specimens
The samples of insulators used this paper are Disc type Glass insulator and Composite insulator. Their dimensions and profile are shown in Table 2 and Fig. 1 (a) and (b) dimensions used are D is the disc diameter, H is the configuration height, and L is the creepage distance of the insulators tested  .
Dimensions of insulators tested PPT Slide
Lager Image
Dimensions of insulators tested PPT Slide
Lager Image
Profiles of insulator tested: (a) Type A; (b) Type B
- 2.3 Preparation of samples
- 2.3.1 Preconditioning of specimens
Before the pollution flashover tests are carried out, the specimens are carefully cleaned to ensure removing of all traces of dirt and grease and then dried naturally  . The specimen surfaces are coated with a very thin layer of dry kieselguhr to destroy the hydrophobicity. After one hour the completion of above procedures, surfaces of the specimens are contaminated with the suspension of sodium chloride and kieselguhr to acquire various equivalent salt deposit density levels 0.0103, 0.0319, 0.054 and 0.076 respectively, which simulate the electric and inert materials  .
- 2.3.2 Artificial pollution of insulators
Solid layer method can be employed to contaminate the surface of the insulators by three methods namely Brushing Method (BM), Dipping Method (DM) and Spraying Method (SM)  .
- 2.4 Testing of samples
Flashover voltages of insulators without any pollution are measured as per reference  and results obtained are shown in Table 3 .
Flashover voltages of the insulators without pollution PPT Slide
Lager Image
Flashover voltages of the insulators without pollution
- 2.4.1 Wetting
After polluting the insulator by one of the solid layer method above said, they are allowed for natural drying.
• The composite long - rod insulators are dried for 24 hours.
• The glass disc insulators are dried for 5 hours.
The insulators are positioned vertically at the centre of the artificial pollution chamber. Then they are completely wetted by spraying purified water with an angle of 45° before applying the test voltage. Now the pollution layer on the insulator surfaces is wetted, a series of flashover tests were carried out using AC supply.
- 2.4.2 Flashover voltage evaluation
The test used here is Even-Rising Voltage Method. Using this method the flashover voltage and the standard deviation are calculated as follows:
Increase the voltage at a random rate till the flashover occurs. Record the flashover voltage, U f . 1-2 minutes later, repeat the above procedure for five times. For the polluted specimen, the series of flashover voltage obtained by above procedure expressed as U f1 , U f2 ,... U fn (n is the number of flashovers). The minimum flashover voltage of the polluted insulator is, PPT Slide
Lager Image
Concerning five effective flashover voltages, average flashover voltage for the each insulator is calculated as: PPT Slide
Lager Image
where, U av – average flashover voltage, U fm – Minimum flashover voltage, N – Total times of effective tests and the standard deviation is calculated as: PPT Slide
Lager Image
where, (σ%) – standard deviation U fm (i) – minimum flashover voltage obtained from the test in the i th time.
- 2.4.3 Measurement of ESDD
After flashover testing, the polluted insulators are dried under bright sunlight and dry pellets of NaCl sticking to the insulator surface are collected by brushing them off with a small paint brush. The collected deposits are then dissolved in 100 ml of distilled water. The conductivity measuring instrument is used to measure the conductivity of each collected salt solution. At the same time, the temperature is also recorded. The conductivities at different temperatures are converted to conductivity at 20℃ using the formula: PPT Slide
Lager Image
where, θ is the solution’s temperature in degree Celsius, σ is the volume conductivity at the temperature θ, σ 20 is the volume conductivity at temperature 20℃ (S/m), and b is a temperature depend factor. (b=0.01905). Finally, the ESDD is determined as: PPT Slide
Lager Image
where, Vol is the volume of the solution in cm 3 , A is the area of the cleaned surface in cm 2 , S a is the salinity of the diluted liquid. It is given as: PPT Slide
Lager Image
Thus the ESDD’s at various pollution degrees are calculated and shown in Table 8 .
Flashover voltages of samples using various polluting methods. PPT Slide
Lager Image
Flashover voltages of samples using various polluting methods.
3. Mathematical Modeling
Dimensional analysis discovers the relations among physical quantities by using their dimensions. This could be done by selecting the number of fundamental dimensions and standard magnitude for these dimensions. Four fundamental physical dimensions namely, length (L), mass (M), time (T) and current (A) are used to develop the relationship among the AC flashover voltage and other parameters that affects the flashover voltage in the outdoor environment. The calculation of FOV primarily depends on the ESDD, conductivity of pollution layer (σ s ), surface area of the insulator (A), resistivity of the pollution layer (r p ) and arc constant. In this paper, four models were developed and four parameters are taken for analysis in each model. In order to utilize the algebraic approach to dimensional analysis, it is convenient to display the dimensions of the respective variables in the matrix format. From the dimensional matrix homogenous algebraic equations are formed and solved. By applying Buckingham’s Π theorem, the final dimension expression of FOV can be obtained.
- 3.1 Model 1:
In this model, FOV, ESDD, conductivity of the pollution layer (σ s ) and arc constant (N 0 ) are taken as the parameters for analysis. Therefore the relation among them can be written as: PPT Slide
Lager Image
The dimensional matrix of the respective variables can be written and shown in Table 4 .
Dimensional matrix of Model – I PPT Slide
Lager Image
Dimensional matrix of Model – I
Homogeneous linear algebraic equations are formed and solved through the above dimensional matrix. The values of k 2 , k 3 and k 4 in terms of k 1 (which represents the FOV) can be expressed as: PPT Slide
Lager Image
By assigning k 1 =1 and Buckingham’s π theorem, the final dimensional expression describing AC FOV is obtained as: PPT Slide
Lager Image
By applying Buckingham’s π theorem, the above equation can be written as: PPT Slide
Lager Image
where dc is the dimensional constant. Therefore the final dimensional expression describing FOV can be written as: PPT Slide
Lager Image
The values of equation (11) is known, the FOV are calculated and shown in Table 10 .
Comparison of calculated results with experimental results of type A Insulator PPT Slide
Lager Image
Comparison of calculated results with experimental results of type A Insulator
- 3.2 Model 2
In this model, FOV, ESDD, surface area of the insulator and arc constant are taken as the parameters for analysis. Therefore the relation among them can be written as: PPT Slide
Lager Image
The dimensional matrix of the respective variables can be written and shown in Table 5 .
Dimensional matrix of Model – II PPT Slide
Lager Image
Dimensional matrix of Model – II
Homogeneous linear algebraic equations are formed and solved through the above dimensional matrix. The values of k 2 , k 3 and k 4 in terms of k 1 (which represents the FOV) can be expressed as: PPT Slide
Lager Image
By assigning k1=1 and Buckingham’s π theorem, the final dimensional expression describing FOV is obtained as: PPT Slide
Lager Image
By applying Buckingham’s π theorem, the above equation can be written as: PPT Slide
Lager Image
The final dimensional expression describing FOV can be written as: PPT Slide
Lager Image
The values of equation (16) is known, the FOV are calculated and shown in Table 10 .
- 3.3 Model 3
In this model, FOV, Pollution severity factor (p), ESDD and arc constant (N 0 ) are taken as the parameters for analysis. Therefore the relation among them can be written as: PPT Slide
Lager Image
The dimensional matrix of the respective variables can be written and shown in Table 6 .
Dimensional matrix of Model – III PPT Slide
Lager Image
Dimensional matrix of Model – III
Homogeneous linear algebraic equations are formed and solved through the above dimensional matrix. The values of k 2 , k 3 and k 4 in terms of k 1 (which represents the FOV) can be expressed as: PPT Slide
Lager Image
By assigning k1=1 and applying Buckingham’s π theorem, the final dimensional expression describing FOV is obtained as: PPT Slide
Lager Image
By applying Buckingham’s π theorem, the above equation can be written as: PPT Slide
Lager Image
The final dimensional expression describing FOV can be written as: PPT Slide
Lager Image
The values of equation (21) is known, the FOV are calculated and shown in Table 10 .
- 3.4 Model 4
In this model, FOV, Resistance of the pollution layer (r p ), Surface area of the insulator (A) and arc constant (N 0 ) are taken as the parameters for analysis. Therefore the relation among them can be written as: PPT Slide
Lager Image
The dimensional matrix of the respective variables can be written and shown in Table 7 .
Dimensional matrix of Model – IV PPT Slide
Lager Image
Dimensional matrix of Model – IV
Homogeneous linear algebraic equations are formed and solved through the above dimensional matrix. The values of k 2 , k 3 and k 4 in terms of k 1 (which represents the FOV) can be expressed as: PPT Slide
Lager Image
By assigning k1=1 and applying Buckingham’s π theorem, the final dimensional expression describing FOV is obtained as: PPT Slide
Lager Image
By applying Buckingham’s π theorem, the above equation can be written as: PPT Slide
Lager Image
The final dimensional expression describing FOV can be written as: PPT Slide
Lager Image
The values of equation (26) is known, the FOV are calculated and shown in Table 10 .
4. Results and Discussion
- 4.1. Experimental analysis
The experimental investigation on the pollution flashover performance of composite long rod insulator and Disc type Glass insulator has been carried out. The average FOV U av (in kV) and standard deviation (%) are shown in Table 8 . for three pollution methods (BM, SM and DM) with different ESDD’s 0.0103, 0.0319, 0.054 and 0.076 mg/cm 2 .
While comparing the FOV of pure insulators and artificially polluted insulator for various ESDD’s, pure insulators shows high withstand ability than polluted insulators shown in Fig. 2 (a) and (b) . FOV of pure Type B insulator is 100 kV and the FOV of polluted Type B insulator with ESDD 0.0103 mg/cm 2 , 0.0319 mg/cm 2 , 0.054 mg/cm 2 , and 0.076 mg/cm 2 is 91.5 kV, 74.8 kV, 61 kV and 45.8 kV respectively. Increase in pollution level on surface of insulators, FOV decreases. The chief cause behind this decrease in FOV is, pollution layer act as a conductor (i.e) it forms a conduction layer on surface of the insulator and it leads to immediate flashover. PPT Slide
Lager Image
Flashover voltage Vs ESDD for insulators Type A & B: (a) Type A insulator; (b) Type B insulator
The average pollution FOV of various insulators decreases with increase in ESDD. For example, BM with the value of ESDD is 0.0103 mg/cm 2 , 0.0319 mg/cm 2 , 0.054 mg/cm 2 , and 0.076 mg/cm 2 , the U av of Type A insulator is 44.5 kV, 33.3 kV, 23.8 kV and 15.5 kV respectively. From the data above we can see that U av decreases by respectively when the ESDD increases.
Average pollution FOV of various insulators is affected by different types of polluting methods. The polluting methods are BM, SM and DM, and the U av values of Type A insulator is 44.5 kV, 41 kV, 39.5 kV respectively for ESDD 0.0103 mg/cm 2 , this inferred that the U av values of BM with various ESDD is higher than that of SM and DM pollution methods. Similarly for Type B insulator U av values of BM is higher than that of SM and DM pollution methods and shown in Fig. 3 (a) and (b) . PPT Slide
Lager Image
Flashover voltage Vs ESDD for different polluting methods: (a) Type A insulator; (b) Type B insulator
Authors witnessed standard deviation σ (%) of U av of three pollution methods BM, SM and DM for Type A and B insulator are 1.08% – 3.62%, 0.99% – 2.91% and 0.99% – 3.02% respectively. In these, standard deviation of BM is large compared to that of SM and DM, because of non-uniform pollution layer on the surface of the insulator and this is the main drawback while applying BM pollution on the surface of the insulator. The advantage is, it’s easy and simple.
Based on Table 8 , relative deviation between BM, SM, and DM pollution methods are calculated and shown in Table 9 . The relative deviation of FOV depends on type of insulator. Deviation range of Δ1 % and Δ2 % vary from 1.308% − 19% and (−2.702%) − (−7.142%) respectively for Type A and B. The Δ1% and Δ2% will vary at different ESDDs for the same type of insulator. For example, Type A insulator ESDD is 0.0103−0.076 mg/cm 2 , the relative deviation of FOV for Δ1% and Δ2 % from are 8.292% − 19% and (−3.571%) − (−7.142%), respectively.
Results of deviations in polluting methods. PPT Slide
Lager Image
Results of deviations in polluting methods.
Maximum standard deviation σ (%) of U av is only 3.62. but the deviations Δ1 % of U av between BM and SM is up to 19% at various ESDD’s for Type A and B insulators. This deviation is not acceptable one. U av between DM and SM is less than −2.702, this deviation is acceptable and it could be negligible.
- 4.2. Mathematical analysis
While comparing the mathematical models (1 to 4) results with experimental results of Type A insulator, mathematical models (1 to 4) results are lesser than experimental results. For example, the experimentally attained FOV of Type A insulator with ESDD 0.0103 mg/cm 2 is 44.5 kV and mathematical models (1 to 4) results for same ESDD is 42.07 kV, 42.25 kV, 42.83 kV and 41.36 kV respectively
Much attention is given to compare the factor that affects the FOV is mathematical models (1 to 4). The experimentally attained FOV of Type A insulator with ESDD 0.0103 mg/cm 2 is 44.5 kV and mathematical models (1 to 4) FOV for same ESDD is, 42.07 kV, 42.25 kV, 42.83 kV and 41.36 kV respectively. The experimentally attained FOV of Type B insulator with ESDD 0.0103 mg/cm 2 is 91.5 kV and mathematical models (1 to 4) FOV for same ESDD is 87.34 kV, 87.7 kV, 89.63 kV and 86.54 kV respectively. It’s inferred that experimentally attained FOV is higher than that of mathematically attained FOV; Model 3 with ESDD 0.076 mg/cm 2 is greater than experimental FOV for Type B insulator, similarly Model 1 with same ESDD and shown in Table 10 . FOV of Model 3 is nearly same with experimental FOV for both type A & B insulator and Fig. 4 (a) and (b) shows comparison of results of mathematical models (1 to 4) with experimental results of Type A insulator. PPT Slide
Lager Image
Mathematical Models (1 to 4) results VS Experimental results of Type A & B insulator Type A insulator; (b) Type B insulator.
Table 11 shows the deviation of mathematical models (1 to 4) FOV from experimental FOV. In this FOV deviation with ESDD 0.0103 mg/cm 2 is (−5.78%), (−5.36%), (−3.89%) and (−7.59%) respectively for Type A insulator. Model 1 for ESDD 0.076 mg/cm 2 is 9.88%, this ESDD show positive deviation and all other ESDD’s shows negative deviation for Type A insulator. Similarly for Type B insulator FOV with ESDD 0.0103 mg/cm 2 is (−4.76%), (−4.33%), (−2.08%) and (−5.73%) respectively. Model 2 for ESDD 0.054 mg/cm 2 is 5.32%, this ESDD show positive deviation and all other ESDD’s shows negative deviation for Type B insulator.
Deviation of mathematical models (1 to 4) FOV from experimental FOV in (%) PPT Slide
Lager Image
Deviation of mathematical models (1 to 4) FOV from experimental FOV in (%)
Accuracy of mathematical models (1 to 4) FOV compared with experimental FOV shown in Table 12 . For Type A insulator, accuracy of FOV for model 3 is high compared to other models. For example, accuracy of model 3 with ESDD 0.0103 mg/cm 2 , 0.0319 mg/cm2, 0.054 mg/cm 2 , and 0.076 mg/cm 2 for Type A insulator is 96.25%, 93.38%, 95.93% and 90.13% respectively, except model 1 with an accuracy of 110.97% for ESDD 0.076 (mg/cm 2 ). For Type B insulator, accuracy of FOV for model 3 is high compared to other models.
Accuracy of mathematical models (1 to 4) FOV from experimental FOV in (%) PPT Slide
Lager Image
Accuracy of mathematical models (1 to 4) FOV from experimental FOV in (%)
5. Conclusion
From the above pollution flashover tests and analysis, some conclusions are obtained as follows:
• 1) FOV decreases with increase in pollution severity in terms of Equivalent Salt Density Deposition (ESDD).
• 2) The AC pollution FOV of polluted insulators using brushing method is greater than using spraying method and dipping method.
• 3) Mathematical model 3 yields better results than other methods for both Type A and B insulators.
• 4) Relative deviation between BM and SM is not acceptable by engineering applications because the deviation is beyond tolerance limit.
• 5) While considering accuracy of the Mathematical models (1 to 4), Model 3 shows higher accuracy compared with other Models in both Type A and B insulators.
BIO N.B.Prakash obtained his B.E Degree by the year 2000 and M.E Degree by the year 2002. Presently he is working as an Associate Professor in EEE Dept. of National Engineering College. His research interest is Digital Signal and Image Processing, Embedded System and High Voltage Engineering. M.Parvathavarthini obtained her B.E Degree by the year 2010 and pursued her M.E Degree in High Voltage Engineering at National Engineering College. Her research interest is Condition monitoring of High voltage apparatus. R.Madavan obtained his B.E Degree by the year 2008 and M.E Degree by the year 2011. Presently he is working as an Assistant Professor in EEE Dept. of National Engineering College. His research area is High Voltage and Insulation Engineering and Condition monitoring of High voltage apparatus.
References