High-Fidelity Battery Model for Model Predictive Control Implemented into a Plug-In Hybrid Electric Vehicle

: Power management strategies have impacts on fuel economy, greenhouse gasses (GHG) emission, as well as effects on the durability of power-train components. This is why different off-line and real-time optimal control approaches are being developed. However, real-time control seems to be more attractive than off-line control because it can be directly implemented for managing power and energy ﬂows inside an actual vehicle. One interesting illustration of these power management strategies is the model predictive control (MPC) based algorithm. Inside a MPC, a cost function is optimized while system constraints are validated in real time. The MPC algorithm relies on dynamic models of the vehicle and the battery. The complexity and accuracy of the battery model are usually neglected to beneﬁt the development of new cost functions or better MPC algorithms. The contribution of this manuscript consists of developing and evaluating a high-ﬁdelity battery model of a plug-in hybrid electric vehicle (PHEV) that has been used for MPC. Via empirical work and simulation, the impact of a high-ﬁdelity battery model has been evaluated and compared to a simpler model in the context of MPC. It is proven that the new battery model reduces the absolute voltage, state of charge ( SoC ), and battery power loss error by a factor of 3.2, 1.9 and 2.1 on average respectively, compared to the simpler battery model.


Introduction
To reduce fuel consumption, greenhouse gasses (GHG) emission or to improve battery lifetime, advanced power management strategies have been studied [1,2]. To achieve these goals, different off-line and real-time optimal control approaches are being developed. Model predictive control (MPC, see Figure 1 [3]) is a promising example. MPC is an advanced control methodology that was originally used in the process industries for power plants and oil refineries since the 1980s [4,5]. MPC is capable of predicting future events, in a finite-horizon, and is able to take controlled actions to optimize a cost function. Proportional, integrator, and derivative (PID) and linear quadratic regulator (LQR) controllers do not have this predictive ability. However, the prediction of MPC is usually computed for a relatively short time horizon in the future due to its computational cost. In recent years, MPC has been used for hybrid electric vehicle (HEV) and plug-in HEV (PHEV) applications, as shown in the comprehensive survey of power management topics [6]. In [6], dynamic programming-based MPC relies on dynamic models of the vehicle, most often linear empirical models obtained by system identification [25]. In the references cited above, the significance of a MPC controller for HEV applications has been highlighted. In this literature review, the authors mainly focus on the influence of a new cost function, but they do not deal with the benefit of a high accuracy battery model. Their models are far simpler than the one introduced in papers [ [26][27][28] for automotive applications. For instance, their impedance battery model is limited to an internal resistance impedance model, as shown in Figure 2. Developing a high-fidelity battery model for a HEV is a key opportunity to improve the global performance of a MPC. This is why the main contribution of this manuscript is developing and evaluating a high-fidelity battery model of a PHEV that has been used for a MPC. In Section 2 of this paper, a Matlab model based on a real PHEV is presented. It provides an overview of the PHEV model and helps to understand the role of the battery model. Then, Section 3 is entirely devoted to the description of the high accuracy battery model for a MPC. It precisely simulates the battery voltage response through one resistance in series with a two resistances/capacitors (RC) network, and replicates the battery aging phenomenon. Finally, in Section 4, the impact of the high-fidelity battery model has been evaluated through empirical work and simulation and compared to the simple internal resistance battery model in the context of MPC.
Batteries 2017, 3, 13 2 of 16 are defined as the primary on-line methods. The former has been of interest in many works. For example, using a MPC, an overall power efficiency is maximized, instead of fuel consumption minimization, for a series HEV [10]. In [11], based on the price of gas and electricity in the United States, a MPC minimizes the cost of the vehicle's energy use for a series PHEV. In another research report, engine transient characteristic is incorporated in a MPC for parallel HEV [12]. In other publications [13,14], the goal of MPC is to reduce the CO2 emissions. To implement a fast MPC, the authors in [14] proposed to compute the entire control law off-line. A global optimization-based MPC is developed, and experimental validation is provided on the test bench in [15]. Stochastic MPC with learning is proposed and validated by simulation in [16]. The MPC-based control strategy is also applied to solve the energy management problem of a series [9,17] and power-split HEVs in [18] and to analyze the potential benefits of integrating ultracapacitors (UC) in the energy storage system (ESS) unit of a power-split HEV in [19][20][21][22]. Using the Pontryagin's minimum principle (PMP), the MPC of fuel cell hybrid vehicles (FCHVs) optimizes battery lifetime while reducing battery energy loss, fuel consumption, and powertrain cost in [23,24]. MPC relies on dynamic models of the vehicle, most often linear empirical models obtained by system identification [25]. In the references cited above, the significance of a MPC controller for HEV applications has been highlighted. In this literature review, the authors mainly focus on the influence of a new cost function, but they do not deal with the benefit of a high accuracy battery model. Their models are far simpler than the one introduced in papers [ [26][27][28] for automotive applications. For instance, their impedance battery model is limited to an internal resistance impedance model, as shown in Figure 2. Developing a high-fidelity battery model for a HEV is a key opportunity to improve the global performance of a MPC. This is why the main contribution of this manuscript is developing and evaluating a high-fidelity battery model of a PHEV that has been used for a MPC. In Section 2 of this paper, a Matlab model based on a real PHEV is presented. It provides an overview of the PHEV model and helps to understand the role of the battery model. Then, Section 3 is entirely devoted to the description of the high accuracy battery model for a MPC. It precisely simulates the battery voltage response through one resistance in series with a two resistances/capacitors (RC) network, and replicates the battery aging phenomenon. Finally, in Section 4, the impact of the highfidelity battery model has been evaluated through empirical work and simulation and compared to the simple internal resistance battery model in the context of MPC.

Vehicle Model
In this section, the architecture of the reference car is shown in Figure 3. It provides an overview of the vehicle model and helps to understand the role of the battery model. The vehicle is a series PHEV. The arrows in Figure 3 indicate the possible direction of the energy flow inside the power-train. An electric generator, an ESS composed of a lithium-ion battery, and an electric motor are connected to a DC bus. During acceleration, the generator or the ESS provide power to the DC bus for the motor while, during regenerative braking, the motor works as a generator and provides power to the ESS through the DC bus. Except for cranking the engine at the very start of a drive, there is no energy transfer from the battery to the engine through the generator while the car is running. Consequently, the energy flow in this situation is minimal and considered as negligible compared to the total energy consumption. This is why the arrows between engine and generator, and between generator and battery are uni-directional. Table 1 summarizes the specifications of the power-train components. The size of the ESS has been selected based on a parametric study on the available size of high energy density modules from A123 (Livona, MI, USA) to minimize fuel consumption [26]. High energy density modules are preferred to high power density modules because the ESS of a PHEV is usually sufficient to meet the power requirements of the vehicle power-train. Therefore, the ESS should be able to store more energy to increase the all-electric range and, as a result, to reduce the fuel consumption of the vehicle. However, the available space inside the concept vehicle test platform (Subaru BRZ 2015) constrains the maximum ESS size. The size of the engine has been selected based on the reduced size engine concept to achieve maximal fuel economy as explained in [29]. A motor unit power has been chosen to keep the acceleration performance of this new PHEV close to its original performance. The unit includes two independent motors for both rear wheels.

Vehicle Model
In this section, the architecture of the reference car is shown in Figure 3. It provides an overview of the vehicle model and helps to understand the role of the battery model. The vehicle is a series PHEV. The arrows in Figure 3 indicate the possible direction of the energy flow inside the powertrain. An electric generator, an ESS composed of a lithium-ion battery, and an electric motor are connected to a DC bus. During acceleration, the generator or the ESS provide power to the DC bus for the motor while, during regenerative braking, the motor works as a generator and provides power to the ESS through the DC bus. Except for cranking the engine at the very start of a drive, there is no energy transfer from the battery to the engine through the generator while the car is running. Consequently, the energy flow in this situation is minimal and considered as negligible compared to the total energy consumption. This is why the arrows between engine and generator, and between generator and battery are uni-directional. Table 1 summarizes the specifications of the power-train components. The size of the ESS has been selected based on a parametric study on the available size of high energy density modules from A123 (Livona, MI, USA) to minimize fuel consumption [26]. High energy density modules are preferred to high power density modules because the ESS of a PHEV is usually sufficient to meet the power requirements of the vehicle power-train. Therefore, the ESS should be able to store more energy to increase the all-electric range and, as a result, to reduce the fuel consumption of the vehicle. However, the available space inside the concept vehicle test platform (Subaru BRZ 2015) constrains the maximum ESS size. The size of the engine has been selected based on the reduced size engine concept to achieve maximal fuel economy as explained in [29]. A motor unit power has been chosen to keep the acceleration performance of this new PHEV close to its original performance. The unit includes two independent motors for both rear wheels. The characteristics of the vehicle model are given in Table 1 and Figure 3. The equations describing the vehicle dynamics, the electric motor and generator, and the engine are included in the Appendix.  The characteristics of the vehicle model are given in Table 1 and Figure 3. The equations describing the vehicle dynamics, the electric motor and generator, and the engine are included in the Appendix A.

Impedance
To capture the battery dynamic accurately, an internal series resistance and a two RC network, as shown in Figure 4, are commonly used to design battery impedance for EV and HEV applications [30][31][32]. The value of each resistance and capacity are based on [28,33]. Indeed, the battery used in the vehicle is the same as the one used in [28,33]: A123 LiFePO 4 prismatic module with its own battery management system (BMS). Only the configuration changes: from 5×22s3p in [28,33] to 7×15s2p in the Subaru BRZ 2015. Table 2 summarizes the characteristics of the battery impedance obtained at laboratory ambient temperature.

Power, Voltage and Current
The power provided or received by the battery is computed as follows: Batteries 2017, 3, 13 4 of 16

Impedance
To capture the battery dynamic accurately, an internal series resistance and a two RC network, as shown in Figure 4, are commonly used to design battery impedance for EV and HEV applications [30][31][32]. The value of each resistance and capacity are based on [28,33]. Indeed, the battery used in the vehicle is the same as the one used in [28,33]: A123 LiFePO4 prismatic module with its own battery management system (BMS). Only the configuration changes: from 5×22s3p in [28,33] to 7×15s2p in the Subaru BRZ 2015. Table 2 summarizes the characteristics of the battery impedance obtained at laboratory ambient temperature.

Power, Voltage and Current
The power provided or received by the battery is computed as follows: In Equations (1) and (2), the battery is charging when Pbatt is negative and is discharging when Pbatt is positive. The battery current is then computed as follows: Then, new voltage is computed using the impedance model of the battery in a discrete time domain as follows: In Equations (1) and (2), the battery is charging when P batt is negative and is discharging when P batt is positive. The battery current is then computed as follows: Then, new voltage is computed using the impedance model of the battery in a discrete time domain as follows: The open circuit voltage of the battery, V OC , is an eight-order polynomial equation varying over state of charge (SoC) representing the SoC-V OC curve shown in Figure 5. The data and polynomial order are extracted from [34]. The open circuit voltage of the battery, VOC, is an eight-order polynomial equation varying over state of charge (SoC) representing the SoC-VOC curve shown in Figure 5. The data and polynomial order are extracted from [34]. Here, it is noted that Ibatt(t) is computed by using Vbatt(t − 1) whereas, in theory, Vbatt(t) is necessary. It is carried out this way because, at this stage, Vbatt(t) is unknown. The error induced by this approximation should be small considering that the voltage fluctuates slowly. However, this error exists and a simple recursive algorithmic method has been developed to minimize it. Indeed, once Vbatt(t) is computed by Equations (4)-(6), it is possible to recompute the value of Ibatt(t) using Equation (3) by replacing Vbatt(t − 1) with Vbatt(t). Then, a new value of Vbatt(t) is also recomputed using Equation (4)-(6) again. Then, the same set of instructions can be repeated as many times as desired. The purpose of this algorithm is to reduce the error due to the approximation implied by Equation ε with  a value as small as wanted. Figure 6 summarizes the current and voltage computation recursive algorithm.

State of Charge
The SoC of the battery is computed by the coulomb counting Equation (7) in the discrete time domain.
1 ∆ (7) Here, it is noted that I batt (t) is computed by using V batt (t − 1) whereas, in theory, V batt (t) is necessary. It is carried out this way because, at this stage, V batt (t) is unknown. The error induced by this approximation should be small considering that the voltage fluctuates slowly. However, this error exists and a simple recursive algorithmic method has been developed to minimize it. Indeed, once V batt (t) is computed by Equations (4)-(6), it is possible to recompute the value of I batt (t) using Equation (3) by replacing V batt (t − 1) with V batt (t). Then, a new value of V batt (t) is also recomputed using Equation (4)-(6) again. Then, the same set of instructions can be repeated as many times as desired. The purpose of this algorithm is to reduce the error due to the approximation implied by Equation (3). This recursive algorithm only stops when ≤ ε with ε a value as small as wanted. Figure 6 summarizes the current and voltage computation recursive algorithm. The open circuit voltage of the battery, VOC, is an eight-order polynomial equation varying over state of charge (SoC) representing the SoC-VOC curve shown in Figure 5. The data and polynomial order are extracted from [34]. Here, it is noted that Ibatt(t) is computed by using Vbatt(t − 1) whereas, in theory, Vbatt(t) is necessary. It is carried out this way because, at this stage, Vbatt(t) is unknown. The error induced by this approximation should be small considering that the voltage fluctuates slowly. However, this error exists and a simple recursive algorithmic method has been developed to minimize it. Indeed, once Vbatt(t) is computed by Equations (4)-(6), it is possible to recompute the value of Ibatt(t) using Equation (3) by replacing Vbatt(t − 1) with Vbatt(t). Then, a new value of Vbatt(t) is also recomputed using Equation (4)-(6) again. Then, the same set of instructions can be repeated as many times as desired. The purpose of this algorithm is to reduce the error due to the approximation implied by Equation ε with  a value as small as wanted. Figure 6 summarizes the current and voltage computation recursive algorithm.

State of Charge
The SoC of the battery is computed by the coulomb counting Equation (7) in the discrete time domain.

State of Charge
The SoC of the battery is computed by the coulomb counting Equation (7) in the discrete time domain.

State of Health
The battery state of health (SoH) is represented by the capacity fade of the battery over cycling. Equations (8)- (13) simulate this phenomenon as a function of the cumulative SoC variation, C-rate, and temperature. Those formulas are extracted from [34] because their tested batteries have the same chemistry and come from the same manufacturer as the battery used in this work, which is the A123 LiFePO 4 battery. This is why it is assumed that the aging characteristic of the large-format battery pack is similar to the cells tested in [35] despite the fact that the battery pack is composed of seven modules of 15 prismatic cells in series and two in parallel and not one cylindrical cell. For example, this assumption was accepted as a correct battery aging model for a causal optimal control-based energy management strategy for a parallel HEV using a lithium ion battery pack in [35]. The following equations explain how the SoH is computed in a discrete domain time.
According to Table 3 of Wang's paper [34], Equation (8) should normally be written this way: Also, according to Equations (3) and (7) of the same paper [34], Equation (11) can be written and z is equal to 0.55. Moreover, B values are a function of C rate , defined as the ratio of the absolute value of the current to the cell capacity (Equation (12)). Wang's paper [34] only provides B values for four different C rate : C/2, 2C, 6C, and 10C. In this paper, by default, the B value is linearly interpolated in function of C rate using data provided in Table 3 in [34], which leads to Equation (10).
By Equation (14) provided in Table 3 in [34], A h is defined as a SoC variation. Indeed, the SoC variation of a battery can be concluded with the following Equation (15): Since this battery is managed by its own BMS, which includes cells balancing, it allows us to hypothesize the following: the aging of battery cells should be roughly homogenous throughout the pack. From this hypothesis, it can be concluded that two of the same batteries, with different capacity, subjected to the same charge and discharge conditions (C-rate, temperature, number of cycles, depth of discharge (DOD), etc.) should be impacted the same way by the capacity fade. In other words, the results published in paper [34] would have been the same if they tested cells with different capacity but with the same chemistry (LiFePO 4 ) and manufacturer (A123).

Model Validation
On the dynamometer, the vehicle has been driven over an HWFET drive cycle from full battery charge (96% of SoC) to its complete depletion (5% of SoC). The goal of this test was to validate the accuracy of the battery model. This is why, during this experiment, the engine and electric generator was not used. The effects of the road slope and the wind speed have not been considered, and so those parameters have been set to zero. The air density has been set to a constant value. At the time of the test, the coolant system, supposed to keep the battery temperature optimal, was not present. However, the internal BMS of the A123 LiFePO 4 prismatic module records the temperature of the battery and provides an alarm response when it goes over the limits. During the experiment, it goes from 25 • C to 38.5 • C. Figure 7 summarizes the experiment.

Electric Motor
During the test, the speed of the car was controlled by a human driver operating an accelerator and brake pedal. The driver tried to follow the HWFET drive cycle as closely as possible. However, pedal sensitivity limitation and driver reactivity cause some inaccuracy that might explain the peaks in error between experimental and simulated motor speed and torque. This is why the vehicle speed profile during this experiment was recorded and used as an input for the simulation.
The goal of those preliminary results is to check the accuracy of the vehicle dynamics and electric motor models, described in Appendix. In fact, a poor design would prevent the evaluation of the performance of the high-fidelity battery model. However, results in Figures 8 and 9 show high accuracy of the vehicle dynamics and motor models with a low average absolute error for the motor speed and good precision for the motor torque (see Table 3).

Electric Motor
During the test, the speed of the car was controlled by a human driver operating an accelerator and brake pedal. The driver tried to follow the HWFET drive cycle as closely as possible. However, pedal sensitivity limitation and driver reactivity cause some inaccuracy that might explain the peaks in error between experimental and simulated motor speed and torque. This is why the vehicle speed profile during this experiment was recorded and used as an input for the simulation.
The goal of those preliminary results is to check the accuracy of the vehicle dynamics and electric motor models, described in Appendix A. In fact, a poor design would prevent the evaluation of the performance of the high-fidelity battery model. However, results in Figures 8 and 9 show high accuracy of the vehicle dynamics and motor models with a low average absolute error for the motor speed and good precision for the motor torque (see Table 3). profile during this experiment was recorded and used as an input for the simulation. The goal of those preliminary results is to check the accuracy of the vehicle dynamics and electric motor models, described in Appendix. In fact, a poor design would prevent the evaluation of the performance of the high-fidelity battery model. However, results in Figures 8 and 9 show high accuracy of the vehicle dynamics and motor models with a low average absolute error for the motor speed and good precision for the motor torque (see Table 3).

Battery
The results of those experiments show that the recursive algorithm does not bring any additional precision to the battery model compared to the two RC network battery model by itself. Nevertheless, they do prove that the two RC network battery model improves the accuracy of the estimation of the voltage and SoC of the battery, as shown in Figures 10 and 11 and Tables 4 and 5, compared to the simple internal resistance battery model. However, the impact of the battery model on the current is minimal, as shown in Figure 12 and Table 6. Those observations have a theoretical explanation. We know that Ibatt = Pbatt/Vbatt. As Vbatt fluctuations are small compared to Pbatt, Ibatt mainly depends on Pbatt. In a high-level abstraction, we could say that Ibatt = f(Pbatt) as Vbatt is quite constant. This is why precision of the battery model will not have much effect on Ibatt. We can also deduce that error on Ibatt

Battery
The results of those experiments show that the recursive algorithm does not bring any additional precision to the battery model compared to the two RC network battery model by itself. Nevertheless, they do prove that the two RC network battery model improves the accuracy of the estimation of the voltage and SoC of the battery, as shown in Figures 10 and 11 and Tables 4 and 5, compared to the simple internal resistance battery model. However, the impact of the battery model on the current is minimal, as shown in Figure 12 and Table 6. Those observations have a theoretical explanation. We know that I batt = P batt /V batt . As V batt fluctuations are small compared to P batt , I batt mainly depends on P batt . In a high-level abstraction, we could say that I batt = f(P batt ) as V batt is quite constant. This is why precision of the battery model will not have much effect on I batt . We can also deduce that error on I batt mainly depends on the vehicle dynamics and electric motor models because P batt clearly depends on those equations (see Equations (1), (2) and (A1)-(A16)). However, the precision of the battery model will have an impact on V batt . Indeed, this battery model is made to simulate the battery response voltage due to current changes. This is why V batt is influenced by the model. Concerning the SoC, it is computed by integrating I batt . Even if the battery model has a locally small impact on I batt , there is a difference between both calculated currents (see Table 5). By integrating those differences, the error, regarding the SoC estimation, accumulates. This explains why the SoC of the two RC network and the internal resistance battery model are slowly diverging from each other while discharging, as shown in Figure 12.          Figure 12. Battery current comparison between the experimental and simulated battery model. For experimental data, V OC has not been measured during the test. However, as the V OC used with the two RC Network battery model is originally derived from experimental data (see Figure 5), the same V OC has been used for computing the battery loss in experimental data. The results, shown in Figures 13 and 14 and Table 7, prove that the simulated losses from the two RC Network battery model are closer to reality than the one simulated by the internal resistance battery model. This observation suggests that the improvement of the battery accuracy may impact the MPC controller. Indeed, as the simulated battery power loss is closer to reality, the power requested to the engine will change and be closer to the real need of the vehicle. Besides, it is noted that the two RC network battery model estimates the bus voltage better. Therefore, when this model is used for the battery in the power-train model of a MPC algorithm, the variables such as voltage, current, torque, and speed of other components, such as engine, generator, and traction motor, will change. Thus, it can be predicted that, in a power-train model of a MPC where a two RC network battery model is deployed, a better estimation of actual variables of the system can be accomplished. change and be closer to the real need of the vehicle. Besides, it is noted that the two RC network battery model estimates the bus voltage better. Therefore, when this model is used for the battery in the power-train model of a MPC algorithm, the variables such as voltage, current, torque, and speed of other components, such as engine, generator, and traction motor, will change. Thus, it can be predicted that, in a power-train model of a MPC where a two RC network battery model is deployed, a better estimation of actual variables of the system can be accomplished.  The employed method to compute both battery power and cumulative losses is described by Equations (16) and (17): Figure 13. Power loss comparison between the experimental data and simulated battery model.
The employed method to compute both battery power and cumulative losses is described by Equations (16) and (17): Eb loss =

Conclusions
In this manuscript, a Matlab-based model has been developed with a high accuracy battery model for a hybridized Subaru BRZ 2015. This manuscript has two contributions. The first contribution is the development, for the first time in the field of a MPC controller, of a high accuracy

Conclusions
In this manuscript, a Matlab-based model has been developed with a high accuracy battery model for a hybridized Subaru BRZ 2015. This manuscript has two contributions. The first contribution is the development, for the first time in the field of a MPC controller, of a high accuracy battery for a PHEV Matlab-based model. Indeed, the new battery model reduces the absolute voltage, SoC and battery power loss error respectively by a factor of 3.2, 1.9 and 2.1 on average compared to the simpler battery model. Those observations highly infer that the performance of MPC will be impacted by the improvement of the high-fidelity battery model. It is likely that the high-fidelity battery model helps to reach more accurate estimation of the behaviors of other power-train components, and thus helps MPC to manage more precisely fuel consumption, battery lifetime, vehicle efficiency, or any other implemented cost function. However, there is no concrete proof that the improvement of the battery model fidelity is necessary to obtain an optimal control solution from the MPC. This is why the second contribution is the question raised in this manuscript: what is the optimal battery model fidelity required to optimize the performance of a MPC controller in the case of automotive application? To answer this question, the first sensitivity analysis of the battery model for MPC is planned to be conducted in a separate work. Following the control loop of Figure 15, this analysis should provide a procedure and results to determine the optimal tradeoff between the complexity and accuracy of the battery model of a PHEV for a MPC controller.

Conclusions
In this manuscript, a Matlab-based model has been developed with a high accuracy battery model for a hybridized Subaru BRZ 2015. This manuscript has two contributions. The first contribution is the development, for the first time in the field of a MPC controller, of a high accuracy battery for a PHEV Matlab-based model. Indeed, the new battery model reduces the absolute voltage, SoC and battery power loss error respectively by a factor of 3.2, 1.9 and 2.1 on average compared to the simpler battery model. Those observations highly infer that the performance of MPC will be impacted by the improvement of the high-fidelity battery model. It is likely that the high-fidelity battery model helps to reach more accurate estimation of the behaviors of other power-train components, and thus helps MPC to manage more precisely fuel consumption, battery lifetime, vehicle efficiency, or any other implemented cost function. However, there is no concrete proof that the improvement of the battery model fidelity is necessary to obtain an optimal control solution from the MPC. This is why the second contribution is the question raised in this manuscript: what is the optimal battery model fidelity required to optimize the performance of a MPC controller in the case of automotive application? To answer this question, the first sensitivity analysis of the battery model for MPC is planned to be conducted in a separate work. Following the control loop of Figure 15, this analysis should provide a procedure and results to determine the optimal tradeoff between the complexity and accuracy of the battery model of a PHEV for a MPC controller.

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

Symbol
Name Units t The discrete time s ∆t The step time s F r The rolling resistance N F g The grading resistance N F w The aerodynamic drag resistance N F a The acceleration resistance N F c The sum of every resistant force applied to the car N β The lateral angle between the wind and car direction rads ϑ The vertical angle between the wind and car direction rads α The road elevation rads The power provided/received by the battery W P m_elec The electric power requested/supply by the electric motor W P g_elec The power provide by the generator W P a The constant power consumed by the auxiliary electric system W P batt The power provided/received by the battery W P c The power requested/provided by the car W P The current provided/received by the battery A SoC The battery state of charge -C batt The battery nominal capacity Ah Q loss The battery capacity fades -Pb loss The battery power loss W Eb loss The battery energy loss Wh ∆SoC The cumulative SoC variation for a given C-rate -DOD The depth of discharge n cycle The number of battery cycle -C r0 The car static rolling coefficient - The car dynamic rolling coefficient - The car drag coefficient ε 0 The gear reduction rd The wheels radius m A f The car front size surface m 2 ρ a The air density kg/m 3 T ext The ambient absolute temperature K P r The ambient pressure Pa R air The specific gas constant for dry air J/(kg·K) R The perfect gas constant J/mol·K E a The activation energy J/mol z The power law factor -B The pre-exponent factor -Appendix A.

Appendix A.1. Vehicle Dynamics
This section describes the vehicle dynamic equations [36][37][38][39]. The power that the car's power-train has to provide can be computed as follows: The rolling resistance of the vehicle is defined by Equation (A2): Mg cos(α(t)) (A2) The static and dynamic rolling coefficients have been experimentally determined. The grading resistance of the vehicle is defined by Equation (A3): F g (t) = Mg sin(α(t)) (A3) The aerodynamics drag resistance is defined as follows: The acceleration resistance of the vehicle is defined by Equation (A7):