Volume 9 - Year 2022 - Pages 101-105
DOI: 10.11159/jffhmt.2022.013

Numerical Study of Electric Motors Cooling Using an Axial Air Flow

Ahmed M. Teamah1, Mohamed S. Hamed1

1 McMaster University, Department of Mechanical Engineering
Thermal Processing Laboratory
1280 Main Street West, Hamilton, Ontario, Canada L8S-4K1
teamaha@mcmaster.ca

Abstract - A numerical study has been carried out considering the heat transfer inside a four-pole synchronous electric motor. The computational work has been carried out using the commercial package ANSYS-CFX 2021 R1. The focus of this study is on the effect of an axial air flow passing through the gap between the stator and the rotor. The rate of cooling of the axial flow in terms of the average Nusselt number has been investigated at different rotational speeds and air flows. All surfaces were considered smooth. The source of heat generation was considered only within the rotor from the electrical windage resistance while the stator was considered insulated. The effect of the rotational speed and axial flow has been represented by a rotational and an axial Reynolds number which were varied in the ranges of 1750- 27000 and 2140- 6425, respectively. The numerical results have been validated using published experimental data with an acceptable deviation. Results showed that the average Nusselt number increases with both the rotational and the axial Reynolds numbers. However, the axial Reynolds number has a more dominant effect on the rate of cooling than the rotational Reynolds number.

Keywords: Heat transfer, Taylor Couette Poiseuille flow, Rotating concentric cylinders.

© Copyright 2022 Authors - This is an Open Access article published under the Creative Commons Attribution License terms Creative Commons Attribution License terms. Unrestricted use, distribution, and reproduction in any medium are permitted, provided the original work is properly cited.

Date Received: 2022-09-11
Date Accepted: 2022-09-15
Date Published: 2022-09-16

1. Introduction

Any electrical motor generates heat from the losses inside its electric windage. Hence, electric motors require proper cooling to avoid overheating. The cooling of a synchronous electric motor can be achieved by passing an axial air flow through the slotted gap between the rotor and the stator. Figure 1 shows the shape and location of such gab in a four-pole electrical motor. The flow rate of the axial air flow required to provide sufficient cooling depends on the rotational speed. Hence, the effect of the rotational speed and the rate of axial air flow on the rate of heat transfer inside the motor gap must be investigated.


Figure 1. Air gap inside the synchronous motors

2. Literature Review

Cooling of electric motors has been studied by many researchers as it is very critical in electric motors operation. Studies like [1-4] focused on the fluid flow and heat transfer of the Taylor-Couette flow inside electric motors, which represents flow between two concentric cylinders without an axial flow.

A review of flow between concentric cylinders was carried out by [5]. The review covered the presence of an axial flow with the Taylor-Couette flow. The problem of interest which involves a Taylor-Couette flow with an axial flow has been referred to as the Taylor-Couette-Poiseuille flow. The Taylor-Couette flow might experience flow instabilities, especially if the inner cylinder is rotating at high speeds while the outer cylinder is fixed. [6] investigated Taylor-Couette flow instabilities.

Other authors studied Taylor-Couette flow (TCF) considering various configurations. [7] considered a multiphase TCF. [8] studied TCF in a narrow spherical gap. [9] considered TCF in a wider gap. [10] considered a laminar-and turbulent TCF. [11] considered a highly turbulent TCF. [12] studied a TCF between two concentric cylinders while the inner cylinder is having ribs. [13] investigated heat transfer characteristics of a TCF in a gap while the outer cylinder has a set of axially distributed slits.

[14] studied the effect of the gab entrance region on TCF with inner rotating cylinder and fixed outer cylinder. They showed that the entrance region in their case is similar to the entrance region of a pipe or a duct. [15] added a set of slots to the motor stator and investigated their effect on the rate of heat transfer of the flow within the annular gap. They showed that these slots have almost no effect on the rate of heat transfer.

To the best of the authors’ knowledge, there has not been any studies that considered the effect of axial flow within a slotted rotor under operating conditions leading to a high rotational Reynolds number reaching 27000. These conditions are of prime interest in real practical applications.

3. Problem Definition

The problem of interest involves the cooling of a synchronous 4-pole electric motor using an axial air flow passing inside the gap between the rotor and the stator. The study resolves the flow and thermal fields inside the annular gap in the presence of the axial air flow through the rotor slots, as shown in Figure 2.

The inner rotating cylinder (i.e., the rotor) is subjected to a uniform heat flux of 500 W/m2, which is the amount of heat generated due to electric losses within the motor windage. The outer cylinder (i.e., the stator) is assumed here as an insulated surface, hence no heat is being transferred to the surrounding air. All the walls were assumed smooth. Thermal radiation has been neglected. The radius ratio, (Rri/Rro), of the considered configuration is 0.65, which is similar to commercial synchronous motor available in market. The length of the rotor is 0.5 m.


Figure 2. Schematic diagram of the present problem

4. Mathematical Model

The fluid flow and heat transfer problem has been simulated using ANSYS-CFX 2021 R1 commercial package. All the governing equations including mass, momentum and energy equations have been solved as referenced in [16].

The turbulence model used in this work is the Shear Stress Transport (SST) turbulence model to resolve the turbulence in the moving air flow. The SST turbulence model was chosen also in order to resolve the turbulent flow near the walls. The grid used in the current simulations has 106 elements. A grid dependence test has been carried out and the current grid has been selected so that the numerical results do not depend on the number of elements within 2%.

The boundary conditions are: at the insulated fixed stator, a zero heat flux is used. At the rotor: a fixed rotating speed (ω) is assumed. Its value was varied according to the value of the rotational Reynolds number (Rer) considered in each simulation. The rotor surface was assumed smooth and subjected to a constant uniformly distributed heat flux of 500 W/m2.

The air velocity at the entrance is constant and its value corresponds to the value of the axial Reynolds number (Rea). Entering air temperature is assumed 25℃. At the exit: air pressure is equal to the atmospheric pressure.

The definitions of the axial Reynolds number (Rea), the rotational Reynolds number (Rer), the Nusselt number (Nu) and the hydraulic diameter (Dh) are provided in the following equations:

(1)

(2)

(3)

(4)

The axial velocity (Va) is assumed equal to the average air velocity through the slots at the entrance. The kinematic viscosity (νair) and the thermal conductivity (Kair) of the air are calculated at the average fluid temperature across the gap. The heat transfer coefficient of air inside the gab is referred to by (h). The angular velocity of the rotor is identified by (ω). Gab width (P) is defined as shown in Figure 2.

4. 1. Validation of the Numerical Results

Validation of the current numerical results has been carried out using experimental results published in [17]. [17] considered heat transfer from an axial air flow inside a slotted rotor of a synchronous electrical motor.

The hydraulic diameter used in [17] was 15.8 mm, the radius ratio was kept fixed at 0.75 and the value of the axial Reynolds number was 4280. The value of the rotational Reynolds number was varied from 1750 to 10000. The heat generated inside the rotor is constant of 500 W/m2. The maximum deviation between the present numerical results and the experimental data reported in [17] is about 16.8%, as illustrated in Figure 3.


Figure 3. Comparison of current numerical results and experimental data reported in [17] at Rea = 4280

5. Results and Discussions

The effect of the parameters investigated in the present study on the rate of heat transfer within the motor gap are: the rotational Reynolds number (Rer) and the axial Reynolds number (Rea) in the range of 1750 to 27000 and 2140 to 6425, respectively.

1. Effect of Rer on Nuav

Figure 4 shows the effect of the rotational Reynolds number on the average Nusselt number inside the motor gap. The axial Reynolds number is kept constant at 2140. Results indicated that, increasing the rotational speed enhanced the rate of heat transfer within the gap. The results showed that the rate of increase in Nusselt number decreased somewhat at Rer about 8000. The slope in Figure 4 at Rer = 8000 has decreased by about 60% compared to the slope at Rer = 1750 (Green Line in Figure 4). The rate of increase in Nusselt number decreased after Rer = 8000 due to the increase in air circulation and the formation of large eddies inside the four longitudinal slots. Figure 5 shows the increase in air circulations and eddies inside one of the slots at Rer=27000 compared to that at Rer=6000.


Figure 4. Effect of rotational Reynolds number on average Nusselt number at Rea = 2140

Figure 5. Streamlines at different Rer at constant Rea of 2140

5.2. Effect of Rea on Nuav

The effect of the axial and rotational Reynolds number on the average Nusselt number is shown in Figure 6. Increasing the axial Reynolds number (Rea) from 2140 to 6425 resulted in an increase of the average Nusselt number inside the motor gap, at all values of Rer considered in this study. Also increasing the rotational Reynolds number at any value of the axial Reynolds number resulted in improvement in the heat transfer inside the gap.

However, results showed that the effect of Rea on the rate of heat transfer is more dominant than the effect of Rer. Therefore, motor cooling is significantly improved by increasing the axial air velocity through the gap between the rotor and the stator. It is clear from Figure 6 that the rate of enhancement in Nu is reduced at Rer equals to about 8000. Increasing Rer above 8000 resulted in stronger eddies of hot air forming inside the gap (i.e., the slots), which decreased the rate of enhancement in Nusselt number, as shown in Figure 6.


Figure 6. Effect of the Rer and Rea on Nuav.

5. Summary and Conclusions

The rate of cooling of a slotted electric motor using an axial air flow has been investigated numerically considering the effects of the rotational and the axial Reynolds numbers. Results showed that increasing the axial air flow or the motor rotational speed improves motor cooling. However, the effect of the axial air flow is more pronounced than the effect of the rotational speed. Results also showed that better cooling is achieved when the rotational speed corresponds to a rotational Reynolds number value below 8000, as the rate of enhancement in the average Nusselt number is decreased above this value.

6. Acknowledgement

The authors would like to acknowledge financial support of this work received from Mathematics of Information Technology and Complex Systems (Mitacs).


References

[1] A. Fasquelle et al., "Coupled electromagnetic acoustic and thermal-flow modeling of an induction motor of railway traction," Applied Thermal Engineering, vol. 30, pp. 2788-2795, 12/01 2010. View Article

[2] Y. Huai, R. Melnik, and P. Thøgersen, "Computational analysis of temperature rise phenomena in electric induction motors," Applied Thermal Engineering - APPL THERM ENG, vol. 23, 05/01 2003. View Article

[3] Z. Kolondzovski, A. Belahcen, and A. Arkkio, "Multiphysics thermal design of a high-speed permanent-magnet machine," Applied Thermal Engineering - APPL THERM ENG, vol. 29, pp. 2693-2700, 09/01 2009. View Article

[4] C.-H. Huang and H.-C. Lo, "A three-dimensional inverse problem in estimating the internal heat flux of housing for high speed motors," Applied Thermal Engineering - APPL THERM ENG, vol. 26, pp. 1515-1529, 10/01 2006. View Article

[5] M. Fénot, Y. Bertin, E. Dorignac, and G. Lalizel, "A review of heat transfer between concentric rotating cylinders with or without axial flow," International Journal of Thermal Sciences - INT J THERM SCI, vol. 50, pp. 1138-1155, 07/01 2011. View Article

[6] G. I. Taylor, "Stability of Viscous Liquid Contained Between Two Rotating Cylinders," Philosophical Transactions of the Royal Society of London, vol. 223, pp. 289-343, 01/01 1923. View Article

[7] D. Maillet, Y. Jarny, and D. Petit, "Problèmes inverses en diffusion thermique - Formulation et résolution du problème des moindres carrés," Ref: TIP201WEB - "Physique énergétique," Jul. 10, 2018. View Article

[8] S. Abbas and A. Shah, "Simulation of different flow regimes in a narrow-gap spherical Couette flow," Applied Mathematics and Computation, vol. 421, p. 126929, 2022/05/15/ 2022. View Article

[9] N. A. Manikandan, K. Pakshirajan, and G. Pugazhenthi, "A novel rotating wide gap annular bioreactor (Taylor-Couette type flow) for polyhydroxybutyrate production by Ralstonia eutropha using carob pod extract," Journal of Environmental Management, vol. 299, p. 113591, 2021/12/01/ 2021. View Article

[10] L. Jirkovsky and L. M. Bo-ot, "Laminar-turbulent transition in Taylor-Couette flow from a molecule dependent transport equation," Physics Letters A, vol. 408, p. 127481, 2021/08/27/ 2021. View Article

[11] G. Luo and Z. Yao, "Decoupling tests on axial heat-transfer in highly turbulent Taylor-Couette flow using thermal waves," Experimental Thermal and Fluid Science, vol. 128, p. 110439, 2021/10/01/ 2021. View Article

[12] M. Matsumoto et al., "Enzymatic starch hydrolysis performance of Taylor-Couette flow reactor with ribbed inner cylinder," Chemical Engineering Science, vol. 231, p. 116270, 2021/02/15/ 2021. View Article

[13] S.-l. Sun, D. Liu, Y.-Z. Wang, S. M. R. S. Naqvi, and H.-B. Kim, "Heat transfer characteristics of Taylor-Couette flow with axially distributed slits using field synergy principle and entropy generation analysis," International Communications in Heat and Mass Transfer, vol. 129, p. 105699, 2021/12/01/ 2021. View Article

[14] M. Molki, K. N. Astill, and E. Leal, "Convective heat-mass transfer in the entrance region of a concentric annulus having a rotating inner cylinder," International Journal of Heat and Fluid Flow, vol. 11, pp. 120-128, 06/01 1990. View Article

[15] C. Gazley, "Heat transfer characteristics of the rotational and axial flow between concentric cylinders," Trans. ASME, vol. 80, pp. 79-90, 01/01 1958. View Article

[16] C. J. S. t. g. ANSYS, ANSYS, "Release 12.0 User's guide," 2009.

[17] M. Fénot, E. Dorignac, A. Giret, and G. Lalizel, "Convective heat transfer in the entry region of an annular channel with slotted rotating inner cylinder," Applied Thermal Engineering, vol. 54, pp. 345-358, 05/14 2013. View Article