Thermoelectric Modeling and Online SOC Estimation of Li-Ion Battery for Plug-In Hybrid Electric Vehicles

The increasing oil price, energy demand, and environmental concern are leading to a global switch towards Plug-In Hybrid Electric Vehicles (PHEVs). In a PHEV, Li-ion battery is considered as the primary propelling source. Therefore, an accurate battery model is required to predict the I-V characteristic and dynamic behavior of a battery.This paper presents a highly effective thermoelectric model of Li-ion battery developed in Simulink. An algorithm is proposed for estimation of state of charge (SOC) and open circuit voltage (OCV) adaptively to notify the exact SOC level for better utilization of battery power and optimal vehicle performance. Thermal behavior of Li-ion battery is investigated for wide temperature range and its effect on resistance, capacity, and OCV is recorded.Theminimum SOC level to which battery can get depleted is calculated using gradient method.The proposed simulation results are analyzed with those of earlier models and found to be better.


Introduction
Internal Combustion Engine (ICE) based automobiles have been causing toxic emissions, global warming, and environmental and ecological danger.Petroleum as a finite fuel and increasing prices of crude oil were motivation to find alternative approaches to propel vehicles.The world is now moving towards hybrid vehicles, which contain an alternative power source along with ICE to reduce the liquid fuel consumption without affecting vehicle performance.
Hybrid Electric Vehicles (HEVs) and PHEVs are now available in the market with reduced petroleum consumption.In PHEV, the battery is used as primary power source and ICE as secondary power source.It provides a longer driving range and easy refueling (recharging) with reduced liquid fuel consumption and toxic emissions.Thus, PHEV is a means to reduce the energy demand and replacing the liquid fuel consumption by storing electrical energy in large on-board rechargeable batteries with high fuel economy and better energy efficiency [1].
Depending upon the application and specification different types of batteries are available in market like leadacid, Nickel, and Li-ion based batteries.Among these, Li-ion batteries are preferred in portable electronics today because they provide higher specific energy, higher specific density, extra durability, and lower self-discharge rate with better safety issues [2,3].
Since batteries play a vital role in PHEVs, it is essential to study their behavior.A holistic understanding of the same would result in better performance of the vehicle.The battery can be modeled using electrochemical, mathematical, analytical and stochastic, impedance based, and electrical circuit based models.The electrical circuit based model depicts the battery behavior easily [4].Electric models are of varying degrees of complexity to capture battery performance with respect to a set of parameters which are explained in detail in the next section.In this paper, a second-order electrical battery model considering diffusion and double layer effects and self-discharge current is proposed.It reflects the effect of temperature on various parameters.And a modified SOC estimation is considered here which incorporates both open circuit voltage method and ampere hour counting.To inquire about instantaneous SOC and OCV during battery usage, an adaptive estimation algorithm is developed.To use the battery, effectively without deteriorating battery health, the threshold SOC value of the battery is also computed.

Modelling and Simulation in Engineering
This paper is organized with different sections; unfolding existing battery models in literature is in Section 2. Section 3 contains the proposed SOC estimation method with the detailed description of weighting factor and correction factor.Section 4 describes the development of the proposed 2RC thermoelectric model.Section 5 discusses the simulation results.Section 6 determines the threshold SOC and Section 7 finally concludes the paper.

Overview of Existing Models
Sean proposes a PSPICE macromodel showing the voltage dependency on SOC, discharge current, resistance, and capacity variation with respect to temperature [5].This is further used by [6] to propose a discrete-time model which is capable of battery lifetime estimation.In 1994, The National Renewable Energy Laboratory (NREL) modeled Liion battery with a voltage source and internal resistance as a function of SOC, temperature, and current flow direction in ADVISOR.Saft America developed the high-power Liion cells and implemented 2-capacitance battery model in PSPICE.It shows a slightly better performance in comparison to NREL's model [7].Chen et al. proposed a model to be used with an equilibrium potential and two internal resistances  1 and  2 where  1 is a function of discharge current, temperature, and life cycle and  2 is a function of state of discharge and temperature but did not count for transient response of the battery [8].
Gao et al. demonstrated a dynamic model of Li-ion battery which depicts a capacity variation on the basis of C-rate and temperature change, equilibrium potential, and transient response of the battery [9].It consists of equilibrium potential, internal resistance (with two components), and a capacitor (transient response of charge double layer) but may result in a better performance by involving self-discharge current and diffusion effect between electrodes.Tremblay et al. presented a battery model, for dynamic simulation software, and added the same in the SimPowerSystems library MATLAB/Simulink.It consists of an internal resistance and a voltage source which is a nonlinear function of battery SOC.It does not account for Peukert, self-discharge, memory effect, and temperature variation [10].
Lee et al. used Li-ion battery model with internal resistance, one RC combination, and a voltage source (OCV as a function of SOC).They estimated SOC using ampere-hour counting and capacity estimation, neglecting coulomb based counting which is also required for accurate SOC estimation [11].Using lumped model, the SOC estimation algorithm is developed at varying temperatures [12,13].It consists of one resistance with two components (series and charge transfer), one RC ladder (diffusion resistance and diffusion capacitance), and voltage source (OCV) but does not account for self-discharge current and transient behavior of a battery.Bhide and Shim developed a circuit based Li-ion battery model using AMESim and also represented the temperature rise in core and the crust.But the model considers only particular discharge rate and different temperature and discharge rate factor functions for different rates [14].Reference [15] combined the electric model developed in [10] and thermal model developed in [14] to derive a thermoelectric analytical model.This model can inspect the behavioral change of battery due to temperature variation, but contains others lacking in [10].
For online SOC estimation, [16] proposes the model with internal resistance, OCV, and two RC circuit combinations.For SOC estimation they considered only coulomb counting and neglected voltage based counting which is required for better accuracy.The model developed in [17] is a blend of previous models and overcomes few of their limitations.It predicts runtime, steady state, and transient response accurately by capturing all the dynamic electrical characteristics of batteries.The RC network is modeled to account the effect of self-discharge losses due to long time storage and also includes transient response but it does not include thermal effects.Reference [18] used a model developed in [17] and incorporated temperature and capacity fading effect, to propose a dynamic model of Li-ion battery using MATLAB/Simulink.To determine online SOC of Li-ion battery, [19] estimated electrical parameters with temperature variation.Determination of battery SOC using a secondorder model is introduced in [16,20].
Kroeze and Krein proposed two models for predicting SOC, terminal voltage, and power losses.These are the following: (1) SOC can be predicted when temperature and cycle number are given and (2) transient behavior of terminal voltage can be figured out where each parameter is a function of SOC [21].Zhang and Chow constructed an equivalent circuit of battery cell which is based on Thevenin's theorem.It describes the SOC variation with current and considers the battery relaxation effect too but lacks online-parameter variation which is important in HEV/PHEV applications [22].This also neglects self-discharge current which results in an approximate error of 3% between experimental and estimated SOC. Randles' model [23] developed for lead-acid batteries is remapped by Gould et al. [24].They implemented an equivalent circuit model to determine state of function of Li-ion battery but did not consider the temperature effect.Based on experimental results, 2RC battery model is proposed and mathematical modeling is performed in [25], but self-discharge current is not discussed.References [26][27][28] also proposed 2RC battery model, but no discussion of self-discharge current is performed.Reference [29] has also used 2RC battery model and temperature effect is also incorporated, but again self-discharge current is not discussed.
A large number of researchers attempted to calculate exact SOC of the battery.Pang et al. [30] used OCV method to calculate the SOC (voltage based SOC, i.e., SOC V ) of the battery.The OCV based SOC estimation technique is advantageous in various aspects as follows: (i) OCV versus SOC characteristic is independent of the age of the Li-ion battery [31] and (ii) this is very accurate but requires some rest time [32].Ampere-hour counting method (current based SOC, i.e., SOC  ) is a suitable method to estimate SOC of the battery as it is easy, direct, and easily implementable.If the current measurement is accurate, then the method is also reliable.But it may have some initial value or accumulated error problems [33].To overcome the shortcomings of both and to utilize the added advantages, these two methods can be combined together.References [19,34] identified the contribution of both SOC V and SOC  together to estimate accurate SOC of the battery but do not include the effects of temperature.

Proposed Model
The proposed model considers the effects of temperature as an independent variable and incorporates self-discharge current as well in calculations.The objective is to model a Li-ion battery to represent its actual characteristics to achieve high accuracy and robustness in run-time SOC estimation.The model is aimed to simulate the dynamic behavior of a Li-ion battery as a second-order equivalent circuit in SIMULINK.All the parameters in the proposed model are multivariable functions of the SOC, current, and temperature.

SOC Estimation.
It is appreciable to characterize the Liion battery to dynamically compute the SOC even in case of temperature variation.Li-ion battery has a very eminent effect of temperature on its performance and various parameters.Under optimal temperature range batteries behave as prescribed, but outside battery cell experiences severe loss of capacity.To characterize the battery performance under the influence of temperature, thermal effect during modeling is deemed.Temperature dependent modeling provides the pertinent information about the parameters under temperature variation.
The vehicle performance is characterized by SOC of battery defined as the ratio of remaining capacity to fully charged capacity: SOC 0 is the initial SOC level of battery.
3.1.1.SOC V Calculation.Cell voltage under reversible conditions, that is, all the reactions are balanced, is called equilibrium voltage which is occasionally referred to as OCV or rest voltage.With this OCV, voltages based SOC (SOC V ) can be estimated using where  0 is battery terminal voltage when SOC = 0% and  1 is battery terminal voltage when SOC = 100%.But, due to change in temperature, equilibrium voltage of battery at any temperature  gets changed as / is temperature coefficient and is constant for the considered temperature range.So the consideration of this OCV with temperature effect will lead to modifying SOC V and will contribute in the final SOC calculation.

SOC 𝑖 Calculation.
The coulomb counting method involves the current integration flowing through the battery to get SOC  : is battery capacity in Ah.Due to change in temperature, the cell reaction rate gets changed which has been depicted here using ( 5), (6), and (7).From the Arrhenius equation, the reaction rate is given as 0 is reaction constant,  is gas constant,   is activation energy, and  is operating temperature.During the electron transfer reaction, electrons require the additional amount of energy to surmount the energy barrier called the activation energy (  = J⋅mol −1 ) which depends on temperature.As for every 10 ∘ C temperature increase, current gets doubled so for Δ temperature change, reaction rate ratio is articulated as is the reaction rate (mole/s) and can be expressed as current.Suppose  1 is the reaction rate at a temperature ( + Δ), that is,  1 = ( + Δ), and  2 at temperature , that is,  2 = ().To represent the effect of temperature on   , ( 5) and ( 6) are equated as (7) and hence   is presented as (8): Li-ion batteries exhibit self-discharge phenomenon even at moderate oxidation levels.It is primarily due to loses occurring at the negative electrode, which results from several side reactions, each with their own activation energy and rate constant.From Arrhenius equation self-discharge current can also be modeled as (5).The battery capacity also affects it; hence self-discharge current is sculpted as (9) considering the effect of temperature on activation energy from (8): 3.1.3.Modified SOC Calculation.Temperature dependent SOC can be deduced as (10) by combining SOC V and SOC  with a weighting factor . Charging and discharging efficiency influence battery dynamics to a great extent; weighting factor allied with SOC  should govern the combined SOC.
Correction factor (CF)  is integrated here, which is a function of SOC as (11) which helps in getting exact SOC during discharging:

Weighting Factor Calculation.
To calculate the weighting factor, the value of OCV and time required to get steady OCV should be considered because under steady condition SOC V has higher accuracy.The entire OCV range is divided into 100 sections; each section weighting value is calculated as According to experimental data, average time taken by OCV to get steady condition is   = 1278 s and time between two samples is  = 192.78s.As  <   , weighting factor  is deduced as (13).Table 1 lists the values of weighting factors and correction factor at various SOC levels: Table 2 brings together the expected SOC values to be generated during discharging and values obtained from simulation.The modified SOC estimation method calculates (very nearby) expected values from the range of 100% to 50%.The error is very less in the range of 0-0.5%.From 50-36% SOC ranges, the error between the expected and obtained SOC is lesser than 2%, but error starts increasing 35% onward and continues to empty state.

Correction Factor Determination.
Determining CF is an essential and more intricate process.Since the battery provides electrical energy to the load (vehicle), more precise SOC estimation is a very important consideration under dynamic state of charge/discharge for the driver and vehicle controller.Over the time sphere, under various operational conditions, the current transfer rate affects the CF so it is closely related to the charge/discharge rate.Battery performance is expected to get varied with the load; and CF gets modified for the varying conditions.CF at the fixed temperature is a function of SOC.Reference [35] substantiates variation of CF with SOC and temperature and justifies it with the mathematical expression derived.When battery alone is considered, then CF can be estimated as (11) but as it gets associated with the load it varies significantly and can alter the performance of the vehicle.
At high SOC levels, the kinetic rates are high in the presence of more unreacted active mass and existence of higher level of electrolyte concentration which leads to rapid material conversion.At low SOC levels, the kinetic rates are relatively low and lead to 100% charge efficiency [36].Change in temperature directly influences the rate of reaction and hence governs CF.So, CF is a function of temperature  and SOC and demonstrated by the simulation results.CF derived can be stored or calculated "on-the-fly" and applied incrementally to correct the SOC.The CF variation with the SOC and temperature is shown in Figure 1.Exploration of CF for the range of SOCs with wide temperature series is performed, that is, 0 ∘ C, 25 ∘ C, 50 ∘ C, and 60 ∘ C (very low, room, and high temperatures), which helps to find a healthy operating range of the battery.The pattern achieved is nearby Cauchy's distribution.At room temperature (25 ∘ C), CF is higher at lower SOCs and gets low down with the higher SOC values.The same pattern is observed for the 0 ∘ C, 50 ∘ C, and 60 ∘ C; only the quantitative measures are changed.Estimate the CF at room temperature and store the values in look-up table.According to the developed relation as given in (14), CFs are calculated "onthe-fly" at various SOCs and temperatures onwards. is the error tolerance, applicable in the cases of temperature variation: 3.4.Proposed SOC Estimation Algorithm.Continuous recording of SOC during vehicle propulsion is important to update power management system.The power management system, then, takes decision to toggle between battery and engine for minimum fuel consumption and better energy efficiency.Change in temperature due to any reason (environmental

Development of Proposed Thermoelectric Model
In a Li-ion battery, Li-ions move from anode to cathode through an organic electrolyte during discharging.During charging the reverse process is followed.Anode and cathode have layered structure and Li-ions move in between layers.
The governing chemical reactions are as follows: The time varying response of battery can be represented by a combination of diffusion and double layer resistance and capacitances [37,38].Proposed battery model is shown in Figure 3, where  is the voltage source which is a function of charge remaining in the battery,  int is internal resistance, and  1 ,  2 ,  1 , and  2 are diffusion and double layer resistances, and capacitances, respectively.The output voltage of the battery is lower than no load voltage due to voltage drop and is represented as  int .
Temperature effect and self-discharge current are included for analyzing run-time characteristics.A number of Li-ion cells are connected in series and in parallel to simulate a powerful battery to be used in hybrid vehicles.The output voltage is calculated as (15) where () is time dependent OCV, () is the time dependent current, and multiplier term associated with () represents the impedance offered.
From (2), temperature dependent voltage based SOC, that is, SOC V (), can be intended as (16) and from (4) temperature dependent current based SOC, that is, SOC  (), can be calculated as (17): SOC at (any temperature) operating temperature  can be calculated using (10), (16), and (17).Electric component values are simulated and determined with varied SOC as shown in Figure 4.The model is largely practicable because all the extracted RC parameters are approximately constant over 20%-100% SOC and change exponentially within the lower limits of SOC of the battery.
Deduce SOC v (T) using (15) New SOC(t) at (t + 1) Assuming initial level of Discharging the battery with user defined load Is SOC = min threshold?
New SOC is considered as charge remaining in battery.
Display to driver

Stop Yes No
Tracking V(t) and I(t) during discharge and I (s-d) as in ( 5) and ( 9) Effect of temperature on I(t)

Simulation Results
Concentration of Li-ions of solid phase with respect to intercalating material (Li-ion) represents OCV and gets affected by changes in chemical composition, pressure, and temperature.The proposed model is simulated in wide temperature range of 25-70 ∘ C and variation of OCV is intrigued with respect to SOC as shown in Figure 5. Increment in OCV is recorded as the temperature increases from ambient and vice versa for the decrease in temperature but fluctuates roughly with temperature variation.At lower temperature, battery cell is inefficient due to more impedance, but, at higher temperature, efficiency improves because the rate of chemical reactions increases due to lesser impedance offered.With the increase in temperature, resistance decreases which is shown in Figure 7.
Johnson locates an increment of three times in resistance as the temperature drops from ambient to 0 ∘ C [40].Johnson et al. articulate that resistance can get varied by eight times depending upon discharge rate [41].For proposed model, it increases by 5.5 times than the ambient level as the temperature drops to the 0 ∘ C. Battery capacity extensively varies with operating conditions and strongly depends on its internal impedance.With changing temperature, impedance of a battery varies a lot, hence capacity.At lower temperatures, cell reaction rate decreases, hence low current flows and vice versa at higher temperatures as shown in Figure 7. Therefore, usable capacity decreases at low temperatures and increases at high temperatures.At 0 ∘ C, capacity is found to be 17.11% of its value at ambient.At higher temperatures, current drawn from battery increases rapidly, and battery depletes swiftly.At 40 ∘ C the usable capacity is 5.6% more than the value at ambient.These results are better than other reported works.
For the fixed temperature, battery resistance variations are small for a complete SOC range from 100% to 10% as shown in Figure 8(a) and results are comparable with [13].To accommodate all temperature data, a larger scale range is chosen; that is why the plots seem to be constant, but variation is observed when plotted for individual cases as shown in Figure 8  change varying from 0.005% to 0.26% on 273 K temperature is observed.Almost the same pattern with little bit variation is observed for other temperatures as well.Battery operation at enormously high temperatures causes an increase in rate of thermal reactions leading to even higher heat generation.Many electrochemical side reactions are damaged in the process.An efficient convective cooling system can avoid such high temperature excursions during repeated cycling.
Simulation of HEV is performed over the ECE EUDC driving cycle and is presented separately in [30].The authors also observed that the proposed SOC estimation method gives better fuel efficiency with 2RC battery model than the default SOC estimation method as shown in Table 3.

Threshold SOC Determination
Li-ion batteries are used extensively in PHEVs to propel the vehicle as primary power source; hence it is essential to monitor continuously its charge status in terms of SOC.Overcharging and depletion below a specified level deteriorates battery health.Since batteries are expensive, it is compulsory to take care of their health for durability.Battery mainly operates in two modes: (1) charge depletion and (2) charge sustaining.In Hybrid Electric Vehicles (HEVs) both modes together are used in such a way that initial and final SOC at the end of the trip are the same.In HEVs, batteries get charged either by the engine (indirectly by fuel) or by regenerative braking.
As far as PHEVs are concerned, the batteries are charged through main supply at home or parking lots; therefore it can be depleted completely at the end of the trip.However, depleting a battery below a threshold level deteriorates its health and hence performance.In available literatures, authors do not find proper justification to choose the lowest SOC level up to which battery should be allowed to discharge.Based on data collected during discharge from 100% to 0% SOC at various temperatures, the authors propose determination of the threshold SOC level using gradient method.The curves shown in Figure 5 where  1 = −6.9991+ 007,   2. Proposed threshold level of SOC is 35% and operating under this is not recommended as exposed in Figure 9. Operating a battery at higher temperatures, threshold occurs sooner and at lower temperatures it occurs later as compared to the room temperature case, which can also be surveyed in Figure 5.

Conclusion
Today, the major challenge is to utilize battery power efficiently without affecting its health to minimize the fuel consumption in PHEVs.The authors proposed a 2RC battery model, for HEV/PHEV applications.They have included temperature and self-discharge current both while modeling in Simulink and developing mathematical expressions.Along with this, the authors used both open circuit voltage and coulomb counting method together to estimate battery SOC.This SOC estimation was further improved using a correction factor.Since battery performance greatly depends on temperature, a generalized thermoelectric model is proposed to characterize its dynamics.
Effect of temperature on various parameters like current, OCV, capacity, and resistance is determined.It is observed that, from 0 ∘ C to 40 ∘ C, OCV values do not vary much, but beyond this they get amplified because of transient effects.Higher resistance at and below room temperature and a lesser amount of offered resistance beyond room temperature are demonstrated.Hence, capacity becomes very high at higher temperature, but at the same time current drawn from the battery turns out to be high; consequently, battery depletes rapidly.The online estimation of SOC is carried out to use the battery power judiciously; in addition, a threshold SOC, that    is, 35%, is gauged and suggests not to use the battery below this threshold due to health reasons.

Figure 1 :
Figure 1: Correction factor variation with temperature and SOC.
From 20 ∘ C to 40 ∘ C curves are very close to each other, but beyond 40 ∘ C OCV increases, which depicts the true property of the battery.While comparing with [39], at 0 ∘ C and 25 ∘ C not much variation in OCV is found, but a significant change at 60 ∘ C is observed as shown in Figures 6(a), 6(b), and 6(c).

Figure 7 :
Figure 7: Resistance and capacity variation with temperature.

Figure 8 :
Figure 8: Resistance variation with varying SOC; (a) all temperatures in same plot and (b) plots for individual temperature.

Table 2 :
Comparison between expected and obtained values of SOCs.

Table 3 :
Fuel economy comparison chart for SOC methods.