Page 166 - Lindens Handbook of Batteries
P. 166

MATHEMATICAL MODELING OF BATTERIES        6.19


                                                                                      e
                                         Specific energy (kWh/kg)          a          c  d
                                            –1
                                          10









                                          10 –2                               b
                                                    10 –2         10 –1         10 0          10 1
                                                              Specific power (kW/kg)
                                        FiguRE  6.12  Simulated  Ragone  plots  constructed  using  a  mechanistic  model.  (a)  Model
                                        predictions show that the original design cell delivers very low specific energies at high power
                                        applications. (b) Increasing the electronic conductivity of the cathode matrix (e.g., by addition
                                        of conductive carbon) shows some improvement. (c) Further enhancement in electronic conduc-
                                        tivity shows little change in cell performance. (d) Diffusion limitations are relaxed by increasing
                                        the  solid-phase  diffusion  coefficient  at  the  cathode  (e.g.,  by  doping).  (e)  Further  limitations
                                        within the solid phase are eliminated by reducing the particle size to a few nanometers.


                                The left side of Eq. (6.47) denotes the rate at which energy is consumed or generated per unit volume.
                                The first term on the right side represents heat transfer by conduction following Fourier’s law. The
                                term q refers to heat generated or consumed due to the reactions that take place during the operation
                                of the battery. Typically, for an electrochemical reaction, this term is expressed as follows: 19

                                                       i ∂         ∂ E    ∂φ  ∂φ
                                                   q =  2  -  φ  - φ    E - T  j   -  i  1  - i i  2     (6.48)
                                                                       
                                                       x ∂    1  2    j  ∂ T   1  x ∂  2  x ∂
                                                                       
                                The first component of q represents the heat generated from the charge-transfer reaction, the second
                                term is the Joule heat generated due to the current flow across the solid matrix, and the last term is the
                                corresponding value for current flow in the electrolyte. Additional complexities, such as heat transfer
                                due to differential phase changes, radiation effects, etc., can be treated by incorporating the amount of
                                                                              ∂ E
                                heat generated from such phenomena in Eq. (6.48). The term  T  j   corrects the open circuit voltage
                                                                               T ∂
                                for changes in entropy with temperature. Open circuit voltages measured at different temperatures
                                may be used to evaluate this term empirically. Changes in other properties such as the diffusivity or
                                conductivity with temperature are often approximated by the Arrhenius equation

                                                                    E   1  1   
                                                         Φ  Φ =  ref  exp  -  a    -            (6.49)
                                                                       
                                                                     RT  T ref    
                                where Φ may represent D , κ , D , etc., Φ  represents the corresponding property measured at
                                                                 ref
                                                       eff
                                                   ef f
                                                           s
                                the reference temperature (T ), and E  is the activation energy.
                                                            a
                                                     ref
                                   Figure 6.13 illustrates the change in battery performance with different degrees of convective
                                cooling during a 3C-rate discharge of a lithium-ion cell. For the adiabatic case, increase in the cell
                                temperature favors higher reaction rates and enhances transport within the electrolyte following
                                [Eq. (6.49)]. The model predicts that without a cooling system in place, the difference in cell
   161   162   163   164   165   166   167   168   169   170   171