First Seminar (14.10.2024) - 2
First Seminar (14.10.2024) - 2
ABSTRACT
60076-7 standard and IEEE C57.91 Standard are used for calculating thermal
characteristics and the loss of life of distribution transformer.
1. INTRODUCTION
1.1. Background
However, these transformers generate heat due to electrical losses, which can
affect their efficiency and lifespan. Managing the internal temperature of transformers
is crucial, as excessive heat can lead to insulation degradation, hot spots, and even
transformer failure. With the rising concern over global climate change and increasing
ambient temperatures, transformers are exposed to higher thermal stresses than
anticipated in their original design. This poses challenges in design and maintaining
their reliability and extending their operational life.
The existing designs of many distribution transformers do not adequately account for
the current and projected increases in ambient temperatures. Moreover, traditional
thermal management techniques may not suffice in maintaining optimal thermal
performance under extreme environmental conditions. Therefore, there is a pressing
need to develop improved thermal analysis models and design strategies that take into
account the rise in ambient temperatures to enhance transformer resilience, efficiency,
and longevity.
The main aim of this thesis is to study the Design Improvement of Distribution
Transformer (1000KVA) with Thermal Analysis under Ambient Temperature Rise
Conditions. It also aims to analyze the thermal behavior of this transformer under
different ambient temperature conditions and to develop constructional design
improvements that enhance their thermal performance, reliability, and lifespan.
2. LITERATURE REVIEW
The history of transformers begins in the early 19th century with important
discoveries in electromagnetism. Hans Christian Ørsted showed that electric currents
create magnetic fields, and Michael Faraday's work in 1831 led to the principle of
electromagnetic induction, which is the basis for transformer operation. These early
findings laid the groundwork for later technological developments.
In recent decades, the rise of modern thermal analysis techniques has allowed
for more precise modeling of heat flow in transformers. Advanced materials, such as
amorphous metals, and better cooling technologies have further reduced thermal stress.
Additionally, as climate change leads to higher ambient temperatures, thermal analysis
plays an increasingly important role in ensuring transformers can withstand rising
environmental heat.
innovation was its use of a closed iron core, which greatly enhanced the efficiency of
magnetic induction between the two coils.
At the same time, William Stanley, working for George Westinghouse in the
United States, created a transformer based on similar principles. His design further
refined the technology, making it more practical for commercial use. Stanley's
transformer was crucial in the development of AC power distribution systems, as it
allowed power to be transmitted over long distances with minimal losses and then
stepped down for safe local use.
Figure 2.1. William Stanley's First Transformer built in 1885, Single Phase AC Power
7
Figure 2.2. William Stanley's First Transformer which was used in the Electrification
of Great Barrington
2.3. Introduction of Distribution Transformer
2.3.1. Operating Principle
A distribution transformer is a crucial component in the electrical power
distribution system, used to step down high-voltage electricity to a lower voltage
suitable for use in homes and businesses. The operating principle of a transformer is
based on electromagnetic induction. When an alternating current (AC) passes through
the primary winding, it creates a magnetic field in the transformer’s core. This varying
magnetic field induces a voltage in the secondary winding, and depending on the
winding ratio, the voltage is either increased (step-up transformer) or decreased (step-
down transformer). In a distribution transformer, the goal is to reduce the high
transmission voltage to the desired level for distribution.
2.3.2 Construction
The basic construction of a distribution transformer consists of the following parts:
2. Windings: There are two types of windings: primary and secondary. The
primary winding receives the high voltage from the transmission lines, while
the secondary winding delivers the stepped-down voltage to the distribution
network.
9
4. Main Tank: Tanks are made up of fabricated rolled steel plates. They are
provided with lifting hooks and inbuilt cooling tubes. In order to minimize the
weight and stray losses, aluminum sheets are also being used instead of Steel
plates. However, due to its lightweight property, nowadays aluminum tank is
more familiar and costly than a steel tank. The main tank is the robust part of
the transformer that serves mainly two purposes:
10
a. It protects the core and windings from the external environment and
provides housing for them.
b. It is used as a container for transformer oil and provides support for
all other external accessories of the transformer.
5. Terminals and Bushings: Terminals and bushings are also important parts of
transformer that are used to connect incoming and outgoing cables of supply
and load. These are connected with the ends of the windings conductor.
Bushings are mainly insulators made up of porcelain or epoxy resins. They are
mounted over the tank and form a barrier between the terminals and the tank.
They provide safe passage for the conductor connecting terminals to the
windings. As windings are of two types and so bushings are also of two types
as named below:
a) High-Voltage Bushings
b) Low-Voltage Bushings
11
to the variation in the ambient temperature of transformer oil inside the main
tank. It is a cylindrical drum-type structure installed on the top of the
main tank of the transformer. It is connected to the main tank through a
pipe and a Buchholz relay mounted on the pipe.
10. Breather: The breather is a cylindrical container filled with silica gel and
directly connected to the conservator tank of the transformer. The main function
of the breather is to supply moisture-free fresh air to the conservator tank during
the expansion and contraction of transformer oil. This is because the transformer
oil when reacting with moisture can affect the insulation and cause an internal
fault in a transformer.
11. Explosion Vent: An explosion vent is a metallic pipe with a diaphragm at one
end and installed on the main tank slightly above than conservator tank. It is
available only in high-rated power transformers. The main function of the
explosion vent is to protect the power transformer against explosion during
excessive pressure build-up in the main tank due to severe internal faults. It acts
as an emergency exit for oil and hot air gases inside the main tank of the
transformer. The explosion vent works on the same principle as the safety valve
works in the pressure cooker. Hence In other words we can also call the
explosion vent as safety valve of the transformer.
12. Radiator and Fans: Since Power losses in the transformer are dissipated in the
form of heat. So a cooling arrangement is required for the power transformer.
14
Losses in a Transformer
Iron (core) losses consists of hysteresis and eddy current losses, both of which
occur in the magnetic core due to the alternating magnetic field. Hysteresis loss in the
core is related to the material's magnetic properties and can be estimated using
following equation is:
= × × .
× (2.1)
Where;
Is hysteresis losses (W)
Is Constant specific to the core material
Is frequency of the AC (Hz)
.
Is Maximum magnetic flux density (T)
Is Volume of the core ( m3)
Eddy current loss is due to induced currents in the core and can be expressed as:
= × × × × (2.2)
Where;
Is Eddy current losses (W)
Is constant related to core material and conductivity
Is thickness of the core laminations (m)
Winding (core) losses are due to the current flowing through the transformer
windings and are load-dependent. Copper losses increase with load, so as the
transformer experiences higher load currents, the heat generated in the windings rises.
The copper loss can be represented by the equation is:
= × (2.3)
Where;
Is Copper losses (W)
Is current through the windings (A)
Is Resistance of the windings (Ω)
Combining both copper and iron losses, the total heat generation can be
expressed as:
= + (2.4)
18
Where;
Is Conduction heat transfer rate (W)
Is thermal conductivity of the material (W/m·K)
" Is cross-sectional area through which heat is being conducted (m²)
#
Is temperature gradient (K/m)
of these temperatures help prevent insulation degradation, reduce the risk of failures,
and ensure long-term reliable operation.
The top-oil temperature rise refers to the difference between the temperature
of the oil at the top of the transformer tank (the hottest part of the oil) and the ambient
temperature. Top-oil temperature rise can be expressed as:
The top-oil temperature rise depends on the load current, losses, cooling
method, and transformer design. The IEEE Loading Guide (IEEE C57.91) provides
equations to estimate the top-oil temperature rise as a function of transformer load and
cooling conditions. For a specific load condition, Governing Equations for Top-Oil
Temperature Rise ∆'#1 can be estimated as:
53×(6( ))7
∆0#1 ( ) = ∆0#1,3 × 4 8 (2.10)
53
Where;
∆0#1 ( ) Is the top-oil temperature rise at time (t)(℃).
∆0#1,3 Is the rated top-oil temperature rise (℃)
Is the ratio of load losses to no-load losses at rated load
:( ) Is per-unit load at time (t)
; Is the oil exponent (typically 0.8 to 1.0)
The hot-spot temperature rise is the temperature rise of the hottest point
within the transformer winding relative to the surrounding oil. This "hot-spot" typically
occurs deep within the windings, where heat dissipation is most limited. The hot-spot
temperature is the most critical factor affecting the life expectancy of the insulation. It
can be expressed as:
22
Hot-spot temperature rise is the primary factor affecting the aging rate of
transformer insulation. Insulating materials, such as paper and pressboard, degrade
faster at higher temperatures. The relationship between hot-spot temperature and
insulation aging is exponential. The aging acceleration factor (?@@ ) is given by:
DEF GHHI
ABB = 24 J
8
(2.13)
Where;
ABB Is the aging acceleration factor
0<= Is the hot-spot temperature
110°C Is the reference temperature at which insulation aging is considered to
proceed at a normal rate.
When the hot-spot temperature exceeds the design limits, the insulation’s life
expectancy decreases significantly. Managing hot-spot temperature is crucial for
maintaining the transformer's health and ensuring long-term reliability.
23
The total hot-spot temperature rise is the sum of the ambient temperature, the
top-oil temperature rise, and the hot-spot temperature rise:
Where;
0( ) Is the transformer temperature at time (t)(℃).
0N Is Final steady-state temperature(℃).
K Is thermal time constant (hours)
2.6.6 Loss of Life Calculation
The loss of life is calculated based on how long the transformer operates at
different temperatures relative to its designed life at the reference temperature (usually
110°C). The key equation used is based on the aging acceleration factor (FAA), which
depends on the transformer’s hot-spot temperature.
DEF (X)GDYZ[
W ]
:OPP O :R S (T. U. ) = V1 S \ ^ (2.16)
Where;
0<= ( ) Is hot-spot temperature at time (t)(℃).
24
They have been used in the industry and will be explained in more details.
25
This type of cooling is most common used in the practice. The ONAN cooling
transformers have their windings and core cooled by oil naturally circulation
(thermosiphon effect). The oil is then cooled by air which is naturally cooled.
Moreover, Oil is kept in circulation by the gravitational buoyancy in the closed-loop
cooling system as shown in Figure 2.13.
As the transformer losses increases, the number and size of the radiators that are
required to cool the oil must increase. Eventually, a point is reached where air and
26
natural convection are not adequate to remove the heat and air must be forced through
the radiators by motor- driven fans (John J. Winder, Jr, 2002), This type of cooling is
termed as ONAF (Oil Natural and Air Forced) cooling, as shown in Figure 2.14. These
fans generate forced air flow in the radiators to increase the oil circulation rate from
tank to radiators. This type of cooling improves the efficiency of the transformer to
increase the capability to operate at a high load, but this flow rate is relatively low.
Because of this, the heat carrying (or dissipating) capacity of the oil is low. The heat
carrying capacity can be defined as
= _. `a (' −' ) (2.17)
Where Q is heat flow in W, m is mass flow rate in kg/s, Cp is specific heat in J/ (kg °
C) and temperatures Tout and Tin are in °C. For the given transformer oil inlet (Tin) and
top oil (Tout) temperatures.
Although the cooling capacity is greatly increased by the use of forced air, increasing
the loading to take advantage of the increased capacity will increase the temperature
gradients within the transformers. A point is reached where the internal temperature
gradient limit the ability to increase load any further. The solution is to increase the oil
velocity by pumping oil as well as forcing air through the radiators. This type of cooling
called OFAF (Oil Forced and Air Forced) cooling. The usual pump placement is at the
27
bottom of the radiators as shown in Figure 2.15, forcing oil from the radiator outlets
into the bottom of the transformer tank in the Sam direction as natural circulation but
at a much higher velocity (John J. Winder, Jr, 2002). By direction the flow of oil within
the transformer windings, reduce considerably the temperature between the top and the
bottom of the radiators (M.T. Ishak, 2009).
In OFAF cooling type the oil forced from the pump into the transformer, its
flow is governed by the least resistance path as well as the buoyancy. Hence, part of the
oil may not enter either windings or core, and may from parallel path outside these two.
Thus, the top oil temperature may reduce because of the mixture of hot oil coming from
the windings and the cool oil coming from the pump. This in turn reduces the
effectiveness of radiators. The solution is to improve the heat dissipation rate by the oil
forced and directed in the windings through the paths. This type of cooling called
ODAF (Oil Direct and Air Forced) cooling, as shown in Figure 2.16, (S.V.Kulkani and
S.A. Khaparde, 2004).
28
Key
A Top-oil temperature derive as the average of the tank outlet oil temperature and
the tank oil pocket temperature
B Mixed oil temperature in the tank at the top of the winding (often assumed to be
the same temperature as A)
C Temperature of the average oil in the tank
D Oil temperature at the inlet of the tank (assumed to be the same as at the bottom
of the winding)
E Bottom of the tank
gr Average winding to average oil (in tank) temperature gradient at rated current
H Hot-spot factor
P Hot-spot temperature
I Average winding temperature determined by resistance measurement
x-axis Temperature
y-axis Relative positions
measured point; calculated point
a) The oil temperature inside the tank increases linearly from bottom to top,
whatever the cooling mode.
b) As a first approximation, the temperature rise of the conductor at any position
up the winding is assumed to increase linearly, parallel to the oil temperature
rise, with a constant difference gr between the two straight lines (gr being the
difference between the winding average temperature rise by resistance and the
average oil temperature rise in the tank).
c) The hot-spot temperature rise is higher than the temperature rise of the
conductor at the top of the winding because allowance has to be made for the
increase in stray losses, for differences in local oil flows and for possible
additional paper on the conductor. To take into account these non-linearity’s,
the difference in temperature between the hot-spot and the top-oil in tank is
made equal to H×gr, that is, ∆θhr = H×gr.
30
Transformer losses are broadly classified as no-load and load losses. The no-
load losses are composed of two components: (1) the hysteresis loss and (2) eddy
current loss. The hysteresis loss is the energy loss when the magnetic material goes
through a cycling state. The eddy current loss is caused when the flux lines pass through
the core. Load losses occur when the output is connected to a load so that current flows
through the transformer from input to output terminals. Losses of energy appear as heat
which is the major cause of transformer temperature rise and its aging. Among the
several factors that affect the transformer performance, the location of the maximum
temperature of the solid insulation named ‘hot spot temperature’ (HST) and its
temperature level have been identified as the main reasons for transformer aging and
failures. In fact, high temperatures in a transformer will drastically shorten the life of
insulating materials used in windings and structures. For every 8 degrees Celsius (°C)
temperature rise, the transformer’s life is cut by one-half; therefore, maintaining cooling
systems is critical. In worst cases, the high HST may lead to transformer failure and
explosion. When the internal temperature of a transformer reaches 150–300 °C under
abnormal conditions, the mineral oils produce hydrogen and methane gases due to
chemical decomposition. At temperatures above 300 °C, ethylene is formed, and above
700°C, large amounts of hydrogen and ethylene are produced. These gases tend to
dissolve partially or entirely in the mineral oil and can form combustible/flammable
mixtures if they escape the transformer oil compartment, leading to unexpected
fire/explosion accidents. Consequently, excellent thermal design is paramount in
transformer design since it determines insulation aging and life.
The thermal modeling of the transformer can be divided into two general
categories. The first category is physical and semi-physical modeling describing a
transformer’s complex heat transfer modes with simple differential equations such as
correlation, lumped, and equivalent circuit models. The other category is the
computational fluid dynamic (CFD) modeling approach based on numerical methods
such as finite element and finite volume methods, which is very accurate but time-
consuming and requires the design detail of the transformer for modeling.
31
Correlation models are analytical tools used to estimate thermal behavior and
performance characteristics of transformers based on empirical relationships and
statistical data. The correlations consist of some thermal parameters which should be
obtained by heat run tests. IEEE C57.91-2011 and IEC 60076-7 suggest well-known
correlation models for estimating the hot spot temperature of transformers. In these
models, the hot spot temperature is assumed to consist of three components: the average
ambient temperature, the top oil rise over ambient temperature, and the winding hot
spot rise over top-oil temperature. The standards present a set of correlations to
calculate these components based on some parameters with a related estimation table
for each. Transformers usually operate on a load cycle that repeats every 24 h. Hence,
the proposed models try to estimate the hot spot temperature of the transformer
considering load variation. Based on the proposed thermal model in, the hot spot
temperature is:
0<= = 0 - + ∆0#b + ∆0<= (2.18)
Where;
0<= Is winding hottest-spot temperature, (°C)
0 - Is average ambient temperature, (°C)
∆0#b Is top-oil rise over ambient temperature, (°C)
∆0<= Is winding hottest-spot rise over top-oil temperature, (°C)
The top-oil temperature rise is calculated by the following expression:
X
∆0#b = c∆0#b,d − ∆0#b, e W1 − S ghi ] + ∆0#b, (2.19)
Where;
∆0#b Is top-oil rise over ambient temperature, (°C)
∆0#b,d Is ultimate top-oil rise over ambient temperature, (°C)
∆0#b, Is initial top-oil rise over ambient temperature, (°C)
S Is the base of natural logarithm
K#b Is the oil time constant of transformer
The initial and ultimate top oil rise are calculated by:
kl7 35
∆0#b, = ∆0#b,3 × j 35
m (2.20)
32
kn7 35
∆0#b, = ∆0#b,3 × j m (2.21)
35
Where;
o Is the ratio of initial load to rated load, per unit.
o Is the ratio of ultimate load to rated load, per unit.
Is the ratio of load loss at rated load to no-load loss on the tap position.
The value of n has been chosen for each cooling model based on Table 2.1.
cooling transformer oil are its heat capacity and thermal conductivity. The electrical
and corresponding thermal laws of resistance and capacitance are as follows:
$
p= . ×R And R=`.× (2.22)
L
0= × q And q=` × (2.23)
G. Swift et al. present a simple equivalent circuit to represent the thermal heat flow
equations for power transformers. The proposed equivalent circuit is shown in Figure
2.18., which consists of two winding-to-oil and oil-to-ambient models. In order to
represent the effect of air or oil cooling convection currents, a nonlinear resistor is
assumed. Also, the heat generated by power losses and ambient temperature are
represented by current and voltage sources, respectively. This model considers three
inputs, including the iron loss (qfe), the copper loss (qcu), and the ambient temperature.
The oil thermal resistance (Roil) is a combination of the oil and the interface resistance
with the moving air layers. In addition, the hot spot resistance (Rhs) consists of the
insulation resistance and the interface resistance with the moving oil next to the
insulation. It should be noted that, for the lumped capacitance, it is assumed that the
temperature distribution within the medium (oil) is uniform at any given instant, which
is a reasonable assumption.
34
The differential equation related to Figure 2.20 for the oil-to-ambient model is:
H
Lslt
q +q =` . × +3 × v0 . −0 - wx (2.25)
slt,u
Where Coil is the thermal capacitance of the oil, and ORy, is the thermal resistance of
the oil at rated conditions. If we define as the ratio of qzU and q S at rated load, then
(2.25) reduces to:
7 35 H H
{|n Lslt
35
× }∆0 .,3 ~
x
=K . × + v0 . −0 - wx (2.26)
Where ∆0ORy, is the difference of 0ORy and 0•_€ at rated load, and KORy is the oil time
constant which can be obtained as follows:
K . = .,3 ×` . (2.27)
Hence, the difference equation corresponding to (2.26) is:
7 35 H H
‚ {|n
•0 . =M × ƒj 35
m × }∆0 .,3 ~
x
− v0 . −0 - wx „ (2.28)
slt
Where D is the difference operator, indicating a slight change in the associated variable.
As is shown in Figure 2.18, the winding-to-oil model is analogous to the oil-to-air
model. So, the equation for calculating the 0ℎP would be similar to (2.28).
In brief, network modeling reduces the complex pattern of the oil flow inside a
transformer winding down to a matrix of simple hydraulic channel approximations
interconnected by junction points or nodes. Figure 2.21, for example, shows the
geometry approximated for a 2D network model of windings of a transformer and its
corresponded THN model. The detailed THN model of windings comprises two
hydraulic network and thermal network models. The hydraulic part involves a network
of ducts and junctions describing the oil flow distribution. In this modeling, the oil flow
and pressure drop correspond to the electrical current and voltage, respectively. The
thermal model considering all types of heat transfer mechanisms is coupled with the
hydraulic models to evaluate the thermal performance of winding and obtain the HST
of the transformer.
36
enough, especially for describing the oil flow velocity. Therefore, the authors proposed
new equations based on CFD simulation data for the Nusselt number and friction
coefficient. In fact, although the thermal models based analytical methods can offer
acceptable accuracy in some cases, in most practical engineering applications, various
assumptions and simplifications need to be made to enable the analytical solution of
the differential equations representing the physical solution. This limits the
applicability of these methods to simple type problems or limits the validity of the
solutions if too many assumptions and simplifications are made. However, the advent
of digital computers has offered a solution to solve differential equations by numerical
methods, which can deal with these drawbacks. For several years’ significant effort has
been devoted to describing the thermal behavior of solid/fluid substances based on
numerical methods, resulting in computational fluid dynamics (CFD) modeling.
‡ˆ
+ ∇. (…† ) = 0 (2.29)
‡
Where † is the velocity vector of fluid (m/s) and ∇. Is known as the divergence operator
and operates as follows:
‡B ‡B ‡B
∇. " = + + (2.30)
‡ ‡ ‡‹
Where T is the pressure, Œ is the dynamic viscosity of the fluid, and A is external body
force.
39
The first law of thermodynamics states that the sources for the variation of the
total energy are the work of forces acting on the system plus the heat transmitted to this
system. Hence, with regard to the sources of energy variations in a fluid system and
according to Fourier’s law of heat conduction for transferring heat in a solid substance,
the conservative form of the energy equation is:
‡(ˆ•| #)
+ ∇. c…`a †'e = ∇. ( ∇') + Ž• (2.32)
‡
Inputs:
• Transformer load
• Top oil temperature
• Ambient temperature
• Tap position
• Cooling operation
Outputs:
• Dynamic rating
• Hot spot temperature
• Loss of life
40
The insulation paper that determines aging rates can be calculated in two
conditions where thermally upgraded can retain a higher percentage in tensile and
strength of insulation paper than non-thermally upgraded.
The DTR algorithm sends periodic reports to asset managers on the actual loss
of life of insulation paper. Further, this algorithm can calculate various pre loading and
overloading scenarios. Moreover, multiple technical challenges may occur in this
algorithm such as applying it to all power transformers in the network and the accuracy
of the calculated or measured transformer hot spot temperature. Many tests are applied
to above models help in obtaining precise results. Even analytical method works in
obtaining the probabilistic distribution of hot spot temperature of the transformer.
Monte Carlo technique is a developed simulation calculates complicity of a model
parameters for estimating the hot spot temperature from various loads. Monte Carlo
technique is used due to its mathematical simplicity and its ability to include more
probabilistic parameters for ambient temperature and transformer load. Further,
through the above technique, the lifetime of the transformer can be addressed.
Arrhenius chemical reaction rate law is used to calculate the relationship for thermal
aging properties of insulation material. This method states logarithm of time for
physical property of the insulation material. Therefore, the aging insulation is governed
by chemical process and its reactions vary with time. According to the Arrhenius
chemical reaction can be presented as in
š
yO™ 1 = "+ (2.34)
L‘›X 5 œ•
Several thermal models are used for estimating hot spot temperature behavior.
41
Start
Specification and Constraints, Specify max; and min; limits of Design Parameter
Is Analysis
Acceptable? No
Were Results
No Valid?
Reassess Design and Modeling
Yes
Conclusion and Discussion
References
[2] T. Wildi, “Electrical Machines & Drives,” Electrical Machines & Drives.
Prentice Hall, Columbus, 1985, doi: 10.1016/c2009-0-15969-x.
[6] "IEEE Recommended Practice for Performing Temperature Rise Tests on Oil-
Immersed Power Transformers at Loads Beyond Nameplate Ratings," IEEE Std
C57.119-2001, p. 0_1, 2002.
[8] J.-M. Mufuta and E. van den Bulck, “Modelling of the mixed convection in the
windings of a disc-type power transformer,” Appl. Therm. Eng., vol. 20, no. 5,
pp. 417–437, Apr. 2000, doi: 10.1016/S1359-4311(99)00034-4.