Nonlinear Temperature-Dependent State Model of Cylindrical LiFePO 4 Battery for Open-Circuit Voltage, Terminal Voltage and State-of-Charge Estimation with Extended Kalman Filter

: Ambient temperature affects the performance of a battery power system and its accuracy in state-of-charge ( SOC ) estimation for electric vehicles and smart grid systems. This paper proposes a battery model that considered ambient temperature, cell temperature, hysteresis voltage and thermal aging on capacity due to multiple charging and discharging. The SOC is then estimated using an extended Kalman ﬁlter. Several forms of validation were tested on an actual cell battery under speciﬁc ambient temperatures to verify the battery cell model, terminal voltage and SOC estimation performance. The SOC estimation results show an improvement in root-mean-squared error as compared to Extended Kalman Filter (EKF) without considering the temperature dependency. The proposed battery temperature-dependent model gave a smaller root-mean square error in SOC and terminal voltage at 5 ◦ C, 15 ◦ C and 45 ◦ C.


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 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 3.2775 × 10 −4 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 (T amb ) temperature refers to the ambient temperature where the battery cell operates while the cell temperature (T cell ) is the surface temperature of the cell. The terms are incorporated into U temp (T cell , T amb ) 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. From the ECM model in Figure 1a, the equations corresponding to the ECM can be expressed as follows. S . OC = − I 3600C n (T amb ) (1) .
. U hs = −β|I|[0.0755(1 − SOC)sign(I) + U hs ] Energies 2018, 11, 2467 5 of 28 U t = U oc (SOC, T cell ) − U 1 − R 0 (SOC, T amb )I + U hs + U temp (T cell , T amb ) (4) where the internal resistance is represented by R 0 in Ω, R 1 and C 1 are the polarization resistance (in Ω) and polarization capacitance (in F) as a function of T amb (in • C) and state-of-charge, SOC ∈ (0, 1). The cell capacity C n is a function of ambient temperature (T amb ). The voltage across the single RC circuit is indicated as U 1 in V. The hysteresis voltage dynamics [63] at equilibrium state is denoted by U hs that is a function of SOC and current I in A. The function sign(I) = 1 for I 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 U t . The temperature dependency term U temp (T) is included in the terminal voltage equation to compensate for the temperature variation on the output voltage. The open-circuit voltage U oc is a function of SOC and cell temperature T cell (in • C). The following will explain how the parameters used in (1)-(4) are modelled.

Cell SOC-OCV Relationship
The SOC in (1) was obtained by the Ampere-counting method. However, the cell capacity C n can be a function of ambient temperature (T amb ) 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 C n (T amb ) used in (1) can be written as follows.
where the internal resistance is represented by 0 in Ω, 1 and 1 are the polarization resistance (in Ω) and polarization capacitance (in F) as a function of (in °C ) and state-of-charge, (0,1) SOC  . The cell capacity is a function of ambient temperature ( ). The voltage across the single RC circuit is indicated as 1 in V. The hysteresis voltage dynamics [63] at equilibrium state is denoted by ℎ that is a function of SOC and current in A. The function sign( ) = 1 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 SOC and cell temperature (in °C ). The following will explain how the parameters used in (1)-(4) are modelled.

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.   [51,[64][65][66] can be seen in Figure 3. The PDT conducted in Figure 3 generated a series of The mapping of open circuit voltage (OCV) to SOC is represented by U oc as follows.
where SOC ∈ (0, 1) and T cell is the cell temperature in • C. The plot of open-circuit voltage, U oc 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 Energies 2018, 11, 2467 6 of 28 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.
Energies 2018, x, x FOR PEER REVIEW 6 of 29 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.    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.     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.

DC
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.  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. 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.

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. 6 ( , ) 0.0492 6.8 10 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 m 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 −1 m −2 K 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.

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 U temp (T cell , T amb ) is included to compensate for the differences on the terminal voltage equation.
where ∆T = T cell − T amb , T cell is obtained from the general energy balance for battery thermal model proposed by Bernardi et al. [69] with the total heat generated (Q gen ) in the battery cell as defined by

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. 6 ( , ) 0.0492 6.8 10 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 m 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 −1 m −2 K 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.  (8) where m is the mass of the battery cell in kg, c p is the specific heat capacity of the cell in J/kgK, h is the convective heat transfer coefficient, ∆T = T cell − T amb , ε is the emissivity of heat term, σ is the Stefan-Boltzmann constant term in J·s −1 m −2 K 4 , dU oc /dT cell is the reversible heat term or entropy coefficient and R c is the contact resistance in Ω. The terminal voltage U t and open-circuit voltage U oc can be seen in (4) and (6), respectively. The values of the parameters related to the battery cell can be found in Table 1.  [63] It is known that RC values depend on the SOC and T amb . It is essential to include the dependent terms SOC and T amb in (2) and (4) as seen in the following Figures 8-10. The regression equations R 0 , R 1 , C 1 have a function of SOC. As compared with the ambient temperature counterpart, the behaviors of R 0 , R 1 , C 1 at different ambient temperature are difficult to model accurately. Hence, the dynamics of R 0 , R 1 , C 1 use in (2) and (4) are embedded into a series of lookup tables such that the values at specific SOC and T amb can be obtained. The relationship between R 0 , R 1 , C 1 and SOC can be expressed as follows. For example, the equations for R 0 , R 1 , C 1 in (9) to (11) were derived by curve-fitting the experimental data at ambient temperature of 25 • C.
where SOC ∈ (0, 1).  [63] 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 Figures 8-10. The regression equations 0 , 1 , 1 have a function of SOC. As compared with the ambient temperature counterpart, the behaviors of 0 , 1 , 1 at different ambient temperature are difficult to model accurately. Hence, the dynamics of 0 , 1 , 1 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 0 , 1 , 1 and SOC can be expressed as follows. For example, the equations for 0 , 1 , 1 in (9) to (11) were derived by curvefitting the experimental data at ambient temperature of 25 °C. 3 where ∈ (0,1).
(a) 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 terminal voltage is also affected by the changes in T cell instead of T amb . 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 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.  The time response, T cell 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. 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. Figure 13. Cell temperature at different ambient temperature (for constant 2 A load current) for 3000 s. Figure 13. Cell temperature at different ambient temperature (for constant 2 A load current) for 3000 s.

Cell Terminal Voltage Model
With U oc (SOC, T cell ) and U temp (T cell , T amb ) obtained, the terminal voltage, U t of the battery cell in (4) can be numerically computed using the simulation block diagram as shown in Figure 14. As seen in Figures 15-19, the terminal voltage at different ambient temperatures such as 5 • C, 15 • C, 25 • C, 35 • C, and 45 • C is plotted. It can be observed that the simulated terminal voltage matches experimental data near to the end of the discharging cycle or lower SOC value. The root means square errors (RMSEs) at each ambient temperature are tabulated in Table 2. The maximum error in the terminal voltage as compared with the actual data is around 0.0663 V. Note that the RMSE can be further improved by tuning the parameters in U oc and U temp in (6) and (7), respectively.

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).
Then, ̇1 ( ) in (2) can be discretized via the following steps.
where t is the current time, 0 The current remains unchanged within one sampling time T that is ( ) ( T) T tk     and solving the integration with lower and upper limit of 0 and T, respectively gives where 11 ( 1) [ and t  is the sampling period of 1 s.
Next, the solution for () hs Ut in (3) can be written as follows.

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).

SOC(t) in
Then, (2) can be discretized via the following steps.
where t is the current time, t 0 is initial time, t 0 = kT and t = (k + 1)T, k = 0, 1, 2, 3 . . . n. The current remains unchanged within one sampling time T that is I(τ) = I(KT) equals to a constant. Defining dt = −dτ, t = (k + 1)T − τ and solving the integration with lower and upper limit of 0 and T, respectively gives where U 1 (k + 1) = U 1 [(k + 1)T], I(k) = I(kT) and ∆t is the sampling period of 1 s.
Next, the solution for U hs (t) 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 X k = SOC k U 1 (k) U hs (k) T , k is the time index, I(k) is the current at k, u k = I k 1 is the input vector, w k is the process noises and v k is the measurement noise. The process and measurement noise satisfy the following equations.
where δ kj = 0 if k = j 1 if k = j is Kronecker delta, Q(k) and R(k) are symmetric positive definite matrix.
The matrices used in (18) are defined as: where ∆t is the sampling period of 1 s. Note that the function of T amb and SOC in R 0 , R 1 , C 1 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 X 0 = 0.9 0.12 0.002 T and the EKF estimator isX(0) = 0.1X(0). The tuning parameters of the EKF are taken as: Q = I 3 , R = 0.0005 and P 0 = 0.1I 3 where I 3 is 3 × 3 identity matrix. The matrix Q, P and R 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.
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 (0) = [0.9 0.12 0.002] T and the EKF estimator is ̂( 0) = 0.1 (0). The tuning parameters of the EKF are taken as: = , = 0.0005 and = 0.1 3 where 3 is 3 × 3 identity matrix. The matrix Q , P and R 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 Figures 22 and 23.  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 Figures 22 and 23. 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 Figures 22 and 23.        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 ( 0 , 1 , 1 ), 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 nontemperature-dependent model adopted the fixed RC parameters [51,[64][65][66]: 0 = 0.0823, 1 = 0.0042, 1 = 52.2. 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].

Conclusions
The effects of the ambient temperature, cell temperature, hysteresis voltage and thermal aging on capacity due to multiple charging and discharging was considered in the battery cell model of the cylindrical lithium iron phosphate battery cells (ANR26650M1-B). The temperature-dependent 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, U temp and RC parameters (R 0 , R 1 , C 1 ), gives a less RMSE as compared to the non-temperature-dependent model. Both model with and without ambient temperature term used the same U oc in (6). The non-temperature-dependent model adopted the fixed RC parameters [51,[64][65][66]: R 0 = 0.0823, R 1 = 0.0042, C 1 = 52.2. 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].

Conclusions
The effects of the ambient temperature, cell temperature, hysteresis voltage and thermal aging on capacity due to multiple charging and discharging was considered in the battery cell model of the cylindrical lithium iron phosphate battery cells (ANR26650M1-B). The temperature-dependent extended Kalman filter (EKF) algorithm was used to estimate the state-of-charge (SOC) of the battery cell. The effectiveness of the proposed scheme was validated experimentally on actual battery cell under different ambient temperatures. As compared with cell model without the ambient temperature term in the SOC model, the proposed battery temperature-dependent model gave a less root-mean-square error at low (i.e., 5 • C and 15 • C) and high ambient temperature (45 • C). For future works, the battery packs will be simulated and tested. The reliability of the SOC estimation under varying conditions of the battery pack will be examined. The prismatic or pouch cell will be examined.
Author Contributions: C.S.C. defined the flow of the paper, theoretical concepts of the unified battery cell model and model simulation and verification that drove this research with Z.G. and J.H.K.C. on cell temperature measurements. C.Z. provided the academic support and advice during the project. All authors discussed and provided comments at all stages.
Funding: This research received no external funding.