DOI QR코드

DOI QR Code

Thermal Analysis of IPMSM with Water Cooling Jacket for Railway Vehicles

  • Park, Chan-Bae (Railroad Propulsion System Research Team, Korea Railroad Research Institute)
  • 투고 : 2013.09.24
  • 심사 : 2013.11.28
  • 발행 : 2014.05.01

초록

In this paper, the water cooling method among the forced coolant cooling methods is considered to be applied to the 110kW-class IPMSM for railway vehicles. First, basic thermal property analysis of the IPMSM is conducted using the three-dimensional thermal equivalent network method. Then, based on the results of the basic thermal property analysis, some design requirements for the water cooling jacket are deduced and a basic design of the water cooling jacket is carried out. Finally, thermal equivalent circuit of the water cooling jacket is attached to the IPMSM's 3D thermal equivalent network and then, the basic thermal and effectiveness analysis are conducted for the case of applying the water cooling jacket to the IPMSM. In the future, the thermal variation trends inside the IPMSM by the application of the water cooling jacket is expected to be quickly and easily predicted even at the design step of the railway traction motor.

키워드

Interior permanent magnet synchronous motor;IPMSM;Water cooling jacket;Thermal analysis;Thermal equivalent network

1. Introduction

The traction motor for the railway vehicle is installed on the axles of the bogie which is located on the lower part of the vehicle. Due to influence of the rail track environment, there is a large amount of dust blown during vehicle operations. Because of such operating environments, a totally-enclosed motor is primarily applied which can fundamentally block entry of any foreign substance from outside into the traction motor. However, the totally-enclosed motor, which has a problem for cooling, is not a suitable structure for the Interior Permanent Magnet Synchronous Motor (IPMSM) which has a permanent magnet embedded in the rotor. Especially, in recent market, where the downsizing of the traction motor both in size and weight is greatly in demand, research on a separate cooling device for the IPMSM is need to have a competitive edge.

In general, there are three cooling methods for the traction motor: natural air-cooling, forced air-cooling, and forced coolant cooling [1]. Of these, forced coolant cooling has the best cooling capacity. The smallest-sized traction motor is able to be designed if forced coolant cooling is applied to the IPMSM. However, for the application to be possible, research on the following is essential: 1) prediction of heat generated from inside the IPMSM which is designed to have high current density, 2) design of cooling channel and device to efficiently discharge the heat generated inside the IPMSM to the outside, and 3) a heat designed cooling channel and device [1]. Fig. 1 shows recent trends of R&D on the railway vehicle’s traction motor. Whereas there was research on the induction motor using natural air-cooling in the past, the direction of research has now changed into the IPMSM with water-cooling jacket. Therefore, in this paper, the water cooling method among the forced coolant cooling methods is considered to be applied to the 110kW-class IPMSM. First, basic thermal property analysis of the IPMSM is conducted using the three-dimensional (3D) thermal equivalent network method (TENM). Then, based on the results of the basic thermal property analysis, some design requirements for the water cooling jacket are deduced and a basic design of the water cooling jacket is carried out. Finally, thermal equivalent circuit of the water cooling jacket is attached to the IPMSM’s 3D thermal equivalent network (TEN). Then, the basic thermal and effectiveness analysis are conducted for the case of applying the water cooling jacket to the IPMSM.

Fig. 1.Recent R&D trend of traction motor for railway vehicles

Table 1.Main specifications of 110 kW-class IPMSM

 

2. Basic Thermal Property Analysis of the IPMSM

Thermal analysis is conducted for the 110kW-class IPMSM model by organizing 3D TEN. For the thermal analysis of the IPMSM, the overall heat distribution and flow of the IPMSM are identified by employing the lumped parameter method. When deducing heat sources such as thermal losses at each part of the heat equivalent circuit, a hybrid process using the distributed parameter method is applied [2]. Fig. 2 shows the thermal analysis concept of the IPMSM and 3D TENM proposed in this paper [3]. As shown in Fig. 2, 3D TEN considered a small control volume of the IPMSM is derived. Based on shapes and materials of each control volume of the IPMSM, the main thermal transfer coefficients (conduction/convection coefficients) are derived. The results are summarized in Table 2 [4]. Fig. 3 shows the transient thermal analysis results of the IPMSM at rated operation (2400rpm). As shown in Fig. 3, even with a relatively short time of operation, the temperature at each part of the IPMSM increases drastically. This is because the IPMSM is a totally-enclosed model without any cooling method, and the channel to forcefully discharge the heat, which has been generated inside the IPMSM, is not considered.

Fig. 2.Concept of thermal analysis of IPMSM and 3-D thermal equivalent network model

Table 2.Main thermal transfer coefficients of 110kW-class IPMSM

Fig. 3.Temperature rising properties of 110 kW-class IPMSM model (2400rpm)

 

3. Deduction of Specifications for a Water Cooling Jacket and Heat Transfer Properties

The cooling device for the IPMSM under consideration in this paper is a solenoid-typed water cooling jacket which covers the exterior of the stator core. According to the structure of the inner pipe, both 1-channel-type and 2-channel-type water cooling jackets are reviewed. Fig. 4 shows the structures of the water cooling jacket for the IPMSM. The following are the boundary conditions for analyzing the fluid movement and heat transfer properties inside the pipe of the water cooling jacket [5, 6]. i) TFor the moment, this paper only considers a circular shape as the cross-section form of the pipe for water cooling. In the fluid inside the pipe, Vm is the average speed of the fluid. While fluid is flowing inside the water cooling jacket pipe with length L, the pressure drop is as shown in Eq. (1). Here, f is the friction coefficient and D ishe surrounding temperature of the pipe inside the water cooling jacket is uniform. ii) The inner surface of the pipe is smooth. iii) The flow speed inside the pipe is uniform regardless of position.

Fig. 4.Water cooling jackets for 110kW-class IPMSM

For the moment, this paper only considers a circular shape as the cross-section form of the pipe for water cooling. In the fluid inside the pipe, Vm is the average speed of the fluid. While fluid is flowing inside the water cooling jacket pipe with length L, the pressure drop is as shown in Eq. (1). Here, f is the friction coefficient and D is the diameter of the pipe. The pump power required to overcome a certain pressure drop of ΔP can be calculated by using Eq. (2); here, V is the volumetric flow rate that passes through the pipe. The friction coefficient and the heat transfer coefficient are higher at the pipe entry where the thickness of the boundary layer is 0, and gradually decrease inside the pipe. Therefore, the pressure drop and heat flux are higher at the pipe entry and improve the average friction coefficient and heat transfer coefficient throughout the entire pipe due to entry area effect. However, the entry area effect may be overlooked for long pipes. In circular pipes, the average Nu is given as shown in Eq. (3) for laminar flows. Here, Nu is the Nusselt number, Re is the Reynolds number, Pr is the Prandtl number, and μb and μs are the dynamic viscosity of fluid calculated on the average temperature and the surface temperature, respectively. The interaction formulas between friction factor, f, and Nu are given in Equations (4) and (5) in case of the completely developed turbulent flows in the smooth circular pipe.

Table 3 shows the main specifications of the water cooling jacket for the 110kW-class IPMSM. The pipe material is copper with high heat conductivity. The total length of the pipe is 8.8(m) and both the 1-channel and 2-channel models are considered. Fig. 5 shows the analysis results of each characteristic parameter following changes in flow velocity of coolant inside the pipe of the water cooling jacket. Fig. 5 (a) shows changes to the Nusselt number. As shown in Fig. 5 (a), turbulence is generated even in slow flow velocity as the diameter of the pipe increases. On the other hand, if the pipe diameter is uniform regardless of pipe length, then the Nusselt number is consistent throughout. Fig. 5 (b) shows the exit temperature variations of the coolant. Since the surrounding temperature of the pipe is 75(℃), it would help in improving the thermal absorption of the coolant if the water cooling jacket is designed such that the exit temperature of the coolant would reach around 75(℃). For such a design, the case in which the flow inside the pipe is turbulent would be most effective, but the required thermal absorption of the pipe and the capacity of the coolant circulation pump would also have to be considered at the same time. Fig. 5 (c) shows changes in NTU. NTU is known as the net transfer unit and it is the index of a heat transfer system performance. In general, for a system with overcome the heat transfer threshold [7, 8]. As shown in Fig. 5 (c), the shorter the pipe length or the bigger the pipe diameter, it is easier to design a system with NTU over 5. Fig. 5 (d) shows changes in the thermal absorption rate of the coolant. As shown in Fig. 5 (d), the larger the pipe diameter, the higher the thermal absorption rate of the coolant, and using one pipe with a long channel yields a higher thermal absorption rate than using two pipes with short channels. Additionally, as the pipe diameter gets larger, the variation of the thermal absorption rate increases even with small changes in the flow velocity of the coolant. In areas where the coolant flow is turbulent, the thermal absorption rate is seen to drastically increase. Moreover, since the total sum of electromagnetic losses of the IPMSM, excluding mechanical loss, is about 5.83(kW), the thermal absorption of the water cooling jacket should be limited to 6(kW). Such limitation is also related to the pump capacity required for coolant circulation; thus, economic aspects should be considered as well [9, 10]. The design specifications of the water cooling jacket for the 110kW-class IPMSM are summarized in Table 4.

Table 3.Main specifications of Water-cooling Jacket

Fig. 5.Analysis results of characteristic parameters by flow velocity increasing on the inside of the pipe

Table 4.Design Results of Water-cooling Jacket for 110 kW-class IPMSM

 

4. Thermal Analysis of 110kW-class IPMSM with the Water Cooing Jacket

4.1 Thermal analysis method of the IPMSM with the water cooling jacket using TENM

In this paper, TENM is employed to analyze the thermal property variations on the IPMSM when the water cooling jacket is applied. More specifically, variations in temperature and current density at the IPMSM coil by heat absorption of the water cooling jacket were analyzed.

The thermal analysis method using TEN for the 110kW-class IPMSM with the water cooling jacket is shown in Fig. 6. As shown in Fig. 6, the coolant flowing through the pipe inside the water cooling jacket acts as an absorbent of the heat surrounding the pipe. Thus, the absorption rate can be equalized to a heat source on the TEN. However, unlike the previous IPMSM TEN in which the losses generated from the stator, rotor, permanent magnet, and coil are input as (+) heat sources, the water cooling pipe of the water cooling jacket should be input as a (−) heat source since the pipe absorbs surrounding heat. To analyze thermal properties of the IPMSM with the water cooling jacket using TENM, it must be assumed that the thermal absorption rate of the coolant flowing through the pipe inside the water cooling jacket is consistent throughout. In Table 4, the total thermal absorption capacity of the water cooling jacket has been limited to 6(kW). Therefore, as shown in Fig. 6, if the water cooling jacket is divided into four control volume models in axial directions, then the jacket can be equalized to four (−) heat sources. Consequently, the thermal absorption capacity of each control volume reaches a maximum of approximately 1.25 (kW).

Fig. 6.Analysis methods by thermal equivalent network when water-cooling jacket is applied to 110kWclass IPMSM

4.2 Thermal analysis results by applying the water cooling jacket

TEN of the 110kW-class IPMSM with the water cooling jacket is created by inserting the water cooling jacket model at the stator core exterior of the IPMSM. The flow velocity of the coolant is increased from 0 to 0.35(m/s) and the temperature variation property inside the IPMSM is predicted. First, as shown in Fig. 3, after 150 minutes of operation without the application of the water cooling jacket, the temperatures of the IPMSM coil, the frame and the permanent magnet reached about 180(℃), 121(℃), 187(℃), respectively. Fig. 7 shows the temperature variation effects on the IPMSM model due to the increase in the coolant flow velocity upon application of the water-cooling jacket. As shown in ℃. 7, the temperature at all parts inside the IPMSM drops as the flow velocity of the water cooling jacket is increased from 0 to 0.35(m/s). Especially, the temperatures at the stator and coil experience a much greater drop than the rotor since these parts are in direct contact with the water cooling jacket. Additionally, when the flow velocity exceeds 0.2(m/s), the flow inside the pipe changes from laminar to turbulent, and the thermal absorption rate and temperature drop rate increase due to coolant. Fig. 8 shows the property variations on the coil of the IPMSM model due to the increase in coolant flow velocity on the water cooling jacket. Fig. 8 (a) shows the property variation between the thermal absorption rate of the water cooling jacket and the temperature of the coil due to the increase of the coolant flow velocity. When the flow velocity is increased up to 0.35(m/s), the coil temperature drops from 180(℃) to 95.1(℃). This is approximately an 85(℃) drop. Fig. 8 (b) and (c) show the property variation among the thermal absorption rate, phase resistance and phase current on the coil. When the flow velocity is increased up to 0.35(m/s), the phase resistance drops from 5.95(mΩ) at 180(℃) to 4.78(mΩ) by approximately 1.17(mΩ) by the temperature drop of the coil. Consequently, the phase current margin of about 35.4(A) from 307(A) at 180(℃) to 342.4(A) can be obtained. Fig. 8 (d) shows the property variation between the thermal absorption rate and the current density of the coil. When the flow velocity is increased up to 0.35(m/s), it can be predicted that the current density of the coil is increased by approximately 0.52(A/mm2) from 4.5(A/mm2) to 5.02(A/mm2). According to the analysis results, the phase resistance decrease of approximately 32(%) is generated because of the temperature drop of the coil caused by the thermal absorption of the water cooling jacket. The decrease in the phase resistance in turn leads to an increase in the coil current density by approximately 21(%).

Fig. 7.Temperature variation properties of 110kW-class IPMSM model due to the increase of the coolant flow velocity on the water-cooling jacket

Fig. 8.Properties variations on the winding of 110kW IPMSM-class model due to the increase of the coolant flow velocity on the water-cooling jacket

In conclusion, the thermal analysis method of the 110 kW-class IPMSM model with the thermal model of the water cooling jacket using the 3D TENM is presented in this paper, and the temperature variation trends inside the IPMSM by the application of the water cooling jacket can be quickly and easily predicted even at the design step of the railway traction motor.

 

5. Conclusion

Water cooling jackets have yet to be applied to any traction motors for railway vehicles in Korea. Thus, this paper handled to conduct basic design of a water cooling jackets and also investigate heat transfer properties. The 3D TENM for the 110 kW-class IPMSM is first employed to conduct a basic thermal property analysis. After analyzing the basic thermal properties of the IPMSM, the necessary specifications for the water cooling jacket were deduced based on the thermal analysis results of the IPMSM. Next, heat transfer properties were analyzed by using an equivalent model of the solenoid-typed pipe which is to be installed within the water cooling jacket of the IPMSM. A basic design model for the 6 kW-class water cooling jacket was deduced. Finally, thermal equivalent circuit of the water cooling jacket was attached to the IPMSM’s 3D TEN, and the basic thermal and effectiveness analysis were conducted for the case of applying the water cooling jacket to the IPMSM. In future studies authors intend to make a 6 kW-class water cooling jacket and a 110 kW-class IPMSM and to validate the proposed design and thermal analysis method of the water cooling jacket with the IPMSM.

참고문헌

  1. A. F. Armor, and M. V. K. Chari, "Heat flow in the stator core of large turbine-generators by the method of three dimensional Finite Elements (Part II: Temperature distribution in the stator iron)," IEEE Transactions on PAS, PAS-95(5), pp. 1657-1668, 1976.
  2. Y. Bayazitoglu and M. N. Ozisik, "Element of Heat Transfer." New York : McGraw-Hill, 1988.
  3. C. B Park, H.W. Lee, B. S. Lee, "A Study on the Reduction of the Stator Iron Loss on Permanent Magnet Synchronous Motor for light Railway Transit Propulsion System," Journal of the Korean Society for Railway, Vol.15, No.4, pp. 376-380, 2012. https://doi.org/10.7782/JKSR.2012.15.4.376
  4. C. B Park, H.W. Lee, B. S. Lee, "Investigation of a thermal analysis method for IPMSM in railway vehicles," Journal of the Korean Society for Railway, Vol.16, No.2, pp. 99-103, 2013. https://doi.org/10.7782/JKSR.2013.16.2.099
  5. C. B. Park, "A Study on the Thermal Characteristics of 110kW-class IPMSM for Light Railway Transit using the 3-Dimensional Thermal Equivalent Network considering Heat Source by Iron Loss Density Distributions," The transactions of the Korean Institute of Electrical Engineers, Vol. 62, No. 7, pp. 1038-1044, 2013. https://doi.org/10.5370/KIEE.2013.62.7.1038
  6. F. Kreith and M.S. Bohn, "Principles of Heat Transfer: 5th edition," St. Paul, MN : West Publishing, 1993.
  7. T. Jokinen and J. Saari, "Modeling of the coolant flow with heat flow controlled temperature sources in thermal networks," IEE Proceeding of Electrical Power Applications, Vol.144, No.5, 1997.
  8. J.H. Choi, S.M. Jang, C. Han, J.H. Ahn, K.H. Sim, "Electromagnetic loss and thermal analysis on high speed PM motor for turbo compressor," KIEE Summer Conference 2012, pp. 618-619, 2012.
  9. J. Driesen, R. J. M. Belmans, and K. Hameyer, "Finite-element modeling of thermal contact resistances and insulation layers in electrical machines," IEEE Transactions on Industry Applications, Vol.37, No.1, pp. 15-20, 2001. https://doi.org/10.1109/28.903121
  10. Batchelor, G. K., "An Introduction to Fluid Dynamics," Cambridge University Press, 1967.

피인용 문헌

  1. The Equilibrium Design of Radial Magnetic Force for Reduction of Vibration in IPM type BLDC Motor vol.11, pp.2, 2016, https://doi.org/10.5370/JEET.2016.11.2.377
  2. Design criteria for water cooled systems of induction machines vol.114, 2017, https://doi.org/10.1016/j.applthermaleng.2016.12.031
  3. Wafter design for totally enclosed electric machines vol.129, 2018, https://doi.org/10.1016/j.applthermaleng.2017.10.011