Battery Energy Storage
Battery Energy Storage
99
100 Electric and Hybrid Vehicles: Design Fundamentals
Table 4.1
Nominal Energy Densities of Sources
Nominal Specific Energy
Energy Source (W h/kg)
Gasoline 12,500
Diesel 12,000
Biodiesel 10,900
Natural gas 9,350
Ethanol 8,300
Methanol 6,050
Hydrogen 33,000
Coal (bituminous) 8,200
Lead-acid battery 35
Li-polymer battery 200
Flywheel (carbon-fiber) 200
• Nickel-metal-hydride (NiMH)
• Lithium-ion (Li-ion)
• Lithium-polymer (Li-poly)
• Sodium-sulfur
The Li-ion battery technology is the most promising among the four battery
chemistries mentioned above. It must also be stated that there are several
different types of Li-ion battery technologies that are being developed for
electric and hybrid vehicles.
The development of batteries is directed toward overcoming the signifi-
cant practical and manufacturing difficulties. The theoretical predictions are
difficult to match in manufactured products due to practical limitations. The
theoretical and practical specific energies of several batteries are given in
Table 4.2 for comparison.
The battery technology has gone through extensive research and develop-
ment efforts over the past 30 years, yet there is currently no battery that can
deliver an acceptable combination of power, energy, and life cycle at a reason-
able cost for high-volume production urban usage electric vehicles. However,
the extensive research and interest in alternative vehicles have resulted in
several promising battery technologies. Pure electric vehicles are available
commercially, but the cost is prohibitive since these require large capacity
batteries. Hybrid electric vehicles minimize the battery capacity through
using the combination of IC engine and electric machines. Even though
Battery Energy Storage 103
TABLE 4.2
Specific Energies of Batteries
Specific Energy (W h/kg)
Battery Theoretical Practical
Lead-acid 108 50
Nickel-cadmium 20–30
Nickel-zinc 90
Nickel-iron 60
Zinc-chlorine 90
Zinc-bromide 70
Silver-zinc 500 100
Sodium-sulfur 770 150–300
Aluminum-air 300
Nickel metal 70
Hydride 150
Li-ion
Finally in this chapter, battery-pack design for electric and hybrid vehicles’
applications will be addressed.
Vcell
+ – Negative
electrode
Separator
Positive Cell
electrode container
Negative electrode
(a) (b) Electrolyte
FIGURE 4.1
Components of a battery cell. (a) Cell circuit symbol; (b) cell cross-section.
Battery Energy Storage 105
releases electrons to the external circuit. For battery charging, a source with
voltage higher than the battery terminal voltage has to be applied so that
current can flow into the battery in the opposite direction. During charging,
electrons are released at the positive electrode and consumed at the negative
electrode; consequently, the positive electrode is oxidized and negative elec-
trode is reduced.
Regardless of the battery cell chemistry, redox reactions take place at the
electrodes during both cell charging and discharging for the release and
absorption of electrons at the terminals. The generalized redox reactions are
given by [1]
Charge
aA ← → cC + nE+ + ne − (4.1)
Discharge
Charge
bB + nE + + ne − ← → dD (4.2)
Discharge
Charge
aA + bB ← → cC + dD (4.3)
Discharge
Chemical reactions 4.1 and 4.2 illustrate that electrons are released and
absorbed during any redox reaction. The positive electrode reaction 4.1
shows that during cell charging, species A within the electrode is oxidized
and becomes energized species C, releasing electron(s) into the external cir-
cuit and positive ion(s) into the electrolyte. Similarly, the negative electrode
reaction 4.2 shows that species B at the electrode combines with positive
ion(s) from the electrolyte and electron(s) from the external circuit to form
energized species D. The converse is true at the two electrodes during cell
discharging. The coefficients a, b, c, and d represent the numbers of moles
associated with the species in the reactions; the coefficient n represents the
number of electrons and ions involved in the redox reactions.
In electric traction applications, battery cell operation is that of cell dis-
charging when the energy is supplied from the battery to the electric motor
for propulsion power and of cell charging when energy is supplied from
an external source to store energy in the battery. In conventional vehicles,
battery cells supply power to electrical accessories while discharging, and
accept energy from an external device to replenish the stored energy during
charging. We will next review the redox reactions during cell charging and
discharging in a few battery chemistries, starting with the lead-acid battery
Battery Energy Storage 107
RL
Electron flow
+ –
e– e– e– e–
H+
PbO2(s) H+
H+ SO42– Pb(s)
H+
SO42–
Water
PbSO4(s) PbSO4(s)
2H2O
FIGURE 4.2
Lead-acid battery: cell discharge operation.
cell. Lead-acid is still the battery choice for powering electrical accessories in
conventional, electric, and hybrid electric vehicles.
Figure 4.2 shows the cell discharge operation of a lead-acid battery cell into
a passive resistive element. The positive electrode made of lead-oxide (PbO2)
is reduced by consuming electrons and ions. The electron supply is through
the external circuit which originates at the negative electrode. The current
flow is therefore out of the positive electrode into the electrical load with the
battery acting as the source. The positive electrode reaction is given by
Electric current
Electron flow
+
e– e– e– e–
H2O
PbO2(s) PbSO4(s) H2O
PbSO4(s) Pb(s)
SO42–
H+
H+
H+
SO42–
FIGURE 4.3
Lead-acid battery: cell charge operation.
The cell charge operation is the reverse of the cell discharge operation. An
external electrical source supplies current into the battery to reverse the
chemical reactions as shown in Figure 4.3. During cell charging, the lead
sulfate is converted back to the reactant states of lead and lead oxide. The
electrons are consumed from the external source at the negative electrode,
while the positive electrode produces the electrons. The current flows into
the positive electrode from the external source, thereby delivering electrical
energy into the cell where it gets converted into chemical energy. The posi-
tive electrode is oxidized, releasing electrons during cell charging as follows:
Discharge
NiOOH + H 2O + e − ← → Ni ( OH ) + OH −
Charge 2
Discharge
Cd + 2OH − ← → Cd(OH)2 + 2e −
Charge
Discharge
Cd + 2NiOOH + 2H 2O ← → 2Ni(OH)2 + Cd(OH)2
Charge
In NiMH batteries, the positive electrode is a nickel oxide similar to that used
in a NiCd battery, while the negative electrode has been replaced by a metal
hydride which stores hydrogen atoms. The concept of NiMH batteries is
based on the fact that fine particles of certain metallic alloys, when exposed to
hydrogen at certain pressures and temperatures, absorb large quantities of the
gas to form the metal hydride compounds. Furthermore, the metal hydrides
are able to absorb and release hydrogen many times without deterioration.
The two electrode chemical reactions in a NiMH battery are as follows:
At the positive electrode,
Discharge
NiOOH + H 2O + e − ← → Ni(OH)2 + OH −
Charge
Discharge
MH x + OH − ← → MH x −1 + H 2O + e −
Charge
where
M stands for metallic alloy, which takes up hydrogen at ambient tempera-
ture to form the metal hydride MHx
x is the number of hydrogen atoms absorbed
110 Electric and Hybrid Vehicles: Design Fundamentals
Discharge
NiOOH + MH x ← → Ni(OH)2 + MH x −1
Charge
4.3.1 Battery Capacity
The amount of charge released by the energized material at an electrode
associated with complete discharge of a battery is called the battery capacity.
The capacity is measured in A h (1 A h = 3600 C or Coulomb, where 1 C is the
charge transferred in 1 s by 1 A current in the SI unit of charge).
The theoretical capacity of a battery can be obtained by Faraday’s law of
electrolysis, which states that the mass of the elemental material altered at
an electrode is directly proportional to the element’s equivalent weight for
a given quantity of electrical charge. The equivalent weight of the elemen-
tal material is given by the molar mass divided by the number of electrons
transferred per ion for the reaction undergone by the material. This number
is known as the valency number of ions for the substance. Mathematically,
Faraday’s law can be written as
Q Mm
mR = (4.4)
F n
where
mR is the mass of the limiting reactant material altered at an electrode
Q is the total amount of electric charge passing through the material
F is the Faraday number or Faraday constant
Mm is the molar mass
n is the number of electrons per ion produced at an electrode
Mm/n is the equivalent weight of the reactant substance. The Faraday num-
ber is given by the amount of electric charge carried by one mole of electrons.
The number of molecules or atoms in a mole is given by the Avogadro num-
ber NA which is equal to 6.022045 × 1023 mol−1. The amount of charge in one
Battery Energy Storage 111
QT = xnF C (4.5)
mR
x=
Mm
Here
mR is the mass of the reactant material in kg
Mm is the molar mass of that material in g/mol
mR n
QT = 0.278 F Ah (4.6)
Mm
The cells in a battery are connected in series and the capacity of the battery
is dictated by the smallest cell capacity. Therefore, QTbattery = QTcell. Six battery
cells connected in series to form a battery are shown in Figure 4.4.
– 12.9 V +
2.15 V
– +
Q Q Q Q Q Q
FIGURE 4.4
Battery cells connected in series.
112 Electric and Hybrid Vehicles: Design Fundamentals
Battery Load
+ Ev
Ri I = constant
SoD (to) = 0
± Ev Vt RL SoD (td) = QT
_
SoD
(a) (b) QT
FIGURE 4.5
(a) Steady-state battery equivalent circuit. (b) Battery open circuit voltage characteristics.
open circuit voltage characteristics are shown in Figure 4.5b. As the battery
is gradually discharged, the internal voltage decreases, while the internal
resistance increases. The open circuit voltage characteristics have a fairly
extended plateau of linear characteristics with a slope close to zero. The open
circuit voltage is not a good indicator of the state of charge; state of charge
of a battery pack needs to be calculated considering discharge current char-
acteristics, battery chemistry, temperature effects, and number of charge/
discharge cycles. Once the battery is completely discharged, the open circuit
voltage decreases sharply with more discharge.
Vt
VFC
Vcut
0 QP SoD
FIGURE 4.6
Battery terminal voltage.
Battery Energy Storage 113
i
A
t=0
Vt RL
+ –
Battery
FIGURE 4.7
Battery capacity measurement.
CP =
∫ i(t)dt
to
(4.7)
where
to is the time at which battery is fully charged
tcut is the time at which battery terminal voltage is at Vcut
Vt
I2
I1
FIGURE 4.8
Constant current discharge curves.
C 100 A h
rate is = 20 A
5 5h
and
100 A h
2C rate is = 200 A
0.5 h
dq
i(t) =
dt
where q is the charge moving through the circuit. The instantaneous theo-
retical state of charge SoCT(t) is the amount of equivalent charge remaining
at the positive electrode and ready to be released by the energized material.
If the state of charge is QT at the initial time to, then SoCT(to) = QT. For a time
interval dt
dSoCT = − dq
= −i(i)dt
Battery Energy Storage 115
Integrating from the initial time to to the final time t, the expression for
instantaneous state of charge is obtained as
t
∫
SoCT (t ) = QT − i ( τ ) dτ
to
(4.8)
SoCT (t) =
QT −
∫ i(τ)dτ × 100%
to
QT
The state of charge will be increasing when a battery is being charged. If the
state of charge is zero initially at t = 0, the state of charge at time t expressed
in percentage form is given by
SoCT (t) =
∫ i(τ)dτ × 100%
0
QT
QT − SoCT (t)
DoD(t) = × 100%
QT
t
=
∫ i(τ)dτ × 100%
to
(4.10)
QT
116 Electric and Hybrid Vehicles: Design Fundamentals
4.3.9 Battery Energy
Energy of a battery is measured in terms of the capacity and the discharge
voltage. To calculate the energy, the capacity of the battery must be expressed
in coulombs. 1 A h is equivalent to 3600 C, while 1 V refers to 1 J (J for joule)
of work required to move 1 C charge from the negative to positive electrode.
Therefore, the stored electrical potential energy in a 12 V, 100 A h battery is
(12)(3.6 × 105) J = 4.32 MJ. In general, the theoretical stored energy
ET = VbatQT
where
Vbat is the nominal no load terminal voltage
QT is the theoretical capacity in C units
1000 Fn nmR
ET = mR Vbat = 9.65 × 107 Vbatt J (4.11)
Mm Mm
Ep =
∫ v(t) i(t) dt W h
to
(4.12)
where
to is the time at which battery is fully charged
tcut is the time in hours at which battery terminal voltage is at Vcut
v is the battery terminal voltage
i is the battery discharge current
Discharge energy E
SE = =
Total battery mass MB
Battery Energy Storage 117
The unit for specific energy is W h/kg. The theoretical specific energy of a
battery using Equation 4.9 is
nVbat mR
SET = 9.65 × 107 × W h/kg (4.13)
Mm mB
If the mass of the battery MB is proportional to the mass of the limiting reac-
tant of the battery mR, then SET is independent of mass. The specific energy of
lead-acid battery is 35–50 W h/kg at C/3 rate. Since practical energy EP varies
with discharge rate, the practical specific energy SEP is also variable.
The term energy density is also used in the literature to quantify the qual-
ity of a battery or other energy sources. The term energy density refers to the
energy per unit volume of a battery. The unit for energy density is W h/L.
4.3.11 Battery Power
The instantaneous battery terminal power is
where
Vbat is the battery terminal voltage
i is the battery discharge current
Using Kirchhoff’s voltage law for the battery equivalent circuit of Figure 4.5a,
Vt = Ev − Rii (4.15)
Power versus current characteristic is shown in Figure 4.9. Using the maxi-
mum power transfer theorem in electric circuits, the battery can deliver max-
imum power to a DC load when the load impedance matches the battery
internal impedance. The maximum power is
Ev2
Pmax = (4.17)
4Ri
Since Ev and Ri vary with the state of charge, Pmax also varies accordingly.
118 Electric and Hybrid Vehicles: Design Fundamentals
Power
Pmax
iPmax Current
FIGURE 4.9
Battery power characteristics.
Maximum power output is needed from the battery in fast discharge con-
ditions in vehicle applications, which occur when the electric motor is heav-
ily loaded. Acceleration on a slope is such a condition, when the motor draws
a lot of current to deliver maximum power required for traction.
The performance of batteries to meet acceleration and hill climbing require-
ments can be evaluated with the help of rated power specifications, which
are based on the ability of the battery to dissipate heat. The rated continuous
power is the maximum power that the battery can deliver over prolonged
discharge intervals without damage to the battery. These do not necessarily
correspond to Pmax on p–i curve of battery characteristics. The rated instanta-
neous power is the maximum power that the battery can deliver over a very
short discharge interval without damage to the battery.
P
SP = (units : W/kg ) (4.18)
MB
where
P is the power delivered by battery
MB is the mass of battery
1000
Gas turbine
Sodium sulfur
Nickel-cadmium
100
Lead-acid
Fuel cell
10
1 10 100 1000
Specific energy (W h/kg)
FIGURE 4.10
Specific power vs. specific energy (Ragone plots) of several batteries, a gasoline engine, and a
fuel cell.
Ragone relationship and is described by Ragone plots. Ragone plots are usually
obtained from constant power discharge tests or constant current discharge
plots. Let us consider the experiment of Figure 4.6 again, but this time the
current i is adjusted by varying RL such that the power output at the battery
terminals is kept constant. The experiment stops when the battery terminal
voltage reaches the cut-off voltage, i.e., Vt = Vcut. We assume that the battery
is fully charged at t = 0. The experiment is performed at several power levels
and the following data are recorded: (1) power p(t) = Vti = P, (2) time to cut-off
voltage tcut, and (3) practical energy EP = Ptcut. The plot of SP vs. SE on a log-
log scale is known as the Ragone plot. The Ragone plots of several batteries
along with alternative energy sources and IC engines are given in Figure
4.10 to give an idea about the relative power and energy capacities of these
different devices.
voltaic cells since volts or potentials are established due to the spontaneous
chemical reactions. The electrolytic cells are those in which reactions are
nonspontaneous, and effected by the imposition of external voltage greater
than the open circuit voltage of the cell. Electrical energy is expended by
the external source and work is done on the electrochemical cell. When a
battery cell is being charged to restore chemical energy within, it functions
as an electrolytic cell. The battery cells will be referred to as electrochemical
cells in general, since both discharging and charging are associated with the
secondary battery cells used in electric and hybrid vehicles.
The fundamentals of a battery cell, or for that matter, of any other electro-
chemical cell, is in the physics and chemistry involved in the energy con-
version process. In this section, we will analyze the fundamentals to gain
sufficient understanding of the process so that electrochemical cell models
can be developed. The models are useful not only for analyzing the electro-
chemical cells, but also for evaluating systems that include the electrochemi-
cal cell-based components for energy storage.
The fundamentals governing the principles and operation of an electro-
chemical cell are chemical thermodynamics, electrochemical reaction rates,
electrode kinetics, and mass transport [1–7]. Each one of these mechanisms
influences the energy conversion from stored chemical energy into electrical
energy in an electrochemical cell. The mechanisms are treated in detail in
this section.
H = U + pV
where
U is the internal energy
p is the pressure
V is the volume of the system
The product of pressure and volume (pV) is related to the work of expansion.
In galvanic devices, stable substances react spontaneously to form new
substances, which mean that the reaction can only go from a higher energy
state to a lower one according to the first law of thermodynamics. The lower
energy state that the substances will assume depends on another measure
known as the entropy of a substance. Entropy is a measure of the disorder
level for the particles (molecules, atoms, etc.) that make up a substance.
Entropy is a property that is specified for every equilibrium state of a sub-
stance. Since entropy is a property, the change in entropy in going from one
state to another is the same for all processes. The SI unit for entropy is J/K.
The second law of thermodynamics states that a spontaneous reaction in
a closed environment proceeds until the maximum entropy is reached for
the substances involved. Based on the two laws of thermodynamics, part
of the enthalpy change has to be reflected in the form of entropy change
(TΔS), where T is the temperature in Kelvin and ΔS is the change of entropy.
Therefore, the maximum available work from a chemical reaction or any
process is [8,9]
∆G = ∆H − T ∆S
∆G = ∆H o − T ∆So (4.19)
122 Electric and Hybrid Vehicles: Design Fundamentals
where
∆H o = ∑ H (products) − ∑ H (reactants)
o
f
o
f
∆G = ∆G of (products) − ∆G of (reactants)
The amount of heat content liberated from the chemical reaction between
metallic zinc and silver chloride solution when mixed at standard-state con-
ditions is −233 kJ/mol of Zn reacted, i.e., ΔHo = −233 kJ/mol. For the chemical
reaction, TΔS = −43 kJ/mol. The Gibbs-free energy in this example is
The example illustrates that all of the enthalpy change cannot be converted
to work; a minimum amount of energy has to be consumed by the entropy
change which is usually reflected in the generated heat during chemical
reactions.
Gibbs-free energy is also useful in determining the maximum thermody-
namic conversion efficiency of a galvanic cell, which is given by the ratio of
Gibbs-free energy and the total enthalpy change:
∆G T ∆S
ηEC = = 1− (4.20)
∆H ∆H
Battery Energy Storage 123
For the given example of silver-zinc battery cell, ηEC = 81.5%. The inherent high
conversion efficiency compared to the thermodynamic upper limit of efficien-
cies in heat engines is an advantage of galvanic devices. However, just like the
heat engines, the practical efficiency of galvanic devices during normal opera-
tion is much lower than the theoretical efficiency. The decrease in efficiency is
directly related to the practical currents required in practical systems.
Gibbs-free energy released in the chemical process of the electrochemical
cell imposes the theoretical limit on the maximum work that can be done by
the cell. This work is the charge transferred per mole under the force of the
open-circuit voltage (OCV) of the cell. With E being the electrode potential
difference at equilibrium, the work done by the cell can be expressed as
= −nFEo (4.21)
The theoretical upper limit of a galvanic cell potential can be obtained from
the above as
∆G
Eo = − (4.22)
nF
We will use the chemical reaction in a NiCd battery cell to calculate the
equilibrium cell potential. The chemical reaction in a NiCd battery cell is
given by
∆G = ∆G of (products) − ∆G of (reactants)
= −293.05 kJ/mol
For the NiCd battery cell, two electrons are involved in the chemical reac-
tion. The theoretical cell potential is thus
−293, 050
Eo = − = 1.52 V
2 × 96, 485
4.4.3 through 4.4.7. The nominal practical cell voltage in a NiCd cell is 1.3 V.
This is true for all electrochemical cells. The thermodynamic cell potential
only gives the theoretical upper value of cell potential for a battery chemistry.
The potential of an electrode is the potential difference between the elec-
trode and the electrolyte that it is in contact with. The electrode potential is
determined with respect to a reference electrode, since an absolute potential
value cannot be obtained. Both chemical and electrical processes contribute
to the electrode potential difference [6]. The environment at the vicinity of
an electrode is changed due to chemical activities between the electrode and
the electrolyte regardless of the electric potential difference at the solid-liq-
uid phase boundary. The measure of the work done to bring a particle to its
assumed potential is the chemical potential. Again, regardless of the changes
in the chemical environment, the transfer across the electric potential is
accomplished by electric work done in its original sense. Although one cannot
separate these two components for single species experimentally, the differ-
ences in the scales of the two environments make it possible to separate them
mathematically [7,8]. The resultant potential for these two kinds of energy
change is the electrochemical potential or simply the electric potential.
Q
x= (4.23)
nF
dQ
i= (4.24)
dt
dx i
= mol/s
dt nF
The reaction rate is typically expressed in mol/s per unit area since elec-
trode reactions are a heterogeneous process occurring only at the electrode–
electrolyte interface. The heterogeneous reaction rate depends on the mass
transfer to the electrodes and various surface effects in addition to the elec-
trode kinetics. The reaction rate per unit area is given by
dx i j
= = mol/s cm 2 (4.25)
dt nFA nF
where
A is the area in cm2
j is the current density in A/cm2
kf
O + ne −
R
kr
where kf and kr are the forward and reverse rate constants, respectively. The
rate constants are the proportionality factors linking the concentration of the
species to the reaction rates. The concentration of species undergoing oxida-
tion at a distance x from the surface and at time t will be denoted as CO(x,t);
hence, the surface concentration is CO(0,t). Similarly, the surface concentra-
tion for the species undergoing reduction is CR(0,t).
The reaction rate obtained from the product of rate constant and species
concentration can be equated to the reaction rate given by Equation 4.25 to
establish the relationship between the species concentrations and faradaic
current. Therefore, for the forward and reverse currents, we have
if
k f CO (0, t) =
nFA
and
ir
k r C R (0 , t ) =
nFA
Battery Energy Storage 127
The net current flow at the electrode is the difference between the forward
and reverse currents
For the sake of simplicity in analyzing the electrode process, we will assume
a single electron transfer (i.e., n = 1) at the electrode–electrolyte interface. In
this case, the rate constants can be related to the electrical potential across the
electrode–electrolyte interface using free energy considerations [8]. For the
standard potential E0, the forward and reverse rate constants are equal; this
constant is known as the standard rate constant and is given the symbol k0.
The rate constants at other potentials are given in terms of the standard rate
constant as
0
k f = k 0 e( F/RT )( − α )( E − E )
0
k r = k 0 e( F/RT )(1−α )( E−E )
Inserting these relations into Equation 4.26 gives the complete current–volt-
age characteristics at the electrode–electrolyte interface
( 0
i = FAk 0 CO (0, t)e( F/RT )( − α )( E−E ) − CR (0, t)e( F/RT )(1−α )( E−E
0
)
) (4.27)
RT CO (0, t)
E = Eo + ln (4.28)
nF CR (0, t)
RT CO*
Eeq = Eo + ln (4.29)
F CR*
where
Eeq is the equilibrium potential
CO* and CR* are the bulk concentrations of the oxidation and reduction
reactants, respectively
Another approach for estimating the terminal potential is the Tafel solu-
tion which assumes that the current is large in one direction or the other.
The approximation means that one of the two exponential terms in the
Butler–Volmer expression of Equation 4.27 is negligible. The Tafel solution
is given by
RT RT
E(t) = E0 + ln(i0 ) − ln (i(t)) (4.30)
αnF αnF
i0 = i f = ir = FAk 0C0*e
− α ( F/RT )(1− α )( Eeq − E0 )
(4.31)
∂C( x , t) ∂ 2C( x , t)
=D
∂t ∂x 2
where
C(x,t) is the active species concentration
D is the diffusion coefficient of the electrochemical cell
x and t are the space and time variables, respectively
δ
Cd (t) = Cbulk − i(t) (4.32)
nFAD
where
δ is the diffusion layer thickness
A is the surface area
element (CPE) instead of Fick’s second law [10]. For a CPE, the phase angle
of the frequency response remains the same for all frequencies. The CPE
transfer function used to represent the overall diffusion process is given by
Cd (s) K
H (s) = = q , 0<q<1
i(s) s
The time response of the CPE is easier to solve for simpler operations such as
constant current discharge. That time response is
The diffusion process described serves the purpose of representing both the
energy storage and the impedance. The energy within a battery is stored or
spatially distributed in the electrolyte in terms of the concentration of the
active material. The movement of active material during cell reactions is con-
trolled by the inherent impedance of the electrolyte. Both mechanisms have
been represented by the diffusion process and are represented as a coupled
mechanism from the electrical perspective. However, for certain analysis, it
is desirable to separate the source and impedance. This separation is desir-
able for certain applications. For example, the fuel cell-type electrochemi-
cal device does not store any energy; the materials for chemical reactions
in a fuel cell are supplied from external fuel. The electrochemical process
is more accurately modeled by separated energy source (the fuel) and the
impedance to the source. Another important need is for the prediction of
how much energy is left in the battery, which is essentially calculating the
SoC of the battery.
The separation of energy storage and impedance enables an improved
modeling of pulsed discharge characteristics of batteries, which is essen-
tially what takes place in electric and hybrid vehicles. For the complicated
discharge currents in such applications, it is difficult to obtain analytical
solutions of the CPE, and often one has to resort to numerical solutions.
RT Celectrode , i
EC = ln (4.34)
nF Cbulk , i
where Celectrode,i and Cbulk,i are the concentrations of the inert ions at the elec-
trode and the bulk solution, respectively. Polarization, regardless of its ori-
gin, is reflected at cell terminals as a voltage reduction from the open circuit
voltage.
hybrid vehicle systems) and for the design of these systems. For example, a
simplified battery model can be used for the dynamic simulation of a hybrid
vehicle to predict the powertrain characteristics as well as the range on elec-
tric-only operating mode. Depending on the simulation objective, the mod-
els can be represented by a set of electrical circuit components or by a simple
empirical relationship of two parameters. These two types of models are the
electric circuit models and empirical models, which are presented in this sec-
tion. The electric circuit-based models are somewhat more complex than the
empirical models, but are extremely useful for vehicle system level analysis.
On the other hand, the empirical models allow a quick evaluation on the
range of a vehicle based on the capacity or energy density.
The energy storage device models presented in this section are useful not
just for electric and hybrid vehicle applications, but also for utility power sys-
tem applications. Distributed power systems require energy storage devices
with similar features as those required for electric and hybrid vehicles.
where
E0 is the initial voltage of the cell
RΩ is the series resistance
k1 is a constant parameter
( )
E(t) = E0 + RΩi(t) + k1 ln ( 1 + i ) sgn(i) + k 2 ln 1 + qd sgn(i) + k 3 ln(1 − qs )
The constants E0, RΩ, k1, k2, and k3 depend on the properties of the elec-
trochemical cell and can be determined from experimental data. While the
Hartley model gives a mathematical representation of the terminal voltage, it
is often convenient to find an equivalent electric circuit model for simulation
and analysis of a battery cell. In the following, several such electric circuit
models representing an electrochemical cell are discussed, starting with a
basic model derived from the Hartley model.
dvd (t) 1 1
= i(t) − vd (t)
dt Cd Cd Rd
Battery Energy Storage 135
Cs Cd
RΩ
– vs(t) + – vd(t) +
i(t)
Rsd Rd
E0
FIGURE 4.11
Electric equivalent circuit battery model.
where vd(t) is the voltage dropped across the RdCd parallel circuit that is pro-
portional to the diffusion charge qd(t). Additional RC circuit elements can be
added to represent the diffusion charge, but we will keep it as a single RC
time constant for our simple model shown in Figure 4.11.
Another key cell dynamic that needs to be modeled is the effect of SoC on
the terminal voltage of the cell. Figure 4.7 showed how the battery terminal
voltage decreases as the cell is being discharged. In the middle of the charac-
teristics, the terminal voltage decrease is approximately linear which can be
modeled by a series capacitor Cs to represent the stored charge in the cell. The
voltage across this storage capacitor Cs is proportional to the stored charge
qs(t). As the SoC of the cell increases or decreases during charging or dis-
charging, the voltage across the capacitor will increase or decrease, respec-
tively. Additionally, an electrochemical cell loses charge while it is at rest. A
resistor can be added in parallel to the storage capacitor to account for this
loss of charge. This resistor Rsd represents the self-discharge of the cell. The
CsRsd circuit elements representing the storage capacitor and self-discharge
resistor are shown in Figure 4.11 in series with the diffusion parameters. The
mathematical representation of this segment of the circuit model in relation
to the terminal current is
dqs (t) 1
= i(t) − qs (t)
dt Rsd
The two other parameters that need to be added to complete the electro-
chemical cell equivalent circuit is a voltage source in series with a resistor
representing the ohmic resistance drop described in Section 4.4.6. The voltage
source is taken to be the open circuit voltage of the cell E0, and RΩ is the ohmic
resistance, both of which are shown in Figure 4.11 in series with the storage
and diffusion parameters. This completes the simple equivalent circuit model
of an electrochemical cell. The values of these circuit elements can be deter-
mined experimentally by applying a step change in battery current. The pro-
cedure for obtaining the parameters of this cell is given in Example 4.1.
136 Electric and Hybrid Vehicles: Design Fundamentals
Example 4.1
Solution
ΔVd, ΔVCs, ΔVRΩ are the voltage differences that need to be calculated from the
test data to obtain the diffusion, storage, and series resistance parameters, respec-
tively. The time to reach 63% of ΔVd is 100 s. Neglecting the self-discharge of the
cells, calculate the battery equivalent circuit parameters.
Let us first calculate the equivalent series resistance of the battery. The volt-
age drop for the series resistance shows up in the output voltage characteristics
as an instantaneous increase or decrease of the terminal voltage due to the step
change in current. The voltage increase due to the 15 A step change in current is
ΔVRΩ = 5.775 − 5.58 = 0.195 V. Therefore, the series resistance value is
∆VRΩ 0.195
RΩ = = = 0.013 Ω
∆I 15
∆Vd 0.098
Rd = = = 0.00653 Ω
∆I 15
The RC time constant for the diffusion parameters is 100 s. Therefore, the diffusion
capacitor Cd can be calculated as
FIGURE 4.12
Test data for a battery to calculate equivalent circuit parameters.
Battery Energy Storage 137
100
Cd = = 15, 306 F
0.00653
The storage capacitor Cs can be calculated from the voltage change due to the
constant current discharge ΔVCs and the change in stored charge. This is calculated
as follows:
∆Q 15 ( 4370 − 3150)
Cs = = = 871, 428.6 F
∆VCs 5.894 − 5.873
Ccapacity = 3600 ⋅ QC ⋅ k1 ⋅ k 2
E0 (VSoC)
Rsd
+ Cslow Cfast Vt
–
Ibat
–
FIGURE 4.13
Run-time battery model proposed by Mora et al. [13].
138 Electric and Hybrid Vehicles: Design Fundamentals
Idischarge (A)
Short time
OCV (V)
Vt (V)
constant
Idischarge
FIGURE 4.14
Example experimental curves to obtain run-time model parameters: (a) SoC vs. open circuit
voltage characteristics; (b) discharge plot for calculating RC time constants.
where
Qc is the battery capacity in A h
k1 an k2 are cycle number and temperature-dependent correction param-
eters, respectively
R Lseries
ZAC
+
E0 (SoC)
+
–
Vt
Ibatt
–
FIGURE 4.15
Impedance-based equivalent electric circuit battery model.
Battery Energy Storage 139
In the equivalent circuit, the diffusion process has been described by a gen-
eral CPE; open circuit voltage and concentration polarization have been
represented by Nernst equations; charge transfer polarization has been
represented by Tafel equation; ohmic voltage drop has been represented by
resistance RΩ; electric double layer has been represented by capacitance Cdbl;
and resistance Rsd represents self-discharge of the cell.
The first principle model is construction and chemistry independent. The
parameters of the model can be derived from experimental response data of
the device, which eliminates the need for the knowledge of electrochemical
if Ec R Vt
V1
Ect +–
Zero initial i
conditions
id
K if
Sq Ce +
– E0 Cdbl
Rsd Co
FIGURE 4.16
First principle battery model with constant current source. (From Xia, L., Behavioral modeling
and analysis of galvanic devices, PhD dissertation, University of Akron, Akron, OH, 2000.)
140 Electric and Hybrid Vehicles: Design Fundamentals
properties and the design details. The model is based on the fact that before
the discharge of any current, the internal voltage E0, the double-layer capaci-
tor voltage V1, and the terminal voltage Vt (the variables are shown in Figure
4.16) are all the same. The charge transfer potential and the concentration
polarization potential are zeroes for this condition. When a load is connected
to the terminals, initially the discharge current is almost entirely supplied by
the double-layer capacitor. As the double-layer capacitor discharges and V1
decreases, the charge transfer potential is established and the faradaic cur-
rent if starts to increase. When current if increases to a point where Ect does
not change appreciably, id becomes minimal. In this situation, the faradaic
current if supplies the majority of the load current. The potential drop in this
region is primarily due to the ohmic resistance.
As an example, the parameters for a generic battery cell are given below [7]:
of constant current discharge experiments give the I vs. tcut data for differ-
ent constant current levels; tcut is the time when the terminal voltage reaches
the cut-off voltage limit Vcut during constant current discharge. The data
obtained is used to fit Peukert’s equation to develop the model as
I ntcut = λ (4.35)
where
I is the constant discharge current
n and λ are curve fitting constants of a particular battery
n is a number between 1 and 2 with the value approaching 1 for smaller cur-
rents, but tends toward 2 for larger currents. The model does not specify the
initial capacity, nor does it model the voltage variation or temperature, and
aging factors. Peukert’s model does not give any terminal i–v information.
Example 4.2
Find the curve fitting constants n and λ for Peukert’s equation for the two mea-
surements available from a constant current discharge experiment of a battery: (1)
(t1,I1) = (10,18) and (2) (t2,I2) = (1,110).
Solution
Equation 4.35 is the Peukert’s empirical formula using the constant current dis-
charge approach. Taking logarithm of both sides of Equation 4.35
( )
Log10 I n × tcut = Log10 ( λ )
1 1
=> Log10 (I ) = Log10 (tcut ) + Log10 ( λ )
n n
Comparing with the equation for a straight line, y = mx + b; I versus tcut curve is
linear on a log-log plot, as shown in Figure 4.17.
The slope of the straight line is
log(t1 /t2 )
Therefore, n = − .
log(I1 /I2 )
−1
For the graph shown, n = = 1.27 [∵ t1 = 10t 2 ]
18 /110
The other constant can now be calculated from Peukert’s equation
λ = 1101.27 × 1 = 391.4 A h
142 Electric and Hybrid Vehicles: Design Fundamentals
I
(Log axis)
100 A (t1,l1)
y = mx + b
y = log(I )
x = log(tcut)
–1/n
18 A (t2,l2)
10 t (Log axis)
cut
FIGURE 4.17
Plot of Peukert’s equation using constant current discharge.
or
λ = 181.27 × 10 = 392.8 A h
λ
E = E0 − iR − Ki
λ − it
where
E is the battery voltage
i is the current
t is the time
The parameters of the model are E0, R, K, and λ representing battery refer-
ence voltage, internal resistance, polarization constant, and reference capac-
ity, respectively. The model parameters have some physical meaning and
relate the electrochemical behavior with the terminal i–v characteristics of
the battery. The output response is expressed as a function of time in this
model; however, the model is difficult to use for discharge patterns other
than constant current discharge.
battery can be used to predict the range of an electric vehicle. Using Peukert’s
equation, we can establish the relationship between Q and I. The practical
capacity of a battery is
Q = I × tcut
Q
=> tcut =
I
Q
In = λ
I
λ
=> Q = n −1
I
∫
SoD = i(τ)dτ
and
SoD
DoD =
Q(i)
SoD is the amount of charge that the battery generates to the circuit. Assume
that at t = t0, the battery is fully charged. Let us consider a small interval of
time dt. Therefore,
d(SoD)
d(DoD) = , where d(SoD) = i(t)dt
Q(i)
We know that Q = λ/In−1 for constant current discharge. Let, Q = λ/in−1 for
time-varying current as well, for the lack of anything better.
Therefore,
idt in
d(DoD) = n −1
= dt
λ/i λ
144 Electric and Hybrid Vehicles: Design Fundamentals
Integrating, we obtain,
t t
in
∫
t0
d(DoD) =
∫
t0
λ
dt
t
in
=> DoD(t) − DoD(t0 ) =
∫
t0
λ
dt
t n
i
DoD(t) =
λ
0t
∫
dt × 100%
(4.36)
Example 4.3
The constant current discharge characteristics of the battery pack used in an elec-
tric vehicle are
The current drawn from the battery during test drives of the electric vehicle for the
SAE schedule J227a has the profile shown in Figure 4.18. The current magnitudes
for the three SAE schedules are given in Table 4.3.
Find the range of the electric vehicle for each of the three schedules.
Solution
Apply the FDM (Equation 4.36) to find the number of driving cycles for DoD = 100%.
From FDM
Ia
Current
Ib
ta tb tc
Time (s)
FIGURE 4.18
Pattern of current drawn from the battery.
Battery Energy Storage 145
Therefore,
t100%
i1.35
1=
∫
0
645 × 3600
dt
72
i1.35
DoD for 1cycle => fcyc =
∫ 645 × 3600 dt
0
( )
72 19 1.35 38
i1.35
∫ dt = 4.31× 10 −7 100t
∫ ∫ (35) dt
1.35
=> fcyc = dt +
645 × 3600
0 0 19 19
1 2.35
= 4.31× 10 −7 9.41
19 + 121.5 (38 − 19)
2 . 35
=> fcyc = 2.74 × 10 −3
1
∴ 1 = N × fcyc => N =
fcyc
1
∴N = = 365 cycles
2.74 × 10 −3
Example 4.4
The data given in Table 4.4 is collected from an experiment on a battery with
mass 15 kg. The data is used to draw the Ragone plot shown in Figure 4.19. Using
the data points (8,110) and (67.5,10), calculate the constants of Peukert’s equation
n and λ.
Given a battery terminal power profile p(t), the specific power SP(t) profile can be
obtained by diving the power profile p(t) by the total vehicle mass mV (Figure 4.20).
The battery is assumed to be fully charged at t = 0.
Let, fr(t) = fraction of available energy provided by battery from 0 to t, where
fr(0) = 0, since SoD(0) = 0. Now, consider the time interval dt over which a fraction
of available energy dfr is provided by the battery
dE dE /mV d (SE)
dfr = = = .
Eavail dEavail /mV SEavail
TABLE 4.4
Data from Constant Power Discharge Test
P (W) tcut (h) EP (W h) SP (W/kg) SE (W h/kg)
(Measured) (Measured) (Calculated) (Calculated) (Calculated)
150 6.75 (150)(6.75) = 1013 150/15 = 10 1013/15 = 67.5
450 0.85 381 30 25.4
900 0.23 206 60 13.7
1650 0.073 120 110 8
Battery Energy Storage 147
SP
(W/kg)
(8,110)
100
60
(67.5,10)
30
SE (W h/kg)
8 13.7 25.4 67.5 100
FIGURE 4.19
Ragone plot for Example 4.4.
0 t 0 t
FIGURE 4.20
Power and specific power.
If dE is the energy provided by battery to the electrical circuit over dt and Eavail is
the total available energy, then
dE = pdt
d (SE) = (SP)dt
Therefore,
SEavail = f (SP)
We will use Peukert’s equation to relate specific power and specific energy as
follows:
(SP)n * SEavail = λ
Therefore,
SP (SP)n+1
dfr = dt = dt
λ /(SP)n λ
148 Electric and Hybrid Vehicles: Design Fundamentals
Integrating,
fr (t ) t
(SP)n+1
∫
fr ( 0 )
dfr =
∫
0
λ
dτ
t
(SP)n+1
=> fr (t ) =
∫
0
λ
dτ (4.38)
Equation 4.38 is the FDM using power density approach. If t = time at which x% of
available energy has been used, then
t
x (SP)n+1
100
=
∫
0
λ
dτ
(SP)n +1
t100%
Note that 1 =
0∫ λ
dτ
At t100%, 100% i.e., all the available energy has been used by the system.
TABLE 4.5
Properties of Electric and Hybrid Electric Vehicles Batteries
Specific
Energy Specific Energy
Battery Type (W h/kg) Power (W/kg) Efficiency (%) Cycle Life
Lead-acid 35–50 150–400 80 500–1000
Nickel-cadmium 30–50 100–150 75 1000–2000
Nickel-metal hydride 60–80 200–400 70 1000
Aluminum-air 200–300 100 <50 Not available
Zinc-air 100–220 30–80 60 500
Sodium-sulfur 150–240 230 85 1000
Sodium-nickel-chloride 90–120 130–160 80 1000
Li-polymer 150–200 350 Not available 1000
Li-ion 90–160 200–350 >90 >1000
4.6.1 Lead–Acid Battery
The lead-acid batteries have been the most popular choice of batteries for
electric vehicles during the initial development stages. The lead-acid battery
has a long history that dates back to the middle of the nineteenth century
and is currently a very mature technology. The first lead-acid battery was
produced as early as in 1859. In the early 1980s, over 100 million lead-acid
batteries were produced per year. The long existence of the lead acid battery
is due to
Cell 1 Cell 2
+ + + + – – – – –
– – – – – + + + +
FIGURE 4.21
Schematic diagram of lead-acid battery showing through-partition connection.
Battery Energy Storage 151
FIGURE 4.22
A lead-acid battery grid.
The lower practical cell voltage between 1.2 and 1.3 V means that more cells
have to be connected in series to get the desired voltage. The specific energy
of NiCd batteries is 30–50 W h/kg, which is similar to that of lead-acid
batteries.
The biggest drawbacks of NiCd batteries are the high cost and the toxicity
contained in cadmium. The environmental concerns may be overcome in
the long run through efficient recycling, but the insufficient power delivered
by the NiCd batteries is another important reason for not considering these
batteries for electric and hybrid electric vehicles applications. The draw-
backs of the NiCd batteries led to the rapid development of NiMH batter-
ies, which are deemed more suitable for electric and hybrid electric vehicle
applications.
4.6.4 Li–Ion Battery
The lithium metal has high electrochemical reduction potential relative to
that of hydrogen (3.045 V) and the lowest atomic mass (6.94), which shows
promise for a battery of 3 V cell potential when combined with a suitable
positive electrode. The interest in secondary lithium cells soared soon
after the advent of lithium primary cells in the 1970s, but the major dif-
ficulty was the highly reactive nature of the lithium metal with moisture
that restricted the use of liquid electrolytes. The discovery in late 1970s by
researchers at Oxford University that lithium can be intercalated (absorbed)
into the crystal lattice of cobalt or nickel to form LiCoO2 or LiNiO2 paved the
way toward the development of Li-ion batteries [17]. The use of metallic-Li is
bypassed in Li-ion batteries by using lithium intercalated (absorbed) carbons
(LixC) in the form of graphite or coke as the negative electrode along with the
lithium metallic oxides as the positive electrode. The graphite is capable of
hosting lithium up to a composition of LiC6. The majority of the Li-ion bat-
teries use either a layered oxide or iron phosphates of lithium as the positive
electrode. The layered positive electrodes of cobalt oxide are expensive, but
proved to be the most satisfactory. Nickel-oxide LiNiO2, which costs less,
can also be used, but is structurally more complex. The performance is simi-
lar to that of cobalt-oxide electrodes. The manganese oxide-based positive
electrodes (LiMn2O4 or LiMnO2) are also used since manganese is cheaper,
widely available, and less toxic. Alternative positive electrode material is the
lithium-iron-phosphate (LiFePO4) which can deliver stable and good perfor-
mance at lower costs.
Battery Energy Storage 153
RL
Electron flow
(discharge)
+ –
e– Li Li+ e–
Li+
LiCoO2 Carbon
Li Li+ Li+
e– e–
Li Li+ Li+
Electrolyte
e– Li Li+ Li+ e–
FIGURE 4.23
Lithium-ion cell. (Courtesy of Research Studies Press Ltd.)
Discharge
Li x C6 ← → 6C + xLi + + xe − where 0 < x < 1
Charge
Discharge
xLi + + xe − + Li(1− x )CoO 2 ← → LiCoO 2
Charge
During cell charge operation, the lithium ions move in the opposite direction
from the positive electrode to the negative electrode. The nominal cell volt-
age for a Li-ion battery is 3.6 V, which is equivalent to three NiMH or NiCd
battery cells.
The lithium ion batteries have high specific energy, high specific power,
high energy efficiency, good high-temperature performance, and low self-
discharge. The components of Li-ion batteries are also recyclable. These
characteristics make Li-ion batteries highly suitable for electric and hybrid
vehicles and other applications of rechargeable batteries. The main draw-
back of Li-ion batteries is that these are very sensitive to overvoltages and
154 Electric and Hybrid Vehicles: Design Fundamentals
4.6.5 Li–Polymer Battery
The Li-polymer evolved out of the development of solid-state electrolytes,
i.e., solids capable of conducting ions but are electron insulators. The solid-
state electrolytes resulted from the research in the 1970s on ionic conduction
in polymers. These batteries are considered solid-state batteries since their
electrolytes are solids. The most common polymer electrolyte is the polyeth-
ylene oxide compounded with an appropriate electrolyte salt.
The most promising positive electrode material for Li-polymer batteries
is vanadium oxide V6O13 [16]. This oxide interlaces up to 8 lithium atoms per
oxide molecule with the following positive electrode reaction:
Discharge
Li x + V6 O13 + xe − ← → Li x V6 O13 where 0 < x < 8
Charge
The Li-polymer batteries have the potentials for the highest specific energy
and power. The solid polymers, replacing the more flammable liquid elec-
trolytes in other type of batteries, can conduct ions at temperatures above
60°C. The use of solid polymers also has a great safety advantage in case of
electric and hybrid electric vehicles accidents. Since the lithium is interca-
lated into carbon electrodes, the lithium is in ionic form and is less reactive
than pure lithium metal. The thin Li-polymer cell gives the added advan-
tage of forming a battery of any size or shape to suit the available space
within the electric and hybrid electric vehicles chassis. The main disadvan-
tage of the Li-polymer battery is the need to operate the battery cell in the
temperature range of 80°C–120°C. Li-polymer batteries with high specific
energy, initially developed for electric vehicle applications, also have the
potential to provide high specific power for hybrid electric vehicle applica-
tions. The other key characteristics of the Li-polymer are good cycle and
calendar life.
100 W/kg at 80% DoD (see Sections 4.3.8 and 4.3.12 for definition of depth of
discharge and specific power, respectively). With the present-day technol-
ogy, the range of zinc-air batteries can be between 300 and 600 km between
recharges.
Other metal air systems have also been investigated but the work has been
discontinued due to severe drawbacks in the technologies. These batteries
include iron-air and aluminum-air batteries where iron and aluminum are
respectively used as the mechanically recyclable negative electrode.
The practical metal-air batteries have two very attractive positive features:
(1) The positive electrode can be optimized for discharge characteristics,
since the batteries are recharged outside the battery and (2) the recharging
time is rapid with a suitable infrastructure.
– +
NaAlCl4
Beta-alumina separator
Liquid sodium
Cell case
FIGURE 4.24
A sodium-nickel-chloride cell.
is similar to that of NaS battery, but the positive sulfur electrode is replaced
by nickel chloride (NiCl2) or a mixture of nickel chloride and ferrous chlo-
ride (FeCl2). The negative electrode and the electrolyte are the same as in
NaS battery. The schematic diagram of a NaNiCl2 cell is shown in Figure
4.24. In order to provide good ionic contact between the positive elec-
trode and the electrolyte, both of which are solids, a second electrolyte of
sodium chloraluminate (NaAlCl4) is introduced in a layer between NiCl2
and beta-alumina. The NaAlCl4 electrolyte is a vital component of the bat-
tery, although it reduces the specific energy of the battery by about 10%
[17]. The operating temperature is again high, similar to that of NaS battery.
The basic cell reactions for the nickel chloride and ferrous chloride positive
electrodes are
Discharge
NiCl 2 + 2Na ← → Ni + 2NaCl (2.58 V)
Charge
Discharge
FeCl 2 + 2Na ← → Fe + 2NaCl (2.35 V)
Charge
The sodium chloride batteries are commonly known as the ZEBRA bat-
teries, which originally resulted from a research collaboration between sci-
entists from the United Kingdom and South Africa in the early 1980s. The
ZEBRA batteries have been shown to be safe under all conditions of use. The
ZEBRA batteries have high potentials for being used as batteries for electric
vehicles and hybrid electric vehicles. There are several test programs per-
formed with the ZEBRA batteries.
TABLE 4.6
USABC Objectives for EV Advanced Battery Packs
Minimum Goals for
Long-Term Long-Term
Parameter Commercialization Goals
Specific energy (W h/kg) 150 200
(C/3 discharge rate)
Specific power (W/kg) 300 400
(80% DoD per 30 s)
Specific power (W/kg), 150 200
Regen. (20% DoD per 10 s)
Recharge time (h) 4–6 3–6
(20% → 100% SoC)
Cost, U.S. $/kW h 150 100
158 Electric and Hybrid Vehicles: Design Fundamentals
TABLE 4.7
USABC Goals for HEV Advanced Energy Storage Systems
Power Assist Power Assist
Parameter (Minimum) (Maximum)
Pulse discharge power, 10 s (kW) 25 40
Peak regenerative pulse power, 10 s (kW) 20 35
Total available energy (kW h) 0.3 at C/1 rate 0.5 at C/1 rate
Maximum weight (kg) 40 60
Maximum volume (L) 32 45
Cost, @1,000,000 units/year (U.S. $) 500 800
The USABC has also set goals for hybrid electric vehicles at two levels of
power-assist, one at the 25 kW level and the other for the 40 kW level. A sub-
set of the goals set by the USABC for power-assist hybrid electric vehicle
energy storage system is listed in Table 4.7. The calendar life for these bat-
teries is targeted for 15 years, while cycle life has been set for 300,000 cycles
for specified SoC increments. The roundtrip energy efficiency has been set to
90% for both 25 and 40 kW power-assist hybrid electric vehicles.
The USABC has also specified goals for two main PHEV battery types:
a high power/energy ratio battery providing 10 mi of all-electric range
(PHEV-10), and a low power/energy ratio battery providing 40 mi of
all-electric range (PHEV-40). PHEV-10 goals are set for a “crossover utility
vehicle” weighing 1950 kg and the PHEV-40 goals are set for a midsize sedan
weighing 1600 kg. Few of these important goals set for PHEV development
are listed in Table 4.8. The calendar life for these batteries is also targeted
for 15 years, and roundtrip energy efficiency has been set to 90% for both
PHEV-10 and PHEV-40. All the specified goals for energy storage systems for
TABLE 4.8
USABC Goals for PHEV Energy Storage Systems
Parameter PHEV-10 PHEV-40
Pulse discharge power, 10 s (kW) 45 38
Peak regenerative pulse power, 10 s (kW) 3,025 35
Available energy for charge depleting 3.4 11.6
mode (kW h)
Available energy for charge sustaining 0.3 0.5
mode (kW h)
Charge depleting life/discharge 5,000/17 5,000/58
throughput (cycles/MW h)
Charge sustaining life cycle (cycles) 300,000 300,000
Maximum weight (kg) 60 120
Maximum volume (L) 40 80
Cost, @100,000 units/year (U.S. $) 1,700 3,400
Battery Energy Storage 159
electric and hybrid electric vehicles are listed in the USABC Web site under
USCAR at www.uscar.org.
FIGURE 4.25
Parameter estimations and pack management in a BMS.
BMS for an energy storage system is designed to have all or some of the fol-
lowing features:
and estimated cell parameters with references or neighboring cells. The volt-
ages and SoC anomalies of one cell compared to the nearby cells is indicative
of poor health of that cell. Similarly, the excessive self-discharge of the pack
raises a flag, and is compared with preset limits to estimate the pack SoH.
The SoH information can be used to replace damaged cells in a pack instead
of replacing the entire pack.
The temperature is the primary environmental factor that affects the SoC
of an energy storage system. Imbalances in temperature among the various
cells in a pack will result in imbalances in the SoC. Additionally, tempera-
ture affects the self-discharge rates. The thermal management in a pack is
part of the cooling system design for the pack, but the temperature infor-
mation of the cells should be effectively utilized for protection and health
monitoring of the cells.
The maximum power available from the battery at a given time is calculated
in the BMS based on the SoC and terminal voltages, ensuring that operating
voltage, current, SoC, and other design limits are not violated. The BMS sets
the power limits during charging and discharging for battery protection.
Batteries could get severely damaged due to inappropriate charging. The
limits are reported to the supervisory controller for powertrain controls in
electric and hybrid vehicles.
There are three levels of management systems: pack-level management,
modular-pack-level management, and cell-level management. Pack-level
management is the most basic one where overall pack voltage and SoC are
monitored, whereas the most complete cell equalization and balancing is
possible when individual cell parameters of voltage, current, and tempera-
ture are monitored. The charging and discharging power managements at
the pack level leaves individual cells vulnerable to damage. In modular-pack-
level management, groups of cells are treated as a module for cell balancing
and equalization; the BMS algorithms depend on group voltage, current,
and temperature measurements rather than pack or individual cell measure-
ments. For packs employing cell equalizer circuits, the BMS generates com-
mand signals for cell voltage equalization based on its measurements and
estimations. The circuits act on these signals to balance the cells or groups
of cells.
Data logging is another important function of the energy storage manage-
ment system. The data for voltage, current, temperature, SoC, and number of
charge/discharge cycles could be stored as a function of time for SoH moni-
toring, diagnostics, and fault analysis.
∫ SoC
SoC
estimation
i(t)
+
Battery Electric
Vbatt
pack drivetrain
–
FIGURE 4.26
Battery SoC measurement.
SoC(t) = CP − SoD(t)
states that both charge and voltage need to be balanced in a series string to
maximize the output of a pack. When a series-string of electrochemical cells
is charged as a pack, slight parameter mismatches in individual cells and
temperature differences result in charge and voltage imbalances. The imbal-
ances adversely affect the vehicle performance by reducing the throughput
of the battery pack.
The chemical reactions in an electrochemical cell depend on the tempera-
ture and pressure. The temperature differences among the cells change the
self-discharge rates, causing imbalances in the charge of the cells. A low cell
temperature reduces chemical activity which increases the cell’s internal
impedance. The increased internal resistance reduces the terminal voltage,
and thus the cell capacity. In addition, manufacturing differences and differ-
ent aging characteristics result in parameter mismatches among individual
cells, which cause voltage and capacity imbalances [20,21].
The charge imbalance also shows up as voltage differences. The imbalances
tend to grow as the pack goes through repeated charge/discharge cycles.
The weaker cells tend to charge slower and the stronger cells charge faster.
The process shortens the pack life and reduces its utilization. The number of
charge/discharge cycles affects some battery chemistry more than the oth-
ers. For example, Li-ion batteries are highly sensitive to overvoltages and
undervoltages. Li-ion batteries are recommended to limit the charge/dis-
charge rates to no more than 2 C, and also to keep the cells charged to at least
40% SoC to minimize aging.
The maximum throughput of the pack can be ensured by balancing the
voltage and charge of individual cells. The cell balancing methods utilize
electronic circuits and control to even out the voltages and SoC of a series-
string of electrochemical cells. The simplest strategy adopted for charging a
series-string of cells is to monitor the cell voltages and discontinue charging
when one of the cells (strongest cell) reaches the voltage limit for individual
cells. Extended charging is another option where charging is continued even
after the strongest cell has reached its capacity to bring the weaker cells up
to capacity. When charging continues to bring the weaker cells to the maxi-
mum voltage, overvoltage results in the stronger cells. Overcharging is not
at all an option with certain battery chemistry, while in others the process
vents hydrogen gas (known as gassing) and removes water from the over-
charged cells. Increased gassing in the cell at elevated temperatures shortens
the cell life.
The overcharging in the stronger cells can be avoided if there is a path
to shunt the charging currents once they reach the voltage limit. Similarly,
the simplest protection during discharging of a pack is to shut down when the
first cell reaches the minimum voltage limit. This cell is consequently the
weakest cell in the series-string and is limiting the capacity of the pack. If
discharging is continued to extract energy from the stronger cells, then the
weaker cell voltage will fall below the minimum voltage level, possibly caus-
ing damage to the cells.
164 Electric and Hybrid Vehicles: Design Fundamentals
4.7.4 Battery Charging
The charging of secondary batteries is accomplished in several phases
using different charging currents. The phases are structured based on the
battery chemistry to minimize the damage on the cells. The initial charg-
ing phase is the bulk charging phase when the cells are charged with the
maximum current to replenish most of the charge lost during discharge.
The last few percentages of SoC are replenished with the absorption charg-
ing phase. The charging current in this phase is kept low to prevent any
damage to the cells. An equalization phase can also be used to fully charge
and balance all the cells in the battery pack. The float charge phase starts
once the battery is fully charged to compensate for energy lost over time
due to self-discharge. Microprocessor controllers are used to set the charg-
ing profile based on an algorithm to tune the charging for a particular type
of battery chemistry.
The battery charging circuits can apply either a constant current or a con-
stant voltage or any combination of the two to design the charge profile. In
the constant current charging method, known as I-charging, a current regu-
lator in the battery charger maintains the set current level. The charging
current levels are adjusted by the current regulator for the different phases
of charging. The charging current can also be applied in the form of pulses
by pulse width modulation (PWM) control of the output voltage. The charg-
ing rate is controlled by adjusting the pulse width. The short durations
between the pulses allow the chemical reactions within the cells to stabilize.
Excessive chemical reactions that could lead to gassing are avoided by using
pulse charging. An example of multistep charging profile with I-charging
and PWM control is shown in Figure 4.27 [22].
For constant voltage charging, a voltage greater than the battery upper
limit voltage is applied by the charger for bulk charging. A constant volt-
age charge is also usually applied during the absorption charging phase.
During the float charging phase, the charger applies a DC voltage slightly
lower than the battery upper limit voltage across the battery. A slight drop
in the battery voltage results in charge being replenished in the battery. This
form of charging is also known as trickle charging used to compensate for
self-discharge in the cells.
Battery Energy Storage 165
ICH
Bulk
charging
Absorption
1.00 charging
Equalization Float
0.75
charging charging
0.50
0.25
FIGURE 4.27
Multistep I-charging with PWM control.
Problems
4.1 Estimate the weight of a 12 V, 100 A h lead-acid battery. Do this by calcu-
lating the reactant masses participating in the overall chemical reaction.
Also, assume that the mass of H2O in the electrolyte is twice the mass
of H2SO4. Neglect battery casing mass, electrode grid mass, separator
mass, and current bus mass. (Note that n = 2 for Pb and PbO2 and n = 1
for H2SO4.)
4.2 In the nickel-cadmium cell, nickel oxyhydroxide, NiOOH is the active
material in the charged positive plate. During discharge it reduces to
the lower valence state, nickel hydroxide Ni(OH)2, by accepting electrons
from the external circuit:
Discharge
2NiOOH + 2H 2O + 2e − ← → 2Ni(OH)2 + 2OH − (00.49 V)
Charge
Discharge
Cd + 2OH − ← → Cd(OH)2 + 2e − (0.809 V)
Charge
Discharge
Cd + 2NiOOH + 2H 2O ← → 2Ni(OH)2 + Cd(OH)2− (1.299 V)
Charge
Estimate the weight of a 11.7 V, 100 A h Ni-Cd battery. Neglect the mass
KOH component of the electrolyte.
4.3 A 12 V battery is connected to a series RL load, as shown in Figure P4.3.
The battery has a rated capacity of 80 A h. At t = 0, the switch is closed
and the battery begins to discharge.
t=0 R
i(t) 90 mH
Battery +
12 V –
FIGURE P4.3
(a) Calculate and plot the battery discharge current, i(t), if the steady-
state discharge rate is C/2. Neglect battery voltage drop.
(b) Calculate and plot SoD(t) in A h for 0 < t < 2 h.
(c) Calculate and plot SoC(t) assuming that at t = 0, the battery is
charged to rated capacity. Assume also that the rated capacity is
the practical capacity.
(d) Calculate the time corresponding to 80% DoD.
SP (W/kg) SE (W h/kg)
10 67.5
110 8
100 A
0.5 m 1 ms t
80 A
FIGURE P4.5
I
32 100
v (km/h)
35
0 19 38 42 4 t (s) 0 19 38 42 4 t (s)
(a) (b)
FIGURE P4.6
168 Electric and Hybrid Vehicles: Design Fundamentals
References
1. D. Linden, Handbook of Batteries, 2nd edn., McGraw-Hill Inc., New York, 1995.
2. C. Mantell, Batteries and Energy Systems, McGraw-Hill Inc., New York, 1983.
3. K. Kordesch and G. Simadar, Fuel Cells and Their Applications, John Wiley & Sons
Inc., Chichester, U.K., 1996.
4. J.O.M. Bockris, Fuel Cells: Their Electrochemistry, McGraw-Hill Inc., New York,
1969.
5. J.S. Newman, Electrochemical Systems, Prentice Hall, Englewood Cliffs, NJ, 1991.
6. T. Erdey-Gruz, Kinetics of Electrode Process, John Wiley & Sons Inc., New York,
1972.
7. L. Xia, Behavioral modeling and analysis of galvanic devices, PhD dissertation,
University of Akron, Akron, OH, 2000.
8. A.J. Bard and L.R. Faulkner, Electroanalytical Methods: Fundamentals and
Applications, John Wiley & Sons Inc., New York, 1996.
9. J.R. Elliott and C.T. Lira, Introductory Chemical Engineering Thermodynamics,
Prentice Hall, Upper Saddle River, NJ, 1999.
10. T.T. Hartley, Insights into the fractional order initial value problem via semi-
infinite systems, NASA TM-1998-208407, November 1998.
11. R.L. Hartman II, An Aging Model for Lithium-ion Cells, PhD dissertation,
University of Akron, Akron, OH, 2008.
12. T.T. Hartley and A.G. Jannette, A first principles model for nickel-hydrogen
batteries, AIAA (American Institute of Aeronautics and Astronautics), in
3rd International Energy Conversion Engineering Conference, San Francisco, CA,
August 2005.
13. M. Chen and G.A. Rincon-Mora, Accurate electrical battery model capable of
predicting runtime and I–V performance, IEEE Transactions on Energy Conversion,
21(2), 504–511, June 2006.
14. S. Abu-Sharkh and D. Doerffel, Rapid test and non-linear model characteriza-
tion of solid-state lithium-ion batteries, Journal of Power Sources, 130, 266–274,
2004.
15. C.M. Shephard, Design of primary and secondary cells, Journal of the
Electrochemical Society, 112, 657–664, July 1965.
16. D.A.J. Rand, R. Woods, and R.M. Dell, Batteries for Electric Vehicles, John Wiley &
Sons Inc., New York, 1998.
17. R.M. Dell and D.A.J. Rand, Understanding Batteries, Royal Society of Chemistry,
Cambridge, U.K., 2001.
18. M.J. Isaacson, R.P. Hollandsworth, P.J. Giampaoli, F.A. Linkowsky, A. Salim,
and V.L. Teofilo, Advanced lithium ion battery charger, in The Fifteenth Annual
IEEE Battery Conference on Applications and Advances, Long Beach, CA, January
2000, pp. 193–198.
Battery Energy Storage 169
19. G.L. Plett, Extended Kalman filtering for battery managements systems of LiPB-
based HEV battery packs: Part 1. Background, Journal of Power Sources, 134,
252–261, 2004.
20. W.F. Bentley, Cell balancing considerations for lithium-ion battery systems, in
Twelfth Annual IEEE Battery Conference on Applications and Advances, New York,
January 1997, pp. 223–226.
21. R.M. Laidig and W.J. Wurst, Technology implementation of stationary battery
failure prediction, in Proceedings of the Ninth Annual IEEE Battery Conference on
Applications and Advances, Long Beach, CA, January 1994, pp. 168–172.
22. N.H. Kutkut and M.D. Divan, Dynamic equalization techniques for series bat-
tery stacks, in IEEE Telecommunications Energy Conference, Boston, MA, October
1996, pp. 514–521.
5
Alternative Energy Storage
171