Modeling, Simulation, and Analysis of Lithium-Ion Batteries for Grid-Scale Applications

Lithium-ion batteries have become universally present in daily life, being used across a wide range of portable consumer electronics. These batteries are advantageous compared to other forms of energy storage due to their high energy density and long cycle life. These characteristics make lithium-ion batteries advantageous for many new and developing applications that require large scale energy storage such as electric vehicles and the utility grid. Typical uses for lithium-ion batteries require consistent cycling patterns that are predictable and easy to approximate across all uses, but new large scale applications will have much more dynamic demands. The cycling patterns for electric vehicles will vary based on each individuals driving patterns and batteries used for energy storage in the grid must be flexible enough to account for continuous fluctuations in demand and generation with little advanced notice. Along with these requirements, large scale applications do not want to sacrifice on cycle life and need to know that adding batteries will make operational and economic sense in specific cases. It is not possible to experimentally validate every possible driving pattern or grid storage need because of the great expense of these large systems and the long timescale required for testing. Therefore modeling of these systems is advantageous to help study specific application constraints and understand how lithium-ion batteries operate under those constraints. A systems level model is developed to study lithium-ion battery systems for use with solar energy (in a solar-battery hybrid system) and electric vehicles. Electrochemical based battery models are used as a component within larger systems. To facilitate fast simulation a single step perturbation and switch method is outlined for increasing the speed and robustness of solving the systems of DAEs that result from the systems level model. Operational characteristics are studied for lithium-ion batteries used to store solar energy within the electric grid. Different grid demands are tested against the system model to better understand the best uses for the solar-battery hybrid system. Both generic site studies and site specific studies were conducted. Solar irradiance data from 2010-2014 was obtained from 10 US based sites and used as an input to the system model to understand how the same system will operate differently at various locations. Technological benefits such as system autonomy were simulated for each site as well as economic benefits based on a time-of-use pricing scenario. These models included the growth of the solid-electrolyte interface layer on the battery electrodes to measure capacity fade during operation. This capacity fade mechanism allowed tracking of the site specific effects on battery life. A systems level model for an electric vehicle was also developed to simulate the growth of the SEI layer caused from different types of driving cycles and charging patterns. Results from both system models are presented along with an optimization method for the solar-battery hybrid model. In addition to modeling, experimental tests of LiFePO4 lithium-ion battery cells were conducted to measure capacity fade associated with different types of cycling throughout a batterys life. Cycling protocols were tested to study traditional capacity fade and also to focus on increasing a cells lifetime benefit through application switching

Global optimization of hybrid systems

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 2006.Includes bibliographical references (p. 339-353).Systems that exhibit both discrete state and continuous state dynamics are called hybrid systems. In most nontrivial cases, these two aspects of system behavior interact to such a significant extent that they cannot be decoupled effectively by any kind of abstraction and must be analyzed simultaneously. Hybrid system models are important in many areas of science and engineering, including flip-flops and latching relays, manufacturing systems, air-traffic management systems, controller synthesis, switched systems, chemical process systems, signaling and decision making mechanisms in (biological) cells, robotic systems, safety interlock systems, and embedded systems. The primary focus of this thesis is to explore deterministic methods for the global optimization of hybrid systems. While the areas of modeling, simulation and sensitivity analysis of hybrid systems have received much attention, there are many challenging difficulties associated with the optimization of such systems. The contents of this thesis represent the first steps toward deterministic global optimization of hybrid systems in the continuous time domain. There are various reasons for wanting to solve optimization problems globally.(cont.) In particular, there are many applications which demand that the global solution be found, for example, formal safety verification problems and parameter estimation problems. In the former case, a suboptimal local solution could falsely indicate that all safety specifications are met, leading to disastrous consequences if, in actuality, a global solution exists which provides a counter example that violates some safety specification. In the latter case, a suboptimal local solution could falsely indicate that a proposed model structure did not match experimental data in a statistically significant manner, leading to the false rejection of a valid model structure. In addition, for many optimization problems in engineering, the presence of nonlinear equality constraints makes the optimization problem nonconvex such that local optimization methods can often fail to produce a single feasible point, even though the problem is indeed feasible. The control parameterization framework is employed for the solution of the optimization problem with continuous time hybrid systems embedded. A major difficulty of such a framework lies in the fact that the mode sequence of the embedded hybrid system changes in the space of the optimization decision variables for most nontrivial problems.(cont.) This makes the resulting optimization problem nonsmooth because the parametric sensitivities of the hybrid system do not exist everywhere, thus invalidating efficient gradient based optimization solvers. In this thesis, the general optimization problem is decomposed into three subproblems, and tackled individually: (a) when the mode sequence is fixed, and the transition times are fixed; (b) when the mode sequence is allowed to vary, and the transition times are fixed; and (c) when the mode sequence is fixed, and the transition times are allowed to vary. Because even these subproblems are nontrivial to solve, this thesis focuses on hybrid systems with linear time varying ordinary differential equations describing the continuous dynamics, and proposes various methods to exploit the linear structure. However, in the course of solving the last subproblem, a convexity theory for general, nonlinear hybrid systems is developed, which can be easily extended for general, nonlinear hybrid systems. Subproblem (a) is the easiest to solve. A convexity theory is presented that allows convex relaxations of general, nonconvex Bolza type functions to be constructed for the optimization problem. This allows a deterministic branch-and-bound framework to be employed for the global optimization of the subproblem.(cont.) Subproblems (b) and (c) are much more difficult to solve, and require the exploitation of structure. For subproblem (b), a hybrid superstructure is proposed that enables the linear structure to be retained. A branch-and-cut framework with a novel dynamic bounds tightening heuristic is proposed, and it is shown that the generation of cuts from dynamic bounds tightening can have a dramatic impact on the solution of the problem. For subproblem (c), a time transformation is employed to transform the problem into one with fixed transition times, but nonlinear dynamics. A convexity theory is developed for constructing convex relaxations of general, nonconvex Bolza type functions with the nonlinear hybrid system embedded, along with the development of bounding methods, based on the theory of differential inequalities. A novel bounding technique that exploits the time transformation is also introduced, which can provide much tighter bounds than that furnished utilizing differential inequalities.by Cha Kun Lee.Ph.D