Next Article in Journal / Special Issue
The Noisy and Marvelous Molecular World of Biology
Previous Article in Journal
Internet of Things in Smart Grid: Architecture, Applications, Services, Key Technologies, and Challenges
Previous Article in Special Issue
Full Statistics of Conjugated Thermodynamic Ensembles in Chains of Bistable Units
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Thermodynamic Model for Lithium-Ion Battery Degradation: Application of the Degradation-Entropy Generation Theorem

Mechanical Engineering Department, University of Texas at Austin, Austin, TX 78712, USA
*
Author to whom correspondence should be addressed.
Inventions 2019, 4(2), 23; https://doi.org/10.3390/inventions4020023
Submission received: 19 February 2019 / Revised: 25 March 2019 / Accepted: 27 March 2019 / Published: 3 April 2019
(This article belongs to the Special Issue Thermodynamics in the 21st Century)

Abstract

:
Presented is a lithium-ion battery degradation model, based on irreversible thermodynamics, which was experimentally verified, using commonly measured operational parameters. The methodology, applicable to all lithium-ion batteries of all chemistries and composition, combined fundamental thermodynamic principles, with the Degradation–Entropy Generation theorem, to relate instantaneous capacity fade (loss of useful charge-holding capacity) in the lithium-ion battery, to the irreversible entropy generated via the underlying dissipative physical processes responsible for battery degradation. Equations relating capacity fade—aging—to battery cycling were also formulated and verified. To show the robustness of the approach, nonlinear data from abusive and inconsistent battery cycling was measured and used to verify formulations. A near 100% agreement between the thermodynamic battery model and measurements was achieved. The model also gave rise to new material and design parameters to characterize all lithium-ion batteries.

Graphical Abstract

1. Introduction

Lithium-ion batteries are rechargeable and exhibit high-energy density, minimal maintenance, low self-discharge, and a long cycle life [1,2]. Lithium-ion battery issues include self-discharge, aging and thermal instability. Models to predict performance over time, use electrical parameters and theories, others combine electrical and chemical phenomena [1,2,3,4]. Although significant advances have been made for battery electrode and electrolyte materials [4,5,6,7] for high-energy density lithium-ion batteries, critical issues concerning safety, reliability, and continuous availability [5,8,9,10,11] limit their use. To mitigate loss of charge-storing capacity with cycling, and to prevent battery overheating, leading to catastrophic thermal events, manufacturers of lithium-ion battery-operated devices and electric vehicles often rely on computationally intensive monitoring algorithms [11,12], which are unreliable at predicting sudden battery instabilities/failures. Due to the complex internal mechanisms involved in the lithium-ion battery operation, battery manufacturers and researchers use expensive empirically fitted mathematical models, to estimate the battery’s response to loading [13,14,15,16]. Needed is a robust and universal model to accurately and consistently characterize the battery and predict degradation under all operating conditions, irrespective of active mechanisms or interactions. Presented is an entirely physics-based battery degradation model that instantaneously characterizes the battery’s response, using commonly measured battery parameters of voltage, current, and temperature.
The chemistries at the battery’s cathode [7]
Li 1 x MO 2 + xLi + + xe LiMO 2
and anode [7]
Li x C 6 C 6 + xLi + + xe
give an overall cell reaction [9]
Li 1 x MO 2 + Li x C 6 LiMO 2 + C 6
involving a nominal cell voltage of approximately 3.6–4.2 V, depending on the transition metal M used. Here, 0 ≤ x < 1. Transition metals include cobalt, manganese, and nickel, among others. In a typical Li-ion cell, the electrodes are active materials Li 1 x MO 2 and Li x C 6 , bonded to current collectors by the electrolyte, usually liquid, gel, or lithium metal polymer, which facilitates transport of lithium ions ( Li + ), between electrodes. During charging, Li + are deintercalated (extracted) at the cathode and the active material is oxidized, whereas the anode active material is reduced and Li + extracted from the cathode are intercalated (inserted) into the anode. The phenomena reverse for discharge [1,4,7].
Full charging applies a constant current which energizes the battery to its rated voltage, followed by a constant-voltage process during which charge current decreases to 1–3% of the battery’s rated current. Cycling a battery dissipates a variable amount of heat, depending on the discharge/charge rates and operating conditions. End of discharge voltage for typical commercial Li-ion batteries is 2.7 V/cell, to avoid damage from deep discharge. An intermediate step, settling, allows the transport phenomena and reaction kinetics to stabilize. Manufacturers specify nominal and typical electrical charge-holding capacity (discharge capacity) of Li-ion batteries at a given discharge rate (or current) in Ampere-hours (Ah).
Battery life, typically 500 to 1200 discharge–charge cycles, depends on many factors [8,9,10]. Even with proper maintenance, the capacity of a Li-ion cell fades with cycling. Several discharge–charge cycles produce more dissipative losses which generate heat [1]. Other degradation mechanisms include lithium corrosion and plating on the anode, and excessive growth or disconnection of the solid–electrolyte interface (SEI) layer on or from the anode, respectively, resulting in a loss of contact [10]. On the cathode, a passivation layer forms, which grows during cycling and reduces capacity over time. Structural changes in the electrodes and irreversible decomposition of the electrolyte, over time, also limit intercalation and diffusion of Li + .
Li-ion batteries are very sensitive to charge rates and Depth of Discharge DoD. Improper charging overheats batteries, causing catastrophic failure. Overcharging facilitates migration of Li + from the layered structure, building up metallic lithium on the anode, and releasing excess oxygen at the cathode. As this continues, pressure in the battery increases and more heat is released, which can eventually cause an explosion. Overdischarging causes similar irreversible damage. Battery capacity and life cycle also depend on design and operating conditions.

Background/Literature Review

In [16,17], lithium-ion battery degradation in electric vehicle application is presented. In both studies, the authors used empirically fitted Thevenin (for battery characterization) and DoD-dependent (for battery aging) models. Using representative real-world drive cycles, a close fit between model and experimental results was obtained [16], while noting the effects of regenerative braking and ambient temperature on model error. In [17], consistent cycling capacity fade test of a 200 Ah lithium-ion battery yielded a model accuracy greater than 96.5%, for all 400 cycles. Karnopp [18] coupled the chemical and electrical interactions of a battery via the Gibbs relation A   J = µ i N ˙ i = G ˙ where A is de Donder’s affinity, J is reaction rate, µ i are chemical potentials, Ni are numbers of moles, and G is Gibbs free energy. The dot above quantities represents differentiation with time t. The entropy produced S’ from the energy dissipated via the Ohmic or chemical reaction, was evaluated from Prigogine’s minimum entropy generation theorem [19] as.
S ˙ = A   J T = V I T
where T is temperature, V is battery voltage, I is battery current, and the apostrophe denotes entropy production. In [20], process irreversibilities in lithium-ion batteries were presented through overpotentials (charge transfer and diffusion). The authors gave a thermodynamic model using Gibbs energy change,
d G = S r e v d T + V d P + V d q
where S r e v = G T P = q V T was the reversible entropy from thermal energy change. The irreversible entropy produced by the Ohmic work dissipation was S ˙ = V I T and the entropy transfer out was S ˙ e x t = h A T T T . Menard et al. [3] gave a similar bond-graph model for one electrode of the lithium-ion battery.
While measurements of temperature-dependent degradation are readily available [10,15,21,22,23,24,25,26,27], models do not adequately estimate the effect of temperature on battery performance. Some use the Arrhenius equation to model temperature-dependent degradation at moderate to high temperatures, and empirical models for lower temperatures [15,22,23,24,25,26,27]. Battery degradation models, not valid for all instantaneous conditions, often fail under unsteady operation, deep discharge, or other nonlinear system interactions; often cannot accurately predict useful life; cannot adequately account for battery aging or parasitic losses; and cannot be easily adapted to other battery types. This study used thermodynamic formulations (of irreversible entropies produced as a consequence of energy dissipation) and the DEG theorem [28], to analyze nonlinear, uncontrolled and cyclically inconsistent lithium-ion battery degradation. Since these entropies in turn depend on those phenomenological parameters (or variables) used by the battery industry to describe degradation, the result was a battery degradation model, based on fundamentals, but expressed in terms of convenient industry parameters. The analyses were then experimentally verified. Definitions of relevant battery parameters were adapted from [1].

2. Degradation-Entropy Generation Theorem Review

Rayleigh [29] through his dissipation function was the first to characterize dissipative forces. In classical irreversible thermodynamics, Onsager [30] developed his famous reciprocity theorem [19]. A quantitative study of degradation of systems through dissipative processes, by Bryant et al. [28], formulated the Degradation-Entropy Generation (DEG) theorem, that established a direct relationship between degradation measure w (which in this article is battery charge content or capacity) and rates of entropy generation S i , using irreversible thermodynamics. Entropy measures disorganization in materials. Since degradation is advanced and permanent disorganization, entropy generation is fundamental to degradation.

2.1. Statement

Given an irreversible material transformation consisting of i = 1, 2, …, n dissipative processes p i , which could describe an energy, work, or heat characteristic of the process, and assuming the mechanism can be described by a parameter or state variable w that measures the effects of the transformation or degradation, such that
w = w p i = w p 1 , p 2 ,     ,   p n ,   i = 1 , 2 , , n
is monotonic in each p i . Then the rate of change of the degradation
w ˙ = i B i S ˙ i
is a linear combination of rates of entropies S ˙ i generated by the dissipative processes p i , where the transform process coefficients
B i = w S i p i
are slopes of degradation w with respect to entropy generation S i ; the | p i notation refers to the process p i being active.
The theorem was stated and proved in [28]. Integrating Equation (7) over time t, composed of cycles wherein B i is constant, yields the accumulated total degradation
Δ w = i B i S i
which is a linear combination of the accumulated entropies S i generated by the dissipative processes p i .

2.2. Generalized Degradation Analysis Procedure

Bryant et al. [28] presented a structured approach to degradation analysis which embeds the physics of the dissipative processes into the energies p i = p i ζ i j ,   j = 1 , 2 , , m ; derives entropy generation S ˙ i as a function of the p i and expresses the rate of degradation w ˙ , as a linear combination of the entropy generation terms, Equation (7). p i is the energy dissipated by internal processes, lost work, transferred heat, a change in the thermodynamic energy (internal energy, enthalpy, Helmholtz or Gibbs free energy), or some other functional form of energy. ζ i j are the time-dependent phenomenological variables used in practice and associated with p i . The degradation/transform coefficients B i must be measured using Equation (8). This approach has the following characteristics:
  • It identifies the degradation measure w, dissipative process energies p i , and phenomenological variables ζ i j ;
  • It finds the entropy generation S i = S i ( p i ) caused by the dissipative processes p i ;
  • It evaluates the coefficients B i by measuring increments, accumulation or rates of degradation versus increments, accumulation or rates of entropy generation;
  • It relates degradation measure to entropy generation, via Equations (7) or (9).
This approach can formulate degradation models consisting of one or many dissipative processes p i , acting simultaneously. The DEG theorem has been previously applied to friction and wear processes [31,32,33,34] and metal fatigue [35,36,37].

3. Formulations

3.1. Fundamental Thermodynamic Formulations

This section reviews the laws of thermodynamics [19,38,39,40,41,42,43,44].

3.1.1. First Law—Energy Conservation

The first law of thermodynamics
d U = δ Q δ W + µ k d N k
for a stationary thermodynamic system, neglecting gravity, dU balances the change in internal energy, δQ is the heat exchange across the system boundary, δW is sum of the work transfer across the system boundary, and µ k d N k is sum of the compositional energy changes. Inexact differential δ indicates path-dependent variables.

3.1.2. Second Law and Entropy Balance—Irreversible Entropy Generation

The Clausius inequality is the change in entropy for a closed system:
d S δ Q T
The inequality indicates the system’s entropy change dS equals or exceeds the measured entropy transfer by heat δQ/T, across the system boundary. For a reversible process, the entropy change
d S = d S r e v = δ Q r e v T
which approximates a quasi-static (very slow) process, in which total entropy change occurs via a reversible heat transfer δ Q r e v .
Expressed as the equality d S = d S e + δ S by Prigogine [19,42], the second law of thermodynamics equates dS to the measured entropy flow d S e across the system boundaries from heat transfer or mass transfer (for open systems), plus any entropy δ S produced through dissipative processes, within the system boundaries. Entropy generation δ S measures the structural disorganization from permanent changes to the system [19,40,42,45], when it evolves. For a closed system (no mass transfer across system boundaries), the second law becomes [19,42]
d S = d S i r r = δ Q T + δ S
where d S i r r is entropy change along an irreversible (real) path, δQ/T is entropy flow by heat transfer, which might be positive or negative, and T is the temperature of the boundary where the energy/entropy transfer takes place. The second law, via Equations (11) and (13), also asserts that entropy generation δ S 0 .

3.1.3. Combining First and Second Laws

Reversible Transformation: For a system undergoing quasi-static heat transfer and compression work, Equation (10), replacing δ Q with δ Q r e v = T d S r e v from Equation (12), gives the combined form of the first and second laws [44]
d U = d U r e v = T d S r e v P d V + µ k d N k
Here P is pressure, V is volume, and T is temperature. Equation (14) is the reversible change in internal energy between two equilibrium states, according to the first law, and is valid for all systems.
Irreversible Transformation: Eliminating δQ from Equation (10) via Equation (13) gives, for compression work P d V , [19,33,34,41,42]
d U = d U i r r = T d S e T δ S P d V + µ k d N k
the irreversible combined first and second laws. Equation (15) applies to irreversible processes. Comparing to the prior equation, reversible entropy change d S r e v was replaced by entropy flow d S e , and entropy generation δ S was included.

3.2. Li-Ion Battery Analysis

Battery degradation observed electrically via capacity fade, occurs chemically through electrode corrosion and evolution of gases, and thermally through hot surroundings and joule heating.
Engineering Model: The thermodynamic system boundary only encloses the battery; the system is closed (battery mass stays in the battery); heat transfers between battery and surroundings; and the system is at equilibrium before and after charging or discharging.

3.2.1. Combining First and Second Laws with Gibbs Potential

Electrochemical energy storage devices are conveniently characterized using the Gibbs free energy [19,40,42,45]
G = U + P V T S
which measures process-initiating energy changes in a thermodynamic system. Differentiating Equation (16) and substituting Equation (14) for dU into the result, give the Gibbs fundamental relation
d G = d G r e v = S r e v d T + V d P + µ k d N k 0
the quasi-static (reversible) change in Gibbs energy of the system between two equilibrium states according to the first law, valid for all systems. At constant pressure dP = 0,
d G = d G r e v = S r e v d T + µ k d N k 0
where d G = d G r e v is the free energy change in the battery, via the reversible (rev) path, maximum for energy transfer out of the battery (discharge), and minimum for energy transfer into the battery (charge). For reactions such as phase transitions and chemical formation/decomposition of substances, change in the Gibbs energy via reversible and irreversible paths, between states, can be used to calculate entropy change in the process.
Embedded in the compositional energy change term µ k d N k are the battery’s electrochemical kinetics, which include mass change due to chemical reactions d N r k , mass flow across system boundary d N e k , and ionic diffusion of mass within the system boundary d N d k , i.e., d N k = d N r k + d N e k + d N d k , where N k is the number of moles of material species k with N r k , N e k , and N d k the reactive, externally transferred and internally diffusive species, respectively. For a battery (closed system), d N e k = 0 gives d N k = d N r k + d N d k , and
µ k d N k = µ k r d N r k + µ h i g h µ l o w d N d k
accounts for chemical reactions r and diffusion d (note superscripts). Variation in chemical potential within the battery induces diffusion, hence µ h i g h and µ l o w are chemical potentials in the high- and low-potential regions, respectively [19,42], and d N d k is the amount of active species transported between both regions. Considering the source of the battery’s electrical energy is chemical and ionic, the reaction and diffusion terms in Equation (19) are replaced by the directly coupled electrical boundary work V d q , such that
µ k d N k = V d q
where V is the terminal voltage and d q = I d t is the charge transferred by battery current I over time dt. The derivations that convert Equation (19) to (20), which couple the electrochemical kinetics to the easily measured V and I, are given in Appendix A. The change in Gibbs free energy (Equation (18)) becomes, for the discharge/charge process,
d G = d G r e v = S r e v d T + V d q 0
Differentiating Gibbs energy Equation (16) with the pressure constant and substituting Equation (15) for dU into the result, give the irreversible form of the Gibbs relation
d G = δ G i r r = S d T + V d q T δ S 0
where d G = δ G i r r is the free electrochemical energy change in the battery, via the irreversible (irr) path, also maximum for energy transfer out of the battery (discharge) and minimum for energy transfer into the battery (charge). In Equation (22), by singling out the first two terms, the phenomenological Gibbs energy change can be defined as
δ G p h e n = S d T + V d q
due, only, to the changes in the physically observable and measurable intensive and extensive state variables ( d T and d q are easily measured). Note the overall change dG can be evaluated via a constant ideal change, d G r e v (Equation (21)), for a reversible process, or via a nonlinear, irreversible, and path-dependent real change, d G i r r (Equation (22))—both assess the total change in Gibbs free energy. With dG between beginning and end states known via d G r e v , Equations (21) and (22) can be combined to introduce the battery’s fundamental Gibbs-based entropy generation or production at constant pressure
δ S = S d T T + V d q T d G r e v T = δ G p h e n T d G r e v T 0
which satisfies the first and second laws. Equation (24) evaluates entropy generation δ S via the difference between Gibbs energy changes between two thermodynamic states, along reversible and phenomenological paths. During discharge (energy extraction), d T 0 , d q 0 and d G r e v 0 , rendering δ S 0 . For charge (energy addition), d T 0 ,   d q 0 and d G r e v 0 , reversing the signs of the middle terms in Equation (24) to preserve accordance with the second law: δ S 0 [42]. In batteries, d G r e v can be evaluated via
d G r e v = V O C d q r e v
where V O C is the battery’s standard potential (or open-circuit voltage) and d q r e v is reversible charge transfer.
Equation (24) defines entropy production as the difference between phenomenological Gibbs entropy d S p h e n = δ G p h e n T = S d T T + V d q T and reversible (ideal) Gibbs entropy d S r e v = d G r e v T , restated as
δ S = δ S p h e n d S r e v 0
where ( d S r e v δ S p h e n ) < 0 for discharge and 0 < ( d S r e v δ S p h e n ) for charge, making Equation (26), a statement of the second law, universally consistent with the thermodynamic process directional signs (IUPAC convention of positive energy into a system and negative work out).
Comparing Equations (12) to (13), and (21) to (22), and recalling the thermodynamic State Principle, shows that the changes in Gibbs entropy and energy between the two states are path-independent, whether the process path is reversible or irreversible, i.e.,
d S = d S r e v = δ S i r r = δ S p h e n δ S ,   d G = d G r e v = δ G i r r = δ G p h e n T δ S
Entropy and energy changes along an irreversible path must sum the instantaneously observed/measured phenomenological entropy and energy changes to the internally generated entropy and energy loss. Most battery analyses use reversible Gibbs energy change d G = d G r e v and entropy change d S = d S r e v [46,47]—only end-state measurements of system variables (before and after a discharge or charge step) are required, unlike d G = δ G i r r = δ G p h e n T δ S and d S = δ S i r r = δ S p h e n δ S , which require instantaneous accounting of all active processes. Since d G and dS during a process can be negative or positive, depending on heat flow TdSe or entropy flow dSe into or out of the battery, neither dG nor dS measures the permanent changes in the battery. On the other hand, entropy generation, Equation (24) or (26), evolves monotonically as per the second law. With δ S = 0 indicating an ideal (reversible) battery-process interaction, Equation (26) also indicates that a portion of any real battery’s energy is never available for external work, δ S > 0 . In Equation (26), the entropy generated by the battery’s internal irreversibilities alone, is in accordance with experience, similar to the Gouy-Stodola theorem used in availability (exergy) analysis [43,45,48,49]. Note since a discharged battery cannot ‘revert’ to a charged state without external recharging, all batteries are thermodynamically irreversible.

3.2.2. Relaxation/Settling and Self Discharge

During active discharging/charging, any heat generated and not instantaneously transferred out builds up. Upon load removal, the battery settles and the heat transfers out, as the Gibbs potential proceeds to a new equilibrium. During settling, the cell voltage relaxes and the battery transfers entropy to the surroundings. To apply Equation (24) to the settling process, the external Ohmic interaction Vdq, which replaced the compositional energy changes in Equations (17) and (18), is replaced by the internal ionic transport process energy (see Equation (19) and Equation (A8) in Appendix A). Phenomenological entropy production during settling becomes
δ S p h e n = S d T T k B R T 2 η M C x d t ln m L i + l o w m L i + h i g h
in terms of gas constant R, Li+ concentration MC, Boltzmann constant kB, distance x and dynamic friction coefficient η (inverse of electrical mobility). Here, S d T , positive for decreasing temperature dT ≤ 0, represents both voltage and thermal relaxation, and the last right-hand side term represents diffusion during settling, all of which proceed spontaneously and significantly slower than the active Ohmic processes [18,19]. Equation (24) also applies to self-discharge during storage, with the external Ohmic term Vdq representing spontaneous charge leakage [2].
With voltage relaxing in opposite directions, for discharge and charge, and often canceling during cycling, entropy production during settling proceeds at the rates of spontaneous cooling and diffusion of the charge species, which, relative to the preceding and subsequent active Ohmic energy and entropy interactions (charge/discharge step), is negligible. Relaxation equilibrium is approached asymptotically, over hours to days, rendering relaxation entropy production miniscule, compared to that during charging and discharging. Additionally, many batteries continue to supply power during charge, removing settling from the cycling schedule. Considering the above, settling effects and its entropy is to be neglected.
Equations (23) and (24), over time increments dt, give rates G ˙ = S T ˙ + V I and S ˙ = S T ˙ T + V I T G ˙ r e v T . Integrating from t0 to t gives the total Gibbs energy and entropy generation
Δ G = t 0 t S T ˙ d t + t 0 t V I d t
S = t 0 t S T ˙ T d t + t 0 t V I T d t t 0 t G ˙ r e v T d t

3.3. Entropy Content S and Internal Free Energy Dissipation −SdT

In Equations (17), (21), and (22), the term −SdT, which represents free energy dissipated internally and not instantaneously transferred out during active work interaction, includes Ohmic dissipation and can include heat from an external source. The temperature change dT is driven by the entropy content S of the system. Resolving entropy into components, S = S µ , T depends on chemical potentials µ and temperature T, for a chemically reactive system. The Gibbs–Duhem equation [19,39,40,41,42,43,44]
S d T V d P + N k d µ k = 0
at constant pressure yields
S d T = N k d µ k
where N k d µ k = N d µ for one active species (e.g., Li+). The molar Gibbs energy at constant temperature and pressure gives the following, via Equation (17),
µ = G N T , P
Substituting, from Equation (25), G ˙ r e v = V O C I r e v and, from Equation (32), S T ˙ = N µ ˙ into Equations (29) and (30) yields
Δ G = t 0 t N µ ˙ d t + t 0 t V I d t
S = t 0 t N µ ˙ T d t + t 0 t V I T d t t 0 t V O C I r e v T d t
Integrating Appendix A Equation (A14), Faraday’s first law, yields
N = C n F
where C is the charge in the battery at time t given as
C t = C 0 ± C t Δ t 0
C 0 is the battery’s initial charge content at t = t 0 , and
C t Δ t = t 0 t I t d t
is the total charge transferred at time t. Combining Equation (20) with Equation (32) gives, in rate form,
S T ˙ = N µ ˙ = C V ˙
Substituting Equation (39) into Equations (34) and (35) yields
Δ G = t 0 t C V ˙ d t + t 0 t V I d t
S = t 0 t C V ˙ T d t + t 0 t V I T d t t 0 t V O C I r e v T d t
Equations (40) and (41) give Gibbs energy change and entropy production as functions of instantaneous voltage V, current I, and temperature T.
The Gibbs–Duhem constant-pressure formulation in rate form, Equation (39), shows interdependence of battery voltage and temperature—temperature rise drops voltage, and voltage drops as charge exits raises temperature (typical for active discharge), see Equations (17), (21) and (22), wherein −SdT reduces the Gibbs energy for positive S and dT. Being dependent on voltage change, charge content and temperature, the batter’s internal accumulation of free energy dissipation −SdT is more appropriately named Electro-Chemico-Thermal (ECT) energy change, which for a reactive system depends on the number of moles, chemical potential, heat capacity, and changes in temperatures of reactive components [19,39,40,41,42,43,44]. Variables V ˙ and I, negative during discharge and positive during charge, establish similar directional signs for the Gibbs energy change and entropy generation terms. From Equation (41), the Ohmic entropy S Ω = V I T d t accounts for the loss in output/input work, the reversible Gibbs entropy S r e v = V O C I r e v T d t measures the minimum possible loss, while the ECT entropy S V T = C V ˙ T d t accounts for loss due to rise in internal entropy content and temperature, as battery potential and charge content decrease. Similarly, from Equation (40), Ohmic work G Ω = t 0 t I V d t and ECT energy G V T = t 0 t C V ˙ d t .
During relaxation/settling, entropy generation in Equation (41) is obtained by applying the Gibbs–Duhem formulation (Equation (31) at constant P) to Equation (28). Without any active charge transfer during settling, Equation (41) indicates that ECT entropy, the first right-hand side term representing entropy accumulation from both voltage and thermal relaxation via Equation (39), is the most significant entropy generation component. Similarly, entropy generation during battery storage, including effects of self-discharge, can be evaluated via Equation (41).

3.4. Degradation-Entropy Generation (DEG) and Capacity Fade in Batteries

In primary cells, degradation in the form of capacity loss is the difference between initial capacity at initial time t 0 and capacity at a latter time t , as the battery discharges monotonically. In secondary cells, the charge step reverses this ‘loss’, making the prior definition unsuitable for describing permanent loss of charge capacity. However, over time, secondary cells degrade and lose their ability to hold charge, resulting in capacity fade. This capacity fade is estimated as
Δ C C C = C 1 | C N |
the difference between the first cycle’s charge content or capacity C 1 and the Nth cycle’s charge content C N Δ t (capacity C N Δ t can be replaced by the total charge transferred C t , measured after a discharge step), sometimes expressed as a percentage of C 1 [15,50]. Using Coulomb/Charge Counting (subscript CC), Equation (42) requires a consistent cycling schedule and low constant discharge rate for cycles 1 and N, hence, is not applicable to irregular and inconsistent cycling schedule, with deep discharges and undercharges. Similar to Section 2.2, the above Gibbs-based formulations are combined with the DEG theorem, as follows:
  • Let available battery capacity or charge content C be a DEG transformation measure and capacity fade (lost discharge/charge capacity) Δ C be the observed/measured degradation, the DEG Equation (9), with Δ C replacing Δ w becomes,
    Δ C = i B i S i
  • From Equation (41), entropy generation S’ = S’ V t , I t , T t , suggesting via Equation (43) that C = C V t , I t , T t . Substituting the entropy generation terms of Equation (41) into Equation (43) gives
    Δ C = B V T t 0 t C V ˙ T d t + B Ω t 0 t V I T d t B G t 0 t V O C I r e v T d t
  • Via Equation (8), with C replacing w, DEG coefficients
    B V T = C S V T ,   B Ω = C S Ω , B G = C S r e v
    can be evaluated via measurements, as slopes of charge C versus ECT entropy S V T , Ohmic entropy S Ω and reversible Gibbs entropy S r e v . As demonstrated in Section 4, the DEG capacity fade, Equation (44), is applicable to all forms of cycling conditions.

4. Experiments

Data sets arose from degradation of four same-model Lithium-ion batteries, under inconsistent and abusive cycling. Parameters were uncontrolled—cycles discharged to different depths and recharged to different heights, with some discharge steps extending significantly below the manufacturer-specified full-discharge depth of 2.7 V, to induce rapid degradation and severely nonlinear behavior, as this would test the robustness of the DEG degradation model. Cycle 1 discharged to 0.8 V and cycles 27 and 29 discharged down to 0.6 V. A constant current of 3 A was used for the charge steps. Figure 1 shows all four batteries during cycling. Measured cyclic data sets are in Appendix B.

4.1. Apparatus

Each setup (Figure 1) had a single-cell 3.7 V 11.5 Ah lithium-ion polymer battery with graphite anode and aluminum current-collectors, model PL-9059156, manufactured by Batteryspace; a Hitech X48 Multi-charger, powered by a DC power supply for the charge cycle; a set of Dale RH-50 1 Ω , 50 W resistors for a standardized uniform resistive load; lead wires with known gages; two duplex-insulated OMEGA TT-K-24-25 Type K thermocouples (one to measure battery temperature and the other to measure ambient temperature); a current, voltage, and resistance meter (or multimeter) (to measure pre- and post-cycling voltages and resistances), a National Instrument controller/data acquisition system NI CompactRIO with an analog input module NI 9215 (to monitor battery and resistor voltages), a digital output module NI 9401 to automate the cycling process; and a thermocouple module NI 9211 (with internal cold-junction compensation and auto-zero correction for measurement accuracy) to monitor ambient and battery temperatures.

4.2. Setup and Procedure

Figure 2 depicts the battery circuit. Monitored during cycling were battery terminal voltage V B dependent on the switch state (set to node “c” during charge, set to node “d” during discharge and opened to measure battery open-circuit voltage V O C ), voltage drop across load resistor V R , battery temperature T, and ambient temperature T near the battery. To obtain sufficient data for statistical significance and establish repeatability, four Li-ion batteries were cycled at two different unsteady discharge rates.
Battery tests were conducted in a well-ventilated area. All tests followed the same procedure.

4.2.1. Setup and Initial Measurements

For one battery, two 1 Ω resistors were used in parallel (actual resistance including wiring R l o a d = 0.6 Ω, measured via a multimeter); for the remaining three batteries, three 1 Ω resistors were used in parallel ( R l o a d = 0.45 Ω); the sensing terminal of one thermocouple was attached to the battery surface and the ambient thermocouple was placed in air, near the battery, with the sensing terminal not in contact with any surface; the discharge–charge cycling circuit was set up, with the resistor loads for cell discharging and current measurements, and a battery charger for charging (Figure 2); the initial battery voltage and temperature, and resistance of the load resistor network were measured.

4.2.2. Cycling

With the battery voltage V B 3.7 V and recording transient data at 0.1 Hz, the resistor network was connected to the battery to begin the discharge step; the resistors were disconnected from the battery after V B fell to 2.7 V (full discharge) or below (over-discharge); the battery was allowed to settle for 20–25 min (for voltage and thermal relaxation); then charged at 3 A, varying the charge depths for subsequent cycles; and allowed to settle for about 30 min. This was repeated until the battery degraded, observed in one of two ways:
  • the battery’s capacity fell to less than two-thirds the initial capacity or
  • the battery began to inflate in geometric volume (close monitoring of Li-ion batteries was required during cycling).

5. Results, Analysis, and Discussion

Monitored parameters changed with time at unsteady rates. Tables listed in this paper contain data for discharge (left side of table) and charge (right side of table). Signs indicated either a parametric decrease/increase or the direction of the process; for example, ECT energy/entropy and Ohmic work/entropy are negative for discharge, and positive for charge. Path-dependent integrals were evaluated using the trapezoidal rule on data over time increment Δ t . With data sampled at 0.1 Hz, Δ t = 10 s. Plots show the direction of accumulation and the rates, starting with the discharge step and ending with the charge step.
Using equations for estimating charge content/transfer, Gibbs energy, and entropy (Section 3), data from the lithium-ion battery cycling experiments are presented in Table 1 and Table 2, for a randomly selected battery #2. Similar trends were observed in the other same-model batteries tested. Table 1 presents the Gibbs energy-entropy data and Table 2 presents the DEG data. In both tables, column 1 contains variable N numbers the discharge–charge cycles. To add to cyclic inconsistencies, a “Half cycle” (Table 1 and Table 2 and Table A1) is a charge-only cycle, which occurs when two charge steps follow a discharge step. For example, in cycle 2, the battery was fully discharged and re-charged to about 50% nominal capacity, rested for a few hours, re-charged again in cycle 3, before discharging in cycle 4; hence, cycle 3, having no discharge step, is a half cycle. In accordance with battery industry, Ah, Wh and Wh/K are used for charge, energy, and entropy, respectively. (1 Ah = 3600 As = 3600 C; 1 Wh = 3600 Ws = 3600 J; 1 Wh/K = 3600 J/K), giving B coefficients units of Ah K/Wh, to differentiate them from the reciprocal of the voltage-temperature coefficient which has units of V/K (i.e., 1 Ah K/Wh = 1 K/V).
Figure 3 plots measured discharge and charge current I, battery voltage V, and temperature T versus time, during cycle 6. With a constant resistive load R l o a d , V and I trend similarly during discharge. In Figure 3, battery temperature rose during discharge but remained relatively steady during charge, trending with ambient temperature.

5.1. Gibbs Energy and Entropy Components

Columns 2 and 7 of Table 1 list the discharge and charge Ohmic work G Ω , and columns 3 and 8 show the ECT energy G V T . Figure 4 plots Ohmic S Ω (columns 4 and 9 of Table 1), ECT S V T (columns 5 and 10) and reversible Gibbs S r e v (columns 6, 11) entropies, versus charge transferred. Continuous curves represent discharge and dashed curves represent charge profiles, both processes proceeding in opposite directions. Measured charge transfer (Figure 4) appears linear vs. Ohmic entropy S Ω (blue plots) and curvilinear vs. ECT entropy S V T (red plots). Voltage relaxation after discharge typically reduces charge S V T . In other words, when the external load is disconnected from the battery, the battery’s potential V immediately rises, indicating voltage elasticity and causing the voltage change during active charge to be less than that during discharge. This contributes to higher Coulombic efficiency of the charge step in battery cycling, and verifies that charging is more thermodynamically reversible than discharging, and hence generates less entropy. Note that a reversible (ideal) step is one that proceeds at a constant current and voltage. Figure 4 shows Ohmic entropy to be several times that of the ECT entropy, seemingly justifying the latter’s absence in many battery analyses. However, the significance of the ECT entropy is underscored by the need to avoid overdischarging and overcharging, and keep batteries cool during operation, for better and longer performance.
Figure 5 plots rates of entropy generation components versus time. Ohmic entropy rate starts relatively high and decreases during discharge; vice versa during charge. ECT entropy generation rate fluctuates significantly with a slight overall change, during both discharge and charge. These fluctuations have higher amplitudes during discharge than charge. Rate of reversible Gibbs entropy, steady for both discharge and charge, is higher for discharge via higher discharge rate.
Figure 6 plots the primary components of entropy generation—phenomenological entropy S p h e n (purple plots) and reversible entropy S r e v (green plots)—during discharge and charge steps of cycle 6. The region between S r e v and S p h e n curves represents the entropy generated S . During discharge and charge, S r e v < S p h e n as anticipated by Equation (26), in accordance with the second law. Entropy generation is significantly higher during discharge with higher values of current and deep discharge.
Similar results were observed for all cycles and all four batteries.
The final ‘Summary/Total’ row of Table 1 shows that the total Ohmic entropy S Ω is about 10 times the total ECT entropy S V T for all the discharge steps combined, whereas all the charge steps collectively accumulated ratio S Ω : S V T at the much more thermodynamically reversible 46:1. Discharge S V T is 4.4 times the charge S V T (via a higher entropy generation rate from discharge rate >> charge rate and deep discharging), while charge S Ω is higher than the discharge S Ω (via the overall charge duration >> discharge duration). Total reversible entropy S r e v is less than the total phenomenological entropy for both discharge and charge, in accordance with thermodynamic formulations.

5.2. DEG—Capacity Versus Entropy

By associating data for entropies and battery charge at various time instants, Figure 7 plots charge C (measured via integration of battery current) versus accumulated entropy components S Ω (Ohmic) and S V T (ECT). Figure 7a shows the three-dimensional plots for discharge, curve ‘D’ on the red plane, and charge, curve ‘C’ on the green plane. End views in Figure 7b show a near perfect coincidence of all data points onto these planes. A goodness of fit R 2 1 was obtained for all measured data, for all cycles of all four batteries tested. The near perfect coincidences onto 2D planes suggest a linear dependence of battery charge on both Ohmic and ECT entropies at every instant of the discharge/charge process, consistent with the predicted linearity of Equation (43). The dimensions of the horizontal axes and their maximum values are determined by the accumulation of the entropy generation components at the end of the discharge/charge step.
In Figure 7, the measured data points for both discharge and charge form distinct paths––deemed Degradation-Entropy Generation (DEG) trajectories which lie on (coincide with) planar surface(s), deemed DEG plane(s)/surface(s). The 3D space of the DEG surface(s), the battery’s DEG domain, appears to characterize the allowable regime in which the battery can operate.
Degradation Coefficients: The red and green DEG planes in Figure 7 have orientations in the 3D DEG domain that can be specified in terms of partial derivatives of charge with respect to Ohmic and ECT entropies respectively. Via Equation (45), these orientations (or partial derivatives) define the Ohmic and ECT degradation coefficients B Ω and B V T . Data from a reference cycle (cycle 1 in this study) was used to obtain the battery’s characteristic DEG coefficients and predict the rest of the battery’s operational life. From cycle 1 data, B Ω = 76.6 Ah K/Wh and B V T = 113 Ah K/Wh for discharge, and for charge, B Ω = 75.5 Ah K/Wh and B V T = 28.3 Ah K/Wh.

5.2.1. Phenomenological Charge, Measured Charge, and Reversible Charge

The Carnot limitation—a corollary constraint of the second law which governs availability of a system’s energy for work—was expressed by Burghardt [45] for a heat source as
E n e r g y   a d d e d   =   A v a i l a b l e   e n e r g y   +   U n a v a i l a b l e   e n e r g y
Defining phenomenological charge transfer C p h e n , the battery’s actual physically observable path, which can be correlated directly with phenomenological entropy via the first two right-hand side terms in Equation (41), as,
C p h e n = C t + Δ C p h e n = B V T t 0 t C V ˙ T d t + B Ω t 0 t V I T d t
where the phenomenological capacity fade Δ C p h e n , a portion of C p h e n not available for external/boundary work during cycling, due to instantaneous dissipation from battery heating and loss of charge and potential, is the difference between the phenomenological charge transfer C p h e n and measured charge transfer C t = t 0 t I t d t . Similarly, defining reversible charge transfer C r e v that establishes the battery’s ideal path and correlates with reversible entropy, the last term of Equation (41),
C r e v = C t + Δ C r e v = B G t 0 t V O C I r e v T d t
where reversible charge fade Δ C r e v is the difference between reversible (ideal) charge transfer
C r e v = t 0 t I r e v d t = I r e v Δ t
and measured charge transfer C t (or accumulated discharge/charge). Here C r e v is maximum during discharge (maximum discharge capacity) and minimum during charge (minimum charge capacity); Δ C r e v represents the portion of C r e v unavailable due to previous permanent degradation and instantaneous dissipation; and I r e v is the constant reversible current, maximum during discharge and minimum during charge.
Following Equation (26), the difference between Equations (47) and (48) derives Equation (44), the DEG capacity fade model, which quantifies actual capacity fade from degradation
Δ C D E G Δ t = C p h e n C r e v = B V T t 0 t C V ˙ T d t + B Ω t 0 t V I T d t B G t 0 t V O C I r e v T d t
as the difference between phenomenological and reversible charge transfer components. Inspection of Equations (48) and (49) indicates that B G = d C r e v d S r e v = T V O C . Substituting I r e v Δ t (from Equation (49) with I r e v constant) for the last term in Equation (50) gives
Δ C D E G Δ t = B V T t 0 t C V ˙ T d t + B Ω t 0 t V I T d t I r e v Δ t
While C t and C r e v are determined from currents I and I r e v , and C p h e n is not directly measurable. This makes Equation (51), the simplified form of the DEG model in Equation (44)—which requires only measurements of V, I and T—convenient for instantaneous evaluation of operational capacity fade, irrespective of discharge rate, depth of discharge, and cycling inconsistencies.
During discharge, C r e v in Equations (48) and (49) represents the overall maximum charge available in the battery, only obtained from new batteries (at t = t 0 , C r e v = C p h e n ) or if the battery operates as a perfect energy source or sink (wherein no output or input power converts to heat or degrades the battery). As the battery ages, C p h e n diminishes (i.e., at t > t 0 , C p h e n < C r e v ). As degradation continues, the amount of energy required to restore the battery’s original charged state increases. Hence, during charge, C r e v is the overall minimum charge required to restore the battery to its initial state, realizable in new batteries (i.e., at t = t 0 , C r e v = C p h e n ), while C p h e n is the actual increasing charge from the external charger, required to restore the battery to its original state (i.e., at t > t 0 , C p h e n > C r e v ). This implies Δ C > 0 consistently represents the capacity fade during both discharge and charge in accordance with positive entropy generation S’ > 0, Equation (26).
Table 2 presents the components of the DEG capacity fade Δ C D E G and Coulomb-Counted capacity fade Δ C C C for the discharge and charge steps of all 32 cycles monitored.

5.2.2. Evaluating Capacity Fade—Battery Cycle Life Model

  • Cyclic values of DEG capacity fade Δ C D E G and Coulomb-Counted capacity fade Δ C C C presented in Table 2 for cycles 1 to 32, were obtained as follows:
  • Evaluate phenomenological charge C p h e n in Table 2, columns 2 (discharge) and 7 (charge), from Equation (47), by combining reference DEG coefficients with each cycle’s phenomenological entropy components S Ω and S V T given in Table 1. For example, the discharge step of cycle 6, row 6 of Table 1, has S Ω = −0.08 Wh/K, and S VT = –0.005 Wh/K. When combined with B Ω = 76.6 Ah K/Wh and B V T = 113 Ah K/Wh from cycle 1, Equation (47) gives C phen = (76.6 × −0.08) + (113 × −0.005) = −6.7 Ah.
  • Evaluate reversible charge C r e v (Table 2, columns 3 and 8) from Equation (49). For example, the cycle 1 starting current I 1 t 1 = −5.2 A gives the cycle 6 discharge I r e v = −5.2 A which, with the cycle 6 discharge duration Δ t = 1.53 h (Table A1 in Appendix B), gives C rev = −5.2 × 1.53= −8.0 Ah.
  • Evaluate DEG capacity fade (Table 2, columns 4 and 9) from Equation (51), Δ C D E G = C p h e n C r e v . For cycle 6 discharge, Δ C DEG = −6.7 − (−8.0) = 1.3 Ah. Similarly, for cycle 6 charge, Δ C DEG = 0.3 Ah.
  • Evaluate cyclic Coulomb-Counted charge transfer (Table 2, columns 5 and 10), C CC = t 0 t I t dt where I(t) is instantaneously measured unsteady current. For cycle 6 discharge and charge, row 6 of Table 2, C CC = −7.2 Ah and 4.4 Ah respectively.
  • Evaluate Coulomb-Counted capacity fade (Table 2, column 6) from Equation (42), Δ C C C = C 1 | C C C | , where C 1 = −6.1 Ah is Coulomb-Counted charge transfer during cycle 1 discharge. For cycle 6, Δ C CC = |−6.1| − |−7.2| = −1.1 Ah.
From the above procedure, substituting cycle 1 discharge’s pair of DEG coefficients (given in Section 4.2) and reversible current—obtained as cycle 1’s starting discharge current I r e v = I 1 t 1 —into Equation (51), the discharge (subscript d) capacity fade model for the Li-ion battery used in this study becomes,
Δ C D E G d = 113 S V T + 76.6 S Ω I 1 t 1 Δ t
Similarly, for the charge steps, with B Ω = 75.5 Ah K/Wh, B V T = 28.3 Ah K/Wh and I r e v = I 1 t n , obtained at the end of cycle 1′s charge step t n , the battery’s charge c capacity fade model is
Δ C D E G c = 28.3 S V T + 75.5 S Ω I 1 t n Δ t
The battery cycle life models, Equations (52) and (53), are in the form of first-order partial differential equations with constant coefficients. The slight difference in discharge and charge B Ω values is primarily attributed to the difference in discharge and charge currents, while the significant difference in discharge and charge B V T is due to the different instantaneous responses to the different operating conditions. Figure 8, plots of all 32 cyclic phenomenological charges vs. cyclic Ohmic and ECT entropies, shows each discharge and charge step’s characteristic operational conditions (depth of discharge, etc.) and response to those conditions, via the scatter in the cyclic C p h e n data. However, as predicted by Equations (52) and (53), and more easily observed in Figure 8b, all cyclic responses by the battery to the given operational conditions lay on the same DEG planes, giving a consistent prediction of the battery’s cycle life.
Cycle 1 characterized the battery via the B coefficients and Equations (52) and (53), or Equation (44). Data from subsequent cycles were used to monitor the battery’s degradation with cycling. Table 1 and Table 2 show that steps with high electro-chemico-thermal (ECT) entropy (relative to Ohmic entropy) had high DEG capacity fade Δ C D E G (discharge steps 27 and 29, both of which had the lowest discharge cut-off voltage of 0.6 V, indicating the deepest discharge—see Table A1 in Appendix B). With each charge step proceeding at the steady rate of 3 A, the charge entropy generation was determined primarily by the duration of charge, hence steps with more accumulated charge tended to have higher entropy generation and hence, higher charge capacity fade (charge steps 1, 7, 25, 31).
Due to the severely inconsistent cycling of varying discharge and charge depths, Coulomb-Counted capacity fade Δ C C C (Equation (42)) which depended only on current, was not applicable to the measured data sets. Discharge values are presented in Table 2, column 6, only for comparison to the DEG capacity fade Δ C D E G . The effects of the irregular discharge durations are easily observed in Δ C C C values in Table 2, with most cycles falsely indicating improvement (increased charge capacity) over cycle 1, even with the insufficient recharge in the preceding charge steps and low discharge depths—down to < 1 V, in several cases. The DEG capacity fade is a function of each cycle’s entropy generation, which is in turn a function of each cycle’s phenomenological state variables (V, I and T). Table 2’s ‘Summary/Total’ row shows the total DEG discharge capacity fade of 39.3 Ah, relative to an ideal/reversible total of 245 Ah discharge obtained from the battery (16.0%), and charge capacity fade of 8.2 Ah relative to a total of 176.8 Ah charge supplied to the battery (4.6%), suggesting that most of the battery’s loss of usable charge-holding capacity occurred during discharge (particularly deep discharge), in accordance with experience. Based on the battery’s nominal capacity rating of 11.5 Ah, the battery’s nominal DEG capacity fade was estimated via
Δ C D E G Δ t = C p h e n C r e v C r e v C n o r m
as (0.16 × 11.5=) 1.83 Ah discharge fade and (0.046 × 11.5=) 0.53 Ah charge fade.
Figure 9 plots capacity fade components—measured charge C t (blue plot), phenomenological charge C p h e n (purple plot), and reversible charge C r e v (green plot)—during discharge and charge steps of cycle 6, with capacity fade Δ C as the region between C r e v and C p h e n curves. During both discharge and charge, C r e v < C p h e n ; note that using the magnitude only, during discharge, C p h e n < C r e v . These consistent trends were anticipated by formulations in accordance with the second law, see previous sections. As with entropy generation in Figure 6, discharge capacity fade was significantly higher than the charge capacity fade, due to a faster and deeper discharge than charge. With relatively low ECT entropy, C p h e n C t (Figure 9).

6. Discussion

Prigogine [42] introduced a universally non-positive (for macroscopic systems undergoing spontaneous processes), interaction-specific “local potential”, which is analogous to the ECT entropy introduced in this study. This article derived and experimentally verified a universally consistent system-based, time-dependent entropy generation. It was shown that,
  • phenomenological entropy generation S p h e n is the sum of Ohmic entropy S Ω and electro-chemico-thermal ECT entropy SVT;
  • entropy generation is the difference between phenomenological S p h e n and reversible S r e v Gibbs entropies, at every instant;
  • entropy generation is always non-negative, in accordance with the second law, whereas components S p h e n and S r e v are directional—positive during charge and negative during discharge. This implies S p h e n S r e v during charge and S p h e n S r e v during discharge, in accordance with experience and thermodynamic laws. This article demonstrated the significance of the previously neglected reversible S r e v and ECT S V T entropies in evaluating entropy generation in batteries.

6.1. Features of the DEG Theorem and Coefficients

The DEG methods can accurately describe the battery’s charge levels (or state of charge SOC) within a discharge–charge cycle and the capacity fade over multiple discharge–charge cycles, using entropy generation components. DEG coefficients relate operational capacity fade to entropy generation in a rechargeable battery at any point in the battery’s life, by quantifying the dissipative processes’ individual contributions to the battery’s degradation. Ohmic coefficient B Ω > 0 via Equation (44), where phenomenological entropy components and charge transfer are negative during discharge and positive during charge, while ECT coefficient B V T has a varying sign. This is implied in Equation (47), rewritten as
B V T = 1 S V T C p h e n B Ω S Ω
where both C p h e n and B Ω S Ω are negative during discharge and positive during charge.

6.1.1. DEG Trajectories, Surfaces, and Domains

The DEG theorem converts degradation failure modeling and concomitant design into a multi-dimensional geometry problem, as Figure 7 and Figure 8 suggest. The volume spanned by entropy trajectories defines the operating and aging region, as well as the consistent parameters for identifying desired characteristics from batteries of all configurations. Phenomenological entropy components serve as axes in the multi-dimensional space of C versus S i , where i numbers the active processes. DEG trajectories characterize the discharge–charge cycle conditions; inclination and location of DEG surfaces appear to characterize a battery’s discharge/charge rates; and the DEG domain (here Charge versus Ohmic and ECT Entropies) seems to characterize the battery life, for all cycles and all rates. A battery having a domain with large accumulated charge dimension and small ECT entropy dimension (relative to Ohmic entropy dimension) delivers power more efficiently.

7. Summary and Conclusions

This study applied irreversible thermodynamics to instantaneous characterization of lithium-ion battery degradation. The DEG theorem was successfully applied to nonlinear lithium-ion battery cycling, using non-intrusive measurements of temperature, voltage, and current. Detailed thermodynamic breakdown of the active electrochemical processes in lithium-ion batteries during cycling was presented, and the kinetics (see Appendix A) were represented with universally consistent and easily measurable battery operational parameters to determine the components of energy change and entropy generation. Via entropy generation evaluation, battery capacity fade model was formulated from the DEG theorem and experimentally verified, for four same-model lithium-ion batteries. The electro-chemico-thermal ECT entropy was introduced and its contribution to capacity fade was demonstrated. In the formulations of this article, operational capacity fade depended on the usual phenomenological variables used by the battery industry, but through the entropy generated by the underlying dissipative processes. A thermodynamic potential—Gibbs free energy—replaced the steady state assumptions of previous DEG applications, and employed the instantaneous applicability of the first and second laws of thermodynamics. The DEG methodology could directly compare the same-model and different-model lithium-ion batteries, as well as technologies, designs, and materials used in manufacturing the battery.

Author Contributions

Conceptualization, J.A.O. and M.D.B.; Data curation, J.A.O.; Formal analysis, J.A.O.; Funding acquisition, M.D.B.; Investigation, J.A.O.; Methodology, J.A.O. and Michael Bryant; Resources, J.A.O. and M.D.B.; Software, J.A.O.; Validation, J.A.O.; Visualization, J.A.O.; Writing—original draft, J.A.O.; Writing—review & editing, J.A.O. and M.D.B..

Acknowledgments

The authors thank the University of Texas at Austin Mechanical Engineering Department’s Fall 2015 semester Mechatronics I students for their assistance in the battery measurements.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

NomenclatureNameUnit
A chemical affinityJ/mol
BDEG coefficientAh K/Wh
C charge, charge transfer or capacityAh
Δ C charge or capacity fadeAh
FFaraday’s constantC/mol
GGibbs energyWh
Idischarge/charge current or rateA
kBBoltzmann constantJ/K
mmasskg
n’number of charge species
Ncycle numberKg/mol
N, Nknumber of moles of substancemol
pdissipative process energyJ
PpressurePa
qchargeAh
QheatJ
Rgas constantJ/mol·K
Sentropy or entropy contentWh/K
S’entropy generation or productionWh/K
ttimesec
TtemperaturedegC or K
Uinternal energyJ
VvoltageV
V volumem3
wdegradation measure
WworkJ
Symbols
μchemical potential
ζphenomenological variable
Subscripts & acronyms
ΩOhmic
0initial
ccharge
ddischarge
ECT, VTElectro-Chemico-Thermal
ttime
revreversible
irrirreversible
phenphenomenological
CCCoulomb-Counted
DEGDegradation-Entropy Generation

Appendix A

Appendix A.1. Electrochemical Kinetics

Reversible reduction and oxidation reactions at the cathode and anode, and transport of ions through the electrolyte underlie lithium-ion battery operation. In Equation (18), these kinetics are embedded in the compositional change term µ k d N k . This section formulates equations governing electrochemical kinetics of charge formation (ionization reaction) and charge transport [19,38,42].

Appendix A.1.1. Charge Intercalation (Absorption) and Deintercalation (Release)

For the reversible chemical reactions at the cathode and anode expressed respectively in Equations (1) and (2), changes in mole numbers of active species can be related stoichiometrically, via
d N L i 1 x M O 2 1 = d N L i + x = d N e x = d N L i M O 2 1 = d ξ c a d N L i x C 6 1 = d N C 6 1 = d N L i + x = d N e x = d ξ a n
where d ξ c a = d ξ a n = d ξ is the reaction extent. Subscripts ca and an represent the cathode and anode, respectively. Combining the above renders an equation involving de Donder’s reaction affinity A, which drives the reaction [1], that is
µ k d N k = A d ξ .

Appendix A.1.2. Reaction Rates

Combining forward reaction rate K f = l f a L i 1 x M O 2 a L i x C 6 applicable to battery charging, reverse reaction rate K r = l r a L i M O 2 a C 6 applicable to battery discharging, affinity A = R T ln K f K r [19], and reaction extent d ξ = ( K f K r ) , gives
µ k d N k = R T ( K f K r ) ln K f K r .
Here l f , l r are the reaction rate constants, often evaluated from the Arrhenius equation l = l 0 e E A / R T , a represents the activity of a species given by the subscript, and R is the gas constant.

Appendix A.1.3. Charge Transport

In a typical reaction, consumption of an active species slows the reaction, hence the chemical potential depends on the amount of the species. Similarly, in a continuous system containing two regions of an active species, diffusion will occur if the amount of the species in one region (and consequently the region’s chemical potential) is lower than in the other region. Diffusion tends to equalize the chemical potentials of both regions, hence the thermodynamic force that drives this flow, the diffusion affinity µ h i g h µ l o w , is the difference in species chemical potentials in regions inside the system [19], e.g., battery. In terms of the reaction extent, diffusion work is given as µ h i g h µ l o w d ξ . As done for chemical affinity (using the reaction rates) above,
µ h i g h µ l o w = R T ln m L i + l o w m L i + h i g h
expressed in terms of molalities m L i + h i g h , m L i + l o w (mol/kg) of Li+ in both regions. Activity a k l o w could also be used for convenience.
Fick’s laws of diffusion and ion conservation [1,19] govern the diffusion work in the presence of concentration gradients. The diffusion current N ˙ d k ( x k ,t) (mol/cm2·s) in terms of Li+ concentration Mc ( x k ,t) (mol/cm3) and the diffusion coefficient D (cm2/s) can be written as
N ˙ d k = D M C x k
Using the position-dependent form of the chemical potential µ k x = µ k + R T l n x k , a consistent thermodynamic formulation for diffusion current is given as [19]
N ˙ d k = R L M C M C x k
which by comparison to Equation (A4) yields L = D M C R . Values of D are widely available [51]. x k measures position.
In addition to the concentration gradients, the presence of an electric field facilitates ion transport between the two regions. Based on a position- and field potential-dependent chemical potential µ k x τ = µ k + R T l n x k + τ k to consider the resulting ionic drift ( τ k measures field potential), Einstein related ionic mobility Γ to the diffusion coefficient D, which in terms of dynamic friction coefficient η (inverse of electrical mobility) and the Boltzmann constant kB gives the Stokes-Einstein equation for electrical mobility
D = k B T η .
Combining diffusion affinity Equation (A3), Fick’s law of diffusion Equation (A4), and Stokes-Einstein relation Equation (A6), gives the diffusion-drift work
µ h i g h µ l o w d N d k = k B R T 2 η M C x d t ln m L i + l o w m L i + h i g h .
Adding chemical reaction (Equation (A2)) and ion diffusion-drift (Equation (A7)), gives compositional energy change in the battery
µ k d N k = R T [ ( K f K r ) ln K f K r k B T η M C x d t ln m L i + l o w m L i + h i g h ]
which accounts for all primary chemical (reaction) and electrical (transport) kinetics.

Appendix A.2. Coupling Reaction and Transport Kinetics with the Electrochemical Potential

Equation (A9) is not convenient for most experimental measurements and performance analyses of a battery in operation. To account for the chemical reactions at the electrodes and the resulting electrical forces in electrolytes, electrochemical potential µ ˜ k = µ k + n F ϕ (a combination of the chemical potential µ k of the active species k and an electrical field potential ϕ) is typically used to characterize the battery’s useful energy, where n is number of species, e.g., electrons, involved in charge transfer (x for lithium-ion batteries, Equations (1) and (2)), F = 96,485 C/mol is Faraday’s constant. Defining electrochemical reaction affinity [19,42] for each electrode half-reaction,
A ˜ c a = µ L i 1 x M O 2 c a + x µ L i + c a + x µ e c a µ L i M O 2 c a n F ϕ c a A ˜ a n = µ L i x C 6 a n µ C 6 a n x µ L i + a n x µ e a n + n F ϕ a n .
Summing both half-reactions and rearranging to obtain the overall electrochemical affinity,
A ˜ = µ L i 1 x M O 2 c a + µ L i x C 6 a n µ L i M O 2 c a µ C 6 a n + x µ L i + c a µ L i + a n + x µ e c a µ e a n n F ϕ c a ϕ a n .
For the directly-coupled electrode reactions, the charge potentials cancel out, i.e., µ L i + c a = µ L i + a n and µ e c a = µ e a n . Rearranging with A ˜ = 0 (equilibrium condition) gives the electromotive force or potential difference between both electrodes or the battery’s terminal voltage
1 n F µ L i 1 x M O 2 c a + µ L i x C 6 a n µ L i M O 2 c a µ C 6 a n = ϕ c a ϕ a n = V
Comparison of the terms in the bracket on the left-hand side of Equation (A12), to the overall reaction in Equation (3), indicates that Equation (A12) could have been derived by directly considering the overall reaction. If µ k = µ L i 1 x M O 2 c a + µ L i x C 6 a n µ L i M O 2 c a µ C 6 a n Equation (A12) becomes
1 n F µ k = V .
Equation (A1) and Faraday’s law of electrolysis give
d N k = d q n F
which combines with Equation (A13) to give,
µ k d N k = V d q ,
allowing the convenient use of the easily measured terminal voltage V and charge transfer dq or current I = d q / d t .

Appendix B

Experimental Results

A battery cycle consisted of discharge followed by charge. The discharge currents were measured as I = V R / R l o a d (see Figure 2). Constant-charge current of 3 A was used for the charge steps. Using equations for estimating battery capacity and Ohmic work (Joule dissipation in the load resistors), Table A1 presents data from the lithium-ion battery cycling experiments—discharge (left side) and charge (right side). N numbers the discharge–charge cycles. Δ t is duration, C t the accumulated/transferred charge. Other column variables are Ohmic work W = IV; battery temperature change Δ T = T d , c T 0 (subscripts denote end of discharge d, end of charge c, and start of charge/discharge 0); average ambient temperature T during the discharge/charge; discharge cut-off voltage V e n d ; Thevenin internal impedance Z = V O C V C C I determined via open-circuit voltage V O C (measured as V B with switch open), closed-circuit voltage V C C (measured as V B with resistive load R l o a d connected to battery), and current I, see Figure 2; and State of Health SoH. Due to varying depths of discharge DoD from cycle to cycle, values at 30 min of discharge are used for Z and SoH in Table A1. Similar trends were observed for all batteries tested [38]. As in Table 1 and Table 2, a “Half cycle” is a cycle comprising of a charge step only.
Table A1. Monitored and processed parameters for the Li-ion battery #2 (initial discharge current: ~5 A, charge current: 3 A).
Table A1. Monitored and processed parameters for the Li-ion battery #2 (initial discharge current: ~5 A, charge current: 3 A).
DischargeCharge
N Δ t
h
C t
Ah
W
Wh
Δ T
degC
T
degC
V e n d
V
Z 0.5 h r
Ω
SoH Δ t
hr
C t
Ah
W
Wh
Δ T
degC
T
degC
11.47−6.10−19.96.525.90.770.0421.003.4910.4741.40.227.6
21.79−8.37−28.610.224.82.430.0250.981.725.1420.20.826.7
3Half cycle1.333.9716.18.125.7
42.00−9.52−32.96.227.82.660.0310.981.283.8315.1−2.727.7
50.64−3.00−10.28.624.62.690.0810.951.815.4221.8−0.326.5
61.53−7.23−24.67.226.12.610.0260.971.474.4117.40.426.7
71.40−6.71−22.96.824.82.670.0320.992.738.1732.80.025.3
82.03−9.74−33.48.824.62.640.0281.00Missing data
91.36−6.39−21.56.324.92.580.0260.983.6010.8043.32.627.1
101.77−8.70−30.37.126.82.870.0240.990.611.817.0−1.727.7
111.09−5.18−17.59.223.52.420.0270.991.624.8619.20.926.2
121.82−8.52−29.06.926.61.880.0320.991.805.4021.52.826.3
13Half cycle2.136.3725.3−0.728.5
142.12−10.54−37.06.329.72.320.0281.030.992.9611.5−0.729.5
151.33−6.24−21.08.824.72.420.0400.981.534.5818.1−0.126.7
161.73−8.27−28.17.924.82.140.0251.001.464.3817.3−0.725.8
171.67−7.97−27.08.925.42.310.0291.001.815.4221.5−0.527.5
182.13−10.01−33.67.326.21.460.0271.021.283.8415.00.727.5
191.36−6.42−21.45.625.62.180.0260.991.474.4017.4−0.626.1
202.00−8.95−29.28.025.80.950.0271.011.695.0620.10.128.2
212.16−9.45−30.49.027.21.100.0281.021.273.7914.9−0.126.7
221.50−6.78−21.87.425.61.580.0310.981.414.2216.70.925.8
231.74−7.85−25.85.925.31.930.0210.951.504.5017.82.026.0
241.80−8.19−26.711.124.01.680.0250.971.735.1720.5−0.226.5
250.68−3.35−11.73.425.63.440.0250.974.1612.4750.92.126.3
262.69−12.37−40.68.025.01.560.0361.001.484.4317.52.426.8
272.44−8.84−25.76.924.90.570.0340.952.768.2733.13.126.7
28Half cycle1.775.3021.13.126.3
291.91−6.41−18.18.425.50.580.0450.933.089.2237.22.325.9
300.91−3.75−11.413.025.21.150.0740.963.339.9740.42.726.3
311.31−5.76−17.78.325.81.800.0510.953.219.6339.27.226.1
320.75−3.19−9.712.325.41.510.1290.941.695.0520.00.925.7

References

  1. Rahn, C.D.; Wang, C.Y. Battery Systems Engineering, 1st ed.; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2013. [Google Scholar] [CrossRef]
  2. Huggins, R.A. Energy Storage: Fundamentals, Materials and Applications, 2nd ed.; Springer: New York, NY, USA, 2010. [Google Scholar] [CrossRef]
  3. Ménard, L.; Fontès, G.; Astier, S. Dynamic energy model of a lithium-ion battery. Math. Comput. Simul. 2010, 81, 327–339. [Google Scholar] [CrossRef]
  4. Bresser, D.; Paillard, E.; Passerini, S. Lithium-ion batteries (LIBs) for medium- and large-scale energy storage: Current cell materials and components. In Advances in Batteries for Medium and Large-Scale Energy Storage; Elsevier Ltd.: Amsterdam, The Netherlands, 2014; Volume 1, pp. 125–211. [Google Scholar] [CrossRef]
  5. Goodenough, J.B. Energy storage materials: A perspective. Energy Storage Mater. 2015, 1, 158–161. [Google Scholar] [CrossRef]
  6. Bresser, D.; Paillard, E.; Passerini, S. Lithium-ion batteries (LIBs) for medium- and large-scale energy storage: Emerging cell materials and components. In Advances in Batteries for Medium and Large-Scale Energy Storage; Menictas, C., Skyllas-Kazacos, M., Lim, T.M., Eds.; Elsevier Ltd.: Amsterdam, The Netherlands, 2014; Volume 1, pp. 213–289. [Google Scholar] [CrossRef]
  7. Kurzweil, P. Lithium battery energy storage: State of the art including lithium-air and lithium-sulfur systems. In Electrochemical Energy Storage for Renewable Sources and Grid Balancing, 1st ed.; Moseley, P.T., Garche, G., Eds.; Elsevier B.V.: Amsterdam, The Netherlands, 2014; pp. 269–307. [Google Scholar] [CrossRef]
  8. Goodenough, J.B.; Kim, Y. Challenges for rechargeable batteries. J. Power Sources 2011, 196, 6688–6694. [Google Scholar] [CrossRef]
  9. Scrosati, B.; Garche, J. Lithium batteries: Status, prospects and future. J. Power Sources 2010, 195, 2419–2430. [Google Scholar] [CrossRef]
  10. Vetter, J.; Novák, P.; Wagner, M.R.; Veit, C.; Möller, K.C.; Besenhard, J.O.; Winter, M.; Wohlfahrt-Mehrens, M.; Vogler, C.; Hammouche, A. Ageing mechanisms in lithium-ion batteries. J. Power Sources 2005, 147, 269–281. [Google Scholar] [CrossRef]
  11. Jiang, J.; Zhang, C. Fundamentals and Applications of Lithium-Ion Batteries in Electric Drive Vehicles; John Wiley & Sons Singapore Pte Ltd.: Singapore, 2015. [Google Scholar] [CrossRef]
  12. Barsukov, Y. Challenges and Solutions in Battery Fuel Gauging; 2004. Available online: http://www.ti.com/lit/ml/slyp086/slyp086.pdf (accessed on 29 March 2019).
  13. Miranda, D.; Costa, C.M.; Lanceros-Mendez, S. Lithium ion rechargeable batteries: State of the art and future needs of microscopic theoretical models and simulations. J. Electroanal. Chem. 2015, 739, 97–110. [Google Scholar] [CrossRef]
  14. Ecker, M.; Gerschler, J.B.; Vogel, J.; Käbitz, S.; Hust, F.; Dechent, P.; Sauer, D.U. Development of a lifetime prediction model for lithium-ion batteries based on extended accelerated aging test data. J. Power Sources 2012, 215, 248–257. [Google Scholar] [CrossRef]
  15. Cordoba-Arenas, A.; Onori, S.; Guezennec, Y.; Rizzoni, G. Capacity and power fade cycle-life model for plug-in hybrid electric vehicle lithium-ion battery cells containing blended spinel and layered-oxide positive electrodes. J. Power Sources 2015, 278, 473–483. [Google Scholar] [CrossRef] [Green Version]
  16. Panchal, S.; McGrory, J.; Kong, J.; Dincer, I.; Agelin-Chaab, M.; Fraser, R.; Fowler, M. Cycling degradation testing and analysis of a LiFePO4 battery at actual conditions. Int. J. Energy Res. 2017, 41, 2565–2575. [Google Scholar] [CrossRef]
  17. Panchal, S.; Rashid, M.; Long, F.; Mathew, M.; Fraser, R.; Fowler, M. Degradation Testing and Modeling of 200 Ah LiFePO4 Battery; SAE Technical Paper; SAE International: Detroit, MI, USA, 2018. [Google Scholar] [CrossRef]
  18. Karnopp, D. Bond graph models for electrochemical energy storage: Electrical, chemical and thermal effects. J. Frankl. Inst. 1990, 327, 983–992. [Google Scholar] [CrossRef]
  19. Kondepudi, D.; Prigogine, I. Modern Thermodynamics: From Heat Engines to Dissipative Structures; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 1998. [Google Scholar]
  20. Esperilla, J.J.; Félez, J.; Romero, G.; Carretero, A. A full model for simulation of electrochemical cells including complex behavior. J. Power Sources 2007, 165, 436–445. [Google Scholar] [CrossRef]
  21. Ramadass, P.; Haran, B.; White, R.; Popov, B.N. Capacity fade of Sony 18650 cells cycled at elevated temperatures: Part II. Capacity fade analysis. J. Power Sources 2002, 112, 614–620. [Google Scholar] [CrossRef]
  22. Feng, X.; Gooi, H.B.; Chen, S. Capacity fade-based energy management for lithium-ion batteries used in PV systems. Electr. Power Syst. Res. 2015, 129, 150–159. [Google Scholar] [CrossRef]
  23. Ashwin, T.R.; Chung, Y.M.; Wang, J. Capacity fade modelling of lithium-ion battery under cyclic loading conditions. J. Power Sources 2016, 328, 586–598. [Google Scholar] [CrossRef] [Green Version]
  24. Long, L.; Wang, S.; Xiao, M.; Meng, Y. Polymer electrolytes for lithium polymer batteries. J. Mater. Chem. A 2016, 4, 10038–10069. [Google Scholar] [CrossRef]
  25. Schmalstieg, J.; Käbitz, S.; Ecker, M.; Sauer, D.U. A holistic aging model for Li(NiMnCo)O2 based 18650 lithium-ion batteries. J. Power Sources 2014, 257, 325–334. [Google Scholar] [CrossRef]
  26. Ecker, M.; Nieto, N.; Käbitz, S.; Schmalstieg, J.; Blanke, H.; Warnecke, A.; Sauer, D.U. Calendar and cycle life study of Li(NiMnCo)O2-based 18650 lithium-ion batteries. J. Power Sources 2014, 248, 839–851. [Google Scholar] [CrossRef]
  27. Käbitz, S.; Gerschler, J.B.; Ecker, M.; Yurdagel, Y.; Emmermacher, B.; André, D.; Mitsch, T.; Sauer, D.U. Cycle and calendar life study of a graphite|LiNi1/3Mn 1/3Co1/3O2 Li-ion high energy system. Part A: Full cell characterization. J. Power Sources 2013, 239, 572–583. [Google Scholar] [CrossRef]
  28. Bryant, M.D.; Khonsari, M.M.; Ling, F.F. On the thermodynamics of degradation. Proc. R. Soc. A: Math. Phys. Eng. Sci. 2008, 464, 2001–2014. [Google Scholar] [CrossRef] [Green Version]
  29. Strutt, J.W.; Rayleigh, B. The Theory of Sound. Nature 1877, 2. [Google Scholar] [CrossRef]
  30. Onsager, L. Reciprocal Relations in Irreversible processes. Phys. Rev. 1931, 37, 405. [Google Scholar] [CrossRef]
  31. Doelling, K.L.; Ling, F.F.; Bryant, M.D.; Heilman, B.P. An experimental study of the correlation between wear and entropy flow in machinery components. J. Appl. Phys. 2000, 88, 2999–3003. [Google Scholar] [CrossRef]
  32. Ling, F.F.; Bryant, M.D.; Doelling, K.L. On irreversible thermodynamics for wear prediction. Wear 2002, 253, 1165–1172. [Google Scholar] [CrossRef]
  33. Bryant, M.D. Entropy and Dissipative Processes of Friction and Wear. FME Trans. 2009, 37, 55–60. [Google Scholar]
  34. Bryant, M.D. On Constitutive Relations for Friction from Thermodynamics and Dynamics. J. Tribol. 2016, 138, 041603. [Google Scholar] [CrossRef]
  35. Naderi, M.; Khonsari, M.M. An experimental approach to low-cycle fatigue damage based on thermodynamic entropy. Int. J. Solids Struct. 2010, 47, 875–880. [Google Scholar] [CrossRef] [Green Version]
  36. Amiri, M.; Modarres, M. An entropy-based damage characterization. Entropy 2014, 16, 6434–6463. [Google Scholar] [CrossRef]
  37. Amiri, M.; Naderi, M.; Khonsari, M.M. An experimental approach to evaluate the critical damage. Int. J. Damage Mech. 2011, 20, 89–112. [Google Scholar] [CrossRef]
  38. Osara, J.A. The Thermodynamics of Degradation. Ph.D. Thesis, The University of Texas at Austin, Austin, TX, USA, May 2017. [Google Scholar]
  39. DeHoff, R.T. Thermodynamics in Material Science, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2006. [Google Scholar]
  40. Callen, H.B. Thermodynamics and an Introduction to Thermostatistics; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 1985. [Google Scholar]
  41. De Groot, S.R. Thermodynamics of Irreversible Processes; North-Holland Publishing Company: Amsterdam, The Netherland, 1951. [Google Scholar] [CrossRef]
  42. Prigogine, I. Introduction to Thermodynamics of Irreversible Processes; Charles C Thomas: Springfield, IL, USA, 1955. [Google Scholar]
  43. Bejan, A. Advanced Engineering Thermodynamics, 3rd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1997. [Google Scholar]
  44. Moran, M.J.; Shapiro, H.N. Fundamentals of Engineering Thermodynamics, 5th ed.; Wiley: Hoboken, NJ, USA, 2004. [Google Scholar]
  45. Burghardt, M.D.; Harbach, J.A. Engineering Thermodynamics, 4th ed.; Harper Collins College Publishers: Glenview, IL, USA, 1993. [Google Scholar]
  46. Reynier, Y.; Yazami, R.; Fultz, B. The entropy and enthalpy of lithium intercalation into graphite. J. Power Sources 2003, 119–121, 850–855. [Google Scholar] [CrossRef]
  47. Yazami, R.; Maher, K. Thermodynamics of lithium-ion batteries. In Lithium-Ion Batteries: Advances and Applications; Pistoia, G., Ed.; Elsevier: Amsterdam, The Netherlands, 2014. [Google Scholar] [CrossRef]
  48. Bejan, A. The method of entropy generation minimization. In Energy and the Environment; Sage Publications: Thousand Oaks, CA, USA, 1990; pp. 11–22. [Google Scholar]
  49. Pal, R. Demystification of the Gouy-Stodola theorem of thermodynamics for closed systems. Int. J. Mech. Eng. Educ. 2017, 45, 142–153. [Google Scholar] [CrossRef]
  50. Fernández, I.J.; Calvillo, C.F.; Sánchez-Miralles, A.; Boal, J. Capacity fade and aging models for electric batteries and optimal charging strategy for electric vehicles. Energy 2013, 60, 35–43. [Google Scholar] [CrossRef]
  51. Lide, D.R. CRC Handbook of Chemistry and Physics, 75th ed.; CRC Press: Boca Raton, FL, USA, 1994. [Google Scholar]
Figure 1. Photo of a portion of actual experimental set-up used in Li-ion battery measurements, showing four Li-ion batteries, a commercial charger, high-power discharge resistors, connector cables, and thermocouple wires.
Figure 1. Photo of a portion of actual experimental set-up used in Li-ion battery measurements, showing four Li-ion batteries, a commercial charger, high-power discharge resistors, connector cables, and thermocouple wires.
Inventions 04 00023 g001
Figure 2. Discharge–charge cycle circuit for the lithium-ion batteries.
Figure 2. Discharge–charge cycle circuit for the lithium-ion batteries.
Inventions 04 00023 g002
Figure 3. Monitored parameters during cycling (cycle 6 data) showing 1.5 h discharge starting at ~5 A, followed by 0.4 h settling and 1.5 h constant-current charge at 3 A.
Figure 3. Monitored parameters during cycling (cycle 6 data) showing 1.5 h discharge starting at ~5 A, followed by 0.4 h settling and 1.5 h constant-current charge at 3 A.
Inventions 04 00023 g003
Figure 4. Phenomenological entropy components (Ohmic and Electro-Chemico-Thermal (ECT)) vs. measured charge transfer during discharge (solid curves) and charge (dashed curves).
Figure 4. Phenomenological entropy components (Ohmic and Electro-Chemico-Thermal (ECT)) vs. measured charge transfer during discharge (solid curves) and charge (dashed curves).
Inventions 04 00023 g004
Figure 5. Entropy generation component rates (Ohmic, ECT, and reversible Gibbs) over time, during discharge and charge.
Figure 5. Entropy generation component rates (Ohmic, ECT, and reversible Gibbs) over time, during discharge and charge.
Inventions 04 00023 g005
Figure 6. Entropy generation S and components during cycle 6′s discharge and charge steps. The region between reversible entropy and phenomenological entropy curves represents entropy generation.
Figure 6. Entropy generation S and components during cycle 6′s discharge and charge steps. The region between reversible entropy and phenomenological entropy curves represents entropy generation.
Inventions 04 00023 g006
Figure 7. (a) 3D plots and linear surface fits of Charge (vertical axes) vs. Ohmic and ECT entropies (horizontal axes) during discharge (red plane, curve ‘D’ starts from upper right corner) and charge (green plane, curve ‘C’ starts from lower left corner) steps of cycle 6, indicating a linear dependence on two active processes. (b) The end projection of the planes and curves of (a), and a visual of the goodness of fit. Axes are not to scale.
Figure 7. (a) 3D plots and linear surface fits of Charge (vertical axes) vs. Ohmic and ECT entropies (horizontal axes) during discharge (red plane, curve ‘D’ starts from upper right corner) and charge (green plane, curve ‘C’ starts from lower left corner) steps of cycle 6, indicating a linear dependence on two active processes. (b) The end projection of the planes and curves of (a), and a visual of the goodness of fit. Axes are not to scale.
Inventions 04 00023 g007
Figure 8. (a) 3D scatter of cyclic phenomenological C p h e n (vertical axis) vs. Ohmic and ECT entropies (horizontal axes) for all discharge steps (purple dots on red plane) and charge steps (white dots on green plane) shows that all 32 cycles lay on cycle 1’s DEG planes, a visual of Equations (52) and (53); (b) end projection of (a) shows a linear relationship. Axes are not to scale.
Figure 8. (a) 3D scatter of cyclic phenomenological C p h e n (vertical axis) vs. Ohmic and ECT entropies (horizontal axes) for all discharge steps (purple dots on red plane) and charge steps (white dots on green plane) shows that all 32 cycles lay on cycle 1’s DEG planes, a visual of Equations (52) and (53); (b) end projection of (a) shows a linear relationship. Axes are not to scale.
Inventions 04 00023 g008
Figure 9. Capacity fade and components for cycle 6′s discharge and charge. The region between reversible and phenomenological charge transfer curves is capacity fade Δ C .
Figure 9. Capacity fade and components for cycle 6′s discharge and charge. The region between reversible and phenomenological charge transfer curves is capacity fade Δ C .
Inventions 04 00023 g009
Table 1. Processed Gibbs energy and entropy parameters for Li-ion battery #2 (Discharge rate: ~5 A, Charge rate: 3 A). Cycle 6 (in bold) is used in the breakdown in this section.
Table 1. Processed Gibbs energy and entropy parameters for Li-ion battery #2 (Discharge rate: ~5 A, Charge rate: 3 A). Cycle 6 (in bold) is used in the breakdown in this section.
DischargeCharge
N G Ω
Wh
G V T
Wh
S Ω
Wh/K
S V T
Wh/K
S r e v
Wh/K
G Ω
Wh
G V T
Wh
S Ω
Wh/K
S V T
Wh/K
S r e v
Wh/K
1−19.93−3.00−0.07−0.010−0.1141.401.000.140.0030.13
2−28.62−2.22−0.09−0.007−0.1320.150.550.070.0020.07
3Half cycle16.140.080.050.0000.05
4−32.85−2.72−0.11−0.009−0.1515.080.430.050.0010.05
5−10.20−0.73−0.03−0.002−0.0521.750.380.070.0010.07
6−24.57−1.64−0.08−0.005−0.1117.440.500.060.0020.06
7−22.92−1.73−0.08−0.006−0.1032.750.800.110.0030.11
8−33.43−2.63−0.11−0.009−0.16Missing data
9−21.52−1.39−0.07−0.005−0.1043.260.710.140.0020.14
10−30.25−1.55−0.10−0.005−0.136.990.160.020.0010.02
11−17.48−1.13−0.06−0.004−0.0719.230.560.060.0020.06
12−29.00−2.89−0.10−0.010−0.1321.490.620.070.0020.07
13Half cycle25.340.870.080.0030.08
14−37.01−3.11−0.12−0.010−0.1611.490.360.040.0010.04
15−20.97−2.00−0.07−0.007−0.1018.070.570.060.0020.06
16−28.13−1.95−0.09−0.006−0.1217.290.510.060.0020.06
17−26.99−2.28−0.09−0.008−0.1221.450.650.070.0020.07
18−33.60−3.65−0.11−0.012−0.1615.000.460.050.0020.05
19−21.43−1.62−0.07−0.005−0.1017.430.420.060.0010.06
20−29.16−3.74−0.10−0.012−0.1520.050.650.070.0020.06
21−30.39−4.53−0.10−0.015−0.1614.900.240.050.0010.05
22−21.75−2.56−0.07−0.008−0.1116.700.270.060.0010.05
23−25.83−1.72−0.09−0.006−0.1317.790.510.060.0020.06
24−26.72−2.70−0.09−0.009−0.1320.470.370.070.0010.07
25−11.73−0.31−0.04−0.001−0.0550.850.560.170.0020.16
26−40.64−5.32−0.13−0.017−0.2117.520.500.060.0020.06
27−25.68−5.24−0.08−0.017−0.1833.060.510.110.0020.11
28Half cycle21.130.260.070.0010.07
29−18.09−3.55−0.06−0.012−0.1437.180.490.120.0020.12
30−11.39−1.53−0.04−0.005−0.0740.440.440.130.0010.13
31−17.74−2.32−0.06−0.008−0.0939.200.540.130.0020.13
32−9.67−1.21−0.03−0.004−0.0619.960.740.070.0020.07
SUMMARY/TOTAL
−2.34−0.234−3.47 2.430.0532.39
Table 2. Processed Degradation-Entropy Generation (DEG) and Coulomb-Counted capacity fade parameters for Li-ion battery #2 (Discharge rate: ~5 A, Charge rate: 3 A). Cycle 6 (in bold) is used in the Capacity Fade discussion.
Table 2. Processed Degradation-Entropy Generation (DEG) and Coulomb-Counted capacity fade parameters for Li-ion battery #2 (Discharge rate: ~5 A, Charge rate: 3 A). Cycle 6 (in bold) is used in the Capacity Fade discussion.
DischargeCharge
N C p h e n
Ah
C r e v
Ah
Δ C D E G
Ah
C C C
Ah
Δ C C C
Ah
C p h e n
Ah
C r e v
Ah
Δ C D E G
Ah
C C C
Ah
1−6.5−7.61.1−6.10.010.710.10.610.5
2−7.7−9.31.6−8.4−2.35.35.00.35.2
3Half cycle3.83.80.04.0
4−9.4−10.41.0−9.5−3.43.83.70.13.8
5−2.5−3.30.8−3.03.15.35.20.15.4
6−6.7−8.01.3−7.2−1.14.64.30.34.4
7−6.8−7.30.5−6.7−0.68.47.90.58.2
8−9.4−10.61.2−9.7−3.6Missing data
9−5.9−7.11.2−6.4−0.310.610.40.210.8
10−8.2−9.21.0−8.7−2.61.51.50.01.6
11−5.0−5.70.7−5.20.94.64.60.04.8
12−8.8−9.50.7−8.5−2.45.35.20.15.4
13Half cycle6.16.10.06.3
14−10.3−11.00.7−10.5−4.43.02.90.13.0
15−6.2−6.90.7−6.2−0.14.64.40.24.6
16−7.6−9.01.4−8.3−2.24.64.20.44.4
17−7.8−8.70.9−8.0−1.95.35.20.15.4
18−9.8−11.11.3−10.0−3.93.83.70.13.8
19−5.9−7.11.2−6.4−0.34.64.30.34.4
20−9.0−10.41.4−9.0−2.95.34.90.45.1
21−9.4−11.21.8−9.5−3.43.83.70.13.8
22−6.3−7.81.5−6.8−0.74.64.10.54.2
23−7.6−9.01.4−7.9−1.84.64.30.34.5
24−7.9−9.41.5−8.2−2.15.35.00.35.2
25−3.2−3.50.3−3.42.812.912.00.912.5
26−11.9−14.02.1−12.4−6.34.64.30.34.4
27−8.0−12.74.7−8.8−2.78.48.00.48.3
28Half cycle5.35.10.25.3
29−6.0−9.93.9−6.4−0.39.18.90.29.2
30−3.6−4.71.1−3.82.49.89.60.210.0
31−5.5−6.81.3−5.80.39.99.30.69.6
32−2.8−3.91.1−3.22.95.34.90.45.1
SUMMARY/TOTAL
−205.7−245.039.3−213.8 185.0176.88.2183.2

Share and Cite

MDPI and ACS Style

Osara, J.A.; Bryant, M.D. A Thermodynamic Model for Lithium-Ion Battery Degradation: Application of the Degradation-Entropy Generation Theorem. Inventions 2019, 4, 23. https://doi.org/10.3390/inventions4020023

AMA Style

Osara JA, Bryant MD. A Thermodynamic Model for Lithium-Ion Battery Degradation: Application of the Degradation-Entropy Generation Theorem. Inventions. 2019; 4(2):23. https://doi.org/10.3390/inventions4020023

Chicago/Turabian Style

Osara, Jude A., and Michael D. Bryant. 2019. "A Thermodynamic Model for Lithium-Ion Battery Degradation: Application of the Degradation-Entropy Generation Theorem" Inventions 4, no. 2: 23. https://doi.org/10.3390/inventions4020023

APA Style

Osara, J. A., & Bryant, M. D. (2019). A Thermodynamic Model for Lithium-Ion Battery Degradation: Application of the Degradation-Entropy Generation Theorem. Inventions, 4(2), 23. https://doi.org/10.3390/inventions4020023

Article Metrics

Back to TopTop