1. Introduction
Lithium iron phosphate (LiFePO
4) batteries have become popular for renewable energy storage devices, electric vehicles [
1,
2,
3] and smart grids [
4,
5,
6]. However, the batteries are quite vulnerable to temperature variation and unforeseen operating conditions such as being overly charged or discharged, which affects its performance and reduces its lifespan in an electric vehicle. There have been many efforts in recent years to enhance the accuracy of the
SOC estimation in a battery management system (BMS) [
7]. A detailed literature review can be found in the following references [
8,
9,
10,
11]. The common Coulomb counting or Ampere-hour method [
12] uses the current reading of the battery over the operating period to calculate
SOC values.
Hence, an accurate and precise physical modeling on electrochemical batteries for
SOC estimation [
13] is needed. The electrochemical model [
14] is one of the first such models using the Lithium-ion, electrolyte-phase concentration and reaction current density in partial differential equations (PDEs). The electrochemical model can provide accurate estimations of the cell behavior over a range of operating conditions at the expense of high computational complexity due to PDEs that are unsuitable for online application. To circumvent the problem, different model reduction techniques were applied. The most common reduced-order model can be achieved by neglecting the electrolyte-phase potential and assuming a constant electrolyte-phase concentration [
15]. Such assumptions are allowed for low charging rates. Hence, the use of model [
16] that preserves some of the physical insights to estimate the
SOC at higher discharge rates. A similar model can be seen in Reference [
17] for higher charge and discharge rates. The reduced electrochemical model of composite electrode Lithium-ion battery [
18] was used for dual-nonlinear observers that could estimate five different cells’
SOC over time. An alternative approach, such as the equivalent circuit model (ECM) that models the internal battery characteristics using resistors and capacitors are commonly applied for
SOC estimation due to their low complexity and simplicity. Some of the early works on
SOC estimation include an adaptive filter such as an extended Kalman filter [
19] on a 1-RC Lithium-ion battery cell. Although the aging effect was included with ±5% of the values estimated by Ah counting, this publication did not consider the ambient temperature as varying model parameter for estimation. Extension to 2-RC Lithium-ion battery cell via EKF [
20] was studied. But the results did not consider the influence of ambient temperature, and the aging [
19] in the
SOC estimation.
As the ambient temperature affects the
SOC estimation, the 1D (one-dimensional) electrochemical with the thermal model was used to understand the heat distribution of surface and internal temperature on a pouch-type lithium-ion battery cell [
21] using thermodynamics and kinetics. The effects on the
SOC estimation was not discussed. In addition, the ECM was used to simulate the electrochemical, thermal model [
22] at an optimal cycling period. The results help to identify the battery cell parameters. However, the influence of ambient temperature and its impact on
SOC estimation were not studied. As a result, a fully coupled electro-thermal model [
23] with a fast charging strategy was modeled as a linear-time-varying model predictive control problem, while the unmeasurable battery internal states such as the
SOC and core temperature were estimated via a nonlinear observer. But the works were surrounded on control to balance the time and core temperature increase instead of
SOC identification.
More works were then performed to validate the lithium-ion battery at a different temperature such as 25 °C [
24] and a low temperature of −10 °C [
25] via a pseudo-two dimensional electrochemical thermal model [
26]. Although the paper [
27] shows that the macroscopic models are strongly controlled by the battery operating temperature conditions, the works were mainly on the proposed model that could affect the battery operating or internal temperature. The results of pulse-relaxation discharge in the electrochemical thermal model [
28] describe the dynamic change of Lithium-ion concentration distribution in phases that were useful for analyzing the polarization of the battery. Another pseudo-two dimensional electrochemical thermal model [
29] showing the battery pack application was then developed for the electric vehicle battery pack and thermal management system to maintain the temperature uniformity and decrease the maximum temperature of the cells. The
SOC was not undertaken to further examine the impact of
SOC estimation. In addition, other than modeling and analyzing the relationship between the temperature (−20 °C to 60 °C) and battery model parameters, the analyzed results did not show how the electrochemical-thermal battery model can be adaptively tuned for the different ambient temperature in the case of the multiple-cell battery. The research was further performed to model the electro-thermal behavior [
30,
31] of the battery of different size of cylindrical cells with the application of the battery pack on an actual electric vehicle under different drive tests. An electro-thermal model [
32] of lithium-ion batteries based on the ECM was proposed. The model parameters of the cell were optimized and validated for a wide range of temperatures and
SOC that led to the use of ECM again. However, the tuning of the battery model at different ambient temperature was not mentioned.
A few nonlinear observer design methods were applied to derive the ECM-based nonlinear
SOC estimators such as sliding mode observer [
33], nonlinear observer [
34] and Lyapunov-based observer methods [
35]. The ambient temperature was not considered during the
SOC estimation. Nevertheless, an online parameter identification [
36] was recently used to estimate the battery parameters and nonlinear Kalman filter for
SOC estimation at different ambient temperatures. Another approach considering the battery temperature using particle swarm optimization and Gauss-Newton algorithm [
37] to estimate the battery parameters was used. The temperature considered was either the battery’s internal or surface temperature.
An alternate non-model based approach including neural networks [
38,
39,
40], fuzzy logic [
41], neural network-fuzzy [
42], support vector machine (SVM) [
43,
44] and extreme learning machine (ELM) [
45,
46,
47,
48,
49,
50,
51] methods was developed to predict the
SOC. These machine learning methods require sufficient large dataset and computation time for the training and validating the
SOC value. They are not quite suitable for the operating battery power system as significant training time is required (except for ELM that uses regularized least squares to compute faster than the conventional quadratic programming approach used in gradient method). A battery model [
52,
53,
54,
55,
56,
57] is briefly discussed that could affect the
SOC estimation. In addition, the ECM is likely to have higher precision in practice due to its simplicity in the parameter identification. Hence, an adaptable ECM based Lithium-ion battery with consideration of the ambient temperature variation for more accurate
SOC estimation is required. In contrast to works done in battery parameter estimation [
58,
59,
60,
61,
62], the proposed method is based on the temperature-dependent cell model.
In summary, the contributions are as follows. A unified lithium-ion battery model that includes the ambient temperature, cell temperature, thermal aging effect on capacity, hysteresis voltage and its impact on state-of-charge (SOC)-open-circuit voltage (OCV) relationship, hysteresis voltage dynamics and terminal voltage is proposed and validated by an actual experimental test.
The paper is organized as follows: The new nonlinear battery model including ambient temperature is derived in a few subsections starting from
Section 2,
Section 3,
Section 4 and
Section 5.
Section 6 represents the methodology for
SOC estimation using the extended Kalman filter.
Section 7 provides the conclusions.
2. Cell Nonlinear Temperature-Dependent State Model
An ambient and cell temperature-dependent
equivalent circuit model (ECM) in
Figure 1a is used. As shown in the circuit model, 1RC is used instead of 2RC. The comparative curves between the 1RC and 2RC can be seen in
Figure 1b. The root mean square error of 1RC and 2RC is around
V (i.e., quite small). Hence, 1RC is used for simplicity as fewer parameters are required for computing terminal voltage and
SOC value.
The ambient
temperature refers to the ambient temperature where the battery cell operates while the cell temperature
is the surface temperature of the cell. The terms are incorporated into
to represent the temperature dependency in obtaining the terminal voltage (instead of an external disturbance for model validation and robustness test). In addition, the proposed battery cell model provides the current-voltage behavior compensated by the hysteresis effect that gives a balance between simplicity and accuracy in modeling the battery cell model. Note that the experimental data used for the subsequent comparison was obtained from the cylindrical lithium iron phosphate battery cells (ANR26650M1-B) in
Table 1.
From the ECM model in
Figure 1a, the equations corresponding to the ECM can be expressed as follows.
where the internal resistance is represented by
in Ω,
and
are the polarization resistance (in Ω) and polarization capacitance (in F) as a function of
(in °C) and state-of-charge,
. The cell capacity
is a function of ambient temperature (
). The voltage across the single RC circuit is indicated as
in V. The hysteresis voltage dynamics [
63] at equilibrium state is denoted by
that is a function of
and current
in A. The function
for
greater or equal to 0 and −1 otherwise. The constant
has a value of 0.0245. The terminal or the output voltage of the battery cell is denoted by
. The temperature dependency term
is included in the terminal voltage equation to compensate for the temperature variation on the output voltage. The open-circuit voltage
is a function of
and cell temperature
(in °C). The following will explain how the parameters used in (1)–(4) are modelled.
3. Cell SOC-OCV Relationship
The
SOC in (1) was obtained by the Ampere-counting method. However, the cell capacity
can be a function of ambient temperature (
) as shown in
Figure 2 where the capacity varies at different ambient temperature. The static capacity curve is obtained by curve-fitting the experimental data. The capacity model
used in (1) can be written as follows.
The map** of open circuit voltage (OCV) to
SOC is represented by
Uoc as follows.
where
and
is the cell temperature in °C. The plot of open-circuit voltage,
to
SOC relationship in (6) as compared to the experimental OCV-SOC obtained from pulse discharge test (PDT) in [
51,
64,
65,
66] can be seen in
Figure 3. The PDT conducted in
Figure 3 generated a series of the discharge current of 2 A for decreasing
SOC values in
Figure 4. In the test, the battery was fully charged with constant current-constant voltage (CCCV) mode under 25 °C ambient temperature to 3.6 V. Then battery cell was then discharged at a pulse current 2 A with 450 s discharging time and 2700 s relaxation time until the terminal voltage reaches the cut-off voltage of 2.5 V as seen in
Figure 5. The OCV was obtained when the battery’s terminal voltage reached equilibrium state during the relaxation period. As seen in
Figure 6, the values of the
SOC were then recorded. The cycle was repeated at other ambient temperatures. The battery cell was tested under specific ambient temperature (i.e., 5 °C, 15 °C, 25 °C, 35 °C, and 45 °C) in a temperature chamber. The reason for using the selected temperature is due to the limitation of the temperature chamber at the time of the testing, although the battery cell can go up to 55 °C (and in some application such as an i3 battery for BMW that could go beyond 70 °C) and reached a minimum of 5 °C only. In this study, the maximum limit of the operating ambient temperature is set to 45 °C with the use of heat absorbing material such as phase change material (PCM). The PCM can allow higher operating ambient temperature if needed.
As shown in
Figure 6, the simulation of
SOC using (1) is compared with the experimental data obtained from Ah counting method (without the ambient temperature in the static capacity). The results show slight differences in
SOC estimation. The simulation result using (1) gives a higher value as compared to the Ah counting method. It shows that it is essential to include an ambient temperature in
SOC estimation. On the other hand, the OCV-SOC relationship proposed in (6) exhibits a close match to the experimental data as seen in
Figure 7a. The constants used in the equation were tuned iteratively to match the experimental OCV-SOC data.
The OCV-SOC at various ambient temperature can be simulated in MATLAB/Simulink as shown in
Figure 7b. It is used to compare the simulated results with the experimental data as shown in
Figure 7c. The root mean square errors (RMSEs) of OCV at different ambient temperature as compared to the proposed model can be seen in
Figure 7d. The maximum RMSE of OCV is around 0.33 V at an ambient temperature of 5 °C followed by 15 °C with RMSE of 0.28 V.
4. Cell Temperature Model
In addition to the cell temperature affecting the OCV-SOC relationship, the ambient temperature [
67,
68] also affects the accuracy of terminal voltage and subsequent
SOC estimation. To account for the influence, the term
is included to compensate for the differences on the terminal voltage equation.
where
,
is obtained from the general energy balance for battery thermal model proposed by Bernardi et al. [
69] with the total heat generated
in the battery cell as defined by
where
is the mass of the battery cell in kg,
is the specific heat capacity of the cell in J/kgK,
is the convective heat transfer coefficient,
,
is the emissivity of heat term,
is the Stefan-Boltzmann constant term in J·s
−1m
−2K
4,
is the reversible heat term or entropy coefficient and
is the contact resistance in
. The terminal voltage
and open-circuit voltage
can be seen in (4) and (6), respectively. The values of the parameters related to the battery cell can be found in
Table 1.
It is known that RC values depend on the
SOC and
. It is essential to include the dependent terms
SOC and
in (2) and (4) as seen in the following
Figure 8,
Figure 9 and
Figure 10. The regression equations
have a function of
SOC. As compared with the ambient temperature counterpart, the behaviors of
at different ambient temperature are difficult to model accurately. Hence, the dynamics of
use in (2) and (4) are embedded into a series of lookup tables such that the values at specific
SOC and
can be obtained. The relationship between
and
SOC can be expressed as follows. For example, the equations for
in (9) to (11) were derived by curve-fitting the experimental data at ambient temperature of 25 °C.
where
.
The terminal voltage is also affected by the changes in
instead of
. The schematic of the experimental setup that utilizes infrared thermal imaging was used to determine the cell temperature as seen in
Figure 11. The laboratory environment was first kept approximately 25 °C and discharged to a programmable electronic load. A thermal infrared camera was used to capture the temperature distribution of the battery cell which was wrapped in non-glossy black tape with a known emissivity of 0.95. The battery cell was discharged under a constant current rate (2 A) at the same ambient laboratory temperature. The experiment was performed on three different battery cells and the average cell temperature was taken. The stop condition for the discharge was set to a cut-off voltage of 2.5 V to prevent damage to the cell. For clarity, only one of the image captured can be seen in
Figure 12. The surface temperature can reach around 29 °C as seen in the thermal image. The same test was then repeated at other ambient temperatures.
The time response,
obtained from the experimental is compared with the simulated result by Bernardi et al. [
69]. It can be observed that there exist some errors. The errors are found to be around 3 °C at the steady-state. The high time constant resulted in the sluggish response as seen in
Figure 13. The same phenomena can be observed at other ambient temperature. As time increases, the cell temperature will increase due to the continuous discharging the cell at 2 A until it reaches the cut-off voltage of 2.5 V. The entire simulation block diagram by MATLAB/Simulink can be seen in
Figure 14. The block diagram was used to simulate the OCV-SOC relationship, cell temperature, terminal voltage and
SOC estimate of the battery cell at different ambient temperature. The input load current can be changed.
6. SOC Estimation by Extended Kalman Filter Method
To facilitate the
SOC estimation process by EKF, the above equations from (1) to (3) have to be discretized [
70]. Equation (6) and (7) will be included in the output equation in (4).
First,
in (1) can be discretized as follows.
Then,
in (2) can be discretized via the following steps.
where
is the current time,
is initial time,
and
.
The current remains unchanged within one sampling time T that is
equals to a constant. Defining
,
and solving the integration with lower and upper limit of 0 and T, respectively gives
where
,
and
is the sampling period of 1 s.
Next, the solution for
in (3) can be written as follows.
Equation (16) can be discretized as
The proposed temperature-dependent battery cell model that includes the ambient temperature, cell temperature, hysteresis voltage dynamics, the thermal aging effect on capacity and terminal voltage can be written as a nonlinear system equation.
where
,
is the time index,
is the current at
,
is the input vector,
is the process noises and
is the measurement noise. The process and measurement noise satisfy the following equations.
where
is Kronecker delta,
and
are symmetric positive definite matrix.
The matrices used in (18) are defined as:
where
is the sampling period of 1 s. Note that the function of
and
in
are omitted for clarity.
Based on the nonlinear battery cell model obtained in the previous sections, an extended Kalman filter (EKF) can be designed to estimate the
SOC of the cell. The EKF remains quite appealing for most practitioners even with the inherent linearization error and noises [
63]. The detrimental effects due to the linearization-induced error on the
SOC estimation are not significant as the battery cell model does not have strong nonlinearity. The state equation is quite linear in the state and the nonlinearity in the output equation has a bounded gradient [
63]. The model mismatch due to linearization can be resolved by setting the initial
SOC value close to the true value obtained from OCV that can give a reasonable estimate of
SOC value. The EKF design [
19] for the state estimation is straightforward and is thus omitted in this paper. The initial conditions of the battery cell model are taken as
and the EKF estimator is
. The tuning parameters of the EKF are taken as:
where
is
identity matrix. The matrix
,
and
are symmetric positive definite. The EKF is quite robust to changes of the tuning parameters, and the process of tuning them is straightforward.
As shown in
Figure 20, with the EKF and proposed cell model use the same values of model parameters, the
SOC values converge to the reference values obtained by Ah counting method. The maximum RMSE of the
SOC estimate is around 0.09. The
SOC converged to the actual value at time increases. When compared at different ambient temperature,
Figure 21 shows the response of the
SOC estimate when 5 °C, 15 °C, 35 °C, and 45 °C. It can be observed at lower ambient temperature; the
SOC value is lowered than the ambient temperature of 25 °C and above.
The battery cell model that includes the ambient temperature, cell temperature, hysteresis effect and thermal aging effect on capacity, OCV-SOC and terminal voltage was simulated. As ambient temperature varies, the response of
SOC deviates from the reference or experimental data (by Ah counting method). However, the proposed cell model with temperature-dependent terms is capable of estimating
SOC despite the variation in ambient temperature as seen in
Figure 22 and
Figure 23.
Similar phenomena can be seen in
Figure 24 and
Figure 25 at different ambient temperatures for the terminal voltage. Without the temperature-dependent terms in the battery cell model, the terminal voltage is higher than the actual terminal voltage (that does not represent the characteristic of the battery cell where the terminal voltage can exhibit a lower value).
The results of the
SOC estimation and terminal voltage using the proposed model with EKF based
SOC estimation method (with and without considering the ambient temperature) is shown in
Table 3. The model without the ambient temperature-dependent term in static capacity,
and RC parameters (
, gives a less RMSE as compared to the non-temperature-dependent model. Both model with and without ambient temperature term used the same
in (6). The non-temperature-dependent model adopted the fixed RC parameters [
51,
64,
65,
66]:
. The similar phenomena can be seen in the terminal voltage where lower RSME can be observed (except for 25 °C and 35 °C). It may be due to the modelling was mostly performed at 25 °C, and thus improvement is not prominent. Nevertheless, the increased in RMSE is not more than 5%. In addition, there exists at least 57% improvement in
SOC estimation at a lower temperature. For the case of terminal voltage, the improvement is around 56% at a lower temperature of 5 °C. There is no improvement at ambient temperature of 25 °C.
In summary, the proposed cell model can estimate the
SOC and terminal voltage under the different and aging thermal effect on the battery cell at a low (i.e., 5 °C and 15 °C) and high ambient temperature (45 °C). Note that the results depend on the type of battery cell used. However, a similar approach can be applied in various applications [
71,
72].