An Optimal Fast-Charging Strategy for Lithium-Ion Batteries via an Electrochemical–Thermal Model with Intercalation-Induced Stresses and Film Growth

Optimal fast charging is an important factor in battery management systems (BMS). Traditional charging strategies for lithium-ion batteries, such as the constant current–constant voltage (CC–CV) pattern, do not take capacity aging mechanisms into account, which are not only disadvantageous in the life-time usage of the batteries, but also unsafe. In this paper, we employ the dynamic optimization (DP) method to achieve the optimal charging current curve for a lithium-ion battery by introducing limits on the intercalation-induced stresses and the solid–liquid interface film growth based on an electrochemical–thermal model. Furthermore, the backstepping technique is utilized to control the temperature to avoid overheating. This paper concentrates on solving the issue of minimizing charging time in a given target State of Charge (SoC), while limiting the capacity loss caused by intercalation-induced stresses and film formation. The results indicate that the proposed optimal charging method in this paper offers a good compromise between the charging time and battery aging.


Introduction
Lithium-ion batteries have been used in many electronic products due to their high cell voltage, high energy density, high power density, convenient operating temperature range, lack of memory property, and high cycle life [1]. When operating a lithium-ion battery efficiently and safely during charging, long charging time, capacity degradation, capacity wastage, and overheating are the main difficulties that need to be overcome.
In recent years, many researchers have made efforts to optimize lithium battery charging. Many results are based on traditional charging patterns, such as the constant-current (CC) pattern and the constant current-constant voltage (CC-CV) pattern, without considering the aging process. For example, Liu and Luo [2] proposed a Taguchi-based algorithm and adopted orthogonal arrays to obtain the optimal rapid-charging strategy for a piecewise CC charging approach, which can charge a battery cell from 0% to 75% within 40 min and increase the cycle life by more than 60%. On the basis of the open circuit voltage (OCV)-resistance equivalent circuit model, Abdollahi et al. [3] presented a closed-form solution for optimally charging a lithium-ion battery; the target function is established through a combination of two consumption functions: time-to-charge (TTC) and energy losses (EL). Here, the CC-CV pattern is selected as the optimal charging scheme, where the current in the CC stage is a function of the ratio of weighting on TTC and EL. Then, Abdollahi et al. [4] proposed the objective function consisting of TTC, EL, and a temperature rise index (TRI). Then, the value of the current in the CC stage is also considered as a function of the ratio of weighting on TTC and EL, 2 of 16 and finally the analytical solution for the optimal problem is derived. Monem et al. [5] studied the influence of three charging strategies including CC, CC-CV, and constant current-constant voltage with negative pulse (CC-CVNP) on the battery's cycle life. The results show that the CC-CVNP pattern with low amplitude and less negative pulses is more efficient than the CC and CC-CV patterns. Liu et al. [6] firstly put forward a triple-objective function for optimal battery charging on the basis of a coupled thermoelectric model. Then, the CC-CV charge strategy is optimized, which offers the best compromise among three significant performance indexes consisting of charging time, energy loss, and temperature rise. Fang et al. [7] permit users to specify charging objectives and reach them by dynamic optimal control for the first time and proposed two charging methods based on the linear quadratic control theory without real-time constrained optimization computation. Compared to the conventional open-loop regulation of fast charging, the close-loop optimal method can be used to accurately control the specific parameters, such as temperature, current, and voltage. Patnaik et al. [8] came up with a constant temperature-constant voltage (CT-CV) charging algorithm that considers battery temperature as a key feedback variable. Then, a simple and easy-to-realize proportional-integral-derivative (PID) controller is employed to implement this close-loop method and the results indicate that the proposed approach achieved a 20% faster charging rate with an identical total temperature increase as compared to the constant current-constant voltage (CC-CV) technique. Klein et al. [9] paid attention to minimum charging time and proposed a simple one-step predictive control algorithm that is capable of solving the time-optimal solution and satisfying the real-time requirement. On the basis of the electrochemical battery model, Pramanik et al. [10] introduced a novel method for optimally charging the lithium-ion battery cell, which establishes the objective function aiming to minimize the charging cost. The result indicates that the optimal charging method presented in the paper [10] can decrease the charging time of a lithium-ion cell, meanwhile guaranteeing the temperature limit when compared with the traditional constant current charging. Considering the electrolyte and thermal dynamics based on a single particle model, the Legendre-Gauss-Radau (LGR) pseudo-spectral approach is used to solve the problem of nonlinear multistate optimal control, and the minimum time charge strategies are analyzed minutely while taking the solid and electrolyte phase concentration limits and temperature constraint into account [11].
Generally speaking, fast charging can accelerate the battery aging processes. For reducing the aging rate, some researchers consider the aging process when optimizing the fast-charging strategy and some good results have been obtained. For example, through coupling incremental capacity (IC) and IC peak area analysis with the mechanistic model, Ansean et al. [12] quantified the mechanism of degradation that leads to the aging of the battery cell. In addition, the results show that aging is caused by a loss of lithium ions and a lower level of loss of active material on the negative electrode. On the basis of cycle-life testing (up to 4500 cycles), Ansean et al. [13] proposed a multistage fast-charging algorithm which allows a full recharging of the cell (0% to 100% SoC) within 20 min (indeed after 4500 cycles are reached) and does not cause any remarkable degradation to the battery cell. Considering the influence of intercalation-induced stress on aging, Suthar et al. [14] used dynamic optimization to achieve the optimal current profile to fast charge a lithium-ion battery through a single-particle model while coupling this with the intercalation-induced stress generation model. In addition, this was the first time protocols for optimally charging batteries while ensuring a minimal mechanical cost to the electrode particles during intercalation were demonstrated. Torchio et al. [15] used the first-principles pseudo-two-dimensional (P2D) model together with the capacity fade mechanisms that work when the battery is operating. Then, the model predictive control (MPC) method based on a linearized model of the P2D model was proposed to approach a target value of the state of charge (SoC) while considering the degradation process of the system as well as the thermal and voltage limits. To analyze the effect of static and dynamic fast-charging current strategies on the degradation properties of lithium-ion batteries, Monem et al. [16] applied both the static and dynamic fast-charging current profiles to a lithium-ion battery cell. After 1700 cycles, the result shows that the dynamic fast-charging current profile had a more outstanding role in reducing the aging rate and the charging time of Energies 2020, 13, 2388 3 of 16 cells than the static fast-charging profile. To reduce the influence of fast charging on battery aging, Ali et al. [17] proposed a temperature control method based on fuzzy logic that protects the batteries from overvoltage and overheating. The result illustrates that the proposed fast-charging pattern spends 9.76% less time during full battery recharging than the conventional CC-CV method, and the approach does not bring significant degradation. As a compromise between the three objectives of high safety, longer lifetime, and a lower charging time, Zou et al. [18] proposed a fast-charging method on the basis of the electrochemical model and MPC theory. Here, the battery optimal charging issue is described in a linear time-varying model for the implementation of the MPC algorithm. Similarly, Lina et al. [19] presented an electrolyte enhanced single particle model with aging mechanisms which considers the effect of electrolyte dynamics, then the dynamic programming DP approach is adopted to obtain the optimal charging profiles to reduce the charge time and battery aging. On the basis of the electrochemical-thermal capacity fade model, Xu et al. [20] used the DP optimization algorithm to minimize capacity fade, temperature rise, and charging time. Although there are many papers considering the aging process in terms of optimizing the fast charge curve, the aging models used in the above-mentioned papers seldom involve intrinsic aging mechanisms and never consider the effects of intercalation-induced stress and film growth together on the aging process.
In this paper, the electrochemical-thermal model is employed to obtain the optimal charging profile and control the temperature of cell. Controlling the temperature of the hottest point inside the cell based on the backstepping method can help avoid overheating. During fast charging, intercalation-induced stress will cause particle fracture which can accelerate the aging process. Hence, restricting intercalation-induced stress to a given range is significant to reduce the degradation of cells when seeking the optimal fast-charging profile. Furthermore, the main aging cause is the growth of film on the surface of particles, and this film is a compound containing lithium which cannot be reused. Therefore, confining the growth rate of surficial film to a proper range can maximize the available lithium, which is another contribution of this paper. To sum up, on the basis of maintaining a constant temperature using the backstepping control method and minimizing the charging time while limiting intercalation-induced stress and the film growth to an appropriate range are the main contributions of this paper.

SPM-Electrolyte-T Model
The SPM-Electrolyte-T Model is a simplified version of the first principle electrochemical model of batteries composed of a single particle model, electrolyte dynamics, and thermal dynamics. Considering the effects of electrolyte and thermal dynamics can help us investigate the fast-charging scheme to reduce battery aging to prolong the lifetime of lithium-ion batteries. The single particle model can be described by the following partial differential equation [21]: with the following conditions: The definitions of the relevant electrochemical parameters are listed in the attached table.
Most studies relevant to the electrochemical model of lithium-ion batteries ignore the influence of electrolyte dynamics for inexpensive computation. However, when the charging or discharging current is high, there would be a remarkable electrolyte concentration potential between the cathode and the anode, which has a significant impact on the prediction of the terminal voltage. Specifically, the process of fast charging requires more than a 2 C current, which has enough power to produce an Energies 2020, 13, 2388 4 of 16 obvious concentration potential. Hence, during the operation of fast charging, it is wise to consider electrolyte dynamics.
The electrolyte dynamics are governed by where D e f f e (c e , T avg ) = D e (c e , T avg )(c e ) · (ε e ) brug is the effective electrolyte diffusivity, which means that the diffusion rate of lithium ions is related to its concentration and temperature; t 0 c is the transference number, which is deemed to be constant; ε + e , ε − e , ε sep e are the volume fractions of the electrolyte in the cathode, anode, and separator, respectively, and L + , L − , L sep are the lengths of the cathode, anode, and separator. The effects of electrolyte dynamics on the terminal voltage equation is presented as follows, considering: where Most model parameters of battery cells are relative to the temperature and are expressed by the following equation: where P 0 is the parameter when the temperature is equal to T 0 ; and E a P is the activation energy of parameter P. The average temperature T avg can be obtained by the following equations: with the boundary conditions Energies 2020, 13, 2388 where r ∈ [0, R b ] is the radial coordinate of a cylindrical battery; α = k/(ρC p ) is proportional to the average thermal conductivity k; ρ is the average mass density and C p is the average specific thermal capacity; . q is the volumetric heat generated rate; and T a is the ambient temperature.
. q can be computed by the following equation: where U o (t) is the open circuit voltage, which is a function of the boundary concentration; U ter = V(t) is the terminal voltage; and ∂U o (t)/∂T(t) is the entropic heat generation which is too small to neglect. The open circuit voltage U o (t) can be given by the following equation: Then, the volumetric heat generated rate . q can be rewritten as follows: .

Intercalation-Induced Stress
The intercalation-induced stress generated by lithium-ion intercalation and deintercalation influences the diffusion of lithium and even causes electrode particle fracture. This electrode particle fracture will accelerate aging. Considering the effects of intercalation-induced stress helps us to investigate the fast-charging scheme to reduce the electrode particle fracture and thus prolong the lifetime of the lithium-ion battery. The intercalation-induced stresses are composed of radial stress σ r and tangential stress σ t , which are dependent of the lithium concentration [22]: . Ω ± , E ± , and ν ± are the partial molar volume, Young's Modulus, and Poisson's ratio of the electrode material in the cathode and anode, respectively, which have different sensitivities to temperature. c ± is the concentration change from the stress-free value. The resultant stress is a weighted sum of σ r and σ t : From [23], we know that the maximum absolute values of radial and tangential stress are located at the center and the surface of the particle, respectively. Taking the anode as an example, we have Energies 2020, 13, 2388 6 of 16 According to Equation (19), we find that the maximum value of the resultant stress is located at the point where the lithium concentration is the minimum value of the electrode particle. That is Figure 1 shows the input current under Urban Dynamometer Driving Schedule (UDDS), which is applied to a battery cell. Then, the maximum absolute values of radial stress, tangential stress, and the resultant stresses are investigated. According to Equation (21), the intercalation-induced stresses are caused by the lithium-ion concentration nonuniformity in space, which is related to the values of the charging and discharging current. When the current is positive, the lithium-ion concentration in the surface of the particle is smaller than in the center of the particle, that is, the maximum radical stress is positive and the maximum tangential stress is negative, and vice versa. From the input current under UDDS, we know that the current is rapidly changing, which can cause the oscillation of lithium-ion concentration and enlarge the imbalance of lithium-ion concentration in space. Hence, the maximum absolute values of radial stress, tangential stress, and the resultant stresses change rapidly over time, as shown in Figures 2-4. From the three figures, we can see that when the input current is negative, the absolute values of the three stresses will increase gradually, and vice versa. This is because the negative current increases the lithium-ion concentration nonuniformity in space and the positive current decreases this nonuniformity. Besides that, we can find that the absolute values of the three stresses have the same changing curve, which agrees with Equations (20) and (21). Another finding is that the value of the maximum tangential stress is larger than the maximum resultant stress, and the maximum resultant stress is always larger than the maximum radial stress. Hence, by adjusting the input current, the radial stress, tangential stress, and the resultant stress can be controlled to decrease the aging of a battery cell.
Energies 2020, 13, 03 FOR PEER REVIEW 7 of 18 resultant stress, and the maximum resultant stress is always larger than the maximum radial stress. Hence, by adjusting the input current, the radial stress, tangential stress, and the resultant stress can be controlled to decrease the aging of a battery cell.   Energies 2020, 13, 03 FOR PEER REVIEW 7 of 18 resultant stress, and the maximum resultant stress is always larger than the maximum radial stress. Hence, by adjusting the input current, the radial stress, tangential stress, and the resultant stress can be controlled to decrease the aging of a battery cell.

Film Growth Model
Ramadass et al. [24] studied a resistive film formed on the anode electrode/electrolyte interphase, which is the main cause of capacity loss. The veracious chemical side reaction generating the resistive film lies on the chemical component of the anode electrode and electrolyte. Furthermore, a simplified and universal model for capacity reduction based on film growth originated from a single particle model is presented here: S Li e P + − + + → (22) where S is the type of solvent, and P is the reaction product.
On account of this nonreversible side reaction, the reaction product builds up a film at the interface of electrode/electrolyte with a thickness varying over time film ( ) t δ . The irreversibly formed membrane together with the initial solid interphase resistance SEI R forms the total resistance of the electrode/electrolyte interface as follows: where p κ is the conductivity of the formed film.
The key state related to the increasing rate of interface film thickness, because of the unwanted lithium loss, is described by where, n a , p M , p ρ , and F are the formed film's specific surface area, molecular weight, mass density, and Faraday's constant, respectively. The variable ( ) s J t is the local volume current density used for the side reaction and described by Butler-Volmer kinetics. Assuming that the lithium reduction reaction is nonreversible and the concentration of lithium ions in the solution changes little, it is possible to approximate ( ) s J t through the following equation:

Film Growth Model
Ramadass et al. [24] studied a resistive film formed on the anode electrode/electrolyte interphase, which is the main cause of capacity loss. The veracious chemical side reaction generating the resistive film lies on the chemical component of the anode electrode and electrolyte. Furthermore, a simplified and universal model for capacity reduction based on film growth originated from a single particle model is presented here: where S is the type of solvent, and P is the reaction product. On account of this nonreversible side reaction, the reaction product builds up a film at the interface of electrode/electrolyte with a thickness varying over time δ f ilm (t). The irreversibly formed membrane together with the initial solid interphase resistance R SEI forms the total resistance of the electrode/electrolyte interface as follows: where κ p is the conductivity of the formed film.
The key state related to the increasing rate of interface film thickness, because of the unwanted lithium loss, is described by where, a n , M p , ρ p , and F are the formed film's specific surface area, molecular weight, mass density, and Faraday's constant, respectively. The variable J s (t) is the local volume current density used for the side reaction and described by Butler-Volmer kinetics. Assuming that the lithium reduction reaction where i 0,s is current density occupied by the lithium-loss side reaction; R is the universal gas constant; and T avg is the average temperature of the battery cell. The component η s (t) denotes the overpotential generated by the side reaction, which is expressed by the following equation: where ∆φ(t) denotes the difference between the solid phase potential and electrolyte potential. The variable U s,re f is the reference equilibrium potential, which is considered as a constant value. The total intercalation current J tot is used for the anode-side solution, which is presented by the sum of the current between the solid and solution J 1 , and the solution decreasing reaction and solution J s . Under three different charging currents (1 C, 2 C, and 3 C), the resistance film growth rates are computed according to Equations (23)-(27) and are presented in Figure 5. From this figure, we see that the film growth rate increases as the charging current or the time increases, owing to the augmentation of the reduction reaction rate and the rise in ∆φ(t) according to Equation (6). This means that a continuous high current and a high SoC is detrimental to reducing the lithium loss.

Temperature Control Based on the Backstepping Technique
The core temperature of a battery cell is higher than its surface temperature, hence, the control of the core temperature being more significant. Confining the maximum temperature of the battery to a certain value can guarantee both the safe operation and efficient charging. However, the core temperature of the battery is not accessible and cannot be measured by sensor directly. Therefore, we have to estimate the core temperature of a battery cell to make sure that the control goal is achieved. H. E. Perez et al. utilized the thermal model based on the equivalent circuit model to obtain the core temperature of a battery cell, which consists of the heat conduction resistance c R , convection resistance u R , core heat capacity c C , and surface heat capacity s C [25]. The thermal model used in this paper is described by the partial differential equation, which can accurately present the evolution of spatial temperature and is more accurate than the model employed by H. E. Perez et al.

Temperature Control Based on the Backstepping Technique
The core temperature of a battery cell is higher than its surface temperature, hence, the control of the core temperature being more significant. Confining the maximum temperature of the battery to a certain value can guarantee both the safe operation and efficient charging. However, the core temperature of the battery is not accessible and cannot be measured by sensor directly. Therefore, we have to estimate the core temperature of a battery cell to make sure that the control goal is achieved. H. E. Perez et al. utilized the thermal model based on the equivalent circuit model to obtain the core temperature of a battery cell, which consists of the heat conduction resistance R c , convection resistance R u , core heat capacity C c , and surface heat capacity C s [25]. The thermal model used in this paper is described by the partial differential equation, which can accurately present the evolution of spatial temperature and is more accurate than the model employed by H. E. Perez et al. To maintain the temperature of the battery cell where we want it to be, we can adjust the temperature or flow velocity of the cooling fluid. Herein, a backstepping method is utilized to design the control law for a prespecified time-invariant reference temperature T re f (r 1 ).
The following error system is thus introduced: where the expected temperature distribution T re f (r 1 ) is a quadratic parabola, which is in accordance with the true case; T sre f and T cre f are the reference surface temperature and core temperature, respectively. The time derivative of T e (r 1 , t) can be obtained according to Equations (12) and (28).
with the conditions ∂T e (0, t) where U is the boundary control law, which is applied to the boundary condition of Equation (12) for restricting T(r 1 , t) to around T re f (r 1 ).
Here, the following invertible backstepping transformation is adopted: which maps (29) into the stable system: where c > 0.5. K(r 1 , ρ) can be solved easily by the following kernel function: with the following conditions: The detailed solution process can be seen in [26]. According to the boundary condition w r 1 (R b , t) = 0 and Equation (29), we have Then, the control law is obtained: Energies 2020, 13, 2388 10 of 16 The control law requires access to the full states T(ρ, t). However, the core temperature cannot be measured by sensor directly. For that, the following boundary temperature observer is proposed: with boundary conditions The estimation error T(r 1 , t) = T(r 1 , t) −T(r 1 , t) is introduced, then substituting Equation (37) for Equation (12), one has with boundary conditions ∂ T(0, t) Using the following backstepping transform, which is similar with (30), we have where w(r 1 , t) satisfies the following stable system: where d > 0.5 as well. M(r 1 , ρ) satisfies the following kernel function: with the boundary conditions: The solution procedure of M(r 1 , ρ) is similar to K(r 1 , ρ).
To verify the performance of the backstepping-based observer and controller, a constant current of 2 C is applied to a battery cell; the curve of the temperature rising is presented in Figure 6. r = 0 and r = 1 represent the center and the surface of a battery cell, respectively. The two horizontal lines are the target temperature of the surface and the center of a battery cell. From Figure 6, we can see that the temperature estimation of both the surface and the center is capable of tracking the true temperature and approaching the target temperature gradually. That is to say that we can control the temperature to a safe range using the proposed temperature controller.
Energies 2020, 13, 2388 11 of 16 of 2 C is applied to a battery cell; the curve of the temperature rising is presented in Figure 6. 0 r = and 1 r = represent the center and the surface of a battery cell, respectively. The two horizontal lines are the target temperature of the surface and the center of a battery cell. From Figure 6, we can see that the temperature estimation of both the surface and the center is capable of tracking the true temperature and approaching the target temperature gradually. That is to say that we can control the temperature to a safe range using the proposed temperature controller.

Problem Formulation
This paper concentrates on the minimization of the charging time in a given target SoC with limits involving current, voltage, intercalation-induced stress, and film formation. The simultaneous nonlinear programming approach is used in this paper. Considering the optimal charging profile with predetermined final SoC under the objective of minimization of charging time, the objective function J is given by

Problem Formulation
This paper concentrates on the minimization of the charging time in a given target SoC with limits involving current, voltage, intercalation-induced stress, and film formation. The simultaneous nonlinear programming approach is used in this paper. Considering the optimal charging profile with predetermined final SoC under the objective of minimization of charging time, the objective function J is given by where I(t) is the input current (A); t f is the final time (s); V(t) is the terminal voltage (V); σ − r,upper (t) and σ − t,upper are the radial and tangential stress upper bounds, respectively; and ∂ . δ f ilm_max (t) is the upper bound of the film growth rate.

Simulation and Results
In order to verify the performance of the proposed optimal fast-charging method, CC-CV, the most common pattern of fast charging, was used for the comparative study, including four aspects: charging time, maximum tangential stress, maximum radial stress, and film growth rate.
In the process of fast charging a lithium battery, the temperature of the battery cell grows quickly and the aging rate increases. To avoid overheating, the backstepping method is used to keep the battery cell at a constant temperature, as shown in Figure 7. A constant temperature in the battery cell not only helps avoid overheating, but also helps compare the performance of the proposed fast-charging method with the CC-CV method. In order to find a compromise between the charging time and battery aging, the proposed optimal fast-charging method in this paper considers two aging factors: intercalation-induced stress and film growth. In the comparative study, the CC-CV pattern adopts three kinds of maximum charging current: 1 C, 2 C, and 3 C. This means that the effect of the choice of the maximum charging current on the comparative results is avoided.
In this paper, the upper limits of the radial and tangential stress are predetermined and optimal charging curves were obtained using the dynamic optimization approach. As we know, when charging occurs, the maximum tangential stress is located at the surface of the particle, and the maximum radial stress is located at the center of the particle. For an anode made of graphite, neither the maximum tangential stress nor the maximum radial stress is not allowed to exceed the yielding stress for any length of time to reduce the risk of anode fracture. The following case is discussed: charging a battery from 0% SoC to 100% SoC. battery cell at a constant temperature, as shown in Figure 7. A constant temperature in the battery cell not only helps avoid overheating, but also helps compare the performance of the proposed fastcharging method with the CC-CV method. In order to find a compromise between the charging time and battery aging, the proposed optimal fast-charging method in this paper considers two aging factors: intercalation-induced stress and film growth. In the comparative study, the CC-CV pattern adopts three kinds of maximum charging current: 1 C, 2 C, and 3 C. This means that the effect of the choice of the maximum charging current on the comparative results is avoided.
In this paper, the upper limits of the radial and tangential stress are predetermined and optimal charging curves were obtained using the dynamic optimization approach. As we know, when charging occurs, the maximum tangential stress is located at the surface of the particle, and the maximum radial stress is located at the center of the particle. For an anode made of graphite, neither the maximum tangential stress nor the maximum radial stress is not allowed to exceed the yielding stress for any length of time to reduce the risk of anode fracture. The following case is discussed: charging a battery from 0% SoC to 100% SoC.  Figure 7 presents the balance of temperature distribution. If this balance is broken, the temperature controller would work to drag the temperature to the target position using various methods, such as changing the fluid temperature and adjusting the fluid flow velocity. Thus, a stable cell temperature can be achieved for safe operation. Figures 8 and 9 show the current and voltage curves of charging a battery from 0% SoC to 100% SoC by utilizing the CC-CV pattern under three maximum charging currents. At the initial charging stage, the input current maintains at the maximum and the terminal voltage climbs fast. When the output voltage reaches the upper value, this voltage is held until the battery cell is charged fully. With the increase in the maximum charging current, the time cost for full charging becomes less and less. Figure 10 indicates the optimal charging current obtained by the dynamic optimization method. The upper value and lower values of the constraints are listed in Table 1.   Figure 7 presents the balance of temperature distribution. If this balance is broken, the temperature controller would work to drag the temperature to the target position using various methods, such as changing the fluid temperature and adjusting the fluid flow velocity. Thus, a stable cell temperature can be achieved for safe operation. Figures 8 and 9 show the current and voltage curves of charging a battery from 0% SoC to 100% SoC by utilizing the CC-CV pattern under three maximum charging currents. At the initial charging stage, the input current maintains at the maximum and the terminal voltage climbs fast. When the output voltage reaches the upper value, this voltage is held until the battery cell is charged fully. With the increase in the maximum charging current, the time cost for full charging becomes less and less. Figure 10 indicates the optimal charging current obtained by the dynamic optimization method. The upper value and lower values of the constraints are listed in Table 1.    From Figure 11, we see that the area under the film growth curve is smaller than those of the other three curves, which means the least lithium loss during full charging. At the same time, the charging time of the optimal charging profile just takes 200 s more than the time cost under CC-CV with a 3 C current, which is fast enough and acceptable.

( )
The yielding stress of the radial stress and tangential stress in a graphite-based anode are both around 30 MPa. If the radial stress or tangential stress is more than 30 Mpa for a long time, it is possible to cause fatigue failure and accelerate the aging process. From Figure 12, the tangential stress generated by optimal charging profile becomes greater than its yielding stress within 1400 s, while the tangential stress under CC-CV with 1 C, 2 C, and 3 C is larger than 30 MPa for more than 1700 s,     From Figure 11, we see that the area under the film growth curve is smaller than those of the other three curves, which means the least lithium loss during full charging. At the same time, the charging time of the optimal charging profile just takes 200 s more than the time cost under CC-CV with a 3 C current, which is fast enough and acceptable.
The yielding stress of the radial stress and tangential stress in a graphite-based anode are both around 30 MPa. If the radial stress or tangential stress is more than 30 Mpa for a long time, it is possible to cause fatigue failure and accelerate the aging process. From Figure 12, the tangential stress generated by optimal charging profile becomes greater than its yielding stress within 1400 s, while the tangential stress under CC-CV with 1 C, 2 C, and 3 C is larger than 30 MPa for more than 1700 s,   Range From Figure 11, we see that the area under the film growth curve is smaller than those of the other three curves, which means the least lithium loss during full charging. At the same time, the charging time of the optimal charging profile just takes 200 s more than the time cost under CC-CV with a 3 C current, which is fast enough and acceptable.
Energies 2020, 13, 03 FOR PEER REVIEW 15 of 18 2200 s, and 3700 s, respectively. From the maximum values of radial stress as shown in Figure 13, we can see that the optimal charging profile causes a shorter period of overstress, which represents less electrode fatigue damage. Hence, the result shows that the optimal charging profile obtained by the dynamic optimization method has less probability of causing stress than the CC-CV profiles and so can slow the aging process and extend the battery life. Therefore, while optimizing the charging speed, the proposed optimal charging method offers a good compromise between charging time and battery aging.   The yielding stress of the radial stress and tangential stress in a graphite-based anode are both around 30 MPa. If the radial stress or tangential stress is more than 30 Mpa for a long time, it is possible to cause fatigue failure and accelerate the aging process. From Figure 12, the tangential stress generated by optimal charging profile becomes greater than its yielding stress within 1400 s, while the tangential stress under CC-CV with 1 C, 2 C, and 3 C is larger than 30 MPa for more than 1700 s, 2200 s, and 3700 s, respectively. From the maximum values of radial stress as shown in Figure 13, we can see that the optimal charging profile causes a shorter period of overstress, which represents less electrode fatigue damage. Hence, the result shows that the optimal charging profile obtained by the dynamic optimization method has less probability of causing stress than the CC-CV profiles and so can slow the aging process and extend the battery life.
dynamic optimization method has less probability of causing stress than the CC-CV profiles and so can slow the aging process and extend the battery life.
Therefore, while optimizing the charging speed, the proposed optimal charging method offers a good compromise between charging time and battery aging.

Conclusions
In this paper, the electrochemical-thermal model was employed to obtain the optimal charging profile and control the temperature of the cell. A temperature controller-based backstepping method can slow the aging process and extend the battery life.
Therefore, while optimizing the charging speed, the proposed optimal charging method offers a good compromise between charging time and battery aging.

Conclusions
In this paper, the electrochemical-thermal model was employed to obtain the optimal charging profile and control the temperature of the cell. A temperature controller-based backstepping method Therefore, while optimizing the charging speed, the proposed optimal charging method offers a good compromise between charging time and battery aging.

Conclusions
In this paper, the electrochemical-thermal model was employed to obtain the optimal charging profile and control the temperature of the cell. A temperature controller-based backstepping method was proposed to keep a relatively constant cell temperature to avoid overheating. Then, the effects of intercalation-induced stress were considered, because these cause particle fracture, which can accelerate the aging process. Furthermore, the growth of film on the surface of particles was also taken into account due to its ability to reduce the reused lithium ions. After that, the optimization objective was established, which minimizes the time cost during charging, while confining intercalation-induced stress and the growth rate of surficial film to a given range. Finally, the simulation was implemented and the results show that the film growth of the optimal charging curve is smaller than that of the other three charging curves, thus demonstrating the least lithium loss during full charging. Besides that, the tangential stress and the radial stress generated by the optimal charging profile were both greater than their yielding stresses for less time during charging as compared with the CC-CV pattern, which means there is less risk of electrode fatigue fracture. Furthermore, the charging time of the optimal charging profile just takes 200 s more than the time cost under CC-CV with a 3 C current, which is fast enough and acceptable. Hence, the proposed optimal charging method offers a good compromise between charging time and battery aging.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: BMS