Detection Method for Soft Internal Short Circuit in Lithium-Ion Battery Pack by Extracting Open Circuit Voltage of Faulted Cell

Early detection of internal short circuit which is main cause of thermal runaway in a lithium-ion battery is necessary to ensure battery safety for users. As a promising fault index, internal short circuit resistance can directly represent degree of the fault because it describes self-discharge phenomenon caused by the internal short circuit clearly. However, when voltages of individual cells in a lithium-ion battery pack are not provided, the effect of internal short circuit in the battery pack is not readily observed in whole terminal voltage of the pack, leading to difficulty in estimating accurate internal short circuit resistance. In this paper, estimating the resistance with the whole terminal voltages and the load currents of the pack, a detection method for the soft internal short circuit in the pack is proposed. Open circuit voltage of a faulted cell in the pack is extracted to reflect the self-discharge phenomenon obviously; this process yields accurate estimates of the resistance. The proposed method is verified with various soft short conditions in both simulations and experiments. The error of estimated resistance does not exceed 31.2% in the experiment, thereby enabling the battery management system to detect the internal short circuit early.


Introduction
Lithium-ion batteries are widely used as a power source in electric devices and electric vehicles [1,2], due to their high power density, high energy efficiency and excellent cycle stability [3,4].The demand for them is expected to rise continuously in the coming years [5,6].However, safety concerns related to the lithium-ion batteries still remain [7][8][9] because hazardous incidents such as fire accidents in the Boeing 787-8 aircraft [10] and battery failures in the Samsung Note7 [11] have been frequently reported by media [12].The main cause of these two events is internal short circuit (ISCr) in the lithium-ion batteries.The ISCr can be caused by manufacturing defects [13,14] and various types of abuse such as overcharge [15,16] and overdischarge [17].Furthermore, when a magnitude of ISCr resistance (R ISCr ) is lower than a particular value [18], a temperature of the battery exceeds a certain point due to the ISCr [19][20][21].Then, the battery may experience thermal runaway with a fire and can even explode [22][23][24][25].Therefore, the detection of soft ISCr, which has a large magnitude of the R ISCr , is more necessary than the detection of hard ISCr with a small magnitude of the R ISCr for user safety to prevent the lithium-ion battery from causing hazardous events such as the thermal runaway.
Recently, for these reasons, methods for detecting the ISCr have been suggested [26,27].When the ISCr happens in the battery cell, terminal voltage of the cell decreases rapidly and then temperature of the cell increases drastically; these two characteristics can be used as pre-determined thresholds for detecting the ISCr [26].However, to obtain the threshold values, prior ISCr-test with the batteries were required.Besides the threshold-based method, a detection algorithm based on an equivalent circuit model of the battery has been introduced [27].Using variation of estimated parameters in both the equivalent circuit model and the energy balance equation, the method for detecting the ISCr was verified with large magnitude of the R ISCr , which was larger than 10 Ω [27]; the 10 Ω is considered to be the minimum value which must be detected early before occurrence of the thermal runaway [23,27].However, this algorithm was verified with only one type of current profile and the similar parameter variation may not be obtained depending on various other current profiles due to absence of the R ISCr in the equivalent circuit model.
The R ISCr is regarded as a promising fault index for detecting the ISCr because it can describe self-discharge phenomenon caused by the ISCr in the equivalent circuit model with R ISCr [28,29] and represent heat generated by the ISCr [30].Therefore, the ISCr detection methods, which directly use the R ISCr as the fault index, have been introduced [29,[31][32][33].With measurement data of the self-discharge current and the terminal voltage of the battery with ISCr in the particular experiment configuration, the R ISCr can be calculated correctly [29,31].However, these experimental methods cannot be used as on-board ISCr detection when the restricted experiment system is not configured in actual application environment.Therefore, after analyzing the self-discharge phenomenon caused by the ISCr, the equation for calculating the R ISCr was derived and then the R ISCr (1∼20 Ω) was estimated to detect the ISCr [32].However, the accuracy of the R ISCr estimated from two different load current profiles was low, because the R ISCr in the equivalent circuit model was not used to estimate open circuit voltage (OCV) of the battery with ISCr.To overcome this error, the previously estimated R ISCr in the model was used to update the model iteratively and to estimate the OCV, and then the next estimated R ISCr (5∼50 Ω) became accurate enough to detect the ISCr [33].
For a large capacity and a high power, a lithium-ion battery pack, where many battery cells are connected, is used in actual applications such as electric vehicles and energy storage system for the grid [34,35].Studies for detecting the ISCr in the battery pack have been suggested [36][37][38].The ISCr of the battery pack was detected based on the correlation coefficient of terminal voltages of individual cells [36], and the deviations of both state of charge (SOC) and heat generation power were used to detect the ISCr in the pack [37].However, these two methods were verified with only the hard ISCr having a magnitude of R ISCr (0.36 Ω for [36] and 0.35∼2.4Ω for [37], respectively), which is small enough to cause the dangerous incidents such as the thermal runaway in the cell [30].Especially, temperature data of the individual cells were needed to detect the hard ISCr correctly [37].Magnitude, differential value and fluctuation of estimated model parameters were acquired to detect the ISCr with a wide range of R ISCr (1 ∼100 Ω) after estimating the parameters in the mean-difference model of the battery pack [38].These three detection methods for the battery pack have a common constraint: all terminal voltages of the individual cells in the battery pack must be provided.
If many battery cells are connected in series in the battery pack, many channels with high accuracy and high sample rate for measuring all the individual cell voltages increase the cost of battery management system (BMS).In addition, a data storage unit, needed to save and monitor the measurement data, can lead to increasing the cost of BMS [39].Furthermore, depending on the applications of the battery pack and the BMS with various topologies, the individual cell voltages may not be provided with high precision and high sample rate and may not be saved due to the cost of the BMS [40,41].When the data of individual cell voltages are not provided from the BMS, the ISCr detection methods for the pack [36][37][38] have trouble in deriving the properties of the faulted battery from the battery pack, resulting in problem of detecting the ISCr of the cell in the pack.Therefore, an algorithm for detecting the soft ISCr with load currents and whole terminal voltages of the battery pack is necessary definitely.Moreover, considering data computation, using the whole terminal voltage for diagnosing the ISCr in the battery pack is more efficient than using all individual cell voltages.
This paper proposes a method for detecting the soft ISCr in the lithium-ion battery pack, where normal batteries and a battery with ISCr are connected in series, with the load currents and the whole terminal voltages of the pack.To reflect the effect of ISCr in the battery pack clearly, the OCV of the faulted cell is extracted from the pack with ISCr.The proposed algorithm estimates the SOC of the battery pack with the extended Kalman filter (EKF) to increase accuracy of the pack SOC.Then, the SOC of the normal cell is obtained with the Coulomb counting method and a stable initial value, which is determined from the estimated SOC of the pack.Using the SOC estimates of the pack and the normal cell, the OCV of the battery cell with ISCr in the pack can be obtained; this shows the self-discharge phenomenon caused by the ISCr obviously.As a result, the R ISCr can be estimated accurately, and the soft ISCr in the battery pack can be detected.To verify the proposed method, various soft ISCr cases were configured for simulation and experiment, and two load current profiles: dynamic stress test (DST) and urban dynamometer driving schedules (UDDS) were used.
The remainder of this paper is organized as follows: the proposed algorithm is carried out in Section 2; the configurations of simulation and experiment are introduced in Section 3; the results of the proposed algorithm and the discussions are presented in Section 4; the conclusions of this study and the outline of future work are summarized in Section 5.

Overall Scheme for ISCr Detection Algorithm
To detect the soft ISCr in the battery pack, the proposed method estimates the R ISCr of the faulted cell using the whole terminal voltages and the load currents of the pack.As a fault index, the estimated R ISCr can directly inform the user of degree of the ISCr in the pack.If the soft ISCr having a large magnitude of the R ISCr is detected, the BMS can give enough time to cope with the ISCr fault.The overall scheme of the proposed algorithm is depicted in Figure 1, which comprises of four estimation steps: estimating pack SOC (SOC p ); normal cell SOC (SOC n ); faulted cell SOC (SOC f ) and R ISCr .When the ISCr occurs in a single battery cell, abnormal properties such as the decrease in terminal voltage and the increase in battery temperature are easily observed [42,43].Hence, measured terminal voltage, current and temperature of the battery cell can be used to detect the ISCr.However, when an ISCr occurs in the battery pack, which only provides the load currents and the whole terminal voltages of the battery pack, the decrease in the whole terminal voltage caused by the ISCr in the battery pack is not obeserved conspicuously.Thus, additional decrease in the OCV f , induced by the self-discharge current in the faulted cell, should be extracted to ensure the high accuracy of R ISCr estimates.First of all, using the equivalent circuit model of the battery pack with ISCr and the EKF algorithm, the SOC p is estimated.As a second step, the SOC n can be calculated by the Coulomb counting method with a initial value of the SOC n which is obtained with the estimated SOC p .Then, using the estimated SOCs of both the battery pack and the normal cells, the OCV of the faulted cell OCV f and the SOC f are obtained at the third and fourth estimation steps.Once the difference between the initial and present SOC f estimates is more than a certain value, which is determined as 0.1 (10% of the total capacity of the cell) and discussed in Section 4.6, the R ISCr of the faulted cell in the battery pack can be estimated and used to detect the soft ISCr in the pack as the fault index.

Equivalent Circuit Model of Battery Pack with ISCr
In Figure 2, the battery pack consists of several normal battery cells and one faulted battery cell, and the whole terminal voltage V t,p and the load current I L are described, where m is the number of cells connected in series in the pack.The normal cell is represented by a simple equivalent circuit model [29] composed of the OCV (V OC ) and an internal resistance R. The sum of V OC s and Rs of normal cells described in Equations ( 1) and ( 2) is used to express the model of normal cells with discretization step in Equation ( 3), and the sum of the terminal voltages of the normal cells V t,n can be induced by Ohm's law in Equation ( 3) [44].The cell with ISCr is represented by the simple model with the R ISCr , which is connected with the model of the normal cell in parallel [29,38].The subscript f is used particularly in parameters associated with the faulted cell.The I L is divided into two currents which are the self-discharge current I 1L flowing through the R ISCr and residual current I 2L (Equation ( 4)), and V t, f is the terminal voltage of the faulted cell.The faulted cell model is described in Equations ( 5) and ( 6) [33,38].To represent the model of the battery pack with ISCr, the V t,p is obtained by adding the V t,n and the V t, f , as shown in Equation (7).

Estimation of Pack SOC Using EKF
The EKF algorithm is a common method to estimate accurate SOC because the estimates are not affected by measurement noise dominantly due to properties of the battery system reflected in the state space model [45][46][47].In this paper, the EKF algorithm is used to estimate the SOC p correctly.After assuming that R ISCr R f in the model of battery pack with ISCr (Equation ( 7)), the simplified V t,p can be expressed with battery pack OCV (V OC,p ) and whole internal resistance R p which is the sum of R n and R f in Equation ( 8).The recursive least squares (RLS) algorithm is used to identify the parameter of the normal battery model [48].Using the RLS algorithm, the R p is obtained from the V t,p and the I L of the normal battery pack.Then, the estimated R p is used to configure the state space model of battery pack with ISCr in the EKF algorithm.When the I L is positive during battery pack charging and negative during discharging, in recursive discrete-time form, the SOC p is calculated with the Coulomb counting method [44] in Equation ( 9), where η is the charging and discharging efficiency, ∆t is the sample period and C n is the nominal capacity of the normal battery pack.In this study, the charging and discharging efficiency is defined as 1.
To estimate the SOC p using the EKF algorithm, the corresponding equations are listed in Table 1, where x k is the SOC p , y k is the V t,p , u k is the I L , k is the sample index, and w k and v k are the zero mean Gaussian noise with covariance of Q and T.

Description Equation Step
State space model Initial assumed values

Estimation of Normal Cell SOC
After starting to estimate the SOC p , the imprecise SOC p can be obtained for several seconds because of the incorrectly assumed initial value of the SOC p .Thus, the estimated SOC p obtained from the stable point p st , where P k is lower than a certain small value, is used to estimate the SOC n , and the p st is the sample index.By assuming that the SOC p (p st ) equals to the SOC n (p st ) at the stable point, the next SOC n can be calculated with the I L and the Coulomb counting method described in Equation (9).

Estimation of OCV and SOC for Faulted Cell
The relationship between OCV and SOC of the normal battery pack (Figure 3) obtained from the prior test [49] is essential for conducting the proposed algorithm.The fully charged battery pack is rested for 3600 s to obtain the OCV, which is equal to the terminal voltage, at 100% SOC.Then the battery pack is discharged with 0.5C (1 A) for 720 s to set 90% SOC, where the C (C-rate) is defined as the charge and discharge current of the battery, and then rested for 3600 s to get the OCV at 90% SOC.By repeating the process, the relationship between OCV and SOC can be obtaind.Using the relationship, the OCV p and the OCV n can be acquired and are expressed in Equations ( 10) and ( 11), where h is the function representing the relationship.By subtracting OCV n s from the OCV p , the OCV f is calculated.In sequence, the SOC f is obtained from the inverse function h −1 of h (Equation ( 12)).In the process of using the h −1 to obtain the SOC f , the accuracy of SOC f is dependent on the slope from the relationship between OCV and SOC; i.e., in the range from 100% SOC to 50% SOC with the steep slope from the relationship a small SOC f error is caused by error of the V OC, f , while a large SOC f error can be induced by the same error of the V OC, f in the range from 50% SOC to 10% SOC with the gradual slope.To avoid this problem, the data of V t,p and I L were used to conduct the proposed algorithm until the estimated SOC f ≤ 0.55.

Calculation of R ISCr
Using the self-discharge phenomenon [32,33], the R ISCr can be estimated with the estimated SOC f .To formulate Equation ( 16) for calculating the R ISCr , the Coulomb counting method is used with respect to the SOC f .I L − I 1L is used in Equation ( 13), because the faulted cell is discharged from the I 2L instead of the I L .
To delete the unknown terms, SOC f (0) and I 1L , the k − 1th equation in Equation ( 13) is subtracted from the kth equation in Equation ( 13), and then the I 1L is replaced with . In Equation ( 14), is used instead of the V t, f because the V t, f is not provided measurement data from the battery pack with ISCr.
The estimated R ISCr from Equation ( 14) using a short interval between k and k − 1 is vulnerable to errors of the estimated SOC f because the slight variation of the estimated SOC f in a short interval cannot reflect dominant self-discharge phenomenon from the ISCr.Therefore, the interval must be extended by adding the k − 1, k − 2, k − 3, • • • , p st + 1th equations in Equation ( 14) to the kth equations in Equation ( 14) (p st + 1 < k); i.e., the R ISCr is estimated once the difference between the initial estimated SOC (SOC f (p st )) and present estimated SOC (SOC f (k)) is more than 0.1, which is 10% of the total capacity of the cell.
The self-discharge phenomenon from the ISCr can be explained with the last term of Equation ( 15).For the normal battery cell with R ISCr ∼ = ∞, the last term can be approximated to zero.However, when the ISCr occurs in the cell, the non-zero last term represents an additional decline in the SOC of the faulted cell due to the self-discharge current.Consequently, the estimated R ISCr ( R ISCr ) can be obtained with Equation ( 16) after Equation ( 15) is rearranged.

Parallel Processing of Proposed Algorithm
Once the SOC f (p st ) is obtained, the SOC f (p st ) is used in Equation ( 16) to calculate the R ISCr continuously.If error exists in the SOC f (p st ), the error can affect the accuracy of R ISCr .Therefore, it is necessary to estimate R ISCr s from the various p st positions and apply these estimated R ISCr s to the ISCr fault index.Based on the estimates of SOC f (p st ), if the difference between SOC f (p st ) and SOC f (k) is more than or equal to 0.01, 0.02, 0.03, • • • , the proposed methods are carried out sequentially and implemented in parallel with previously executed method to diversify the stable point of estimated SOC of the faulted cell.In this paper, four proposed methods were executed additionally in parallel.

Simulation Configuration
The simulation model of the battery pack with ISCr was configured with MATLAB/Simulink (MATLAB R2017b, MathWorks, Natick, MA, USA) [50].In both simulation and experiment configuration, the battery pack was composed of four normal battery cells and a faulted battery cell with ISCr.A first-order RC model [51], where the RC represents a parallel resistor-capacitor sub-circuit, was utilized to describe the normal cell in the battery pack.The parameters for the first-order RC model were estimated with the RLS algorithm [48] and the experimental data of terminal voltage and load current, measured when the normal cell was discharged by load current profiles.The normal cell is same with the cell in the battery pack of the experimental environment.In this study, two load current profiles such as the DST and the UDDS were used to verify the proposed method in simulation and experiment.Prior characteristic tests for obtaining the capacity and the relationship between OCV and SOC of the normal cell were conducted, and these two data were also used for the simulation configuration.The simulation model of the ISCr-faulted cell was represented by connecting the R ISCr with the normal model in parallel.In this study, various resistance values such as 50 Ω, 30 Ω, 20 Ω, 10 Ω and 5 Ω were used to represent diverse ISCr fault conditions.

Experimental Configuration
Figure 4 shows the experimental set-up for the ISCr tests with the battery packs.Two identical battery packs, A and B, were used to acquire the experimental data.The battery packs were tested in a thermal chamber, and the ambient temperature was maintained at 25 ± 1 • C. The important specifications of the cell are shown in Table 2.The prior tests for battery packs were conducted to obtain the capacities and the relationship between OCV and SOC.Based on the nominal capacity in Table 2, for capacity test, the packs were charged with the constant-current constant-voltage (CC-CV) protocol.For all experiments including the prior tests and the ISCr fault tests, when the batteries were charged, the CC-CV protocol was used.Charge-current was 0.5C (1 A) with upper cutoff-voltage 4.2 V in the CC mode, and cutoff-current was 0.05C (0.1 A) in the CV mode.Then, the packs were discharged with 0.5C as CC discharging [52].The discharged capacities were regarded as true capacities.The measured true capacities of the battery packs, 2.1974 Ah and 2.1949 Ah for pack A and B, respectively, were used for acquiring the correct relationship between OCV and SOC.To check the distribution of capacities of cells in the pack, the capacity test was conducted for individual cells.As a result, the mean and standard deviation are 2.1944 Ah and 0.0062, respectively, for pack A; and 2.1959 Ah and 0.0071, respectively, for pack B. The experiments of the battery packs were conducted using a battery test device (Regenerative Battery Pack Test System 17020, Chroma, Taoyuan, Taiwan) with the sample period of 0.1 s.The five 10 Ω resistances, which have ±5% tolerance, were used to make various true R ISCr s.Their measured true values were 49.91 Ω, 29.98 Ω, 20.00 Ω, 10.02 Ω and 5.00 Ω, and were used to calculate relative errors of the fault index.These resistances were connected with one of the cells in the pack in parallel to represent the ISCr, and a switch was used to initiate the ISCr faults in the packs when the load current profiles were applied to the packs.For the ISCr experiments, the initial SOC of the pack was set to 90% after the pack was charged and rested for 3600 s.Subsequently, the DST current profile was used to discharge the pack A until its SOC reaches 10% of its total capacity to prevent the pack from being over-discharged, while the UDDS current profile was used for the pack B.

Estimation Results of SOCs for Pack and Faulted Cell
As illustrated in Figure 6a, the initial values of SOC p estimates were 0.5 due to the initially assumed value of the x 0 , and the estimated SOC p s of the battery pack with ISCr decreased faster than that of the normal battery pack.However, the difference between the normal SOC p and others with various ISCr fault cases did not reflect significantly the effect of self-discharge from the ISCr in the battery pack, consequently leading to inaccurate estimates of the R ISCr .Thus, deriving the properties of the faulted cell from the pack with ISCr was necessary to observe the effect of ISCr clearly.The estimated SOC f is described in Figure 6b with the specific SOC f range from 90% SOC to 55% SOC.There were no estimated values of SOC f before the p st was reached, because the SOC n was calculated after the p st was reached.Please note that by extracting the V OC, f with these V OC,n and V OC,p from the pack with ISCr, the effect of the self-discharge phenomenon, which is not dominantly observed in the SOC p (Figure 6a), becomes enlarged noticeably in the SOC f (Figure 6b).

Estimated R ISCr s from Parallel Processing
When the battery pack A was connected with true R ISCr (49.91 Ω) and discharged with the DST current profile, the R ISCr s were estimated from different four stable points of estimated SOC of the faulted cell (SOC f (p st )) using the parallel processing (Figure 7).Although the four proposed methods were executed additionally in parallel with the firstly implemented method to diversify the SOC f (p st )s, the number of SOC f (p st )s can be different in accordance with the ISCr fault conditions in the pack; i.e., even though the four proposed algorithms are added sequentially in all ISCr fault cases, new SOC f (p st )s may not be extracted because of the condition, where the SOC f (p st ) was obtained if the P k was lower than the certain small value (1.4 × 10 −6 ) in Section 2.4.In the case of the experiment with the true R ISCr (49.91 Ω), the 1st stable point were obtained from the firstly executed method, while the three stable points were extracted from the the four added algorithms. Figure 7 shows the slightly different R ISCr s for the true R ISCr (49.91 Ω).Although the R ISCr,4 obtained from the 4th stable point was most accurate among them in this case, the order of accuracy of R ISCr s from the diverse stable points was changed depending on various ISCr fault cases.Therefore, to reflect all R ISCr s in an ISCr fault condition, mean value (R ISCr ) of them was used as the fault index in simulation and experiment results.

Estimation Results of R ISCr in Simulation
Figure 8a shows the R ISCr s with different ISCr fault conditions when the DST current profile was applied to the pack with ISCr, and Figure 8b depicts the R ISCr s obtained from the UDDS current profile.The R ISCr s in the specific SOC f range from 90% SOC to 55% SOC are described in Figure 8. Depending on the various ISCr fault cases, the SOC f (p st ) and the 1st stable point, where the R ISCr,1 starts to be obtained, were different.Thus, the different start points for estimating the R ISCr s are aligned to zero in Figure 8, and the time of x-axis represents the total time used to obtain the R ISCr with the R ISCr s.Estimation results from two different current profiles were similar, and the values of R ISCr slightly fluctuated because they were affected by waveforms of the current profiles.
The R ISCr in Equation ( 16) is estimated with difference between the SOC f (p st ) and the SOC f (k), which reflects the self-discharge phenomenon caused by the ISCr.At the same time, to discharge the pack, the difference between the SOC f (p st ) and the SOC f (k) becomes large as the magnitude of R ISCr is small (Figure 6b).In addition, the large difference can be obtained from long discharge time.The large difference between the SOC f (p st ) and the SOC f (k) is insensitive to errors of estimated SOC f , leading to accurate estimates of the R ISCr .In cases of ISCr 50 Ω, 30 Ω and 20 Ω, compared to cases of ISCr 10 Ω and 5 Ω, large errors of the R ISCr , caused by errors of estimated SOC f , occur in the early stage for estimating the R ISCr because of subtle difference between the SOC f (p st ) and the SOC f (k) in short discharge time.However, the R ISCr s in cases of ISCr 50 Ω, 30 Ω and 20 Ω gradually approach true R ISCr s because the difference between the SOC f (p st ) and the SOC f (k) becomes large gradually as the pack is discharged, resulting in the dominant effect of the self-discharge phenomenon compared with that of the estimation errors of SOC f .It is also the reason that the R ISCr is used as the fault index, which was calculated with the R ISCr s for the different stable points in the parallel processing considering errors caused by the position of SOC f (p st ).Although the large difference was obtained for each cases ISCr 10 Ω and 5 Ω, the R ISCr still had errors, because using the V t,p m instead of the V t, f in Equation ( 16) greatly affected the R ISCr estimates compared to the cases of ISCr 50 Ω, 30 Ω and 20 Ω; i.e., while the V t,p m was similar with the V t, f for ISCr 50 Ω, 30 Ω and 20 Ω, the difference between the V t,p m and the V t, f became large in cases of ISCr 10 Ω and 5 Ω, leading to large errors of the R ISCr .
To evaluate the accuracy of the R ISCr , relative errors of the estimates were calculated with Equation (17).The relative errors of the final values in R ISCr s (final relative error) with the various ISCr fault conditions are shown in Table 3.It should be noted that the final relative errors are less than or equal to 10% except for the ISCr 5 Ω case.Even though the final relative errors of ISCr 5 Ω in the DST and the UDDS were about 26%, the ISCr fault (5 Ω) can be detected with the R ISCr .In addition, although the magnitude of R ISCr (10 Ω) is regarded as the minimum value which must be detected early [23,27], the R ISCr (5 Ω) in the cell cannot sufficiently increase the temperature of the battery which is too low to cause the thermal runaway [30,37].The reason that the magnitude of the R ISCr is 5 Ω and more for verification is to show that the proposed algorithm can be applied to detection of the soft ISCr.If the proposed method detects various soft ISCr faults, we can conclude that there is sufficient time to provide against the thermal runaway in the battery pack.

Estimation Results of R ISCr in Experiment
Figure 8c,d show the estimation results of the R ISCr in the experiment with the DST and the UDDS current profiles.The tendency of the obtained R ISCr in experimental results was similar to that of the simulation due to the reasons explained in Section 4.4.In addition, the simplified model of the battery pack with ISCr, which was induced by assuming R ISCr R f (Equation ( 8)), was validated, because all estimated R ISCr s were indeed much greater than the R f in both simulation and experiment.
In particular, the main difference between simulation results and experimental results was that the relative errors in the experiment increased (Table 4) because in all ISCr faults cases, the R ISCr s in the experiment were generally under-estimated compared to those of the simulation.Contrary to the configuration of the simulation, the characteristics of individual cells in the experiment, such as capacity, internal resistance, and relationship between OCV and SOC, were not identical.Due to both the model difference and measurement noise, the errors of SOC p estimates of the experiment increased compared to that of the simulation, leading to large errors of the estimated SOC f in the experiment.This large errors of SOC f estimates increased the values of the denominator in Equation ( 16); this was main cause to incur the under-estimation of R ISCr in the experiment.Meanwhile, in case of ISCr 5 Ω, the final relative errors in the experiments decreased compared to the results of simulation, because the R ISCr s for ISCr 5 Ω with large relative errors in the simulation became close to the true value of R ISCr by the under-estimation in the experiment.Although the relative errors of the experiment increased compared to those of the simulation, the obtained R ISCr was accurate enough to be used as the fault index to detect the soft ISCr before the thermal runaway occurs in the battery pack, and to classify the various ISCr fault conditions.

Other Discussions
For the normal battery pack, the self-discharge current I 1L cannot flow through R ISCr and the faulted cell is discharged by the load current I L .Therefore, the SOC f is represented as Equation (18) with the I L , and the equation for obtaining the R ISCr for the normal battery pack ( R ISCr,n ) can simply be expressed with both the measured V t,p and the (Equation ( 19)) which is the difference value between the estimated SOC f and the true SOC f .
When the SOC f is estimated with the DST load currents and the whole terminal voltages of the normal battery pack in the specific region, the maximum value and relative error of the were 0.0197 and 3.6%, respectively, for the simulation; and were 0.0502 and 8.4%, respectively, for the experiment.Actually, in Section 4.5 as the reason for the under-estimation of R ISCr in the experiment, it was checked that the error of estimated SOC f for the experiment was larger than it for the simulation.These two maximum errors were used to calculate the R ISCr,n for the normal cases in both simulation and experiment.Due to the non-zero value of , the obtained R ISCr,n monotonically increased.In Figure 9, the dotted lines show R ISCr,n s of the normal battery pack in the specific region and represent reliable maximum values of R ISCr s for all ISCr fault cases.If the obtained R ISCr exists in the region above the dotted line, the R ISCr are determined as unreliable estimation values, while the region under the dotted line is defined as the reliable estimation region of the R ISCr .From the stable point p st , sufficient time to obtain reliable R ISCr is necessary.If the R ISCr is obtained when the k is very close to the p st in Equation ( 16), the R ISCr can be located in the region above the dotted line, because the values of the dotted line are small in the early stage shown in Figure 9 and the obtained R ISCr for ISCr fault cases can be affected by both the estimation errors of SOC f and measurement noise in the V t,p and the I L .Thus, in the proposed method, the R ISCr was estimated once the difference between the SOC f (p st ) and the SOC f (k) was more than or equal to 0.1; this condition guaranteed that the R ISCr s in various ISCr faults of simulation and experiment exist in the reliable estimation region.
When a battery pack is manufactured with used lithium-ion cells, a variation in the characteristics of individual cells in the pack becomes large [41], resulting in large error of R ISCr estimates.However, the battery pack made with used cells can be operated with balanced voltage and SOC due to the proper screening process in configuring the battery pack [53].Therefore, the proposed method can be applied in both fresh and reused battery pack.
If an ambient temperature in operation environment of the battery varies and the battery model which does not reflect thermal properties of the battery is used to estimate the SOC, errors of estimated SOC become large [54][55][56][57], leading to large error of R ISCr estimates.Although in this study the proposed method focused on detecting the soft ISCr at constant temperature, depending on real applications the ambient temperature can be changed [58].Therefore, detection of ISCr in the battery under varing ambient temperature is a maningful and interesting subject of reaserch as a future work.

Conclusions
In this paper, a method for detecting the ISCr early in the lithium-ion battery pack was introduced.The battery pack with ISCr was represented with the equivalent circuit model with the R ISCr and the EKF algorithm was used to estimate SOC p accurately.The OCV of the faulted cell was derived from the battery pack with ISCr to reflect the self-discharge phenomenon caused by the ISCr in the battery pack clearly, because the effect of ISCr in the battery pack was not observed in both the V t,p and the SOC p obviously.Using the Coulomb counting method and the stable initial value of the SOC n , obtained from the estimated SOC p , the SOC n was calculated.The OCV p and the OCV n were acquired from the relationship between OCV and SOC of the normal battery pack, and then the OCV f was calculated with these two OCV values.Subsequently, the R ISCr (5∼50 Ω) of the battery pack with ISCr was estimated accurately using the self-discharge phenomenon in the SOC f .The proposed algorithm was verified for various soft ISCr fault conditions such as diverse magnitudes of true R ISCr and two load current profiles in both the simulation and the experiment.In addition, through estimating the R ISCr from the normal battery pack and analyzing it, it was checked that estimated R ISCr s in the various fault cases were reliable.Using the proposed algorithm, the R ISCr was estimated with high accuracy, and the soft ISCr in the battery pack can be detected using the R ISCr as the fault index.Our future research will focus on increasing the accuracy of the R ISCr estimates and extending the availability of our proposed algorithm to both the aged battery pack and the battery under varing ambient temperature.

Figure 1 .
Figure 1.The scheme of the proposed algorithm.

Figure 2 .
Figure 2. Equivalent circuit model of the lithium-ion battery pack with internal short circuit (ISCr).

Figure 3 .
Figure 3. Relationship between open circuit voltage (OCV) and state of charge (SOC) of the normal battery pack.

Figure 5 Figure 5 .
Figure 5 depicts the terminal voltages of battery pack A and B in the experiments depending on the magnitudes of R ISCr .The terminal voltages of the packs were measured until the pack SOC reached 10%.As the magnitude of R ISCr is small, the terminal voltages decreased rapidly compared with the voltages of the normal battery pack, leading to rapid termination of experiments for the battery packs with ISCr.The additional decline in terminal voltages caused by the self-discharge phenomenon of the ISCr was not observed clearly.In addition, because the terminal voltages of the packs were affected by waveforms of load current profiles in common, the terminal voltages fluctuated and the slight difference of voltages between the normal pack and the pack with ISCr was not monitored readily.Therefore, it was difficult to detect the soft ISCr directly with only the measurment data of terminal voltages of the packs.

Figure 6 .
Figure 6.Estimated SOCs with different R ISCr s in experiment with DST current profile: (a) battery pack and (b) faulted battery cell.

Figure 7 .
Figure 7.Estimated R ISCr s from different stable points in experiment with true R ISCr 49.91 Ω and DST current profile.

Figure 8 .
Figure 8.Estimated R ISCr s from the various ISCr fault cases in simulations: (a) DST current profile and (b) UDDS current profile and in experiments: (c) DST current profile and (d) UDDS current profile.

Figure 9 .
Figure 9. Reliable maximum R ISCr,n in simulation and experiment with normal battery pack and DST current profile: (a) simulation and (b) experiment.

Table 3 .
Final relative errors (%) of R ISCr in simulation depending on the ISCr faults.

Table 4 .
Final relative errors (%) of R ISCr in experiment depending on the ISCr faults.