Numerical Analysis of Wick-Type Two-Phase Mechanically Pumped Fluid Loop for Thermal Control of Electric Aircraft Motors

: The development of thermal control systems has become an important issue in next-generation electric aircraft design due to the increase in heat exhausted with electriﬁcation. In this paper, a wick-type two-phase mechanically pumped ﬂuid loop system for future electric aircraft was proposed through the investigation of current two-phase ﬂow cooling technology. Taking the experimental electric aircraft X-57 as an example, a wick-type two-phase mechanically pumped ﬂuid loop with four evaporators for transporting 12 kW of waste heat within an 80 ◦ C temperature limit was proposed and its feasibility was conﬁrmed. A numerical model was constructed and validated to predict the operating characteristics of a two-phase mechanically pumped ﬂuid loop. The optimal pump outputs under-even and uneven heat load conditions and was investigated for the ﬁrst time by considering the vapor-liquid separation conditions in each ﬂow path and the power consumption of the pump. Under the optimal pump output condition, the operating characteristics of the wick-type two-phase mechanically pumped ﬂuid loop system were calculated. The calculation results indicate that the proposed wick-type two-phase mechanically pumped ﬂuid loop is suitable as the thermal control system for an X-57 electric aircraft motor, as the calculation results satisﬁed the operational requirements of the motor.


Introduction
Recently, the electrification of next-generation aircraft has become more important due to the need to reduce the emission of carbon dioxide.Several small electric aircraft have been proposed, and a flight experiment has also been conducted [1][2][3][4].Based on the knowledge obtained from small electric aircraft, a large passenger aircraft was proposed [5][6][7][8], and the application was planned to commence before 2035 to meet the requirements of the International Air Transport Association (IATA) roadmap [9].However, with electrification, the heat generation density of the battery, electric motors, and avionics is expected to increase rapidly due to the miniaturization of the aircraft body.As a result, it is difficult for the existing air-cooling method to exhaust all the waste heat to the surroundings.Currently, due to the temperature limit, the flight period of an experimental electric aircraft is confined to approximately one hour.Therefore, advanced thermal control systems for electrification aircraft are required [10][11][12][13][14].
A similar thermal management issue also exists in satellites and spacecraft.A robust, capillary-pumped, two-phase heat transfer loop heat pipe (LHP) was proposed for cooling the electronic components installed in satellites and spacecraft.This LHP does not require electric power and features a large heat transfer distance, high heat transfer capacity, and flexible configuration [15,16].In addition, the ground uses of LHPs for cooling engines and central processing units (CPU) have been proposed [17][18][19].However, it is difficult for the LHP to cool intense heat sources, such as the electric aircraft's motor, as a result of the insufficient driving force.To solve this problem, an additional driving force is required.
Instead of the capillary pump, a mechanical pump can be employed to cool components with high heat generation, owing to its effectiveness in other applications.A single-phase, mechanically pumped fluid loop was used as the heat transfer device in the Mars rover, Curiosity, to cool the waste heat from a radioisotope thermoelectric generator (nuclear battery) [20].The two-phase, mechanically pumped fluid loop (2PMPFL) was proposed as a future heat transfer device owing to the following characteristics: (1) better heat transfer performance compared to single-phase flow; and (2) lower power consumption due to the low mass flow rate in the 2PMPFL.
Two types of 2PMPFLs were mentioned by Daimaru, T. et al. [21].The first is the wickless-type 2PMPFL, which uses a mechanical pump to circulate the working fluid.Existing research on wickless-type 2PMPFLs includes studies conducted on the ground and on the Kibo module of the international space station, to investigate the two-phase flow boiling phenomena in the evaporator [22][23][24][25][26][27][28].In the wickless-type 2PMPFL, the two-phase flow boiling behavior in the evaporator is difficult to predict, given the lack of vapor-liquid separation.Poor flow stability results in an uneven temperature distribution in the evaporator surface.A more reliable 2PMPFL system is necessary due to an even temperature distribution across the evaporator is desired.
The second type of 2PMPFL, mentioned by Daimaru, T. et al. [21], is the wick-type 2PMPFL.In this type, a porous structure called wick is employed in the evaporator to achieve vapor-liquid separation.Compared to the wickless-type 2PMPFL, the wick-type 2PMPFL offers the following advantages.(1) The wick enables vapor-liquid separation in the evaporator, which enhances flow stability, especially in the wick-type 2PMPFL system including multiple evaporators.As a result, the operating characteristics are easier to predict.(2) The wick forms a uniform vapor-liquid boundary, creating an even temperature distribution in the evaporator.(3) The wick can provide additional capillary pump capability, alongside the mechanical pump capability, which improves the operating performance of the wick-type 2PMPFL.
Several experimental approaches have been conducted to investigate the heat transfer performance of wick-type 2PMPFLs with one evaporator [29][30][31][32][33][34][35], two evaporators [36,37], and six evaporators [38].However, a detailed design method and numerical model for the wick-type 2PMPFL was not established as a result of the complicated operating behavior.
In this paper, a wick-type 2PMPFL with multiple evaporators was proposed as the potential thermal control system for cooling multiple heat sources in future electric aircraft.Firstly, a one-dimensional steady-state numerical model for the proposed wick-type 2PMPFL cooling system with four evaporators was established to confirm whether it satisfies the thermal control requirements of the electric aircraft motor.The four heat sources (generator, converter, inverter, and rotor) in the motor were assumed to be cooled through four evaporators, and the surface of the wing was assumed to be the heat dissipation area.In the numerical model, for the first time, the optimal output condition of the pump (i.e., the optimal pump capacity) was deduced by considering the vapor-liquid separation in each flow path and the power consumption.Under the optimal conditions of the pump, the following attributes were investigated; the operating temperature of the evaporators and accumulators, and the mass flow rate in each flow path.

Wick-Type Two-Phase Mechanically Pumped Fluid loop 2.1. Operating Characteristics
Figure 1 presents an image of the proposed wick-type 2PMPFL thermal management system with four evaporators.The wick-type 2PMPFL absorbs waste heat from the powered machine and transfers the waste heat to a cold surrounding on the condenser side through the phase change of the inner working fluid (i.e., vapor to liquid).A mechanical pump is used to circulate the working fluid, which efficiently transports heat from the evaporator to the condenser.through the phase change of the inner working fluid (i.e., vapor to liquid).A mechanical pump is used to circulate the working fluid, which efficiently transports heat from the evaporator to the condenser.Evaporators with isothermal temperature variations were developed in the authors' previous research [39][40][41].A porous structure, called a wick, is presented in each evaporator.As a uniform vapor-liquid boundary can be generated in the wick, the temperature variation in the evaporating surface is small.As a result, the heat source can be cooled with an even temperature distribution heat transfer surface.Furthermore, the wick also contributes toward the flow stability of the working fluid in the evaporators by separating the vapor and liquid flows; this makes predicting the operating characteristics of the 2PMPFL less complicated.
There are two flow paths in the wick-type 2PMPFL.One is the main flow path, where the working fluid flows from the liquid line's divergence point C through the primary wick in each evaporator, the vapor line junction point A, the condenser, and then returns to point C. The other is the bypass flow path, where the working fluid flows from the liquid line's divergence point C through each compensation chamber (CC), the bypass line junction point B, and then returns to point C.An accumulator was installed near the mechanical pump to regulate the amount of working fluid in the fluid loop when the heat load varies.
During the operation of the wick-type 2PMPFL, the mass flow rate throughout the wick-type 2PMPFL is controlled by the mechanical pump.In the main flow path, the mass flow rate can be estimated using heat loads.After calculating the mass flow throughout the wick-type 2PMPFL and the main flow path, the mass flow rate in the bypass line can be calculated.Evaporators with isothermal temperature variations were developed in the authors' previous research [39][40][41].A porous structure, called a wick, is presented in each evaporator.As a uniform vapor-liquid boundary can be generated in the wick, the temperature variation in the evaporating surface is small.As a result, the heat source can be cooled with an even temperature distribution heat transfer surface.Furthermore, the wick also contributes toward the flow stability of the working fluid in the evaporators by separating the vapor and liquid flows; this makes predicting the operating characteristics of the 2PMPFL less complicated.
There are two flow paths in the wick-type 2PMPFL.One is the main flow path, where the working fluid flows from the liquid line's divergence point C through the primary wick in each evaporator, the vapor line junction point A, the condenser, and then returns to point C. The other is the bypass flow path, where the working fluid flows from the liquid line's divergence point C through each compensation chamber (CC), the bypass line junction point B, and then returns to point C.An accumulator was installed near the mechanical pump to regulate the amount of working fluid in the fluid loop when the heat load varies.
During the operation of the wick-type 2PMPFL, the mass flow rate throughout the wicktype 2PMPFL is controlled by the mechanical pump.In the main flow path, the mass flow rate can be estimated using heat loads.After calculating the mass flow throughout the wick-type 2PMPFL and the main flow path, the mass flow rate in the bypass line can be calculated.

Design Method
In this work, as the first step, a wick-type 2PMPFL was proposed to realize the thermal control of the Maxwell X-57 aircraft's cruise motor, which was developed by the ESAero corporation in 2017 [42].As the existing thermal control system is incomplete and the steady-state temperature is beyond the operating temperature limit, the current flight period of the X-57 is limited to one hour before the temperature exceeded the operating limit.
The power of the X-57 cruise motor is 60 kW, and the efficiency is assumed to be 80%, considering the design margin.Therefore, supposing all the losses are transformed to heat, the amount of heat generated is 12 kW.The heat generation system consists of four components: a generator, converter, inverter, and rotor.The maximum amount of heat generated by each component was estimated to be 3 kW.A part of the wing area (0.53 m × 3 m) at the top and bottom surfaces was used as the heat-dissipating surface.During the flight period, the air temperature was assumed to be 20 • C, and the flight speed was 76.9 m/s.The temperature limit of the cruise motor was set to 100 • C. Considering the design margin, the temperature limit of the evaporator was set to 80 • C.
Table 1 shows the configurations of the wick-type 2PMPFL based on Figure 1, where a 12.7 mm tube was used as the vapor line and the condensing line with an inner diameter of 10.9 mm, whereas a 6.35 mm tube was used as the liquid line with an inner diameter of 4.6 mm.The length of the condensing line in the main flow path/bypass flow path was 12.1 m/2.3 m.The length of the liquid line in the main flow path and bypass flow path was 1.83 m/1.5 m.The volume of the accumulator was 2.55 L. R245fa was selected as the working fluid considering its desirable properties: high critical temperature; a low global warming potential (more environmentally friendly); and excellent thermal and chemical stability.The charging amount was calculated to be 4.058 kg by considering the vapor-liquid distribution.

Modeling System
The mathematical model was established by referring to the LHP numerical model [43].For the wick-type 2PMPFL, zero-dimensional numerical models were used for the evaporator and CC accumulator, and one-dimensional models were used for the vapor line, condenser line, liquid line, and bypass line to predict the temperature and pressure of each component.
Equations ( 1) and ( 2) show the energy balance conditions in each evaporator.T evap and T ev represent the temperatures of the heating surface and the evaporating surface in wick, respectively.
where Q load is the total heat load, Q evap_amb is heat lost to the surroundings, and G evap_ev T evap − T ev is the heat transferred to the wick's evaporating surface.
where G wick (T ev − T cc ) is the heat lost to the CC side, and .mλ is the heat absorbed through evaporation of the working fluid.Equation (5) shows the energy balance condition in each CC.Q cc_amb represents the heat leak from the CC to the surroundings, and T ll_out is the liquid line outlet temperature.
The temperature variation and pressure variation in the vapor, condenser, liquid, and bypass line were calculated using Equations ( 6) and (7), respectively.T f and T a represent the working fluid's temperature and the temperature of the surroundings, respectively; G f_a represents the thermal conductance between the working fluid and the surroundings.
Müller-Steinhagen and Heck's equations (Equations ( 8) and ( 9)) were used to calculate the pressure drop in the two-phase flow region, and Fujii's equation (Equations ( 10)-( 17)) was used to calculate the Nusselt number over the inner surface of the condensing tube in the two-phase flow region [44,45].
Nu F = 0.018Re 0.9 Nu FG_G = 0.725 Energies 2022, 15, 1800 6 of 15 In Equations ( 8) and ( 9), dP dL f,2ph represents the gradient of the two-phase flow pressure drop, dP dL f,v represents the gradient of the pressure drop assuming that the entire working fluid is in the vapor phase, and dP dL f,l represents the gradient of pressure drop assuming that the entire working fluid is in the liquid phase.
In Equations ( 10)-( 17), Nu F is the Nusselt number in the forced convection condensation region, and Nu FG_G is the Nusselt number from the co-existing forced and natural convection condensation region to the end of the natural convection condensation region.The larger value between Nu F and Nu FG_G was used as the average Nusselt number for the heat transfer between the working fluid and the inner wall of the condensing tube.
Equation ( 18) was employed to calculate the heat transfer coefficient between the heat dissipation surface (aircraft wing) and the surroundings.u ∞ represents the flight speed.The wing was assumed to be cooled through turbulent airflow.

Calculation Procedure
Figure 2 shows a flowchart describing the operation of the control system.First, the accumulator temperature T acc1 and heat load conditions at each evaporator were inputted.Unlike an LHP, in the wick-type 2PMPFL, the heat leaked to the accumulator is negligible; therefore, the accumulator is at the surrounding temperature during operation.Assuming that the working fluid in the accumulator and the condensing region operate in a vaporliquid saturation state, the temperature and pressure in the accumulator and the condensing region satisfy the Clausius-Clapeyron equation (Equation ( 19)).If the pressure difference between the accumulator and the condenser is small, the temperature difference of the working fluid between the accumulator and the condenser is also small; this results in incomplete condensation and ceases the operation.To ensure complete condensation, the accumulator must be heated to increase the temperature of the working fluid in the condenser.Therefore, a sufficient temperature difference between the working fluid and condenser wall is generated for condensation.The calculation results present the wick-type 2PMPFL's operating characteristics when the accumulator is operated at the lowest permitted temperature for completing condensation under each corresponding heat load condition.
After inputting the accumulator temperature T acc1 and consequently the heat load conditions at each evaporator, the initial value of the pump output, evaporating temperature in the evaporator, and heat leak from the evaporator side to the CC side were set.Next, the mass flow rate in the main and bypass flow paths and the CC temperature T cc1 were calculated using the initial and input values.Thereafter, the temperature and pressure values in the vapor, condenser, liquid, and bypass lines were sequentially ordered along the flow direction of the working fluid.The accumulator is saturated with vapor-liquid, which allows for the calculation of a second temperature, T acc2 , using the saturation vapor curve, provided the pressure value of the accumulator is available.The calculation result for the accumulator's temperature is acquired by adjusting the initial value, T acc1 , until the difference between T acc1 and T acc2 is less than 0.1 • C.After inputting the accumulator temperature T acc and consequently the heat load conditions at each evaporator, the initial value of the pump output, evaporating temperature in the evaporator, and heat leak from the evaporator side to the CC side were set.Next, the mass flow rate in the main and bypass flow paths and the CC temperature T cc were calculated using the initial and input values.Thereafter, the temperature and pressure values in the vapor, condenser, liquid, and bypass lines were sequentially ordered along the flow direction of the working fluid.The accumulator is saturated with vaporliquid, which allows for the calculation of a second temperature, T acc , using the saturation vapor curve, provided the pressure value of the accumulator is available.The calculation result for the accumulator's temperature is acquired by adjusting the initial value, T acc , until the difference between T acc and T acc is less than 0.1 C.
Subsequently, the energy equation for the CC (equation ( 5)) is used to calculate another CC temperature, T cc .The calculation result for the CC's temperature T cc is acquired by following the same steps as those for the accumulator's temperature.Subsequently, the energy equation for the CC (Equation ( 5)) is used to calculate another CC temperature, T cc2 .The calculation result for the CC's temperature T cc2 is acquired by following the same steps as those for the accumulator's temperature.
The mass flow rate in the vapor line was calculated using the corresponding heat load condition, and the mass flow rate in the bypass line was calculated using the pump output and the mass flow rate in the vapor line.The optimal output of the pump was calculated using the power consumption and the vapor-liquid separation condition.The conditions for the optimal output are satisfied when only the vapor exists in the vapor line and only liquid exists in the bypass line.The mass flow rate in the main and bypass flow paths should follow the relationship that the pressure drop of the main flow path is equal to the pressure drop of the bypass flow path.In the numerical model, the optimal pump output can be calculated by adjusting the initial value of the pump output until the pressure difference between the main and bypass flow paths is less than 10 Pa.Under optimal operation, the pump consumes the minimum power required, circulates the working fluid, and realizes the separation of vapor and liquid in the evaporator.If the pressure drop of the main flow path is greater than that in the bypass flow path, a part of the vapor penetrates the CC and flows through the bypass flow path, resulting in a high operating temperature and a deteriorating operating state.By contrast, if the pressure drop of the bypass flow path is greater, some of the liquid will flow out of the wick with the vapor.If this occurs, the pump consumes more power than required, and the separation of the vapor and liquid flows in the evaporator cannot be realized.
According to Figure 1, assuming that the mass flow rates of the returning liquid to each CC, in the main flow path of each evaporator, and in the bypass flow path of each evaporator are  20)-( 22), provided all the evaporators are operated under even heat load conditions.Under uneven heat load conditions, assuming that evaporator 1 is heated with the highest heat load and the initial values of the temperature and mass flow rate, the initial value of the evaporating temperature in the other evaporators is adjusted until the difference in the pressure calculation results at point A through the flow paths with . .

Calculation Conditions
The flow of the working fluid was approximated to be a one-dimensional incompressible viscous flow.Referring to the National Institute of Standards and Technology (NIST) reference fluid properties [46], the thermophysical properties of the working fluid were used in the functions of temperature in the form of a polynomial.Uniform temperature distribution was assumed at the heating surface of the evaporator and the heat dissipation surface of the condenser.The only liquid flow was considered to exist in the CC and bypass flow path.

Calculation Results
The validation of the numerical model for two types of two-evaporators LHP was conducted at first.When comparing the operating characteristics of two-evaporators LHP with the wick-type 2PMPFL with four evaporators, the only difference is that in the twoevaporators LHP, the separation of vapor and liquid flow is realized in the evaporator as a result of the wick, while in the wick-type 2PMPFL with four evaporators, the vapor and liquid flow cannot be realized unless the optimal pump capacity is attained.Therefore, under the optimal pump capacity condition (i.e., the vapor-liquid separation condition in the wick-type 2PMPFL system's evaporator), the operating characteristics of the twoevaporators LHP and the wick-type 2PMPFL with four evaporators are identical.If the numerical model is validated in the two-evaporators LHP, the validation of the numerical model for wick-type 2PMPFL with four evaporators under the optimal pump capacity condition will also be realized.
Figure 3 presents the comparison between calculation and the experimental results of one LHP with two evaporators (wick material: polytetrafluoroethylene, working fluid: acetone).Figure 4 presents the comparison between the calculations and the experimental results of another LHP with two evaporators (wick material: stainless steel 316L, working fluid: ammonia).The detailed configurations of both the two-evaporators LHPs were described in the references [47,48].In Figures 3 and 4, the x-axis shows the heat load at each evaporator, while the y-axis shows the calculations and the experimental results of each Energies 2022, 15, 1800 9 of 15 evaporator case's steady-state temperature.Since the temperature difference between the calculations and the experimental results was small in both types of two-evaporators LHP, the validation of the numerical model for two-evaporators LHP was realized.Therefore, under the optimal pump capacity condition, the validation of the wick-type 2PMPFL was realized.
Figure 3 presents the comparison between calculation and the experimental results of one LHP with two evaporators (wick material: polytetrafluoroethylene, working fluid: acetone).Figure 4 presents the comparison between the calculations and the experimental results of another LHP with two evaporators (wick material: stainless steel 316L, working fluid: ammonia).The detailed configurations of both the two-evaporators LHPs were described in the references [47,48].In Figures 3 and 4, the x-axis shows the heat load at each evaporator, while the y-axis shows the calculations and the experimental results of each evaporator case's steady-state temperature.Since the temperature difference between the calculations and the experimental results was small in both types of two-evaporators LHP, the validation of the numerical model for two-evaporators LHP was realized.Therefore, under the optimal pump capacity condition, the validation of the wick-type 2PMPFL was realized.

Figure 3.
Calculations and experiment results of the two-evaporators type LHP's each evaporator case's temperature according to the reference [47] (wick material: polytetrafluoroethylene, working fluid: acetone).The following calculations of the wick-type 2PMPFL with four evaporators were performed under the assumption that the pump operated at its optimal output.The factors investigated were the mass flow rate in each flow path, the minimum operating temperature of the accumulator, and the operating behavior of the evaporators and CCs. Figure 5 shows the calculation results for the evaporator temperature and the lowest accumulator temperature under even heat load conditions.The heat load for each evaporator was increased from 0.3 kW to 3.6 kW with an increment of 0.3 kW.The wick-type 2PMPFL operated until the heat loads were increased to 3.6 kW, where the operating temperature of each evaporator was higher than the 80 C limit.The wick-type 2PMPFL transferred a total of 13.2 kW (3.3 kW for each evaporator), which satisfies the design requirement of a total heat transfer capability of 12 kW (3 kW for each evaporator).The following calculations of the wick-type 2PMPFL with four evaporators were performed under the assumption that the pump operated at its optimal output.The factors investigated were the mass flow rate in each flow path, the minimum operating temperature of the accumulator, and the operating behavior of the evaporators and CCs. Figure 5 shows the calculation results for the evaporator temperature and the lowest accumulator temperature under even heat load conditions.The heat load for each evaporator was increased from 0.3 kW to 3.6 kW with an increment of 0.3 kW.The wick-type 2PMPFL operated until the heat loads were increased to 3.6 kW, where the operating temperature of each evaporator was higher than the 80 • C limit.The wick-type 2PMPFL transferred a total of 13.2 kW (3.3 kW for each evaporator), which satisfies the design requirement of a total heat transfer capability of 12 kW (3 kW for each evaporator).
temperature under even heat load conditions.The heat load for each evaporator was increased from 0.3 kW to 3.6 kW with an increment of 0.3 kW.The wick-type 2PMPFL operated until the heat loads were increased to 3.6 kW, where the operating temperature of each evaporator was higher than the 80 C limit.The wick-type 2PMPFL transferred a total of 13.2 kW (3.3 kW for each evaporator), which satisfies the design requirement of a total heat transfer capability of 12 kW (3 kW for each evaporator).Following the same load conditions as those in Figures 6 and 7 shows the ratio of the mass flow rate in the main to the bypass flow path, as well as the subcooled length in the condenser.As the four evaporators received the same heat loads, the ratio of the mass flow rate in the main flow path to the mass flow rate in the bypass flow path was the same for each evaporator, as shown in Equation (23).Furthermore, as indicated in Figure 6, the ratio of the mass flow rate in the main flow loop to the mass flow rate in the bypass flow loop increased with the heat loads, whereas the subcooled region length decreased, except in the fluctuating regions.This was attributed to the rough calculation with 5 • C temperature variations in each step in the accumulator temperature.With a more accurate calculation of the accumulator's temperature (i.e., smaller temperature variations in each step), the fluctuating regions will disappear.m b4 (23) Figure 7 shows the head (Equation ( 24)) and capacity (Equation ( 25)) under each corresponding heat load condition when the optimal pump output was realized.The mechanical pump, which satisfies both requirements, will be selected in the future.
Capacity = 60 .m ρ × 1000 (25) ature variations in each step in the accumulator temperature.With a more accurate calculation of the accumulator's temperature (i.e., smaller temperature variations in each step), the fluctuating regions will disappear.
Figure 6.Calculation results of mass flow ratio and subcooled region length under even heat loads conditions.
Figure 7 shows the head (Equation ( 24)) and capacity (Equation ( 25)) under each corresponding heat load condition when the optimal pump output was realized.The mechanical pump, which satisfies both requirements, will be selected in the future.Table 2 shows each evaporator's temperature and the percentage of mass flow rate in each returning liquid line (m c , m c , m c , m c ), each vapor line (m a , m a , m a , m a ), and each bypass line (m b , m b , m b , m b ) under the uneven heat load conditions, where the heat loads were assumed to be 3 kW, 2.7 kW, 2.4 kW, and 2.1 kW for evaporators 1, 2, 3, and 4, respectively.Therefore, temperature control for each evaporator was feasible, even though all the evaporators operated with different heat loads under the optimal pump output conditions.As shown in Table 2, the mass flow rate of the returning liquid CC 1, m c , was the highest, contrary to m b in the bypass flow path with the smallest value.Table 2 shows each evaporator's temperature and the percentage of mass flow rate in each returning liquid line (m c1 , m c2 , m c3 , m c4 ), each vapor line (m a1 , m a2 , m a3 , m a4 ), and each bypass line (m b1 , m b2 , m b3 , m b4 ) under the uneven heat load conditions, where the heat loads were assumed to be 3 kW, 2.7 kW, 2.4 kW, and 2.1 kW for evaporators 1, 2, 3, and 4, respectively.Therefore, temperature control for each evaporator was feasible, even though all the evaporators operated with different heat loads under the optimal pump output conditions.As shown in Table 2, the mass flow rate of the returning liquid CC 1, .m c1 , was the highest, contrary to .m b1 in the bypass flow path with the smallest value.

Conclusions
In this paper, to solve the thermal control system of future electric aircraft, a wicktype 2PMPFL system with four evaporators was proposed for cooling the cruise motor of the X-57 electric aircraft.A numerical model was established to predict the operating characteristics of the wick-type 2PMPFL system.For the first time, the pump's optimal output condition based on mass flow rate was investigated, which enabled the selection of an appropriate pump according to the mass flow rate calculated from the heat load conditions.Based on the validation of the numerical model for the two-evaporators LHP, the numerical model for wick-type 2PMPFL with four evaporators under the optimal pump output condition was validated.Under the even heat load conditions, the wick-type 2PMPFL system can transfer 13.2 kW (3.3 kW for each evaporator) within the temperature upper limit of 80 • C, which satisfied the design goal of 12 kW (3 kW for each evaporator).Under the uneven heat load conditions with the optimal pump output, it is found that the vapor-liquid separation was realized in four evaporators simultaneously, and the temperature control for each evaporator was feasible.These achievements in the proposal of a wick-type 2PMPFL system will contribute toward establishing an effective future thermal control system for electric aircraft.

Figure 1 .
Figure 1.Schematic of a wick-type 2PMPFL with four evaporators (not to scale).

Figure 1 .
Figure 1.Schematic of a wick-type 2PMPFL with four evaporators (not to scale).
and (4) represent the thermal conductance of the wick.d wick and A wick represent the thickness of the wick and the cross-sectional area of the wick, respectively.

Figure 2 .
Figure 2. Flow chart of control system operation.

Start=
Set evaporating temperature, Heat leak from evaporator to CC, total pump output, decide Calculate mass flow rate, Calculate pressure drop and temperature variation in vapor line, condenser, liquid line and bypass line Calculate from local energy balance equation Yes

Figure 2 .
Figure 2. Flow chart of control system operation.
respectively, the relation between the mass flow rates in each flow path satisfies Equations ( less than 10 Pa.Subsequently, the initial value of the mass flow rate is adjusted on each evaporator's bypass line until the difference in the pressure calculation results at point B through the flow paths with .

Figure 4 .
Figure 4. Calculations and experiment results of the two-evaporators type LHP's each evaporator case's temperature according to the reference [48] (wick material: stainless steel 316L, working fluid: ammonia).

Figure 4 .
Figure 4. Calculations and experiment results of the two-evaporators type LHP's each evaporator case's temperature according to the reference [48] (wick material: stainless steel 316L, working fluid: ammonia).

Figure 5 .Figure 5 .
Figure 5. Calculation results of the evaporator temperature and accumulator's lowest temperature under even heat loads conditions.

Figure 7 .
Figure 7. Calculation results for pump capacity and head requirements under even heat loads conditions.
dP dL f,v Gradient of the pressure drop assuming that the entire working fluid is in the vapor phase, (Pa•m −1 ) dP dL f,2ph Gradient of the two-phase flow pressure drop, (Pa•m −1 ) dP dT sat Gradient of the pressure drop, (Pa•K −1 ) dP total Pressure drop in the whole loop, (Pa) dT Temperature variation, (K) d wick Thickness of the wick, (m) f Darcy friction factor, (-) G Parameter for calculation, (Pa•m −1 ) Ga D Galilei number, (-) g Gravitational acceleration, (m•s −2 ) G evap_ev Thermal conductance between evaporator case and evaporating surface, (W•K −1 ) G f_a Thermal conductance between working fluid and surroundings, (W•K −1 ) G wick Thermal conductance of wick, (W•K −1 ) h Heat transfer coefficient, (W•m −2 •K −1 ) H L Condensation number, (-) k air Thermal conductivity of air, (W•m −1 •K −1 ) k wick Thermal conductivity of wick, (W•m −1 •K −1 ) .m Mass flow rate, (kg•s −1 ) .m a1 Mass flow rate in the main flow path of evaporator 1, (kg•s −1 ) .m a2 Mass flow rate in the main flow path of evaporator 2, (kg•s −1 ) .m a3 Mass flow rate in the main flow path of evaporator 3, (kg•s −1 ) .m a4 Mass flow rate in the main flow path of evaporator 4, (kg•s −1 ) .m b1 Mass flow rate in the bypass flow path of evaporator 1, (kg•s −1 ) .m b2 Mass flow rate in the bypass flow path of evaporator 2, (kg•s −1 ) .m b3 Mass flow rate in the bypass flow path of evaporator 3, (kg•s −1 ) .m b4 Mass flow rate in the bypass flow path of evaporator 4, (kg•s −1 ) .m c1 Mass flow rate of the returning liquid to CC 1, (kg•s −1 ) .m c2 Mass flow rate of the returning liquid to CC 2, (kg•s −1 ) .m c3 Mass flow rate of the returning liquid to CC 3, (kg•s −1 ) .m c4 Mass flow rate of the returning liquid to CC 4, (kg•s −1 ) .m L Mass flow rate of liquid, (kg•s −1 ) Nu F Nusselt number in the forced convection condensation region, (-) Nu FG_G Nusselt number from the co-existing forced and natural convection condensation region to the end of the natural convection condensation region, (-) P acc Pressure of accumulator, (Pa) P con Pressure of condenser, (Pa) Pr Prandtl number, (-) Pr L Prandtl number of liquid, (-) Q cc_amb Heat exchange amount from CC to surroundings, (W) Q evap_amb Heat exchange amount from evaporator to surroundings, line outlet, (K) T w Temperature of wing, (K) u Velocity, (m•s −1 ) u ∞ Flight speed, (m•s −1 ) x Vapor quality, (-) η Parameter for calculation, (-) λ Latent heat, (J•kg −1 )

Table 2 .
Calculation results of operating temperature and the mass flow rate in each flow path for each evaporator.

Table 2 .
Calculation results of operating temperature and the mass flow rate in each flow path for each evaporator.