## Abstract

To enhance worldwide environmental conditions, the air transport industry must drastically reduce carbon dioxide emissions. Electrification of aircraft propulsion systems is one way to meet this demand. In particular, the focus is on obtaining single-aisle aircraft with partial turbo-electric propulsion and approximately 150 passenger seats by the 2030s. To develop a single-aisle aircraft with partial turbo-electric propulsion, an air-cooled interior permanent magnet (IPM) motor with an output of 2 MW is desired. One of the most difficult problems in air cooling is that air-cooling performance decreases with increasing altitude because the air density decreases. To investigate the effect of altitude on air-cooling performance in the IPM motor, the authors formulated mathematical system equations to describe heat transfer inside the target air-cooled IPM motor, and mathematical analytical solutions were obtained. The most severe condition is the top-of-climb condition. For this condition, a designer should choose cooling air mass flow rates that keep the temperature of the permanent magnets below the maximum temperature limit of 100 °C and the temperature of the coils below the maximum temperature limit of 250 °C. Here, the sizes of the air-cooling channels strongly affect air cooling with the IPM motor. In this paper, the authors briefly review the mathematical formulations and their solutions, investigate the effect of channel size on air-cooling performance in an IPM motor, and explore the optimum configuration and settings for the air cooling channels.

## 1 Introduction

The Air Transport Action Group published a vision “Waypoint 2050X” [1], which stated the commercial aviation industry plans to decarbonize and achieve worldwide complete net-zero emissions by 2050. The International Air Transport Association including aircraft manufacturers, engine manufacturers, airports, and air traffic management, declared adoption of this goal [2,3] to support the United Nations Climate Change Conference (COP26). One of the innovations to decarbonize technology is electrification of aircraft propulsion systems. Felder at the National Aeronautics and Space Administration presented six promising types of aircraft electric propulsion systems [4].

Type 1: All-electric, in which batteries drive electric motor-powered fans.

Type 2: Parallel hybrid, in which battery-powered electric motors assist turbofans.

Type 3: Turbo-electric, in which gas turbines generate electricity and their electricity drives electric motor-powered fans.

Type 4: Series hybrid, in which gas turbines generate electricity, batteries provide electricity, and both electricity sources drive electric motor-powered fans.

Type 5: Partial turbo-electric, in which turbofans generate electricity and their electricity drives additional electric motor-powered fans.

Type 6: Series/partial hybrid, in which turbofans generate electricity, batteries provide electricity, and both sources drive additional electric motor-powered fans.

Currently, many companies are developing short-range aircraft with all-electric propulsion systems, such as those in Refs. [5–7], because all-electric systems are the simplest propulsion devices. The power density of cutting-edge batteries was 400 Wh/kg in 2021 [8], which is less than the 12,000 Wh/kg of kerosene. Thus, all-electric propulsion devices can be used only for short-range aircraft due to weight limitations for flight; in 2021, there are no aircraft with all-electric propulsion flying regional or longer ranges. Even by the 2030s, it would be difficult to realize single-aisle aircraft with all-electric propulsion and approximately 150 passenger seats.

However, single-aisle aircraft with partial turbo-electric propulsion will be possible in the 2030s, so there are many concepts for developing single-aisle aircraft with partial turbo-electric propulsion by the 2030s. One of the promising single-aisle aircraft with partial turbo-electric propulsion is the STARC-ABL [9]. Figure 1 shows the targeted single-aisle tube- and wing-configuration with attached generators and with a ducted, electrically driven, boundary layer-ingesting tailcone fan. Boundary layer ingestion realizes a high propulsive efficiency *η*_{pro}. Here, an electric motor used to drive the tailcone fan must output approximately 2 MW. In addition, this electric motor must be a lightweight and high-power density (= output/mass, W/kg) motor to enable flight. Furthermore, the electric motor must be air-cooled because there is no oil-lubricated gas turbine engine near the tail, unlike cases with attached generators. We have two choices in selecting the electric motor: superconductive motors and normal conductive motors. Although a superconductive motor is extremely lightweight and can operate at a gradually higher temperature than before [10], cryogenic cooling is required, and air cooling is not sufficient. Thus, superconductive motors look promising, but unfortunately, they are still a future technology at this time in the early 2020s. Among normal conductive motors, an interior permanent magnet (IPM) motor is presumed to be the best choice because of its simpler structure, higher efficiency, and higher power density.

Here, when the temperature of a permanent magnet exceeds its own maximum temperature limit, its magnetic flux density gradually decreases and does not recover. It is called “thermal demagnetization.” Aircraft fly at a high altitude of 13,000 m (42,650 ft), where the ambient pressure decreases to 16 kPa (2.3 psi) in an intercontinental cruise, and aircraft can take off from an airport in a desert during summer at a high temperature of 50 °C (122 °F). Under any conditions, a reliable design method for thermal management of an air-cooled electric motor is essential to prevent permanent demagnetization.

Many researchers have studied thermal management of electric motors with experimental approaches [11], with numerical approaches [12,13], and with both [14]. These approaches are very important and effective, but they are overly time-consuming for the early design phase; instead, theoretical approaches are very convenient and provide results quickly with changing design parameters.

In early designs of electric motors, thermal analyses using a thermal circuit model were most common, such as in Ref. [15]. On the other hand, in early designs for heat exchangers, several differential equations were formulated for heat balancing with all fluid flows, this system of differential equations was mathematically integrated, and the results (i.e., temperature distributions, heat flow rates, pressure losses, etc.) were obtained as in Refs. [16–20]. Tao et al. also applied the latter approach to perform thermal management of entire electric automobiles, including an electric motor, but they used matlab/simulink and amesim to solve this system of differential equations [21].

The authors and colleagues have advanced the latter approach using differential equations to analyze an air-cooled electric IPM motor driving an aviation fan and obtained transient solutions of the mathematical analytical equations instead of numerical solutions [22]; the mathematical analytical equations instantaneously derive the solution values without iterations, unlike numerical solutions. Then, the authors and colleagues obtained steady solutions with numerical time integration of the transient solutions. These solutions can be useful in allowing engineers to design air-cooled IPM motors early in design because steady solutions can be obtained quickly by using Microsoft Excel on a personal computer. Using this method, the authors discussed the capability and limitations of air cooling for the target IPM motor. From the results, the following predictions were made [22]: (1) internal recirculating air and external air should be separated inside and outside of the motor casing so electric parts in the casing also remain safe. (2) The internal air should move heat from the rotor to the outer heat exchanger, the heat exchanger should collect heat from the stator, and the external air should move all heat from the heat exchanger to the surroundings. (3) In the heat exchanger, the external air should flow between a stator and the internal air, and the flow direction of the external air should be opposite that of the internal air. (4) The internal air mass flowrate should be chosen to maintain turbulent internal flow through the heat exchanger. (5) There is the optimum flowrate of the internal air to minimize the rotor temperature. (6) The rotor temperature increases with decreasing internal air mass flowrate below the optimum value. (7) The rotor temperature increases with increasing internal air mass flowrate beyond the optimum value. (8) The rotor and stator temperatures are linearly proportional to the ambient temperature. (9) The rotor temperature increases with decreasing ambient pressure. (10) An increase in altitude causes a decrease in power consumption for both external air and internal air. (11) The most severe condition is the top-of-climb condition, and external and internal air mass flow rates should be chosen to keep the temperature of the permanent magnets below the maximum temperature limit of 100 °C and the temperature of the coils below the maximum temperature limit of 250 °C.

The above predictions were obtained for constant air-cooling channel sizes; however, the authors noticed that the sizes of air-cooling channel strongly affected air-cooling performance in the IPM motor. A designer has several choices for air-cooling channel sizes in the heat exchanger for internal and external air and choices for inner air-cooling channel sizes in the rotor. In this paper, the authors first briefly review the mathematical formulations and their solutions, then investigate the effect of channel size on air-cooling performance in the IPM motor, and finally explore the optimum configuration and settings for the air-cooling channels.

## 2 Targeted Air-Cooled Electric Motors Driving Aviation Fans

This section describes a targeted air-cooled IPM motor driving an aviation fan. Assumptions used to derive mathematical formulations and their solutions are defined.

### 2.1 Interior Permanent Magnet Motors Using Normal Conductors.

The IPM motor is an alternating current synchronous electric motor. Figure 2 shows a schematic of the IPM motor. In general, a rotor is located inside a stator with a cylindrical shape. In the rotor, permanent magnets are placed inside a yoke made of a pile of magnetic steel sheets, from which the name “IPM” is derived. Permanent magnets do not require electricity. This method has many merits: there is no Joule heating loss, no electrical wiring, and no brush, and the system features a simple structure. The yoke concentrates its magnetic force on the surface of the rotor. Therefore, the rotor efficiency of the IPM motor is increased, and the rotor can be very compact. In addition, compared with surface permanent magnetic motors, the IPM motor can operate at a higher rotational rate because the permanent magnets are held tight by the yoke against a centrifugal force. Furthermore, the IPM motor rotor rotates synchronously with the frequency of a rotating magnetic field composed of electric coils under a wide range of loads. Thus, the rotation rate is highly controllable.

As a permanent magnet, a “neodymium magnet” made of Nd_{2}Fe_{14}B is the best choice because it has the strongest magnetic flux of approximately 1.3 T, which is the value for a mass-produced neodymium magnet. The most advanced neodymium magnet with an EH or VH suffix can be used at a temperature of approximately 200 °C (392 °F), which is given as the catalog value [23]. However, this catalog value corresponds to the maximum allowable temperature when a neodymium magnet is positioned without a magnetic field, but the maximum allowable temperature must be lower to prevent demagnetization when a neodymium magnet is placed in a strong magnetic field such as a rotor in an IPM motor. Therefore, we set the maximum operating temperature of the neodymium magnet to 100 °C based on the experiences of our colleagues. In the rotor, heat is generated by hysteresis loss and eddy-current loss. The generated heat from the hysteresis loss and the eddy-current loss is relatively small compared with the Joule heating loss in the target condition, although these losses depend on operating conditions such as the rotational speed, power, voltage, current, and shape of each component.

In the stator, there are electromagnetic coils and a yoke, and heat is generated by Joule heating loss, hysteresis loss, and eddy-current loss. The Joule heating loss is dominant. To maintain the material health of the thermal insulators, we set the maximum operating temperature of the coils below a temperature of 250 °C (482 °F).

Consequently, more heat is generated in the stator than in the rotor of an IPM motor. Thus, the stator temperature is higher than the rotor temperature in general. In addition, the annular air layer between the rotor and stator is narrow due to electromagnetic performance, and direct heat transfer (DHT) via the annular air layer should count. Thermal barrier coatings on the outer surface of the rotor and the inner surface of the stator are effective in reducing direct heat transfer via the annular air layer.

For a targeted IPM motor, we presume the following factors:

Assumption 1. The output of the IPM motor is 2 MW. In a practical flight, the required output of the motor at the top-of-climb could be less than that at the takeoff. However, this design allows pilots to access the maximum output power at the top-of-climb; this constitutes safe redundancy.

Assumption 2. The rotational rate is 8000 rpm (i.e., 838 rad/s). In motor-driven fans, the rotation rate could be determined to suit the target thrust and flying velocity different from gas turbine-driven fans; thus, the rotation rate at the maximum output power was selected.

Assumption 3. The outer diameter of the rotor,

*d*_{RL}, is 0.25 m, the outer diameter of the permanent magnets,*d*_{PML}, is 0.245 m, the inner diameter of the permanent magnets,*d*_{PMM}, is 0.20 m, the inner diameter,*d*_{RM}, of the rotor's yoke is 0.15 m, and the shaft diameter,*d*_{A}, is 0.10 m. The yoke of the rotor is made of magnetic steel, and its thermal conductivity,*λ*_{steel}, is 50 W/(m·K). A thermal barrier coating (TBC) with a thickness,*δ*_{TBC,R}, of 0.0005 m is placed on the outer surface of the rotor, a TBC with a thickness,*δ*_{TBC,S}, of 0.0005 m is placed on the inner surfaces of the stator, and their thermal conductivities,*λ*_{TBC,S}and*λ*_{TBC,R}, are 0.025 W/(m·K).Assumption 4. An air layer with a thickness,

*δ*_{Air}, of 0.001 m is between the TBCs. The inner diameter of the stator,*d*_{SM}, is*d*_{RL}+ 2(*δ*_{TBC,R}+*δ*_{Air}+*δ*_{TBC,S}), the outer diameter of the coils' section,*d*_{SCL}, is 0.30 m, and the outer diameter of the stator's yoke,*d*_{SL}, is 0.40 m. The length of both the stator and rotor,*ℓ*, is 0.50 m. The yoke of the stator is made of magnetic steel with*λ*_{steel}of 50 W/(m·K).Assumption 5. A heat exchanger is placed outside of the stator, which has a two-story structure. External air flows through the channels on the floor closer to the stator, and internal air flows through the channels on the floor farther from the stator. Both channels of the heat exchanger, for internal and external air, have square shapes of

*a*×*b*mm^{2}, and the number of channels is*n*. The channels are made of aluminum alloy, and the thermal conductivity of the heat exchanger,*λ*_{Al}, is 200 W/(m·K). The thickness,*δ*_{f}, of the floor of the heat exchanger attached on the outer surface of the stator is 0.005 m, and the thickness,*δ*_{b}, of the bulkhead between the external and internal air of the heat exchanger is 0.0005 m.Assumption 6. The efficiency is 99% or higher, so the rotor produces a heat rate,

*Q*_{R}, of 1 kW and the stator produces a heat rate,*Q*_{S}, of 15 kW. These values were derived by electromagnetic analyses corresponding to the maximum output and the worst case scenario. Other cases indicated lesser values of*Q*_{R}and*Q*_{S}even when the efficiency decreased because the output decreased more. With direct heat transfer,*Q*_{S→R}, from the stator via the annular air layer to the rotor, cooling air flows should remove*Q*_{R}+*Q*_{S→R}from the rotor, and cooling air flows should remove*Q*_{S}−*Q*_{S→R}from the stator.Assumption 7. The temperature of the rotor embedding permanent magnets never exceeds 100 °C, and the temperature of the stator does not exceed 250 °C.

### 2.2 Air-Cooling Configuration Schemes for Interior Permanent Magnet Motors.

The rotor is isolated from the motor casing, so internal air should cool the rotor. The outer shape of the rotor is fully cylindrical in a fixed stator. Many conventional motors undergo air cooling, in which air flow between the rotor and the stator cools the rotor. Such configurations that flow onto the outer surface of the rotor are Taylor–Couette or Taylor–Couette–Poiseuille flows [24]. However, there are several problems, which are described below.

Problem 1: The cylindrical channel between the rotor and stator is generally narrow to ensure electromagnetic performance; thus, the pressure loss is relatively large, and it is difficult to achieve a sufficiently large mass flowrate for the internal air flow through the cylindrical channel. A motor driving an aviation fan should have a small aspect ratio of its diameter/length to keep an airflow channel downstream of the fan; thus, the temperature of the internal air gradually increases in the axial direction until it reaches the rotor temperature, which occurs at the halfway point along the axial direction. At this point and beyond, the internal air does not cool the rotor because the air temperature is the same or higher than the rotor temperature.

Problem 2: In the targeted IPM motors, the stator generates heat, *Q*_{S}, larger than the heat, *Q*_{R}, generated by the rotor; thus, the temperature on the inner surface of the stator may be larger than the temperature on the outer surface of the rotor. In this case, the internal air no longer cools the rotor.

Instead, the inner surface of the rotor is useful in cooling, and the authors presume this air-cooling method for the rotor. The authors determined which air-cooling configuration is best for an IPM motor in a previous paper [22]. Based on these results, the air-cooling configuration shown in Fig. 3 is best; first, internal recirculating air moves heat from the rotor to the heat exchanger attached outside of the stator, then the heat exchanger collects heat from the stator; next, the heat exchanger moves all heat to the external air, and finally, the external air moves all heat to the surroundings.

### 2.3 Air Physical Properties at an Arbitrary Temperature and Pressure.

Aircraft operate in all places across the Earth and must continue to operate even in severe conditions, such as in taking off from a desert in summer, which corresponds to an inlet temperature of 50 °C, or cruising at a high altitude of 13,000 m, where the air pressure is 16 kPa.

In this article, we presume that air is a perfect gas with a molar mass *m*_{mol} of 28.8 × 10^{−3 }kg/mol and a specific heat ratio of 1.4. In the pressure range 16–110 kPa (corresponding to an altitude range of 0–13,000 m) at which almost all civil aircraft operate, the effects of pressure on viscosity and thermal conductivity are negligible, so viscosity *μ*_{air} (Pa·s) and thermal conductivity *λ*_{air} (W/(m K)) are estimated with Sutherland's equations [25]. Additionally, we presume that air is incompressible over each air channel in the rotor or heat exchanger; namely, the physical properties of air are the same as values at the inlet of each air channel, and thus we presume:

Assumption 8. For internal air, the inlet pressure to the rotor is kept at ambient pressure. The outlet pressure from the rotor (corresponding to the inlet pressure to the internal compressor) is less than the ambient pressure due to pressure loss. The internal compressor compresses the internal air, and the outlet pressure from the internal compressor (corresponding to the inlet pressure to the heat exchanger) is the highest pressure in the system. The outlet pressure from the heat exchanger (corresponding to the inlet pressure to the rotor) returns to the ambient pressure.

Assumption 9. The fan and external compressor compress the external air, and the outlet pressure from the external compressor (corresponding to the inlet pressure to the heat exchanger) reaches the highest pressure. The outlet pressure from the heat exchanger decreases to ambient pressure due to pressure loss by the heat exchanger.

## 3 Mathematical Analytical Solutions for Temperatures of the Air-Cooled Interior Permanent Magnet Motor

This section derives the mathematical formulations for each heat-transfer process in the air-cooled IPM motor.

### 3.1 Direct Heat Transfer From the Stator Coils to the Rotor Permanent Magnets Via the Annular Air Layer.

Direct heat transfer, from the hot coils in the stator, via the inner surface of the stator, the thermal barrier coatings on the inner surface of the stator, the annular air layer, the thermal barrier coatings on the outer surface of the rotor, the outer surface of the rotor, and the electromagnetic steel to the cold permanent magnets in the rotor cannot be ignored. The heat transfer consists of conduction, convection, and radiation. In cases of the target temperature range, the effect of the radiation is less than the conduction and convection (see Appendix A1); here, the direct heat transfer by the conduction and convection is estimated.

_{CL}

*ρ*

_{air}and

*μ*

_{air}are the density and viscosity of air, respectively, which are the average values of

*T*

_{S}and

*T*

_{R}

*. ω*

_{R}is the rotational rate of the rotor, and

*δ*

_{C}is the difference in diameter between coaxial cylinders of

*d*

_{SM}− 2

*δ*

_{TBC,}

_{S}− 2

*δ*

_{TBC,}

_{R}−

*d*

_{RL}. Nu

_{CL}is defined as

where *λ*_{air} is the thermal conductivity at the average temperature of *T*_{S} and *T*_{R}*. K*_{air} is the coefficient for overall heat transfer via the annular air layer between two thermal barrier coatings. *K*_{air} is defined at *d*_{RL} + 2*δ*_{TBC,}_{R} of the outer surface of the rotor coating.

*K*

_{S}

_{→}

_{R}from the coils, via the thermal barrier coating, the annular air layer, the thermal barrier coating, and the yoke, to the permanent magnet is estimated.

*K*

_{S}

_{→}

_{R}is defined at

*d*

_{RL}of the outer surface of the rotor

*Q*

_{S}

_{→}

_{R}between the coils and the permanent magnets is

### 3.2 Heat Transfer Between the Internal Air and the Rotor.

The rotor and shaft are mechanically connected, so they both rotate at the same rate *ω* (rad/s), as shown in Fig. 4. The rotor generates heat of *Q*_{R} (W), and the direct heat *Q*_{S→R} is added; the internal air recirculating in a motor case removes all heat via the inner surface of the rotor. This type of heat transfer was experimentally investigated by Seghir-Ouali et al. [27].

*P*

_{IR}(m) normal to the axis is

where *d*_{RM} (m) is the diameter of the inner surface of the rotor.

*U*

_{A}through the rotor and tangential velocity

*U*

_{R}of the inner surface of the rotor are

where *ω* (rad/s) and *X* (rpm) are rotational rates of the rotor.

_{A}and the rotational Reynolds number Re

_{R}of the internal air through the rotor are

^{5}≤ Re

_{R}), the experimental correlation between Nusselt and rotational Reynolds numbers proposed in Ref. [27] is

^{3}< Re

_{R}< 2.77 × 10

^{5}and 0 < Re

_{R}< 3 × 10

^{4}), another experimental correlation was proposed in Ref. [27]

*h*

_{R}is

*λ*

_{IR}(W/(m K)) is the thermal conductivity of the internal air.

*K*

_{R}is the coefficient for overall heat transfer via the yoke and heat resistance on the inner surface of the rotor because the permanent magnets are placed inside the yoke.

*K*

_{R}is defined at

*d*

_{RM}of the inner surface of the rotor

*T*

_{R,t}(K). On the other hand, the temperature

*T*

_{IR}

_{x}(K) of the internal air is increased by heat received from the inner surface of the rotor. At position

*x*(m) in the axial direction from the inlet of the rotor, the local heat flux

*q*

_{IR}

_{x}(W/m

^{2}) for an infinitesimal region in the axial direction with a length of

*dx*(m) is

*C*

_{p}(J/(kg K)) is the specific heat of the internal air. Then

*c*

_{0}is an integration constant. By applying the boundary condition

*T*

_{IR}

_{x}=

*T*

_{IRin}at

*x*=

*0*

*T*

_{IR}

_{x}=

*T*

_{IRout}is applied at

*x*=

*ℓ*

_{R}; finally, the outlet temperature

*T*

_{IRout}of the internal air is solved analytically

where *m*_{R} (kg) is the mass of the rotor, *C*_{P,R} (J/(kg·K)) is the specific heat of the rotor, and *T*_{R,t} (K) is the temperature of the permanent magnets at time *t* (s). *T*_{R,t} can be calculated over time; thus

Case 1: In cases of

*Q*_{R}+*Q*_{S→R}>*Q*_{I←R},*T*_{R,t}increases with time.Case 2: In cases of

*Q*_{R}+*Q*_{S→R}<*Q*_{I←R},*T*_{R,t}decreases with time.Case 3: In cases of

*Q*_{R}+*Q*_{S→R}=*Q*_{I←R},*T*_{R,t}is steady.

### 3.3 Heat Transfer Among the Stator, External Air, and Internal Air in the Heat Exchanger Attached Outside of the Stator.

Figure 6 shows a schematic for the configuration with a stator-external air-internal air heat exchanger. This heat exchanger is directly attached on the outer surface of the stator yoke. This heat exchanger consists of two-story structures: channels for the external air are located on the floor closer to the stator, and channels for the internal air are located further away on the floor. The stator generates heat of *Q*_{S} (W), and the direct heat of *Q*_{S→R} is subtracted. Solid structures of the stator yoke and the heat exchanger directly conduct heat generated by the stator to the external air from the bottom side. At the same time, the bulkhead of the heat exchanger passes heat from the internal air to the external air from the top side.

*V*

_{I}and

*V*

_{E}of internal and external air passing through the heat exchanger are

where $m\u02d9I$ and $m\u02d9E$ (kg/s) are mass flow rates, *ρ*_{I} and *ρ*_{E} (m^{3}/kg) are densities, *n*_{I} and *n*_{E} (m) are the numbers of channels, *a*_{I} and *a*_{E} are the heights, and *b*_{I} and *b*_{E} (m) are widths of the internal and external air channel, respectively.

*d*

_{I}and

*d*

_{E}(m) are

*S*

_{I}and

*S*

_{E}(m

^{2}) are the cross section and

*P*

_{I}and

*P*

_{E}(m) are the peripheral lengths of the internal and external air channel, respectively. Pipe Reynolds numbers Re

_{P,I}and Re

_{P,E}of the internal and external air through each channel in the heat exchanger are

_{I}and Pr

_{E}of the internal and external air passing through each air channel in the heat exchanger are

_{P,I}of the internal air with fully developed temperature fields in laminar flow (Re

_{P,I}< 3000) or turbulent flow (3000 ≤ Re

_{P,I}< 10

^{7}) in a pipe with a constant wall temperature is

_{P,E}of the external air is

*h*

_{I}and

*h*

_{E}are

*λ*

_{I}and

*λ*

_{E}(W/(m K)) are the thermal conductivities of the internal and external air, respectively. The heat transfer coefficient between the stator and the external air is

*h*

_{E}, and the overall heat transfer coefficient between the external air and the coil via the yoke is

*K*

_{S}(W/(m

^{2}K)). Here,

*K*

_{S}is defined at

*d*

_{SL}+ 2

*δ*

_{F}for the floor closer to the stator of the heat exchanger

*K*

_{B}(W/(m

^{2}K)), is

where *δ*_{B} (m) is the thickness and *λ*_{Al} (W/(m K)) is the thermal conductivity of the bulkhead between the internal and external air.

*x*(m) in the axial direction from the left end of the heat exchanger, for an infinitesimal region in the axial direction with a length of

*dx*(m), local heat flux

*dq*

_{S}

_{→}

_{Ex}(W/m) is exchanged between the stator and the external air, and local heat flux

*dq*

_{I}

_{→}

_{Ex}(W/m) is exchanged between the internal air and the external air. The coils are made of copper; thus, the thermal conductivity in the axial direction is much greater than that in air, and the temperature of the heat exchanger can be presumed to be constant in the axial direction at

*T*

_{S,t}(K). The internal air flows in the opposite direction against the external air flow because this configuration increases the heat exchange efficiency. Namely, the internal air flows from the right end to the left end in Fig. 7. The temperature

*T*

_{Ix}(K) of the internal air decreases as heat

*dq*

_{I}

_{→}

_{Ex}is released to the external air via the bulkhead. On the other hand, the external air receives

*dq*

_{S}

_{→}

_{Ex}from the stator and

*dq*

_{I}

_{→}

_{Ex}from the internal air via the bulkhead at the same time. Thus, the temperature

*T*

_{Ex}(K) of the external air increases. Based on these characteristics, the energy balances of the internal and external air are

*T*

_{Ex}and

*T*

_{Ix}; thus, we apply boundary conditions of

*T*

_{I}

_{0}=

*T*

_{I}

_{out}and

*T*

_{E}

_{0}=

*T*

_{E}

_{in}at

*x*=

*0 and*

*T*

_{IℓR}=

*T*

_{I}

_{in}and

*T*

_{EℓR}=

*T*

_{E}

_{out}at

*x*=

*ℓ*

_{R}(see Ref. [22] for details). Then, the analytical solutions are obtained from [22]

*Q*

_{I→E}transferred from internal to external air is equal and opposite to the heat

*Q*

_{I}received by the internal air; thus

*Q*

_{I→E}and exchange heat

*Q*

_{S→E}from the stator to external air equals the heat

*Q*

_{E}received by the external air; thus

*Q*

_{S→E}and heat generation rate

*Q*

_{S}in the stator changes the stator temperature

*T*

_{S,t}; thus

where *m*_{S} (kg) is the mass of the stator, *C*_{P,S} (J/(kg K)) is the specific heat of the stator, and *T*_{S,t} (K) is the temperature of the coils at time *t* (s). Therefore, *T*_{S,t} can be calculated over time, so

Case 1: In cases of

*Q*_{S}−*Q*_{S→R}>*Q*_{I}+*Q*_{E},*T*_{S,t}increases with time.Case 2: In cases of

*Q*_{S}−*Q*_{S→R}<*Q*_{I}+*Q*_{E},*T*_{S,t}decreases with time.Case 3: In cases of

*Q*_{S}−*Q*_{S→R}=*Q*_{I}+*Q*_{E},*T*_{S,t}is steady.

## 4 Pressure Loss and Power Consumption of Cooling Air in the Air-Cooled Interior Permanent Magnet Motor

This section estimates the pressure difference for each air-cooling channel contained in the targeted air-cooled IPM motor.

### 4.1 Pressure Loss and Power Consumption of Internal Air Passing Through the Rotor.

*f*

_{IR}with coaxial Reynolds number Re

_{C}and rotational Reynolds number Re

_{R}[28]. The definition of Re

_{R}is the same as in Eq. (8), but coaxial Reynolds number Re

_{d,C}is newly defined. Here, the inner diameter of the rotor is

*d*

_{RM}, the outer diameter of the shaft is

*d*

_{A}, and the difference between

*d*

_{RM}and

*d*

_{A}is

*δ*

_{C}. By considering a hydraulic diameter

*d*

_{C}(m) of the channel between the coaxial cylinders

*S*(m

^{2}) is the cross section of the channel and

*P*(m) is the peripheral length of the channel. Namely, the hydraulic diameter

*d*

_{C}of the channel between the coaxial cylinders equals the difference

*δ*

_{C}between

*d*

_{RM}and

*d*

_{A}. Thus, the coaxial Reynolds number Re

_{C}for the internal air passing through the rotor is

*f*

_{IR}is almost the same as

*f*

_{IR}for no rotation even in laminar and turbulent flows in cases where Re

_{C}is not very small and both cylinders rotate at the same rate [28]. Therefore, the pressure loss coefficient

*f*

_{IR}is

for turbulent flows.

_{C}changes with Re

_{R}as laminar flow changes to turbulent flow [28]; however, we presume laminar flows in cases of 0 < Re

_{C}< 2000 and turbulent flows in cases of 2000 ≤ Re

_{C}based on experimental results reported by Yamada and Watanabe [28]. Here, the total pressure loss Δ

*p*

_{lossIR}(Pa) is

### 4.2 Pressure Loss and Power Consumption of the Internal and External Air Passing Through the Heat Exchanger Attached Outside of the Stator.

As shown in Fig. 6, the channels for internal air and external air are rectangular with a height of *a* and width of *b*; thus, hydraulic diameters *d*_{I} and *d*_{E} and pipe Reynolds numbers Re_{P,I} and Re_{P,E} were defined in Eq. (22). The total pressure loss is estimated with the Hargen–Poiseuille and Blasius equations in association with the Darcy–Weisbach equation.

*f*

_{I}and

*f*

_{E}are

_{P,I}< 3000 or Re

_{P,E}< 3000), and

for turbulent flows (3000 ≤ Re_{P,I} or 3000 ≤ Re_{P,E}).

*p*

_{lossIH}and Δ

*p*

_{lossEH}(Pa) are

### 4.3 Effect on the Energy Balance of Power Consumption by Compressors for Internal and External Air.

*W*

_{E}

*ρ*

_{E}

_{in}is the density at the fan inlet and

*η*

_{comp}

_{E}is the adiabatic efficiency of the fan and external compressor. Power

*W*

_{E}increases the temperature

*T*

_{E}

_{in}at the inlet of the heat exchanger

*W*

_{I}consumed to force the internal air through the rotor and the heat exchanger is

*ρ*

_{IRout}is the density at the inlet of the internal compressor and

*η*

_{comp}

_{I}is the adiabatic efficiency of the internal compressor. Power

*W*

_{I}increases the temperature

*T*

_{I}

_{in}at the inlet of the heat exchanger and also increases the heat

*Q*

_{I→E}exchanged from the internal to external air because this power

*W*

_{I}is ultimately converted into heat; thus

Equation (51) indicates that the power consumption *W*_{I} affects the rotor temperature *T*_{R}.

## 5 Effect of Direct Heat Transfer and Heat Resistance on Air-Cooling Performance

### 5.1 Difference Between the Current Model and the Authors' Conventional Model.

To obtain mathematical analytical solutions that maintain the maximum temperature limits for permanent magnets in the rotor and coils in the stator, the authors and colleagues introduced a mathematical model that describes heat transfer inside the targeted air-cooled IPM motor [22]. The solutions indicated many characteristics and provided considerable knowledge regarding the targeted IPM motor under severe operational conditions, such as top-of-climb.

However, as shown in Table 1, the temperature in the entire rotor at the permanent magnets and the inner heat-transfer surface were assumed to be constant in the conventional model. The temperature over the entire stator at the coils and the outer heat-transfer surface were also assumed to be constant. In practical situations, temperatures at the permanent magnets and coils must be higher than that in the isothermal stator and rotor because both the rotor and stator have internal heat resistance. In addition, in the model in Ref. [22], the rate of direct heat transfer from the stator coils via the annular air layer to the rotor permanent magnets was implied for the constant heat generated from the stator and rotor. In practical situations, the rate of direct heat transfer must be altered based on the temperature difference between the coils and permanent magnets.

Physics | Current model | Conventional model in Ref. [22] |
---|---|---|

Heat resistance of electromagnetic steel between the rotor inner surface and the permanent magnets when calculating the heat exchange rate between the permanent magnets and the internal air | Considered | Not considered |

Heat resistance of electromagnetic steel and aluminum between the coils and the air channel surfaces of the heat exchanger when calculating the heat exchange rate between the coils and the external air | Considered | Not considered |

DHT between the coils and the permanent magnets in the rotor | Considered | Not considered |

Heat resistance of the thermal barrier coating on the stator inner surface when calculating DHT | Considered | Not considered |

Heat resistance of the thermal barrier coating on the rotor outer surface when calculating DHT | Considered | Not considered |

Physics | Current model | Conventional model in Ref. [22] |
---|---|---|

Heat resistance of electromagnetic steel between the rotor inner surface and the permanent magnets when calculating the heat exchange rate between the permanent magnets and the internal air | Considered | Not considered |

Heat resistance of electromagnetic steel and aluminum between the coils and the air channel surfaces of the heat exchanger when calculating the heat exchange rate between the coils and the external air | Considered | Not considered |

DHT between the coils and the permanent magnets in the rotor | Considered | Not considered |

Heat resistance of the thermal barrier coating on the stator inner surface when calculating DHT | Considered | Not considered |

Heat resistance of the thermal barrier coating on the rotor outer surface when calculating DHT | Considered | Not considered |

### 5.2 Effect of the Air Mass Flow Rate on Temperatures of the Permanent Magnet and Coil.

Figure 8 compares the results estimated with the current model (solid curves) and with the conventional model (dashed curves). The left figures show the temperature changes in permanent magnets *T*_{R} (in blue) and coils *T*_{S} (in red) with changes in the mass flowrate of a focused air flow in the horizontal axis, and the right figure shows the power consumption needed to drive the external air *W*_{E} (in green) and the internal air *W*_{I} (in orange) corresponding to *T*_{R} in the vertical axis. Refer to Ref. [22], where the detailed tendencies were explained. Differences between the results of the current model and the conventional model are considered in Sec. 5.

Figure 8(a) shows the results for cases with a constant external air mass flowrate of 1.4 kg/s. The coil temperature *T*_{S} (red solid curve) estimated by the current model is higher than the *T*_{S} (red dashed curve) estimated by the conventional model because of heat resistance from the coils to the external air via the magnetic steel and the heat exchanger. In the region where $m\u02d9I\u22650.07$ kg/s (solid curve) or $m\u02d9I\u22650.045$ kg/s (dashed curve), the internal air flow is turbulent in the heat exchanger, and heat is exchanged well. Thus, as shown in Fig. 3, the inlet temperature *T*_{IRin} of the rotor sufficiently decreased, and the permanent magnets were well cooled, so *T*_{R} was always lower than *T*_{S}. Instead, in the region where $m\u02d9I<0.07$ kg/s (solid curve) or $m\u02d9I<0.045$ kg/s (dashed curve), the internal air flow is laminar in the heat exchanger, and heat is not exchanged sufficiently. Thus, as shown in Fig. 3, the inlet temperature *T*_{IRin} of the rotor is relatively high, and the permanent magnets are not sufficiently cooled, so *T*_{R} has a gap within this region of $m\u02d9I$.

Direct heat transfer *Q*_{S→R} flows from the hotter part to the colder part. In the turbulent region of $m\u02d9I$, the coil temperature *T*_{S} is higher than the permanent magnet temperature *T*_{R}; thus, *Q*_{S→R} is positive, and the temperature difference of *T*_{S} (red solid curve) − *T*_{R} (blue solid curve) estimated by the current model is greater than that of *T*_{S} (red dashed curve) − *T*_{R} (blue dashed curve) estimated by the conventional model. On the other hand, in the laminar region of $m\u02d9I$, the coil temperature *T*_{S} is lower than the permanent magnet temperature *T*_{R}; thus, *Q*_{S→R} is negative, and the temperature difference of *T*_{R} (blue solid curve) − *T*_{S} (red solid curve) estimated by the current model is smaller than that of *T*_{R} (blue dashed curve) − *T*_{S} (red dashed curve) estimated by the conventional model.

An increase in $m\u02d9I$ causes a decrease in the permanent magnet temperature *T*_{R} in a smaller $m\u02d9I$ region because the pipe Reynolds number Re_{P,I} and Nusselt number Nu_{P,I} increase and the internal air is cooled better through the heat exchanger. Then, in a larger $m\u02d9I$ region, an increase in $m\u02d9I$ causes a rebound in the permanent magnet temperature *T*_{R} because the pressure loss used to recirculate the internal flow increases and the power consumption *W*_{I} increases. *W*_{I} is finally converted into heat in the internal air, and this heat should be transferred from the internal air to the external air. Therefore, the permanent magnet temperature *T*_{R} gradually increases.

Figure 8(b) shows the results for cases with a constant internal air mass flowrate of 0.28 kg/s. Under this condition, the flow of internal air through the heat exchanger is always turbulent so the coil temperature *T*_{S} is higher than the permanent magnet temperature *T*_{R}; thus, *Q*_{S→R} is positive, and *T*_{S} (red solid curve) estimated by the current model is higher than *T*_{S} (red dashed curve) estimated by the conventional model.

Although the internal air mass flowrate $m\u02d9I$ is constant, the power *W*_{I} (orange curves) consumed by the internal air decreases with increasing $m\u02d9E$ (corresponding to decreasing *T*_{R}). This is because an increase in $m\u02d9E$ causes a decrease in the temperature of the entire system; thus, the internal air is cooled better, its density increases, its velocity decreases and *W*_{I} decreases.

### 5.3 Effect of Ambient Temperature on Air-Cooling Performance.

Figure 9 shows the effect of the ambient temperature *T*_{amb} on air-cooling performance, i.e., the results for *T*_{R} (light gray curves) and *T*_{S} (dark gray curves) and the power consumed by internal air flow *W*_{I} (light gray curves) and external air flow *W*_{E} (dark gray curves) with *T*_{amb} in the case of constant volume flow rates $V\u02d9E$ of 1.4 m^{3}/s and $V\u02d9I$ of 0.32 m^{3}/s, which correspond to $m\u02d9E$ = 1.7 kg/s and $m\u02d9I$ = 0.34 kg/s at the ambient temperature *T*_{amb} of 15 °C and the ambient pressure *p*_{amb} of 101.3 kPa. The reason why volume flow rates $V\u02d9E$ and $V\u02d9I$ are fixed is that a compressor operating at a constant rotational rate generally discharges a constant volume flowrate when the inlet temperature and pressure change. For practical flight, the temperature at takeoff has a wide range, from −40 °C (at an airport in a polar region in winter) to 50 °C (at an airport in a desert region in summer), so we estimate *T*_{R} and *T*_{S} in this temperature range. As expected, a higher *T*_{amb} results in higher *T*_{R} and *T*_{S}. Both *T*_{R} and *T*_{S} are linearly proportional to *T*_{amb}.

### 5.4 Effect of Ambient Pressure on Air-Cooling Performance.

Figure 10 shows the effects of ambient pressure *p*_{amb} on *T*_{R}, *T*_{S}, *W*_{E}, and *W*_{I} in the case of constant $V\u02d9E$ and $V\u02d9I$. Civilian aircraft usually fly at altitudes ranging within 0 m and 13,000 m, so ambient pressures *p*_{amb} range between 101.3 kPa and 16 kPa based on the ISA. We estimate *T*_{R} and *T*_{S} in this pressure range at a constant temperature of 15 °C because we first want to know the effect of pressure alone on air-cooling performance.

*T*_{S} increases with decreasing *p*_{amb} because *ρ*_{E} decreases and $m\u02d9E$ decreases at a constant volume flowrate $V\u02d9E$ of 1.4 m^{3}/s.

*T*_{R} increases with decreasing *p*_{amb}. This is because the density *ρ*_{IR}, rotating Reynolds number Re_{R} in Eq. (8), and pipe Nusselt number Nu_{P,R} in Eq. (9) between the rotor and the internal air decrease. This is also because the density *ρ*_{I}, pipe Reynolds number Re_{P,I} in Eq. (22), and pipe Nusselt number Nu_{P,I} in Eq. (25) of the internal air passing through the heat exchanger decrease. When *p*_{amb} values are between 56 kPa and 54 kPa, the flow of the internal air through the rotor is switched from the rotation-dominant flow to the rotation-inertia-affected flow. Thus, Nu_{P,R} is switched from in Eq. (9) to in Eq. (10), and there is a gap in *T*_{R} when *p*_{amb} values are between 56 kPa and 54 kPa. In addition, when *p*_{amb} values are between 28 kPa and 26 kPa, the flow of the internal air through the heat exchanger is switched from turbulent flow to laminar flow. Thus, Nu_{P,I} is switched from in Eq. (25) to in Eq. (24), and there is a gap in *T*_{R} when *p*_{amb} values are between 28 kPa and 26 kPa.

Here, we consider the effect of pressure *p*_{amb} on power consumption *W*_{E}. Power consumption *W*_{E} is proportional to pressure loss Δ*p*_{lossEH}. From Eq. (46), Δ*p*_{lossEH} is proportional to *ρ*_{E}*f*_{E}. From Eq. (44), *f*_{E} is inversely proportional to Re_{P,E} in a laminar flow, and from Eq. (45), *f*_{E} is inversely proportional to Re_{P,E} to the power of 1/4 in a turbulent flow. From Eq. (22), Re_{P,E} is proportional to $m\u02d9E$, and $m\u02d9E$ is proportional to *ρ*_{E} because of the constant volume flowrate $V\u02d9E$. From air characteristics, the density *ρ*_{E} is proportional to *p*_{amb} under a constant temperature. Consequently, *W*_{E} is almost independent of *p*_{amb} in a laminar flow region, and *W*_{E} is almost proportional to *p*_{amb} to the power of 3/4 in a turbulent flow under a constant volume flowrate $V\u02d9E$. *W*_{I} also has the same tendency.

### 5.5 Effect of Altitude on Air-Cooling Performance.

Figure 11 shows results for *T*_{E}_{in} (gray dashed dotted curve), *p*_{E}_{in} (dashed double-dotted curve), *T*_{R} (light gray curves), *T*_{S} (dark gray curves), the power consumption of the internal air flow *W*_{I} (light gray curves) and the external air flow *W*_{E} (dark gray curves) with changing altitude under conditions based on ISA standard air in the case with a $V\u02d9E$ of 1.4 m^{3}/s and a $V\u02d9I$ of 0.32 m^{3}/s. *T*_{amb} decreases by 6.5 °C/1,000 m below *A *=* *11,000 m, but *T*_{amb} becomes constant at −56.5 °C above *A *=* *11,000 m. *p*_{amb} gradually decreases with increasing *A*.

*T*_{S} is affected by both *T*_{amb} and *p*_{amb}. The effect of *T*_{amb} prevails over the effect of *p*_{amb} until approximately *A *=* *2,000 m; thus, *T*_{S} gradually decreases. *T*_{S} rebounds beyond that point because the effect of *p*_{amb} gradually prevails over the effect of *T*_{amb} with increasing *A*. Above *A *=* *11,000 m, only the effect of *p*_{E}_{out} operates, and *T*_{S} rapidly increases.

*T*_{R} also has the same tendency as *T*_{S}, although *T*_{R} has the minimum value at approximately *A *=* *6,500 m. There is a gap in *T*_{R} when *A* values are between 6,250 m and 6,500 m because the flow of the internal air through the rotor is the rotation-dominant flow at *A* less than 6,250 m, but it is the rotation-inertia-affected flow at *A* greater than 6,500 m. There is another gap in *T*_{R} when *A* values are between 11,250 m and 11,500 m because the flow of the internal air through the heat exchanger is turbulent at *A* less than 11,250 m, but it is laminar at *A* greater than 11,500 m.

Decreasing *T*_{E}_{in} increases the external air density *ρ*_{E}, and decreasing *p*_{E}_{in} decreases *ρ*_{E}. With increasing altitude, the pressure effect prevails over the temperature effect; thus, *ρ*_{E}, $m\u02d9E$, and *W*_{E} decrease with increasing altitude *A* at a constant $V\u02d9E$. As with *ρ*_{E}, the internal air density *ρ*_{I}, $m\u02d9I$, and *W*_{I} decrease with increasing altitude *A* for a constant $V\u02d9I$.

Civilian aircraft usually fly within an altitude range of 0 m to 13,000 m, so the most severe condition is at *A *=* *13,000 m in Fig. 11. Although this result does not satisfy the maximum temperature limit of 100 °C for the permanent magnets, it was found that the top-of-climb condition was the most severe, which derives the highest temperatures of permanent magnets and coils.

## 6 Effect of Air-Cooling Channel Sizes on Air-Cooling Performance for the Top-of-Climb Condition

This section optimizes the sizes of air-cooling channels to reach the minimum value of the sum of power consumption of the external and internal air flows with a setting that keeps permanent magnet temperatures below the maximum temperature limit of 100 °C and coil temperatures below the maximum temperature limit of 250 °C.

### 6.1 Effect of Clearance Between Rotor and Shaft.

As shown in Fig. 4, there is an annular air layer between the inner surface of the rotor and the outer surface of the shaft. Figure 12 shows the coil temperature *T*_{S} and permanent magnet temperature *T*_{R} as a function of changing rotor inner diameter *d*_{R}. There is a gap in *T*_{R} when the *d*_{R} values are between 0.11 m and 0.1125 m. Decreasing *d*_{R} increases the pressure loss through the rotor, increases *W*_{I}, increases the inlet temperature of the heat exchanger *T*_{I}_{in}, decreases *ρ*_{I}_{in}, decreases Re_{P,I}, and switches the internal air flow through the heat exchanger from turbulent to laminar at a *d*_{R} of 0.11 m. A much larger *d*_{R} (beyond the range of Fig. 12) would switch the internal air flow through the rotor from the rotation-inertia-affected flow to the rotation-dominant flow, and another gap in *T*_{R} would appear.

### 6.2 Effect of Clearance Between Rotor and Stator.

There is an annular air layer in the clearance space between the rotor and stator. The thickness of the annular air layer affects the rate of direct heat transfer between the hotter coils and the colder permanent magnets, and thus, a larger thickness of the annular air layer is preferable for decreasing the direct heat transfer rate in general. On the other hand, a smaller thickness of the annular air layer is preferable for strengthening the electromagnetic force. This paradox may be resolved by applying thermal barrier coatings on the inner surface of the stator and the outer surface of the rotor. Figure 13 shows the effects of the thickness of each thermal barrier coating on the direct heat transfer rate *Q*_{S→R}, the permanent magnet temperature *T*_{R}, and the coil temperature *T*_{S} when the total thickness of three layers, the annular air layer and two thermal barrier coating layers, is constant. Increasing the thickness of each thermal barrier coating drastically reduces *Q*_{S→R} and slightly decreases *T*_{R}.

In Fig. 14, the effect of the annular air layer thickness is evaluated for the case with a constant thickness for each thermal barrier coating layer. Changing the thickness of the annular air layer has little effect on *Q*_{S→R}, *T*_{S} and *T*_{R}. This is because the heat resistance of the thermal barrier coating layers is much larger than that of the annular air layer; thus, the overall heat transfer coefficient between the coils and permanent magnets *K*_{S→R} is almost decided by the values of the heat resistances *λ*_{TBC,}_{S}/*δ*_{TBC,}_{S} and *λ*_{TBC,}_{R}/*δ*_{TBC,}_{R} in Eq. (4).

### 6.3 Effect of the External-Air-Channel Height and Width of the Heat Exchanger.

Figure 15 shows the coil temperature *T*_{S}, the permanent magnet temperature *T*_{R}, the power consumption by the compressors to drive internal air flow *W*_{I} and external air flow *W*_{E} with changing external-air-channel sizes *a*_{E} = *b*_{E} and the same volume flow rates $V\u02d9E$ = 1.4 m^{3}/s and $V\u02d9I$ = 0.32 m^{3}/s. In this figure, *a*_{E} and *b*_{E} change together and maintain the same value. In a range of smaller values *a*_{E} = *b*_{E}, decreasing the values of *a*_{E} = *b*_{E} causes a larger pressure loss because of a faster *V*_{E}; thus, *W*_{E} increases. At the same time, to compress the external air to resist the larger pressure loss, the external-air inlet temperature of the heat exchanger *T*_{E}_{in} is larger due to adiabatic compression at the fan and the external compressor, and thus, it raises *T*_{E}_{in} and *T*_{S}, resulting in a higher internal-air temperature and a higher *T*_{R}. When *a*_{E} = *b*_{E} is less than 0.004 m, *T*_{I}_{in} is sufficiently high, and a high *T*_{I}_{in} causes a small density *ρ*_{I}_{in} and a small Re_{P,I}; thus, the internal flow converts to laminar flow, and there is a gap in *T*_{R} when *a*_{E} = *b*_{E} is 0.004 m. Instead, in a range with larger values of *a*_{E} = *b*_{E}, increasing the values of *a*_{E} and *b*_{E} results in a higher *T*_{S} and *T*_{R} because of a smaller Reynolds number Re_{P,E} in Eq. (22); thus, the external air flow removes less heat through the heat exchanger. Therefore, it causes a higher *T*_{S}, a higher internal-air temperature, and a higher *T*_{R}. Here, the pressure losses of the external air decrease because of a slower *V*_{E} in Eq. (46); thus, a smaller *W*_{E} is achieved. When *a*_{E} = *b*_{E} is greater than 0.0185 m, Re_{P,I} is less than 3000; thus, the internal air flow converts to laminar flow. Laminar flow has a much lower heat transfer coefficient with air channels, as shown in Eq. (24); therefore, there is a gap in *T*_{R} when *a*_{E} = *b*_{E} is 0.0185 m. With a smaller *T*_{S}, *a*_{E} = *b*_{E} = 0.005 m is the optimum value. With a smaller *T*_{R}, *a*_{E} = *b*_{E} = 0.0075 m is the optimum value. With a smaller *W*_{I}, *a*_{E} = *b*_{E} = 0.008 m is the optimum value. With a smaller *W*_{E}, the largest *a*_{E} = *b*_{E} is the optimum value. One should choose the more stringent condition among *T*_{S}, *T*_{R}, *W*_{I}, and *W*_{E} in considering the overall balance.

Figure 16 shows *T*_{S}, *T*_{R}, *W*_{I}, and *W*_{E} with changing an external-air-channel width *b*_{E} under various constant values of an external-air-channel height *a*_{E} at the volume flow rates $V\u02d9E$ = 1.4 m^{3}/s and $V\u02d9I$ = 0.32 m^{3}/s. In the left figure of Fig. 16, *T*_{S} (red curves) exhibits the minimum value at an *a*_{E} of 0.005 m as in Fig. 15. In a range of *a*_{E} less than 0.005 m, a smaller *a*_{E} derives a larger *T*_{S} because of a larger pressure loss. In a range of *a*_{E} greater than 0.005 m, a larger *a*_{E} derives a larger *T*_{S} because of a smaller Re_{P,E} in Eq. (22). Also, *T*_{R} (blue curves) exhibits the minimum value at an *a*_{E} of 0.0075 m as in Fig. 15. The tendency of *T*_{R} is the same as that of *T*_{S}. In the right figure of Fig. 16, a smaller *a*_{E} derives a larger *W*_{E} because of a larger pressure loss. Instead, *W*_{I} is almost independent of *a*_{E} because *T*_{R} and *T*_{S} do not affect *W*_{I} in the same manner as in Fig. 9.

### 6.4 Effect of Internal-Air-Channel Height and Width of the Heat Exchanger.

Figure 17 shows *T*_{S}, *T*_{R}, *W*_{I}, and *W*_{E} with changing internal-air-channel sizes *a*_{I} = *b*_{I} at the same volume flow rates $V\u02d9E$ = 1.4 m^{3}/s and $V\u02d9I$ = 0.32 m^{3}/s. In this figure, *a*_{I} and *b*_{I} change together while maintaining the same values. In the left figure of Fig. 17, for a range of smaller values of *a*_{I} = *b*_{I} less than 0.004 m, decreasing *a*_{I} = *b*_{I} causes larger pressure loss resulting in a higher temperature of the internal air flow and a higher *T*_{R}. In a range of *a*_{I} = *b*_{I} greater than 0.004 m, larger *a*_{I} = *b*_{I} derives a larger *T*_{R} because of a smaller Re_{P,I} in Eq. (22); thus, the internal air flow radiates less heat through the heat exchanger, and it causes a higher internal-air temperature as well as a higher *T*_{R}. In the right figure of Fig. 17, smaller *a*_{I} = *b*_{I} derives a larger *W*_{I} because of a larger pressure loss. Instead, *W*_{E} is almost independent of *a*_{I} = *b*_{I} because *T*_{R} and *T*_{S} do not affect *W*_{E} the same as Fig. 9. With a smaller *T*_{S}, *a*_{I} = *b*_{I} = 0.0045 m is the optimum value. With a smaller *T*_{R}, *a*_{I} = *b*_{I} = 0.004 m is the optimum value. With a smaller *W*_{I}, the largest *a*_{I} = *b*_{I} is the optimum value. *W*_{E} is constant at any *a*_{I} = *b*_{I}.

Figure 18 shows *T*_{S}, *T*_{R}, *W*_{I}, and *W*_{E} with changing an internal-air-channel width *b*_{I} under various constant values of an internal-air-channel height *a*_{I} at the volume flow rates $V\u02d9E$ = 1.4 m^{3}/s and $V\u02d9I$ = 0.32 m^{3}/s. In the left figure of Fig. 18, *T*_{R} (blue curves) exhibits the minimum value at *a*_{I} of 0.002 m. In a range of *a*_{I} less than 0.002 m, a smaller *a*_{I} derives a larger *T*_{R} because of a larger pressure loss. In a range of *a*_{I} greater than 0.004 m, a larger *a*_{I} derives a larger *T*_{R} because of a smaller Re_{P,I} in Eq. (22). Additionally, there is a gap in *T*_{R} due to the transition from turbulent to laminar flow of the internal air in the heat exchanger caused by a smaller Re_{P,I}. In the right figure of Fig. 18, a smaller *a*_{I} derives a larger *W*_{I} because of a larger pressure loss. Instead, *W*_{E} is almost independent of *a*_{I} because *T*_{R} and *T*_{S} do not affect *W*_{E} in the same manner as Fig. 9.

### 6.5 Air-Cooling Performance at the Optimum Setting.

Figure 19 shows the results for *T*_{E}_{in}, *p*_{E}_{in}, *T*_{R}, *T*_{S}, *W*_{I}, and *W*_{E} with changing altitude under conditions based on the ISA standard for air at a $V\u02d9E$ of 1.4 m^{3}/s and a $V\u02d9I$ of 0.32 m^{3}/s with optimum settings of the air-cooling channel sizes: *d*_{R} = 15 mm, *d*_{A} = 10 mm, *δ*_{TBC} = 0.5 mm, *δ*_{AIR} = 1 mm, *a*_{E} = 7.5 mm, *b*_{E} = 14 mm, *a*_{I} = 4 mm, and *b*_{I} = 14 mm. The entire tendency seen in Fig. 19 is similar to the results estimated with the current model (by solid curves) shown in Fig. 11, although the air-cooling-channel sizes are optimized in Fig. 19, and the air-cooling-channel sizes are not optimized in Fig. 11. The results shown in Fig. 19 satisfy the maximum temperature limit for the permanent magnets and coils, i.e., *T*_{R} should be less than *T*_{R}_{limit} of 100 °C and *T*_{S} should be less than *T*_{S}_{limit} of 250 °C with the maximum power consumption of 5.35 kW (*W*_{I} = 189 W and *W*_{E} = 5.16 kW). The current model results shown in Fig. 11 do not satisfy the maximum temperature limits, and the maximum power consumption is 3.38 kW (*W*_{I} = 41 W and *W*_{E} = 3.34 kW).

Therefore, optimization of air-cooling channel sizes is very important in keeping the temperature limits and minimizing power consumption to enhance the total efficiency of the electric IPM motor. In a practical flight, the required output of the motor at the top-of-climb could be less than that at the takeoff. However, this optimum setting allows pilots to access the maximum output power at the top-of-climb; this constitutes safe redundancy. Of course, the proposed theoretical analysis method is effective at variable outputs and rotation rates along with a practical flight path.

## 7 Conclusion

In this paper, the effect of altitude on air-cooling performance in an air-cooled 2 MW IPM motor driving an aviation fan was theoretically investigated. Compressor-driven internal air cools permanent magnets (PMs) in a rotor to satisfy the PM maximum temperature limit, and it carries heat from the rotor to the heat exchanger placed outside of a stator. Compressor-driven external air cools coils in the stator via the heat exchanger to satisfy the coil maximum temperature limit, and it carries all heat to the surroundings.

In cases of constant volume flowrates of both airs, with increasing the ambient temperature, the temperatures of the PM and coils linearly increase. The power driving the compressors is almost constant despite the changing ambient temperature. Additionally, with decreasing the ambient pressure, temperatures of the PM and coils rapidly increase because the Reynolds number decreases due to decreasing air density; therefore, the flow mode switches between the rotation-dominant flow and rotation-inertia-affected flow of internal air through the rotor, and the flow mode switches between the laminar flow and turbulent flow of internal air through the heat exchanger. Especially, to satisfy both maximum temperature limits, keeping the turbulent flow mode of internal air through the heat exchanger is necessary. With changing ambient pressure, the power driving the compressor is almost constant in the laminar flow mode of internal air through the heat exchanger, and the power is almost proportional to the ambient pressure to the power of 3/4 in the turbulent flow mode. Furthermore, the effect of altitude on temperatures of the PM and coils was investigated. The most severe condition is the top-of-climb condition. Under this condition, the optimized air-cooling channel sizes satisfy the maximum temperature limits for the PM and coils although the nonoptimized air-cooling channel sizes do not satisfy the maximum temperature limits at the same volume flow rates of internal and external cooling airs.

How and where to switch each flow mode depends on each physical model and correlation used. If each physical model and correlation would be more improved, this approach would derive a more accurate estimation. This approach is very useful for early thermal designs of air-cooled electric motors driving aviation fans.

## Acknowledgment

The authors acknowledge Professor H. Ohsaki, and Professor Y. Terao, (at The University of Tokyo) and Professor S. Shuchi, and Professor Y. Kataoka, (at Akita Prefectural University) for their technical advice. The author sincerely appreciates IHI Corporation, which supported this study technically. The authors are also grateful to Mr. H. Oyori, and Mr. K. Takahashi, (at IHI Corporation) and Ms. A. Mamada, (at The University of Tokyo) for their support.

## Funding Data

New Energy and Industrial Technology Development Organization (NEDO) (Project No. JPNP15005; Funder ID: 10.13039/501100003051).

## Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

## Nomenclature

*a*=height of each channel for each air type (m)

*b*=width of each channel for each air type (m)

*c, C*=constants

*C*_{P},*C*_{V}=isobaric and isochoric specific heats (J/(kg·K))

*d*=diameter or hydraulic diameter (m)

*D*=drag force exerted on a flying airplane (N)

*E*_{pro},*E*_{kin}=propulsion energy, kinetic energy (W)

*F*=thrust generated by the propulsion system (N)

*h*=heat transfer coefficient (W/(m

^{2}·K))*K*=overall heat transfer coefficient (W/(m

^{2}·K))- $\u2113$ =
length of heat exchanger or rotor in the axial direction (m)

*m*=mass (kg)

- $m\u02d9$ =
mass flow rate (kg/s)

*n*=number of channels for each air type

*N*=number of heat transfer units

- Nu
_{CL}, Nu_{P}= clearance Nusselt number, pipe Nusselt number

*p*=pressure (Pa)

*P*=peripheral length of each channel for each air type (m)

- Pr =
Prandtl number

*q*=heat flux (W/m

^{2})*Q*=heat flow rate (W)

*R*, $\u211c$ =gas constant (J/(kg·K)) and universal gas constant (J/(mol·K))

- Re
_{A}, Re_{P}, Re_{R}= axial, pipe, and rotational Reynolds numbers

*S*=cross section of each channel for each air type (m

^{2})*t*=time (s)

*T*=temperature (K)

*u*=average jet velocity from the propulsion system under the main wings (m/s)

*U*=tangential velocity of the inner surface of the rotor (m/s)

*v*=flight velocity of an airplane (m/s)

*V*=axial velocity (m/s)

- $V\u02d9$ =
volume flow rate (m

^{3}/s)*w*=local air velocity after a flying airplane (m/s)

*W*=power consumption (W)

*x*=position (m)

*X*=rotation rate in rotations per minute (rpm)

*α, β, γ, ζ, ξ, σ, ϕ, ψ, θ*=combined variables

*γ*=specific heat ratio

*δ*=_{B}thickness of the bulkhead (m)

*δ*=_{C}diameter difference between coaxial cylinders (m)

*δ*_{TBC,}_{R},*δ*_{TBC,}_{S}=thermal barrier coating thickness for the rotor and the stator (m)

*ε*=ratio of heat capacity flow rates

*η*_{comp}=compressor efficiency

*η*_{pro}=propulsion efficiency

*λ*=thermal conductivity (W/(m·K))

*μ*=viscosity (Pa·s)

*ρ*=density (kg/m

^{3})*ω*=rotation rate in radians per second (rad/s)

### Subscripts

- A =
shaft

- Al =
aluminum

- B =
bulkhead between internal air and external air

- E =
external air

- H =
heat exchanger

- I =
internal air

- in =
inlet of heat exchanger

- L =
laterad (outside)

- M =
medial (inside)

- mol =
per one mol

- out =
outlet of heat exchanger

- PM =
permanent magnet

- R =
rotor

- S =
stator

- SC =
coil in stator

- t =
time

- TBC =
thermal barrier coating

- x =
position from the left end of each heat exchanger or rotor

### Function

### Appendix A

##### A1. Radiative Heat Transfer Rate.

*T*

_{S}= 250 °C, a PM temperature

*T*

_{R}= 100 °C, a temperature at the inner surface of the TBCs on the coil

*T*

_{TBC,}

_{S}= 190 °C, and a temperature at the outer surface of TBC on the rotor

*T*

_{TBC,}

_{R}= 160 °C, the radiative heat transfer

*Q*

_{rad,}

_{S}from the coils and that

*Q*

_{rad,}

_{R}from the rotor via the annular air layer are as follows, respectively:

*σ*

_{SB}is Stefan–Boltzmann constant of 5.67 × 10

^{−8}(W/(m

^{2}K

^{4})),

*d*

_{SM}is the inner diameter of the coils,

*d*

_{RL}and

*ℓ*

_{R}are the outer diameter and the length of the rotor, and

*ε*

_{S}and

*ε*

_{R}are emissivities of TBCs on the coils and rotor which are both assumed to be 0.5. Solving the two equations simultaneously

##### A2. Validation of the Proposed Method.

The proposed theoretical estimation used several models, such as a heat transfer coefficient for each flow scheme in a channel with solid walls, which all were already validated experimentally in the original manuscripts. In addition, a design process for an air-to-air heat exchanger is very common in the thermal engineering field. However, it is important for the proposed theoretical estimation to be validated experimentally or numerically. The authors and colleagues are now preparing validation experiments for the targeted IPM motor with the air cooling designed by the proposed theoretical estimation described in this paper; however, some more time is still needed. Instead, for validation of the proposed theoretical estimation, numerical results of an air-cooled IPM motor with the configuration shown in Table 2 by Jercic et al. [29] were used. They showed the steady-state temperature distribution of PMs and coils as shown in Table 3 when the minimum cooling-air flowrate is 0.096 m^{3}/s at a rotation speed of 1700 rpm.

Length of rotor and stator | 327 mm | Inner dia. of external-air channels | 228 mm |

Shaft diameter | 59 mm | Outer dia. of external-air channels | 242.5 mm |

Inner dia. of PMs | 91 mm | Outer dia. of heat exchanger | 262.5 mm |

Outer dia. of PM | 139 mm | Height of external-air channels | 7.5 mm |

Inner dia. of coils | 142 mm | Width of external-air channels | 7.5 mm |

Outer dia. of coils | 190 mm | Number of external-air channels | 36 |

Inner dia. of heat exchanger | 222.5 mm | All values were read and estimated from figures in Ref. [29] by the authors |

Length of rotor and stator | 327 mm | Inner dia. of external-air channels | 228 mm |

Shaft diameter | 59 mm | Outer dia. of external-air channels | 242.5 mm |

Inner dia. of PMs | 91 mm | Outer dia. of heat exchanger | 262.5 mm |

Outer dia. of PM | 139 mm | Height of external-air channels | 7.5 mm |

Inner dia. of coils | 142 mm | Width of external-air channels | 7.5 mm |

Outer dia. of coils | 190 mm | Number of external-air channels | 36 |

Inner dia. of heat exchanger | 222.5 mm | All values were read and estimated from figures in Ref. [29] by the authors |

Power | 118.5 kW | Rotation rate | 1700 rpm |

Efficiency | 97.05% | Flowrate of external-air | 0.096 m^{3}/s |

Power | 118.5 kW | Rotation rate | 1700 rpm |

Efficiency | 97.05% | Flowrate of external-air | 0.096 m^{3}/s |

The heat generation rate should be 3495 W when the power is 118.5 kW and the efficiency is 97.05%. The authors assumed the heat generation rate from the PM is 218.5 W and that from the coils is 3277.25 W. The motor of Jercic et al. has a simple air-cooling scheme different from the authors; therefore, the authors assumed several settings to allow the proposed theoretical estimation to simulate the motor of Jercic et al. First, internal air is always aspirated from the surroundings instead of recirculation; thus, inlet temperatures at both the heat exchanger and the rotor are fixed at 295 K. Second, the authors set the flowrate of the internal air at 0.00177 m^{3}/s to keep the temperature of the PM at 147.6 °C because Jercic et al. did not describe how they cooled the rotor in Ref. [29]. Third, the authors set a heat transfer coefficient at contact surfaces of internal air in the outer heat exchanger at 5 W/(m^{2}⋅K) to simulate a natural convection cooling by the surrounding air (Table 4).

Heat generation rate from PM | 218.5 W | Heat generation rate from coils | 3277.25 W | ||||

Flowrate of internal air | 0.00177 m^{3}/s | Flowrate of external air | 0.096 m^{3}/s | ||||

Inlet temp. at the heat exchanger of internal air | 295 K | Inlet temp. at the heat exchanger of external air | 295 K | ||||

Inlet temp. at the rotor of internal air | 295 K | Heat transfer coefficient on contact surfaces of internal air | 5 W/m^{2} K |

Heat generation rate from PM | 218.5 W | Heat generation rate from coils | 3277.25 W | ||||

Flowrate of internal air | 0.00177 m^{3}/s | Flowrate of external air | 0.096 m^{3}/s | ||||

Inlet temp. at the heat exchanger of internal air | 295 K | Inlet temp. at the heat exchanger of external air | 295 K | ||||

Inlet temp. at the rotor of internal air | 295 K | Heat transfer coefficient on contact surfaces of internal air | 5 W/m^{2} K |

Using the above assumptions, the proposed theoretical method in this paper estimates the temperature of coils *T*_{S} = 138.2 °C in Table 5. The numerical analysis by Jercic et al. in Fig. 8 in Ref. [29] showed temperatures of the coils at seven points, whose average value was *T*_{S} = 141.1 °C

Jercic et al. [29] | Proposed theoretical estimation | |
---|---|---|

Temperature of PM | T_{R} = 147.6 °C | T_{R} = 147.6 °C |

Temperature of coils | 116.7 °C, 154.5 °C, 141.2 °C, 141.5 °C, 145.8 °C, 133.6 °C, 154.5 °C, average T_{S} = 141.1 °C | T_{S} = 138.2 °C |

Jercic et al. [29] | Proposed theoretical estimation | |
---|---|---|

Temperature of PM | T_{R} = 147.6 °C | T_{R} = 147.6 °C |

Temperature of coils | 116.7 °C, 154.5 °C, 141.2 °C, 141.5 °C, 145.8 °C, 133.6 °C, 154.5 °C, average T_{S} = 141.1 °C | T_{S} = 138.2 °C |

The proposed theoretical method can estimate a reasonable result in the above case, although the authors must validate this theoretical method by some experiments in the near future.

## References

**2**(2), pp.

**12**(1).10.4271/14-12-01-0005