Detection of Internal Short Circuit in Lithium Ion Battery Using Model-Based Switching Model Method

Minhwan Seo 1, Taedong Goh 2, Minjun Park 1, Gyogwon Koo 1 and Sang Woo Kim 1,* 1 Department of Electrical Engineering, Pohang University of Science and Technology, 77 Cheongam-Ro, Nam-Gu, Pohang 37673, Korea; mhseo09@postech.edu (M.S.); parkmj@postech.edu (M.P.); ggkoo99@postech.edu (G.K.) 2 Department of Creative IT Excellence Engineering and Future IT Innovation Laboratory, Pohang University of Science and Technology, 77 Cheongam-Ro, Nam-Gu, Pohang 37673, Korea; ehd1116@postech.edu * Correspondence: swkim@postech.edu; Tel.: +82-54-279-2237; Fax: +82-54-279-2903


Introduction
Li-ion batteries have high power density, high energy efficiency and a long cycle life [1], and are therefore used as electric energy storage and power sources for electric devices and electric-drive vehicles.However, the Li-ion battery can develop dangerous malfunctions [2,3] such as internal short circuit (ISCr) [4,5] and cell reversal [6]; the main causes of these phenomena are overcharge [7] and overdischarge [8].The ISCr may cause thermal runaway when the temperature rise by the ISCr in the battery exceeds a certain point [5] or the ISCr resistance R ISC f is lower than a certain value [9].Then a fire and an explosion can occur by the thermal runaway [10][11][12].The ISCr is the main cause of battery fire accidents in Boeing 787-8 aircraft [13].Therefore, a method to detect the ISCr is necessary before the thermal runaway happens in the Li-ion battery.
For these reasons, studies to detect the ISCr have been presented [14][15][16][17].The ISCr can be detected by determining certain thresholds such as reduction of terminal voltage and increase of batteries temperature [14], but to obtain the thresholds, this method requires prior ISCr tests with batteries.Therefore, model-based algorithms have been presented to detect the ISCr by identifying variations of parameters in the model [15,16].Using equivalent circuit models of a normal battery and a battery with ISCr as thresholds, characteristic parameters are obtained to detect ISCr in a battery pack [15], but this method can be used only when the battery with ISCr is connected to several normal batteries Energies 2017, 10, 76 2 of 13 in series, and the terminal voltages of both the battery with ISCr and the normal batteries are provided.The ISCr can be detected by using variation of estimated parameters in the equivalent circuit model and the energy balance equation [16], but this method must be verified with other load current profiles to check whether the estimated parameters show similar variation.When the ISCr occurs in the Li-ion battery, the terminal voltage increases once the battery is recharged, but the voltage reaches a stable value [18]; R ISC f can be calculated using the charging current and the terminal voltage.However, calculation of R ISC f by this method requires knowledge of the specific charging current that makes the terminal voltage reach the stable value in the battery with ISCr.
An early version [17] of the algorithm proposed in this paper estimated R ISC f by using self-discharge from the ISCr to detect it, but the accuracy of R ISC f estimates was low; to solve this problem, this paper presents a model-based switching model method (SMM).As a fault index, R ISC f is estimated accurately by using SMM to detect the ISCr.To verify the proposed algorithm, environments of simulation and experiment are configured and two load current profiles: dynamic stress test (DST) and urban dynamometer driving schedule (UDDS) are used.The proposed algorithm is explained in Section 2, the environments of simulation and experiment are introduced in Section 3. The results of the simulation and the experiment are presented and discussed in Section 4. Finally, the conclusion and outline of future work are presented in Section 5.

Switching Model Method
In accordance with the estimated state of charge (SOC) defined as the present capacity of the battery as a proportion of its total capacity, the model of Li-ion battery with ISCr is switched to the updated model of Li-ion battery with ISCr.If the variation between initial estimated SOC and current estimated SOC is ≥0.2, the R ISC f is estimated.Estimates of R ISC f fluctuate due to variation in load currents, so as a fault index, the mean R ISC f of estimated R ISC f s is used to detect the ISCr instead of the estimated R ISC f .R ISC f is calculated from the previous estimated R ISC f s and the current estimated R ISC f , then used to change the model with ISCr to the updated model with ISCr.Then, using the updated model improves the accuracy of open circuit voltage (OCV) estimates, so R ISC f is estimated accurately.Iteratively, the next R ISC f is calculated using the next estimated R ISC f and used to update the model again.We call this process the SMM (Figure 1).
normal batteries in series, and the terminal voltages of both the battery with ISCr and the normal batteries are provided.The ISCr can be detected by using variation of estimated parameters in the equivalent circuit model and the energy balance equation [16], but this method must be verified with other load current profiles to check whether the estimated parameters show similar variation.When the ISCr occurs in the Li-ion battery, the terminal voltage increases once the battery is recharged, but the voltage reaches a stable value [18];   can be calculated using the charging current and the terminal voltage.However, calculation of   by this method requires knowledge of the specific charging current that makes the terminal voltage reach the stable value in the battery with ISCr.
An early version [17] of the algorithm proposed in this paper estimated   by using selfdischarge from the ISCr to detect it, but the accuracy of   estimates was low; to solve this problem, this paper presents a model-based switching model method (SMM).As a fault index,   is estimated accurately by using SMM to detect the ISCr.To verify the proposed algorithm, environments of simulation and experiment are configured and two load current profiles: dynamic stress test (DST) and urban dynamometer driving schedule (UDDS) are used.The proposed algorithm is explained in Section 2, the environments of simulation and experiment are introduced in Section 3. The results of the simulation and the experiment are presented and discussed in Section 4. Finally, the conclusion and outline of future work are presented in Section 5.

Switching Model Method
In accordance with the estimated state of charge (SOC) defined as the present capacity of the battery as a proportion of its total capacity, the model of Li-ion battery with ISCr is switched to the updated model of Li-ion battery with ISCr.If the variation between initial estimated SOC and current estimated SOC is ≥0.2, the   is estimated.Estimates of   fluctuate due to variation in load currents, so as a fault index, the mean   ̅̅̅̅̅̅̅ of estimated   s is used to detect the ISCr instead of the estimated   .  ̅̅̅̅̅̅̅ is calculated from the previous estimated   s and the current estimated   , then used to change the model with ISCr to the updated model with ISCr.Then, using the updated model improves the accuracy of open circuit voltage (OCV) estimates, so   is estimated accurately.Iteratively, the next   ̅̅̅̅̅̅̅ is calculated using the next estimated   and used to update the model again.We call this process the SMM (Figure 1).

Equivalent Circuit Models
A normal Li-ion battery can be represented by an equivalent circuit model (Figure 2a) [15,17,19] that consists of OCV   , internal resistance , load current  and terminal voltage   .The Li-ion Energies 2017, 10, 76 3 of 13

Equivalent Circuit Models
A normal Li-ion battery can be represented by an equivalent circuit model (Figure 2a) [15,17,19] that consists of OCV V oc , internal resistance R, load current I and terminal voltage V t .The Li-ion battery with ISCr can be represented by a similar equivalent circuit model (Figure 2b) [15,17,20] where I 1 f is the current that flows within the battery, and I 2 f is the current that flows through the R ISC f .
Energies 2017, 10, 76 3 of 13 battery with ISCr can be represented by a similar equivalent circuit model (Figure 2b) [15,17,20] where  1 is the current that flows within the battery, and  2 is the current that flows through the   .
(a) (b) This model with ISCr has been verified to mimic the ISCr [17].Especially, subscript  (fault) is used to distinguish between parameters related to the model with ISCr and parameters of the normal battery model.The model with ISCr is described in Equations ( 1) and ( 2) by Ohm's law with a discretization step [15]:

OCV Estimation
The recursive least squares (RLS) algorithm is usually used to estimate parameters in the normal battery model [19].In this paper, the RLS algorithm was used to estimate model parameters such as   and   (Figure 2b).The initial value of the covariance matrix was [500 −250; −250 210], and initial values of the parameter vector   that contains estimated parameters were the terminal voltage measured at the first time and 0.05.The forgetting factor is typically a value between 0.95 and 1 to get finely estimated parameters.In this study, the forgetting factor was set to 0.9995.To distinguish the OCV estimates from the two models with ISCr,  , is the estimated OCV from the model with ISCr, and  , is the estimated OCV from the updated model with ISCr.The equation used in the RLS algorithm is  This model with ISCr has been verified to mimic the ISCr [17].Especially, subscript f (fault) is used to distinguish between parameters related to the model with ISCr and parameters of the normal battery model.The model with ISCr is described in Equations ( 1) and ( 2) by Ohm's law with a discretization step [15]:

OCV Estimation
The recursive least squares (RLS) algorithm is usually used to estimate parameters in the normal battery model [19].In this paper, the RLS algorithm was used to estimate model parameters such as V oc f and R f (Figure 2b).The initial value of the covariance matrix was [500 −250; −250 210], and initial values of the parameter vector θ f that contains estimated parameters were the terminal voltage measured at the first time and 0.05.The forgetting factor is typically a value between 0.95 and 1 to get finely estimated parameters.In this study, the forgetting factor was set to 0.9995.To distinguish the OCV estimates from the two models with ISCr, V OC f ,pre is the estimated OCV from the model with ISCr, and V OC f ,upd is the estimated OCV from the updated model with ISCr.The equation used in the RLS algorithm is where y f is a measurable quantity and ∅ f is a vector of known quantities.The RLS algorithm cannot estimate R ISC f directly, because the number of unknown parameters is bigger than the number of known parameters.In Equation (3), θ f has two estimated parameters, , which are combined with three unknown parameters, V OC f ,pre , R f and R ISC f .Therefore, using the assumption [18] that the first parameter pre can be estimated from the RLS algorithm.
Energies 2017, 10, 76 4 of 13 However, this approximation-assumption causes error that reduces the accuracy of R ISC f estimates.Therefore, the updated model described in Equations ( 4) and ( 5) must be used to avoid this assumption.
If the initial R ISC f is estimated at the point, at which the model with ISCr is switched to the updated model with ISCr, the unknown parameter R ISC f is assumed to be the initial estimated R ISC f .The V OC f ,pre is used to estimate the initial R ISC f , and the method to estimate SOC and R ISC f will be explained in Section 2.3.Then the current I 1 f can be calculated using Equation ( 4) and used as input data of the RLS algorithm in Equation ( 5).Once the input data changes from I f to I 1 f , the estimated parameters of θ f also change and the V OC f ,upd can be estimated directly without the approximation-assumption.Once the next R ISC f is estimated using V OC f ,upd , R ISC f is calculated and then used to update the model and estimate V OC f ,upd iteratively.
After switching the model, the estimated OCVs with SMM began to be more accurate than the estimated OCVs without SMM (Figure 3) because of elimination of the error from the approximation-assumption.The true OCVs used for verification were calculated using the true value of R ISC f , the coulomb counting method and the relation between OCV and SOC.
Energies 2017, 10, 76 4 of 13 estimates.Therefore, the updated model described in Equations ( 4) and ( 5) must be used to avoid this assumption.
If the initial   is estimated at the point, at which the model with ISCr is switched to the updated model with ISCr, the unknown parameter   is assumed to be the initial estimated   .The  , is used to estimate the initial   , and the method to estimate SOC and   will be explained in Section 2.3.Then the current  1 can be calculated using Equation ( 4) and used as input data of the RLS algorithm in Equation ( 5).Once the input data changes from   to  1 , the estimated parameters of   also change and the  , can be estimated directly without the approximationassumption.Once the next   is estimated using  , ,   ̅̅̅̅̅̅̅ is calculated and then used to update the model and estimate  , iteratively.After switching the model, the estimated OCVs with SMM began to be more accurate than the estimated OCVs without SMM (Figure 3) because of elimination of the error from the approximationassumption.The true OCVs used for verification were calculated using the true value of   , the coulomb counting method and the relation between OCV and SOC.

SOC and 𝑅 𝐼𝑆𝐶𝑓 Estimation
The estimated SOCs can be obtained from the relation between OCV and SOC (Figure 4).Because of additional self-discharge due to  2 flowing through the   , the estimated SOCs decline more in the battery with ISCr than in the normal battery.Furthermore, decrease in   represents increase in the severity of the ISCr in the Li-ion battery and in the decline of estimated SOCs.Using the self-discharge phenomenon,   can be estimated [18].
The coulomb counting method is usually used to calculate true SOCs from the load current, true initial SOC and total capacity [21].This method uses Equation ( 6) with a discretization step, where   is the maximum capacity of the battery, and  is the sampling rate.To eliminate the unknown term   (0) ,   () is subtracted from   ( + 1), and  2 is replaced with  2 =     in Equation (6).

SOC and R ISC f Estimation
The estimated SOCs can be obtained from the relation between OCV and SOC (Figure 4).Because of additional self-discharge due to I 2 f flowing through the R ISC f , the estimated SOCs decline more in the battery with ISCr than in the normal battery.Furthermore, decrease in R ISC f represents increase in the severity of the ISCr in the Li-ion battery and in the decline of estimated SOCs.Using the self-discharge phenomenon, R ISC f can be estimated [18].
The coulomb counting method is usually used to calculate true SOCs from the load current, true initial SOC and total capacity [21].This method uses Equation ( 6) with a discretization step, where C max is the maximum capacity of the battery, and T is the sampling rate.To eliminate the unknown term SOC f (0), SOC f (k) is subtracted from SOC f (k + 1), and Energies 2017, 10, 76 5 of 13 Exact estimation of the   by using the variation of estimated SOCs is difficult in the short interval between k + 1 and k.Therefore, the interval must be increased by adding the k − 1th to p + 1th equations to the kth equation in Equation (7) to clearly show the self-discharge phenomenon where k is a current iteration and p is an initial iteration (p + 2 < k).
( + 1) −   () =      ( + 1) −      ( + 1) (7) The first term on the right side of Equation ( 8) describes the discharge from the load current, and the second term represents the self-discharge from the ISCr.  can be estimated using Equation (9), which is a rearrangement of Equation (8).
The choice of time to estimate the   for switching the model is important, because if ISCr is 50 Ω, 30 Ω or 20 Ω, the self-discharge from the ISCr is not observed dominantly in the early iterations of the process of estimating SOC (Figure 4); i.e., the ratio of decrease in SOC due to self-discharge to total decrease in SOC must be large enough to reduce the effect of errors of the estimated SOCs and to clearly show the effect of self-discharge from ISCr.Accordingly, we determine that the model should be switched when the variation between initial estimated SOC and current estimated SOC is ≥0.2; i.e., 20% of the total capacity of the Li-ion battery.

Load Current Profiles
Two load current profiles were used as input data to the simulation and the experiments (Figure 5).We named these current profiles DST 5 A and UDDS 5 A; both have the minimum value of −5 A. We also used both DST 3 A and UDDS 3 A to verify the proposed algorithm with various data.Exact estimation of the R ISC f by using the variation of estimated SOCs is difficult in the short interval between k + 1 and k.Therefore, the interval must be increased by adding the k − 1th to p + 1th equations to the kth equation in Equation (7) to clearly show the self-discharge phenomenon where k is a current iteration and p is an initial iteration (p The first term on the right side of Equation ( 8) describes the discharge from the load current, and the second term represents the self-discharge from the ISCr.R ISC f can be estimated using Equation (9), which is a rearrangement of Equation (8).
The choice of time to estimate the R ISC f for switching the model is important, because if ISCr is 50 Ω, 30 Ω or 20 Ω, the self-discharge from the ISCr is not observed dominantly in the early iterations of the process of estimating SOC (Figure 4); i.e., the ratio of decrease in SOC due to self-discharge to total decrease in SOC must be large enough to reduce the effect of errors of the estimated SOCs and to clearly show the effect of self-discharge from ISCr.Accordingly, we determine that the model should be switched when the variation between initial estimated SOC and current estimated SOC is ≥0.2; i.e., 20% of the total capacity of the Li-ion battery.

Load Current Profiles
Two load current profiles were used as input data to the simulation and the experiments (Figure 5).We named these current profiles DST 5 A and UDDS 5 A; both have the minimum value of −5 A. We also used both DST 3 A and UDDS 3 A to verify the proposed algorithm with various data.

Configuration of Simulation Environment
In this study, a first-order RC model [22] was used to build a simulation model.To represent the ISCr,   was connected in parallel at the terminal of the first-order RC model (Figure 6).The simulation model was configured using MATLAB/Simulink [23,24].Resistance  0 , resistance  1 and capacitance  1 (Figure 6) were estimated using the RLS algorithm with experimental data of cell A and DST 5 A [19], then used to build the simulation model.

Configuration of Experiment Environment
In this study, two identical cells A and B (Table 1) were used to configure the experiment and to get data in various environments.The temperature was about 18-26 °C when cells A and B were tested.The initial SOC was set to 70% in both cells.The load current profiles were differently applied to each battery.DST 5 A and UDDS 5 A were used in experiments with cell A, and DST 3 A and UDDS 3 A were used for cell B. To prevent batteries from overdischarge, these load current profiles were applied to the batteries until their SOCs reached 10% of total capacity.To make various values of resistance such as 50 Ω, 30 Ω, 20 Ω, 10 Ω and 5 Ω, the five 10 Ω resistances were combined.The tolerance of the 10 Ω resistance was ±5%.The true values of these resistances were measured such as 49.91 Ω, 29.93 Ω, 19.92 Ω, 9.95 Ω and 4.98 Ω respectively, and used to calculate the relative error in the experimental data.When the load current profile was applied to the cell, the switch was used to connect the cell and resistances 50 Ω, 30 Ω, 20 Ω, 10 Ω or 5 Ω in parallel to represent a battery with ISCr.In 10 seconds after the load current was sent, the switch was

Configuration of Simulation Environment
In this study, a first-order RC model [22] was used to build a simulation model.To represent the ISCr, R ISC f was connected in parallel at the terminal of the first-order RC model (Figure 6).The simulation model was configured using MATLAB/Simulink [23,24].Resistance R 0 f , resistance R 1 f and capacitance C 1 f (Figure 6) were estimated using the RLS algorithm with experimental data of cell A and DST 5 A [19], then used to build the simulation model.

Configuration of Simulation Environment
In this study, a first-order RC model [22] was used to build a simulation model.To represent the ISCr,   was connected in parallel at the terminal of the first-order RC model (Figure 6).The simulation model was configured using MATLAB/Simulink [23,24].Resistance  0 , resistance  1 and capacitance  1 (Figure 6) were estimated using the RLS algorithm with experimental data of cell A and DST 5 A [19], then used to build the simulation model.

Configuration of Experiment Environment
In this study, two identical cells A and B (Table 1) were used to configure the experiment and to get data in various environments.The temperature was about 18-26 °C when cells A and B were tested.The initial SOC was set to 70% in both cells.The load current profiles were differently applied to each battery.DST 5 A and UDDS 5 A were used in experiments with cell A, and DST 3 A and UDDS 3 A were used for cell B. To prevent batteries from overdischarge, these load current profiles were applied to the batteries until their SOCs reached 10% of total capacity.To make various values of resistance such as 50 Ω, 30 Ω, 20 Ω, 10 Ω and 5 Ω, the five 10 Ω resistances were combined.The tolerance of the 10 Ω resistance was ±5%.The true values of these resistances were measured such as 49.91 Ω, 29.93 Ω, 19.92 Ω, 9.95 Ω and 4.98 Ω respectively, and used to calculate the relative error in the experimental data.When the load current profile was applied to the cell, the switch was used to connect the cell and resistances 50 Ω, 30 Ω, 20 Ω, 10 Ω or 5 Ω in parallel to represent a battery with ISCr.In 10 seconds after the load current was sent, the switch was

Configuration of Experiment Environment
In this study, two identical cells A and B (Table 1) were used to configure the experiment and to get data in various environments.The temperature was about 18-26 • C when cells A and B were tested.The initial SOC was set to 70% in both cells.The load current profiles were differently applied to each battery.DST 5 A and UDDS 5 A were used in experiments with cell A, and DST 3 A and UDDS 3 A were used for cell B. To prevent batteries from overdischarge, these load current profiles were applied to the batteries until their SOCs reached 10% of total capacity.To make various values of resistance such as 50 Ω, 30 Ω, 20 Ω, 10 Ω and 5 Ω, the five 10 Ω resistances were combined.The tolerance of the 10 Ω resistance was ±5%.The true values of these resistances were measured such as 49.91 Ω, 29.93 Ω, 19.92 Ω, 9.95 Ω and 4.98 Ω respectively, and used to calculate the relative error in the experimental data.When the load current profile was applied to the cell, the switch was used to connect the cell and resistances 50 Ω, 30 Ω, 20 Ω, 10 Ω or 5 Ω in

Comparison between Results with SMM and without SMM
To ensure fair comparison between R ISC f estimates with and without SMM from the same point, the R ISC f s were not estimated in the method without SMM before the model with ISCr was switched to the updated model with ISCr in the method with SMM.After switching the model, R ISC f with SMM began to converge more accurately on the true value of 10 Ω than R ISC f without SMM did (Figure 8).The reason of this superiority is that the estimated OCVs and SOCs became accurate due to the updated model, which removed the error imposed by the approximation-assumption.It is also reason that the R ISC f was estimated with the accurate SOC estimates.When the proposed algorithm without the assumption was used, the accuracy of R ISC f estimates was generally improved and the relative error of the final value of R ISC f decreased greatly (Tables 2 and 3).Because the decrease in the magnitude of true R ISC f represented increase in the error from the assumption, the difference between relative errors with and without SMM increased as the magnitude of true R ISC f decreased.
Energies 2017, 10, 76 7 of 13 turned on.Therefore, the load current and the terminal voltage were measured after 10 s; sample interval was 0.1 s.

Relation between OCV and SOC Test
The relation between OCV and SOC was obtained by a prior test [25] and is necessary to use the proposed algorithm.After the battery had been charged fully, it was rested for 3600 s to obtain a value of OCV that is equal to the terminal voltage.Then the battery was discharged with 0.5 C for 720 s to set the SOC to 90%, then rested for 3600 s to get the value of OCV.The OCV-SOC curve (Figure 7) could be obtained by repeating the process until the value of SOC reached 0%.

Comparison between Results with SMM and without SMM
To ensure fair comparison between   ̅̅̅̅̅̅̅ estimates with and without SMM from the same point, the   s were not estimated in the method without SMM before the model with ISCr was switched to the updated model with ISCr in the method with SMM.After switching the model,   ̅̅̅̅̅̅̅ with SMM began to converge more accurately on the true value of 10 Ω than   ̅̅̅̅̅̅̅ without SMM did (Figure 8).The reason of this superiority is that the estimated OCVs and SOCs became accurate due to the updated model, which removed the error imposed by the approximation-assumption.It is also reason that the   was estimated with the accurate SOC estimates.When the proposed algorithm without the assumption was used, the accuracy of   estimates was generally improved and the relative error of the final value of   ̅̅̅̅̅̅̅ decreased greatly (Tables 2 and 3).Because the decrease in the magnitude of true   represented increase in the error from the assumption, the difference between relative errors with and without SMM increased as the magnitude of true   decreased.

Effect of Magnitude of True R ISC f in the Simulation
Load current profiles DST 5 A and UDDS 5 A were used to execute the simulation model that represented cell A with ISCr.Initial SOC of the simulation model was 70%, like the configuration of the experiment with cells A and B.
In cases ISCr 50 Ω and 30 Ω for DST 5 A and ISCr 50 Ω, 30 Ω and 20 Ω for UDDS 5 A, the R ISC f s fluctuated much more than other ISCr faults (Figure 9) because the effect of self-discharge from ISCr was too small to be represented in the estimated outcomes like the normal battery.

Effect of Magnitude of True 𝑅 𝐼𝑆𝐶𝑓 in the Simulation
Load current profiles DST 5 A and UDDS 5 A were used to execute the simulation model that represented cell A with ISCr.Initial SOC of the simulation model was 70%, like the configuration of the experiment with cells A and B.
In cases ISCr 50 Ω and 30 Ω for DST 5 A and ISCr 50 Ω, 30 Ω and 20 Ω for UDDS 5 A, the   ̅̅̅̅̅̅̅ s fluctuated much more than other ISCr faults (Figure 9) because the effect of self-discharge from ISCr was too small to be represented in the estimated outcomes like the normal battery.Although the fluctuations were high in the early stage, the magnitude of the fluctuations gradually decreased and R ISC f s converged on values near the true R ISC f .The relative error of ISCr 5 and Energies 2017, 10, 76 9 of 13 10 Ω was large although the self-discharge from ISCr 5 Ω and 10 Ω was influential.Because the self-discharge from ISCr 5 Ω and 10 Ω largely increased the decline of SOC in the early stage of interations (Figure 4), the initial estimated SOC had estimation error.This error affected the accuracy of R ISC f estimates for ISCr 5 Ω and 10 Ω.The final value of R ISC f had relative error ≤14.2% (Table 4) and the ISCr could be detected early with high accuracy of R ISC f estimates.R ISC f also fluctuated in experimental results but also decreased and converged on values near the true R ISC f (Figure 10).The main difference between simulation results and experiment results was that the maximum relative error of ISCr 50 Ω and 30 Ω significantly increased from 14.2% in the simulation to 49.3% in cell A (Table 5).
Energies 2017, 10, 76 9 of 13 interations (Figure 4), the initial estimated SOC had estimation error.This error affected the accuracy of   estimates for ISCr 5 Ω and 10 Ω.The final value of   ̅̅̅̅̅̅̅ had relative error ≤14.2% (Table 4) and the ISCr could be detected early with high accuracy of   estimates.10).The main difference between simulation results and experiment results was that the maximum relative error of ISCr 50 Ω and 30 Ω significantly increased from 14.2% in the simulation to 49.3% in cell A (Table 5).Especially, when true   was 50 Ω or 30 Ω, the effect of self-discharge from ISCr relatively decreased because the error of estimated SOCs increased due to noise in the experimental environment.Despite this noise, the ISCr 20 Ω, 10 Ω, and 5 Ω had relative error ≤19.7% because the effect of self-discharge from ISCr was much bigger than the increase of the error in estimated SOCs.Therefore, the ISCr fault could be detected early before thermal runaway occurred in the Li-ion battery with ISCr.However, the low accuracy of estimated   in ISCr 50 Ω and 30 Ω remains a problem; to overcome it, the effect of self-discharge from ISCr must be increased by decreasing the C-rate of the load current profiles to increase the time over which the battery completely discharges.Especially, when true R ISC f was 50 Ω or 30 Ω, the effect of self-discharge from ISCr relatively decreased because the error of estimated SOCs increased due to noise in the experimental environment.Despite this noise, the ISCr 20 Ω, 10 Ω, and 5 Ω had relative error ≤19.7% because the effect of self-discharge from ISCr was much bigger than the increase of the error in estimated SOCs.Therefore, the ISCr fault could be detected early before thermal runaway occurred in the Li-ion battery with ISCr.However, the low accuracy of estimated R ISC f in ISCr 50 Ω and 30 Ω remains a problem; to overcome it, the effect of self-discharge from ISCr must be increased by decreasing the C-rate of the load current profiles to increase the time over which the battery completely discharges.The experiment with cell B was conducted using DST 3 A and UDDS 3 A, which were the load current profiles with low C-rate.When the C-rate of load current profiles decreased, the area between the estimated SOCs of normal cell B and the estimated SOCs of cell B with ISCr 50 Ω increased more than that of cell A (Figure 11).The area represents the decline of SOC due to self-discharge by ISCr.Therefore, this change increased the influence of self-discharge, and the accuracy of the estimated R ISC f in ISCr 50 Ω and 30 Ω was improved (Figure 12).Accordingly, the maximum relative error of ISCr 50 Ω and 30 Ω decreased greatly from 49.3% in cell A to 22.1% in cell B (Table 6), and the relative error was ≤26.1%.The experiment with cell B was conducted using DST 3 A and UDDS 3 A, which were the load current profiles with low C-rate.When the C-rate of load current profiles decreased, the area between the estimated SOCs of normal cell B and the estimated SOCs of cell B with ISCr 50 Ω increased more than that of cell A (Figure 11).The area represents the decline of SOC due to self-discharge by ISCr.Therefore, this change increased the influence of self-discharge, and the accuracy of the estimated   in ISCr 50 Ω and 30 Ω was improved (Figure 12).Accordingly, the maximum relative error of ISCr 50 Ω and 30 Ω decreased greatly from 49.3% in cell A to 22.1% in cell B (Table 6), and the relative error was ≤26.1%.The experiment with cell B was conducted using DST 3 A and UDDS 3 A, which were the load current profiles with low C-rate.When the C-rate of load current profiles decreased, the area between the estimated SOCs of normal cell B and the estimated SOCs of cell B with ISCr 50 Ω increased more than that of cell A (Figure 11).The area represents the decline of SOC due to self-discharge by ISCr.Therefore, this change increased the influence of self-discharge, and the accuracy of the estimated   in ISCr 50 Ω and 30 Ω was improved (Figure 12).Accordingly, the maximum relative error of ISCr 50 Ω and 30 Ω decreased greatly from 49.3% in cell A to 22.1% in cell B (Table 6), and the relative error was ≤26.1%.

Effect of Variation in OCV-SOC Curve
The OCV-SOC curve can be changed because of the capacity fade of Li-ion battery caused by the cycle aging and the calendar aging [3].In this study, the amount of capacity fade of the aged battery (cell A) was 2.7% of total capacity.The OCV-SOC curves for both fresh cell A and aged cell A were almost equal.Therefore, the OCV-SOC curves for fresh cells A and B were used in the proposed algorithm.However, the severe capacity fade can cause significant change in the OCV-SOC curve.This change can generate the considerable error of estimated SOCs which results in the error of   estimates in Equation (9).
In summary, when the load current pofiles with high C-rate were used, the ISCr fault in ISCr 20~5 Ω range could be detected early before thermal runaway happened in the Li-ion battery with ISCr.Furthermore, when the load current profiles with low C-rate were used to increase the effect of self-discharge in ISCr 50~30 Ω, the proposed algorithm could detect the ISCr fault early in ISCr 50~30 Ω with high accuracy of the   estimates.In addition, the study considering the variation in the OCV-SOC curve should be proceeded continuously to improve the accuracy of the   estimates.

Conclusions
In this paper, a model-based SMM is introduced to detect ISCr in the Li-ion battery.Using the equivalent circuit model of the battery with ISCr and the RLS algorithm, the OCV is estimated.The SOC is estimated by using its relationship with OCV.Then   is estimated using the selfdischarge phenomenon of the ISCr.The SMM greatly increased the accuracy of the estimated   .The proposed algorithm was verified in simulations and experiments using two load current profiles.The effect of the magnitude of true   on estimated   and the effect of C-rate of load current on estimated   were analyzed.The   can be estimated with high accuracy using the proposed algorithm, and as fault index, the   ̅̅̅̅̅̅̅ can be used to detect the ISCr early.Our future research will concentrate on extending our proposed algorithm to detection of ISCr in an aged battery and a battery pack.

Effect of Variation in OCV-SOC Curve
The OCV-SOC curve can be changed because of the capacity fade of Li-ion battery caused by the cycle aging and the calendar aging [3].In this study, the amount of capacity fade of the aged battery (cell A) was 2.7% of total capacity.The OCV-SOC curves for both fresh cell A and aged cell A were almost equal.Therefore, the OCV-SOC curves for fresh cells A and B were used in the proposed algorithm.However, the severe capacity fade can cause significant change in the OCV-SOC curve.This change can generate the considerable error of estimated SOCs which results in the error of R ISC f estimates in Equation (9).
In summary, when the load current pofiles with high C-rate were used, the ISCr fault in ISCr 20~5 Ω range could be detected early before thermal runaway happened in the Li-ion battery with ISCr.Furthermore, when the load current profiles with low C-rate were used to increase the effect of self-discharge in ISCr 50~30 Ω, the proposed algorithm could detect the ISCr fault early in ISCr 50~30 Ω with high accuracy of the R ISC f estimates.In addition, the study considering the variation in the OCV-SOC curve should be proceeded continuously to improve the accuracy of the R ISC f estimates.

Conclusions
In this paper, a model-based SMM is introduced to detect ISCr in the Li-ion battery.Using the equivalent circuit model of the battery with ISCr and the RLS algorithm, the OCV is estimated.The SOC is estimated by using its relationship with OCV.Then R ISC f is estimated using the self-discharge phenomenon of the ISCr.The SMM greatly increased the accuracy of the estimated R ISC f .The proposed algorithm was verified in simulations and experiments using two load current profiles.The effect of the magnitude of true R ISC f on estimated R ISC f and the effect of C-rate of load current on estimated R ISC f were analyzed.The R ISC f can be estimated with high accuracy using the proposed algorithm, and as fault index, the R ISC f can be used to detect the ISCr early.Our future research will concentrate on extending our proposed algorithm to detection of ISCr in an aged battery and a battery pack.

Figure 1 .
Figure 1.The scheme of the switching model method.

Figure 1 .
Figure 1.The scheme of the switching model method.

Figure 6 .
Figure 6.Configuration of simulation model: the first-order RC model with   .

Figure 6 .
Figure 6.Configuration of simulation model: the first-order RC model with   .

Figure 6 .
Figure 6.Configuration of simulation model: the first-order RC model with R ISC f .

Figure 8 . 8 .
Figure 8.Comparison of   ̅̅̅̅̅̅̅ with and without SMM (experiment: cell A, DST 5 A and ISCr 10 Ω). 8. Comparison of R ISC f with and without SMM (experiment: cell A, DST 5 A and ISCr 10

Figure 9 .
Figure 9.   ̅̅̅̅̅̅̅ in the simulation model with ISCr 50, 30, 20, 10 and 5 Ω: (a) DST 5 A; (b) UDDS 5 A. Although the fluctuations were high in the early stage, the magnitude of the fluctuations gradually decreased and   ̅̅̅̅̅̅̅ s converged on values near the true   .The relative error of ISCr 5 Ω and 10 Ω was large although the self-discharge from ISCr 5 Ω and 10 Ω was influential.Because the self-discharge from ISCr 5 Ω and 10 Ω largely increased the decline of SOC in the early stage of

Figure 11 .Figure 11 .
Figure 11.Comparison of estimated SOCs of normal cells and of cells with ISCr 50 Ω (experiments: DST 5 A for cell A and DST 3 A for cell B).

Figure 11 .Figure 12 Figure 12 .
Figure 11.Comparison of estimated SOCs of normal cells and of cells with ISCr 50 Ω (experiments: DST 5 A for cell A and DST 3 A for cell B).
, ,   and   . of   can approximate  , because the   >>   , the  , can be estimated from the RLS algorithm.However, this approximation-assumption causes error that reduces the accuracy of

Table 1 .
Key specification of Li-ion battery.

Table 1 .
Key specification of Li-ion battery.

Table 1 .
Key specification of Li-ion battery.

Table 2 .
Relative error (%) of the final value of R ISC f depending on the ISCr faults in the experiment with cell A and DST 5 A.

Table 3 .
Relative error (%) of the final value of R ISC f depending on the ISCr faults in the experiment with cell A and UDDS 5 A.

Table 2 .
Relative error (%) of the final value of   ̅̅̅̅̅̅̅ depending on the ISCr faults in the experiment with cell A and DST 5 A.

Table 3 .
Relative error (%) of the final value of   ̅̅̅̅̅̅̅ depending on the ISCr faults in the experiment with cell A and UDDS 5 A.

Table 4 .
Relative error (%) of the final value of R ISC f depending on the ISCr faults in the simulation model.Effect of Magnitude of True R ISC f in the Experiment with Cell A

Table 4 .
Relative error (%) of the final value of   ̅̅̅̅̅̅̅ depending on the ISCr faults in the simulation model.Effect of Magnitude of True   in the Experiment with Cell A   ̅̅̅̅̅̅̅ also fluctuated in experimental results but also decreased and converged on values near the true   (Figure

Table 5 .
Relative error (%) of the final value of R ISC f depending on the ISCr faults in the experiment with cell A.

Table 5 .
Relative error (%) of the final value of   ̅̅̅̅̅̅̅ depending on the ISCr faults in the experiment with cell A.

Table 5 .
Relative error (%) of the final value of   ̅̅̅̅̅̅̅ depending on the ISCr faults in the experiment with cell A.

Table 6 .
Relative error (%) of the final value of   ̅̅̅̅̅̅̅ depending on the ISCr faults in the experiment with cell B.

Table 6 .
Relative error (%) of the final value of R ISC f depending on the ISCr faults in the experiment with cell B.