1. Introduction
The lifetime of lithium-based battery cells is strongly correlated to the temperature levels to which they are exposed during operation and storage. Heat is dissipated in a battery cell both during charging and discharging and its intensity depends on the chemistry of the cells as well as on the rate of this process, which is generally characterized in terms of the so-called C-rate (i.e., the ratio between the actual charge/discharge current and the theoretical current draw under which the battery would deliver its nominal rated capacity in one hour). To improve battery life, it is necessary to keep its operating temperature within a prescribed optimal range and minimize the temperature gradients within a cell as well as the battery pack. Excessively high or low temperatures can reduce lifetime, increase safety risks and even cause permanent damage to the battery [
1].
Electric vehicles (EVs) are equipped with battery packs consisting of hundreds or even thousands of battery cells connected in a specific serial-parallel configuration. As the cells aging rate depends on the operating temperature, a nonuniform temperature distribution within the battery pack (temperature gradient) may lead over time to a significant discrepancy in the characteristics of the cells, thereby affecting the overall performance of the battery pack. For example, as discussed in Ref. [
1], a single poorly performing cell may lead to insufficient charging and discharging of the whole battery pack.
The optimal operating temperature of battery packs depends on the chemistry of the cells, the required performance and the desired lifetime. In literature, the optimal temperature range indicated for Li-ion battery cells differs. Refs. [
1,
2] indicate a range between 15 °C and 35 °C, while Ref. [
3] between 20 °C and 40 °C, and Refs. [
4,
5] between 20 °C to 50 °C. Regarding the temperatures gradients within a battery pack, there is, instead, general agreement that differences in the operating temperature of the single cells should not exceed 5 °C [
4]. It is, then, apparent that the thermal management of the battery pack is key to ensuring safe operation and avoiding fast degradation of the cell performance.
It is possible to distinguish between three battery thermal management (BTM) system configurations [
6]:
Passive: heat dissipation occurs without energy consumption;
Active: a (gas or liquid) coolant is forced to circulate through the battery pack;
Combined or hybrid: battery cooling is performed by a combination of both passive and active cooling solutions.
In [
7], Huaqiang Liu et al. documented the state-of-the-art developments in recent years concerning various BTM systems for batteries. The authors list the advantages and disadvantages of four types of active cooling systems based on the coolant characteristics as:
Air cooling (forced convection)
Liquid cooling
Evaporative cooling (heat pipe)
Melting of a phase change material
Passive cooling has several advantages, such as no power consumption, no additional control systems required, construction and implementation simplicity, no additional safety risks, low maintenance needs and low-cost potential. However, it is very challenging to fulfil all the requirements in terms of cooling demand and temperature uniformity of battery packs for EVs only by means of passive cooling. One of the most promising passive solutions is the use of a phase changing material (PCM). A PCM stores and releases thermal energy via a melting−freezing cycle: when the material freezes, it releases thermal energy to the surrounding environment and when it melts, an equal amount of thermal energy is absorbed. In theory, the temperature of the PCM remains constant during phase transition. In practice, due to the unavoidable presence of impurities in commercially available PCMs, the solid−liquid phase transition occurs with a temperature glide.
In the literature, the use of a PCM in Li-ion batteries has been investigated experimentally as well as by numerical simulations. The PCM allows the temperature rise in the battery pack to be reduced (below the upper safety limit) and a homogeneous temperature distribution inside the pack and the cell itself to be maintained. Ling et al. [
4] have tested experimentally two different thermal management systems for a module made of 20–18,650 Samsung cells:
The tests showed that, if the battery pack was charged at 1-C rate and discharged at 1.5/2-C rate, the cell temperature rose above 60 °C with the passive BTM solution only after two charging/discharging cycles. A room temperature of 25 °C was maintained during the tests. In the contrast, with the hybrid solution, the temperature of the battery pack did not exceed 45 °C. It must be noted that the thermal capacity of the PCM/EG composite was high since the cells were arranged in the composite with a square pattern at a distance of 12 mm from one another.
Hemery et al. [
8] experimentally assessed the performance of a BTM system adopting a PCM and that of an air cooled BTM in the case of natural and forced convection. The battery pack mock-up was equipped with 27 electric heaters embedded in cylindrical stainless-steel cases to reproduce the heat dissipation occurring in the battery cells. The authors measured the temperature rise within the battery pack for the two BTM concepts during multiple subsequent driving cycles. Large temperature variations were measured with the air-cooled BTM, while with the use of PCM, the increase over time in the average cell wall temperature was significantly lower. As long the PCM was not fully melted, the temperature increase was around 1 °C. Similarly, the temperature gradients within the mock-up were largely reduced by the PCM.
W.Q. Li et al. [
9] studied a BTM concept where the battery cells were packed in a structure made of copper metal foam saturated with paraffin (RT 44HC). The authors performed experiments to assess this BTM concept versus a passively air-cooled BTM system. The tests showed that, with the air-cooled battery pack, the cells located in the center of the module reached the upper temperature limit indicated by the manufacturer (i.e., 65 °C) under a 1-C discharge rate after approximately 50 min. This temperature limit was not exceeded by the outer cells of the battery pack, thanks to the cooling provided by ambient air through natural convection. The temperature difference (DT) between the inner and outer cells was higher than 10 °C. The battery pack structure combining copper metallic foam with paraffin featured, instead, a maximum temperature lower than 50 °C and the temperature difference across the pack reduced to about 4 °C for the same battery discharge rate. If the cells were surrounded only by paraffin, the maximum temperature of the inner cells increased to 53 °C, while that of the outer cells was around 47 °C (DT = 6 °C). It was concluded that the metallic foam increased the thermal conductivity within the PCM, thus improving the performance of the passive BTM system.
Similar results were found by Guiwen Jiang et al. [
10], who investigated the adoption in the battery pack of a composite made of a phase change material and expanded graphite (CPCM). The authors also varied the amount of expanded graphite (EG) in the composite and observed that the leakage of melted PCM occurring during phase transition tended to reduce as the EG mass fraction was increased. Without PCM/CPCM, the surface temperature of a cell increased almost linearly over time and reached a temperature higher than 60 °C in less than one hour. With PCM, the cell temperature started increasing only after 700 s, when the PCM that was in contact with the cells had completely melted. Due to the poor thermal conductivity of the material, the PCM outer layers did not melt. The expanded graphite of the composite material allows the thermal conductivity to be enhanced. The outcome was a more homogeneous temperature within the battery pack and, therefore, a more effective BTM system. The amount of EG in the composite must be, nevertheless, limited since it decreases the overall specific enthalpy variation of the CPCM during phase transition. The authors concluded that the EG fraction has to be limited to around 16–20 wt.%. This represents a trade-off between the reduction in PCM leakage and the decrease in the overall thermal inertia of the composite.
In [
11], Huang et al. studied three different kinds of BTMs. In the baseline solution, 30 cylindrical 18,650 lithium iron phosphate batteries were inserted in a PCM-graphite matrix. The other BTMs made use of flat heat pipes to cool the PCM-graphite matrix. The working fluid of the heat pipes was either air (BTM solution No. 2) or ethyl alcohol (BTM solution No. 3). During a single full discharge test, the melting temperature was not reached for all the BTM configurations, even for a 2-C discharge rate. If the discharge/charging cycle was repeated, the temperature of the PCM in the uncooled configuration progressively increased, while for the BTM systems using heat pipes the average temperature tended to stabilize. With a 3-C discharge/charge rate, the cell temperature was around 50 °C for the BTM system adopting air as the working fluid of the heat pipes and 45 °C in the ethyl alcohol case.
Based on the results reported in the literature, it is possible to conclude that the adoption of a PCM in the battery pack is not sufficient to maintain cell temperature low enough to preserve their performance and lifetime. In this respect, M. Alipour et al. [
5] report that, for the majority of Li-ions batteries, the optimal operating range is between 10–35 °C.
This work presents the design of a combined active/passive BTM solution for the battery pack of an electric van. The research was carried out in the framework of the European funded project Everlasting, that aims at enhancing the range, the lifetime and the safety of Li-ion batteries through innovative battery management. The proposed BTM system combines the use of PCM in the battery pack with air cooling. The cooling air is precooled to maintain the battery temperature below 35 °C under all operating conditions. The air is cooled by a refrigerating loop whose realization requires the adaptation of the heat pump installed aboard the eVan for cabin air conditioning. The design of the heat pump is, however, not documented in this paper, which is structured as follows:
Section 2 provides an overview of the application and of the design requirements of the battery pack.
Section 3 documents the CFD model of the battery pack used to assess the hybrid BTM system and the experimental setup realized to validate the numerical results.
Section 4 compares the performance of the fully passive solution against that of the hybrid BTM system for a prescribed driving cycle of the eVan.
3. Modeling of the Hybrid BTM System
To assess the combined active/passive BTM solution of the battery pack, a combination of modelling and experimental activities were carried out. Computational fluid dynamics (CFD) simulations were performed to calculate the temperature variations of the cells and to simulate the effect of two cooling air channels at the top and the bottom of the battery pack on the temperature gradient among the cells. An experimental mock-up was, then, built to test the passive and the active cooling solutions, and thus validate the numerical calculations. The CFD modelling approach and the experimental setup are discussed in the following.
3.1. CFD Model
The CFD model was developed in ANSYS Fluent. A grid independence study was carried out and showed that a mesh size of 2.1 × 10
6 elements was the most suitable for this study. A portion of the mesh domain is shown in
Figure 7a. The simulation domain encompassed only a portion of a battery pack module, as symmetry boundary conditions could be applied at planes parallel to the main flow direction of the cooling air and passing through the vertical axes of the cells, see lateral boundaries of the CFD model in
Figure 8a. Moreover, it was also possible to identify a symmetry plane with respect to the vertical axes of the cells: only half of the cells required, then, to be modeled.
For the purely passive thermal management solution, the standard pressure-based solver was used. The considered time step was equal to 5 s while the simulation was deemed converged when the scaled energy residual was lower than 1 × 10-11. When air cooling was simulated, the standard k-eps model was used. The time-step of the solution varied from 0.05 s at the beginning of the simulation and it was progressively increased up to 1 s. The adopted convergence criterion required that the scaled residual for the continuity equation decreased to 1 × 10-3.
The PCM-graphite matrix surrounds the cells along 45 mm of their height. The lower and upper ends of the cells extends into the cooling air channels by 10 mm, as shown in
Figure 7b.
Figure 8b shows an exemplary temperature distribution in the battery pack module calculated by the CFD model at the end of a charge/discharge cycle performed at a 1-C rate.
3.2. Experimental Setup
To validate the calculations of the CFD model and to assess the actual performance of the PCM-EG matrix, an experimental setup and a mock-up of the battery pack were built. A schematic of the experimental setup is shown in
Figure 9. The prescribed temperature of the cooling air mass flow rate was attained in all operating conditions via a heat exchanger located downstream of the air compressor. The mockup was instrumented as follows:
Two air flowmeters, one for each cooling channel.
Absolute and differential pressure sensors at the inlet and at the outlet of the mock up.
Forty thermocouples type K to measure:
- o
External temperature of the cells at the center of the PCM—graphite matrix (position 1 in
Figure 9) and in the cooling channels (position 2)
- o
Temperature inside the PCM—graphite matrix (position 3)
- o
Air temperature in the cooling channel at the top of the battery pack (position 4)
Figure 10 shows two photos of the setup. The first one displays one of the cooling air channels highlighting the direction of the airflow (a), while the second one shows the whole mockup with the instrumentation cabling (b).
The battery pack mockup does not feature real cells. As shown in
Figure 11, these were replaced by electric heating elements in an aluminum cylinder. The current of the heating elements could be adjusted to reproduce the heat dissipation of the battery cells. The thermal power that was transferred from the heating elements to the PCM matrix and the cooling air could be accurately estimated by measuring the electrical power. The aluminum casings also featured two grooves along the external surface in contact with the PCM-graphite matrix. Thermocouples were inserted into the grooves to measure the temperature of the battery pack mockup. The heating elements were arranged in the PCM matrix with a triangular pattern. To reduce the overall mass and size of the mockup, the space between two cells/heating elements was limited to 2 mm. The height of each heating element—aluminum cylinder assembly was 65 mm for a total weight of 53 g, 4 g more than the actual electric cells.
Two kinds of tests were performed with the battery pack mock-up. The first consisted of heating up the HSS till the PCM was completely melted. The temperature of the composite was measured at different locations to estimate the melting time and the temperature gradients within the battery pack mock-up. In the second test, the battery pack was cooled by an airflow at 30 °C, starting from an initial condition where the PCM was already completely melted. The solidification time and the temperature gradients in the mock-up were estimated based on the thermocouple measurements.
3.2.1. PCM Melting Time
The PCM melting time was monitored for different values of thermal dissipation of the cells, ranging from 0.55 W to 6.6 W/cell (corresponding to a 0.5-C–4-C discharge rate range). The experimental results and the predictions of the CFD model are shown in
Figure 12 and
Figure 13, respectively.
It is apparent that there is a large deviation in the actual and predicted melting time. This was confirmed also for the other thermal power levels, as reported in
Table 4.
The difference between the experimental and the simulated melting time strongly increases, the lower the thermal power released by the cells/heating elements. Experimental tests were repeated several times and the same trend was observed. The cause thereof lay in the heat losses occurring at the top and bottom surfaces of a cell mock-up which were in direct contact with the PEEK casing that absorbed part of the thermal power released by the heating elements. No clearance between the heating elements and the external casing was considered in the design as preliminary CFD calculations showed an inefficient flow distribution when there was a gap between the top/bottom surfaces of the cells and the external PEEK casing. A possible solution to reduce these heat losses may consist of adding an insulation layer between the upper surface of the cooling air channels and the two ends of a heating element. However, this modification could not be put in place due to the very tight construction of the PEEK casing.
3.2.2. Thermal Gradients in the Battery Pack
Figure 14,
Figure 15 and
Figure 16 show the wall temperatures of the cells when the thermal dissipation is, respectively, equal to 0.55 W, 2.2 W and 6.6 W. The figures report also the maximum temperature difference between the cells over time. This quantity is indicated as
DT max-min in the charts and is deemed representative of the thermal gradients among the cells. It is calculated as the difference between the maximum and minimum temperature of the cells in the battery pack mock-up.
In all the experimental tests reported in the figures above, DT max-min was less than 2.5 °C as long as the PCM was not completely melted in a given position. When the thermal power of the heating elements was low (0.55 W), the thermal energy was evenly transferred by heat conduction throughout the PCM matrix and the value of DT max-min remained almost constant, except when there were inhomogeneities in the state of the PCM over the mock-up: in the zones where PCM melting was completed the temperature of the cells rose rapidly. Once the PCM was fully melted throughout the composite matrix, the value of DT max-min lowered to levels similar to those in the initial phase of the experiment. When the thermal power of the heating elements was higher (6.6 W), the value of DT max-min tended to increase more rapidly if the PCM was completely melted. However, in all test cases, the temperature differences among the cells remained below 5 °C. This is a key feature enabled by the adoption of the EG-PCM matrix in the battery pack.
3.2.3. Thermal Gradients in the Battery Pack when Air Cooling Is Activated
As previously discussed, the amount of PCM of the HSS is not sufficient to maintain the battery pack at temperatures lower than 40 °C in all operating conditions. Thus, a cooling system is needed. In the experiments described in this section, air cooling was activated to cool down the cells from an initial temperature of 44 °C, while the heating elements were turned off. A possible situation where operating conditions similar to those of this experiment may occur is the charging of the eVan battery pack. Future charging stations may, indeed, feature a cooling system to cool down to the battery pack once the charging is completed.
Figure 17,
Figure 18 and
Figure 19 show the measured surface temperature of the different cells in the battery pack mock-up when the airflow speed was 0.5 m/s, 1 m/s and 2 m/s. The cooling air temperature at the battery pack inlet was kept constant and equal to 30 °C. The chart of the figures also report the variation of
DT max-min over time, see the red curve in the charts.
As expected, the higher the air velocity, the shorter the cooling period. At the same time, the temperature differences among the cells tended to increase, though only slightly. The maximum value of
DT max-min never exceeded 4 °C. The experimental tests also showed that melting and solidification of the PCM did not occur over the same temperature range, as already observed in the literature, e.g., in Ref. [
17]. The melting of the PCM considered in this project occurs around 37 °C while its solidification starts at about 34 °C.
For this kind of experiment, the CFD calculations are in good agreement with the measurements.
Figure 20 shows the cell temperature when the battery mock-up was cooled by an airflow characterized by a speed of 1 m/s. The same figure also displays the evolution over time of
DT max-min, that is, the difference between the minimum and maximum temperature measured in the battery pack, when the airflow speeds were 0.5 m/s, 1 m/s and 2 m/s.
The trend of DT max-min determined by CFD simulations is similar to that observed in the experiments. For instance, with reference to the case where the airflow velocity was equal to 1m/s, the maximum values of DT max-min according to the CFD model and the estimated one based on the thermocouples measurements are both close to 3.5 °C. The CFD model tends to predict a more homogeneous temperature distribution than that found experimentally. The value of DT max-min remains above 1 °C during PCM solidification for about 3500 s in the experiment against 2000 s of the CFD simulation. At the same time, the numerical results shows that about 7000 s are required to cool down the mock-up to about 30 °C. This estimate is very close to that observed in the experiment. The time required for PCM solidification is also comparable in the test and the simulation, though the temperature range at which this phenomenon occurred differs by about 6 °C. Similar conclusions can be drawn also for the other two test cases.
3.3. CFD Model Validation
The comparison of the simulation results with the experimental data previously presented revealed that:
- (1)
The melting time experimentally observed was twice as high as that predicted by the CFD model. Given that such a large discrepancy was not observed when the battery pack was cooled down and the PCM solidified, the melting time estimated on the basis of the thermocouple measurements was, arguably, highly influenced by the thermal losses occurring through the mock-up casing as the heating elements were in direct contact with the external walls of the casing. For the same reason, the temperature differences among the cells predicted by the CFD model did not match those experimentally observed during PCM melting.
- (2)
In the contrast, when the battery pack mock-up was cooled down and the heating elements were turned off, the CFD model captured reasonably well the trend over time of the temperature difference among the cells and the time required to cool down the mock-up.
- (3)
The temperature range over which the solidification of the PCM occurred was about 6 °C lower than that in the numerical simulations. The reason thereof is that the PCM, when cooled from a superheated liquid state, reaches a subcooled metastable state before the melting process starts. This nonequilibrium thermodynamic state has been observed for different PCMs during solidification [
11]. The simple thermodynamic model implemented in the CFD model cannot account for that.
Regarding the last consideration,
Figure 21 compares the equivalent specific heat capacity of the composite considered in the model with that estimated based on the experimental data. The latter can be evaluated by dividing the thermal energy provided by the electric heaters over a given time interval, here taken equal to 15 s, by the mass of the composite and the average temperature variation measured during the same period.
Based on the considerations above, it can be argued that the CFD model correctly predicts the heat conduction within the battery pack and the heat transfer between the cooling air and the composite matrix. Thus, it can be considered validated.
5. Conclusions
The study presented in this paper investigated an active/passive solution for the battery pack of an electric van. In the first part of this study and also in
Section 4.2, it was shown that, a purely passive BTM solution employing a PCM-expanded graphite composite matrix would require a large amount of PCM in order to keep the temperature of the battery pack at or below 40 °C in all the driving conditions of the electric van. It was, then, decided to couple this passive solution with an air cooling system that can be activated when necessary. Notice that the cooling system may be necessary also in the case the capacity of the HSS were enough to maintain the temperature of the cells in the prescribed range. The PCM is, indeed, to be cooled down to return into solid phase. This can be accomplished, for example, when the vehicle is parked after a drive.
The assessment of such a hybrid BTM configuration has been performed based on a combination of numerical analysis and experiments in ad-hoc setup. A mock-up of the battery pack was built: a PCM-expanded Graphite matrix houses 40 electric heating elements whose current can be modulated to reproduce the thermal dissipation of the battery cells. Simulations and experimental tests were performed and compared for different cases.
The melting time experimentally observed was higher than that predicted by the CFD model, but this discrepancy was not observed during the solidification of the PCM. This discrepancy could be attributed to the thermal losses occurring through the mock-up casing as the heating elements were in direct contact with the external walls of the casing. Moreover, the temperature range over which the PCM melted was 6 °C lower than that of the numerical simulations. This occurred because the simple thermodynamic model cannot predict the metastable state reached by the liquid phase before the onset of solidification.
The mockup was then used to emulate the heat dissipation of the cells during a reference driving cycle of the eVan and the thermal management solution as designed. Results showed that considering a highway driving cycle, even starting from an initial temperature of the composite of 40 °C the battery pack temperature can be maintained below 40 °C by a cooling air flow characterized by a velocity of 2.5 m/s and a temperature of 30 °C. Moreover, the temperature gradient within the battery pack remains below 5 °C, meeting the design requirements for the application at hand.