Omar Ahmed, Anne Goj Transformer Protector Corp, Humble, TX
Abstract
Power transformers are the cornerstone of the electricity grid and their protection is crucial for all the generation, transmission and distribution companies around the world. Power transformer explosions are common and lead to fire, risk to human life, power outages, environmental pollution, generating huge expenses.
Therefore, a strategy to prevent transformer tank rupture was developed [1]. This strategy is based on the direct mechanical response of up to several fast Depressurization Sets (DS) to the first dynamic pressure peak induced by electrical faults, which we define as the TRANSFORMER PROTECTOR (TP). In the same paper, the operation of this fast depressurization strategy was experimentally validated: if the oil and gases evacuation through the fast Depressurization Set is activated within milliseconds by the first dynamic pressure peak before static pressure increases, the transformer explosion can be prevented.
In order to minimize the costs and the risks of experimental testing, a compressible two-phase flow numerical simulation tool based on a 3D finite volume method was developed to study transformer explosions and possible strategies for their prevention. The current paper shows the applications of this simulation software towards studying the influence of real mechanical transformer protections. In this paper, simulations on a transformer geometry with and without Oil Bushing Cable Box (OBCB) protections will be presented and their results discussed.
1 Introduction
Power transformers are highly dangerous due to the large quantity of oil they contain in combination with high voltage elements. When the oil loses its dielectric properties, low impedance faults may occur, which result in arcing. The first dynamic pressure peak increases at the location of the arc, leading to pressure waves propagating at the speed of sound within the oil, about 1000m/s. These pressure waves reflect off tank walls which build up static pressure within the transformer tank, ultimately resulting in a potentially violent tank rupture and fires when static pressure reaches the tank limit.
Transformers are generally built to withstand a static pressure of 0.7 bar (10.15 psi) gauge pressure. In 2004 the independent Brazilian high voltage laboratory CEPEL conducted a series of experiments on transformer arcing and depressurization [1]. Simulations show that the internal overpressure for the parameters of the CEPEL tests would have resulted in static pressures of 7 bar (101.5 psi) gauge, 10 times greater than the withstand limit of the transformer tanks.
Section VIII of the ASME code defines how to determine the thickness of thin shells under internal pressure. As an example, Division I, UG-27 is used to determine the thickness of a spherical shell. If we make the assumption that the internal pressure is much less than the product of the maximum allowable stress of a material and the efficiency of the bolts/welds, then we can conclude that the transformer tank should have been approximately 10 times thicker to avoid transformer explosion without protection [1]. Such a tank design would induce additional economical and industrial strains. For this reason, an alternative explosion prevention method is investigated.
The strategy reported in this paper involves the rapid evacuation of oil and gases within milliseconds of the arc event by means of a passive mechanical system, called the TRANSFORMER PROTECTOR. A
complete experimental study of this strategy was previously reported in [1]. As live tests are both highly expensive and dangerous, an alternative computational simulation tool has been developed based on the results of the CEPEL tests.
The aim of this paper is to present the results of the application of this computer simulation tool towards understanding the operation of a rapid oil evacuation strategy under a variety of different situations. First, a computational simulation tool is briefly discussed. Then, the parameters for the simulations are presented. Finally, the time evolution of an internal overpressure inside protected and unprotected transformers is detailed. The simulation tool is used to look at the effect of the location and number of fast tank Depressurization Sets, under two arc configurations.
2 Computational Simulation Tool
The described physical phenomena are modeled as a 3D compressible two-phase flow, using a set of partial differential equations based on a 5 equation model developed in [2]. These equations represent the conservation of mass, momentum, and energy, as well as the advection of the volume fraction for each phase. Effects such as gravity, viscosity, and heat conduction are added in the modeling to adhere to physical constraints. One of the major characteristics of the model is its ability to accurately depict the pressure wave propagation inside liquids and gases. A finite volume algorithm is adopted to transform the system of differential equations into algebraic equations, and therefore computationally solve the system. The volumes are defined by an unstructured 3D mesh, to allow a precise description of complex geometries such as transformer tanks.
In the test series, the equations are solved allowing for phase change through cavitation. This properly models extremely low local pressures as a consequence of fast depressurization.
3 Simulation Description
In order to study the fast direct-tank depressurization protection strategy, it is important to accurately model the transformer geometry, understand the initialization of the phenomena, and account for the tank explosion prevention technology.
3.1 Transformer Geometry
As there are no concrete standards for the design of transformer tanks, we chose to focus on the components and dimensions which are characteristic of a real transformer tank. A generic 400 MVA transformer tank geometry has been created for the purpose of simulations, and is shown in Figure 1. This geometry models the transformer with three inner cores but no other internal components (such as support beams and cables). This geometry was used in three primary configurations:
1) No Protection
2) Single Horizontal Depressurization Set (HDS)
3) Single HDS + 3 Oil Bushing Cable Box (OBCB) protections (one each on turrets 1, 2 and 3)
With two arc locations:
1) “Left” – Entrance to turret 1
2) “Right” – Entrance to turret 3
The final meshes consisted of approximately 150,000 tetrahedral elements.
Rather than downstream piping elements from an actual TP, atmospheric boundary conditions are used at the locations of the Depressurization Sets.
3.2 Description of the Electric Arc
Tests on industrial scale transformers [1] provide detailed information on electrical arc phenomenon.
a) Overview: When a low impedance arc occurs in an oil insulated transformer, the liquid oil surrounding the arc location is quickly vaporized and a gas bubble is created. As the liquid oil does not allow for expansion, the gas bubble becomes highly pressurized at the point of vaporization. This pressure difference between the gas bubble and the surrounding oil creates a first dynamic pressure peak, leading to dynamic pressure waves when interacting with the tank structure. These dynamic pressure waves propagate throughout the transformer, at the speed of sound, building static pressure within the
b)After the initial oil vaporization and gas bubble creation, the electrical arc will no longer be in direct contact with the transformer oil, and therefore further vaporization is limited. Instead, the arc energy transferred to the gas is used to heat the gas up, generating plasma. However, the cracking of the vaporized oil does not have as pronounced effect on the pressure wave generation as the initial vaporization process.
c) Study Arcing Conditions: In the “left” arc location, the arc occurs in the transformer tank at the connection to turret 1. In the “right” arc location, the arc occurs in the transformer tank at the connection to The arc is simulated by placing a pressurized volume of gas, or “bubble”, in the transformer. Then, over the duration of the arc, energy is added into the region, causing further rise in pressure due to the equation of state. This method for modeling the arc has been used for previous studies such as the parametrical study on TP configurations [3].
The bubble is roughly cylindrical in shape, with a length and radius of approximately 0.4m and 0.2m, respectively, and under 11 bars of pressure. However, due to mesh variations, the result is not exactly identical.
3.3 Modeling the Protection
Simulating the operation of the depressurization system is a difficult task. This has been accomplished in the computational tool by removing the separation between the oil and the Decompression Chamber as a function of the pressure difference across the barrier. Furthermore, the opening function is constrained to open in progressively from the center first, and then outwards. A comparison of the oil evacuation from live tests (left) and simulations (right) is shown in Figure 2.
3.4 Simulation Parameters
The simulation parameters used for the tests are described in Table 1, Table 2 and Table 3, below. The
oil and gas are modeled with the stiffened gas equation of state.
Density (kg/m3) |
Gamma (unit-less) |
Stiffness (Pa) |
|
Oil |
843.0 |
1.7 |
7.0*108 |
Gas |
1.0 |
1.43 |
0.0 |
Table 1 - Fluid Properties
Alpha |
Pressure (Pa) |
|
Transformer Tank |
1.0*10-5 |
1.0*105 |
Arc “bubble” |
0.99999 |
1.1*106 |
Outlet (when open) |
0.99999 |
1.0*105 |
Table 2 - Initial Conditions
Current (Peak) |
Voltage (Peak) |
Frequency |
Duration |
Total Energy |
2.0*103 A |
80.0*103 V |
50.0 Hz |
0.110 s |
7.1986 MJ |
Table 3 - Arc Parameters
4 Simulation Results
Table 4 shows the peak pressures for the six tests. Table 5 shows the activation times of the individual DS and the Depressurization Time of the transformer tank. The Activation Time is defined as the time from the beginning of the arc to the opening of the DS. The Depressurization Time is defined as the period from the beginning of the arc to the final time the dynamic pressure dips below the static pressure withstand limit of the transformer tank, the point after which tank rupture is no longer possible.
The Activation Times for the left arc location show the OBCB protection is 20% faster than the Activation Time of the HDS alone. Switching to the right arc location the OBCB activation is 70% faster than the only HDS Activation Time.
Tank * |
Turret 1 |
Turret 2 |
Turret 3 |
|
“Left” - No Protection |
19.3 bar |
15.0 bar |
10.7 bar |
18.2 bar |
“Left” - HDS |
16.8 bar |
12.4 bar |
9.4 bar |
15.5 bar |
“Left” - HDS + 3 OBCB |
15.0 bar |
12.1 bar |
7.0 bar |
9.3 bar |
“Right” - No Protection |
18.6 bar |
17.7 bar |
11.9 bar |
15.5 bar |
“Right” - HDS |
17.0 bar |
15.8 bar |
10.1 bar |
12.4 bar |
“Right” - HDS + 3 OBCB |
15.0 bar |
11.5 bar |
6.6 bar |
11.5 bar |
* Peak “Tank” pressure may occur at any point in the transformer, including turrets
Table 4 - Peak Pressures in Transformer
Depressurization Time |
Activation HDS |
Activation OBCB Turret 1 |
Activation OBCB Turret 2 |
Activation OBCB Turret 3 |
|
“Left” - No Protection |
∞ |
n/a |
n/a |
n/a |
n/a |
“Left” – HDS |
400 ms |
2.07 ms |
n/a |
n/a |
n/a |
“Left” - HDS + 3 OBCB |
146 ms |
1.95 ms |
1.56 ms |
2.84 ms |
5.33 ms |
“Right” - No Protection |
∞ |
n/a |
n/a |
n/a |
n/a |
“Right” – HDS |
350 ms |
5.33 ms |
n/a |
n/a |
n/a |
“Right” - HDS + 3 OBCB |
139 ms |
5.32 ms |
5.97 ms |
3.08 ms |
1.65 ms |
Table 5 - Depressurization and Activation Times
“Left” Arc Location |
“Right” Arc Location |
Figure 3 - No Protections |
The Depressurization Times are much faster for the case with the OBCB protections. In the left arc case the OBCB Depressurization Time is 36.5% the HDS alone. The right arc case Depressurization Time is 39.7% the time taken by the HDS alone. These lower Depressurization Times are approximately 2.6 times faster than the HDS only cases.
4.1 Effect of the Arc Location
The results show a large degree of symmetry based on the location of the arc. Figure 3 shows the pressure profiles on the unprotected transformer. From these figures, we can see that the pressure profiles in the two tests are very similar, where the profiles of Turret 2 are nearly the same, and the profiles of Turret 1 and Turret 3 are mirror images of each other when switching the arc location.
We note some important differences between Left and Right locations. Table 6 gives a comparison of the peak pressures at important points in the transformer. Although the pressure profiles are similar, pressure peaks vary by as much as 10%.
Pressure (bars) |
|||
“Left” |
“Right” |
Difference |
|
Main Tank |
19.3 |
18.6 |
- 0.7 |
Far Turret |
18.2 |
17.7 |
- 0.5 |
Near Turret |
15.0 |
15.5 |
+ 0.5 |
Middle Turret |
10.7 |
11.9 |
+ 1.2 |
Table 6 – Peak Pressure Comparisons based on Arc Location – Unprotected Transformers
Finally, we note that with OBCB protections regardless of where the initial fault occurs a DS will likely be close to the arc. In the left arc location the first DS activates
1.25 times earlier. The right case is even faster activating 3 times sooner than the HDS only due to its proximity to the arc.
4.2 Effect of the OBCB protections
We can see a trend that the additional OBCB protections allows for the substantially more quickly activation of a DS, if an arc is sufficiently close to a turret.
More importantly, the results show that the addition of 3 OBCB protections reduces the Depressurization Time as compared to that of a single horizontal DS. As shown in Table 5, the Depressurization Time of the 3+1 configuration was, on average, 2.6 times the Depressurization Time in the single HDS in the simulations performed.
Furthermore, the simulations show that the additional OBCB protections help limit the peak pressures experienced at the pressure sensors located in all three bushing turrets, as seen in Table 4. The pressure evolution of both cases is presented in Figure 4 and 5.
“Left” Arc Location |
“Right” Arc Location |
Figure 4 - Single Horizontal Depression Set |
“Left” Arc Location |
“Right” Arc Location |
Figure 5 - Horizontal Depressurization Set + 3 OBCB Protections |
5 Conclusions
A three dimensional simulation tool based on a five equation hydrodynamic model has been used to investigate a protection strategy against transformer tank rupture. The use of this protection strategy serves as an alternative to designing transformer tanks as pressure vessels, which by application of the ASME standard would require an impractical factor of ten increase in the tank thickness. By applying the simulation tool, the efficiency of a fast tank depressurization method has been studied.
Simulations have been performed on a 400 MVA transformer tank to determine the effect of adding three Oil Bushing Cable Box protections in addition to a single Horizontal Depressurization Set. As modeled, the turrets on which the additional OBCB protections are installed are openly connected to the main tank, without internal obstructions between the main tank and bushing turrets.
These simulations indicate that the additional protections reduce the overall depressurization time of the transformer by a factor of 2.6 by activating up to 3 times earlier. Furthermore, the OBCB protections reduce the peak pressure in the bushings, which are vulnerable to rupture.
6 References
- S. MULLER, R. BRADY, G. de BRESSY, P. MAGNIER, G. PERIGAUD, “Prevention of Transformer Tank Explosion, Part 1: Experimental Tests on Large Transformers”, ASME PVP 2008.
- R. BRADY, S. MULLER, G. de BRESSY, P. MAGNIER, G. PERIGAUD, “Prevention of Transformer Tank Explosion, Part 2: Development and Application of a Numerical Simulation Tool”, ASME PVP 2008
- R. BRADY, S. MULLER, M. PETROVAN- BOIARCIUC, G. PERIGAUD, B.LANDIS, “Prevention of Transformer Tank Explosion, Part 3: Design of Efficient Protections using Numerical Simulations”, ASME PVP 2009
7 Biography
Omar Ahmed is currently working at Transformer Protector Corporation as a Research Engineer, where he is modelling the physics associated with transformer explosion and depressurization strategies. He specializes in optimizing the computational fluid dynamics algorithms to model energy transfer from the arcing event to the transformer oil, and the subsequent pressure wave propagation.
Omar completed his undergraduate degrees in Mathematics and Physics at the University of Texas, Austin in 2006, and completed his Masters in Geophysical Fluid Dynamics at Rice University in 2009.
Anne Goj is currently working at Transformer Protector Corporation as a Research Engineer where she spends so much time calculating quantities with physics that she occasionally wonders how her degrees all have “chemistry” written on them.
Anne studied theoretical and computational physical chemistry at Cornell University before relocating to Texas in 2007.