Advanced
Multi-physics Analysis for Temperature Rise Prediction of Power Transformer
Multi-physics Analysis for Temperature Rise Prediction of Power Transformer
Journal of Electrical Engineering and Technology. 2014. Jan, 9(1): 114-120
Copyright © 2014, The Korean Institute of Electrical Engineers
  • Received : March 01, 2013
  • Accepted : June 25, 2013
  • Published : January 01, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Hyun-Mo Ahn
Dept. of Electrical Engineering, Dong-A University, Korea. (hmahn@donga.ac.kr)
Joong-Kyoung Kim
R&D Center, Power & Industrial System Performance Group, HYOSUNG Corporation, Korea .(jkkim1976@hyosung.com)
Yeon-Ho Oh
Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute, Korea. (yhoh@keri.re.kr)
Ki-Dong Song
Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute, Korea. (yhoh@keri.re.kr)
Sung-Chin Hahn
Corresponding Author: Dept. of Electrical Engineering, Dong-A University, Korea .(schahn@dau.ac.kr)

Abstract
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.
Keywords
1. Introduction
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 [1 , 2] . 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 [3] .
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 [4 , 5] .
Previous work concerning the thermal dynamics of power transformers has been conducted by Marina A. Tsili and Eleftherios I. Amoiralis [4] , who derived a hybrid numerical-analytical technique for power transformer thermal modeling, and by Carlos Ortiz and Adam W. Skorek [5] , 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 [6 , 7] , 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 [6 , 7] . To couple thermal, electromagnetic, and fluid fields, an algorithm employing an UDF [8] is proposed.
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 [9 , 10] 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 [11] 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 Fig. 1(a) , consists of a reduce model of the transformer core and approximate windings, whereas the thermal-fluid component, shown in Fig. 1(b) , 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.
PPT Slide
Lager Image
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; Fig. 2 shows the multi-physics algorithm used to assess the material properties that are dependent on the temperature characteristics.
PPT Slide
Lager Image
Flow chart of multi-physics analysis
As can be seen in the Fig. 2 , 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 [1] . The amount of copper loss incurred can be calculated as a function of the main current running through the primary and secondary windings.
PPT Slide
Lager Image
where Pw is the power loss in the coil winding [W], σe is the electric conductivity [S/m],
PPT Slide
Lager Image
is the current density [A/m 2 ], and dv is the incremental volume element [m 3 ]. The electrical properties of copper are dependent on temperature, with the resistivity changing as temperature increases and the power loss correspondingly increasing:
PPT Slide
Lager Image
where αc is the temperature coefficient, ρ e,0 is the baseline electric resistivity [Ω·m] of copper, and T is the temperature [℃]. Thus, at 20 [℃], the calculated initial copper losses of HV and LV windings are 163,020.6 [W/m 3 ] and 188,305.7 [W/m 3 ], respectively.
Core losses can be grouped into two loss categories [1] : 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 [9 , 10 , 12] .
PPT Slide
Lager Image
where Pc is the core loss [W]; kh , kc and ke are the hysteresis loss coefficient, classical eddy current loss coefficient, and excess eddy current loss coefficient, respectively; f is the frequency [Hz]; and Bm 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 Table 1 ; in this case, the total core loss owing to the magnetic flux density is analytically calculated by using core loss coefficients, as shown in Fig. 3 .
Core Loss Coefficients @ 60Hz (Calculation)
PPT Slide
Lager Image
Core Loss Coefficients @ 60Hz (Calculation)
PPT Slide
Lager Image
Core losses due to the magnetic flux density @60Hz
The total core loss shown in Fig. 3 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 Fig. 4 .
PPT Slide
Lager Image
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 [4 , 8] :
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where,
PPT Slide
Lager Image
is the heat flow [W], k is the thermal conductivity [W/m·℃] of the medium, A is the cross section [m 2 ], ε is the surface emissivity, and σ is the Boltzmann constant [W/m 2 ·K 4 ].
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 [2] .
PPT Slide
Lager Image
PPT Slide
Lager Image
where, Qgen is heat source [J] and m 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 [11] .
PPT Slide
Lager Image
where, Nu is the Nusselt numbers, h is the heat transfer coefficient [W/m 2 ·℃] and δ 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 [5 , 13 , 14] :
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where, ρ is the density [kg/m 3 ],
PPT Slide
Lager Image
is the velocity [m/s],
PPT Slide
Lager Image
is the acceleration of gravity [m/s 2 ], μ is the viscosity [kg/m·s], p is the pressure [N/m 2 ], cp is the specific heat of the fluid at constant pressure [J/kg·℃], Φ is a dissipation function, and q′ is the energy per unit volume [W/m 3 ]. Turbulence phenomena in the cooling oil can be modeled using a k ε 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 [13] .
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
As shown in Fig. 5 , 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.
PPT Slide
Lager Image
Properties of cooling oil
The temperature distributions in a power transformer are shown in Fig. 6 : 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.
PPT Slide
Lager Image
Temperature distribution
As can be seen in Fig. 7 , 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.
PPT Slide
Lager Image
Comparison of temperature rise results
4. Conclusion
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.
Acknowledgements
This work was supported by the Dong-A University research fund.
BIO
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.
References
Kulkarni S. V. , Khaparde S. A. 2004 Transformer Engineering - Design and Practice Marcel Dekker, Inc.
Alegi G. L. , Black W. Z. 1990 “Real-time thermal model for an oil-immersed, forced-air cooled transformer” IEEE Trans. on Pow. Del. 5 (2) 991 - 999    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. 2009 “Hybrid numerical-analytic technique for power transformer thermal modeling” IEEE Trans. Magn. 45 (3) 1408 - 1411    DOI : 10.1109/TMAG.2009.2012647
Ortiz C. , Skorek A.W. , Lavoie M. , Benard P. 2009 “Parallel CFD analysis of conjugate heat transfer in dry-type transformer” IEEE Trans. Ind. Appl. 45 (4) 1530 - 1534    DOI : 10.1109/TIA.2009.2023561
Ahn H. M. , Kim M. S. , Song J. S. , Hahn S. C. 2011 “A study of coupled electromagnetic-thermal field analysis for temperature rise prediction of power transformer” Trans. KIEE 60 (10) 1838 - 1845    DOI : 10.5370/KIEE.2011.60.10.1838
Ahn H. M. , Lee B. J. , Hahn S. C. 2011 “An efficient investigation of coupled electromagnetic-thermal-fluid numerical model for temperature rise prediction of power transformer” in Proceedings of ICEMS 2011 Conference Beijing, China
2011 Fluent Product Documentation, Release 12.1 Ansys Inc.
Yamazaki K. , Fukushima N. 2010 “Iron-loss modeling for rotating machines: comparison between Bertotti's three-term expression and 3-D eddy-current analysis” IEEE Trans. Magn. 46 (8) 3121 - 3124    DOI : 10.1109/TMAG.2010.2044384
Bertotti G. 1988 “General properties of power losses in soft ferromagnetic materials” IEEE Trans. Magn. 24 (4) 621 - 630    DOI : 10.1109/20.43994
Kim J. K. , Hahn S. C. , Park K. Y. , Kim H. K. , Oh Y. H. 2005 “Temperature rise prediction of EHV GIS bus bar by coupled magneto-thermal finite element method” IEEE Trans. Magn. 41 (5) 1636 - 1639    DOI : 10.1109/TMAG.2005.846117
da Luz M. V. F. , Leite J. V. , Benabou A. , Sadowski N. 2010 “Three-phase transformer modeling using a vector hysteresis model and including the eddy current and the anomalous losses” IEEE Trans. Magn. 46 (8) 3201 - 3204    DOI : 10.1109/TMAG.2010.2049006
El Wakil N. , Chereches N.-C. , Padet J. 2006 “Numerical study of parameters affecting the temperature distribution in a disc-type transformer winding” J. Therm. Sci. 45 (6) 615 - 626    DOI : 10.1016/j.ijthermalsci.2005.09.002
Smolka J. , Nowak A. J. 2011 “Shape optimization of coils and cooling ducts in dry-type transformers using computational fluid dynamics and genetic algorithm” IEEE Trans. Magn. 47 (6) 1726 - 1731    DOI : 10.1109/TMAG.2011.2109731