Proposing a Hybrid BTMS Using a Novel Structure of a Microchannel Cold Plate and PCM

: The battery thermal management system (BTMS) for lithium-ion batteries can provide proper operation conditions by implementing metal cold plates containing channels on both sides of the battery cell, making it a more effective cooling system. The efﬁcient design of channels can improve thermal performance without any excessive energy consumption. In addition, utilizing phase change material (PCM) as a passive cooling system enhances BTMS performance, which led to a hybrid cooling system. In this study, a novel design of a microchannel distribution path where each microchannel branched into two channels 40 mm before the outlet port to increase thermal contact between the battery cell and microchannels is proposed. In addition, a hybrid cooling system integrated with PCM in the critical zone of the battery cell is designed. Numerical investigation was performed under a 5C discharge rate, three environmental conditions, and a speciﬁc range of inlet velocity (0.1 m/s to 1 m/s). Results revealed that a branched microchannel can effectively improve thermal contact between the battery cell and microchannel in a hot area of the battery cell around the outlet port of channels. The designed cooling system reduces the maximum temperature of the battery cell by 2.43 ◦ C, while temperature difference reduces by 5.22 ◦ C compared to the straight microchannel. Furthermore, adding PCM led to more uniform temperature distribution inside battery cell without extra energy consumption.


Introduction
Lithium-ion batteries are crucial in terms of voltage, energy density, self-discharge rate, and cycle life when compared to other energy storage batteries. Their unique properties make them essential for electric vehicles (EVs) and energy storage systems [1]. With the advancement of battery materials and structures, the energy density of lithium-ion batteries has been steadily increasing. As of 2020, the average energy density has reached 300 Wh/kg [2]. Improving the energy density of batteries leads to a higher demand for thermal safety. However, temperature dramatically affects the performance and lifespan of lithium-ion batteries. Low temperatures cause a decrease in battery capacity by slowing down the chemical reaction rate and increasing internal resistance [3,4]. Raised temperatures expedite the deterioration of the battery's structural components, decreasing battery performance and lifespan. This also leads to a drop in thermal safety and can even cause battery thermal runaway [5,6]. Implementing a suitable BTMS is critical to ensure the reliable operation of battery-powered electric vehicles. BTMS keeps thermal safety, battery performance, and longevity while preventing thermal runaway [7]. Consequently, developing a BTMS to maintain the battery temperature within the optimal range of 20 • C-50 • C and improve temperature uniformity is essential [8,9]. Lithium-ion batteries use three types of thermal management: active cooling, passive cooling, and hybrid cooling.
The active cooling method involves the consumption of energy to manage temperature, which is achieved through techniques such as water-cooling or air-cooling thermal management systems [10,11]. Conversely, the passive cooling technique involves temperature management without energy expenditure, typically utilizing materials with high thermal storage capacity, such as PCM [12,13]. Many scholars have extensively researched the air-cooling method and battery layout. For example, Sharma et al. [14,15] studied the effect of air inlet placement on cylindrical battery temperature using forced air cooling. Results showed that bi-directional air intake reduced battery temperature differences and decreased fan energy consumption. Wang et al. [16] enhanced the cooling efficiency of the air-cooled battery BTMS in electric vehicles using parallel plates. The optimal model had two parallel plates with an optimal length and height of 1.5 and 30 mm. As a result of three optimization schemes, T max and ∆T max were reduced by 3.37 K and 5.5 K, respectively. Widyantara et al. [17] introduced an air-cooling system to modify the number of cooling fans and the inlet air temperature. The numerical model of 74 V and 2.31 kWh battery packing was simulated using the lattice Boltzmann method. Results demonstrated that three cooling fans with a 25 • C inlet air temperature achieved the best performance with low power required. The optimized configuration kept all battery cells inside the optimal temperature range, even though the maximum temperature difference was still 15 • C. Liquid cooling systems are more efficient for heat transfer than air cooling systems due to the better thermal conductivity of liquids. They can cool batteries directly with a highly thermally conductive and insulating liquid or indirectly via a cooling plate [18,19]. Wang et al. [20] compared the thermal performance of parallel and serial flow types of BTMS. The parallel flow channel cooling method showed a more uniform temperature distribution, with an optimized layout of parallel channels maintaining the maximum battery temperature below 36 • C while keeping the maximum temperature difference within 4.17 • C at a 3C discharge rate. The U-shaped lightweight liquid cooling method developed by Li et al. [21] improves the thermal safety and weight optimization of prismatic battery cells. It led to a 21% decrease in maximum battery temperature and a 45% reduction in cooling plate weight, showing promise for electric vehicle use. Magnini and Thome [22] researched the effect of flow parameters on heat transfer in slug flow boiling in microchannels. They found that a thin liquid film and a short liquid slug enhance the heat transfer coefficient. Li et al. [23] studied the flow patterns of R1234ze(E) in a microchannel tube. The results showed that due to its lower vapor density, the annular flow starts at a lower quality for R1234ze(E) than R32 and R134a. Oil spray cooling for end winding thermal management was studied by Wang et al. [24] The presence of oil spray was effective in maintaining safe winding temperature under higher heat loads, with elevated flow rates and spray temperatures resulting in more uniform temperature distribution. Jayarajan and Azimov [25] proposed a zig-zag serpentine flow pattern design for cold plates, optimizing several parameters for six different designs. The five-channel design with 18 mm channel width and the seven-channel design with 16 mm width were the best configurations, and the highest temperature at the exit region was 330.84 K for the design with three channels. To enhance the thermal performance of battery modules in EVs and reduce power consumption, Wei et al. [26] optimized the heat transfer of liquid cold plates with serpentine channels. The results showed that by maintaining the pressure drop below 1000 Pa, the maximum temperature differences for discharging rates at 1 C, 2 C, 3 C, and 4 C were controlled within 0.29 K, 1.11 K, 2.17 K, and 3.43 K, respectively. To improve the heat transfer between the battery bottom and the cold plate, Dong et al. [27] designed a wavy channel, which reduced the maximum temperature of the battery module by 1.75 • C compared to the straight-channel design. Lee and Garimella [28] presented an experimental investigation of saturated flow boiling heat transfer in microchannel heat sinks. Their results demonstrated that the pressure drop across the microchannels increases rapidly with heat flux, and the local heat transfer coefficient increases almost linearly with heat flux at low to medium heat fluxes. Li et al. [29] presented the heat transfer coefficient and pressure drop of R32 in a 24-port microchannel tube. The results showed that R32 has a Energies 2023, 16, 6238 3 of 20 higher heat transfer coefficient and lower pressure drop than R134a. Lin et al. [30] proposed a composite silica gel (CSG) coupled with a cross-structure mini-channel cold plate as a cooling system for battery modules. Results showed that the CSG-based liquid system can control the temperature below 45 • C and maintain the temperature difference within 2 • C at a 3C discharge rate. Worwood et al. [31] presented a novel graphite-based fin material for cooling large format pouch-type batteries, with an in-plane thermal conductivity five times greater than aluminum with the same weight. Results showed that the new fin can reduce the peak measured temperature and surface temperature gradient by up to 8 • C and 5 • C, respectively, when compared to aluminum fins under an aggressive electric vehicle duty cycle. In contrast to the active TMS, the passive method does not consume energy. This method's most common cooling system type is heat pipe and PCM. PCM absorbs a significant amount of battery heat with a small volume change during phase change by its high latent heat. High thermal conductivity additives like aluminum wire, foam metal, carbon fiber, and graphite are commonly incorporated into composite Phase Change Materials (CPCM) to enhance thermal conductivity [32]. Huang et al. [33] suggested a flexible CPMC that helps reduce contact resistance in a battery thermal management system. The cooling efficiency of the BTMS was enhanced as the battery temperature decreased by 10 • C at an 18C discharge rate. V.G. Choudhari et al. [34] conducted a numerical analysis on a battery pack with a fin structure and PCM. The findings indicated that the fin structure led to an 8.17% reduction in the maximum battery temperature and increased heat transfer efficiency. Rabiei et al. [35] proposed six microchannel configurations composed of a wavy wall and metal foam embedded microchannel. The study indicated that the wavy wall microchannel performed better heat dissipation at a higher inlet velocity, while at a low inlet velocity, the metal foam embedded microchannel can cool down the battery cell with higher performance and lower power consumption. Ki et al. [36] developed a BTMS including a metal foam as a cooling passage in a 24-pouch cell battery module. The study revealed that with a 4.2 times lower flow rate, the system could attain an optimal temperature and temperature deviation of less than 2 • C. In a study conducted by Weng et al. [37], various fin structures such as V, Y, and X shapes were designed and utilized as heat flow channels in a PCM enclosure. The findings indicated that the X-shaped structure, which possessed an outstanding branched configuration, exhibited the most excellent thermal performance among the tested fin structures.
Based on the literature review, several methods have been explored to improve cooling systems. Using the liquid-cooling system, PCM alone has certain limitations, and combining liquid cooling as an active cooling method and PCM as a passive method has become a primary focus of researchers [38]. This paper presents a new liquid cooling system that utilizes a microchannel cold plate to manage the thermal conditions of a pouch Liion battery by an enhanced coolant distribution design. The primary objectives of this cooling system are to achieve a uniform temperature distribution across the batteries and to efficiently dissipate heat, which are crucial requirements for the effective thermal management of Li-ion batteries. The cooling system consists of inlets/outlets and a thermal conductive plate. The objective of reducing energy consumption and attaining the optimal battery temperature is followed through the proposal of an active and passive cooling system composed of a cooling plate and PCM. The cold plate cooling system comprises an inlet/outlet and a thermally conductive plate. To ensure a uniform temperature distribution in critical areas of the battery, where heat tends to accumulate at the end part, the feasibility of utilizing a Branch Micro-channel Cold Plate (BMCP) instead of a State Micro-channel Cold Plate (SMCP) is examined. BMCP improves heat transfer by increasing the surface area between the battery cell and coolant channels. The branching pattern of the microchannels creates a larger surface area for heat transfer, and the BMCP also enhances the mixing of the coolant. This is particularly effective in the critical zone, where the temperature difference between the cell and coolant is most significant [39]. Furthermore, this study analyzes the effect of operating parameters, including coolant inlet temperature and mass flow rate, on the cooling system. The temperature difference across the battery module and the Energies 2023, 16, 6238 4 of 20 maximum temperature inside the battery pack are also examined to ensure optimal thermal performance and prevent safety concerns. The findings of this study represent progress toward developing a BTMS with enhanced controllability. Table 1 compares the cooling system of the present study and recent studies. The following statements summarize the key concepts of this work:

•
The performance of BTMS was analyzed under active and passive cooling techniques using liquid, PCM, and BMCP.

•
To evaluate the effect of PCM on the heat dissipation performance of BTMS, the minimum value of PCM weight and ideal thermal performance are suggested.

•
To prevent thermal runaway, the effect of BMCP with SMCP has been studied as a factor that contributes to maintaining uniform temperature distribution in the cells.

•
To achieve the best design, the cooling capability has been analyzed under various environmental conditions to maintain the battery temperature within the optimal range.

Model Description
A Li-ion battery cell which used in Volkswagen and BMW i3 commercial EVs is selected for the present study. This battery cell is composed of a carbonate-based electrolyte, Cu cathode, and an Al anode with a capacity of 7 Ah. The battery has 45 mm width, 145 mm height, and 9mm thickness. The geometrical and thermodynamic characteristics of the battery cell are prepared in Table 2. This study investigates and analyzes a module consisting of ten battery cells, eleven cooling plates, and PCM ( Figure 1). To reduce the computational cost and time, a half cell with a half cooling plate and microchannels is simulated and a symmetry boundary condition is applied. Information about the cold plate is provided in Table 3. Also, Table 4 presents the related characteristics of RT-27 PCM utilized in the present research. Specifically, a worst-case scenario involving a discharge rate of 5C is simulated to evaluate the thermal management system. To calculate the transient heat generation within the cells, a User-Defined Function (UDF) code was developed.

Governing Equations
The computational region is partitioned into three sub-regions: the cold plate (aluminum), the fluid part (water), and the battery cell. The governing equations for each of these

Governing Equations
The computational region is partitioned into three sub-regions: the cold plate (aluminum), the fluid part (water), and the battery cell. The governing equations for each of these sub-regions are formulated as follows.

Battery Cell
The generation heat of the battery cell fluctuates and depends on various factors, such as the size of the cell, state of charge (SOC), discharge rate, cell temperature, and electrochemical reaction rate. For this study, the relationship between heat production and electrical parameters, suggested by Bernarda et al. [47], was utilized. By utilizing the thermodynamic energy equation on a correlation derived from a complete cell, the following result is obtained: In the aforementioned equation, the variables T, U, . Q, V, and I represent the cell temperature, open-circuit voltage, cell heat generation rate, terminal voltage, and electric current of cell, respectively.
Equation (1) can be restated as shown in previous studies [48,49]: In the above equation, the variables F, ∆S, . q, I, and R i represent Faraday number (96,485 C/mol), entropy change, rate of internal heat generation per unit volume, discharge current of cell per unit volume, and internal equivalent resistance, respectively. The internal equivalent resistance is influenced by the temperature and state of charge of the cell and can be mathematically represented as described in previous studies [48,49]. ∆S is dependent on the state of charge (SOC) and can be formulated as follows: The state of charge is defined as: where C 0 is the battery capacity, I is the discharge current, and t is the discharge duration (simulation starts at t = 0 and SOC = 1). Additionally, the energy conservation equation for a single cell is given by: where C p,b , k b , and ρ b are specific heat capacity, thermal conductivity, and density, respectively.

Coolant
The fluid section of the system is governed by the equations of mass, momentum, and energy conservation, which can be mathematically expressed as: where p is the pressure.

Cold plate
The energy conservation equation for the solid domain cold plate can be mathematically expressed as follows: where k s , ρ s , and C p,s are thermal conductivity, density, and heat capacity, respectively.

PCM
The enthalpy-porosity method is employed to model the phase change that occurs during the melting process of a solid into a liquid. During the phase change process, the energy equation is expressed as follows: The enthalpy of the PCM is denoted by H and can be calculated using the following formula: The latent heat that exists between the solid and liquid phases is determined by utilizing the energy equation, which can be expressed through the following equations: Energies 2023, 16, 6238 (14) In these equations, S denotes the source term of the PCM, H represents the total enthalpy, and T indicates the temperature. The term H 0 means the sensible enthalpy, while ∆H characterizes the latent heat content that changes between the solid and liquid phases.

Initial and Boundary Condition
At the beginning of the simulation, the ambient temperature was assigned to all domain temperatures. Three different ambient conditions (20 • C, 25 • C, and 35 • C) were considered for this study. A symmetric boundary condition was applied to both sides of the computational domain, as well as to all walls except for the inlet and outlet, which were treated as adiabatic walls. The velocity at the inlet was set to specific magnitudes (0.1, 0.3, 0.5, 0.7, and 1 m/s), the pressure at the microchannel outlet was set equal to the ambient pressure, and the fluid inlet temperature was set to the ambient temperature for each simulation.

Numerical Method and Mesh Independency
The ANSYS workbench is employed to generate a hexahedral structured mesh using appropriate processing tools for computational fluid dynamics (CFD) analysis. The governing equation is solved using a pressure-based solver in ANSYS Fluent software. To couple pressure and velocity, the SIMPLE algorithm is employed. The internal heat generation of the batteries is determined using a user-defined function (UDF) based on mathematical models. To ensure grid independence, three different grid sizes (150,000, 228,000, 380,000, and 454,000) were utilized for the computational domain. It was observed that the maximum temperature and PCM liquid fraction at 100 s, variation between the last two grids during discharge process with 5C discharge rate, 25 • C ambient temperature, and mass flow rate of 0.44 kg/hr was less than 5% and 5.7%, respectively. Figure 2 illustrates the variation of the maximum temperature for the various grids.

Validation of Numerical Model
In this study, validation was conducted for the battery cell section, cooling section, and PCM due to the use of different domains. The present study compares the mean temperature of a 100 Ah battery under natural air-cooling conditions at two different discharge rates, starting from an initial temperature of 20 °C with a convection heat transfer

Validation of Numerical Model
In this study, validation was conducted for the battery cell section, cooling section, and PCM due to the use of different domains. The present study compares the mean temperature of a 100 Ah battery under natural air-cooling conditions at two different discharge rates, starting from an initial temperature of 20 • C with a convection heat transfer coefficient of 5 Wm −2 K −1 . The experimental results of Lin et al. [50] are also included in the comparison and illustrated in Figure 3a. Also, a comparison between the results of a microchannel with straight walls and the findings of Fan et al. [51] at varying flow rates during a 600-s discharge process at an ambient temperature of 27 • C was conducted. The comparison result is presented in Table 5. The maximum difference of 2.92% was achieved for the standard temperature deviation (SDT). SDT is calculated as:

Results and Discussion
This section presents a numerical investigation of the two proposed cold plates for a specific pouch LIB at a high discharge current of 5C for the water-cooling method by various inlet velocities of 0.1 m/s to 1 m/s and three ambient conditions of 20 °C, 25 °C, and 35 °C. Figures 4 and 5 show the variation of maximum temperature and temperature difference of a battery cell for straight microchannels and branched type during a 5C discharge rate, coolant temperature of 25 °C, inlet velocity of 0.5 m/s, and ambient temperature of 25 °C. Both considered parameters experienced a lower value in the branched microchannel cold plate. In the straight microchannel cold plate, the battery maximum temperature rose sharply within 94 s after the discharge process started, to 30.6 °C, and between 95 s and 165 s, the maximum temperature slightly decreased and after that increased to reach 36.61 °C. Battery cell temperature behavior in the branched microchannel  In Figure 3b, the fraction of PCM melting is depicted over 18 min, and a comparison is made between the numerical results of this study and the experimental and numerical findings reported by Shmueli et al. [52]. Initially, the PCM starts melting at 22 • C when the adjacent walls' temperature is above the fluid's melting point at 30 • C. The numerical results obtained from this study have demonstrated acceptable accuracy with a deviation of less than 5% compared to the experimental data, as shown in Figure 3b.

Results and Discussion
This section presents a numerical investigation of the two proposed cold plates for a specific pouch LIB at a high discharge current of 5C for the water-cooling method by  Figures 4 and 5 show the variation of maximum temperature and temperature difference of a battery cell for straight microchannels and branched type during a 5C discharge rate, coolant temperature of 25 • C, inlet velocity of 0.5 m/s, and ambient temperature of 25 • C. Both considered parameters experienced a lower value in the branched microchannel cold plate. In the straight microchannel cold plate, the battery maximum temperature rose sharply within 94 s after the discharge process started, to 30.6 • C, and between 95 s and 165 s, the maximum temperature slightly decreased and after that increased to reach 36.61 • C. Battery cell temperature behavior in the branched microchannel cold plate was different compared to the straight microchannel cold plate, owing to enhanced water flow distribution in critical zones. As shown in Figure 4, the battery maximum temperature starts to increase as the discharge process started, and reached 28.4 • C in 95 s, which is 2.2 • C lower compared to the straight microchannel cold plate. Following, the battery temperature after experiencing a constant temperature between 95 s and 165 s start to increase to attain at 34 • C at the end of discharge. These results show that the maximum temperature of the battery cell reduced by 2.61 • C, which led to the battery cell experiencing better working conditions in BMCP compared with SMCP. In addition, temperature difference variation as another important parameter illustrated in Figure 5. As mentioned before, the maximum temperature difference in the battery cell should remain below 5 • C during battery discharging. As shown in Figure 5, the temperature difference in the SMCP cooling system stayed lower than 5 • C until 267 s after the discharge process was initiated and after, a higher value achieved. However, in the BMCP cooling system, the temperature difference does not reach value higher than 5 • C. The temperature difference at 700 s is 9 • C and 4 • C for SMCP and BMCP, respectively. was initiated and after, a higher value achieved. However, in the BMCP cooling system, the temperature difference does not reach value higher than 5 °C. The temperature difference at 700 s is 9 °C and 4 °C for SMCP and BMCP, respectively.    Figure 6 depicts the temperature contour of the battery cell for SMCP and BMCP at a 5C discharge rate, 0.5 m/s inlet velocity of coolant, and 25 • C ambient temperature. It can be observed that the temperature contour and temperature distribution are the same with the zone around the inlet port being cold and the outlet being hot. During discharge, as coolant liquid flows inside microchannels, the generated heat inside battery cells is conducted to the coolant liquid and the water bulk temperature increases. As expected, the maximum temperature of water was achieved at the outlet port surface. In similar working conditions, SMCP utilizes four straight microchannels to dissipate generated heat along the length of the battery cell, and, as shown in Figure 6a, owing to the liquid coolant, the temperature increased and as a result, the convection heat transfer coefficient decreased and maximum temperature occurred in the vicinity of the outlets port. To overcome this issue, four straight microchannels branched to increase liquid channels and battery cell surface thermal contact. As illustrated in Figure 6b, the maximum temperature in the battery decreased significantly; meanwhile, a better temperature distribution was achieved compared to SMCP.

Performance Characteristics of Enhanced BTMS
One of the important LIB challenges faced by researchers to design a sufficient BTMS is low thermal conductivity causing ununiform temperature distribution inside the battery cell. Figure 7 shows the temperature difference variation for both designed cooling systems at a 5C discharge rate, 25 • C ambient temperature, and specific range of inlet velocity: 0.1 to 1 m/s. It can be seen from Figure 7 that SMCP is not capable of keeping the temperature difference below 5 • C in all inlet velocities during discharge. In 1 m/s, as the best-case, the temperature difference is 9 • C. Obviously, increasing inlet velocity in SMCP has a negligible effect on temperature difference owing to the low thermal conductivity of the battery cell. BMCP, by increasing the thermal contact surface between the battery cell and coolant channels, increases the heat transfer coefficient in the critical zone in which the cell temperature has a significant difference from the inlet port surrounding area. As shown in Figure 7, the temperature difference reduced from 7 • C to 3 • C while the inlet velocity increased to 1 m/s. In addition, proper coolant flow distribution led to higher thermal performance, and increasing inlet velocity from 0.1 m/s to 1 m/s caused a 57% lower temperature difference, while SMCP is able to decrease temperature difference by only 20%. maximum temperature of water was achieved at the outlet port surface. In similar working conditions, SMCP utilizes four straight microchannels to dissipate generated heat along the length of the battery cell, and, as shown in Figure 6a, owing to the liquid coolant, the temperature increased and as a result, the convection heat transfer coefficient decreased and maximum temperature occurred in the vicinity of the outlets port. To overcome this issue, four straight microchannels branched to increase liquid channels and battery cell surface thermal contact. As illustrated in Figure 6b, the maximum temperature in the battery decreased significantly; meanwhile, a better temperature distribution was achieved compared to SMCP. One of the important LIB challenges faced by researchers to design a sufficient BTMS is low thermal conductivity causing ununiform temperature distribution inside the battery cell. Figure 7 shows the temperature difference variation for both designed cooling systems at a 5C discharge rate, 25 °C ambient temperature, and specific range of inlet velocity: 0.1 to 1 m/s. It can be seen from Figure 7 that SMCP is not capable of keeping the temperature difference below 5 °C in all inlet velocities during discharge. In 1 m/s, as the best-case, the temperature difference is 9 °C. Obviously, increasing inlet velocity in SMCP

Effect of Operating Conditions
EVs must be capable of operating in various environmental conditions. In this section, the thermal performance of the designed cooling systems at different ambient temperatures of 20 • C, 25 • C, and 35 • C under a 5C discharge rate and specific range of inlet velocities of 0.1 m/s to 1 m/s is investigated. Figure 8 illustrates the battery maximum temperature variation with respect to inlet velocity in the three mentioned ambient conditions for both proposed cooling systems. Water as the coolant liquid enters the microchannel without precooling and the thermal equilibrium is ambient. Obviously, the maximum temperature decreases as inlet velocity increases in all ambient conditions due to the higher convection heat transfer coefficient. At 0.1 m/s inlet velocity, the maximum temperature difference between BMCP and SMCP is slightly higher than 3 • C, while the inlet velocity reaches a 1 m/s temperature, the difference between the two cooling systems reduced to 2 • C because the heat generation rate in the battery cell remains constant while the heat transfer capacity in both systems increased.
cell and coolant channels, increases the heat transfer coefficient in the critical zone in which the cell temperature has a significant difference from the inlet port surrounding area. As shown in Figure 7, the temperature difference reduced from 7 °C to 3 °C while the inlet velocity increased to 1 m/s. In addition, proper coolant flow distribution led to higher thermal performance, and increasing inlet velocity from 0.1 m/s to 1 m/s caused a 57% lower temperature difference, while SMCP is able to decrease temperature difference by only 20%.

Effect of Operating Conditions
EVs must be capable of operating in various environmental conditions. In this section, the thermal performance of the designed cooling systems at different ambient temperatures of 20 °C, 25 °C, and 35 °C under a 5C discharge rate and specific range of inlet velocities of 0.1 m/s to 1 m/s is investigated. Figure 8 illustrates the battery maximum temperature variation with respect to inlet velocity in the three mentioned ambient conditions for both proposed cooling systems. Water as the coolant liquid enters the microchannel without precooling and the thermal equilibrium is ambient. Obviously, the maximum temperature decreases as inlet velocity increases in all ambient conditions due to the higher convection heat transfer coefficient. At 0.1 m/s inlet velocity, the maximum temperature difference between BMCP and SMCP is slightly higher than 3 °C, while the inlet velocity reaches a 1 m/s temperature, the difference between the two cooling systems reduced to 2 °C because the heat generation rate in the battery cell remains constant while the heat transfer capacity in both systems increased. In view of the temperature difference parameters, two critical zones are the inlet port and outlet port, which are the cold and hot areas of battery cells, respectively. Figure 9 shows the temperature difference variation inside the battery cell with respect to inlet velocity at different ambient temperatures of 20 °C, 25 °C, and 35 °C. LIB thermal conductivity in the thickness direction is significantly lower than along the surface. As a result, in a high water temperature, the lower temperature gradient in the thickness direction In view of the temperature difference parameters, two critical zones are the inlet port and outlet port, which are the cold and hot areas of battery cells, respectively. Figure 9 Energies 2023, 16, 6238 14 of 20 shows the temperature difference variation inside the battery cell with respect to inlet velocity at different ambient temperatures of 20 • C, 25 • C, and 35 • C. LIB thermal conductivity in the thickness direction is significantly lower than along the surface. As a result, in a high water temperature, the lower temperature gradient in the thickness direction caused better temperature distribution and the maximum temperature difference decreased. As shown in Figure 9, the maximum temperature difference for both SMCP and BMCP increased when the ambient temperature was reduced from 35 • C to 20 • C. In addition, the inlet velocity effect on temperature differences was investigated. The maximum temperature difference in 0.1 m/s in 20 • C ambient conditions for SMCP and BMSP is 12.56 • C and 10 • C, respectively. By increasing the inlet velocity to 1 m/s, this parameter reduced to 9.79 • C and 3.22 • C, respectively. As is clear, BMSP reduces the temperature difference significantly compared to SMCP in the same working conditions. Meanwhile, increasing velocity led to a reduced temperature difference of 22% and 67.8% in SMCP and BMSP, respectively. . Battery maximum temperature difference with respect to velocity at a 5C discharge rate for SMCP and BMCP cooling systems.
Figures 10-12 depicts battery cell temperature contour at different ambient temperatures of 20 °C, 25 °C, and 35 °C under a 5C discharge rate respectively. As can be observed, BMCP has a significant effect on enhancing uniform temperature distribution inside the battery cell. Outlet ports as a critical zone in the battery has the highest temperature in the SMCP cooling system, and, as is shown, BMCP, by enhanced liquid coolant distribution design, reduces the temperature deference in all ambient conditions. In addition, the inlet velocity variation effect is revealed in this figure. By increasing the inlet velocity from 0.1 m/s to 1 m/s in both cooling systems, the temperature distribution become more uniform. Vin = 0.1 m/s Vin = 0.5 m/s Vin = 1 m/s Figure 9. Battery maximum temperature difference with respect to velocity at a 5C discharge rate for SMCP and BMCP cooling systems.
Figures 10-12 depicts battery cell temperature contour at different ambient temperatures of 20 • C, 25 • C, and 35 • C under a 5C discharge rate respectively. As can be observed, BMCP has a significant effect on enhancing uniform temperature distribution inside the battery cell. Outlet ports as a critical zone in the battery has the highest temperature in the SMCP cooling system, and, as is shown, BMCP, by enhanced liquid coolant distribution design, reduces the temperature deference in all ambient conditions. In addition, the inlet velocity variation effect is revealed in this figure. By increasing the inlet velocity from 0.1 m/s to 1 m/s in both cooling systems, the temperature distribution become more uniform.
tures of 20 °C, 25 °C, and 35 °C under a 5C discharge rate respectively. As can be observed, BMCP has a significant effect on enhancing uniform temperature distribution inside the battery cell. Outlet ports as a critical zone in the battery has the highest temperature in the SMCP cooling system, and, as is shown, BMCP, by enhanced liquid coolant distribution design, reduces the temperature deference in all ambient conditions. In addition, the inlet velocity variation effect is revealed in this figure. By increasing the inlet velocity from 0.1 m/s to 1 m/s in both cooling systems, the temperature distribution become more uniform.

PCM's Role in Thermal Performance Enhancement
PCM is able to dissipate a large amount of heat through latent heat without consuming any energy. In addition, PCM can be applied in regions that liquid cooling or other cooling methods cannot deploy. Figure 13 shows the maximum temperature variation with respect to inlet velocity for the hybrid cooling system and BMCP for a 5C discharge rate at ambient temperatures of 20 °C and 25 °C. As can be seen, the maximum temperature in 0.1 m/s and 1 m/s decreases 1.5 °C and 0.93 °C , respectively, at 20 °C ambient conditions. Adding PCM to the cold plate has an insignificant effect at reducing the maximum

PCM's Role in Thermal Performance Enhancement
PCM is able to dissipate a large amount of heat through latent heat without consuming any energy. In addition, PCM can be applied in regions that liquid cooling or other cooling methods cannot deploy. Figure 13 shows the maximum temperature variation with respect to inlet velocity for the hybrid cooling system and BMCP for a 5C discharge rate at ambient temperatures of 20 °C and 25 °C. As can be seen, the maximum temperature in 0.1 m/s and 1 m/s decreases 1.5 °C and 0.93 °C , respectively, at 20 °C ambient con-

PCM's Role in Thermal Performance Enhancement
PCM is able to dissipate a large amount of heat through latent heat without consuming any energy. In addition, PCM can be applied in regions that liquid cooling or other cooling  Figure 13 shows the maximum temperature variation with respect to inlet velocity for the hybrid cooling system and BMCP for a 5C discharge rate at ambient temperatures of 20 • C and 25 • C. As can be seen, the maximum temperature in 0.1 m/s and 1 m/s decreases 1.5 • C and 0.93 • C, respectively, at 20 • C ambient conditions. Adding PCM to the cold plate has an insignificant effect at reducing the maximum temperature. Meanwhile, as PCMs are deployed in regions where temperature distribution is ununiform, the temperature difference experienced better conditions compared to BMCP. Figure 14 illustrates the temperature difference inside the battery for BMCP and hybrid cooling with respect to inlet velocity. Temperature difference at 20 • C and 0.1 m/s inlet velocity is 10.5 • C and 6 • C for BMCP and hybrid cooling, respectively. Increasing inlet velocity led to higher heat transfer capacity by the liquid microchannel and the PCM plays little role in heat dissipation and a portion of PCM remains unmelted. Figure 15 depicts the liquid fraction of PCM for various inlet velocities (0.1 m/s, 0.5 m/s, and 1 m/s) and an ambient condition of 20 • C. At 0.1 m/s, PCM fully melted in 240 s and after that, just the liquid cooling system performed as the heat dissipation system. Nevertheless, at 0.5 m/s and 1 m/s, 43% and 97% of PCM remained unmelted. These results revealed that the optimum mass flow rate of cooling liquid can improve hybrid cooling system performance to reach the highest value.

Conclusions
In this study, three cooling systems were proposed to dissipate the generated heat from a 7 Ah LIB pouch cell. The first system consists of a cold plate and straight micro-

Conclusions
In this study, three cooling systems were proposed to dissipate the generated heat from a 7 Ah LIB pouch cell. The first system consists of a cold plate and straight microchannel (SMCP), the second one utilizes a novel branched microchannels cold plate (BMCP) to achieve the best flow distribution strategy. The last system adds PCM to the cold plate as a passive cooling system. Numerical simulations were performed under three ambient conditions: 20 • C, 25 • C, and 35 • C, for a wide range of inlet velocities (0.1 m/s, 0.3 m/s, 0.5 m/s, 0.7 m/s, and 1 m/s), and a discharge rate of 5C.
Results showed that SMCP was not able to maintain the maximum temperature of the battery cell below 50 • C at low inlet velocity. So, a higher rate of mass flow rate is required, which consumes large amount of energy. In addition, in most situations, SMCP cannot provide a proper temperature distribution inside the battery cell.
Increasing the number of microchannels or designing various flow paths like serpentine, zig-zag and etc., was deployed in many studies to decrease the battery maximum temperature. These conventional solutions lead to a lower temperature of the battery cell, while it required higher pumping power, and a higher mass of coolant fluid. However, in this study, branched microchannels only focus on increasing the heat transfer rate in critical areas of the battery cell. Results of BMCP revealed that the proposed path of coolant flow rate, which increased thermal contact with the battery cell critical area, reduced the maximum temperature difference from 10.7 • C to 4 • C, which is in acceptable range. In addition, the maximum temperature decreased by 4.12 • C at a 20 • C ambient temperature and 0.5 m/s inlet velocity. By increasing inlet velocity, the maximum temperature decreases as an obvious result. However, owing to the low thermal conductivity of the battery cell, the rate of decrease in temperature for SMCP was lower than for BMCP with an effective water distribution design.
Finally, to achieve the best thermal performance of the designed BTMS, PCM as a passive cooling system added to BMCP. PCM was only added to special zones and the surrounding microchannel outlet port, which are the most critical zones with the highest battery temperature, to provide a better uniform temperature distribution without consuming energy. The hybrid cooling system decreased the maximum temperature of the battery cell by 1.57 • C and the maximum temperature difference was 4 • C compared to BMCP. Increasing inlet velocity reduced PCM's role in heat dissipation and an optimum flow rate of coolant was required to reach the best performance of BTMS. A hybrid cooling system can enhance BTMS thermal performance mainly by deploying passive cooling methods which do not consume energy. However, adding PCM to the cooling system should be restricted with respect to BTMS total weight. The mass of PCM in the current study can be optimized in future works to meet the best thermal performance with optimized weight of PCM.