The stateofcharge (SOC) and stateofhealth (SOH) estimation of batteries play important roles in managing batteries for automotive applications. However, an accurate state estimation of a battery is difficult to achieve because of certain factors, such as measurement noise, highly nonlinear characteristics, strong hysteresis phenomenon, and diffusion effect of batteries. In certain vehicular applications, such as idle stop–start systems (ISSs), significant errors in SOC/SOH estimation may lead to a failure in restarting a combustion engine after the shutoff period of the engine when the vehicle is at rest, such as at a traffic light. In this paper, a dual extended Kalman filter algorithm with a dynamic equivalent circuit model of a lead–acid battery is proposed to deal with this problem. The proposed algorithm adopts a battery model by taking into account the hysteresis phenomenon, diffusion effect, and parameter variations for accurate state estimations of the battery. The validity of the proposed algorithm is verified through experiments by using an absorbed glass mat valveregulated lead–acid battery and a battery sensor cable for commercial ISS vehicles.
I. INTRODUCTION
In most countries, global warming has become a significant issue that requires urgent action in terms of decreasing greenhouse gas emissions, particularly in the automotive sector. Electric vehicles (EVs), such as battery EVs, hybrid EVs (HEVs), and plugin HEVs, are promising solutions that can considerably lower the amount of greenhouse gas emissions and achieve better fuel efficiency than conventional vehicles. Electric cars will undoubtedly play an increasingly large role in many countries in the decades ahead as energy independence and environmental concerns intensify. However, EVs are estimated to gain only a modest ground up to 2020
[1]
. Gasoline and dieselpowered vehicles are improving faster than expected and will continue to dominate the global landscape. In addition, carmakers are trying to meet future emission targets mainly through improvements in internal combustion engines (ICEs). Idle start–stop (ISS) systems are becoming increasingly prevalent in ICEs because of their ability to reduce emissions and fuel consumption in a costeffective manner. An ISS system is an easy and lowcost technology that enables an ICE to power off automatically when the car is stopped and to restart upon demand. As a result, the ICE undergoes far more starting events, and the behavior of the ICE during startup becomes more critical. In this kind of vehicle, valveregulated lead–acid (VRLA) batteries are often employed because of their maintenancefree characteristics and higher resistance to shock and vibrations. Batteries operate under intense conditions because of their frequent charge and discharge cycles and cranking; therefore, accurately estimating the state of charge (SOC) and state of health (SOH) of batteries is crucial to avoid failures of restarting the engine after being shut down. Moreover, battery state information is essential to enhance the efficiency of battery energy utilization, lengthen the lifetime of the battery, and prevent permanent damage to the battery.
Several methods of estimating the SOC and SOH of a battery have been proposed by using coulomb counting, artificial neutron networks (ANNs), fuzzy logic (FL), and extended Kalman filters (EKFs)
[2]

[10]
. The coulomb counting method can be simply implemented by integrating the battery current over time
[11]
. However, the reliability of SOC estimation cannot be maintained when errors accumulate over the operation time and when the initial value of the SOC is wrong. The ANNs and FL approaches can estimate the SOC of a battery with an arbitrary initial SOC value. However, these methods are relatively expensive to implement and require the training data of the battery
[12]
. The EKF is a popularly accepted tool that provides a theoretically welldesigned and timeproven method to filter the measurements of the system input and output to produce an intelligent estimation of the dynamic state of a system
[13]

[15]
. The dual EKF (DEKF) algorithm is a combination of two EKFs, in which the SOC is estimated by the first EKF, and the capacity is estimated by the second EKF
[16]
. The main advantage of this method is that it can provide reliable estimation results while the battery is in operation. As a result, the DEKF algorithm is preferred among the abovementioned methods. However, because the estimation accuracy depends mainly on the adoption of a comprehensive and wellparameterized battery model, robust SOC/SOH estimation is possible only when the model is accurate enough.
Although precise modeling of a battery plays an important role in estimating SOC/SOH accurately, it is a complicated and challenging task because VRLA batteries show nonlinear characteristics, and various factors should be considered. For example, the parameters of the equivalent circuit model of a battery depend mainly on operating conditions, such as the SOC of the battery, charge/discharge current, and temperature
[17]
. The hysteresis phenomenon is a historydependent characteristic caused by internal chemical processes that cause difficulties in estimating the SOC, particularly when the battery state turns over from a partial charge or discharge
[18]
,
[19]
. In addition, the diffusion effect, which has a relatively long time constant, has to be considered to improve the estimation accuracy.
In this paper, an accurate SOC/SOH estimation method using DEKF for absorbed glass mat (AGM) type VRLA batteries is proposed. A practical battery model that considers hysteresis and diffusion phenomenon is introduced. The hysteresis effect is modeled by calculating the normalized integration of the charge throughput during a partial cycle, and the diffusion phenomenon is considered by reconstructing the open circuit voltage (OCV) decay depending on the rest time. The size of the state matrix calculation for the Kalman filter operation can be significantly reduced by incorporating the diffusion effect into the OCV. In addition, the dependency of the battery parameters on various internal and external conditions, such as SOC, temperature, and charge/discharge current and their combinations, are considered based on pretest results from the battery. All the test methods for modeling hysteresis, diffusion, and parameter variations are detailed in the following sections. The proposed algorithm is implemented in a battery sensor cable (BSC) with an ARM7 microprocessor and validated through experimental results.
II. PRETEST METHODS FOR MODELING AN AGM VRLA BATTERY
 A. Popular Model of a Lead–Acid Battery
Fig. 1
shows a popular model of a lead–acid battery that consists of a voltage source, a series resistor, and two sets of parallel resistor–capacitor circuits. In the figure, E
_{eq}
represents the equilibrium voltage of the battery, and R
_{i}
represents the resistance of the contacts, the inter cell connections, and the electrolyte. R
_{ct}
and C
_{dl}
are the chargetransfer resistor and electric doublelayer capacitor, respectively, which represent the charge transfer reaction; R
_{d}
and C
_{d}
are the diffusion resistor and diffusion capacitor, respectively, which represent the diffusion phenomenon caused by the grade of the concentration of the electrolyte near the electrode
[20]
. The phenomenon causes a second overvoltage on the electrode potential called “diffusion overvoltage.” This model is suitable for describing the dynamic behavior of a VRLA battery. However, the model is unsuitable for a Kalman filter operation because using the model would lead to computational overkill.
Popular model of a lead–acid battery.
A lithium–ion battery, a simple method that uses a single average OCV–SOC relationship, can be employed because the hysteresis phenomenon is negligible. Hence, SOC estimation can be achieved with high precision without considering the hysteresis effect
[22]

[24]
. However, modeling an AGM VRLA battery is not feasible without considering the hysteresis phenomenon because a large error can result from estimating the OCV of a battery during a dynamic operation. In this study, the hysteresis is modeled by approximating it with a parallelogram. Then, the hysteresis factor is introduced to track the voltage variations caused by hysteresis depending on the charge throughput.
All the parameters of the model vary sensitively depending on the SOC, temperature, and charge/discharge current, as previously mentioned. Thus, the parameter variations have to be modeled carefully. To model the parameter variations, comprehensive pretests have to be conducted by using current pulse tests. All the pretests required to model an AGM VRLA battery are described in the following sections.
 B. Pretests to Develop the Proposed AGM VRLA Battery Model
In this section, a comprehensive explanation for the pretests required to model an AGM VRLA battery is described in detail. The following three kinds of pretests must be conducted to model the battery:

1) Full cycle test to derive the relationship between the SOC and the OCV and to extract the diffusion overvoltage:

2) Partial cycle test to model the hysteresis loop; and

3) Variable current pulse test to model the parameter variations.
An AGM VRLA battery (Solite AGM70L – DIN, 12 V, 70 Ah) is used for the tests, and all the tests are conducted in a constant temperature chamber. The battery is connected to a bipolar DC power supply (NF BP4610), and a program created in LabView 11.0 is used to perform the tests automatically. The charge/discharge cycle profile created in LabView is transmitted to the external control port of the bipolar DC power supply via the D/A converter of the data acquisition board (NI PCI DAQ 6154), and the bipolar power supply charges and discharges the battery exactly as commanded. A sensing circuit for the voltage and current is used to measure the terminal voltage and the current of the battery and the data acquisition board record them.
1) Full Cycle Test to Derive the Relationship between the SOC and the OCV:
To derive the OCV–SOC relationship, a charge/discharge test that uses variable current pulses is performed, as shown in
Fig. 2
. The magnitude and duration of the current pulse are C/20 and 1 h, respectively. Hence, each pulse is equivalent to 5% of the actual capacity. Twenty of these discharge current pulses are applied to the fully charged battery to fully discharge the battery. Then, the battery is charged again with the same current pulse until it reaches 60% of the SOC. Then, the duration of the current pulse is reduced to 0.5 h, and the current pulse is applied to the battery until it reaches 90% of the SOC. At this point, the battery is fully charged by a constant voltage charge until it reaches 100% of the SOC. The relaxation time between each pulse is 3 h, thereby ensuring that the measurement of the true OCV is possible.
Fig. 3
shows the measured voltage and current of the battery obtained through the discharge test at 25 ℃. The charge/discharge current and the terminal voltage of the battery are recorded at every second in a host computer through a data acquisition board. The full cycle test procedure is then repeated at 18, 0, and 45 ℃ The measured terminal voltages of the battery at 25 ℃ after charge and discharge at different SOC values are presented using threedimensional plots in
Fig. 4
. The OCV–SOC curve with a 3 min relaxation and a 3 h relaxation can be obtained as the projections of the battery voltage at 3 min and 3 h after the charge/discharge pulse current are removed. The OCV during charging is higher than that during discharging at the same SOC because of the hysteresis effect
[21]
.
Full cycle test procedure to obtain the OCV–SOC relationship of the AGM VRLA battery.
Measured current and voltage waveforms of the AGM VRLA battery during the discharge at 25 ℃.
OCV–SOC relationship at each SOC obtained from the full cycle test at 25 ℃.
On the basis of the full cycle test results, the measured OCV curves at 25 ℃ with a 3min relaxation and a 3 h relaxation and the difference between them are depicted in
Fig. 5
. Each of the OCV–SOC curves can be modeled by a fifthorder polynomial function and represented as follows:
(a) OCV–SOC curve with 3min and 3g relaxation at 25 ℃. (b) OCV difference between 3min and 3h relaxation curves.
The coefficients of the polynomial equations are listed in
Table I
. The measured curves at 18, 0, and 45℃ are also modeled with Equ. (1), and the OCV value at a certain temperature can be determined by using an interpolation technique that can construct new data based on the range of the obtained data of the OCV from the full cycle tests
[25]
.
COEFFICIENTS OF THE POLYNOMIAL EQUATIONS FOR EACH OCV–SOC CURVE AT 25 ℃
COEFFICIENTS OF THE POLYNOMIAL EQUATIONS FOR EACH OCV–SOC CURVE AT 25 ℃
The OCV lies between the OCV curve with the 3 min relaxation and the 3 h relaxation after the charge transfer overvoltage vanishes; thus, it can be modeled by using diffusion factor ζ, which can be calculated with time constant τ
_{d}
of the R
_{d}
–C
_{d}
circuit as
During the charge, the OCV can be reconstructed by using diffusion factor ζ as follows:
Thus, during the charging, the OCV can be calculated as Equ. (4) depending on the relaxation time, and the OCV can also be calculated in the same manner as follows:
and
As mentioned in Section II A, the diffusion overvoltage needs to be extracted and merged into the OCV to simplify the battery model.
Fig. 6
(a) shows the terminal voltage variation of the battery voltage after the C/20 charge current is removed at 65% of the SOC. The voltage decay can be divided into two parts by the charge transfer overvoltage and the diffusion overvoltage, except for the immediate voltage drop, because of resistance R
_{i}
. The charge transfer overvoltage is caused by R
_{ct}
C
_{dl}
in the battery model and is dominant in the transient voltage in the first 3 min and vanishes afterwards. Thus, the voltage decay after 3 min can be considered diffusion overvoltage that has a relatively long time constant. The terminal voltage of the battery can be represented by Equ. (6), and the parameters can be extracted simply by using the curve fitting
where y
_{0}
is the dc offset, t
_{rest}
is the relaxation time of the battery after charging or discharging, and A
_{1}
and A
_{2}
are the coefficients of the first and the second exponential terms, respectively. τ
_{1}
and τ
_{2}
represent the time constants of the R
_{ct}
–C
_{dl}
and R
_{d}
–C
_{d}
circuit, respectively. This test is repeated at each SOC, and the average value of τ
_{2}
is used to calculate diffusion factor ζ in Equ. (2).
Voltage variation after a charge current pulse at 65% of SOC at 25 ℃. (a) Measured voltage and fitted results with the model in Fig. 1. (b) Extracted charge transfer overvoltage and diffusion overvoltage.
In these results, the extracted charge transfer overvoltage and diffusion overvoltage can be drawn as
Fig. 6
(b), and the battery model can be simplified by merging the diffusion overvoltage into equilibrium voltage E
_{eq}
. E
_{eq}
is reconstructed to a voltage source, which varies not only with the SOC but also with the relaxation time, as illustrated in
Fig. 7
[26]
.
Proposed battery model with the reconstructed OCV.
2) Partial Cycle Test to Model the Hysteresis Loop:
Obtaining the voltage trajectory when hysteresis occurs is necessary because the hysteresis loop can be modeled by a parallelogram
[17]
,
[19]
,
[21]
. The partial cycle test is conducted to obtain information about the hysteresis loop. The current pulses illustrated in
Fig. 8
are used for the partial cycle tests. In this case, a current pulse equivalent to 4% of the actual capacity of the battery is used for the tests. At first, a fully charged battery is rested for 3 h. The battery is discharged continuously to 66% of the SOC and charged to 90% of the SOC by using six current pulses. Then, the battery is discharged again to 66% of the SOC with the same current pulses, as depicted in
Fig. 8
. This test is repeated at three different SOC ranges to derive the average value of the parameters to model the hysteresis with a parallelogram.
Partial cycle test procedure for the hysteresis loop modeling.
In these results, the hysteresis loop can be modeled with the partial cycle test results by a parallelogram, as shown in
Fig. 9
. The actual OCV during charging and discharging lies somewhere in between the 3 h relaxation charge and discharge curves depending on the charge throughput
[18]
.
(a) Modeling of the hysteresis loop with the partial cycle test results by a parallelogram. (b) Partial cycle tests performed at three different SOC ranges.
Thus, the transition of the OCV can be empirically reconstructed by introducing a hysteresis model. Hysteresis factor α is introduced to model the movement of the OCV between the lower boundary curve (3 h relaxation discharge curve OCV
_{d3h}
) and the upper boundary curve (3 h relaxation charge curve OCV
_{c3h}
) during a partial charge/discharge cycle. The OCV reconstruction considering the hysteresis effect can be performed by
Equ. (7) can be simplified as follows:
The empirical equation for hysteresis factor α corresponding to the charge throughput can be described by Equ. (9), as shown in
Fig. 9
(a):
Hysteresis factor α varies at a different rate according to the charge throughput; thus, the values of α
_{1}
and α
_{2}
must be calculated by Equs. (10) and (11), and updated with the charge throughput at each time index k
For instance, when a battery is being charged, α
_{1}
is calculated first until the charge throughput reaches 4%, and then α
_{2}
is calculated until the charge throughput reaches 24%. When the charge throughput exceeds 24%, the hysteresis effect disappears, and the OCV follows the boundary curve. Similarly, when the battery is being discharged, α
_{1}
is decreased to 0 first, and then α
_{2}
is decreased. The charge throughput can be calculated by
where C
_{b}
(Ah) is the capacity of the battery, and Δt is the sampling time. The
K
value, ΔSOC
_{1max}
, and ΔSOC
_{2max}
can be simply extracted by using least square fitting with the test data and a parallelogram, as shown in
Fig. 9
(a).
Fig. 9
(b) shows three hysteresis test results performed at different SOC values, and the results are used to calculate hysteresis factor α, which exhibits small differences at each SOC value.
Table II
shows the variation of the
K
values at three different SOC ranges for the partial cycle tests. The average value of
K
can be used to simplify the hysteresis voltage calculation because the values of
K
are almost the same regardless of the SOC range.
KVALUES OBTAINED AT THREE DIFFERENT SOC RANGES BY THE PARTIAL CYCLE TESTS
K VALUES OBTAINED AT THREE DIFFERENT SOC RANGES BY THE PARTIAL CYCLE TESTS
3) Variable Current Pulse Tests to Model Parameter Variations:
DEKF is a popular technique that can provide a powerful and intelligent solution for the online state estimation of a battery
[16]
. However, the accuracy of DEKF depends significantly on the accuracy of the parameter values in the battery model. As previously mentioned, battery parameters vary according to the operating conditions of the battery, such as the SOC, temperature, and charge/discharge current value
[27]
,
[28]
. Thus, the parameter variations of the battery model have to be modeled carefully and taken into account for the DEKF operation. To model the parameter variations of the battery, variable charge/discharge current pulses, which have different amplitudes (C/1, C/2, C/5, C/10, C/20, and C/40 [A]) and durations (1, 2, 5, 12, 24, and 51 [min]), are applied to the battery, as shown in
Fig. 10
.
Variable current pulse tests to model the parameter variation.
Tests are conducted at each SOC level (10% SOC step from 100% to 40%) to obtain information about the parameter variations. Thereafter, the tests are repeated at different temperatures (18, 0, 25, and 45 ℃) to investigate the temperature dependency of the battery parameters. After all the data are collected, the curve fitting technique is used to extract battery parameters, such as R
_{i}
, R
_{ct}
, and C
_{dl}
. As a result, the parameter variations caused by changes in the SOC, temperature, and charge/discharge current value can be modeled.
Fig. 11
shows the variation of charge transfer resistance R
_{ct}
at different temperatures. R
_{ct}
varies exponentially depending on the current values.
Fig. 12
shows the variation of R
_{ct}
depending on the charge/discharge current at 25 ℃, and R
_{ct}
is fitted by the exponential function as
Variation of R_{ct} according to the charge/discharge current at each temperature.
Variation of R_{ct} depending on the charge/discharge current at 25 ℃ fitted by the exponential functions.
The coefficients of Equ. (13), namely, m
_{Rct}
and n
_{Rct}
, are given in
Table III
.
COEFFICIENTS OF EQUATIONS (13), (14), AND (15)
COEFFICIENTS OF EQUATIONS (13), (14), AND (15)
Fig. 13
shows the variation of R
_{ct}
depending on the SOC of the battery. In this case, the normalized value of R
_{ct}
with respect to its value at SOC = 70% depending on the SOC is drawn and fitted by the secondorder polynomial equation in
Variation of the normalized R_{ct} value with respect to its value at 70% SOC depending on the SOC fitted by a secondorder polynomial function.
The coefficients of Equ. (14), namely, a
_{Rct}
, b
_{Rct}
, and c
_{Rct}
, are given in
Table III
.
Fig. 14
shows the variation of R
_{ct}
depending on temperature. In this case, the normalized value of R
_{ct}
with respect to its value at 25 ℃ depending on the temperature is drawn and fitted by the exponential function in
Variation of the normalized R_{ct} value with respect to its value at 25 ℃ fitted by an exponential function.
The coefficients of Equ. (15), namely, m’
_{Rct}
, and n’
_{Rct}
, are given in
Table III
.
As investigated in
[29]
, the parameters can be adjusted according to the SOC, temperature, and charge/discharge current value by using the following:
The parameters of the battery model in the DEKF can now be updated by using Equ. (16) on a realtime basis, thereby maintaining the accuracy of the battery model.
III. STATE ESTIMATION TECHNIQUE OF AN AGM VRLA BATTERY BY USING THE DEKF ALGORITHM
 A. DEKF Algorithm for SOC/SOH Estimation
The SOC indicates the relative level of charge presently held by a battery, whereas the SOH refers to the general decline in battery performance with respect to usage or aging. Therefore, estimating the SOC and SOH of a battery is important to predict the remaining useful energy and the remaining service life of the battery. In this paper, the DEFK approach is used to estimate the SOC and SOH of a battery. The DEFK process essentially combines two EKFs that run in parallel, where one is the state filter for estimating the SOC, and the other is the weighting filter for estimating the capacity
[16
,
30]
. At each time interval, the state filter uses the priori value of the weight filter and vice versa. The stepbystep calculation procedure is summarized as follows.
The nonlinear state–space model can be represented as
and
where
x
_{k}
is the state vector of the battery model,
θ
_{k}
is the time varying battery capacity,
u
_{k}
is the exogenous input,
y
_{k}
is the system output, and
w
_{k}
,
v
_{k}
,
r
_{k}
, and
e
_{k}
are independent Gaussian noise processes with covariance matrices
respectively. f(·,·,·) represents a nonlinear transition matrix function, and g(·,·,·) represents a nonlinear measurement matrix function.
The computing procedure for the DEFK can be summarized as follows
[16]
.
Step 1: Initializing at k=0.
where E[.] is the statistical expectation operator.
Step 2: Approximating the nonlinear functions.
Step 3: Updating the time for each filter.
State filter:
Weight filter:
Step 4: Updating the measurement for each filter.
State filter:
Weight filter:
 B. SOC/SOH Estimation using the DEKF Algorithm
The state space equation for the proposed battery model in
Fig. 7
in a discrete form can be represented as
and
The terminal voltage of the battery model can be represented by using a nonlinear function as follows:
The OCV value at each time index
k
can be reconstructed by using Equs. (4), (5), and (6)
Equation (34) can be rewritten in a compact form as
where
While the state equation in Equ. (36), as shown below, does not require an approximation, the nonlinear measurement function g(.,.) requires an approximation because the OCV–SOC relationship is nonlinear:
Thus, the Taylor series approximation described in the following equations is used for linearization:
and
where
The partial derivatives are computed at each time step, whereas the total derivatives are computed recursively with zero as initial values.
IV. EXPERIMENTAL VERIFICATION OF THE PROPOSED ALGORITHM
To verify the validity of the proposed estimation algorithm by using the DEKF, the proposed algorithm is tested with an AGM VRLA battery. A Solite AGM70LDIN (12 V, 70 Ah) battery is used to verify the proposed method. The algorithm is implemented by using C language in the ARM7 TDMI microprocessor integrated into the BSC for the ISS EV application, as shown in
Fig. 15
.
Implementation of the proposed algorithm in theARM7TDMI microprocessor of the battery sensor cable.
The battery is fully charged up to 100% of the SOC and then discharged to 50% of the SOC. Thereafter, the battery is charged and discharged repeatedly with different amounts of charge, and eventually, charging is stopped at 78% of the SOC. In this experiment, coulomb counting is considered the reference value of the estimation for the performance comparison.
The DEKF algorithm initially uses the coulomb counting method to estimate the SOC value. Then, the error between the measured terminal voltage and the estimated terminal voltage is used to correct the priori SOC value that is calculated by the coulomb counting method in a realtime update step. In this way, the DEKF can solve the problem of determining the initial SOC and eliminate the accumulated error caused by coulomb counting over longterm operation. Fifty percent of the SOC is used as an initial value instead of 100% of the SOC to verify that the DEKF converges to the true SOC value with a wrong initial SOC value. As shown in
Fig. 16
, the SOC estimation error at the beginning is 50%. After 1.5 h of operation, the SOC error is reduced to less than 15%, and afterward, it remains less than 5%. The error becomes larger when the battery changes its state from charged to discharged, and vice versa, because of the strong hysteresis effect.
SOC estimation results and errors of the proposed algorithm.
However, the error in the SOC estimation is still less than 5% over the entire test period. On the basis of these results, the proposed battery model, which includes the hysteresis, diffusion effect, and parameter variations, is good enough to reproduce the electric behavior accurately to provide reliable SOC estimation.
To guarantee the operation of the DEKF algorithm for estimating battery capacity, a pulse current cycle test that comprises a sequence of discharge pulses and relaxation, followed by a sequence of charge pulses and relaxation, is conducted. The terminal voltage and the current of the battery are recorded by the host computer and data acquisition board every second during the 450 h of the test. The proposed algorithm is also coded by LabView and tested with the recorded data. Capacity estimation results with a true initial value (0.98 Cn = 68.6 Ah) and wrong initial values, which are smaller and larger than the true value, are shown in
Fig. 17
. Five different initial capacity values (0.8, 0.9, 1.0, 1.1, and 1.2 Cn) are used in this test. As shown in
Fig. 17
, the real capacity of the battery can be estimated successfully by the proposed algorithm with less than 5% error after approximately 300 h in all cases. Results demonstrate the performance of the proposed algorithm in predicting the SOH of the battery; the proposed algorithm is determined to be suitable for ISS EV applications.
Capacity estimation results with the different initial capacity values.
V. CONCLUSION
In this paper, a novel technique for estimating the SOC/SOH of AGM VRLA batteries based on the DEKF algorithm has been proposed, and its performance has been verified by experiments with an AGM VRLA battery. The battery model is simplified by reconstructing the OCV, considering the diffusion and hysteresis effects of the battery for the light computation of the Kalman filter. The parameter variations of the battery model according to the SOC, current, and temperature have been modeled by using pretest results with the battery and applied to the model during operation in a realtime basis. The observed error of the SOC and SOH estimation is less than 5%. The proposed method can be used in any kind of battery management system for vehicle and energy storage applications.
BIO
Van Huan Duong was born in Bac Giang, Vietnam, in 1985. He received his B.S. and M.S. degrees in Electrical Engineering from Hanoi University of Science and Technology, Vietnam, in 2008 and Soongsil University, Korea, in 2013, respectively. He is currently working toward his Ph.D. degree in Electrical Engineering at the University of Wollongong, Australia. His current research interests include battery management system and microcontroller unit/digital signal processor applications.
Ngoc Tham Tran was born in Quang Nam, Vietnam, in 1987. He received his B.S. and M.S. degrees in Electrical Engineering from Danang University of Technology, Vietnam, in 2010 and Soongsil University, Republic of Korea in 2015, respectively. He is currently with the Department of Science and Technology, Danang University of Technology. His current research interests include battery management system for electric vehicles, such as battery modeling and stateofcharge and stateofhealth estimation.
Woojin Choi was born in Seoul, Republic of Korea, in 1967. He received his B.S. and M.S. degrees in Electrical Engineering from Soongsil University, Republic of Korea, in 1990 and 1995, respectively, and his Ph.D. degree, also in Electrical Engineering, from Texas A&M University, USA, in 2004. From 1995 to 1998, he was with Daewoo Heavy Industries as a research engineer. In 2005, he joined the School of Electrical Engineering, Soongsil University. His research interests include modeling and controlling electrochemical energy sources, such as fuel cells, batteries, and supercapacitors; power conditioning technologies in renewable energy systems; and DCDC converters for fuel cells and hybrid electric vehicles.
DaeWook Kim was born in Seoul, Republic of Korea, in 1973. He received his B.A. and Ph.D. degrees in Economics from Yonsei University, Republic of Korea, in 1999 and the University of California in Davis, USA, in 2004, respectively. From 2004 to 2007, he was with the Korea Institute for Industrial Economics and Trade as a research fellow. In 2007, he joined the Department of Economics, Soongsil University. His current interests include energy economics, particularly market structure and competition in energy industries.
Mosquet X.
,
Devineni M.
,
Mezger T.
,
Zablit H.
,
Dinger A.
,
Sticher G.
,
Gerrits M.
,
Russo M.
2011
“Poweing Autos to 2020, The Era of the Electric Car?,”
The Boston Consulting Group
Leksono E.
,
Haq I. N.
,
Iqbal M.
,
Soelami F. X. N.
,
Merthayasa I. N.
“State of charge (SoC) estimation on LiFePO4 battery module using Coulomb counting methods with modified Peukert,”
in Rural Information & Communication Technology and ElectricVehicle Technology (rICT & ICeVT), 2013 Joint International Conference on
2013
1 
4
Shahriari M.
,
Farrokhi M.
2013
“Online stateofhealth estimation of VRLA batteries using state of charge,”
IEEE Trans. Ind. Electron.
60
(1)
191 
202
DOI : 10.1109/TIE.2012.2186771
Lin H.T.
,
Liang T.J.
,
Chen S.M.
2013
“Estimation of battery state of health using probabilistic neural network,”
IEEE Trans. Ind. Informat.
9
(2)
679 
685
DOI : 10.1109/TII.2012.2222650
Feng R.
,
Zhao S.
,
Lu X.
“Online estimation of Dynamic StateofCharge for lead acid battery based on fuzzy logic,”
in Measurement, Information and Control (ICMIC), 2013 International Conference on
2013
447 
451
Zheng C.
,
Fu Y.
,
Mi C. C.
2013
“State of charge estimation of lithiumion batteries in electric drive vehicles using extended Kalman filtering,”
IEEE Trans. Veh. Technol.
62
(3)
1020 
1030
DOI : 10.1109/TVT.2012.2235474
Charkhgard M.
,
Farrokhi M.
2010
“Stateofcharge estimation for lithiumion batteries using neural networks and EKF,”
IEEE Trans. Industrial Electron.
57
(12)
4178 
4187
DOI : 10.1109/TIE.2010.2043035
Tsai M.F.
,
Peng Y.Y.
,
Tseng C.S.
,
Li N.S.
“Modeling and estimation of state of charge for LithiumIon batteries using ANFIS architecture,”
in Industrial Electronics (ISIE), 2012 IEEE International Symposium on
2012
863 
868
Fleischer C.
,
Waag W.
,
Bai Z.
,
Sauer D. U.
2013
“Adaptive online stateofavailablepower prediction of lithiumion batteries,”
Journal of Power Electronics
13
(4)
516 
527
DOI : 10.6113/JPE.2013.13.4.516
Zhao Y.
,
Yun H.
,
Liu S.
,
Jiao H.
,
Wang C.
2014
“Stateofcharge estimation for lithiumion batteries using a multistate closedloop observer,”
Journal of Power Electronics
14
(5)
1038 
1046
DOI : 10.6113/JPE.2014.14.5.1038
Kim T.
,
Qiao W.
,
Qu L.
“Realtime state of charge and electrical impedance estimation for lithiumion batteries based on a hybrid battery model,”
in Applied Power Electronics Conference and Exposition (APEC), 2013 TwentyEighth Annual IEEE
2013
563 
568
Vasebi A.
,
Partovibakhsh M.
,
Bathaee S. M. T.
2007
“A novel combined battery model for stateofcharge estimation in leadacid batteries based on extended Kalman filter for hybrid electric vehicle applications,”
Journal of Power Sources
174
(1)
30 
40
DOI : 10.1016/j.jpowsour.2007.04.011
Rezvanizaniani S. M.
,
Liu Z.
,
Chen Y.
,
Lee J.
2014
“Review and recent advances in battery health monitoring and prognostics technologies for electric vehicle (EV) safety and mobility,”
Journal of Power Sources
256
110 
124
DOI : 10.1016/j.jpowsour.2014.01.085
Plett G. L.
2004
“Extended Kalman filtering for battery management systems of LiPBbased HEV battery packs: Part 1. Background,”
Journal of Power Sources
134
(2)
252 
261
DOI : 10.1016/j.jpowsour.2004.02.031
Welch G.
,
Bishop G.
“An introduction to the Kalman filter,”
University of North Carolina at Chapel Hill Chapel Hill
NC, USA
Technical Report1995
Plett G. L.
2004
“Extended Kalman filtering for battery management systems of LiPBbased HEV battery packs: Part 3. State and parameter estimation,”
Journal of Power Sources
134
(2)
277 
292
DOI : 10.1016/j.jpowsour.2004.02.033
Waag W.
,
Käbitz S.
,
Sauer D. U.
2013
“Experimental investigation of the lithiumion battery impedance characteristic at various conditions and aging states and its influence on the application,”
Applied Energy
102
885 
897
DOI : 10.1016/j.apenergy.2012.09.030
Tang X.
,
Zhang X.
,
Koch B.
,
Frisch D.
“Modeling and estimation of Nickel Metal Hydride battery hysteresis for SOC estimation,”
in Prognostics and Health Management, 2008. PHM 2008. International Conference on
2008
1 
12
Eichi H. R.
,
Chow M.Y.
“Modeling and analysis of battery hysteresis effects,”
in Energy Conversion Congress and Exposition (ECCE), 2012 IEEE
2012
4479 
4486
Wiegman H. L. N.
,
Ph.D. dissertation
1999
“Battery state estimation and control for power buffering applications,”
Dept. Elect. and Comp. Eng., Uni. of Wisconsin  Madison
Ph.D. dissertation
Roscher M. A.
,
Sauer D. U.
2011
“Dynamic electric behavior and opencircuitvoltage modeling of LiFePO4based lithium ion secondary batteries,”
Journal of Power Sources
196
(1)
331 
336
DOI : 10.1016/j.jpowsour.2010.06.098
He H.
,
Xiong R.
,
Zhang X.
,
Sun F.
,
Fan J.
2011
“Stateofcharge estimation of the lithiumion battery using an adaptive extended Kalman filter based on an improved thevenin model,”
IEEE Trans. Veh. Technol.
60
(4)
1461 
1469
DOI : 10.1109/TVT.2011.2132812
Lee S.
,
Kim J.
,
Lee J.
,
Cho B. H.
2008
“Stateofcharge and capacity estimation of lithiumion battery using a new opencircuit voltage versus stateofcharge,”
Journal of Power Sources
185
(2)
1367 
1373
DOI : 10.1016/j.jpowsour.2008.08.103
Dai H.
,
Wei X.
,
Sun Z.
,
Wang J.
,
Gu W.
2012
“Online cell SOC estimation of Liion battery packs using a dual timescale Kalman filtering for EV applications,”
Applied Energy
95
227 
237
DOI : 10.1016/j.apenergy.2012.02.044
Kahaner D.
,
Moler C.
,
Nash S.
1989
Numerical methods and software
PrenticeHall, Inc.
Upper Saddle River, NJ, USA
Chen M.
,
RinconMora G. A.
2006
“Accurate electrical battery model capable of predicting runtime and IV performance,”
IEEE Trans. Energy Convers.
21
(2)
504 
511
DOI : 10.1109/TEC.2006.874229
Liao L.
,
Zuo P.
,
Ma Y.
,
Chen X.
,
An Y.
,
Gao Y.
,
Yin G.
2012
“Effects of temperature on charge/discharge behaviors of LiFePO4 cathode for Liion batteries,”
Electrochimica Acta
60
269 
273
DOI : 10.1016/j.electacta.2011.11.041
Gomez J.
,
Nelson R.
,
Kalu E. E.
,
Weatherspoon M. H.
,
Zheng J. P.
2011
“Equivalent circuit model parameters of a highpower Liion battery: Thermal and state of charge effects,”
Journal of Power Sources
196
(10)
4826 
4831
DOI : 10.1016/j.jpowsour.2010.12.107
Shuo T.
,
Munan H.
,
Minggao O.
2009
“An experimental study and nonlinear modeling of discharge behavior of valveregulated lead acid batteries,”
IEEE Trans. Energy Convers.
24
(2)
452 
458
DOI : 10.1109/TEC.2008.2001447
Kim J.
,
Lee S.
,
Cho B.
2012
“Discharging/charging voltagetemperature pattern recognition for improved SOC/capacity estimation and SOH prediction at various temperatures,”
Journal of Power Electronics
12
(1)
1 
9
DOI : 10.6113/JPE.2012.12.1.1