1. Introduction
Currently, the world faces increasingly severe environmental problems, including the greenhouse effect and industrial pollution [
1]. It is also beset by the shortage of resources, especially fossil fuel energy. To resolve these urgent problems, people consider green cars including battery electric vehicles (BEVs), fuel cell electric vehicles (FCEVs), hybrid EVs (HEVs), and so on, to be critically important [
2]. As key components of green cars, batteries have attracted great attention in recent years. Among the current commercial batteries, lithium-ion batteries (LIBs) are considered as the best choice for automobile application, due to their high power/energy density, long life cycle, safety, and environmental friendliness [
3,
4,
5]. However, safe operating conditions are hard to ensure for LIBs, so the battery management system (BMS) is essential when LIBs are used in vehicles [
6,
7]. To monitor and better manage LIBs in the battery pack, the pack voltage, pack current, cell voltage, cell temperature, cell impedance and so on must be detected.
Battery impedance has been regularly employed to investigate battery parameters and states. Buller et al. [
8] proposed an impedance-based simulation model for supercapacitors and LIBs, and the results showed that the simulated and measured voltage data were in excellent agreement. Xu et al. [
9] analyzed the impedance spectra of LIBs and derived a simplified impedance model for the estimation of battery state of charge (SOC). Based on the different value of the simplified equivalent circuit elements at different SOCs, Westerhoff et al. [
10] developed a novel method to estimate the battery SOC, and this method could be considered as complementary to the ampere-hour integral method. Moreover, battery impedance has been employed to estimate the battery’s internal temperature and determine the thermal model parameters [
11,
12,
13,
14].
Considering the importance of battery impedance, many researchers have studied the impedance characterization during practical applications. Schönleber et al. [
15] derived the impedance of a LIB electrode from the intercalation reaction of lithium and studied its dependency on SOC, and the results showed that the electrolyte diffusion impedance is independent of the SOC. Huang et al. [
16] systematically analyzed how different SOCs of 7.5~93.0% can influence battery impedance, and the results indicated that both charge transfer resistance and diffusion resistance increased dramatically when SOC is below 26.5%. Jiang et al. [
17] investigated the impedance of batteries with different life cycles (up to 520 cycles) and found that when the frequency is higher than 100 Hz, the impedance phase of different, aging batteries would be different to each other. Other than the influence of SOC and battery aging status on battery impedance, battery temperature also has a great impact on impedance. Andre et al. [
18] investigated the influence of temperature and SOC on LIBs and, through the cell impedance, the non-linear correlation between temperature and impedance was found. Raijmakers et al. [
19] studied the impedance characterization of a pristine nickel cobalt aluminum oxide battery (NCA) cell at different temperatures and found that ohmic resistance would decrease as the temperature increases and inductive characteristics barely change. Xiao et al. [
20] found that the battery impedance is sensitive to temperature because the ions’ mobility and permittivity are dependent on temperature.
During operation, heat generation caused by charge/discharge may contribute to the temperature rise, and undoubtedly, the temperature rise is inhomogeneous. Usually, the cell internal temperature will always be higher, and the temperature will decrease from inside to outside, which can be called the temperature gradient. It is a common phenomenon, but very few studies have investigated the influence of the temperature gradient. In references [
21,
22], a lumped-parameter electro-thermal model was proposed and applied. The parameters of the electrical model were modified according to the average battery temperature. However, the core temperature was obviously higher than the surface temperature, and the influence of the temperature gradient was ignored. The influence of the temperature gradient on an LIB’s impedance characterization has not been extensively studied either. Troxler et al. [
23] designed an artificial temperature gradient for the battery and analyzed the performance of LIBs. For a cell with similar average temperature, the larger the temperature gradient, the greater the behavior deviation from a cell with a uniform temperature will be. However, that study only focused on the charge transfer resistance. In reference [
11], the overall resistance of the battery with inhomogeneous temperature distribution was similar to the resistance of the battery with homogeneous temperature distribution, and the homogeneous temperature was the mean value of the inhomogeneous temperature.
The objective of this paper is to investigate battery impedance characterization considering the internal temperature gradient, and then prove that the model considering the temperature gradient will perform better. We will propose a novel discretization model and experiments with the design of artificial temperature gradients for the battery with the help of a heating plate. The remaining part of this paper is organized as follows:
Section 2 will introduce the battery impedance, and a discretization model, based on the non-linear impedance model, is proposed to study the influence of the battery internal temperature gradient on impedance. In
Section 3, several experiments are conducted to investigate the impact on battery impedance, and a heating plate is used to construct the artificial temperature gradient.
Section 4 will provide some experimental results and discuss the influencing factors on battery impedance. The verification results of the proposed discretization model are also shown in order to verify the model performance. Finally, some conclusions are drawn and summations provided in
Section 5.
3. Experiment
EIS experiments are undertaken to investigate the battery’s characteristics. A pouch-type LiFePO4 battery with a nominal capacity of 40 Ah is used in the experiment, and there are 51 electrode layers stacked inside the battery. The size of the battery is 133 mm × 16.1 mm × 184 mm.
Experiments with temperature gradients or not are designed and conducted in this section. In the experiment without a temperature gradient, usually little heat is generated in the EIS experiment, so the environmental chamber temperature is considered as the battery temperature. In the experiment with temperature gradients, a heating plate is used to construct an artificial temperature gradient in order to obtain an accurate and qualitative result. One side of the battery is covered with the heating plate, and the other side is exposed to the environment. Apparently, due to heat convection, the real temperature of the side coated with the heating plate may be lower than the heating plate temperature, and due to heat conduction, the real temperature of the side exposed to the environment may be higher than the environmental temperature. We use thermocouples to obtain the actual temperature and the temperature gradient of the battery. In addition to this, in order to confirm the assumption that the temperature gradient is linearly distributed, two thermocouples are placed inside the battery.
Figure 6 shows the position of the thermocouples inside the battery. T-type thermocouples are used in the experiments, and the diameter of the thermocouples is 0.28 mm. One thermocouple is placed in the halfway point within the battery, and another is situated at the quarter point within the battery. The battery cell under experiment is shown in
Figure 7a.
Figure 7 shows the EIS experimental setup. An electrochemical workstation, as shown in
Figure 7b, including the electrochemical interface and the frequency response analyzer, is used to generate the sinusoidal current and then analyze the frequency response of the cell in order to obtain the battery impedance. A power booster (
Figure 7c) can enlarge the sinusoidal current generated by the electrochemical workstation. A temperature measurement station (
Figure 7d), including the wireless temperature unit and the wireless logging station, is used to obtain the temperature and log the data. The battery cell is placed in an environment chamber, as shown in
Figure 7e.
3.1. Experiments without the Temperature Gradient
To acquire the characterization of battery impedance and the relationship between the parameters of the battery’s non-linear impedance model and temperature, the EIS experiment without the temperature gradient is conducted first.
Figure 8 shows the diagram of battery EIS experiment without the temperature gradient.
Table 1 shows the detailed EIS experiment steps without the temperature gradient. Enough rest time is necessary before the test to improve the accuracy of the experiment. A sinusoidal current with 2 A amplitude and a frequency ranging from 0.01 Hz to 10 kHz is loaded on the battery.
3.2. Experiments with Artificial Temperature Gradients
As mentioned above, a heating plate is used to construct the temperature gradient.
Figure 9 shows the diagram of the battery EIS experiment with temperature gradients. One side of the battery cell is covered with a heating plate. Besides the thermocouples inside the battery, there are thermocouples stuck to the surface of the battery. Although the heating plate temperature and the chamber temperature can be set to specific values, the actual temperature of the battery surface measured by the thermocouples will be used to analyze the temperature gradient. The battery maintains 50% SOC during the experiment with an artificial temperature gradient in order to focus on the correlation of the temperature gradient.
Table 2 shows the detailed EIS experiment steps with a temperature gradient. Enough rest time is also necessary before the test to improve the accuracy of the experiment.
4. Results and Discussion
4.1. Battery EIS Characterization without the Temperature Gradient
Figure 10 shows the EIS experiment results of the battery with a uniform temperature.
Figure 10a shows the EIS of the battery at different temperatures, and from this we can see that, with the increase of temperature, the impedance becomes smaller when the frequency is low. While in the high-frequency region, the EIS curves with different temperatures coincide with each other and this is clearer in the enlarged part of
Figure 10a. This indicates that the temperature mainly influences the battery diffusion and charge transfer reaction processes, while the battery’s inductance and ohmic resistance are not affected by the temperature.
Figure 10b shows the EIS of a battery with the same temperature but different SOCs. From this picture, we can see that all the curves are close to each other in the middle and high frequencies, and in the low-frequency region, the most obvious difference is the slope. With the increase of the SOC, the slope of the straight line becomes larger. This indicates that SOC mainly impacts the battery diffusion process.
The EIS results of the battery with uniform temperature are shown in
Figure 11 with the frequency as a variable;
Figure 11a,b shows the impedance magnitude and phase of the battery under different temperatures, while
Figure 11c,d is the impedance magnitude and phase of the battery at different SOCs.
Figure 11a shows that the battery impedance will increase as the temperature decreases, and the battery impedance at −20 °C is almost 10 times the impedance at 20 °C at the frequency of 0.01 Hz to 0.1 Hz. The impedance phase curves shown in
Figure 11b have a great discrimination between 1 Hz and 10 Hz, and that is why the impedance phase can be used to estimate battery temperature [
12]. In a similar way to
Figure 10b, the impedance magnitude curves in
Figure 11c and the phase curves in
Figure 11d are close to each other in all frequency ranges, and it can be concluded that the SOC has little impact on the battery impedance.
4.2. Battery EIS Characterization with Temperature Gradients
Temperature distribution inside the battery is analyzed before obtaining the EIS experiment results of the battery with the artificial temperature gradients.
Figure 12 shows the temperature distribution inside the battery with the artificial gradients. Due to the battery’s thickness and heat transfer, the actual temperature gradient may be different to the temperature gradient that is set. Six groups of data are shown in this figure, including three groups with a smaller set temperature gradient (10 °C), and three groups with a larger set temperature gradient (35 °C). To get a more intuitive view of these data, we give them a linear fit, and the fitting results are also shown in the figure. From this, we can see that the fitting results are excellent, and the R-Square of all fitting results are greater than 0.98. It can be concluded that the temperature inside the battery is linearly distributed, and the assumption proposed in
Section 2 is reasonable.
As mentioned above, the temperature has a great impact on battery impedance. The next point to study is the influence of the temperature gradient on the impedance.
Figure 13 shows four sets of data, including one experiment at a 10 °C uniform temperature and three experiments with different temperature gradients. In the experiment at a 10 °C uniform temperature, the battery temperature gradient is considered as zero. Although the other three experiments are implemented in different environmental temperature, they have one thing in common, the actual heating plate temperature is about 10 °C (10.59 °C, 10.28 °C, 9.76 °C), so these experiments are considered to have a similar maximum temperature. In the figure, from the black line to the magenta line, the temperature gradient becomes larger.
Figure 13a shows that battery impedance will increase as the temperature gradient becomes greater from 0.01 Hz to 10 Hz, and
Figure 13b shows that the battery impedance phase will decrease as the temperature gradient becomes greater from 1 Hz to 100 Hz. It may be concluded that if the maximum battery temperature is fixed, the smaller the minimum temperature, the greater the battery impedance magnitude.
There is another point worth noting in
Figure 12: the magenta line and the olive line intersect at the half of the battery. According to the assumption that the temperature is linearly distributed, the temperature at the battery’s center can be considered as the average battery temperature. A situation whereby several experiments under different temperature gradients but the same battery average temperature demonstrate similar impedance characterization should be verified and discussed.
We select some experiments with approximate average temperature, including the lower average temperature (1 °C) and the higher average temperature (21 °C). Besides these, we also supplement the experiments with those of 1 °C and 21 °C uniform temperature.
Figure 14 shows the experimental EIS of a battery with approximate average temperature, and in
Figure 14b,d, the impedance phase of a battery with different temperature gradients, but in both cases the similar average temperatures almost coincide with each other. In
Figure 14a,c, the impedance magnitude of each temperature gradient differs from each other in the low frequency range and is unrelated to the temperature. It may be concluded that batteries with similar average temperature possess similar impedance characterizations.
Some summaries and conclusions about the battery impedance characterization can be drawn according to the above experimental results:
- (1)
Temperature mainly influences the diffusion and charge transfer reaction processes, while the battery inductance and the ohmic resistance are not affected. SOC only has an impact on the battery diffusion process. In general, the temperature has a greater impact on battery impedance than SOC.
- (2)
For batteries with the same maximum temperature, battery impedance will increase as the temperature gradient becomes greater, especially in the low-frequency region, and the battery impedance phase will decrease as the temperature gradient become greater in the middle-frequency region.
- (3)
In our experiment with artificial temperature gradients, we find that the temperature distribution can be approximately considered to be linear, so the central temperature is equal to the average temperature. The experiments with similar average temperature possess similar impedance characterization, and this is consistent with the idea of lumped parameter model.
4.3. Determination of the Parameters of a Battery Non-Linear Impedance Model
The results of experiments without the temperature gradient are used to fit the equivalent circuit model, and one of the fitting results is shown in
Figure 15. This indicates that the experimental data and the fitting data almost coincide at all the experimental frequencies, and
Figure 15a shows that in the high-frequency region the EIS is not a straight line perpendicular to the abscissa axis; so, maybe a resistance parallel with inductance
L can improve the fitting. Overall, the fitting results are acceptable and can be used in the discretization model.
Figure 16 shows the fitting parameters of the battery model, and
Ro and
R1 decrease exponentially with temperature increase, while
Q1, and
Q2 increase exponentially when the temperature increases.
We employ a polynomial fit to get the relationship between the battery impedance model parameters and temperature, and the formula can be listed as:
where
t is the temperature in centigrade, and
L,
Ro,
R1,
Q1,
Q2,
α and
β are the model parameters.
The above work has enabled the relationship between the battery impedance model parameters and the temperature to be obtained, and the relationship will be used in the battery discretization model.
4.4. Verification of Battery Impedance Model Considering the Temperature Gradient
In most cases, the battery’s surface temperature or average temperature is used to modify the battery’s impedance. However, the above analysis makes it obvious that the internal temperature gradient has a great influence on battery impedance.
Figure 17 shows the comparison of the experimental results and lumped model impedance modified by the surface temperature. The parameters of the lumped model are obtained according to Equation (5) and battery surface temperature. Plainly, there are large deviations between the experimental results and the lumped model impedance modified by the surface temperature, no matter what the maximum surface or minimum surface temperatures are.
According to the discretization model shown in Equation (4) and the relationship between the impedance model parameters and temperature in Equation (5), the impedance of a battery with artificial temperature gradients can be obtained. Three experiments with a higher temperature gradient (the actual temperature difference between battery surfaces is greater than 20 °C) are chosen to verify the battery discretization model.
Figure 18 shows the verification of the battery discretization model. The discretization model considering the temperature gradient and the lumped model considering only the average temperature are also compared.
Overall, the proposed battery discretization model is closer to the actual experimental data than the lumped model. From
Figure 18b,e,h, we can see that both the proposed model and the lumped model deviate from the actual experimental data from 0.01 Hz to 10 Hz. This may be because that the relationship between temperature and the fractional order of
CPE1 and
CPE2 may not be linear. From
Figure 18c,f,i, in most cases, the proposed model’s impedance phase is closer to the experimental data, which proves the accuracy of the proposed model.
Figure 19 shows the errors of the discretization model and the lumped model. Obviously, the magnitude error of the discretization model is smaller in the whole frequency range, while the phase errors in
Figure 19d–f show that the lumped model performs a little better than the discretization model in the low-frequency region. Apparently, the lower the temperature, the greater the magnitude error. In
Figure 19a–c, the experiment with the lower average temperature has a greater error in the low frequency range. In the high-frequency region, the relationship between the phase error and the frequency is approximately linear and independent of temperature. This may be due to the battery impedance model’s inherent defect that a single inductor may not represent the real characterization of battery impedance in the high-frequency region. In the whole frequency range, the maximum impedance magnitude error of the battery discretization model is 1.46 mOhm, and the maximum impedance phase error is 6.96°.
Table 3 compares the statistical data of the discretization model and the average temperature model. Obviously, the accuracy of the discretization model is higher, especially the impedance magnitude of the proposed model. Among the three experiments with different temperature gradients, the root mean squared error (RMSE) of the proposed model is at least 15% smaller than that of the average temperature model. For the model impedance phase, the proposed model is slightly superior.
Through the above analysis, the performance of the proposed model is proved to be accurate and effective. Compared with the lumped model with the average temperature, the proposed model is identical to the actual battery impedance. With the help of the discretization model, we can obtain a more accurate characterization of the battery.
4.5. Limitations and Future Work
In the proposed discretization model, the temperature distribution in the direction parallel to the battery taps is ignored. In fact, the temperature distribution is inhomogeneous in this direction. The temperature of the part close to the battery taps is higher. Moreover, this model may not be suitable for cylindrical-type batteries because the volume of inner and outer layers is different and the calculation of each layer impedance will become complex.
The proposed discretization model is based on the battery non-linear impedance model. It is worth noting that this model is generic for batteries with similar characterizations like the impedance spectroscopy shown in
Figure 1. For other lithium-ion batteries with different battery chemistries, the battery non-linear impedance model may be different and should be determined first. In that case, the discretization method can be applied to batteries with various chemical materials. Verification of the adaptability of this discretization method to other battery chemistries shall be conducted in future work.
Considering the thickness of the battery used in this paper, the temperature distribution is assumed to be linear. Unlike the artificial temperature gradient in this paper, the core and surface temperatures will construct the temperature gradient in actual applications of batteries. Future work will focus on the temperature distribution of the battery with a larger size. Based on an electrothermal model, the core temperature will be estimated online, and then the model parameters will be modified according to the proposed discretization model. The modified model parameters will eventually improve the precision of the battery electrothermal model.