# Patent

## Title
Method and system for battery state of charge estimation

## Abstract

A powertrain having a Lithium-ion (Li-ion) traction battery including solid active particles is operated according to a state of charge of the battery based on an estimated Li-ion concentration profile of a representative solid active particle of a reduced-order model of the battery. The concentration profile is estimated according to the reduced-order model from the measured voltage and current of the battery.

## Background
TECHNICAL FIELD
The present invention relates to estimating the state of charge (SOC) of a battery.
BACKGROUND
An electric vehicle includes a traction battery. The battery has a state of charge (SOC) indicative of the current charging state of the battery. Knowledge of the battery SOC enables usable battery energy and battery power capability to be calculated.
The SOC of a battery may be estimated indirectly. One approach to estimating the battery SOC involves battery current integration (i.e., coulomb counting). A problem with such battery current integration approaches is that highly accurate current sensors are required to ensure the accuracy of the SOC estimation. Another approach uses battery voltage information to estimate the SOC. A problem with such battery voltage translation approaches is that many of them are either map-based or simplified equivalent circuit model based.
A battery is operated around a given battery SOC set-point in charge sustaining driving mode. The battery SOC estimation provided by the noted existing estimation approaches may be valid in this case as the long time duration dynamics are negligible. However, slow diffusion dynamics are dominant in charge depleting driving mode. The battery SOC estimation provided by the noted existing estimation approaches will not be valid due to the high non-linearity of the slow dynamics and dynamic components along wide frequency ranges.
Another battery SOC estimation approach should be developed which can be used over wide SOC, temperature ranges, and various vehicle driving conditions.
SUMMARY
In embodiments, a powertrain having a Lithium-ion (Li-ion) traction battery pack is operated according to values of performance variables of the battery. The performance variables include a state of charge (SOC) of the battery and available battery power limits. The performance variables are estimated from the state variables, represented by Li-ion concentrations in a reduced-order electrochemical battery model. The battery model is further reduced using uneven discretization of the derived ordinary differential equation of the model. A method to estimate the battery SOC based on the estimated Li-ion concentration profile representing electrochemical dynamics of the positive and negative electrodes in the battery cell is provided. The concentration profile is estimated from the measured voltage and current of the battery using Extended Kalman Filter (EKF). The EKF is constructed based on the reduced-order battery model.
In embodiments, a Lithium-ion (Li-ion) traction battery including solid active particles is operated according to a SOC of the battery based on an estimated Li-ion concentration profile representing electrochemical dynamics of a battery cell or pack. Electrochemical dynamics are modeled using a diffusion equation, which captures the electrochemical dynamics of both electrodes, i.e., a positive electrode and a negative electrode. The reduced-order model adopts the structure of an electrode averaged model, but the model validity is extended to wider ranges of operating conditions by adding medium-to-slow dynamics information to the state variables, which are represented by effective Li-ion concentration profiles. The effective Li-ion concentration profile is estimated according to the model from the measured voltage and current of the battery.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 illustrates a block diagram of a power storage system of an electric vehicle in accordance with embodiments of the present invention;
FIG. 2 illustrates a representation of a full-order electrochemical battery model and a reduced-order electrochemical battery model (i.e., an electrode-averaged battery model) derived from the full-order battery model pursuant to a battery state of charge (SOC) estimation process in accordance with embodiments of the present invention;
FIG. 3A illustrates a representation of an even discretization of the Li-ion concentration profile along the electrode solid particle radius for the reduced-order battery model;
FIG. 3B illustrates a representation of an uneven discretization of the Li-ion concentration profile along the electrode solid particle radius for the reduced-order battery model;
FIG. 4 illustrates a graph of open circuit voltage of the battery versus battery SOC;
FIG. 5 illustrates a block diagram illustrating operation of a method and a system for estimating the SOC of a battery using the reduced-order battery model in accordance with embodiments of the present invention;
FIG. 6A illustrates a graph depicting the battery SOC over time including a comparison between the actual battery SOC and the battery SOC calculated using the reduced-order battery model in accordance with embodiments of the present invention;
FIG. 6B illustrates a graph depicting the concentration profile versus normalized electrode solid particle radius at a time instant including a comparison of the average concentration and the concentration estimated in accordance with embodiments of the present invention; and
FIG. 6C illustrates a graph depicting the concentration profile versus normalized electrode solid particle radius at another time instant including a comparison of the average concentration and the concentration estimated in accordance with embodiments of the present invention.



DETAILED DESCRIPTION
Detailed embodiments of the present invention are disclosed herein; however, it is to be understood that the disclosed embodiments are merely exemplary of the invention that may be embodied in various and alternative forms. The figures are not necessarily to scale; some features may be exaggerated or minimized to show details of particular components. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a representative basis for teaching one skilled in the art to variously employ the present invention.
Referring now to FIG. 1, a block diagram of a power storage system of an electric vehicle 10 such as a hybrid electric vehicle in accordance with embodiments of the present invention is shown. The power storage system includes a traction battery 12. Battery 12 is configured to provide electrical power to an electric machine 14 such as a motor/generator that converts the electrical power to motive power for propelling vehicle 10. Battery 12 includes any number and/or type of energy storage cells. For instance, battery 12 is a Lithium-Ion (Li-ion) battery having a plurality of Li-ion power storage cells.
A controller 16 is configured to monitor parameters associated with battery 12 and control the operation of the battery based on the parameters and other factors. Controller 16 is further configured to estimate the state of charge (SOC) of battery 12 pursuant to a battery SOC estimation process in accordance with embodiments of the present invention. The SOC estimation process includes estimating the SOC of battery 12 from the (measured) battery voltage in real time using a reduced-order electrochemistry model of the battery.
Objectives of the battery SOC estimation process include: 1) generating a reduced-order electrochemical battery model (i.e., an electrode-averaged battery model) which can capture the dynamics of battery 12 with sufficient accuracy; 2) minimizing the number of parameters for tuning the reduced-order battery model; and 3) estimating battery SOC from the measured battery voltage in real-time using the reduced-order battery model.
The derivation of the reduced-order or electrode-averaged model of battery 12 (“reduced-order battery model”) is based on the understanding of the full-order electrochemical model of battery 12. This full-order battery model is manipulated to the reduced-order model with the assumption that each electrode (i.e., cathode and anode) can be represented by the averaged electrode particle. This reduced model is then further reduced with the assumption that the diffusion dynamics can be represented by a representative diffusion dynamics accounting for the combined dynamics of the cathode and anode. This makes the resulting model size to be one-half of a model with two single particles. The resulting diffusion dynamics are the averaged diffusion dynamics of the cathode and the anode.
Separately, instantaneous voltage drops by the electrolyte and the over-potential, which is computed from the Butler-Volmer current density equation, are combined in a single internal resistance term. These assumptions make the reduced-order or electrode-averaged battery model even simpler as open potential and electrolyte potential terms do not need to be computed separately.
In turn, Li-ion concentration profiles, which represent the battery cell dynamics, are estimated from input current profiles of output voltage measurement profiles using Extended Kalman Filter (EKF) based on the reduced-order battery model. The battery SOC is then computed from the estimated Li-ion concentration profiles.
Referring now to FIG. 2, an illustration of the full-order electrochemical battery model 20 and an illustration of the reduced-order electrochemical battery model 30 derived from full-order battery model 20 are shown. Again, reduced-order battery model 30 is derived from full-order battery model 20 as part of a battery SOC estimation process in accordance with embodiments of the present invention. Then, reduced-order battery model 30 is further simplified to estimate battery electrochemical states represented by Li-ion concentration profiles and to estimate the battery SOC. Although the simplification of reduced-order battery model 30 may lose prediction accuracy under high charging and discharging conditions, important diffusion dynamics in the solid particle are described in reduced-order battery model 30.
Again, with reference to FIG. 2, dynamics of the positive and negative electrodes may be represented by the combined dynamics of each electrode. As shown in FIG. 2, full-order battery model 20 includes a plurality of active material, solid particles 22 immersed within an electrolyte solution 24 in the anode and a plurality of active material, solid particles 26 immersed within an electrolyte solution 28 in the cathode. As known to one of ordinary skill, a Li-ion 29 (only one shown for simplicity) moves along a path from solid particles 22 of the anode to a path of solid particles 26 of the cathode to thereby generate electric current I upon an application of a voltage V across the battery.
As further shown in FIG. 2, reduced-order battery model 30 is simplified relative to full-order battery model 20 by including an averaged solid particle 32 immersed in an electrolyte solution 34 in the anode and an averaged solid particle 36 immersed in an electrolyte solution 38 in the cathode. Averaged solid particle 32 is representative of the plurality of solid particles 22 in the anode of full-order battery model 20 and averaged solid particle 36 is representative of the plurality of solid particles 26 in the cathode of full-order battery model 20.
The derivation of reduced-order battery model 30 from full-order battery model 20 will now be described in detail. As known to one of ordinary skill, full-order battery model 20 is represented by the following equations:



  
    
      
        
          
            
              ∇
              →
            
            k
          
          ⁢
          
            X
            eff
          
        
        ⁢
        
          
            
              ∇
              x
            
            →
          
          ⁢
          
            ∅
            e
          
        
      
      +
      
        
          
            
              ∇
              x
            
            →
          
          ⁢
          
            K
            D
            eff
          
        
        ⁢
        
          
            
              ∇
              x
            
            →
          
          ⁢
          1
        
        ⁢
        n
        ⁢
        
          
        
        ⁢
        
          c
          e
        
      
    
    =
    
      -
      
        j
        Li
      
    
  




  
    
      
        
          
            ∇
            →
          
          x
        
        ⁢
        
          σ
          eff
        
      
      ⁢
      
        
          
            ∇
            →
          
          x
        
        ⁢
        
          ∅
          s
        
      
    
    =
    
      j
      Li
    
  




  
    
      
        ∂
        
          ɛ
          
            e
            ⁢
            
              
            
            ⁢
            
              c
              e
            
          
        
      
      
        ∂
        t
      
    
    =
    
      
        
          
            ∇
            →
          
          x
        
        ⁢
        
          (
          
            
              D
              e
              eff
            
            ⁢
            
              
                
                  ∇
                  →
                
                x
              
              ⁢
              
                c
                e
              
            
          
          )
        
      
      +
      
        
          
            1
            -
            
              t
              0
            
          
          F
        
        ⁢
        
          j
          Li
        
      
    
  




  
    
      
        ∂
        
          c
          s
        
      
      
        ∂
        t
      
    
    =
    
      
        
          ∇
          →
        
        r
      
      ⁢
      
        (
        
          
            D
            s
          
          ⁢
          
            
              
                ∇
                →
              
              r
            
            ⁢
            
              c
              s
            
          
        
        )
      
    
  



where φ is the electric potential, c is the Li-ion concentration, subscript s and e represent the active electrode solid particle and the electrolyte respectively, σeff is the effective electrical conductivity of the electrode, κeff is the effective electrical conductivity of the electrolyte, κDeff is the liquid junction potential term, Ds is the diffusion coefficient of Li-ion in the electrode, Deeff is the effective diffusion coefficient of Li-ion in the electrolyte, t0 is the transference number, and F is the Faraday constant.
The Butler-Volmer current density equation is as follows:




  
    j
    Li
  
  =
  
    
      a
      s
    
    ⁢
    
      
        j
        0
      
      ⁡
      
        [
        
          
            exp
            ⁡
            
              (
              
                
                  
                    
                      α
                      a
                    
                    ⁢
                    F
                  
                  RT
                
                ⁢
                η
              
              )
            
          
          -
          
            exp
            ⁡
            
              (
              
                
                  -
                  
                    
                      
                        α
                        c
                      
                      ⁢
                      F
                    
                    RT
                  
                
                ⁢
                η
              
              )
            
          
        
        ]
      
    
  




where, αa is the transfer coefficient for anodic reaction, αc is the transfer coefficient for cathodic reaction, R is the gas constant, T is the temperature, η=φs−φe−U(cse) is the over potential at the solid-electrolyte interface at an active solid particle, and j0=k(ce)αa(cs,max−cse)αa(cse)αc.
The battery terminal voltage V is according to the following equation:




  V
  =
  
    
      
        φ
        s
      
      ⁡
      
        (
        
          x
          =
          L
        
        )
      
    
    -
    
      
        φ
        s
      
      ⁡
      
        (
        
          x
          =
          0
        
        )
      
    
    -
    
      
        
          R
          f
        
        A
      
      ⁢
      I
    
  




Reduced-order battery model 30 is derived from full-order battery model 20 by making the following assumptions as shown in FIG. 2:
∇xØs is 0 (the gradient of the electric potential in the solid with respect to the direction from the anode to the cathode (i.e., the x-direction) is zero);
∇xØe is 0 (the gradient of the electric potential in the electrolyte with respect to the x-direction is zero); and
∇xce is 0 (the gradient of the Li-ion concentration in the electrolyte (ce) with respect to the x-direction is zero).
As a result of these assumptions, reduced-order battery model 30 is represented by the following ordinary differential equation (ODE) of each electrode (1):




  
    
      
        
          
            ∂
            
              c
              s
            
          
          
            ∂
            t
          
        
        =
        
          
            
              ∇
              →
            
            r
          
          ⁢
          
            (
            
              
                D
                s
              
              ⁢
              
                
                  
                    ∇
                    →
                  
                  r
                
                ⁢
                
                  c
                  s
                
              
            
            )
          
        
      
    
    
      
        (
        1
        )
      
    
  




The derived ODE represents the electrochemical dynamics of each electrode.
With the assumption that the electrochemical dynamics of the positive and negative electrodes can be combined and then represented by one ODE, an effective Li-ion concentration profile representing electrochemical dynamics combining both electrodes is estimated by the ODE. In other words, the Li-ion concentration profile for each of averaged solid particles 32 and 36 is considered to be the same for simplification.
The battery terminal voltage is according to the following equation (2):

V=U(cse)−R0i  (2)

where Cse is the effective Li-ion concentration at the solid-electrolyte interface and R0 is the internal resistance term which depends on the open circuit voltage and is a known quantity of battery 12.
The ODE (1) describes the Li-ion diffusion dynamics and other medium-to-slow electrochemical dynamics of battery 12. The states of the equation (1) are Li-ion concentrations along the discretized particle radius and are given by the following matrix equation (3):

cs=[cs,1 . . . cs,Mr-1]T  (3)

The control input is battery input current represented by the following equation (4):

uk=ik  (4)

The (measured) system output (y) is the battery terminal voltage and the battery terminal voltage expression is represented by the following equation (5):

y=U(cse)−R0i  (5)

Equation (5) for the measured battery terminal voltage corresponds to the proposed terminal voltage equation (2) (V=U(cse)−R0i).
As described herein, only one ordinary differential equation (i.e., the state-space equation (1)) is solved for both of the anode and the cathode. Thus, the diffusion dynamics from equation (1) represent the averaged sense of the diffusion dynamics of both the anode and the cathode. The battery terminal voltage is reduced to the terminal voltage equation (2) with the assumption that the over-potential and voltage drops by electrolyte can be represented by the internal resistance R0. This assumption enables to exclude over-potential calculations associated with full-order battery model 20 and thus the resulting reduced-order battery model 30 is even simpler.
Again, the state vector according to the matrix equation (3) (cs=[cs,1 . . . cs,Mr-1]T) represents the Li-ion concentration profiles along the electrode solid particle radius. The averaged Li-ion concentration profile is calculated according to the following equation (other expressions can also be used):




  
    
      c
      _
    
    s
  
  =
  
    
      
        
          c
          
            s
            ,
            1
          
        
        ⁢
        
          r
          1
          3
        
      
      +
      
        
          ∑
          
            i
            =
            1
          
          
            Mr
            -
            1
          
        
        ⁢
        
          
            3
            8
          
          ⁢
          
            (
            
              
                c
                
                  s
                  ,
                  i
                
              
              +
              
                c
                
                  s
                  ,
                  
                    i
                    +
                    1
                  
                
              
            
            )
          
          ⁢
          
            
              (
              
                
                  r
                  i
                
                +
                
                  r
                  
                    i
                    +
                    1
                  
                
              
              )
            
            2
          
          ⁢
          
            (
            
              
                r
                i
              
              -
              
                r
                
                  i
                  +
                  1
                
              
            
            )
          
        
      
    
    
      r
      
        Mr
        -
        1
      
      3
    
  




The SOC of battery 12 may then be estimated according to the following equation:



  
    
      S
      ⁢
      
        O
        ^
      
      ⁢
      C
    
    =
    
      
    
    ⁢
    
      
        
          
            θ
            _
          
          s
        
        -
        
          θ
          
            0
            ⁢
            %
          
        
      
      
        
          θ
          
            100
            ⁢
            %
          
        
        -
        
          θ
          
            0
            ⁢
            %
          
        
      
    
  




  where




  
    
      
        θ
        _
      
      s
    
    =
    
      
        
          c
          _
        
        s
      
      
        c
        
          s
          ,
          
            ma
            ⁢
            
              
            
            ⁢
            x
          
        
      
    
  



with θ0% at 0% SOC; with θ100% at 100% SOC; and cs,max being maximum Li-ion concentration.
Referring now to FIGS. 3A and 3B, with continual reference to FIG. 2, representations of an even discretization 40 and an uneven discretization 42 of the Li-ion concentration profile along the electrode solid particle radius for reduced-order battery model 30 are respectively shown. Objectives of uneven discretization 42 include obtaining the same level of prediction accuracy with a reduced order model while maintaining Li-ion concentration profile observability.
For uneven discretization 42, the state-space, ordinary differential equation (1), which describes the Li-ion concentration diffusion dynamics in the averaged solid particle, converts to the following matrix equation (6):




  
    
      
        
          
            
              c
              .
            
            s
          
          =
          
            
              Ac
              s
            
            +
            Bi
          
        
        ⁢
        
          

        
        ⁢
        where
        ⁢
        
          

        
        ⁢
        
          A
          =
          
            [
            
              
                
                  
                    -
                    
                      
                        2
                        ⁢
                        
                          D
                          
                            s
                            ,
                            k
                          
                        
                      
                      
                        r
                        1
                        2
                      
                    
                  
                
                
                  
                    
                      2
                      ⁢
                      
                        D
                        s
                      
                    
                    
                      r
                      1
                      2
                    
                  
                
                
                  ...
                
                
                  0
                
                
                  0
                
              
              
                
                  0
                
                
                  ⋱
                
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  0
                
              
              
                
                  ⋮
                
                
                  
                    
                      α
                      j
                    
                    ⁡
                    
                      (
                      
                        
                          1
                          
                            Δ
                            ⁢
                            
                              
                            
                            ⁢
                            
                              r
                              
                                j
                                -
                                1
                              
                            
                          
                        
                        -
                        
                          1
                          
                            r
                            j
                          
                        
                      
                      )
                    
                  
                
                
                  
                    -
                    
                      
                        α
                        j
                      
                      ⁡
                      
                        (
                        
                          
                            1
                            
                              Δ
                              ⁢
                              
                                
                              
                              ⁢
                              
                                r
                                j
                              
                            
                          
                          +
                          
                            1
                            
                              Δ
                              ⁢
                              
                                
                              
                              ⁢
                              
                                r
                                
                                  j
                                  -
                                  1
                                
                              
                            
                          
                        
                        )
                      
                    
                  
                
                
                  
                    α
                    ⁡
                    
                      (
                      
                        
                          1
                          
                            Δ
                            ⁢
                            
                              
                            
                            ⁢
                            
                              r
                              j
                            
                          
                        
                        +
                        
                          1
                          
                            r
                            j
                          
                        
                      
                      )
                    
                  
                
                
                  ⋮
                
              
              
                
                  0
                
                
                  
                    
                  
                
                
                  
                    
                  
                
                
                  ⋱
                
                
                  0
                
              
              
                
                  0
                
                
                  0
                
                
                  ...
                
                
                  
                    
                      α
                      
                        Mr
                        -
                        1
                      
                    
                    ⁡
                    
                      (
                      
                        
                          1
                          
                            Δ
                            ⁢
                            
                              
                            
                            ⁢
                            
                              r
                              
                                Mr
                                -
                                2
                              
                            
                          
                        
                        -
                        
                          1
                          
                            r
                            
                              Mr
                              -
                              1
                            
                          
                        
                      
                      )
                    
                  
                
                
                  
                    -
                    
                      
                        α
                        
                          Mr
                          -
                          1
                        
                      
                      ⁡
                      
                        (
                        
                          
                            1
                            
                              Δ
                              ⁢
                              
                                
                              
                              ⁢
                              
                                r
                                
                                  Mr
                                  -
                                  2
                                
                              
                            
                          
                          -
                          
                            1
                            
                              r
                              
                                Mr
                                -
                                1
                              
                            
                          
                        
                        )
                      
                    
                  
                
              
            
            ]
          
        
        ⁢
        
          

        
        ⁢
        
          B
          =
          
            
              [
              
                
                  
                    0
                  
                  
                    ...
                  
                  
                    0
                  
                  
                    
                      
                        -
                        
                          
                            α
                            
                              Mr
                              -
                              1
                            
                          
                          ⁡
                          
                            (
                            
                              1
                              +
                              
                                
                                  Δ
                                  ⁢
                                  
                                    
                                  
                                  ⁢
                                  
                                    r
                                    
                                      Mr
                                      -
                                      1
                                    
                                  
                                
                                
                                  r
                                  
                                    Mr
                                    -
                                    1
                                  
                                
                              
                            
                            )
                          
                        
                      
                      ⁢
                      
                        1
                        
                          
                            Fa
                            s
                          
                          ⁢
                          
                            D
                            s
                          
                        
                      
                    
                  
                
              
              ]
            
            T
          
        
      
    
    
      
        (
        6
        )
      
    
  




The state-space matrix equation (6) (ċs=Acs+Bi) is the continuous expression for the state-space equation around the equilibrium.
The discretized expression for the state-space equation around the equilibrium is represented by the following equation (7):

ck+1=Akck+Bkik  (7)

where
Ak=(I+AΔt)
Bk=BΔt
The battery terminal voltage equation is then given by the following equation (8):




  
    
      
        
          
            y
            k
          
          =
          
            
              
                H
                k
              
              ⁢
              
                c
                s
              
            
            -
            
              
                R
                0
              
              ⁢
              
                i
                k
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
        ⁢
        
          

        
        ⁢
        
          
            H
            k
          
          =
          
            
              
                ∂
                
                  y
                  k
                
              
              
                ∂
                
                  c
                  s
                
              
            
            =
            
              [
              
                
                  
                    0
                  
                  
                    ...
                  
                  
                    0
                  
                  
                    
                      
                        
                          
                            ∂
                            
                              U
                              ⁡
                              
                                (
                                
                                  c
                                  se
                                
                                )
                              
                            
                          
                          
                            ∂
                            
                              c
                              se
                            
                          
                        
                        
                      
                      k
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




and where, with the assumption that cse=cMr-1




  
    
      
        
          ∂
          
            U
            ⁡
            
              (
              
                c
                se
              
              )
            
          
        
        
          ∂
          
            c
            se
          
        
      
      
    
    k
  
  =
  
    
      
        
          ∂
          
            U
            ⁡
            
              (
              
                SOC
                se
              
              )
            
          
        
        
          ∂
          
            SOC
            se
          
        
      
      ⁢
      
        
          ∂
          
            SOC
            se
          
        
        
          ∂
          
            c
            
              se
              ⁢
              
                
              
            
          
        
      
    
    =
    
      
        1
        
          
            c
            
              ma
              ⁢
              
                
              
              ⁢
              x
            
          
          ⁡
          
            (
            
              
                θ
                
                  100
                  ⁢
                  %
                
              
              -
              
                θ
                
                  0
                  ⁢
                  %
                
              
            
            )
          
        
      
      ⁢
      
        
          d
          
            U
            ⁡
            
              (
              
                SOC
                se
              
              )
            
          
        
        
          d
          
            SOC
            se
          
        
      
    
  




In the output matrix Hk, a pre-calculated




  
    d
    
      U
      ⁡
      
        (
        
          SOC
          se
        
        )
      
    
  
  
    d
    
      SOC
      se
    
  




profile with respect to the estimated Li-ion concentration at the solid-electrolyte interface is used.

FIG. 4 illustrates a graph 80 including a plot line 82 of open circuit voltage of battery 12 versus the battery SOC. As shown in graph 80, the pre-calculated




  
    d
    
      U
      ⁡
      
        (
        
          SOC
          se
        
        )
      
    
  
  
    d
    
      SOC
      se
    
  




profile 84 corresponds to the slope 86 of plot line 82 at any given point along plot line 82.

From Ak and Hk (or from A and H) battery 12 is fully observable. (Note that the symbol A or Ak may be replaced herein with the symbol F or Fk.)
The model based Extended Kalman Filter (EKF) design follows the following series of equations.
Predict:

{circumflex over (x)}k|k−1=f({circumflex over (x)}k−1|k−1,uk−1)

Pk|k−1=Fk−1Pk−1|k−1Fk−1T+Qk 

Update:

{tilde over (y)}k=yk−h({circumflex over (x)}k|k−1)

Sk=HkPk|k−1HkT+Rk 

Kk=Pk|k−1HkTSk−1 

{circumflex over (x)}k|k={circumflex over (x)}k|k−1+Kk{tilde over (y)}k 

Pk|k=(I−KkHk)Pk|k−1 

where:



  
    
      
        
          
            
              
                F
                
                  k
                  -
                  1
                
              
              =
              
                
                  ∂
                  f
                
                
                  ∂
                  x
                
              
            
            
          
          
            
              
                x
                ^
              
              
                
                  k
                  -
                  1
                
                |
                
                  k
                  -
                  1
                
              
            
            ,
            
              u
              
                k
                -
                1
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            H
            k
          
          =
          
            
              ∂
              h
            
            
              ∂
              x
            
          
        
      
      
    
    
      
        x
        ^
      
      
        k
        |
        
          k
          -
          1
        
      
    
  




  
    
      P
      
        k
        |
        k
      
    
    =
    
      
        cov
        ⁡
        
          (
          
            
              x
              k
            
            -
            
              
                x
                ^
              
              
                k
                |
                k
              
            
          
          )
        
      
      =
      
        E
        ⁡
        
          (
          
            
              (
              
                
                  x
                  k
                
                -
                
                  
                    x
                    ^
                  
                  
                    k
                    |
                    k
                  
                
              
              )
            
            ⁢
            
              
                (
                
                  
                    x
                    k
                  
                  -
                  
                    
                      x
                      ^
                    
                    
                      k
                      |
                      k
                    
                  
                
                )
              
              T
            
          
          )
        
      
    
  




  
    
      P
      
        k
        |
        
          k
          -
          1
        
      
    
    =
    
      
        cov
        ⁡
        
          (
          
            
              x
              k
            
            -
            
              
                x
                ^
              
              
                k
                |
                
                  k
                  -
                  1
                
              
            
          
          )
        
      
      =
      
        E
        ⁡
        
          (
          
            
              (
              
                
                  x
                  k
                
                -
                
                  
                    x
                    ^
                  
                  
                    k
                    |
                    
                      k
                      -
                      1
                    
                  
                
              
              )
            
            ⁢
            
              
                (
                
                  
                    x
                    k
                  
                  -
                  
                    
                      x
                      ^
                    
                    
                      k
                      |
                      
                        k
                        -
                        1
                      
                    
                  
                
                )
              
              T
            
          
          )
        
      
    
  




  
    
      S
      k
    
    =
    
      
        cov
        ⁡
        
          (
          
            
              y
              k
            
            -
            
              
                y
                ^
              
              k
            
          
          )
        
      
      =
      
        E
        ⁡
        
          (
          
            
              (
              
                
                  y
                  k
                
                -
                
                  
                    y
                    ^
                  
                  k
                
              
              )
            
            ⁢
            
              
                (
                
                  
                    y
                    k
                  
                  -
                  
                    
                      y
                      ^
                    
                    k
                  
                
                )
              
              T
            
          
          )
        
      
    
  




  
    
      Q
      k
    
    =
    
      
        cov
        ⁡
        
          (
          
            ω
            k
          
          )
        
      
      =
      
        E
        ⁡
        
          (
          
            
              ω
              k
            
            ⁢
            
              ω
              k
              T
            
          
          )
        
      
    
  




  
    
      R
      k
    
    =
    
      
        cov
        ⁡
        
          (
          
            v
            k
          
          )
        
      
      =
      
        E
        ⁡
        
          (
          
            
              v
              k
            
            ⁢
            
              v
              k
              T
            
          
          )
        
      
    
  



The overall SOC estimation of battery 12 is summarized as follows.
Using the constructed state-space models and EKF, the Li-ion concentration profiles of battery 12 are estimated using according to the following equations:

ck+1=Fkck+Bkik 

yk=Hkcs−R0ik 

The estimated Li-ion concentration profiles being represented by the following equation:

ĉk+1=ĉk+Kk(yk−ŷk)

The SOC of battery 12 is then estimated from the Li-ion concentration profile information according to the following equation:



  
    
      S
      ⁢
      
        O
        ^
      
      ⁢
      C
    
    =
    
      
        
          
            θ
            _
          
          s
        
        -
        
          θ
          
            0
            ⁢
            %
          
        
      
      
        
          θ
          
            100
            ⁢
            %
          
        
        -
        
          θ
          
            0
            ⁢
            %
          
        
      
    
  




  where




  
    
      
        θ
        _
      
      s
    
    =
    
      
        
          c
          _
        
        s
      
      
        c
        
          s
          ,
          
            ma
            ⁢
            
              
            
            ⁢
            x
          
        
      
    
  




  
    
      
        c
        _
      
      s
    
    =
    
      
        
          
            c
            
              s
              ,
              1
            
          
          ⁢
          
            r
            1
            3
          
        
        +
        
          
            ∑
            
              i
              =
              1
            
            
              Mr
              -
              1
            
          
          ⁢
          
            
              3
              8
            
            ⁢
            
              (
              
                
                  c
                  
                    s
                    ,
                    i
                  
                
                +
                
                  c
                  
                    s
                    ,
                    
                      i
                      +
                      1
                    
                  
                
              
              )
            
            ⁢
            
              
                (
                
                  
                    r
                    i
                  
                  +
                  
                    r
                    
                      i
                      +
                      1
                    
                  
                
                )
              
              2
            
            ⁢
            
              (
              
                
                  r
                  i
                
                -
                
                  r
                  
                    i
                    +
                    1
                  
                
              
              )
            
          
        
      
      
        r
        
          Mr
          -
          1
        
        3
      
    
  



Referring now to FIG. 5, a block diagram 50 illustrating the operation of estimating the battery SOC using reduced-order battery model 30 in accordance with embodiments of the present invention is shown. Again, controller 16 is configured to implement the operation for estimating the battery SOC using reduced-order battery model 30. The operation employs a closed-loop estimator 52 according to reduced-order battery model 30. The design of estimator 52 pursuant to the Extended Kalman Filter (EKF) techniques described herein is characterized by the diffusion dynamics equation ck+1=Fkck+Bkik and the measured battery voltage equation yk=Hkcs−R0ik. As indicated above, the control input to estimator 52 is the current ik−1 54 of battery 12. Estimator 52 generates an estimated battery voltage vk 56 based on the battery current ik−1 54. The actual battery voltage vk 58 is measured. The estimated and measured battery voltages 56 and 58 are input to a comparator 60. Comparator 60 generates a difference Δv 62 indicative of the voltage difference between the estimated and measured battery voltages 56 and 58. The voltage difference Δv 62 is input back to estimator 52, which in turn generates another battery voltage estimate 56. This process repeats to cause the voltage difference Δv 62 to be close to zero. As a result the voltage difference Δv 62 being close to zero, estimator 52 is able to generate a relatively accurate estimated (Li-ion) concentration profile cs,k 64. The estimated concentration profile 64 is averaged by a calculator 66 as described herein to produce averaged concentration cs 68. A battery SOC calculator 70 estimates the battery SOC 72 based on the averaged concentration 68 as described herein.
A method and system for estimating the SOC of a battery from the measured battery voltage using a reduced-order battery model in accordance with embodiments of the present invention includes the following features. The electrochemical process of the battery is captured by reduced-order battery model 30. This means that once model 30 is calibrated, it can be used over a wide range of operating points, SOC, and temperature. The proposed reduced-order battery model 30 is much simpler than typical electrochemical battery models such as full-order battery model 20. This is possible by introducing several assumptions as described herein. Battery model 30 can be tuned by a relative small number of model parameters. The Extended Kalman Filter (EKF) is used for the SOC estimation and its formulation is simple enough to be implemented in the battery management system (BMS) of controller 16.
The method and system for estimating the SOC of a battery in accordance with embodiments of the present invention include other features. The method and system provide an alternative way to estimate the battery SOC with sufficient accuracy only using battery current input and the terminal voltage output information in real time. Reduced-order battery model 30 is relatively simple enough to be used for control design. The model calibration is much easier than the case of equivalent circuit battery model when the prediction accuracy is same. The SOC estimation accuracy is much less sensitive to the current sensor due to the inherent noise rejection capability of the EKF.
Referring now to FIG. 6A, a graph 90 depicting versus time the actual battery SOC and the battery SOC calculated using the reduced-order battery model in accordance with embodiments of the present invention is shown. Graph 90 includes a plot line 92 of the actual battery SOC. Graph 90 further includes a plot line 94 of the battery SOC calculated as described herein using the reduced-order battery model. As shown in graph 90, plot lines 92 and 94 essentially lie on top of one another thereby illustrating that the calculated battery SOC is highly accurate to the actual battery SOC.
FIG. 6B illustrates a graph 100 depicting the Li-ion concentration profile versus the normalized electrode solid particle radius at a certain time instant (namely, with reference to FIG. 6A, 200 seconds at which the battery SOC is about 0.9). Graph 100 includes a plot line 102 of the average concentration. Graph 100 further includes a plot line 104 of the concentration estimated as described herein. As such, graph 100 illustrates a comparison of the average concentration and the concentration estimated in accordance with embodiments of the present invention at a given time instant associated with a given battery SOC value.
FIG. 6C illustrates a graph 110 depicting the Li-ion concentration profile versus the normalized electrode solid particle radius at a certain time instant (namely, with reference to FIG. 6A, 1000 seconds at which the battery SOC is about 0.2). Graph 110 includes plot line 112 of the average concentration and plot line 114 of the estimated concentration estimated. As such, graph 110 illustrates a comparison of the average concentration and the concentration estimated in accordance with embodiments of the present invention at another time instant associated with another battery SOC value.
While exemplary embodiments are described above, it is not intended that these embodiments describe all possible forms of the present invention. Rather, the words used in the specification are words of description rather than limitation, and it is understood that various changes may be made without departing from the spirit and scope of the present invention. Additionally, the features of various implementing embodiments may be combined to form further embodiments of the present invention.

## Claims

1. A method for a powertrain having a Lithium-ion (Li-ion) traction battery, including solid active particles, and a traction motor, the method comprising:
providing electric power from the traction battery to the traction motor for propelling a vehicle; and
operating, by a controller in communication with the traction battery, the traction battery in providing electric power to the traction motor according to a state of charge of the battery based on an estimated Li-ion concentration profile of a representative solid active particle of a reduced-order model of the battery, the concentration profile being estimated according to the model from measured voltage and current of the battery.

2. The method of claim 1 wherein:
the reduced-order model employs an ordinary differential equation representing combined electrochemical dynamics of a positive electrode and a negative electrode of a battery cell of the battery.

3. The method of claim 2 further comprising:
deriving the reduced-order model from a full-order electrochemical model of the battery using uneven discretization of the ordinary differential equation.

4. The method of claim 1 further comprising:
computing the state of charge of the battery from the estimated concentration profile.

5. The method of claim 1 further comprising:
outputting a battery terminal voltage using the reduced-order model.

6. The method of claim 1 further comprising:
estimating the concentration profile using the reduced-order model.

7. The method of claim 5 further comprising:
estimating the concentration profile according to the reduced-order model from the measured voltage and current of the battery using Extended Kalman Filter (EKF) techniques.

8. A vehicle comprising:
a traction motor;
a Lithium-ion (Li-ion) traction battery including solid active particles; and
a controller to operate the battery to provide electricity to the motor according to a state of charge of the battery based on an estimated Li-ion concentration profile of a representative solid active particle of a reduced-order model of the battery and to estimate the concentration profile according to the model from measured voltage and current of the battery.

9. The vehicle of claim 8 wherein:
the controller is further configured to estimate the concentration profile using an Extended Kalman Filter (EKF).

10. The vehicle of claim 8 wherein:
the controller is further configured to output battery terminal voltage estimated according to the reduced-order model.