Experimental and Numerical Investigation on an Integrated Thermal Management System for the Li-Ion Battery Module with Phase Change Material

Lightweight power battery modules with outstanding thermal performance are urgently required given the rapid development of electric vehicles. This study proposes a composite phase change material coupled with forced convection as an integrated thermal management system (ITMS) with the aim to control the temperature’s rising tendency and maintain the temperature distribution uniformly within an appropriate range among the battery modules. The thermal behavior effects of airflow rates on the thermal management system were investigated in detail by combining experiments and numerical simulations. Comparisons were conducted between an air cooling system with an optimum flow rate and the ITMS. Experimental results revealed that the cooling effect of the ITMS was better than that of the forced cooling system at a 3m/s airflow rate. The maximum temperature in the designed battery module was limited to 63.2°C. The maximum temperature difference was limited to 4.8°C at a 4 C discharge rate. This research indicates that the ITMS is an effective and optimized approach to control and balance the temperature among battery modules, thereby providing engineers with design optimization strategies for similar systems.


Introduction
Lithium-ion (Li-ion) batteries have been broadly applied in electric vehicles (EVs) owing to their advantages, such as low discharge rate, high capacity, and long lifespan. Heat could be generated during the chemical reaction and ohmic resistance process and lead to temperature rise. However, the performance of Li-ion batteries is easily affected by heat. Hence, the temperature must be controlled within a safe range. A temperature range of −10°C to 50°C is generally recognized as a tolerable operating condition; the recommended temperature range for optimal Li-ion battery performance is 20°C-40°C [1][2][3]. When heat accumulates to a certain degree inside a battery, the temperature will rapidly rise, and even thermal runaway will occur [4]. The large-scale battery module is restricted by heat dissipation technology in commercial development. Thus, an effective thermal management technology that can guarantee a safe temperature range for battery modules and keep the temperature difference as small as possible must be explored [5].
Many types of battery thermal management systems (BTMS) have been reported in the literature in recent years. The types of BTMS include the air cooling system [6,7], liquid cooling system [8][9][10][11], phase change material (PCM) cooling system [12][13][14], and heat pipe cooling system [15][16][17]. The air cooling system, which is one of the earliest forms of cooling technology, has supposed to be the simplest and cheapest approach, but its heat transfer coefficient and specific heat capacity of air are much lower than those of the other cooling systems; thus, efficiently dissipating heat and suppressing the temperature difference at a high discharge rate is difficult, especially under certain extreme conditions [18]. Wu et al. [19] conducted experiments and simulations and reported that the heat generated inside the battery pack could not be efficiently dissipated through natural convection cooling. By contrast, forced convection can satisfy the requirement for controlling temperature rise for the battery. Park and Jung [20] focused on the arrangement of batteries and parasitic power consumption on the thermal behavior of the air cooling system. The numerical simulation results revealed that a battery module with a small distance for batteries was desirable for the air cooling system. However, the consumed power of the air cooling system for large heat load conditions would become much larger than that of the liquid cooling system for thermal management. An ideal cooling performance of the air cooling system could be obtained by optimizing the design of the heat-dissipating structure. However, a certain degree temperature difference still remains, especially at a high discharge rate. The liquid cooling system has attracted much attention due to its high cooling efficiency and better temperature uniformity than the air cooling system. With regard to thermal battery management, the liquid cooling system can manage the battery temperature in the module within a normal operating range and maintain the temperature difference within 3°C or lower [21]. Vita et al. [22] compared the air and liquid cooling systems, measuring the prismatic Li-ion phosphate battery at different discharge rates under constant current. Irreversible and entropy-based heat generation was also taken into account for mathematical modeling, which was aimed at simulating the temperature evolution and distribution for a battery module. The aforementioned authors further summarized the pros and cons of the two cooling approaches. Huo et al. [23] proposed the cooling system for thermal management of the rectangular Li-ion battery assisted with minichannel cold plates. Various factors, such as flow direction, number of channels, flow rate of coolant mass, and environmental temperature, were investigated in depth by the numerical simulation method. The same research group explored a series of works on the liquid cooling system assisted with minichannel cold plates, which showed a relatively high cooling efficiency [24,25]. In spite of the higher cooling effect of the liquid cooling system, several problems, including high cost, leakage, and difficult sustainment, must be overcome.
The PCM, a novel thermal management system in actual application, has received great attention and has been extensively explored due to its simple structure, no power consumption, and high latent capacity [26]. When the battery is discharging, PCM can absorb the generated heat and store it as latent heat. A large portion of heat is utilized during charge, and a small part of it is transferred to the operating environments. When the battery temperature drops below the melting point or is left to shelve, PCM releases the stored heat to the battery module, thus effectively keeping the battery module at a relatively constant temperature and ensuring an excellent temperature difference. Hallaj and Selman [27] investigated the PCM on a battery module and found that the battery temperature was substantially lower and more uniform than the air cooling system at different discharge rates. The PCM cooling system demonstrated many advantages for the battery module, such as effectively prolonging the cycling life [28], decreasing thermal runaway [29], and sustaining a relatively high operating temperature [30]. Many studies have been conducted to improve the thermal conductivity of BTMS by using various methods, such as adding paraffin (PA) with metal foam/mesh [31][32][33], metal fins/metallic particle [34], graphene [35,36], carbon nanotubes [37], and expanded graphite (EG) [38]. Ling et al. [39] found that the factors of the packing density of the composite PCM and mass fraction PA affected the thermal performance. The result verified that the PA/EG composite displayed a small temperature rise rate, proper operating temperature, and uniform temperature distribution for the battery module. The composite PCM with a paraffin density of 890 kg/m 3 and a mass fraction of 75% could demonstrate optimal thermal management performance. Fathabadi [40] presented a hybrid active-passive system with PA/EG as the composite PCM. The results demonstrated that the temperature rise and difference decrease with the increase in EG mass fraction. Wu et al. [41] designed a 2D thermal model associated with PA/EG composites for battery thermal management. The effects of different EG mass fractions associated with thermal conductivity and different phase change enthalpy were observed at high discharge rates. The investigations revealed that the EG mass fraction of 15%-20% was the optimal proportion to improve the thermal performance of battery modules.
With the increasing heat generation and power of the scale-up battery module, the heat dissipation requirement of batteries for single thermal management has become difficult to achieve. Thus, it is very necessary to apply the integrated thermal management system (ITMS) to address the battery's thermal safety. Several ITMS types, such as the heat pipe cooling/liquid cooling system, PCM/heat pipe cooling system, and air cooling/PCM cooling system, are currently available. Wu et al. [42] proposed a heat pipe-assisted PCM-based cooling system. The experimental results indicated that the maximum temperature of the battery module could be controlled at less than 50°C at a 5 C discharge rate. A low and stable temperature fluctuation could be obtained under several cycling conditions. Jiang et al. [43] designed a tube shell for a battery module with composite PCM coupled with the forced air cooling system. The results indicated that the ITMS could greatly reduce the temperature rise and maintain the temperature difference within 1-2°C. Lazrak et al. [44] studied PCM by combining ITMS with experiments and simulations, which found that the ITMS could uniformly distribute around the battery and reduce the temperature.
The performance of the battery is sensitive to temperature. Meanwhile, the effectiveness of the BTMS can be greatly improved, especially at high discharge rates. Thus, a feasible and simple thermal management system with multiple cooling approaches for large-scale battery modules should be developed. In this work, an ITMS with PCM coupled with forced air cooling has been proposed for 13 (Series) * 4 (Parallel) battery modules composed of 52 cylindrical batteries. By means of coupled air cooling with PCM, the problem of temperature difference in battery modules with air cooling can be effectively solved. At the same time, the PCM has a higher thermal conductivity and greater heat dissipation area, which effectively reduces the external thermal resistance 2 International Journal of Photoenergy of the battery module and greatly improves the heat dissipation efficiency of the battery module. Single forced air cooling for batteries was first investigated through experiments and simulations at 1, 2, and 3 m/s airflow rates. The ITMS was then used for the battery modules with different current rates.
Comparisons between ITMS and forced convection (FC) at the 3 m/s airflow rate were analyzed. The designed ITMS could improve the cooling effectiveness of battery modules, especially on a large scale.

Experiment
In order to study the battery thermal management system with PCM coupled with air cooling, form-stable composite PCM modules were prepared, and then the temperature characteristics of the battery modules with and without PCM were constructed and compared. Finally, experimental data are used to verify the numerical simulation error, so as to improve the accuracy of subsequent experiments. Figure 1 shows the raw materials for the composite PCM preparation. These materials include paraffin, EG, and epoxy resin (ER). The EG (average particle size: 150 μm; fixed carbon content: ≥99%) was obtained from Qingdao Nansha Development Graphite Co., Ltd. Paraffin was purchased from Shanghai Joule Wax Co., Ltd. (carbon content: ≥99.99%; melting point: 50°C; and thermal conductivity: 0.24 W/(m K)). ER, a matrix material for enhancing the machinability of the composite PCM, was provided by Hunan Baxiongdi Co., Ltd. The ER (E-44, 6061) and curing agent (molecular weight: 650) mass ratio was 1 : 1. The expandable graphite was first dried at 60°C for 10 h in a vacuum oven and placed into a muffle furnace. Finally, expandable graphite was subjected to heat treatment at 900°C for 7.5 min to obtain the EG. Predetermined amounts of paraffin were maintained at a constant temperature of 70°C for 45 min in a bath until they completely melted. The paraffin was then mixed with the aforementioned EG under stirring at 2000 rad min -1 for 15 min. Subsequently, ER was added to the mixture under magnetic stirring at 2000 radmin -1 for 15 min. Curing agent B was then placed into a mixture, which was continuously stirred at 2500 rad min -1 for 5 min. Finally, the obtained PA/EG/ER composite PCM was poured into a preprepared mold and then dried for 24 h at a constant temperature of 80°C. The determination of the PCM amount is related to the thermal conductivity, specific heat capacity, heat production, and other thermophysical parameters of the battery module. The heat generation rate of the battery module at a 4 C discharge rate must be determined to calculate the quality required by PCM. The mass proportions of paraffin, EG, and ER were 50%, 3%, and 47%, respectively. The melting point and latent heat of the composite PCM were measured with a differential scanning calorimeter (DSC2910, Texas Instruments Inc., USA) in an N 2 atmosphere. The thermal conductivity was measured using an LFA447 NanoFlash™ system (range: 0.1-2000 W/(m K); accuracy: ±5%; and repeatability: ±3%). Figure 2 displays the detailed assembly processing and the specific design schemes of the battery modules. The PCM block was fashioned through a milling machine to acquire a smoothsurface and standard-dimension PCM block consisting of 52 uniformly distributed holes (Figures 2(a) and 2(b)). Thereafter, the commercial cylindrical 18650 batteries with outstanding electrical-chemical properties consistent with the arranged thermocouples were inserted into the PCM heat dissipation skeleton. The electrical connection was then established with a mash welder with 35 W input power and 15 A working current to avoid unsuccessful welding. Two pairs of Ni sheets were utilized to connect the adjacent batteries to avoid damage from high charge-discharge current or high temperature on the power batteries ( Figure 2(c)). After the electrical connection was completed, a multimeter was employed to measure the voltage of every string and confirm  Table 1 illustrates the parameters of the battery module used in this study. Figure 3 displays a schematic of the experimental system. A battery testing system (BTS-50V120A-NTF, Shenzhen Neware Electronics Co., Ltd., China) was used to conduct the charge and discharge tests for the battery modules. T-type thermocouples (OMEGA type TT-T-30-SLE-1M) were installed on the middle surface of each battery to record the temperature data by using a PC-based data acquisition unit (Agilent 34970A) with 1 s interval. All thermocouples were calibrated with a standard thermometer within the temperature range of 0-100°C. The airflow rates of the FC module were controlled and measured by a DC power supply (YK-AD12015, Yucoo Network Equipment Co., Ltd.) and anemograph (GM816A, Shenzhen BENETECH Instrument Co., Ltd.), respectively. The battery modules were placed in the thermostat (BTH-800, Dongguan GBELL Co., Ltd., China) for constant temperature treatment. The temperature range of the thermostat is from −40°C to 150°C, and its temperature departure is ≤±2°C. The heating and cooling rates are 3 and 1°C min -1 , respectively.

Measurement of Battery
Modules. The ITMS and FC cooling system were tested at 25°C. During the experiments, the battery modules obeyed the following charge protocols: 1 C charge rate with a voltage cutoff limit of 54.6 V in galvanostatic mode until the current dropped to the cutoff limit 0.5 A in potentiostatic mode. The allotted equilibration time was at least 1 h after the charging process. Subsequently, the battery modules were discharged at rates of 1, 2, 3, and 4 C. The energy conservation of the battery is expressed as follows: where k is the thermal conductive coefficient; and ρ and C p represent the average density and specific heat of the battery, respectively. Q is the heat generated by unit battery. This equation provides the basic tool for heat conduction analysis, and from its solution, the temperature distribution of the battery can be obtained as a function of time.
The energy conservation of air is given as follows: where∇ v ! presents the velocity vector, C p,a is the specific heat of the air, ρ a is the density, T a is the temperature, and k a is the thermal conductive coefficient of air. The continuity equation for air is written as follows: The momentum conservation equation for air is expressed as follows: where μ represents the dynamic viscosity of air and P represents the static pressure.

Boundary Conditions and the Grid Independence Test.
The boundary conditions of the outlet and inlet were set as the pressure outlet and velocity inlet. No slip boundary was used for the inside surface of the channels. The inlet air temperature was set to 25°C, which was equal to the ambient temperature. Table 2 shows the parameters of the material properties used in the simulation. The grid independence test is essential to ensure the calculation accuracy of the simulation. A grid number of 1204654 is chosen. Figure 4 illustrates the test result. The 5 International Journal of Photoenergy temperature of the middle surface of the single 18650 cylindrical battery in 600 s at a 4 C discharge rate was utilized to calculate the independence of the grid number. A single battery with insulation cotton packaging was tested at a constant discharge rate of 4 C to obtain the thermal behavior and heat generation through the experimental platform ( Figure 3).
The batteries with thermocouples were placed in the middle of the outer surface. The following charge-discharge operations were conducted: a 1 C charge rate with a voltage cutoff limit of 4.2 V in galvanostatic mode was used until the current dropped to the cutoff limit 0.5 A in potentiostatic mode. The battery was discharged at a 4 C discharge rate until the voltage decreased to 2.8 V. Figure 5 presents the comparison between the simulation and experimental results. The relative error did not exceed 2.4%. The simulated results were in agreement with the experimental data, thereby indicating that the battery heat generation in simulations is reasonable [46]. By setting the battery in an adiabatic environment for discharging at different rates, the data of temperature change was collected, and the corresponding calorific value of the battery was calculated according to the formula. In the process of charging and discharging, the lithium-ion battery has four forms of heat generation, mainly       International Journal of Photoenergy including the following: reaction heat Q r , polarization reaction heat Q p , side reaction heat Q s , and Joule heat Q j . Total heat generation can be expressed by the following formula: The average heat generation rate is expressed as follows: where t is the working time of the discharging process, Q is the total heat generation by the battery, and q is the average heat generation rate of the battery. The heat generation rates of the battery were 23.37 and 218.95 W at 1 and 4 C discharge rates, respectively, according to the calculated internal heat source in the adiabatic environment.

Results and Discussion
In this section, the ITMS and FC cooling system for the largescale battery module were discussed and analyzed in detail. First, the large-scale battery modules with the FC cooling system with various parameters, such as the airflow rate and discharge rate, were investigated through experiments and simulations. The numerical simulation and experimental results are in agreement. The comparison between the ITMS and the FC cooling system for decreasing the maximum temperature and balancing the temperature difference is researched.

Heat Dissipation Evaluation during the Charge-Discharge
Cycle Process. In actual driving operations, the running status of power batteries always changes with the practical road condition. The influence of airflow rates on the heat dissipation performance of the as-prepared FC cooling system during cycling was investigated in detail. Figure 6 exhibits that the maximum temperature gradually increased with the increase in the charge-discharge experiment. The thermal energy heat for the first discharging process was stored in the battery module. If the absorbed heat generated by the battery module cannot be dissipated promptly to the environment due to the low heat transfer efficiency, then the battery module will be exposed to a high temperature and exhibit a rising preliminary temperature during the following cycles. Consequently, the initial temperatures of the abovementioned FC cooling system with a 0 m/s airflow rate reached 46.8°C, 59.4°C, 62.1°C, and 62.9°C in the first, second, third, and fourth cycles, respectively, thereby increasing the maximum temperature for each discharging process. Hence, the battery module with a 0 m/s airflow rate could 7 International Journal of Photoenergy not be controlled within the suitable temperature range. Considering the demand for heat dissipation performance, high airflow rates of 1, 2, and 3 m/s were applied in the FC cooling system. The results indicated that the stable initial and maximum temperatures were maintained at approximately 46.6°C, 45.6°C, and 41.4°C accompanied with 1, 2, and 3 m/s, thereby showing enhanced heat transfer efficiency.
At the airflow rate of 3 m/s for the battery module, the temperature at the beginning of the first, second, third, and fourth cycles reached 41.2°C, 41.2°C, 41.3°C, and 41.3°C, respectively, thereby indicating that the temperature range decreased with time and was always within the safe operating temperature range of the power battery. The test results of the FC cooling system with 1 and 2 m/s showed a similar change trend. These results indicate that FC could benefit the thermal management and thermal safety performance of the battery module. The thermal management system with FC can increase the heat dissipation rate of the battery module and reduce the maximum temperature of the battery. With the increase in wind speed, the cooling effect is further optimized. However, for large-scale modules, the wind speed needs to reach 3 m/s in order to present a relatively obvious cooling efficiency. Figure 7(a) illustrates that when the discharge rate was 1 C, the peak temperature of the FC cooling system reached 45.7, 43.9, 41.9, and 39.6°C at 0, 1, 2, and 3 m/s, respectively. The peak temperature (3 m/s) of the FC cooling system decreased by 15.4% compared with that of the FC cooling system (0 m/s) under the same ambient temperature condition. When the airflow rate was 0 m/s, no air disturbance was observed outside the battery module, except the natural convection heat transfer by air. At this time, the convection heat transfer coefficient was small, which cannot eliminate the heat on the battery surface in a timely manner. Consequently, the battery module temperature was high. The air disturbance around the battery became increasingly strong with the airflow increase. At this time, the battery module surface underwent intense convection heat transfer. The convection heat transfer coefficient was large, and the heat removed by the disturbed air increased, resulting in the low temperature of the battery module. This  International Journal of Photoenergy finding indicates the outstanding controlling temperature capacity and heat dissipation effect. The same tendency was observed at 2, 3, and 4 C discharge rates. The same conclusions could be verified through the numerical simulations ( Figure 8). The heat dissipation performance of the FC cooling system with a 3 m/s airflow rate was superior to that of the other three operating conditions.

Temperature Responses of the Battery Modules at Different Discharge Rates.
The conventional discharge current is applied when the power battery module is used under normal working conditions. At this time, the heat generated by the battery module is small. We can see the law through the simulation diagram. However, the battery module will discharge at a high rate when the power battery module for EV works under extreme conditions, such as climbing hills and accelerating. Consequently, the overall temperature of the battery module becomes much higher than that of the conventional discharge rate. Figures 8(a)-8(d) demonstrate that the maximum temperatures at a 1 C discharge rate were 45.5°C, 43.1°C, 42.4°C, and 39.2°C at 0, 1, 2, and 3 m/s, respectively. Figures 8(e)-8(h) show that the maximum temperature at a 4 C discharge rate was 95.5°C, 94.7°C, 91.4°C, and 89.7°C at 0, 1, 2, and 3 m/s, respectively. This finding shows a declining trend with the increased airflow rate. The peak temperature can be measured in the battery's middle position, whereas the temperature of the FC cooling system was higher than that on the central surface of the battery. In comparison with 0 m/s, the 3 m/s airflow rate shows a decrease of 6.3°C and 5.8°C at 1 and 4 C discharge rates, respectively. The cooling effect worsens with the discharge rate increase because of the heat that accumulated in the center of the module. Such heat cannot easily dissipate to ambient surroundings quickly and in a timely manner. The forced air convection cooling method with a high airflow rate had a thermal management effect, but it is always accompanied with high installation and maintenance costs. Hence, the FC cooling system scheme is not the optimum choice at high discharge rates.
When the battery module is discharged at a low rate, the FC cooling system provides a stable cooling effect. However, when the battery heat generation rate increases, the air cooling cannot effectively cool the entire battery module. Severe heat accumulation occurs inside the battery module, which reduces the temperature consistency of the battery module. Figure 9 presents the relationship among shock, bending, and tensile strength with the increase in the PA content for the as-prepared composite PCM. We tested the different components of PCM and found that the overall mechanical properties of PCM showed a declining trend when the PA content increases within a certain range. The PCM mechanical properties were optimal overall when the PA content was 50%, which is conducive to the practical application of PCM. The main reasons for the aforementioned phenomenon are as follows: (1) the participation of the increased amount of low molecular paraffin The increase in paraffin content could cause a decline in the abovementioned mechanical properties among these composite PCM. This phenomenon should be recognized as the crystals of the paraffin acting as the propagation of stress cracking and defect points for the initiation. The curve of shock strength was slight, which was different from the bending and tensile strength, while the paraffin ratio was greater than 50%. This condition could be attributed to the 3D network structure of the composite PCM with 47% mass ratio epoxy, which can surround the paraffin and create a compact crosslinked structure throughout the composite

10
International Journal of Photoenergy PCM. When the PA, ER, and EG contents are 50%, 47%, and 3%, the three mechanical properties exhibited optimum values.
The PCM with good mechanical properties can provide good protection for the battery. The bending and tensile strength of PCM can be improved by adding ER. Therefore, in order to adapt to the practical application process, it is necessary to adjust the correct proportion of PA and ER to make the composite PCM have good latent heat value and mechanical properties.

Thermal Management Analysis of Battery Modules with
Different Heat Dissipation Schemes. The peak temperature and consistency are regarded as fundamental technical parameters considering heat generation and thermal safety. Figure 10 demonstrates that the peak temperatures of the PCM coupled with the FC system decreased to 40.1°C, 47.4°C, 53.4°C, and 63.2°C, thereby showing recession rates of 12.3%, 21.4%, 31.2%, and 31.3% compared with those of the air forced convection thermal management system when the discharge rates were 1, 2, 3, and 4 C, respectively. The curves of the tempera-ture rise of the PCM coupled with the FC cooling system were smoother and gentler as the discharge current increased compared with those of a single air FC scheme, thereby indicating a low temperature ascending gradient and heat generation rate. When the peak temperature reaches 45°C, the temperature rise rate of the PCM+FC module will decrease. This is precisely because the PCM starts to work, which can absorb the heat of the battery and effectively reduce the temperature rise rate of the battery. The temperature distribution will be discussed and analyzed later.
Compared with the single FC cooling system, the PCM +FC cooling system can achieve more efficient cooling efficiency. Even during the high rate discharge process, the maximum temperature of the PCM+FC battery module is about 25°C lower than that of the FC battery module. This reflects that the cooling effect of the PCM+FC cooling system is more suitable for the conditions of the battery module at a high discharge rate. 11 International Journal of Photoenergy temperature uniformity, the maximum temperature difference of the ITMS (3 m/s) and FC cooling system (3 m/s) was studied at four discharge rates ( Figure 11). The temperature uniformity of the two modules deteriorated with time. The maximum temperature difference of the IMTS with 3 m/s was the minimum among these cases, thereby showing that the IMTS balanced the temperature in the battery more efficiently. The battery module containing the composite PCM with a 3 m/s airflow rate exhibited a small, gradual decrease. This phenomenon is attributed to the composite PCM absorbing a large number of heat generated by the batteries utilizing the autologous latent heat during the liquidsolid phase transition process, controlling the peak temperature, and balancing the temperature difference between the batteries in the modules. The thermal conductivity of the composite PCM largely decreased when it reached or exceeded the thermal saturation state, thereby hindering heat dissipation and causing the temperature difference to ascend slowly during the rest of the period. Consequently, the PCM coupling FC cooling system for thermal management exhibited optimum control and stretched the temperature capacity regardless of the constant current discharge condition or uninterrupted charge-discharge cycling operations, which is beneficial for the improvement of thermal safety performance and further increases service life in the practical industrial process.
For poor temperature consistency applicability of the FC cooling system, the PCM+FC cooling system can effectively solve this problem. The PCM can absorb the heat generated by each battery well, and the flow of air can improve the heat dissipation of PCM. It can be seen that the PCM+FC battery module can maintain good temperature consistency even under a high discharge rate.

Conclusions
In this paper, a novel ITMS with PCM coupled with forced air convection cooling has been proposed for large-scale Liion battery modules. The impact factors, including the airflow rate, different discharge rate, and PCM fraction, have been considered to optimize the temperature of the battery module. The main conclusions can be summarized as follows: (1) The change of the airflow rate can mainly affect the cooling effectiveness of the battery module, and the optimum velocity was 3 m/s (2) The numerical simulation has been validated with the experimental results, which were in good agreement for the BTMS with the air cooling system at 1 and 4 C discharge rates (3) The ITMS, including the PCM coupled with the air cooling system for the battery module, demonstrated excellent thermal management effectiveness at 1, 2, 3, and 4 C. The battery module at a 4 C discharge rate can control the maximum temperature below 63.2°C and maintain the maximum temperature difference within 4.8°C Therefore, it should be concluded that the ITMS could promote the development of the thermal management for battery modules in practical application.
Nomenclature C: Capacity (J kg -1 K -1 ) C p : Specific heat (J kg -1 K -1 ) C p,a : Specific heat of air (J kg -1 K -1 ) k a : Thermal conductive coefficient of air (W m -1 K -1 ) m: Mass (kg) p: Pressure (Pa) T: Temperature (°C) t: Time (s) T a : Temperature of air (°C) Q r : Reaction heat Q p : Polarization reaction heat Q s : Side reaction heat Q j : Joule heat Acronyms ITMS: Integrated thermal management system EVs: Electric vehicles BTMS: Battery thermal management system PCM: Phase change material FC: Forced convection Greek Symbols

Data Availability
The data are available upon requesting the corresponding author.

Conflicts of Interest
The authors declare that they have no conflicts of interest.