Li-Ion Battery Lifetime Model ’ s Influence on the Economic Assessment of a Hybrid Electric Bus ’ s Operation

The present paper is focused on the evaluation of the economic influence of a battery lifetime model upon the optimal sizing and energy management strategy of a dual energy storage system (ESS) composed of Lithium-ion batteries and supercapacitors. The operation of a Hybrid Bus is taken as a case study in order to evaluate the effects of battery lifetime models’ accuracy on ESS sizing and operation in a heavy-duty application. For this purpose, two different lifetime models (a Wöhler-curve-based model and a semi-empirical model) were applied in the multi-objective optimisation of a hybrid electric urban bus. Differences up to ca. 8% on the daily operation costs and ca. 25% on the dual ESS costs were estimated depending on the lifetime model considered for the optimisation.


Introduction
The development of new electrochemical energy storage technologies and the growing maturity of lithium-ion (Li-ion) batteries are promoting the penetration of advanced electro-mobility solutions.In the urban context, one of the most promising alternatives to reduce emissions is the use of hybrid electric buses (HEB), which are considered a more economically viable alternative than full-electric buses in the medium term [1].
The implementation of battery-based energy storage systems (ESSs) entails inherent challenges, since the reduction of fuel consumption is related to a higher use of the batteries.This leads to a great number of battery charge/discharge cycles and thus a potential early degradation [2].To cover the lifetime of a heavy-duty application, e.g., in the case of urban buses, multiple battery replacements might be necessary.Therefore, extending battery lifetime, both in terms of calendar and cycle life, would be crucial to make electric or hybrid buses a cost-competitive alternative to conventional buses [3].To solve this issue, the use of a dual ESS, composed of batteries (BTs) and supercapacitors (SCs), is investigated in this paper.A dual ESS may reduce the power demands on the batteries by using a complementary high power energy storage device (in this case SC) which shaves sudden power peaks.
The sizing of the ESS [4][5][6] is commonly addressed during the design stage of a vehicle.Despite the fact that, in a dual ESS based on BTs and SCs, the function of each energy storage unit during operation is quite defined (energy provided by the BT pack and power peak regulation by the SC pack), a proper methodology is still required to define their suitable sizing and to obtain efficient and economic performance.Thus, it can be noted that the system's sizing and operation are key factors to ensure the fulfillment of the application requirements, the proper harnessing of the onboard energy sources, and the defined economic objectives [7,8].In close-to-market applications, the economic issue becomes in an important decision factor as the solution has to be competitive and profitable both for the manufacturer and the end-users.Therefore, both technical and economic factors have to be evaluated to develop more efficient and cost-competitive HEBs [9].
In this context, the lifetime model of a battery becomes an important tool to evaluate its degradation behaviour during vehicle operation, as the battery lifetime is reported to be typically shorter than that of supercapacitors [10].The proper assessment of battery lifetime [11] may even affect the design and definition of business models, either for vehicle manufacturers or vehicle operators, conditioning how certain economic targets could be achieved [1].Nevertheless, most publications in the literature rarely emphasize lifetime modelling when describing the sizing and operation of an ESS [12][13][14].This paper extends the analysis provided in [15] by analysing the HEB application, in which different operation constraints and ESS demands are implied.
The aim of this paper is to analyse the techno-economic influence of the battery lifetime model on the sizing and operation of a hybrid ESS.A multi-objective (MO) is defined, which covers the total costs of the battery ESS and the supercapacitor ESS and the total fuel costs.Two different battery lifetime models were considered for the targeted MO optimisation: (i) a Wöhler-curve-based lifetime model (simple and cost-efficient but less accurate) and (ii) a semi-empirical lifetime model (costly and time-consuming, but accurate) [16].Results obtained in the optimisation, which involved dual ESS sizing and operation strategy definition, were obtained with each of the two lifetime models.These results are compared in economic terms, and the ageing behaviour estimated with each model is analysed to approach how the dual ESS would perform in a real application.

Scenario Overview
The scenario analysed in this paper is based on a Series Hybrid Electric Bus (SHEB) for urban operation.Figure 1 depicts the powertrain configuration of the SHEB including a dual ESS [17].The main propulsion element is the electric motor (EM) powered from the direct current (DC) bus and mechanically connected to the vehicle transmission (T).The EM operates as generator during the regenerative braking phase.The SHEB includes a dual ESS composed of: (i) a Li-ion battery that provides high energy density; and (ii) a supercapacitor system with high power density, which is conceived to absorb or inject power peaks during regenerative braking or acceleration events, respectively.Among the different possible dual ESS connection topologies [6], a fully active or parallel-connection topology was selected.Such a topology allows for decoupling the operation of batteries and supercapacitors in order to maximise the functionalities of each energy storage technology.The genset considered includes an Internal Combustion engine (ICE) and an electric generator (EG) set.
World Electric Vehicle Journal 2018, 9, x FOR PEER REVIEW 2 of 14 efficient and economic performance.Thus, it can be noted that the system's sizing and operation are key factors to ensure the fulfillment of the application requirements, the proper harnessing of the onboard energy sources, and the defined economic objectives [7,8].In close-to-market applications, the economic issue becomes in an important decision factor as the solution has to be competitive and profitable both for the manufacturer and the end-users.Therefore, both technical and economic factors have to be evaluated to develop more efficient and cost-competitive HEBs [9].In this context, the lifetime model of a battery becomes an important tool to evaluate its degradation behaviour during vehicle operation, as the battery lifetime is reported to be typically shorter than that of supercapacitors [10].The proper assessment of battery lifetime [11] may even affect the design and definition of business models, either for vehicle manufacturers or vehicle operators, conditioning how certain economic targets could be achieved [1].Nevertheless, most publications in the literature rarely emphasize lifetime modelling when describing the sizing and operation of an ESS [12][13][14].This paper extends the analysis provided in [15] by analysing the HEB application, in which different operation constraints and ESS demands are implied.
The aim of this paper is to analyse the techno-economic influence of the battery lifetime model on the sizing and operation of a hybrid ESS.A multi-objective (MO) is defined, which covers the total costs of the battery ESS and the supercapacitor ESS and the total fuel costs.Two different battery lifetime models were considered for the targeted MO optimisation: (i) a Wöhler-curve-based lifetime model (simple and cost-efficient but less accurate) and (ii) a semi-empirical lifetime model (costly and time-consuming, but accurate) [16].Results obtained in the optimisation, which involved dual ESS sizing and operation strategy definition, were obtained with each of the two lifetime models.These results are compared in economic terms, and the ageing behaviour estimated with each model is analysed to approach how the dual ESS would perform in a real application.

Scenario Overview
The scenario analysed in this paper is based on a Series Hybrid Electric Bus (SHEB) for urban operation.Figure 1 depicts the powertrain configuration of the SHEB including a dual ESS [17].The main propulsion element is the electric motor (EM) powered from the direct current (DC) bus and mechanically connected to the vehicle transmission (T).The EM operates as generator during the regenerative braking phase.The SHEB includes a dual ESS composed of: (i) a Li-ion battery that provides high energy density; and (ii) a supercapacitor system with high power density, which is conceived to absorb or inject power peaks during regenerative braking or acceleration events, respectively.Among the different possible dual ESS connection topologies [6], a fully active or parallel-connection topology was selected.Such a topology allows for decoupling the operation of batteries and supercapacitors in order to maximise the functionalities of each energy storage technology.The genset considered includes an Internal Combustion engine (ICE) and an electric generator (EG) set.The route profile considered includes several stops (which were considered as zero-emission zones), a maximum speed of 50 km/h, and a travelling distance of around 27 km (round trip).The considered route also included a full-electric operation zone with a driving section of around 2.4 km.Thus, the considered SHEB is intended to operate both in hybrid mode (powered by the genset and the dual ESS) and in full-electric mode (only propelled by the dual ESS).The speed profile depicted in Figure 2 was repeated several times to complete a daily profile of 16 h of operation.

Aux
World Electric Vehicle Journal 2018, 9, x FOR PEER REVIEW 3 of 14 The route profile considered includes several stops (which were considered as zero-emission zones), a maximum speed of 50 km/h, and a travelling distance of around 27 km (round trip).The considered route also included a full-electric operation zone with a driving section of around 2.4 km.Thus, the considered SHEB is intended to operate both in hybrid mode (powered by the genset and the dual ESS) and in full-electric mode (only propelled by the dual ESS).The speed profile depicted in Figure 2 was repeated several times to complete a daily profile of 16 h of operation.The simulation model was developed in Matlab/Simulink.The powertrain was modelled following a backwards (or effect-cause) approach [17].The power consumed by the bus at each discrete step was calculated following a predefined speed profile by going upstream through the vehicle components [17].

Electrical Modelling
The dual ESS considered is composed of an SC pack and Li-ion BT pack connected on a parallel topology as described in the previous section.The electrical modelling was simplified to limit the computational cost [15].This was performed considering that the implemented simplified model (described below) would be capable of capturing the most relevant part of the transient voltage variation when responding to varying power demands.A more detailed electric model of the batteries would also reproduce the capacitive nature of the voltage response on the batteries.However, this would have a minimal influence on the obtained results, as the Energy Management Strategy (EMS) described below is defined in terms of the energy demands and the available energy on the ESS.
A supercapacitor model consisting of a 3000 F electric double-layer capacitor   (F) in series with a resistance  _ (Ω) was implemented as described in Figure 3.The simulation model was developed in Matlab/Simulink.The powertrain was modelled following a backwards (or effect-cause) approach [17].The power consumed by the bus at each discrete step was calculated following a predefined speed profile by going upstream through the vehicle components [17].

Electrical Modelling
The dual ESS considered is composed of an SC pack and Li-ion BT pack connected on a parallel topology as described in the previous section.The electrical modelling was simplified to limit the computational cost [15].This was performed considering that the implemented simplified model (described below) would be capable of capturing the most relevant part of the transient voltage variation when responding to varying power demands.A more detailed electric model of the batteries would also reproduce the capacitive nature of the voltage response on the batteries.However, this would have a minimal influence on the obtained results, as the Energy Management Strategy (EMS) described below is defined in terms of the energy demands and the available energy on the ESS.
A supercapacitor model consisting of a 3000 F electric double-layer capacitor C SCcell (F) in series with a resistance R int_SCcell (Ω) was implemented as described in Figure 3.
World Electric Vehicle Journal 2018, 9, x FOR PEER REVIEW 3 of 14 The route profile considered includes several stops (which were considered as zero-emission zones), a maximum speed of 50 km/h, and a travelling distance of around 27 km (round trip).The considered route also included a full-electric operation zone with a driving section of around 2.4 km.Thus, the considered SHEB is intended to operate both in hybrid mode (powered by the genset and the dual ESS) and in full-electric mode (only propelled by the dual ESS).The speed profile depicted in Figure 2 was repeated several times to complete a daily profile of 16 h of operation.The simulation model was developed in Matlab/Simulink.The powertrain was modelled following a backwards (or effect-cause) approach [17].The power consumed by the bus at each discrete step was calculated following a predefined speed profile by going upstream through the vehicle components [17].

Electrical Modelling
The dual ESS considered is composed of an SC pack and Li-ion BT pack connected on a parallel topology as described in the previous section.The electrical modelling was simplified to limit the computational cost [15].This was performed considering that the implemented simplified model (described below) would be capable of capturing the most relevant part of the transient voltage variation when responding to varying power demands.A more detailed electric model of the batteries would also reproduce the capacitive nature of the voltage response on the batteries.However, this would have a minimal influence on the obtained results, as the Energy Management Strategy (EMS) described below is defined in terms of the energy demands and the available energy on the ESS.
A supercapacitor model consisting of a 3000 F electric double-layer capacitor   (F) in series with a resistance  _ (Ω) was implemented as described in Figure 3.  Assuming that a string contains S SC supercapacitor cells in series and the SC pack groups B SC strings in parallel, the equations describing the equivalent electrical model for the SC pack were implemented as follows: being C SC (F) and R SC (Ω) the equivalent capacitance and internal resistance at SC pack level, respectively.
In the case of the battery, a model consisting of a voltage source V oc_BTcell (V) in series with a resistance R int_BTcell (Ω) was implemented as presented in Figure 3.In this case, assuming that a string contains S BT battery cells in series and the BT pack groups B BT strings in parallel, the equations for the equivalent electrical model of the BT pack were implemented as follows: being U BT (V) and R BT (Ω) the voltage and equivalent internal resistance for the BT pack, respectively.Table 1 summarises the main parameters of the battery and supercapacitor references considered in this scenario.On the one hand, the maximum C rate for charging and discharging the batteries was limited to 3.5 C to maintain the reliability of the lifetime model considered (described in the next section) and extend the battery lifetime [18].

Lifetime Modelling
With respect to the ageing performance of the supercapacitors, it was considered that these devices can typically withstand a significant amount of cycles, comparatively larger than Li-ion batteries.Therefore, a constant value of 10 6 cycles and a maximum lifetime of 10 years were implemented as the upper limit for the SC lifetime.On the contrary, the lifetime estimation of the Li-ion batteries was one of the core elements of the present study.Two different battery lifetime models were considered to evaluate their impact upon the ESS sizing and the determination of the optimal operation constraints: (i) a Wöhler-curve-based lifetime model and (ii) a semi-empirical lifetime model.These two different approaches represent different levels of complexity and experimental labour costs, but also are unequally accurate for prediction of BT lifetime under real operation conditions.

Wöhler-Curve-Based Lifetime Model
The use of Wöhler-curve-based ageing models has been typically considered in the literature, especially for sizing purposes or to perform economic analysis in applications where a battery-based ESS is integrated.
The idea behind this modelling approach consists of mathematically expressing the number of events i that can occur during the lifetime of a battery until it reaches its End-of-Life (EOL).The mathematical expression is described in Equation ( 5).In this case, i represents a certain Depth of Discharge (DOD) and the model evaluates the effect of the DOD upon the degradation of the battery, typically done using a Rainflow cycle counting algorithm [19,20].LL i represents the lifetime decrease caused by the occurrence of a certain number of i events (i.e., the number of charge-discharge full cycles at a certain DOD).
where N max E i is the maximum number of events i that the battery can withstand and NE i the number of events counted.Thereby, LL i represents the lifetime lost as consequence of the occurrence of a certain number of i events; i.e., as a consequence of a number of charge-discharge cycles at a certain DOD.Similarly, for the whole range of DOD (from 0 to 100%), the total loss of lifetime is expressed according to Equation (6): Accordingly, when the lifetime lost LL equals 1, it is considered that the cell has reached its EOL.This method allows for evaluating the corresponding loss of lifetime for a certain battery State of Charge (SOC) profile, from which the number of events at each DOD range can be defined.
Then, considering the duration of the SOC profile introduced, the total lifetime can be calculated as the inversion of LL, typically measured in years (for a one-year SOC profile): The maximum number of events NE max i at each DOD was deduced from the Wöhler curve of the batteries considered.In this case, for the considered LFP-based cell, the Wöhler curve obtained from experimental testing described in [21] was implemented as displayed in Figure 4.
World Electric Vehicle Journal 2018, 9, x FOR PEER REVIEW 5 of 14 The idea behind this modelling approach consists of mathematically expressing the number of events i that can occur during the lifetime of a battery until it reaches its End-of-Life (EOL).The mathematical expression is described in Equation (5).In this case,  represents a certain Depth of Discharge (DOD) and the model evaluates the effect of the DOD upon the degradation of the battery, typically done using a Rainflow cycle counting algorithm [19,20].  represents the lifetime decrease caused by the occurrence of a certain number of i events (i.e., the number of chargedischarge full cycles at a certain DOD).

𝐿𝐿 𝑖 =
(5) where    is the maximum number of events i that the battery can withstand and   the number of events counted.Thereby,   represents the lifetime lost as consequence of the occurrence of a certain number of i events; i.e., as a consequence of a number of charge-discharge cycles at a certain DOD.Similarly, for the whole range of DOD (from 0 to 100%), the total loss of lifetime is expressed according to Equation ( 6): Accordingly, when the lifetime lost  equals 1, it is considered that the cell has reached its EOL.This method allows for evaluating the corresponding loss of lifetime for a certain battery State of Charge (SOC) profile, from which the number of events at each DOD range can be defined.
Then, considering the duration of the SOC profile introduced, the total lifetime can be calculated as the inversion of , typically measured in years (for a one-year SOC profile): The maximum number of events    at each DOD was deduced from the Wöhler curve of the batteries considered.In this case, for the considered LFP-based cell, the Wöhler curve obtained from experimental testing described in [21] was implemented as displayed in Figure 4.

Semi-Empirical Lifetime Model
A thorough semi-empirical lifetime model was implemented in the study, comprising both cycling and calendar ageing.Such an ageing model allowed for more accurately evaluating the impact of the lifetime model upon the dual ESS sizing [21].The implemented semi-empirical

Semi-Empirical Lifetime Model
A thorough semi-empirical lifetime model was implemented in the study, comprising both cycling and calendar ageing.Such an ageing model allowed for more accurately evaluating the impact of the lifetime model upon the dual ESS sizing [21].The implemented semi-empirical battery lifetime model evaluates the capacity fade experienced over time by superimposing the effect of calendar and cycle operation.
The capacity loss due to calendar life Q loss cal (%) was calculated as described in [22] according to Equation (8).Results were obtained from an extensive testing matrix including cells stored on climatic test chambers at different ambient temperatures and different levels of SOC [22].
where T (K) and SOC (%) are the ambient temperature and the state of charge at which the cell is stored, respectively, t is the time elapsed on storage, and α 1 , α 2 , β 1 , β 2 are fitting coefficients.
The capacity loss due to cycling, Q loss cyc (%), was calculated as described in Ref. [18] according to Equations ( 9) and (10).In this case, the effect of the C rate was limited to 3.5 C and its influence on battery lifetime was neglected as no conclusive effect was concluded from the ageing study performed and described in depth in [5,8].Thereby, only the effects of the DOD and the Ampere hour throughput were considered.Further details about cycle ageing behaviour and specific modelling issues can be found in [16,18].
where Q loss cyc (%) is the equivalent capacity fade caused by cycling; DOD (%) represents the DOD at which the cycles are performed, and Ah the Ampere hours throughput during the considered cycling period.Additionally, γ 1 , γ 2 , γ 3 , α 3 , α 4 , β 3 , and β 4 are constant fitting coefficients and k represents the acceleration factor caused by static cycling ageing tests [18].Finally, the resultant total capacity loss Q loss (%) is calculated as the sum of the capacity loss caused by the calendar life Q loss cal and by the cycle life ageing Q loss cyc [16] as shown in Equation (10).
The precision of this model has been thoroughly validated under different ageing conditions based on the methodology deeply described in [16,21].For all of the cases considered in validation, the root mean square error (RMSE) prediction error calculated for the semi-empirical ageing model was below 1.4%.

Rule-Based Energy Management Strategy
In order to operate the vehicle under the scenarios of hybrid and full-electric performance, and also considering the proposed dual ESS, different operation modes were defined [23].The performance of the Rule-Based Energy Management Strategy (RB-EMS) thus depends on the instantaneous driving mode.

Hybrid Driving Mode
In this operation mode, both the genset and the dual ESS operate together.As depicted in Figure 5, the performance of the hybrid driving mode depends on the State of Charge of the battery (SOC BT ) with two defined control levels (SOC u_ctrl and SOC l_ctrl ) defining the energy mode for the SHEB operation.Thus, within the hybrid driving mode the following energy modes were defined: (i) Depleting mode (DM), (ii) Sustaining mode (SM), and (iii) Charging mode (CM).
In the first two modes, both the genset and the BT pack supply vehicle energy demands.Figure 6a,b depict the proposed strategy.The power demand P DEM (k) at each discrete step k is split between the genset and BT and SC packs and mathematically described according to Equations ( 12) and ( 13): being: P genset (k) (kW) the genset power target.P BT_dch (kW) the BT power target during discharge.P BT_ch (kW) the BT power target during charge.
V nom_BTcell (V) the nominal voltage of the BT cell.I 1C_BTcell (A) the nominal current of the BT cell.
C rate the C-rate limitation for the BT operation.p dch the ratio for BT pack discharging (p dch ∈ [0-1]) with different values for DM, SM, and Electric Mode (EM).p ch the ratio for charging the BT pack (p ch ∈ [0-1]).In the first two modes, both the genset and the BT pack supply vehicle energy demands.Figure 6a,b depict the proposed strategy.The power demand   () at each discrete step  is split between the genset and BT and SC packs and mathematically described according to Equations ( 12) and ( 13): The BT pack in depleting mode provides more energy for traction than in sustaining mode ( ℎ_ >  ℎ_ ).The genset power target   () (kW) is selected to operate between the power range  _2 ≤   () ≤  _1 , which is defined by optimisation (see Section 5).In the first two modes, both the genset and the BT pack supply vehicle energy demands.Figure 6a,b depict the proposed strategy.The power demand   () at each discrete step  is split between the genset and BT and SC packs and mathematically described according to Equations ( 12) and ( 13): The BT pack in depleting mode provides more energy for traction than in sustaining mode ( ℎ_ >  ℎ_ ).The genset power target   () (kW) is selected to operate between the power range  _2 ≤   () ≤  _1 , which is defined by optimisation (see Section 5).Regarding Figure 6: __ℎ (kW) is the maximum allowable power target for charging the SC pack power of the DC/DC). _ ,    (kWh) is the energy provided/stored in the BT pack during a traction and braking phase.
_ ,   (kWh) is the energy provided/stored in the SC pack during a traction and braking phase.The BT pack in depleting mode provides more energy for traction than in sustaining mode (p dch_DM > p dch_SM ).The genset power target P genset (k) (kW) is selected to operate between the power range P genset_2 ≤ P genset (k) ≤ P genset_1 , which is defined by optimisation (see Section 5).
Regarding Figure 6: P SC_max_ch (kW) is the maximum allowable power target for charging the SC pack (maximum power of the DC/DC).E BT_tr , E BT rg (kWh) is the energy provided/stored in the BT pack during a traction and braking phase.
E SC_tr , E SC rg (kWh) is the energy provided/stored in the SC pack during a traction and braking phase.
E CB (kWh) is the energy dissipated in the crowbar.E genset (kWh) is the energy provided by the genset.E BT_ch (kWh) is the energy provided by the genset to charge the BT pack.
In the case of the charging mode, the battery does not provide energy.The genset is in charge of providing energy both to fulfil the vehicle's demand as well as to recharge the BT pack as described by the following equation:

Full-Electric Driving Mode
In this mode, Figure 6c, only the dual ESS operates to propel the vehicle, while the genset is turned off.During a traction phase (P DEM (k) ≥ 0), the power target P BT_dch (kW) (Equation ( 11)) was configured at the maximum optimal power target for discharging (p dch = p dch_EM ).In this case, p dch_EM represents the optimal ratio to configure Equation ( 12) when the vehicle is in charging mode.The remaining energy would then be supplied by the SC pack.
In a regenerative braking phase (P DEM (k) < 0), the same power targets (P BT_ch (kW), P SC_max_ch (kW)) as in hybrid driving mode were applied.

Multi-Objective Optimisation Problem
To apply the MO optimisation to the proposed scenario, the fitness function that quantifies the score of each evaluated solution pursues finding the minimum operation costs for the SHEB as described in Equation ( 14): where C T (€⁄day) represents the total operation cost of the SHEB; BT Tcost (X), SC Tcost (X), and Fuel Tcost represent the economic models (described below) employed to estimate the operation cost of the BT pack, SC pack, and genset, respectively; and X is the vector containing the design variables in the proposed optimisation, Equation ( 15), subject to the constraints of the space of feasible solutions described in Equation ( 16).
X = [SOC u_ctrl , SOC l_ctrl , P genset_1 , P genset_2 , p dch_DM , p dch_SM , p dch_EM , p ch_CM , B BT , B SC World Electric Vehicle Journal 2018, 9, 28 9 of 13 In this scenario, the variables S BT and S SC were assumed to be constant in order to reach a battery pack voltage of 600 V and a supercapacitor voltage of 390 V. Thus, the dual ESS sizing optimisation will only consider the variation of energy capacity (number of branches: B SC , B BT ).The MO problem is solved by means of genetic algorithms (GA).The iterative process carried out by the GA is described in detail in [24].

Operation Cost of the Dual ESS
The total operation cost of the dual ESS was defined according to Equation ( 17) [15].As can be observed, the term u refers to either the battery or the supercapacitor pack.Thus, the cost components corresponding to the battery and supercapacitor pack costs in Equation ( 14) were calculated equally according to the following equations.
being u M y (/year) the annualised maintenance cost of the u pack.The annualised capital cost for the u pack u Ca y (/year) was defined as follows: being C kW dcdc (€/kW) the referential cost of the DC/DC converter; P dcdc u (kW) the power rate of the DC/DC converter; C kWh u (€/kWh) the referential cost of the considered u energy storage technology; Ca u (kWh) the sizing of the u pack; and I (%) and T (years) the banking interest rate and the lifetime of the whole system (SHEB), respectively.The annualised replacement cost (cycling cost) of the u pack, u Re y (€/year), was defined as follows: where ceil(x) rounds its argument x to the higher integer value; and Li f e u (years) is the lifetime estimation for the u pack provided by the corresponding lifetime model (as described in Section 3.2).

Results and Discussion
The optimisation of the operation constraints and the sizing of the dual ESS were performed considering the two different battery lifetime models described in Section 3.2.In the two cases, all the remaining elements of the modelled SHEB were unaltered, including the EMS, route, dual ESS electrical models, and the lifetime model employed to estimate the lifespan of the supercapacitor pack.
Figure 7 shows the three solutions selected from the Pareto front obtained in the optimisation with (a) the Wöhler-curve-based and (b) the semi-empirical lifetime models.The solution number 1 was selected as the most suitable solution for each of the two cases.Table 2 summarises the optimal set of values for the design variables assigned for the optimisations performed with each lifetime model.
In the case of the optimal solution calculated with the Wöhler-curve-based lifetime model, the dual ESS would be composed of a battery pack of 16.56 kWh (182S12P) and a supercapacitor pack of 0.86 kWh (144S2P).The operating costs of the SHEB would be ca.170 €/day, being the daily costs of the batteries and supercapacitors 37.2 and 30.6 €/day, respectively.For the optimal solution calculated with the semi-empirical lifetime model, the selected sizing for the dual ESS would be composed of a battery pack of 13.80 kWh (182S10P) and a supercapacitor pack of 0.86 kWh (144S2P).In that case, the operating costs of the SHEB would be ca.183 €/day, being the daily costs of the batteries and supercapacitors 19.8 and 30.9 €/day, respectively.In the case of the optimal solution calculated with the Wöhler-curve-based lifetime model, the dual ESS would be composed of a battery pack of 16.56 kWh (182S12P) and a supercapacitor pack of 0.86 kWh (144S2P).The operating costs of the SHEB would be ca.170 €/day, being the daily costs of the batteries and supercapacitors 37.2 and 30.6 €/day, respectively.For the optimal solution calculated with the semi-empirical lifetime model, the selected sizing for the dual ESS would be composed of a battery pack of 13.80 kWh (182S10P) and a supercapacitor pack of 0.86 kWh (144S2P).In that case, the operating costs of the SHEB would be ca.183 €/day, being the daily costs of the batteries and supercapacitors 19.8 and 30.9 €/day, respectively.When the semi-empirical lifetime model was employed, the optimisation tended towards a smaller battery sizing.The lifetime of different battery and supercapacitor sizing solutions was evaluated in order to further explore the reasons behind such differences for the optimal case selection.The optimal set of constraints calculated for each selected optimal solution (one for each of the two lifetime models) was maintained while varying only the sizing of the dual ESS.Thereby, Figure 8 shows the variation in the battery lifetime estimation (for each of the two lifetime models) considering the operating conditions defined for each of the two selected optimal solutions.The obtained lifetime values were scaled on a p.u. basis, considering as reference value the lifetime estimated with the Wöhler-curve-based lifetime model (the solution with the greatest lifetime in the leftmost plot in Figure 8 is 1 p.u.).
The results in Figures 7 and 8 revealed a significant influence of the selected lifetime model on the selection of the optimal solution.Indeed, the semi-empirical lifetime model provided significantly larger lifetime estimations.This allowed the optimisation to increase the operation constraints while maintaining a smaller dual ESS sizing compared to the optimal case calculated with the Wöhler-curve-based lifetime model.
As can be observed in Figure 9, the batteries are subjected to a repetitive SOC pattern in which deep BT charges and discharges dominate their degradation.In Figure 9b, it can be observed that in the case of the optimisation performed with the semi-empirical model a better exploitation of the batteries was imposed even when accounting for a smaller dual ESS sizing.On the contrary, in the case of the optimisation performed with the Wöhler-curve-based model, the optimal set of operation constraints forced the battery to operate in a saturated status in which the battery remained most of the time at low SOC levels.Moreover, several micro-cycles were performed in that region.The typical shape of the Wöhler curves provides a large amount of cycles at low DOD ranges, driving the optimisation to that zone of battery operation.However, experimental battery ageing labours at low DOD values are costly and time-consuming, and thus the reliability of any lifetime model in such an area is usually lower.smaller battery sizing.The lifetime of different battery and supercapacitor sizing solutions was evaluated in order to further explore the reasons behind such differences for the optimal case selection.The optimal set of constraints calculated for each selected optimal solution (one for each of the two lifetime models) was maintained while varying only the sizing of the dual ESS.Thereby, Figure 8 shows the variation in the battery lifetime estimation (for each of the two lifetime models) considering the operating conditions defined for each of the two selected optimal solutions.The obtained lifetime values were scaled on a p.u. basis, considering as reference value the lifetime estimated with the Wöhler-curve-based lifetime model (the solution with the greatest lifetime in the leftmost plot in Figure 8 is 1 p.u.).
The results in Figures 7 and 8 revealed a significant influence of the selected lifetime model on the selection of the optimal solution.Indeed, the semi-empirical lifetime model provided significantly larger lifetime estimations.This allowed the optimisation to increase the operation constraints while maintaining a smaller dual ESS sizing compared to the optimal case calculated with the Wöhler-curve-based lifetime model.As can be observed in Figure 9, the batteries are subjected to a repetitive SOC pattern in which deep BT charges and discharges dominate their degradation.In Figure 9b, it can be observed that in the case of the optimisation performed with the semi-empirical model a better exploitation of the batteries was imposed even when accounting for a smaller dual ESS sizing.On the contrary, in the case of the optimisation performed with the Wöhler-curve-based model, the optimal set of operation constraints forced the battery to operate in a saturated status in which the battery remained most of the time at low SOC levels.Moreover, several micro-cycles were performed in that region.The typical shape of the Wöhler curves provides a large amount of cycles at low DOD ranges, driving the optimisation to that zone of battery operation.However, experimental battery ageing labours at low DOD values are costly and time-consuming, and thus the reliability of any lifetime model in such an area is usually lower.On the contrary, the used semi-empirical lifetime model was validated under several operating conditions, including micro-cycles [16].Thus, the lifetime estimations obtained in the low DOD operation area are more accurate than those obtained with the Wöhler-curve-based lifetime model, and the optimisation was driven to other areas of battery operation in which an extended battery lifetime was achieved.
From the figures presented in this section, a significant impact of the lifetime model could be observed in the obtained optimisation results.Thereby, simplifications on the battery lifetime estimation, e.g., when using the widely referenced Wöhler-curve-based lifetime model, may yield inconvenient decisions on the battery sizing and operation strategy.This would certainly have an impact upon the system cost and performance, but most importantly it may even influence the business model definition and the cost structure and revenue streams.

Conclusions
In this paper, a thorough analysis of the influence of the battery lifetime model upon the optimisation results for ESS sizing and operation constraints was presented considering a SHEB as the case study.
The obtained results suggest a significant influence of the lifetime model's accuracy on the costs associated with both the operation of the hybrid bus and the sizing of the BT and SC.A On the contrary, the used semi-empirical lifetime model was validated under several operating conditions, including micro-cycles [16].Thus, the lifetime estimations obtained in the low DOD operation area are more accurate than those obtained with the Wöhler-curve-based lifetime model, and the optimisation was driven to other areas of battery operation in which an extended battery lifetime was achieved.
From the figures presented in this section, a significant impact of the lifetime model could be observed in the obtained optimisation results.Thereby, simplifications on the battery lifetime estimation, e.g., when using the widely referenced Wöhler-curve-based lifetime model, may yield inconvenient decisions on the battery sizing and operation strategy.This would certainly have an impact upon the system cost and performance, but most importantly it may even influence the business model definition and the cost structure and revenue streams.

Conclusions
In this paper, a thorough analysis of the influence of the battery lifetime model upon the optimisation results for ESS sizing and operation constraints was presented considering a SHEB as the case study.
The obtained results suggest a significant influence of the lifetime model's accuracy on the costs associated with both the operation of the hybrid bus and the sizing of the BT and SC.A difference of up to ca. 8% on SHEB daily operation costs was estimated, depending on the employed lifetime model (170 €/day for Wöhler-based optimisation and 183 €/day for semi-empirical-based optimisation).Similarly, when only dual ESS costs were considered, a difference of ca.25% was estimated depending on the employed lifetime model (67.8 €/day for Wöhler-based optimisation and 50.7 €/day for the semi-empirical-based optimisation).
According to the obtained results, and considering the strong impact of the ESS upon total system costs, a detailed and accurate evaluation of battery lifetime is crucial.Therefore, the use of an accurate battery lifetime model is advisable to gain full control of the system and its cost structures.

Figure 4 .
Figure 4. Number of cycles versus Depth of Discharge (DOD) curve considered for the Whöler curve [21].

Figure 4 .
Figure 4. Number of cycles versus Depth of Discharge (DOD) curve considered for the Whöler curve [21].
Depleting and Sustaining modeWorld Electric Vehicle Journal 2018, 9, x FOR PEER REVIEW 8 Full-electric driving mode

Figure 8 .
Figure 8. Lifetime estimated for different battery and supercapacitor sizing values with the set of optimal operation constraints obtained with each of the two lifetime models considered.The lifetime is scaled in p.u. considering the maximum lifetime obtained for the Wöhler-based optimisation.

Figure 8 .Figure 9 .
Figure 8. Lifetime estimated for different battery and supercapacitor sizing values with the set of optimal operation constraints obtained with each of the two lifetime models considered.The lifetime is scaled in p.u. considering the maximum lifetime obtained for the Wöhler-based optimisation.World Electric Vehicle Journal 2018, 9, x FOR PEER REVIEW 12 of 14

Figure 9 .
Figure 9.Comparison of the daily battery operation profile for the optimal solutions calculated: (a) with the Wöhler-curve-based and (b) semi-empirical lifetime models.
E CB (kWh) the energy dissipated in the crowbar.
BT_min SOC Figure 5. Energy Management Strategy (energy modes) in hybrid driving mode.
_ (V) the nominal voltage of the BT cell. 1_ (A) the nominal current of the BT cell.  the C-rate limitation for the BT operation. ℎ the ratio for BT pack discharging ( ℎ ∈ [0-1]) with different values for DM, SM, and Electric Mode (EM). ℎ the ratio for charging the BT pack ( ℎ ∈ [0-1]).  (kWh) the energy dissipated in the crowbar.
() (kW) the genset power target. _ℎ (kW) the BT power target during discharge. _ℎ (kW) the BT power target during charge.
() (kW) the genset power target. _ℎ (kW) the BT power target during discharge. _ℎ (kW) the BT power target during charge. _ (V) the nominal voltage of the BT cell. 1_ (A) the nominal current of the BT cell.  the C-rate limitation for the BT operation. ℎ the ratio for BT pack discharging ( ℎ ∈ [0-1]) with different values for DM, SM, and Electric Mode (EM). ℎ the ratio for charging the BT pack ( ℎ ∈ [0-1]).  (kWh) the energy dissipated in the crowbar.

Table 2 .
Optimal solution selected for the two optimisation cases.