In this paper, a method for multi-physics analysis of the temperature-dependent properties of an oil-immersed transformer is discussed. To couple thermal fields with electromagnetic and fluid fields, an algorithm employing a user defined function (UDF) is proposed. Using electromagnetic analysis, electric power loss dependent on temperature rise is calculated; these are used as input data for multi-physics analysis in order to predict the temperature rise. A heat transfer coefficient is applied only at the outermost boundary between transformer and the atmosphere in order to reduce the analysis region. To verify the validity of the proposed method, the predicted temperature rises in high-voltage (HV) and low-voltage (LV) windings and radiators were compared with the experimental values.
Magnetic circuits and windings are principal sources of loss, and thus of temperature rise, in a transformer. Core and copper losses are mainly responsible for heat generation, and such heat must be dissipated in order to keep the core and winding parts from reaching a temperature that will cause deterioration of the insulation
. As maximum allowable temperatures depend on factors such as the ambient temperature, transformer design, loading condition, and amount of cooling provided, the accurate estimation of the transformer temperature is critical design factor, and limits for temperature rise are specified in International Standard
In general, temperature rise can be predicted by means of a heat transfer analysis using the finite element method (FEM) or by mean of a heat fluid analysis using the finite volume method (FVM). Of these two, heat transfer analysis has the advantage of requiring less computing time, but has poor accuracy for circuits with complex geometry. Thermal fluid analysis, on the other hand, is more accurate over a wide range of circuit geometries, but requires a longer computing time; thus, in order to improve the fluid transfer analysis, it is necessary to develop new methods for reducing the computing time
Previous work concerning the thermal dynamics of power transformers has been conducted by Marina A. Tsili and Eleftherios I. Amoiralis
, who derived a hybrid numerical-analytical technique for power transformer thermal modeling, and by Carlos Ortiz and Adam W. Skorek
, who performed a parallel computational fluid dynamics (CFD) analysis of conjugate heat transfer in a dry-type transformer. Although these studies produced some interesting results, they mainly investigated the thermal problem itself, as opposed to the interactions of the physical processes of the thermal, electromagnetic, or fluid fields; as such, accurate temperature rise prediction was not attempted.
In our previous works
, we dealt with coupled problem such as electromagnetic-thermal field of oil-immersed power transformer. We mainly investigated coupled electromagnetic-thermal problem.
Now, in this paper, a multi-physics analysis that can produce an accurate temperature rise prediction in an oil-immersed transformer is presented to obtain the more accurate results than our previous works
. To couple thermal, electromagnetic, and fluid fields, an algorithm employing an UDF
In order to accurately predict the temperature rise, accurate modeling of heat sources is needed. Electric power losses can be calculated by means of electromagnetic analysis; in order to do so accurately, the copper loss can be considered to be temperature-dependent, whereas the core loss can be accurately estimated using Bertotti’s model
as derived from Steinmetz’s equation. Electric power losses calculated in this way can then be used as input data for multi-physics analysis; in order to predict the temperature rise accurately, the material properties of cooling oil can be considered to be temperature-dependent. Also, we apply heat transfer coefficient obtained by calculating the Nusselt number
at only the outermost boundary between the transformer and the atmosphere this restriction is applied in order to reduce the computing time. To verify the validity of this proposed analysis method, its results in terms of temperature rise in the high-voltage (HV) and low-voltage (LV) windings and radiator can be compared with the experimental results.
2. Multi-Physics Analysis Model
The transformer studied was an oil-immersed, single-phase unit (1MVA, 6600/690V, 60Hz). The multi-physics analysis model used to assess power 1oss in this transformer has an electromagnetic and a thermal-fluid component; to reduce the computing time, the electromagnetic component, shown in
, consists of a reduce model of the transformer core and approximate windings, whereas the thermal-fluid component, shown in
, consists of a model of the transformer parts. To reduce the analysis region, we applied a quarter symmetry model to determine boundary conditions such as the symmetry boundary and the heat transfer coefficient.
Multi-physics analysis models
3. Multi-Physics Analysis
- 3.1 Multi-physics analysis method
In order to obtain an accurate analysis of characteristics of the transformer system, an assessment of the relationships between multi-physical phenomena (multi-physics) is required. For the study described in this paper, we employed a user defined function (UDF) algorithm to solve the multi-physical problem of coupled thermal, electromagnetic and fluid fields. Some of the material properties relevant to both the electromagnetic and the fluid fields are dependent on temperature and thus on the thermal effects;
shows the multi-physics algorithm used to assess the material properties that are dependent on the temperature characteristics.
Flow chart of multi-physics analysis
As can be seen in the
, the process involves calculating electric power losses using electromagnetic analysis and then predicting temperature rise using thermal-fluid analysis, with this sequence of analysis iteratively processed until the temperature reaches a steady state.
- 3.2 Analysis of electric power loss (Heat Source)
Two main types of losses can occur in an electric power transformer: copper loss and core loss
. The amount of copper loss incurred can be calculated as a function of the main current running through the primary and secondary windings.
is the power loss in the coil winding [W],
is the electric conductivity [S/m],
is the current density [A/m
is the incremental volume element [m
]. The electrical properties of copper are dependent on temperature, with the resistivity changing as temperature increases and the power loss correspondingly increasing:
is the temperature coefficient,
is the baseline electric resistivity [Ω·m] of copper, and
is the temperature [℃]. Thus, at 20 [℃], the calculated initial copper losses of HV and LV windings are 163,020.6 [W/m
] and 188,305.7 [W/m
Core losses can be grouped into two loss categories
: hysteresis loss and eddy current loss. Hysteresis loss, which varies as a function of frequency, can be defined as the area within a static magnetization loop. Eddy current losses can be further subdivided into classical eddy current losses and excess eddy current losses. The classical eddy current loss can be derived from Maxwell’s equation and is proportional to the square of the frequency, whereas the excess eddy current loss is a residual core loss mechanism (i.e., one that is not described by hysteresis loss or classical eddy current loss); it is mathematically defined as the total core loss minus the hysteresis and classical eddy current losses. Total core loss can then be separated using Bertotti’s model
is the core loss [W];
are the hysteresis loss coefficient, classical eddy current loss coefficient, and excess eddy current loss coefficient, respectively;
is the frequency [Hz]; and
denotes the maximum flux density [T]. Bertotti’s model can be used to mitigate errors arising from the use of Steinmetz’s equation, which is typically less accurate over a maximum operating magnetic flux density of 1 [T]. In order to solve Bertotti’s model, the three core loss coefficients can be analytically calculated by using core loss data (the B-P curve) shown in
; in this case, the total core loss owing to the magnetic flux density is analytically calculated by using core loss coefficients, as shown in
Core Loss Coefficients @ 60Hz (Calculation)
Core Loss Coefficients @ 60Hz (Calculation)
Core losses due to the magnetic flux density @60Hz
The total core loss shown in
is the sum of the hysteresis loss, the classical eddy current loss, and the excess eddy current loss, and there is a good agreement between the measured loss and the calculated total core loss values over the operating range of a power transformer. The magnetic flux density and core loss distributions in the transformer core can be numerically calculated by using an FEM for a transformer’s magnetic field, as shown in
Magnetic field analysis results of a transformer core
Here, it can be seen that the average magnitude of loss at the transformer core is 1,357.1 [W]; this represents an error of approximately 7.9 [%] compared to the design value. The electric power losses generated by combining the initial copper and core losses in this method can then be used as input data for a multi-physics analysis.
- 3.3 Temperature ri analysis
A heat transfer mechanism requires the existence of a temperature difference; additionally, in all modes of heat transfer, movement of heat proceeds from a higher to a lower temperature medium. Heat transfer in a transformer is effected through conduction, convection, and radiation, given respectively as below
is the heat flow [W],
is the thermal conductivity [W/m·℃] of the medium,
is the cross section [m
is the surface emissivity, and
is the Boltzmann constant [W/m
Heat transfer in an oil immersed transformer mainly occurs by circulation of cooling oil. Heat generated by electric power losses is transmitted to tank and radiator. According to law of energy conservation, the energy balance equations in insulation oil are given as below
is heat source [J] and
is mass [kg]. The heat transfer coefficient can then be defined on the solidfluid boundary (i.e., the tank and the radiator). In order to reduce computing time, the boundary condition defined average heat transfer coefficient can be obtained practically by applying Nusselt number
is the Nusselt numbers,
is the heat transfer coefficient [W/m
is the characteristic length [m].
Typically, the heat transfer in a power transformer is primarily influenced by the circulation of the cooling oil, the fluid dynamics of which can be modeled using CFD. Fluid flow and temperature distribution within a transformer are governed by three conservation equations. These can be solved simultaneously with a semi-implicit method for pressure-linked equations (SIMPLE) algorithm, which uses a combination of continuity and momentum equations in order to derive an equation for pressure correction
is the density [kg/m
is the velocity [m/s],
is the acceleration of gravity [m/s
is the viscosity [kg/m·s],
is the pressure [N/m
is the specific heat of the fluid at constant pressure [J/kg·℃], Φ is a dissipation function, and
is the energy per unit volume [W/m
]. Turbulence phenomena in the cooling oil can be modeled using a
turbulence model having two transport equations one dealing with the kinetic energy of turbulence and the other involving the dissipation rate.
Many properties of the cooling oil are affected by temperature as below
As shown in
, it can be seen that the specific heat increases linearly with temperature, whereas the density and thermal conductivity decrease linearly with temperature and the dynamic viscosity decreases in a nonlinear manner. In fact, except for viscosity, all the cooling oil properties change linearly with temperature.
Properties of cooling oil
The temperature distributions in a power transformer are shown in
: sequential and multi-physics analyses show maximum temperatures of 84.9 [℃] and 111.7 [℃], respectively, at the upper part of the LV winding. This maximum temperature at the upper part is attributed to the fact that while heat flow in the duct region of the inner winding is reduced by the mechanical and insulation structure, and thermal energy in the cooling oil moves upwards owing to buoyancy. Furthermore, as the amount of heat generated in the LV winding is greater than that in the HV winding, the temperature in the region of the LV winding is higher than that in the region of the HV winding. A comparison of the results of the temperature distribution by both the sequential and multi-physics analysis models shows a maximum temperature differential between the two of 26.8 [℃] because heat sources are increased by temperature rise. In order to verify the validity of multi-physics analysis, the temperatures rise at certain points on the HV and the LV windings as well as the radiator were measured using a temperature rise test and then compared with the numerically predicted results.
As can be seen in
, the amount of temperature rise predicted by the multi-physics analysis model showed greater agreement with the experiment than that predicted by the sequential model.
Comparison of temperature rise results
In this paper, a multi-physics analysis used to predict the temperature rise in an oil-immersed transformer was described. In this method, an algorithm using a UDF is performed in order to couple thermal, electromagnetic field, and fluid fields. To accurately predict the temperature rise, accurate modeling of heat sources is required. In order to calculate heat source parameters accurately, copper loss is considered to be temperature dependent, and core loss is calculated using Bertotti’s model. The numerically calculated magnetic flux density values and the analytically calculated core loss coefficients could be used to solve this model. The results from this method agreed well with the design values.
To reduce the region of analysis, we used a quarter symmetry model determined by the boundary conditions such as the symmetry boundary and the heat transfer coefficients. The electric power losses and boundary conditions thus calculated were then used as input data for multi-physics analysis; in order to predict the temperature rise accurately, the material properties of cooling oil considered to be temperature-dependent. In order to verify the validity of the analysis, the predicted temperature rises in the HV and LV windings and radiator were compared with the experimental values. From this, it could be seen that the temperature rise predicted by the multi-physics analysis is closely in accord with the experimentally obtained values; this gives us confidence that the method proposed in this paper will be useful in improving the thermal design of power transformers.
This work was supported by the Dong-A University research fund.
Hyun-Mo Ahn He received the B.S. and M.S. degrees in electrical engineering from Dong-A University, Busan, Korea, in 2009 and 2011, respectively. He is currently working toward the Ph.D. degree. His research interests are multi-physics analysis and numerical analysis of power apparatus and electric machines.
Joong-Kyoung Kim He received the B.S., M.S., and Ph.D. degrees in electrical engineering from Dong-A University, Busan, Korea, in 2002, 2004 and 2008, respectively. He is currently a Chief Researcher with the R&D Center, Power & Industrial Systems Performance Group, HYOSUNG Corporation.
Yeon-Ho Oh He received the B.S. and M.S. degrees in electrical engineering from Dong-A University, Busan, Korea, in 1991 and 1993, respectively. He is currently a Research Engineer with the Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute (KERI).
Ki-Dong Song He received the B.S. and M.S. degrees in electrical engineering from InHa University, Incheon, Korea, in 1988 and 1990, respectively. He got his Ph.D. degree in electrical engineering from Busan National University in 2003. He worked for Korea Electrotechnology Research Institute (KERI) from 1990 and is now the director of the Power Apparatus Research Center, HVDC Research Division. His research interests focus on analysis and design for AC Power Devices and HVDC circuit breakers.
Sung-Chin Hahn He received the B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National University, Seoul, Korea, in 1979, 1981, and 1992, respectively. He is currently a Professor with the Department of Electrical Engineering, Dong-A University, Busan, Korea. His current research covers multi-physics analysis and optimal design of power apparatus and electric machines.
Kulkarni S. V.
Khaparde S. A.
Transformer Engineering - Design and Practice
Marcel Dekker, Inc.
Alegi G. L.
Black W. Z.
“Real-time thermal model for an oil-immersed, forced-air cooled transformer”
IEEE Trans. on Pow. Del.
DOI : 10.1109/61.53112
IEEE C57.91-1995-IEEE Guide for Loading Mineral-Oil-Immersed Transformers
Tsili Marina A.
Amoiralis Eleftherios I.
Kladas Antonios G.
Souflaris Athanassios T.
“Hybrid numerical-analytic technique for power transformer thermal modeling”
IEEE Trans. Magn.
DOI : 10.1109/TMAG.2009.2012647
“Parallel CFD analysis of conjugate heat transfer in dry-type transformer”
IEEE Trans. Ind. Appl.
DOI : 10.1109/TIA.2009.2023561
Ahn H. M.
Kim M. S.
Song J. S.
Hahn S. C.
“A study of coupled electromagnetic-thermal field analysis for temperature rise prediction of power transformer”
DOI : 10.5370/KIEE.2011.60.10.1838
Ahn H. M.
Lee B. J.
Hahn S. C.
“An efficient investigation of coupled electromagnetic-thermal-fluid numerical model for temperature rise prediction of power transformer”
in Proceedings of ICEMS 2011 Conference
Fluent Product Documentation, Release 12.1
“Iron-loss modeling for rotating machines: comparison between Bertotti's three-term expression and 3-D eddy-current analysis”
IEEE Trans. Magn.
DOI : 10.1109/TMAG.2010.2044384
“General properties of power losses in soft ferromagnetic materials”
IEEE Trans. Magn.
DOI : 10.1109/20.43994
Kim J. K.
Hahn S. C.
Park K. Y.
Kim H. K.
Oh Y. H.
“Temperature rise prediction of EHV GIS bus bar by coupled magneto-thermal finite element method”
IEEE Trans. Magn.
DOI : 10.1109/TMAG.2005.846117
da Luz M. V. F.
Leite J. V.
“Three-phase transformer modeling using a vector hysteresis model and including the eddy current and the anomalous losses”
IEEE Trans. Magn.
DOI : 10.1109/TMAG.2010.2049006
El Wakil N.
“Numerical study of parameters affecting the temperature distribution in a disc-type transformer winding”
J. Therm. Sci.
DOI : 10.1016/j.ijthermalsci.2005.09.002
Nowak A. J.
“Shape optimization of coils and cooling ducts in dry-type transformers using computational fluid dynamics and genetic algorithm”
IEEE Trans. Magn.
DOI : 10.1109/TMAG.2011.2109731