Battery Thermal Management Systems: Current Status and Design Approach of Cooling Technologies

: In the current context of transition from the powertrains of cars equipped with internal combustion engines to powertrains based on electricity, there is a need to intensify studies and research related to the command-and-control systems of electric vehicles. One of the important systems in the construction of an electric vehicle is the thermal management system of the battery with the role of optimizing the operation of the battery in terms of performance and life. The article aims to critically analyze the studies and research conducted so far related to the type, design and operating principles of battery thermal management systems (BTMSs) used in the construction of various shaped Li-ion batteries, with focus on cooling technologies. The advantages and disadvantages of the individual components, as well as of the proposed BTM solutions, are extensively investigated, with regard also to the adaptability of these systems to the different Li-ion battery shapes. The information thus synthesized provides the necessary and important information and proposes future directions in research to those interested in this topic to be used to increase the efﬁciency of the thermal management systems of the battery (and with it the global efﬁciency of the electric vehicle).


Introduction
Interest in electric-driven vehicles like hybrid electric vehicles, plug-in hybrid electric vehicles and battery electric vehicles has undoubtedly increased in recent years, following the intensification of environmental regulations on greenhouse gas emissions [1][2][3]. The present lithium-ion (Li-ion) batteries, with their large specific energy, high specific power, low self-discharge rate, high voltage, relatively long life and good recyclability, have been widely used in electric vehicles (EVs) and are regarded as the most suitable energy storage device for electric-driven vehicles [4,5]. However, operating and even storage temperature will affect the performance of Li-ion batteries [6]. Studies have shown that temperature is the main factor influencing battery aging with a negative impact on capacity and internal resistance as well [7]. Additionally, an uneven temperature distribution inside a battery pack can lead to electrically unbalanced cells, and an electric imbalance leads to a capacity loss of the entire battery pack, and to the overcharge of the affected cells during charging, resulting in power losses and increased temperatures [8,9].
Most of the temperature effects are related to chemical reactions occurring in the batteries and also materials used in the batteries [10]. The operational temperature range of EV batteries system commonly varies from −40 • C to 60 • C [11,12]. However, it is generally accepted that to obtain the best performances, the preferred working temperature of EV batteries ranges from 15 • C to 35 • C [6,[13][14][15]. Temperatures outside these limits have a negative impact on the thermal performance of batteries, resulting in significant capacity loss and lifetime reduction. For each temperature degree rise, it is estimated that the battery lifetime shortens by two months [16]. At the same time, it is recommended that the maximum temperature difference inside a battery pack does not exceed 5 • C [17].
• C-rate, charge/discharge current; • State of charge: is in strong connection with the electrochemical reactions and the diffusion of Lithium ions; • Temperature: at high temperatures the electrochemical reactions are intensified and the internal resistance decreases; • Electrochemistry: the active materials influence significantly the amount of heat generated; • State of health: with age the internal resistance increases.
Currently, one important research field is the development of improved cathode and anode materials to enhance the specific power and energy, along with a cost reduction [46]. Safety being the primary concern for both battery and EV manufacturers, it is clear that extreme temperatures must not be reached inside the batteries during the operation of electric-driven vehicles.
In recent years, many researchers attempted to predict the thermal behavior of Li-ion batteries. Mathematical models are developed to estimate the temperature in the core of batteries based on ambient and surface temperatures [47]. Drake et al. [48] determined the generated heat of cylindrical 18,650 cells based on the variation of internal and external temperatures and on the heat flux from its outer surface, by means of a heat flux meter. Other studies [49,50] experimentally investigated the heat generation of similar 18,650 type batteries using an isothermal heat conduction calorimetry (IHC) and an accelerating rate calorimetry (ARC). During the IHC method, the batteries are maintained at a constant temperature by means of a metal or liquid sink, while during the ARC method the battery is allowed to raise its temperature and the heat rejected by the battery is measured [46]. Further experimental investigations [51] present the anisotropic nature of thermal conductance inside the cylindrical batteries, caused by the number of electrode-electrolyte interfaces in radial direction and the absence of these in axial direction. A detailed analysis made by Nazari et al. [45] presents the contribution of reversible and irreversible heat generation, of the active components of batteries (anode, cathode, separator, current collectors), for different Li-ion battery capacities, chemistries and C-rate operations.

Battery Thermal Management Systems
The most important characteristics of a BTMS include a small size, low weight, inexpensiveness, ease of installation, rigidity, reliability and easy maintenance [52]. The main classifications of these systems are the following: active or passive, series or parallel, cooling or heating, air or liquid or phase change material (PCM) or heat pipe (HP) or thermoelectric cooler (TEC) [28,39,53,54]. Active systems consume extra energy to power fans or pumps and are usually implemented in air and liquid cooling systems. Passive systems need specific structures on the surface of batteries to remove heat, such as PCMs or HPs. Generally, active systems have a better heat dissipation capacity, at the cost of increased complexity. PCM-based passive systems are generally used with the aim of improved temperature uniformity, while HP-based systems for their very high thermal conductivity and quick response [55].

Air-Based BTMS
The advantages of these systems are the low cost, simplicity, light extra weight [56,57], and despite their low-conductivity and low-efficiency [58], they are still studied by researchers and used in electric-driven vehicles. A first classification of air-based systems can be made considering whether it is natural or forced convection. Despite its simple structure, the natural method creates large thermal gradients in the battery pack and is highly dependent on ambient conditions [59]. Therefore, a forced-air system is necessary for most applications, despite its higher costs and additional power consumption. A second classification of air-based systems is based on the source of the air, which can be external, cabin pre-conditioned or pre-conditioned by an auxiliary heat exchanger [6].
Regardless of the battery shape, one of the main research directions is the air flow inlet and outlet scheme. Peng et al. [60] made a numerical investigation of different battery configurations, different positions of the air inlet and outlet, and the number of inlets and outlets. After demonstrating that for a battery pack of 20 cylindrical cells the best layout is a 4 × 5 scheme, the authors proposed 18 schemes regarding the inlet and outlet positioning, presented in Figure 1a. Schemes 3, 4, 12 and 13 give the lowest maximum temperature, due to the air flowing through almost all the areas inside the battery pack. An optimization study is made considering a BTMS model with multiple outlets, also shown in Figure 1a, which can slow down the battery heating process. The authors concluded that a longer fluid flow path results in a more uniform temperature distribution, both the minimum temperature and the smallest temperature gradient being obtained in scheme 24. Similar air-cooling configurations are numerically investigated by Chen et al. [61], comparing the thermal performance of asymmetrical and symmetrical parallel air-cooled BTMSs for prismatic batteries. By analyzing various inlet and outlet positions, as shown in Figure 1c, the authors found that the symmetrical systems achieve lower maximum temperature and better temperature uniformity. Various inlet positions and with different number of fans were experimentally investigated by Li et al. [62]. A novel BTMS was proposed by using silica cooling plates with a copper mesh, presented in Figure 1b. The measurements show that for a 5 C discharge rate and an inlet air velocity of 5 m/s, the 1.5 mm thick double silica cooling plate design offers the best results. The authors also found that for the analyzed pouch cells, the side cooling is more effective than the frontal one.
Akinlabi et al. [63] conducted an experimental investigation of an air cooling BTMS with symmetrical vents, at different air flow rates and C-rates. The results point towards the hypothesis that for each current rate there is an optimal air flow rate, above which extra flow will offer little improvement at the cost of increasing energy consumption. The authors suggest that in order to obtain an effective air-cooled BTMS for EVs, bidirectional flow paths are necessary. Wang et al. [64] proposed a reciprocating air flow strategy to reduce the temperature gradient on the surface of a laminated pouch battery. Using infrared images at different depth of discharge (DOD) values, the authors concluded that at the beginning of the discharge process the air flow direction should be forward (towards the terminals) and then change its direction only at a certain DOD. For a 3C discharge rate, the reciprocating process is best to start at 50% DOD. It was found that aperiodic cycles offer the optimum temperature distribution throughout a discharge process, the maximum temperature difference decreasing by 49.5% compared to a one-directional flow, respectively, by 25.5% compared to an isochronous reciprocating cycle.
For prismatic batteries, Z-type, U-type and I-type air flow configurations, shown in Figure 2a-d, are commonly investigated. Zhang et al. [65] developed a transient heat transfer mathematical model that can accurately calculate air velocities and battery temperatures for the three parallel air-flow configurations. The optimization strategy proposes the increase of width around the battery with the highest temperature and the decrease around the one with the lowest temperature, which can enhance temperature uniformity. Considering all the analyzed configurations, the U-type model with optimized channel width and divergence/convergence ducts achieves the best thermal performance and the lowest power consumption. A Z-type air-based BTMS's improvement was proposed by Wang et al. [66] by changing the flow distribution using parallel plates, adjusting the airflow velocity in each cooling channel. An example of adding two parallel plates is presented in Figure 1e. The results show that the ideal solution from a thermal point of view is to put a plate in every cooling channel, however this has a significant negative impact on pressure loss. The acceptable solution found by the authors for the analyzed system is the use of two parallel plates.
Energies 2021, 14, x FOR PEER REVIEW 5 of 33 crease around the one with the lowest temperature, which can enhance temperature uniformity. Considering all the analyzed configurations, the U-type model with optimized channel width and divergence/convergence ducts achieves the best thermal performance and the lowest power consumption. A Z-type air-based BTMS's improvement was proposed by Wang et al. [66] by changing the flow distribution using parallel plates, adjusting the airflow velocity in each cooling channel. An example of adding two parallel plates is presented in Figure 1e. The results show that the ideal solution from a thermal point of view is to put a plate in every cooling channel, however this has a significant negative impact on pressure loss. The acceptable solution found by the authors for the analyzed system is the use of two parallel plates.
(a) (b) (c) Figure 1. (a) Schemes for the battery pack with various inlet and outlet number and position (adapted from source [60]); (b) physical layout of a pouch battery using double silica cooling plates with a cooper mesh (adapted from source [62]); (c) schematics of asymmetrical and symmetrical BTMSs (adapted from source [61]).
For cylindrical batteries, investigations regarding aligned, staggered and cross configurations can be made. Fan et al. [67] performed a series of tests demonstrating that besides occupying the least amount of space from these three arrangements, the aligned configuration also offers the best temperature uniformity and cooling effectiveness, followed by the staggered and lastly the cross battery pack. Regarding the air inlet temperature, although lower values reduce the maximum temperature in the pack, it is shown that an increase from 20 °C to 30 °C at 1 C discharge rate reduces temperature non-uniformity by 12% at an inlet velocity of 0.6 m/s, and by 20% at 1 m/s, respectively. Finally, the study reveals that higher air inlet velocity values lead to an exponential increase in power consumption and a decrease of power efficiency, adding that the aligned configuration requires 23% less energy consumption than the cross pack.  [60]); (b) physical layout of a pouch battery using double silica cooling plates with a cooper mesh (adapted from source [62]); (c) schematics of asymmetrical and symmetrical BTMSs (adapted from source [61]).
For cylindrical batteries, investigations regarding aligned, staggered and cross configurations can be made. Fan et al. [67] performed a series of tests demonstrating that besides occupying the least amount of space from these three arrangements, the aligned configuration also offers the best temperature uniformity and cooling effectiveness, followed by the staggered and lastly the cross battery pack. Regarding the air inlet temperature, although lower values reduce the maximum temperature in the pack, it is shown that an increase from 20 • C to 30 • C at 1 C discharge rate reduces temperature non-uniformity by 12% at an inlet velocity of 0.6 m/s, and by 20% at 1 m/s, respectively. Finally, the study reveals that higher air inlet velocity values lead to an exponential increase in power consumption and a decrease of power efficiency, adding that the aligned configuration requires 23% less energy consumption than the cross pack.  Figure 2. Schematics of parallel air-cooled BTMSs (a) Z-type 3D (adapted from source [65]); (b) Z-type 2D (adapted from source [65]); (c) U-type 2D (adapted from source [65]); (d) I-type 2D (adapted from source [65]); (e) Z-type with two parallel plates (adapted from source [66]).
Zhou et al. [56] proposed a novel cooling strategy based on air distribution pipes for cylindrical batteries and investigated its performance by numerical calculations. The structure of the proposed system is presented in Figure 3. Air enters the pipes from the inlet and then is distributed in-between the batteries. The authors observed that orifices with a 1 mm diameter allow maintaining the maximum temperature below the upper limit, at 38.05 °C. An increase to a diameter of 2 mm further reduced the maximum temperature to 32.55 °C, but also resulted in an increase of temperature gradient inside the analyzed module. Moreover, increasing the rows of orifice also has a positive impact on the thermal performance. The authors also studied the effects of inlet pressure and discharge rate on the heat transfer capacity of the system, highlighting the importance of correlation between inlet pressure, discharge rate and cooling effect.

Liquid-Based BTMS
Liquid-cooled systems are the most used in practical applications [39], due to their high efficiency and compactness. Liquid possesses good heat capacity and heat transfer coefficient, representing a suitable solution for the BTM of EVs [68]. Compared to an aircooling system, 2-3 times less energy is needed to maintain the same average temperature of the batteries [69]. The disadvantages are the increased complexity, added weight, extra costs, rigidity of the structure [55,70]. Moreover, a significant additional power is required Figure 2. Schematics of parallel air-cooled BTMSs (a) Z-type 3D (adapted from source [65]); (b) Z-type 2D (adapted from source [65]); (c) U-type 2D (adapted from source [65]); (d) I-type 2D (adapted from source [65]); (e) Z-type with two parallel plates (adapted from source [66]).
Zhou et al. [56] proposed a novel cooling strategy based on air distribution pipes for cylindrical batteries and investigated its performance by numerical calculations. The structure of the proposed system is presented in Figure 3. Air enters the pipes from the inlet and then is distributed in-between the batteries. The authors observed that orifices with a 1 mm diameter allow maintaining the maximum temperature below the upper limit, at 38.05 • C. An increase to a diameter of 2 mm further reduced the maximum temperature to 32.55 • C, but also resulted in an increase of temperature gradient inside the analyzed module. Moreover, increasing the rows of orifice also has a positive impact on the thermal performance. The authors also studied the effects of inlet pressure and discharge rate on the heat transfer capacity of the system, highlighting the importance of correlation between inlet pressure, discharge rate and cooling effect.  Figure 2. Schematics of parallel air-cooled BTMSs (a) Z-type 3D (adapted from source [65]); (b) Z-type 2D (adapted from source [65]); (c) U-type 2D (adapted from source [65]); (d) I-type 2D (adapted from source [65]); (e) Z-type with two parallel plates (adapted from source [66]).
Zhou et al. [56] proposed a novel cooling strategy based on air distribution pipes for cylindrical batteries and investigated its performance by numerical calculations. The structure of the proposed system is presented in Figure 3. Air enters the pipes from the inlet and then is distributed in-between the batteries. The authors observed that orifices with a 1 mm diameter allow maintaining the maximum temperature below the upper limit, at 38.05 °C. An increase to a diameter of 2 mm further reduced the maximum temperature to 32.55 °C, but also resulted in an increase of temperature gradient inside the analyzed module. Moreover, increasing the rows of orifice also has a positive impact on the thermal performance. The authors also studied the effects of inlet pressure and discharge rate on the heat transfer capacity of the system, highlighting the importance of correlation between inlet pressure, discharge rate and cooling effect.

Liquid-Based BTMS
Liquid-cooled systems are the most used in practical applications [39], due to their high efficiency and compactness. Liquid possesses good heat capacity and heat transfer coefficient, representing a suitable solution for the BTM of EVs [68]. Compared to an aircooling system, 2-3 times less energy is needed to maintain the same average temperature of the batteries [69]. The disadvantages are the increased complexity, added weight, extra costs, rigidity of the structure [55,70]. Moreover, a significant additional power is required

Liquid-Based BTMS
Liquid-cooled systems are the most used in practical applications [39], due to their high efficiency and compactness. Liquid possesses good heat capacity and heat transfer coefficient, representing a suitable solution for the BTM of EVs [68]. Compared to an aircooling system, 2-3 times less energy is needed to maintain the same average temperature of the batteries [69]. The disadvantages are the increased complexity, added weight, extra costs, rigidity of the structure [55,70]. Moreover, a significant additional power is required to circulate the fluid at the desired flow rate, increasing the overall energy consumption of the vehicle [71].
Liquid-cooling is easily applicable for prismatic batteries due their regular shape, allowing a compact arrangement of the BTMS, while cylindrical batteries require special structures to integrate the cooling channels, increasing the system's complexity and mass. The most common solution for prismatic and pouch batteries is the implementation of cooling plates. General investigations analyze the temperature and velocity distribution inside the cooling plates [72], and several studies [73][74][75] investigate by means of orthogonal test the influence of different parameters, such as liquid inlet temperature, cooling plate width, mass flow rate or number of channels. Shang et al. [73] obtained an improvement of 12.61% compared to the maximum temperature without the orthogonal optimization. Ye et al. [74] achieved a lowering of maximum temperature by 5.24% and of the pressure drop by 16.88%. The sensitivity analysis of the orthogonal test also showed that the greatest influence on the pressure loss is caused by the center channel distance and the size of the inlet plenum. The numerical calculations conducted by Monika et al. [75] showed that for cold plates placed between the batteries, a parallel flow arrangement with the inlet placed near the tab region represents the best solution. The authors found that a 5-channeled cold plate of 4 mm width represents the best trade-off between heat transfer and pressure drop. Patil et al. [76] conducted a parametric study to evaluate the thermal performance of different U-turn type microchannel cold plates. The results suggest that a surface area coverage ratio of 0.75 and a hydraulic diameter of the channels of 1.54 mm offer the best cooling capacity. The numerical simulations also show that alternating single inlet and outlet channels decrease the maximum temperature and improve temperature uniformity.
Cai et al. [29] investigated the effect of mass maldistribution on the temperature uniformity of cooling plates, especially from the point of view of channel configuration. The authors concluded that serpentine flow channels offer a better mass distribution than the cavity or parallel configurations. Splitter blades as shown in Figure 4a are used to reduce pressure drop and the results show that narrow channels and longer distances of h 2 offer the better thermal performance for cooling plates. A double-layered cooling channel is proposed by Deng et al. [77], where the coolant enters the cold plate and converges to the center, and then through a vertical connection flows along four diagonal straights channels towards the outlets. The conducted multi-objective optimization shows that the proposed design has better performances than the conventional serpentine and parallel straight channels. Amalesh et al. [78] investigated new mini-channel cold plate designs for the lateral cooling of prismatic batteries, operating at a discharge rate of 3 C. The new proposed designs are represented in Figure 4c. Although all seven designs present better cooling capacity than the conventional design and offer an improved temperature uniformity, the circular slot channel and the zig-zag channel design seem to be the most suitable from a thermal point of view. However, a more detailed study also shows that for almost the same cooling performances, the zig-zag profile causes about 3 times larger pressure drops, proving it not to be the ideal solution to be implemented in the analyzed BTMS. Shen et al. [79] propose a conventional cold plate structure beneath the prismatic batteries, but using refrigerant as coolant medium, and investigated the effects of different flow structures of the cold plate, presented in Figure 4b. The numerical calculations show that although the longer flow paths have a better cooling effect, the increase of the refrigerant's pressure drop is far more significant. Regarding the cross-section of the channels, the decrease of channel thickness causes an increase in collision and friction during the movement of bubbles, which enhances the heat transfer capability of the refrigerant. Moreover, transient simulations show that due to the transition phase of the refrigerant at nearly constant temperature, temperature differences between batteries can be greatly reduced by such a cooling strategy.  Zhang et al. [80] investigated the effect of adding a thermally conductive interface between the batteries and the cold plate, as shown in Figure 5a. The authors prepared a hexagonal boron nitride and silicone rubber (h-BN-SR) composite which presents better mechanical properties and thermal conductivity than pure silicone rubber. At a constant discharge rate of 4 C, the use of the prepared interface decreased the maximum battery temperature by 2.9 °C compared to a 50 µm interface gap. Other studies [81,82] performed comparative analyses between cold plate and other type BTMSs. Dubey et al. [81] compare the performance of immersion cooling and cold plate cooling for cylindrical batteries, the two systems being presented in Figure 5b. The numerical simulations show that at high C-rates the immersion cooling leads to significantly lower battery temperatures. Simultaneously, temperature uniformity across the module is better using a cold plate, because of the higher heat transfer rate and lower thermal characteristics of the dielectric fluid. It Zhang et al. [80] investigated the effect of adding a thermally conductive interface between the batteries and the cold plate, as shown in Figure 5a. The authors prepared a hexagonal boron nitride and silicone rubber (h-BN-SR) composite which presents better mechanical properties and thermal conductivity than pure silicone rubber. At a constant discharge rate of 4 C, the use of the prepared interface decreased the maximum battery temperature by 2.9 • C compared to a 50 µm interface gap. Other studies [81,82] performed comparative analyses between cold plate and other type BTMSs. Dubey et al. [81] compare the performance of immersion cooling and cold plate cooling for cylindrical batteries, the two systems being presented in Figure 5b. The numerical simulations show that at high C-rates the immersion cooling leads to significantly lower battery temperatures. Simultaneously, temperature uniformity across the module is better using a cold plate, because of the higher heat transfer rate and lower thermal characteristics of the dielectric fluid. It was observed that for the same liquid flow rates, immersion-based cooling results in 2.5-3 times better thermal conductance and 15-25 times lower pressure drop compared to the cold plate. Bhattacharjee et al. [82] performed a comparative study between an air-cooled, liquid-cooled with cold plate and immersion-based BTMS. Proving that liquidcooling is more effective than air-cooling, the authors propose a cooling system where the pouch batteries are covered by a nonmetallic sheet and immersed in a dielectric fluid. Such a cooling strategy allows the dissipation of heat in a shorter period, therefore improving the thermal performance of liquid based BTMS.
Energies 2021, 14, x FOR PEER REVIEW 9 of 33 was observed that for the same liquid flow rates, immersion-based cooling results in 2.5-3 times better thermal conductance and 15-25 times lower pressure drop compared to the cold plate. Bhattacharjee et al. [82] performed a comparative study between an air-cooled, liquid-cooled with cold plate and immersion-based BTMS. Proving that liquid-cooling is more effective than air-cooling, the authors propose a cooling system where the pouch batteries are covered by a nonmetallic sheet and immersed in a dielectric fluid. Such a cooling strategy allows the dissipation of heat in a shorter period, therefore improving the thermal performance of liquid based BTMS.
Wang et al. [83] investigated a liquid-based cooling system using another coolant medium than the commonly used water, water/glycol mixture or mineral oil. Indicating the low energy storage density of these coolants, the authors suggest the use of novel coolants made of phase change material (PCM) emulsions, which flow inside 4 mm diameter tubes contained in a heat sink around the cylindrical batteries, as presented in Figure  6a. In the study, two phase change material nano-emulsions containing 10 wt% and 20 wt% OP28E were prepared and investigated from the point of view of their dispersion stability, thermophysical properties or thermal reliability. The studied emulsions showed 1.32 times and 3.65 times improved energy storage capacity compared to water. A different coolant medium is proposed also by Liu et al. [84], who present a BTMS based on liquid metal, which can endure extreme conditions with high ambient temperatures or high discharge rates better than water. Given the shape of the analyzed cylindrical cells, a thermal conductive structure is used with copper pipes as a heat transfer channel, as shown in Figure 6b. An important part of such a study is the power consumption analysis, and the authors found that a variable velocity strategy instead of a traditional constant velocity one can significantly improve the energy efficiency of the BTMS.
A thermally conductive structure around cylindrical batteries combined copper pipes is developed also by Lv et al. [70] and illustrated in Figure 6c. The graphene oxidemodified silica gel has an improved thermal conductivity and transfers the heat efficiently to the fins-enhanced cooper cooling tubes. The analyzed water-cooled battery module presents better cooling capabilities than an air-cooled or a liquid-cooled module without the silica gel. During fast-charging operations with 2 C and 3 C rates, the maximum temperature of the module is at 37.7 °C and 42.0 °C, respectively, while the maximum thermal gradient is kept below 4 °C and 5 °C. Tang et al. [85] proposed a novel liquid-cooling BTMS using microchannel in a heat conducting block with variable contact surface around the cylindrical batteries for a better temperature uniformity, as illustrated in Figure 6d. The authors demonstrated that the bigger the angle difference between contact surfaces, the lower the temperature gradient in the module, while the maximum temperature is also reduced. From a temperature difference of 3.95 °C at an angle difference of 2.5°, the Wang et al. [83] investigated a liquid-based cooling system using another coolant medium than the commonly used water, water/glycol mixture or mineral oil. Indicating the low energy storage density of these coolants, the authors suggest the use of novel coolants made of phase change material (PCM) emulsions, which flow inside 4 mm diameter tubes contained in a heat sink around the cylindrical batteries, as presented in Figure 6a. In the study, two phase change material nano-emulsions containing 10 wt% and 20 wt% OP28E were prepared and investigated from the point of view of their dispersion stability, thermophysical properties or thermal reliability. The studied emulsions showed 1.32 times and 3.65 times improved energy storage capacity compared to water. A different coolant medium is proposed also by Liu et al. [84], who present a BTMS based on liquid metal, which can endure extreme conditions with high ambient temperatures or high discharge rates better than water. Given the shape of the analyzed cylindrical cells, a thermal conductive structure is used with copper pipes as a heat transfer channel, as shown in Figure 6b. An important part of such a study is the power consumption analysis, and the authors found that a variable velocity strategy instead of a traditional constant velocity one can significantly improve the energy efficiency of the BTMS.
A thermally conductive structure around cylindrical batteries combined copper pipes is developed also by Lv et al. [70] and illustrated in Figure 6c. The graphene oxide-modified silica gel has an improved thermal conductivity and transfers the heat efficiently to the fins-enhanced cooper cooling tubes. The analyzed water-cooled battery module presents better cooling capabilities than an air-cooled or a liquid-cooled module without the silica gel. During fast-charging operations with 2 C and 3 C rates, the maximum temperature of the module is at 37.7 • C and 42.0 • C, respectively, while the maximum thermal gradient is kept below 4 • C and 5 • C. Tang et al. [85] proposed a novel liquid-cooling BTMS using microchannel in a heat conducting block with variable contact surface around the cylindrical batteries for a better temperature uniformity, as illustrated in Figure 6d. The authors demonstrated that the bigger the angle difference between contact surfaces, the lower the temperature gradient in the module, while the maximum temperature is also reduced. From a temperature difference of 3.95 • C at an angle difference of 2.5 • , the system descended to 2.58 • C at 15 • . Moreover, the authors demonstrated that the concept with variable contact surface is benefic for temperature uniformity both for cooling and preheating performance.
Energies 2021, 14, x FOR PEER REVIEW 10 of 33 system descended to 2.58 °C at 15°. Moreover, the authors demonstrated that the concept with variable contact surface is benefic for temperature uniformity both for cooling and preheating performance. Figure 6. Liquid cooling BTMSs for cylindrical batteries (a) 3D geometry of the phase change material nano-emulsionbased liquid cooling (adapted from source [83]); (b) structure of liquid-cooled battery pack using liquid metal and copper pipes (adapted from source [84]); (c) schematic diagram and of the graphene oxide-modified silica gel structure with finsenhanced cooling tubes (adapted from source [70]); (d) structure of the battery with variable surface of the heat conducting block (adapted from source [85]).
Wavy channels with liquid flow are another solution for cylindrical batteries and are also adopted by automotive companies like Tesla [6]. The mathematical model developed by Zhao et al. [86] for simulating the thermal behavior of cylindrical batteries cooled by wavy channels with liquid flow show that at the correct fluid inlet conditions (0.5 m/s inlet velocity and 25 °C), the channeled liquid cooling system has the capacity to maintain maximum temperature below 35 °C and maximum temperature difference below 1 °C. Given the risk of increased temperature non-uniformities, especially for the batteries that are close to the outlet of channeled liquid cooling systems, the authors propose heat transfer enhancement methods for the downstream batteries or multi-passage short channels to replace the long wavy single channel. Zhou et al. [87] highlighted the possibility to use half-helical duct structures as a heat exchanger around cylindrical batteries, as shown in Figure 7. Some advantages of such a system are its compact structure, the improvement of heat transfer due to the mixing of fluid particles and the large number of possibilities regarding the number and location of inlet and outlet connections. The pitch and number Figure 6. Liquid cooling BTMSs for cylindrical batteries (a) 3D geometry of the phase change material nano-emulsion-based liquid cooling (adapted from source [83]); (b) structure of liquid-cooled battery pack using liquid metal and copper pipes (adapted from source [84]); (c) schematic diagram and of the graphene oxide-modified silica gel structure with fins-enhanced cooling tubes (adapted from source [70]); (d) structure of the battery with variable surface of the heat conducting block (adapted from source [85]).
Wavy channels with liquid flow are another solution for cylindrical batteries and are also adopted by automotive companies like Tesla [6]. The mathematical model developed by Zhao et al. [86] for simulating the thermal behavior of cylindrical batteries cooled by wavy channels with liquid flow show that at the correct fluid inlet conditions (0.5 m/s inlet velocity and 25 • C), the channeled liquid cooling system has the capacity to maintain maximum temperature below 35 • C and maximum temperature difference below 1 • C. Given the risk of increased temperature non-uniformities, especially for the batteries that are close to the outlet of channeled liquid cooling systems, the authors propose heat transfer enhancement methods for the downstream batteries or multi-passage short channels to replace the long wavy single channel. Zhou et al. [87] highlighted the possibility to use half-helical duct structures as a heat exchanger around cylindrical batteries, as shown in Figure 7. Some advantages of such a system are its compact structure, the improvement of heat transfer due to the mixing of fluid particles and the large number of possibilities regarding the number and location of inlet and outlet connections. The pitch and num- ber of helical ducts did not show any significant differences in the cooling performance, however the concept of alteration of flow direction reduced the maximum temperatures to 30.5 • C and 4.6 • C, respectively, while a larger diameter of the helical ducts also proved to be beneficial. of helical ducts did not show any significant differences in the cooling performance, however the concept of alteration of flow direction reduced the maximum temperatures to 30.5 °C and 4.6 °C, respectively, while a larger diameter of the helical ducts also proved to be beneficial. Cao et al. [88] addressed the problem of numerically investigating a full-size battery pack. It is highlighted that in such an analysis the attention is directed towards the power consumption of the entire system, the coolant flow distribution between different modules and even the electrical circuit networking connection. The study proposes a new design for T-junctions, aiming the reduction of temperature non-uniformity between the modules and enhance flow distribution, with regard to the Reynolds number and pump power requirements.

PCM-Based BTMS
Phase change materials represent an interesting alternative to the conventional forced-air and liquid cooling systems, because it could allow the implementation of a passive BTMS. It is applied especially for cylindrical batteries, because it allows the filling of the unused spaces between the batteries. PCMs have high latent heat capacity and are able to maintain the temperature of the batteries in a small range, close to the melting temperature of the PCM [89]. The melting point is in the 5.5-76 °C temperature range and the phase change process is usually not isotherm [9]. In case of continuous charge and discharge cycles at high rates, there is the risk of the PCM melting completely, leading to thermal instability and leakage of the PCM. Moreover, the melting of the PCM automatically increases its volume, which must be considered when designing a PCM-based BTMS. Another disadvantage of PCMs is their low thermal conductivity, resulting in a difficult heat dissipation.
Numerous studies addressed the problem of low thermal conductivity. Azizi et al. [90] constructed a BTMS made of PCM and aluminum wire mesh plates as a thermal conductivity enhancer. The wire mesh plates with high voidage values are preferred over the aluminum foams for a better filling of the pores during the phase change of the PCM. Zhang et al. [91] study the influence of different concentrations of aluminum nitrite (AlN)based CPCMs. The measurements showed that increasing the AlN mass fraction improved thermal conductivity from 1.48 W/mK to 4.331 W/mK at 50 °C, as well as the tensile strength, bending strength and impact strength. Zhang et al. [92] designed and experimentally investigated a BTMS using copper metal foam saturated PCM. For a 5 C discharge rate the maximum temperature was lowered to 47.86 °C compared to 54.115 °C in the case of pure PCM. The temperature non-uniformity increasing rate also showed smaller values. Jiang et al. [93] addressed the problem of low thermal conductivity of paraffin by absorbing it into expanded graphite (EG) with various mass fractions and then by compressing the resulted powder into cylinder-shaped samples. The authors found that an increased EG mass fraction results in an improved thermal conductivity of 13.85 W/mK, and in a decreased PCM leakage, being at 0.38 wt% loss after 10 phase change Cao et al. [88] addressed the problem of numerically investigating a full-size battery pack. It is highlighted that in such an analysis the attention is directed towards the power consumption of the entire system, the coolant flow distribution between different modules and even the electrical circuit networking connection. The study proposes a new design for T-junctions, aiming the reduction of temperature non-uniformity between the modules and enhance flow distribution, with regard to the Reynolds number and pump power requirements.

PCM-Based BTMS
Phase change materials represent an interesting alternative to the conventional forcedair and liquid cooling systems, because it could allow the implementation of a passive BTMS. It is applied especially for cylindrical batteries, because it allows the filling of the unused spaces between the batteries. PCMs have high latent heat capacity and are able to maintain the temperature of the batteries in a small range, close to the melting temperature of the PCM [89]. The melting point is in the 5.5-76 • C temperature range and the phase change process is usually not isotherm [9]. In case of continuous charge and discharge cycles at high rates, there is the risk of the PCM melting completely, leading to thermal instability and leakage of the PCM. Moreover, the melting of the PCM automatically increases its volume, which must be considered when designing a PCM-based BTMS. Another disadvantage of PCMs is their low thermal conductivity, resulting in a difficult heat dissipation.
Numerous studies addressed the problem of low thermal conductivity. Azizi et al. [90] constructed a BTMS made of PCM and aluminum wire mesh plates as a thermal conductivity enhancer. The wire mesh plates with high voidage values are preferred over the aluminum foams for a better filling of the pores during the phase change of the PCM. Zhang et al. [91] study the influence of different concentrations of aluminum nitrite (AlN)-based CPCMs. The measurements showed that increasing the AlN mass fraction improved thermal conductivity from 1.48 W/mK to 4.331 W/mK at 50 • C, as well as the tensile strength, bending strength and impact strength. Zhang et al. [92] designed and experimentally investigated a BTMS using copper metal foam saturated PCM. For a 5 C discharge rate the maximum temperature was lowered to 47.86 • C compared to 54.115 • C in the case of pure PCM. The temperature non-uniformity increasing rate also showed smaller values. Jiang et al. [93] addressed the problem of low thermal conductivity of paraffin by absorbing it into expanded graphite (EG) with various mass fractions and then by compressing the resulted powder into cylinder-shaped samples. The authors found that an increased EG mass fraction results in an improved thermal conductivity of 13.85 W/mK, and in a decreased PCM leakage, being at 0.38 wt% loss after 10 phase change cycles, without impacting significantly the phase change temperature of the pure PCM. Wu et al. [94] suggest that an EG mass fraction of 15%-20% is recommended in order to increase thermal conductivity. An enhanced CPCM is proposed, made of pyrolytic graphite sheets (PGS), which presents better heat dissipation and temperature uniformity compared to PCM. The propagation of thermal runaway is also analyzed, the PCM/PGS module needing a smaller spacing between batteries to prevent propagation, leading to a reduction of 71.4% in volume compared to the PCM module. Talluri et al. [95] concluded that expanded graphite PCM is the best solution for PCM-based BTMS, enhancing heat transfer rate by 16.4%, and observed no leakage using this solution. During high discharge rates of 25 min, the maximum temperature with PCM increased to 28 • C compared to 35.9 • C without PCM.
Besides low thermal conductivity, the problem of PCM leaking has also been intensely investigated. In liquid phase the homogeneity of the PCM module decreases, and even worse, safety risks rise because of the build-up of leaked, combustible PCM. A solution could be microencapsulation, but high costs [96], the required laborious procedures [97] and the change in volume during temperature variations [98] reduce its applicability. Chen et al. [99] proposed a PCM with high thermal conductivity and a solid-solid phase change, enhanced with fins to further increase thermal performance. Lv et al. [100] developed a paraffin-based CPCM with improved anti-leakage and anti-volume-change properties, obtained by adding nanosilica into the CPCM, conform the schematics from Figure 8a. The nanosilica dispersed in the expanded graphite pores adsorb liquid phase paraffin and restrain the expansion of graphite. This prevents leakage and reduces greatly the volume change of the system. The results show that after 20 h at 55 • C the loss of mass of the prepared samples is 1.15% for CPCM-NS0, 0.54% for the CPCM-NS7 and the lowest 0.39% for the CPCM-NS5. When increasing the temperature to 60 • C, these losses become 4.08%, 1.25% and 0.53%, respectively, showing that a too high participation of nanosilica can lead to an increase of leakage, because of the aggregation of nanosilica resulting in a nonhomogenous distribution of the CPCM's components. The volume change experiments show that the volumetric thermal expansion increment in CPCM-NS0 is 4 times higher than that of the CPCM-NS5.5, which also means a better cooling performance of the proposed CPCM due to the reduced gap between the batteries and the CPCM module. Liu et al. [101] also focused on preparing a low leakage CPCM with high thermal conductivity and high latent heat. The proposed paraffin/expanded graphite/high-density polyethylene/nanosilver CPCM is presented in Figure 8b. Huang et al. [102] developed a flexible form-stable CPCM in order to overcome both the issues of paraffin leakage after melting and of high thermal contact resistance of rigid CPCMs. Flexible polymer is used as a supporting material and expanded graphite is added to enhance thermal conductivity. The proposed solution reduces thermal contact resistance to 461 mm 2 K/W compared to 1818 mm 2 K/W for a conventional CPCM, which means a reduction of temperature by 6.5 • C.
Huang et al. [103] also highlight the need of increasing thermal conductivity of pure PCM, but simultaneously express the necessity of finding a balance between latent heat and thermal conductivity. The authors declare that above 5 W/mK the decrease of maximum temperature becomes less significant and suggest an optimal latent heat of 225 J/g and an optimal porosity of 94% for the best thermal management performance.
Yan et al. [104] conducted an experimental study to determine the optimal phase change temperature of a PCM used in BTMS, highlighting the need to test these systems in dynamic cycling, which is closer to real-life conditions. The authors compared three PCMs with phase change temperatures of 36, 45 and 58 • C. The experiments prove that in dynamic cycling conditions, the laying-aside time between steps is of great importance. The authors recommend the use of the PCM with a phase change temperature of 45 • C. Gscwander et al. [105] analyzed the storage capacity of different PCM emulsions in function of their phase transition behavior and supercooling effect. The melting temperature of the different PCMs was between 16 and 40 • C and a 35 wt% formulation was prepared for the emulsions. Compared to water in the same temperature range, storage capacities increased 2-2.7 times. and high latent heat. The proposed paraffin/expanded graphite/high-density polyethylene/nano-silver CPCM is presented in Figure 8b. Huang et al. [102] developed a flexible form-stable CPCM in order to overcome both the issues of paraffin leakage after melting and of high thermal contact resistance of rigid CPCMs. Flexible polymer is used as a supporting material and expanded graphite is added to enhance thermal conductivity. The proposed solution reduces thermal contact resistance to 461 mm 2 K/W compared to 1818 mm 2 K/W for a conventional CPCM, which means a reduction of temperature by 6.5 °C. Huang et al. [103] also highlight the need of increasing thermal conductivity of pure PCM, but simultaneously express the necessity of finding a balance between latent heat and thermal conductivity. The authors declare that above 5 W/mK the decrease of maximum temperature becomes less significant and suggest an optimal latent heat of 225 J/g and an optimal porosity of 94% for the best thermal management performance.
Yan et al. [104] conducted an experimental study to determine the optimal phase change temperature of a PCM used in BTMS, highlighting the need to test these systems in dynamic cycling, which is closer to real-life conditions. The authors compared three PCMs with phase change temperatures of 36, 45 and 58 °C. The experiments prove that in dynamic cycling conditions, the laying-aside time between steps is of great importance. The authors recommend the use of the PCM with a phase change temperature of 45 °C. Gscwander et al. [105] analyzed the storage capacity of different PCM emulsions in function of their phase transition behavior and supercooling effect. The melting temperature of the different PCMs was between 16 and 40 °C and a 35 wt% formulation was prepared for the emulsions. Compared to water in the same temperature range, storage capacities increased 2-2.7 times.
Weng et al. [106] investigated a single cylindrical battery surrounded by a PCM layer, analyzing the effect of three layer thicknesses, illustrated in Figure 9a. Given that a too thin PCM layer cannot store enough thermal energy, while a too large one decreases the energy density of a battery pack, the authors found that the optimal cooling performance is offered by the 10 mm thick PCM layer. Jilte et al. [107] highlighted the disadvantages of the conventional PCM-based BTMS, consisting of a rectangular PCM-filled battery pack, namely the difference of heat rejection rate between the batteries in the middle and the corner of the pack. A different cooling system is proposed, consisting of embedded PCM flow passages, presented in Figure 9b. The numerical simulations show that a 4 mm thick PCM is sufficient for heat dissipation, the maximum cell temperature rise being below 5 °C under 40 °C ambient air temperature supply. Additionally, the temperature uniformity is kept within 0.05 °C at 2 C rate and below 0.12 °C at 4 C discharge rate. Nasehi et al. [108]  Weng et al. [106] investigated a single cylindrical battery surrounded by a PCM layer, analyzing the effect of three layer thicknesses, illustrated in Figure 9a. Given that a too thin PCM layer cannot store enough thermal energy, while a too large one decreases the energy density of a battery pack, the authors found that the optimal cooling performance is offered by the 10 mm thick PCM layer. Jilte et al. [107] highlighted the disadvantages of the conventional PCM-based BTMS, consisting of a rectangular PCM-filled battery pack, namely the difference of heat rejection rate between the batteries in the middle and the corner of the pack. A different cooling system is proposed, consisting of embedded PCM flow passages, presented in Figure 9b. The numerical simulations show that a 4 mm thick PCM is sufficient for heat dissipation, the maximum cell temperature rise being below 5 • C under 40 • C ambient air temperature supply. Additionally, the temperature uniformity is kept within 0.05 • C at 2 C rate and below 0.12 • C at 4 C discharge rate. Nasehi et al. [108] investigated the use of three PCM layers around a prismatic battery, as shown in Figure 9c, each PCM having different properties, but concluded that the best performances are attained when using a single PCM. Another concept presented by Zhao et al. [109] focuses also on the capacity of preventing thermal runaway propagation in a battery module, by means of a PCM composite board illustrated in Figure 9d. In order to balance heat transfer and thermal isolation capacity simultaneously, the authors propose a structure composed of a thermally conductive shell, PCM to absorb the generated heat and an insulation panel in the middle of the structure, to slow down heat transfer between batteries and improve safety. Ling et al. [110] analyzed the thermal performance of PCMs with different thermal conductivities at low-temperature operation. Two composites are prepared, a 60 wt% RT44HC/expanded graphite with high thermal conductivity, and a 60 wt% RT44HC/fumed silica with low thermal conductivity, and tested in a BTMS at −10 °C and 5 °C. Although the warm battery does not cool down as fast, PCMs with low thermal conductivity are not suited for BTMSs, causing significant temperature differences in the pack (>12 °C). The composite with expanded graphite prevented the overheating of the battery and improved temperature uniformity. Zhong et al. [22] proposed a BTMS able to cool but also preheat the battery pack for insuring normal operating temperatures also in cold climates. The study presents CPCM cooling with an additional resistance wire preheating function, as shown in Figure 10. The preheating system can increase battery temperature to 40 °C in 300 s, using slightly more energy than a pure resistance wire system. Ling et al. [110] analyzed the thermal performance of PCMs with different thermal conductivities at low-temperature operation. Two composites are prepared, a 60 wt% RT44HC/expanded graphite with high thermal conductivity, and a 60 wt% RT44HC/fumed silica with low thermal conductivity, and tested in a BTMS at −10 • C and 5 • C. Although the warm battery does not cool down as fast, PCMs with low thermal conductivity are not suited for BTMSs, causing significant temperature differences in the pack (>12 • C). The composite with expanded graphite prevented the overheating of the battery and improved temperature uniformity. Zhong et al. [22] proposed a BTMS able to cool but also preheat the battery pack for insuring normal operating temperatures also in cold climates. The study presents CPCM cooling with an additional resistance wire preheating function, as shown in Figure 10. The preheating system can increase battery temperature to 40 • C in 300 s, using slightly more energy than a pure resistance wire system. Energies 2021, 14, x FOR PEER REVIEW 15 of 33 Figure 10. Battery module design with CPCM cooling and a resistance wire wrapped around the batteries (adapted from source [22]).

Heat Pipe-Based BTMS
Heat pipes are heat transfer devices with extremely high thermal conductivity, 90 times higher than copper [111], ensuring better heat dissipation and temperature uniformity. They represent a passive, compact, lightweight BTMS, with no need of maintenance and long life cycle [112]. The working principle consists in the phase change of the internal working fluid, which absorbs heat at the evaporation section and releases it at the condensation section. As long as there is a temperature gradient between the two sections, the transfer of heat goes on, without any power consumption. Despite its very high thermal conductivity however, in most applications heat pipe is coupled with forced air or liquid cooling for a more effective heat dissipation at the condensation section.
Heat pipe-based BTMSs are studied mostly for prismatic and pouch batteries due to the adaptability of the shapes and the large available contact areas, while in case of cylindrical batteries special attention needs to addressed to the contact surface between heat pipe and batteries.
Zhang et al. [113] proposed a passive heat pipe-based BTMS for prismatic batteries, using flat heat pipes with fins, as presented in Figure 11a. Compared to natural convection and aluminum plates, the maximum temperature was 73.7% lower and the maximum temperature difference 50.1% lower. Behi et al. [112] designed a smart heat pipe-based BTMS, using the least number of heat pipes based on the multizone analysis of temperature distribution determined by thermal imaging. The high discharge investigation revealed that the most critical zone of a prismatic LTO cell is the center top region. The thermal analysis, illustrated in Figure 11c, showed that a single heat pipe with an optimal placement is sufficient. Kleiner et al. [114] introduced a concept where heat pipes are attached to the terminals of prismatic batteries, as presented in Figure 11b. This reduces the heat flow path and thermal resistance from the heat generation sources, located in the busbar and inside the battery, to the cooling system at the bottom of prismatic batteries. The results show that without any increase in charging time, temperature at the terminals and in the jelly roll are 19 °C, respectively 2 °C lower. Temperature distribution is also improved, due to conduction of 27% of the generated heat.
Liu et al. [115] highlighted the necessity of developing more accurate mathematical models that describe the heat transfer characteristics of heat pipes. The authors propose a segmented thermal model in which the evaporator and condenser region have their own individual equivalent conductivity and thermal parameters, which proves to be more accurate than the traditional non-segmented model. Additionally, the results show that for an effective heat pipe based BTMS it is recommended to have assistance by forced convection. Figure 10. Battery module design with CPCM cooling and a resistance wire wrapped around the batteries (adapted from source [22]).

Heat Pipe-Based BTMS
Heat pipes are heat transfer devices with extremely high thermal conductivity, 90 times higher than copper [111], ensuring better heat dissipation and temperature uniformity. They represent a passive, compact, lightweight BTMS, with no need of maintenance and long life cycle [112]. The working principle consists in the phase change of the internal working fluid, which absorbs heat at the evaporation section and releases it at the condensation section. As long as there is a temperature gradient between the two sections, the transfer of heat goes on, without any power consumption. Despite its very high thermal conductivity however, in most applications heat pipe is coupled with forced air or liquid cooling for a more effective heat dissipation at the condensation section.
Heat pipe-based BTMSs are studied mostly for prismatic and pouch batteries due to the adaptability of the shapes and the large available contact areas, while in case of cylindrical batteries special attention needs to addressed to the contact surface between heat pipe and batteries.
Zhang et al. [113] proposed a passive heat pipe-based BTMS for prismatic batteries, using flat heat pipes with fins, as presented in Figure 11a. Compared to natural convection and aluminum plates, the maximum temperature was 73.7% lower and the maximum temperature difference 50.1% lower. Behi et al. [112] designed a smart heat pipe-based BTMS, using the least number of heat pipes based on the multizone analysis of temperature distribution determined by thermal imaging. The high discharge investigation revealed that the most critical zone of a prismatic LTO cell is the center top region. The thermal analysis, illustrated in Figure 11c, showed that a single heat pipe with an optimal placement is sufficient. Kleiner et al. [114] introduced a concept where heat pipes are attached to the terminals of prismatic batteries, as presented in Figure 11b. This reduces the heat flow path and thermal resistance from the heat generation sources, located in the busbar and inside the battery, to the cooling system at the bottom of prismatic batteries. The results show that without any increase in charging time, temperature at the terminals and in the jelly roll are 19 • C, respectively 2 • C lower. Temperature distribution is also improved, due to conduction of 27% of the generated heat.
Liu et al. [115] highlighted the necessity of developing more accurate mathematical models that describe the heat transfer characteristics of heat pipes. The authors propose a segmented thermal model in which the evaporator and condenser region have their own individual equivalent conductivity and thermal parameters, which proves to be more accurate than the traditional non-segmented model. Additionally, the results show that for an effective heat pipe based BTMS it is recommended to have assistance by forced convection. (c) Figure 11. (a) Passive BTMS using flat heat pipes with fins (adapted from source [113]); (b) thermal cooling of prismatic battery using heat pipe (adapted from source [115]); (c) front and back side of prismatic cell with heat pipe in 8 C discharging rate (adapted from source [112]).
Pulsating or oscillating heat pipes (PHP or OHP) have a simple structure to work during a two-phase flow regime, due to the thermally driven back-and-forth oscillation of bubbles and liquid plugs within its serpentine tubes. Moreover, conventional working fluids can be replaced with other mediums that improve heat transfer performance. Chen et al. [116] proposed an improved pulsating heat pipes using water-based TiO2 nanofluid as working medium. An increase in battery heat generation leads to an increase in the temperature gradient between the condensation and evaporation section of the hat pipe, leading to an accelerated start of the PHP. Simultaneously, an ambient temperature increase to 35 °C is overcame by the proposed system, resulting in a maximum temperature of 42.22 °C on the battery surface and a maximum temperature difference below 2 °C. Wei et al. [117] proposed a plug-in oscillating heat pipe consisting of a flat-plate evaporator and a tube condenser, using as working fluid binary-fluid mixtures of ethanol aqueous solution. The measurements show that the volumetric filling ration and water-ethanol mixing ratio have a significant influence on the thermal performance of oscillating heat pipes, affecting startup evaporator temperature, thermal resistance or air-to-boil temperature. The best heat transfer characteristics were found to be for a 1:1 water-ethanol mixing ratio and a volumetric filling ratio of 30%.

Hybrid BTMS
Hybrid BTMSs usually combine active and passive methods, with the aim of enhancing the heat transfer process of passive systems and offering the possibility to actively control their working process. In general, PCMs are used for a better temperature uniformity and the active system reduces the risk of heat saturation, while HPs are used for their high thermal conductivity which significantly improves local heat transfer and the active system enhances the heat release at the condensation section [55]. Although forced air can be successfully implemented as the active system for lower thermal load systems, liquid cooling offers higher efficiency. Another classification criterion is the combination Figure 11. (a) Passive BTMS using flat heat pipes with fins (adapted from source [113]); (b) thermal cooling of prismatic battery using heat pipe (adapted from source [115]); (c) front and back side of prismatic cell with heat pipe in 8 C discharging rate (adapted from source [112]).
Pulsating or oscillating heat pipes (PHP or OHP) have a simple structure to work during a two-phase flow regime, due to the thermally driven back-and-forth oscillation of bubbles and liquid plugs within its serpentine tubes. Moreover, conventional working fluids can be replaced with other mediums that improve heat transfer performance. Chen et al. [116] proposed an improved pulsating heat pipes using water-based TiO 2 nanofluid as working medium. An increase in battery heat generation leads to an increase in the temperature gradient between the condensation and evaporation section of the hat pipe, leading to an accelerated start of the PHP. Simultaneously, an ambient temperature increase to 35 • C is overcame by the proposed system, resulting in a maximum temperature of 42.22 • C on the battery surface and a maximum temperature difference below 2 • C. Wei et al. [117] proposed a plug-in oscillating heat pipe consisting of a flat-plate evaporator and a tube condenser, using as working fluid binary-fluid mixtures of ethanol aqueous solution. The measurements show that the volumetric filling ration and water-ethanol mixing ratio have a significant influence on the thermal performance of oscillating heat pipes, affecting startup evaporator temperature, thermal resistance or air-to-boil temperature. The best heat transfer characteristics were found to be for a 1:1 water-ethanol mixing ratio and a volumetric filling ratio of 30%.

Hybrid BTMS
Hybrid BTMSs usually combine active and passive methods, with the aim of enhancing the heat transfer process of passive systems and offering the possibility to actively control their working process. In general, PCMs are used for a better temperature uniformity and the active system reduces the risk of heat saturation, while HPs are used for their high thermal conductivity which significantly improves local heat transfer and the active system enhances the heat release at the condensation section [55]. Although forced air can be successfully implemented as the active system for lower thermal load systems, liquid cooling offers higher efficiency. Another classification criterion is the combination method, which can be parallel, series or both. The main disadvantages of hybrid BTMSs regularly are their weight, volume, complexity and energy consumption.
Ling et al. [118] showed that using pure PCM as a passive system, at discharge rates of 1.5-2 C the battery temperature of 60 • C was reached at the third cycle, while using a hybrid system coupled with forced air cooling the maximum temperature was maintained under 46 • C. Another study, performed by Zhao et al. [16] compared various concepts of heat pipe-based BTMSs: natural convection, fan cooling in horizontal and vertical position, and wet cooling by immersion and by using water sprays, respectively. The study does not recommend the cooling of the condenser section in a bath, because of the formation of bubbles that drastically thermal conductivity. The results also show that the inclination angle of ultra-thin heat pipes has no effect on their performances. The water spray strategy proved to be the most efficient, keeping the maximum temperature under 30 • C with only four sprays for 3000 s.
A hybrid BTMS combining PCM and water-cooling is proposed by Molaeimanesh et al. [119], for prismatic cells in different configurations, presented in Figure 12a-d. The numerical simulations show that the PCM stay in solid phase for a longer period and that the surface temperature of the battery is better controlled with the parallel/series configuration. However, this solution offers poor thermal uniformity since both the PCM and the cooling water are in direct contact with the battery. The most reduced temperature differences are recorded for the series configuration. Another hybrid BTMS for prismatic batteries, including a high thermal conductivity and low leakage CPCM and water-cooling is presented by Liu et al. [101]. By coupling the CPCM with a single water channel during a 3 C discharge -1 C charge test, the maximum temperature was reduced from 92 to 59.5 • C and the maximum temperature difference from 16.95 to 6.1 • C, demonstrating a better cooling effect than the CPCM-only system. Further, the authors proposed several configurations by changing the hose connections, presented in Figure 12e. The lowest temperature and temperature difference are obtained using the scheme 6, with four cold-water inlets, which also means a more reduced pressure loss. The results show that in case discharge-charge cycles it is more effective the change the water inlet temperature than the flow rate. From an energy consumption point of view, when the battery pack is moderately heated, varying only the pump power for regulating the flow rate is sufficient, while during severe heat generation, both the pump and the compressor should be adjusted simultaneously. Bai et al. [120] proposed the use of a water-cooling plate placed near the electrodes of pouch batteries, where most of the heat is generated, while the PCM is used to dissipate the generated heat of other parts of the battery and ensure temperature uniformity, as illustrated in Figure 12f. The results showed that increasing mass flow rate was more efficient than other methods in improving the thermal performance, at the cost of an increased energy consumption. The authors state that such a system does not need the increase of thermal conductivity of PCM by additives, and rather suggest the choice of a PCM with a higher melting point or the decrease of water inlet temperature.
For cylindrical cells, Song et al. [121] developed a numerical model using PCM between the batteries and a cold plate placed under, and as supplementary heat transfer elements, a heat spreading plate and thermal columns were added. The results showed that compared to the PCM-only system, the rate of temperature increase is significantly reduced after the start of PCM melting. Zhao et al. [122] proposed a hybrid BTMS for cylindrical batteries based on CPCM coupled with liquid cooling tubes, presented in Figure 12g. The addition of copper foam to the pure PCM reduces battery surface temperatures by 14 • C and offers great temperature uniformity.  Figure 12. Hybrid BTMSs using PCM coupled with liquid cooling (a) parallel configuration (adapted from source [119]); (b) series 1 configuration (adapted from source [119]); (c) series 2 configuration (adapted from source [119]); (d) parallel/series configuration (adapted from source [119]); (e) different hose connection schemes (adapted from source [101]); (f) cooling plate placed near the electrodes (adapted from source [120]); (g) CPCM with copper foam (adapted from source [122]).
Feng et al. [123] proposed a cooling device for the thermal and strain management of cylindrical batteries using a design that combines heat pipes and fins, presented in Figure  13a. The optimal BTMS for the presented concept is the one with forced convection, using a fan in the center of the battery pack. At a 1 C discharge rate, the maximum temperature is 15 °C lower than with natural convection and the rate of temperature increase is lowered. Ye et al. [124] used micro heat pipes that have grooves and fins on the inner walls, Figure 12. Hybrid BTMSs using PCM coupled with liquid cooling (a) parallel configuration (adapted from source [119]); (b) series 1 configuration (adapted from source [119]); (c) series 2 configuration (adapted from source [119]); (d) parallel/series configuration (adapted from source [119]); (e) different hose connection schemes (adapted from source [101]); (f) cooling plate placed near the electrodes (adapted from source [120]); (g) CPCM with copper foam (adapted from source [122]).
Feng et al. [123] proposed a cooling device for the thermal and strain management of cylindrical batteries using a design that combines heat pipes and fins, presented in Figure 13a. The optimal BTMS for the presented concept is the one with forced convection, using a fan in the center of the battery pack. At a 1 C discharge rate, the maximum temperature is 15 • C lower than with natural convection and the rate of temperature increase is lowered. Ye et al. [124] used micro heat pipes that have grooves and fins on the inner walls, increasing the heat transfer capacity. Experimental measurements were made adding fins at the end of the micro heat pipes and using forced air-convection. The study of Liu et al. [115] analyzes besides the enhanced thermal model also the cooling performance of ultra-thin micro heat pipes combined with forced air-convection. The results show that for an effective heat pipe based BTMS it is recommended to have assistance by forced convection. The coupling of the two technologies can reduce maximum temperature by 7.1 • C at a 2 C discharge rate, compared to no heat pipes. Another system combining heat pipe with forced air cooling was developed by Tran et al. [125], where thermally conductive resin fills the spacing between the cylindrical batteries and the heat pipe with heat sink are placed on the lateral shell of the module, as shown in Figure 13b. Compared to using only the heat sink, thermal resistance was reduced by 30% under natural convection and by 20% under forced air convection. A similar conduction element between cylindrical batteries and heat pipes, was designed also by Wang et al. [126], highlighting the importance of the contact surface. A sensitivity analysis is performed, and it is found that the height of the conduction element is the most influential on temperature distribution, followed by the circumference angle of the contact surface, while the variation of thickness of the conduction element and battery spacing are insignificant regarding thermal homogeneity. increasing the heat transfer capacity. Experimental measurements were made adding fins at the end of the micro heat pipes and using forced air-convection. The study of Liu et al. [115] analyzes besides the enhanced thermal model also the cooling performance of ultrathin micro heat pipes combined with forced air-convection. The results show that for an effective heat pipe based BTMS it is recommended to have assistance by forced convection. The coupling of the two technologies can reduce maximum temperature by 7.1 °C at a 2 C discharge rate, compared to no heat pipes. Another system combining heat pipe with forced air cooling was developed by Tran et al. [125], where thermally conductive resin fills the spacing between the cylindrical batteries and the heat pipe with heat sink are placed on the lateral shell of the module, as shown in Figure 13b. Compared to using only the heat sink, thermal resistance was reduced by 30% under natural convection and by 20% under forced air convection. A similar conduction element between cylindrical batteries and heat pipes, was designed also by Wang et al. [126], highlighting the importance of the contact surface. A sensitivity analysis is performed, and it is found that the height of the conduction element is the most influential on temperature distribution, followed by the circumference angle of the contact surface, while the variation of thickness of the conduction element and battery spacing are insignificant regarding thermal homogeneity.
(a) (b) Figure 13. Hybrid BTMSs using heat pipe coupled with air cooling (a) schematic of heat pipe and fin BTMS for cylindrical batteries (adapted from source [123]); (b) battery module with resin matrix and heat pipe (adapted from source [125]).
A hybrid BTMS concept consisting of L-shaped heat pipes and a cooling plate is also considered by researchers. Yuan et al. [127] proposed heat pipe-copper plate structures for prismatic batteries, as presented in Figure 14a, while Gan et al. [128] designed wavy aluminum sleeves and circular heat pipes to fit easily inside a module of cylindrical batteries, as illustrated in Figure 14b. Although the latter study found that the length of the Figure 13. Hybrid BTMSs using heat pipe coupled with air cooling (a) schematic of heat pipe and fin BTMS for cylindrical batteries (adapted from source [123]); (b) battery module with resin matrix and heat pipe (adapted from source [125]).
A hybrid BTMS concept consisting of L-shaped heat pipes and a cooling plate is also considered by researchers. Yuan et al. [127] proposed heat pipe-copper plate structures for prismatic batteries, as presented in Figure 14a, while Gan et al. [128] designed wavy aluminum sleeves and circular heat pipes to fit easily inside a module of cylindrical bat- teries, as illustrated in Figure 14b. Although the latter study found that the length of the heat pipe's condenser section lowers maximum temperature and improves temperature uniformity, the former shows that the size of the evaporation section has a greater influence than the size of the condensation section. In both studies, coolant flow rate is shown to have an important influence on maximum temperature, but a smaller one on temperature distribution. heat pipe's condenser section lowers maximum temperature and improves temperature uniformity, the former shows that the size of the evaporation section has a greater influence than the size of the condensation section. In both studies, coolant flow rate is shown to have an important influence on maximum temperature, but a smaller one on temperature distribution.
(a) (b) Figure 14. Hybrid BTMSs using heat pipe coupled with cooling plate (a) L-shaped heat pipe-copper plates and liquid cooling plate (adapted from source [127]); (b) hybrid BTMS consisting L-shaped heat pipes with aluminum sleeves and liquid cooling plate (adapted from source [128]).
Using a single heat pipe at the most critical zone of a battery, Behi et al. [112] performed an analysis with cold plates at the lateral of the prismatic cells. It is found that compared to natural air cooling, the hybrid BTMS reduces maximum temperature by 32.6%. Alihosseini et al. [129] presented a BTMS where a water channel flows through the condenser section of the heat pipe, while the evaporation section is placed inside an aluminum plate between two battery rows for a better temperature distribution, as shown in Figure 15a. The presented system is able to maintain a constant surface temperature after three cycles, at the cost of 10% extra weight, compared to the situation without heat pipes.
The coupling of heat pipes and refrigerant lines is also analyzed, with the advantage of high heat transfer coefficient during phase change. The study of Liang et al. [130], performed on a BTMS for prismatic batteries, presented in Figure 15b, shows that although a reduced coolant temperature results in lower temperatures in the battery pack, it also leads to higher temperature differences. At the end of discharge process, lower coolant temperatures cause a significantly less uniform Li + concentration within the cells, this being the main cause of capacity loss. Reducing coolant temperature by 10 °C also lowers the battery module voltage, resulting in a capacity loss of 0.88%-1.17% at 5 C discharge. A similar system for prismatic batteries proposed by Yao et al. [131] can maintain the battery maximum temperature at preset values of even 25 °C, but with the cost of 38.41% lower coefficient of performance compared to a preset temperature of 35 °C. In this case also, the authors highlight the necessity to optimize the structure of refrigerant pipeline in order to reduce temperature non-uniformity in the battery module.
Other studies combine the high latent heat of PCMs and the advantage of temperature uniformity with the very high thermal conductivity of heat pipes, acting as a heat transfer enhancer of the PCM. The experimental study conducted by Zhao et al. [132] on the BTMS presented in Figure 16a shows that heat pipes lengthen the time needed for the melting of PCM during dynamic cycling and can improve temperature uniformity. The maximum temperature difference is 28.9% lower than in the case of a PCM-only BTMS. Figure 14. Hybrid BTMSs using heat pipe coupled with cooling plate (a) L-shaped heat pipe-copper plates and liquid cooling plate (adapted from source [127]); (b) hybrid BTMS consisting L-shaped heat pipes with aluminum sleeves and liquid cooling plate (adapted from source [128]).
Using a single heat pipe at the most critical zone of a battery, Behi et al. [112] performed an analysis with cold plates at the lateral of the prismatic cells. It is found that compared to natural air cooling, the hybrid BTMS reduces maximum temperature by 32.6%. Alihosseini et al. [129] presented a BTMS where a water channel flows through the condenser section of the heat pipe, while the evaporation section is placed inside an aluminum plate between two battery rows for a better temperature distribution, as shown in Figure 15a. The presented system is able to maintain a constant surface temperature after three cycles, at the cost of 10% extra weight, compared to the situation without heat pipes. Huang et al. [133] compared a system assisted by PCM around cylindrical batteries and coupled with air-cooling, a system assisted by PCM and coupled with liquid-cooling and a pure PCM BTMS, respectively. The high thermal conductivity of heat pipes coupled with liquid cooling allows the control of temperature difference below 3 °C, while the other systems, especially the pure PCM BTMS cannot transport heat to ambient quickly enough.
Another concept of hybrid BTMS coupling PCM and HP in series configuration is the positioning of heat pipes in contact with the batteries. Zhang et al. [134] proposed a separation type BTMS, as illustrated in Figure 16b. The generated heat is transferred outside the battery pack, stored in the PCM and dissipated to the environment. The proposed model can maintain maximum temperature below 45 °C at 4 C discharge rate, without an Figure 15. Hybrid BTMSs using heat pipe coupled with liquid cooling (a) water channel passing through the condenser section of heat pipes (adapted from source [129]); (b) condenser section in contact with coolant channel (adapted from source [130]).
The coupling of heat pipes and refrigerant lines is also analyzed, with the advantage of high heat transfer coefficient during phase change. The study of Liang et al. [130], performed on a BTMS for prismatic batteries, presented in Figure 15b, shows that although a reduced coolant temperature results in lower temperatures in the battery pack, it also leads to higher temperature differences. At the end of discharge process, lower coolant temperatures cause a significantly less uniform Li + concentration within the cells, this being the main cause of capacity loss. Reducing coolant temperature by 10 • C also lowers the battery module voltage, resulting in a capacity loss of 0.88%-1.17% at 5 C discharge. A similar system for prismatic batteries proposed by Yao et al. [131] can maintain the battery maximum temperature at preset values of even 25 • C, but with the cost of 38.41% lower coefficient of performance compared to a preset temperature of 35 • C. In this case also, the authors highlight the necessity to optimize the structure of refrigerant pipeline in order to reduce temperature non-uniformity in the battery module.
Other studies combine the high latent heat of PCMs and the advantage of temperature uniformity with the very high thermal conductivity of heat pipes, acting as a heat transfer enhancer of the PCM. The experimental study conducted by Zhao et al. [132] on the BTMS presented in Figure 16a shows that heat pipes lengthen the time needed for the melting of PCM during dynamic cycling and can improve temperature uniformity. The maximum temperature difference is 28.9% lower than in the case of a PCM-only BTMS.
Huang et al. [133] compared a system assisted by PCM around cylindrical batteries and coupled with air-cooling, a system assisted by PCM and coupled with liquid-cooling and a pure PCM BTMS, respectively. The high thermal conductivity of heat pipes coupled with liquid cooling allows the control of temperature difference below 3 • C, while the other systems, especially the pure PCM BTMS cannot transport heat to ambient quickly enough.
(a) (b) Figure 15. Hybrid BTMSs using heat pipe coupled with liquid cooling (a) water channel passing through the condenser section of heat pipes (adapted from source [129]); (b) condenser section in contact with coolant channel (adapted from source [130]).
Huang et al. [133] compared a system assisted by PCM around cylindrical batteries and coupled with air-cooling, a system assisted by PCM and coupled with liquid-cooling and a pure PCM BTMS, respectively. The high thermal conductivity of heat pipes coupled with liquid cooling allows the control of temperature difference below 3 °C, while the other systems, especially the pure PCM BTMS cannot transport heat to ambient quickly enough.
Another concept of hybrid BTMS coupling PCM and HP in series configuration is the positioning of heat pipes in contact with the batteries. Zhang et al. [134] proposed a separation type BTMS, as illustrated in Figure 16b. The generated heat is transferred outside the battery pack, stored in the PCM and dissipated to the environment. The proposed model can maintain maximum temperature below 45 °C at 4 C discharge rate, without an auxiliary fan, which is only recommended when longtime high discharge rates are necessary. Jiang et al. [89] suggest the choice of a melting point for the PCM at least 3 °C higher than the ambient temperature for a full recovery of the PCM, and a heat transfer coefficient at the condensation section of the heat pipe between 30 and 60 W/m 2 K for a sustainable operation. Zhou et al. [135] presented a BTMS where the heat pipes are concentrated in a radiator, enhanced with aluminum fins, and the batteries are immersed in a flame retardant, insulated phase change liquid, proving better thermal performance compared to forced air-cooling.
(a) (b) Figure 16. Hybrid BTMSs using PCM coupled with HP in series configuration (a) PCM assisted by HP and forced convection (adapted from source [132]); (b) HP assisted by PCM and forced convection (adapted from source [134]). Figure 16. Hybrid BTMSs using PCM coupled with HP in series configuration (a) PCM assisted by HP and forced convection (adapted from source [132]); (b) HP assisted by PCM and forced convection (adapted from source [134]).
Another concept of hybrid BTMS coupling PCM and HP in series configuration is the positioning of heat pipes in contact with the batteries. Zhang et al. [134] proposed a separation type BTMS, as illustrated in Figure 16b. The generated heat is transferred outside the battery pack, stored in the PCM and dissipated to the environment. The proposed model can maintain maximum temperature below 45 • C at 4 C discharge rate, without an auxiliary fan, which is only recommended when longtime high discharge rates are necessary. Jiang et al. [89] suggest the choice of a melting point for the PCM at least 3 • C higher than the ambient temperature for a full recovery of the PCM, and a heat transfer coefficient at the condensation section of the heat pipe between 30 and 60 W/m 2 K for a sustainable operation. Zhou et al. [135] presented a BTMS where the heat pipes are concentrated in a radiator, enhanced with aluminum fins, and the batteries are immersed in a flame retardant, insulated phase change liquid, proving better thermal performance compared to forced air-cooling.
Parallel configurations of PCM and HP are also analyzed. Wang et al. [136] proposed a structure that utilizes completely the spaces between the cylindrical batteries and does not increase the volume of the battery pack, while cost is maintained relatively low, by using aluminum both at the heat pipes and fins. The experimental measurements show that the heat pipe improves thermal conductivity inside the battery and lowers the maximum temperature, while the PCM is more effective at reducing the temperature non-uniformity. Another study [137] proposes the use of a pulsating heat pipe between two prismatic batteries and filling the spaces inside the serpentine channels with PCM, while immersing the condenser section in 25 • C water. The study recommends the PCM melting temperature to be higher than the saturation temperature of the PHP's working fluid, and the positioning of the adiabatic section outside the region of the terminals.
Another possibility is the use of TECs as the active part, which are a relatively new technology for the use in BTMSs and are usually used as a support device for the system, coupled with other methods, with the aim of increasing the efficiency of existing BTMSs. The working principle of TECs is the conversion of voltage to a temperature difference, and one its greatest advantages is the possibility of using it for both cooling and heating purposes [138]. Moreover, its cold side can be lowered below ambient temperature [139]. In the design of TEC-based BTMSs, the control of the feeding electric current must be chosen in conformity with the thermal load of the system. Higher currents are an inefficient method of absorbing larger amounts of heat, because of the drop of the coefficient of performance [140].
Li et al. [141] also proposed a TEC-based hybrid BTMS, presented in Figure 17a and compared three systems: only TEC, TEC starting at 40 • C coupled with forced air cooling, and TEC coupled with forced air cooling both starting at 40 • C. The less complex pure TEC system presents very good thermal performance compared to the hybrid systems as well, but with larger energy waste. Considering cooling performance and energy consumption, the hybrid system with only the TEC starting at 40 • C is the most suitable. Zhang et al. [142] proposed the positioning of the TEC modules between the batteries and water-cooling plates and investigated the best cooling structure and water inlet velocity. Lyu et al. [143] proposed an advanced hybrid BTMS design, illustrated in Figure 17b. A liquid coolant removes heat from the cylindrical batteries via a copper holder interface for safety reasons. The cold side of the TEC is in contact with the liquid, while heat is removed from the hot side of the TEC using forced air cooling. The study shows that the proposed system has great potential, experimental measurements presenting a 43 • C temperature drop, from 55 to 12 • C, on the battery surface, at 12 V supply of the TEC.
Parallel configurations of PCM and HP are also analyzed. Wang et al. [136] proposed a structure that utilizes completely the spaces between the cylindrical batteries and does not increase the volume of the battery pack, while cost is maintained relatively low, by using aluminum both at the heat pipes and fins. The experimental measurements show that the heat pipe improves thermal conductivity inside the battery and lowers the maximum temperature, while the PCM is more effective at reducing the temperature non-uniformity. Another study [137] proposes the use of a pulsating heat pipe between two prismatic batteries and filling the spaces inside the serpentine channels with PCM, while immersing the condenser section in 25 °C water. The study recommends the PCM melting temperature to be higher than the saturation temperature of the PHP's working fluid, and the positioning of the adiabatic section outside the region of the terminals.
Another possibility is the use of TECs as the active part, which are a relatively new technology for the use in BTMSs and are usually used as a support device for the system, coupled with other methods, with the aim of increasing the efficiency of existing BTMSs. The working principle of TECs is the conversion of voltage to a temperature difference, and one its greatest advantages is the possibility of using it for both cooling and heating purposes [138]. Moreover, its cold side can be lowered below ambient temperature [139]. In the design of TEC-based BTMSs, the control of the feeding electric current must be chosen in conformity with the thermal load of the system. Higher currents are an inefficient method of absorbing larger amounts of heat, because of the drop of the coefficient of performance [140].
Li et al. [141] also proposed a TEC-based hybrid BTMS, presented in Figure 17a and compared three systems: only TEC, TEC starting at 40 °C coupled with forced air cooling, and TEC coupled with forced air cooling both starting at 40 °C. The less complex pure TEC system presents very good thermal performance compared to the hybrid systems as well, but with larger energy waste. Considering cooling performance and energy consumption, the hybrid system with only the TEC starting at 40 °C is the most suitable. Zhang et al. [142] proposed the positioning of the TEC modules between the batteries and water-cooling plates and investigated the best cooling structure and water inlet velocity. Lyu et al. [143] proposed an advanced hybrid BTMS design, illustrated in Figure 17b. A liquid coolant removes heat from the cylindrical batteries via a copper holder interface for safety reasons. The cold side of the TEC is in contact with the liquid, while heat is removed from the hot side of the TEC using forced air cooling. The study shows that the proposed system has great potential, experimental measurements presenting a 43 °C temperature drop, from 55 to 12 °C, on the battery surface, at 12 V supply of the TEC.  Another TEC-based hybrid BTMS concept consists of coupling PCM with TEC modules, which can heat or cool the PCM so that it regains its capacity to retain thermal energy. Song et al. [144] considered a conventional PCM module surrounding the batteries and recommend positioning the TECs on the facing sides that are closer to each other. Jiang et al. [145] added liquid cooling to this concept, as presented in Figure 18a. The study shows that for a battery heat generation of 6 W, the temperature of 50 • C is reached within al. [145] added liquid cooling to this concept, as presented in Figure 18a. The study shows that for a battery heat generation of 6 W, the temperature of 50 °C is reached within 930 s under natural convection, 1275 s under liquid cooling and 5335 s using the proposed hybrid system.
Further, Zhang et al. [146] designed a BTMS combining a heat pipe with TEC, presented in Figure 18b. In this system, the heat on the hot side of the TEC is rapidly conducted to the heat sink region of the BTMS, reducing the surface temperature rising rate. The investigation shows that for the analyzed battery, at a 1 C discharge rate no BTMS is necessary, at 1.5 C a heat pipe based BTMS is needed, at 2 and 2.5 C the combination of heat pipe and TEC is required, while at 3 C discharge duration must be limited.

Discussion and Conclusions
This paper presents comprehensive, detailed information of the design, functionality, performances and applicability of basic and hybrid BTMSs to different Li-ion battery types. The main classification of basic systems is divided into four groups (air, liquid, phase change material, heat pipe), while for the hybrid systems numerous combinations are presented, the most important ones being as follows: • PCM + air/liquid; • HP + air/liquid; • PCM + HP (+ air/liquid); • TEC + others.
A first meaningful discussion is regarding the advantages and disadvantages of each individual element that can be used in BTMSs, being it passive or active, listed in Table 1. Further, the characteristics regarding their implementation in BTMSs and adaptability to the different battery shapes are detailed in Table 2. The disadvantages of basic BTMSs are reduced through combining them with other systems into hybrid BTMSs, presented in Table 3 and Table 4, or various other methods. Further, Zhang et al. [146] designed a BTMS combining a heat pipe with TEC, presented in Figure 18b. In this system, the heat on the hot side of the TEC is rapidly conducted to the heat sink region of the BTMS, reducing the surface temperature rising rate. The investigation shows that for the analyzed battery, at a 1 C discharge rate no BTMS is necessary, at 1.5 C a heat pipe based BTMS is needed, at 2 and 2.5 C the combination of heat pipe and TEC is required, while at 3 C discharge duration must be limited.

Discussion and Conclusions
This paper presents comprehensive, detailed information of the design, functionality, performances and applicability of basic and hybrid BTMSs to different Li-ion battery types. The main classification of basic systems is divided into four groups (air, liquid, phase change material, heat pipe), while for the hybrid systems numerous combinations are presented, the most important ones being as follows: • PCM + air/liquid; • HP + air/liquid; • PCM + HP (+ air/liquid); • TEC + others.
A first meaningful discussion is regarding the advantages and disadvantages of each individual element that can be used in BTMSs, being it passive or active, listed in Table 1. Further, the characteristics regarding their implementation in BTMSs and adaptability to the different battery shapes are detailed in Table 2. The disadvantages of basic BTMSs are reduced through combining them with other systems into hybrid BTMSs, presented in Tables 3 and 4, or various other methods.     Air cooling represents the most easily applicable solution for the thermal management of batteries. However, the low specific heat of air still is a major concern when considering its application in full-size battery packs. A possible solution could be the use of active cooling and heating, using an evaporator or heater core. This would allow a reduced dependency on ambient conditions, an increased control over the inlet temperature in the battery pack and therefore a more efficient management of the temperatures inside the pack. Although such a method would increase the heat transfer, the occupied large volume still represents a significant disadvantage. The inlet and outlet scheme and the air flow configuration are also promising research directions, with numerous possible combinations. However, the possibility of implementing these schemes in a full-size battery pack should also be investigated, considering the need to implement the required inlet and outlet channels in the tight package of a vehicle. Regarding the uneven temperature distribution of air cooling, this can be reduced by using air flow channels with variable widths or by distribution pipes.
The large use of prismatic and pouch batteries in commercial vehicles means that liquid cooling plates are one of the most implemented cooling systems in electric vehicles, due to their good heat capacity and compact arrangement. Most research studies focus on various parameters regarding heat transfer and pressure drop, but only a few consider the uneven heat generation in the batteries, which is even more significant for the larger prismatic and pouch batteries. Therefore, more attention in future studies should be accorded to the possibility of implementing bidirectional or reciprocating liquid flows, considering the local heat generation, therefore improving the temperature distribution. In the presented study, examples of improving the temperature uniformity offered by cooling plates are presented using double-layered plates and a better mass distribution in the channels. The heat transfer capacity of conventional liquid cooling systems can further be enhanced using PCM emulsions or liquid metal. For cylindrical batteries, the poor adaptability to indirect liquid cooling is represented by the use of wavy channels, half helical ducts or conduction elements, which significantly increase complexity and add weight to the system.
Adapting the existing refrigerant circuit of the vehicles to serve also for the thermal management of batteries is a concept with a relatively high probability of practical application. The presented studies clearly show the potential of such systems, having a high heat transfer capacity, a capability to cool below ambient temperature and the possibility of heating. Despite these advantages, there are no mentioned studies to analyze the sizing of such a system for a full-size battery pack, nor the necessary changes to the existing air conditioning system to cope with the additional thermal load of the battery. More effort should be directed towards the development of a large-scale system, with precise control algorithms for cabin comfort, battery maximum temperature and temperature uniformity.
The main disadvantages of pure PCMs, namely the low thermal conductivity and leakage problems can be tackled by adding other materials such as EG, AlN, nanosilica, copper mesh foam and others, to form CPCMs. Although the enhancement of pure PCM with additional materials in different concentrations is a promising research direction, the preparation processes are very laborious, as well as time and energy consuming. It could be of interest by how much the preparation of large quantities of CPCM increases the overall costs of full-size battery packs. Moreover, in the case of pouring hot, melted CPCM directly around the batteries, further investigations are needed regarding the possible negative effects on Li-ion batteries.
Another disadvantage, the limited capacity of storing thermal energy as latent heat, represented by the complete melting of PCM, can be solved by coupling with air or liquid cooling for a better heat dissipation and recovery of the PCM. At present, the main reasons that stop the large-scale applicability of PCMs are the leakage problems, with a negative impact on the safety and reliability of the system and the high additional mass.
Heat pipes coupled with air, water, or refrigerant cooling present good heat transfer efficiency due to the excellent thermal conductivity of HPs, but to also obtain temperature uniformity, it is possible to further add PCM. Research regarding alternative materials in the construction of heat pipes should be closely followed. The possible replacement of copper with aluminum would significantly reduce the costs and could also save some weight, but more investigations in this direction are necessary to determine by how much would the thermal conductivity be affected.
The low COP of TECs means that they are usually used as a supporting device to existing BTMSs, with the advantages of high cooling capacity and possibility of heating. A BTMS using only TEC could be applicable after a significant increase of the energy efficiency of these. In order to be economically competitive with other systems, the thermoelectric figure of merit should increase approximately three times, which for the moment is only possible by using sophisticated and expensive nanostructures.
Both HPs and TECs increase the overall costs of the system, slowing down their large-scale applicability in electric-driven vehicles. However, the use of cheaper alternative materials in the construction of these components could result in an intensified interest in HP-based and TEC-based BTMSs.
The results and good practices obtained so far in the design, configuration and development of BTMSs for Li-ion batteries can further become a starting point for studies related to the thermal management of fuel cells. It is well known that at least at European level the massive use of hydrogen for the supply of electricity is desired in future years, with the help of technologies based on fuel cells for electricity generation.
The general conclusion that can be issued is that considering the great variety of types, design and approaches regarding BTMSs, researchers in this field must pay attention to future research and developments in this domain. Applying and integrating the proposed novelty systems and models could improve the thermal and energy efficiency of full-size battery packs.