The field to which the disclosure generally relates to includes electrochemical device estimators and methods of making and using the same.
Vehicles having an electrochemical device may use a electrochemical device property estimator.
A number of illustrative variations may include a product comprising: a electrochemical device state estimator constructed and arranged to obtain terminal voltage data of a electrochemical device; determine an equivalent circuit model which operates in a manner approximating the terminal voltage data comprising an open circuit voltage VOC comprising a constant voltage source V0 and a voltage Vs across a capacitor Cs, in series with an overpotential circuit comprising an overpotential voltage V1; determining at least one of the power capabilities of the equivalent circuit model; and, estimating at least one of the power capabilities of the electrochemical device based upon the determined power capabilities of the equivalent circuit model.
A number of illustrative variations may include a product comprising: a electrochemical device state estimator constructed and arranged to obtain terminal voltage data of a electrochemical device; determine an equivalent circuit model which operates in a manner approximating the terminal voltage data comprising an open circuit voltage VOC comprising a constant voltage source V0 and a voltage Vs across a capacitor Cs, in series with an overpotential circuit comprising an overpotential voltage V1; determining at least one of the power capabilities of the equivalent circuit model; and, estimating at least one of the power capabilities of the electrochemical device based upon the determined power capabilities of the equivalent circuit model.
Other illustrative variations within the scope of the invention will become apparent from the detailed description provided hereinafter. It should be understood that the detailed description and specific examples, while disclosing variations within the scope of the invention, are intended for purposes of illustration only and are not intended to limit the scope of the invention.
Select examples of variations within the scope of the invention will become more fully understood from the detailed description and the accompanying drawings, wherein:
The following description of the variations is merely illustrative in nature and is in no way intended to limit the scope of the invention, its application, or uses.
Techniques and technologies may be described herein in terms of functional and/or logical block components and various processing steps. It should be appreciated that such block components may be realized by any number of hardware, software, and/or firmware components configured to perform the specified functions. For example, an embodiment of a system or a component may employ various integrated circuit components, e.g., memory elements, digital signal processing elements, logic elements, look-up tables, or the like, which may carry out a variety of functions under the control of one or more microprocessors or other control devices. In addition, those skilled in the art will appreciate that variations may be practiced in conjunction with any number of data transmission protocols.
For the sake of brevity, conventional techniques related to hybrid and electric vehicle operation, electrochemical device operation, electrochemical device diagnostics, vehicle computing modules, and other functional aspects of the systems (and the individual operating components of the systems) may not be described in detail herein. Furthermore, the connecting lines shown in the various figures contained herein are intended to represent example functional relationships and/or physical couplings between the various elements. It should be noted that many alternative or additional functional relationships or physical connections may be present in a number of variations of the subject matter.
As used herein, a “node” means any internal or external reference point, connection point, junction, signal line, conductive element, or the like, at which a given signal, logic level, voltage, data pattern, current, or quantity is present. Furthermore, two or more nodes may be realized by one physical element (and two or more signals can be multiplexed, modulated, or otherwise distinguished even though received or output at a common mode). The description herein refers to elements or nodes or features being “connected” or “coupled” together. As used herein, unless expressly stated otherwise, “connected” means that one element/node/feature is directly joined to (or directly communicates with) another element/node/feature, and not necessarily mechanically. Likewise, unless expressly stated otherwise, “coupled” means that one element/node/feature is directly or indirectly joined to (or directly or indirectly communicates with) another element/node/feature, and not necessarily mechanically.
In a number of variations, the controller 16 may comprise any type of control module or vehicle controller known in the art, and can be equipped with nonvolatile memory, random access memory (RAM), discrete and analog input/output (I/O), a central processing unit, and/or communications interfaces for networking within an automotive communications network. In a number of variations the controller 16 and possibly other illustrative blocks, modules, processing logic, and circuits described in connection with the embodiments disclosed herein may be implemented or performed with a general purpose processor, a content addressable memory, a digital signal processor, an application specific integrated circuit, a field programmable gate array, any suitable programmable logic device, discrete gate or transistor logic, discrete hardware components, or any combination thereof, designed to perform the functions described herein. A processor may be realized as a microprocessor, a controller, a microcontroller, or a state machine. A processor may also be implemented as a combination of computing devices, e.g., a combination of a digital signal processor and a microprocessor, a plurality of microprocessors, one or more microprocessors in conjunction with a digital signal processor core, or any other such configuration. Moreover, the steps of a method or algorithm described in connection with the variations disclosed herein may be embodied directly in hardware, in firmware, in a software module executed by a processor, or in any practical combination thereof. In a number of variations, a software module may reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. In this regard, an exemplary storage medium may be coupled to a processor such that the processor can read information from, and write information to, the storage medium. In a number of variations, the storage medium may be integral to the processor.
In a number of variations, the subject matter and certain aspects of the embodiments thereof may be described in the general context of computer-executable instructions, such as program modules, executed by one or more computing modules, controllers, or other devices. In a number of variations, program modules may include routines, programs, objects, components, data structures, and/or other elements that perform particular tasks or implement particular abstract data types. In a number of variations, the functionality of the program modules may be combined or distributed as desired in various variations. In a number of variations, a computing device that performs computer-executable instructions may include at least some form of computer readable media. Computer readable media may be any available media that can be accessed by a computing device and/or by applications executed by a computing device.
In a number of illustrative variations the product 10 may include any number of sensors. In a number of variations, the sensors may be may be capable of detecting one or more conditions which may include but are not limited to sound, pressure, temperature, acceleration, state of electrochemical device charge, state of electrochemical device power, current, voltage or magnetism and may be capable of producing at least one of sensor data or sensor signals and may sense and be at least one of polled or read by a control system or controller 16. In such variations the control system or controller 16 and any number of sensors may be used to obtain impedance, resistance, current, or terminal voltage data from the electrochemical device 12 at a number of temperatures and states of charge. Based at least upon obtained voltage data, an equivalent circuit model, such as a finite impulse response model circuit (FIR), R+N(R∥C) model circuit or R∥(R+C)N model circuit, or may be another type, may be determined, which may operate in a manner approximating the obtained electrochemical device impedance, resistance, current, terminal voltage data, or other condition of the electrochemical device 12. In a number of variations, the equivalent circuit model may be constructed by first determining a relation of electrochemical device 12 current to equivalent circuit model voltage over a period of time, and solving for a necessary number and value of each equivalent circuit component in adherence with a predetermined current voltage relation. In a number of variations, the equivalent circuit model may be constructed by first determining a relation of electrochemical device 12 current to equivalent circuit model current over a period of time, and solving for a necessary number and value of each equivalent circuit component in adherence with a predetermined current relation. In a number of variations, the equivalent circuit model may be constructed by first determining a relation of electrochemical device 12 voltage to equivalent circuit model current over a period of time, and solving for a necessary number and value of each equivalent circuit component in adherence with a predetermined current voltage relation. In such variations, the control system or controller 16 may be used to determine at least one of the power capabilities of the equivalent circuit model by use of mathematical algorithms including, but not limited to, differential equations, domain matrix exponentials, Laplace transform(s), Fourier transform(s), Fourier series, or any method of integrating a system of ordinary differential equations. Lastly, the control system may be used to determine at least one of the power capabilities of the electrochemical device 12 based upon at least one of the determined power capabilities of the equivalent circuit model.
In a number of variations, the SOC of the electrochemical device 12 may be dynamically tracked to determine when to charge the electrochemical device 12. The controller 16 may be suitably configured to control the SOC between 0% to 100%. Accurate knowledge of internal electrochemical device parameters may be important for purposes of controlling the electrochemical device systems of hybrid electric (and purely electric) vehicles. The methodology described herein adaptively extracts at least one internal parameters for electrochemical devices 12. The advantages of the proposed methodology over existing techniques are better adaptation to variations in the environment, electrochemical device state, and driving conditions, higher computational efficiency, and lower implementation cost including but not limited to, more accurate power prediction models for better power management in the electrochemical device 12 and the product 10 through the controller 16. In a number of variations, this new methodology may provide better power prediction for better power management in an electrochemical device 12 which may benefit range, performance and life of the electrochemical device 12.
To better control the electrochemical devices 12, controllers 16 may determine and process electrochemical devices 12 parameters such as, but not limited to, open-circuit voltage (VOC), ohmic resistance, terminal voltage, terminal current, the electrochemical device capacitance, or may be another condition or parameter. For non-limiting example, the VOC may be used to estimate the electrochemical device state of charge (SOC), which is an index associated with electrochemical device condition. In a number of variations, electrochemical device 12 terminal voltage and current may be used to estimate VOC and other electrochemical device internal parameters that may not be directly measurable during operation.
In a number of variations, a controller 16 (see
In a number of illustrative variations, voltage data for any value of voltage in the methodology discussed herein may be processed in real time using any number of linear regression methods which may include but are not limited to the use of a Kalman filter, weighted recursive least squares (WRLS) analysis, or any other method known in the art. In such variations, the equivalent circuit model may be constructed to operate in a manner approximating the processed data.
In a number of illustrative variations, the electrochemical device 12 yielding the processed data upon which the equivalent circuit model is based may have voltage and current limits. In such variations, for the sake of avoiding damage to the electrochemical device 12, power predictions for the electrochemical device 12 may be made by holding the equivalent circuit model current at an extreme constant value and determining whether the resultant circuit voltage will remain within the voltage limits of the electrochemical device 12. If it is determined that the circuit voltage will remain within the voltage limits of the electrochemical device 12, a current-limited power may be predicted for the electrochemical device 12 based on the extreme constant current. If, it is determined that the circuit voltage will not remain within the voltage limits of the electrochemical device 12, then the circuit voltage may be held at an extreme constant value within the electrochemical device 12 voltage limits, and the current corresponding to the extreme constant voltage may then be determined. A voltage-limited power may then be predicted for the electrochemical device 12 based on the extreme circuit voltage.
In a number of illustrative variations, the current or voltage of the system may held at a constant extreme, and a mathematical algorithm comprising at least one of a numerical integration of a differential equation system, Fourier series, a Fourier transform, Laplace transform, a closed-form expression for a FIR circuit model, a closed-form expression for a R+N(R∥C) model circuit or R∥(R+C)N model may be used in conjunction with a predetermined current voltage relationship of the equivalent circuit model to solve for the electrochemical device 12 power at time t.
In a number of variations, the equivalent circuit model 200 also includes a series resistance 210 disposed between a third node 212 and a fourth node 218. The series resistance 210 is labeled R in
In a number of variations, the equivalent circuit model 200 may also include an electrochemical device voltage source or open circuit voltage source 220 (labeled V0) disposed between the fourth node 218 and the second node 204. In a number of variations, the equivalent circuit model 200 may also include a capacitor Cs disposed between the fourth node 218 and the second node 2014. In a number of variations, the capacitor Cs may have zero initial voltage and include its contribution to terminal voltage as current flows through the electrochemical device. In a number of variations, Cs may depend on SOC at the start of a pulse or may get its value from a look-up table for the VOC. An open circuit voltage 224 (labeled VOC) is defined between fourth node 218 and second node 204. Referring to
In a number of variations, the measurable or terminal voltage V in
V(t)=VOC+V1(t) where VOC=V0+Vs(t) (1)
Power(t)=I(t)(VOC+V1(t)) (2)
In a number of variations, the controller 16 may continually estimate VOC and the parameters of V1 from the observations of terminal voltage and current time history. In a number of variations, the open circuit voltage VOC can be modeled as a source of constant voltage V0 in series with a capacitor Cs in equation 3 below where Q(t)=∫t
V
OC(t)=V0(t0)+Q(t)/CS (3)
In a number of variations the voltage Vs across the capacitor CS may be modeled by the differential equation shown in equation 4.
In a number of variations, the voltage Vs across the capacitor CS may be dependent on temperature T. In a number of variations, the voltage Vs across the capacitor CS may be determined for a number of different temperatures and modeled according to equation 4 below based on the different values of T, which can be determined via look-up tables. In a
number of variations, the value of Cs derives from the slope of
at a current temperature T, however, this may change with varying temperature.
V
0(SOC(t))=V0(SOC0)+k(SOC(t)−SOC0) (5)
In equation 5,
where SOC is the state of charge of the electrochemical device 12 and
where QT is the total charge capacity of the electrochemical device 12. State of Charge can be defined according to equation 6 below where Qa, is the total charge on the anode, and QT is the total charge that the electrochemical device 12 accepts as the anode goes from empty to full. SOC may be regressed in real time based on look-up tables from measured parameters monitored by the controller 16.
where it is assumed that the anode is empty at time t0.
In a number of variations, power according to equation 2 may be estimated or determined in two cases: 1) where a constant extreme value of current I(t) could be found; or 2) where a constant extreme value of voltage V(t) could be found and a mathematical algorithm would solve for power as a function of time t along with finding overpotential voltage V1(t).
Under case 1), when current is constant I(t)=I0, then the equation for voltage across the capacitor Vs can be written in equation 7 below.
In a number of variations, under case 1), it is necessary to solve for V1(t) in equation 1 with the constant current value I0, and known initial conditions for V0, CS, and I0. This may be done via a mathematical algorithm. In a number of variations, the mathematical algorithm may include one of: (A) using a numerical method to approximately solve the governing differential equations to obtain the approximate time evolution of the overpotential V1(t), (B) using a Laplace transform of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t), (C) using a matrix exponential of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t), (D) using a closed form expression for a FIR circuit model of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t), (E) using a closed form expression for a circuit model of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t).
For cases (A), (B), and (C) these may be known techniques known in the art. In a number of variations, for case (A), the equivalent circuit model 200 can be linear or non-linear, and it may be solved using well-known numerical integration methods such as but not limited to the Runge-Kutta method, the Adams-Bashforth method, the Euler method or may be another method.
In a number of variations for cases (A), (B), and (C), the capacitor Cs may have zero initial voltage and may depend on the state of charge (SOC) at the start of the pulse. In a number of variations, this can be determined from a look-up table for the VOC curve. An example VOC curve vs. SOC is shown in
In a number of variations, for a non-constant Cs and/or a non-linear model for V1, one can calculate the power P using ordinary differential equation system as shown in equations 9
where V is a vector of state variables, f is a vector-valued function, g is a scalar valued function, and Φ represents the parameters, e.g., resistance and capacitance values, of the model. The exact form of the functions f and g depends on which equivalent circuit model employed.
In a number of variations, this model may be incorporated to show VOC in Equation (1). In a number of variations, Cs may be a function of temperature where the value of Cs may be found from the slope of the VOC curve interpolated from a look up table. In a number of variations, the look-up table may be used internally in the controller 16. In a number of variations, the look-up table may be measured experimentally in the electrochemical device 12 and stored for use in the controller 16. In a number of variations, the constant voltage V0 and the parameters of overpotential V1 may be regressed in real time using the controller 16. In a number of variations, V1 can be modeled with a method given the kinetics of the model. In a number of variations, this case will provide a showing of V1(t).
In a number of variations, the inclusion of the capacitor Cs may better predict power by accounting for changes in open-circuit voltage during a power pulse. This model of the open-circuit voltage may be combined with different modeling methods used in electrochemical device 12 controllers 16 for the overpotential V1 to determine the overall power given equation (1) including, but not limited to, finite impulse response (FIR)(pulse at 1, 2 or 10 seconds by non-limiting example), R+N(R∥C), R+C+N(R∥C), or may be another type. In a number of variations, in the case of a R+N(R∥C) circuit model, the new circuit model becomes a Cs+R+N(R∥C) circuit model, which is naturally in R+(N+1)(R∥C) form wherein the new Cs element can be considered as an R∥C pair with R=∞. In a number of variations, in a FIR estimate circuit, the power of charge and discharge can be sustained by assuming VOC stays constant during the power pulse. For longer power pulses, enough charge may be passed in or out of the electrochemical device 12 so that VOC moves to a substantially new value during the pulse. In a number of variations, this movement is greater for operation in regions where the VOC curve is steeper, as determined by look up tables. In a number of variations, equation (3) above can better model the change in VOC value as the linear approximation of the VOC around the initial VOC value which shows linear relation to the capacitor Cs circuit element. In a number of variations, if the modeling method for overpotential V1 is linear then the overall (C+(FIR), C+R+N(R∥C), C+R∥(R+C)N, C+Any Linear RC circuit) including the capacitance is linear, and then the state of power can be calculated using closed-form formulas.
In a number of variations, when the ODE model of the overpotential is linear, as it will be for any equivalent circuit model that is a combination of linear resistors and capacitors, the compound circuit built by inserting the capacitor Cs in series with the overpotential circuit remains linear. Because it is linear, this compound circuit can be solved using (B) Laplace transform methods or (C) matrix exponentials. These techniques are well known to those schooled in the art, so we do not further elaborate here.
For case (D), in a number of variations, V1 can be modeled in a FIR model where V1 is calculated from equation 8 below.
where V0(SOC(t))=V0(SOC0)+k(SOC(t)−SOC(0)) where
and where SOC is the state of charge of the electrochemical device and
where Q is the total state of charge of the electrochemical device wherein the overpotential V1 is determined via equation 11 below:
V
1(t)=Σn=0LHn I0(t−nΔt) (11)
wherein the currents values I(t−nΔt) before t=0 are preceding measured values, the current values after t=0 are constant and equal to the above mentioned constant input current, Δt is the sampling period, and wherein {H(nΔt); n=0,1, . . . L} are the coefficients of the FIR filter, {H(nΔt)} can also be viewed as a digitized impulse response.
In a number of variations, in the case of an FIR model, the solution when current is held at a constant value takes a simple form which may be advantageous for implementation.
For case (E), in a number of illustrative variations, and as illustrated by
In a number of illustrative variations, for the purpose of determining the power capabilities of an R+N(R∥C) circuit, it may be assumed that input current, I is held constant at an extreme (allowable, insofar as the device 12 is not damaged) value, say I0, for a chosen interval, t seconds. In such variations, assuming N R∥C pairs, voltage across capacitor i, Ci at time, t may be solved for according to equation 12 below
and power at time, t may be predicted according to equation 13
In a number of variations, once one has an equation for V1(t) has been determined via the cases (A)-(E) above, power can be calculated from equation 14 below:
In a number of variations, for case (2), where a constant extreme value of voltage V(t) could be found and a mathematical algorithm would solve for power as a function of time t, then voltage is held constant at Vconst=V0+Vs+V1, however, neither Vs nor V1 is constant and current now varies with time I(t). Therefore, in a number of variations, the capacitor Cs is considered in series with the overpotential circuit 265 and modified into a modified overpotential circuit 267. The power calculation is now shown in equation 15 below:
P(t)=VconstI(t) (15)
In a number of variations, it is necessary to solve for I(t) in equation (15) with the constant current value Vconst. This may be done via a mathematical algorithm. In a number of variations, the mathematical algorithm may include one of: (A) using a numerical integration of an ordinary differential equation system for the combined circuit where it may be linear or non-linear for the time evolution of the current I(t), (B) using a Laplace transform of the circuit voltage to formulate an equation for the time evolution of the current I(t), (C) using a matrix exponential of the circuit voltage to formulate an equation for the time evolution of the current I(t), or (D) using a closed form expression for a FIR circuit model of the circuit voltage to formulate an equation for the time evolution of the overpotential V1(t).
For cases (A), (B), and (C) these may be known techniques known in the art. In a number of variations, for case (A), the equivalent circuit model 200 can be linear or non-linear.
In a number of variations, one can calculate the power P using ordinary differential equation system, which in a general form is as shown in equations 16 below.
where V is a vector of state variables, f is a vector-valued function, g is a scalar valued function, and t represents the parameters, e.g., resistance and capacitance values, of the model. The exact form of the functions f and g depends on which equivalent circuit model employed. To solve system (7) in the most general case, one recognizes that Vconst=g(V, I) is an algebraic equation and therefore known numerical techniques for solving differential-algebraic equations (that is, a DAE solver) can be employed.
In a less general case, I appears in g(V,I) in such a way that the relationship can be inverted to find I=h(Vconst,V). That is, h(Vconst,V) is a function such that Vconst=g(V, h(Vconst, V)) holds true identically for all possible values of Vconst and the vector V. In such a case, I can be eliminated from the DAE system (7) to obtain the ODE system with equation 17.
which can then be solved using a numerical method for ODE systems such as but not limited to the Runge-Kutta method, the Adams-Bashforth method, and the Euler method. Once V(t) has been so obtained, the power is evaluated as shown in equation 18 below.
P(t)=Vconsth(Vconst,V(t) (18)
Furthermore, when in (7) the functions f(V,I,T; Φ) and g(V,I) are linear in V and I, then general methods for solving linear ODE system, such as the Laplace transform technique or matrix exponentials, can be applied to evaluate V(t) and I(t), so that power can be evaluated as P(t)=VconstI(t). As these techniques are well known to those schooled in the art, we do not further elaborate here.
For case (D), in a number of variations, V1 can be modeled in a FIR model where V1 is calculated from equation 8 below.
where
and where SOC is the state of charge of the electrochemical device and
where Q is the total state of charge of the electrochemical device wherein the overpotential V1 is determined via equation 20 below
V
1(t)=Σn=0LHnI(t−nΔt) (20)
wherein the currents values I(t−nΔt) before t 32 0 are preceding measured values, the current values after t 32 0 are constant and equal to the above mentioned constant input current, Δt is the sampling period, and wherein {H(nΔt); n=0,1, . . . L} are the coefficients of the FIR filter, {H(nΔt)} can also be viewed as a digitized impulse response.
A number of variations may include a method including using a state of power predictor comprising an equivalent which may be modeled based on terminal voltage from an electrochemical device 12 such as but not limited to a electrochemical device, supercapacitor or other electrochemical device and processing the runtime values of that model circuit using any number of known real-time linear regression algorithms including, but not limited, to a weighted recursive least squares (WRLS), Kalman filter or other means. In a number of variations, the constant voltage V0 and the parameters of overpotential V1 may be regressed in real time using the controller 16 using these methods. The method may also include a controller 16 constructed and arranged to receive input from the state of power predictor, compare the input from the predictor with predetermined values and take action such as send a signal representative of the predicted state of power or take other action when the input from the predictor is within a predetermined range of the predetermined values. In a number of variations the controller 16 may be constructed and arranged to prevent a particular usage of an electrochemical device based upon the state of power prediction. In a number of variations, the controller may be constructed and arranged to use at least a segment of voltage-based electrochemical device power capability estimation data, and using real-time linear regression, which may be used as a method of estimating future behavior of a system based on current and previous data points, to provide a robust state of power predictor. Linear regression may be performed by forming an RC circuit which may be equivalent to electrochemical impedance spectroscopy data and processing the runtime values of that RC circuit using any number of known real-time linear regression algorithms which may include, but are not limited to, a weighted recursive least squares (WRLS) algorithm, Kalman filter algorithm or other means. In a number of variations the controller 16 and any number of sensors may be used to obtain impedance data from a electrochemical device at a number of electrochemical device temperatures and electrochemical device states of charge; building an equivalent R+N(R∥C) or R∥(R+C)N circuit which operates in a manner approximating the obtained voltage data; determining at least one of the power capabilities of the equivalent circuit model by use of domain matrix exponentials, a Laplace transform, a Fourier transform, a Fourier series, or any other method of integrating a system of ordinary differential equations; and, estimating at least one of the power capabilities of the electrochemical device based upon at least one of the determined power capabilities of the equivalent circuit model.
The following description of variants is only illustrative of components, elements, acts, products and methods considered to be within the scope of the invention and are not in any way intended to limit such scope by what is specifically disclosed or not expressly set forth. The components, elements, acts, products and methods as described herein may be combined and rearranged other than as expressly described herein and still are considered to be within the scope of the invention.
Variation 1 may include a method comprising: obtaining terminal voltage data of a electrochemical device; determining an equivalent circuit model which operates in a manner approximating the terminal voltage data comprising an open circuit voltage VOC comprising a constant voltage source V0 and a voltage Vs across a capacitor Cs, in series with an overpotential circuit comprising an overpotential voltage V1; determining at least one of the power capabilities of the equivalent circuit model; and, estimating at least one of the power capabilities of the electrochemical device based upon the determined power capabilities of the equivalent circuit model.
Variation 2 may include a method as set forth in variation 1 wherein the terminal voltage is obtained at a number of temperatures and modeled as a function of temperature within the equivalent circuit model.
Variation 3 may include a method as set forth in variations 1-2 wherein the value of the voltage Vs across the inserted capacitor Cs is a function of temperature and modeled as a function of temperature within the equivalent circuit model and where
where QT is the total capacity of the electrochemical device,
and SOC is the state of charge of the electrochemical device.
Variation 4 may include a method as set forth in variations 1 wherein the open circuit constant voltage V0 and overpotential V1 are regressed in real time.
Variation 5 may include a method as set forth in variations 1-4 wherein the overpotential circuit model is a linear circuit.
Variation 6 may include a method as set forth in variations 1-4 wherein the overpotential circuit is a non-linear circuit.
Variation 7 may include a method as set forth in variations 1-6 wherein the overpotential circuit comprises at least one capacitor.
Variation 8 may include a method as set forth in variations 1-7 wherein the overpotential circuit comprises at least one resistor.
Variation 9 may include a method as set forth in variations 1-8 wherein estimating at least one of the power capabilities of the electrochemical device based upon the power capabilities of the equivalent circuit comprises:
imposing a constant input current, I0 upon the equivalent circuit model;
solving for the overpotential voltage, V1(t) according to V(t)=V(0)+Vs+V1(t), where
and initial conditions V(0), and Cs are known, and solving for V1(t) through a mathematical algorithm;
predicting an equivalent circuit model power at time, t according to
and, correlating the equivalent circuit model power at time, t to the power of the electrochemical device at time t where V0(SOC(t))=V0(SOC0)+k(SOC(t)−SOC(0)) where
and where SOC is the state of charge of the electrochemical device.
Variation 10 may include a method as set forth in variations 1-8 wherein estimating at least one of the power capabilities of the electrochemical device based upon the power capabilities of the equivalent circuit comprises:
imposing a an extreme constant input voltage, V upon the equivalent circuit model such that V=(V0+Vs+V1);
solving for an equivalent circuit model current at time t by assuming a constant overpotential for the equivalent circuit model through a mathematical algorithm;
predicting for an equivalent circuit model power at time t via the equation:
Power(t)=VI(t); and,
correlating the equivalent circuit model power at time, t to the power of the electrochemical device at time t where V0(SOC(t))=V0(SOC0)+k(SOC(t)−SOC(0)) where
and where SOC is the state of charge of the electrochemical device.
Variation 11 may include a method as set forth in variation 9 wherein the mathematical algorithm comprises using a numerical integration of the ordinary differential equation of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t).
Variation 12 may include a method as set forth in variation 9 wherein the mathematical algorithm comprises using a Laplace transform of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t).
Variation 13 may include a method as set forth in variations 9 wherein the mathematical algorithm comprises using a matrix exponential of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t).
Variation 14 may include a method as set forth in variation 9 wherein the mathematical algorithm comprises using a close form expression for a FIR circuit model of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t).
Variation 15 may include a method as set forth in variation 9 wherein the mathematical algorithm comprises using a close form expression for a R+N(R∥C) circuit model of the overpotential circuit voltage to formulate an equation for the time evolution of the overpotential V1(t).
Variation 16 may include a method as set forth in variation 10 wherein the mathematical algorithm comprises using a numerical integration of the ordinary differential equation of the overpotential circuit current to formulate an equation for the time evolution of the current overpotential I(t).
Variation 17 may include a method as set forth in variation 10 wherein the mathematical algorithm comprises using a Laplace transform of the overpotential circuit current to formulate an equation for the time evolution of the overpotential current I(t).
Variation 18 may include a method as set forth in variation 10 wherein the mathematical algorithm comprises using a matrix exponential of the overpotential circuit current to formulate an equation for the time evolution of the overpotential I(t).
Variation 19 may include a method as set forth in variation 10 wherein the mathematical algorithm comprises using a close form expression for a FIR circuit model of the overpotential circuit current to formulate an equation for the time evolution of the overpotential I(t).
Variation 20 may include product comprising: a electrochemical device state estimator constructed and arranged to obtain terminal voltage data of a electrochemical device; determine an equivalent circuit model which operates in a manner approximating the terminal voltage data comprising an open circuit voltage VOC comprising a constant voltage source V0 and a voltage Vs across a capacitor Cs, in series with an overpotential circuit comprising an overpotential voltage V1; determining at least one of the power capabilities of the equivalent circuit model; and, estimating at least one of the power capabilities of the electrochemical device based upon the determined power capabilities of the equivalent circuit model.
The above description of select variations within the scope of the invention is merely illustrative in nature and, thus, variations or variants thereof are not to be regarded as a departure from the spirit and scope of the invention.