#### Abstract

This paper presents a methodology for the sizing of a heavy-duty fuel cell commercial vehicle. The parameters scanning model and the long-term stochastic drive cycle are adopted for this proposed sizing framework. The dynamic programming algorithm is employed as the energy management strategy to assess the performance of sizing. The efficacy of this framework is evaluated, and a detailed analysis for the hydrogen consumption is given in the results. In addition, a prediction analysis based on the price performance of the next decade is also given in this work.

#### 1. Introduction

The energy and environment are two limitations for the transportation industry [1, 2]. The heavy-duty commercial vehicle is imperative to turn to new energy for its high fuel consumption and emission. As a fuel-renewable and environment-friendly on-board propulsion system, hydrogen fuel cells are of great significance for alleviating the current energy and environment dilemma.

Restricted by the dynamic properties of fuel cell, fuel cell vehicles (FCVs) are always equipped with at least two energy sources: a fuel cell stack and a battery pack or a supercapacitor, which are also known as fuel cell hybrid electric vehicles (FCHEVs) [3]. Therefore, FCHEV is generally treated as a particular electric vehicle which equipped with two or three power sources [4]. Many approaches are presented in the literature to optimize the vehicle operating cost, which can be divided into two levels: the optimization level and the design level. The optimal operation aims to find the optimal operating points for fuel cell to achieve the minimum hydrogen cost, which is determined by the energy management strategy (EMS). The optimal design aims to determine the optimal sizes of vehicle components [5, 6].

The key role of EMS for a FCHEV is to decide the power splitting for energy sources at each instant, while at the same time, the drivability constraints of vehicle must be satisfied [7]. There have been many researches on EMS to improve the fuel economy. In general, rule-based EMS and optimization-based EMS are two approaches widely studied [8, 9]. Under current technology, rule-based EMS has been widely used in real vehicles’ controller for its advantages in real time, simplicity, and cost. However, subjectivity and uncertainty of this strategy may lead the control rules far to the optimal rules [8]. On the other hand, thanks to the increasing computing power of on-board processors, more complex EMS has the potential to be applied in real vehicles; optimization-based EMS has become the focus of current academic research. Optimization-based EMS is often represented as optimization problems for constrained systems in finite-time domain; the common optimization algorithms include dynamic programming (DP) [10], Pontryagin’s minimum principle (PMP) [11], and genetic algorithm (GA) [12]. In addition, two optimization frameworks (equivalent consumption minimization strategy (ECMS) [13] and model prediction control (MPC) [14–16]) are also applied to this strategy to improve the control effect.

In addition to the EMS, many researchers focused their attention on the sizing of energy sources to improve vehicle performance. The main purpose of sizing of energy sources is to optimize the dynamic and fuel economy with fewer cost; reasonable parameters can maximize the potential of a vehicle’s EMS [17]. With the development of electric vehicle, numbers of algorithms have been adopted to solve this problem. As early as 2004, Lukic and Emadi [18] verified the influence of drivetrain hybridization on energy consumption and dynamic performance. In their work, the hybridization factor was defined, and the EMS was implemented as charge sustaining. Kim and Peng [19] proposed a combined optimization method of EMS and parameters sizing for FCHEV. Different subsystem-scaling models were adopted to predict the characteristics of components on different sizes, and a near-optimal EMS was designed with the inspiration of stochastic dynamic programming results. In [20], the scholars constructed a stochastic drive cycle to simulate the randomness of real-world driving conditions, and ECMS was used as an online EMS; the equivalent factors of ECMS were adjusted with the proposed sizing methodology. In [21], a parameters’ sizing methodology based on ordinal optimization and dynamic programming was proposed. Hu et al. [22] gave a compared analysis for two optimization-based EMS, and the influence of downsizing the battery was also studied. In their further study, the influence of driving pattern on sizing was studied in [23]; the comparison of three different electrochemical energy storage systems and the sizing framework for hybrid electric vehicle are also given [24, 25]. In addition, Karaoğlan et al. [26] studied the influence of gear ratio on fuel economy and emissions.

For the EMS, optimization-based strategy has been the focus of current research; and for the parameters sizing, systematic algorithms are used to find the optimal/near-optimal sizes of components for FCHEV; rule-based EMS or optimization-based EMS are adopted as an associated problem for these algorithms. Although existing research has proposed a number of approaches for parameters sizing, these approaches are based on the few drive cycles or only one standard drive cycle, and parameters’ sizing for hybrid heavy-duty commercial vehicles is little studied.

In this paper, a sizing methodology for a fuel-cell/battery commercial vehicle is presented. It is based on the DP-based EMS and parameters scanning model, and a long-term stochastic drive cycle is adopted for the sizing framework. The main contributions are as follows: (a) A long-term stochastic drive cycle that is based on the cluster and stochastic procedure is adopted, and the processes of generate stochastic drive cycle are also available. (b) A detailed analysis on oversizing the fuel cell stack and battery is given. (c) A predictive analysis is also given with considering the price performance in the next decade.

The remainder of this paper is organized as follows. In Section 2, the model of a heavy-duty fuel cell commercial vehicle is established. In Section 3, the design process of the long-term stochastic drive cycle is proposed. In Section 4, the DP-based EMS is introduced. In Section 5, the sizing methodology of propulsion system is presented. Section 6, provides results and discussion, and in Section 7, the conclusion of this paper is given.

#### 2. Mathematical Model of an FCHEV

In general, for an FCHEV, fuel cell stack is employed as the main power source, and battery is adopted as an energy storage system, which is also known as the auxiliary source. The structure of selected FCHEV is shown in Figure 1. The electric motor is the only conversion equipment of electric energy and kinetic energy, fuel cell stack and battery pack composed of the propulsion system to provide power to permanent magnet synchronous motor and the motor’s inverter (DC-AC), and a DC-DC converter is used to boost the voltage from fuel cell stack to DC bus. In this work, the powertrain system is divided into a combined propulsion subsystem and a transmission subsystem. The parameters sizing of propulsion system is the main work of this paper.

##### 2.1. Vehicle and Drivetrain

The mathematical model building of FCHEV is reversed with the target of tracking drive cycle, and the impacts of lateral dynamics and rotating mass are ignored. At every moment of vehicle’s operation, the torque and speed of wheels are given bywhere is the rolling radius of the wheel, the FCHEV’s mass with full load; the vehicle’s velocity at the current moment; the air density; the frontal area; the aerodynamic drag’s coefficient; the rolling resistance coefficient; the road’s inclination angle. The motor’s torque , speed , and power are formulated as follows:where and are the efficiency and gear ratio of the final drive, respectively.

For this work, the motor power is defined as positive during traction while negative during breaking. The propulsion system consists of fuel cell stack and battery pack outputs the power demand of the motor, as shown in the following formula: where and are the power of fuel cell stack and battery pack. Parameters of the heavy-duty fuel cell commercial vehicle are listed in Table 1.

##### 2.2. Fuel Cell Stack

Fuel cell is an electrical device that converts chemical energy into electricity. Proton exchange membrane fuel cell (PEMFC) is the most common on-board propulsion system for fuel cell vehicles. A complete on-board PEMFC system includes fuel cell stack, hydrogen storage system, hydrogen circuit, air circuit, water circuit, and coolant circuit [3]. The models for these systems are very complicated. This work is mainly concerned with the power-split of the bus but not the detailed conversion process of fuel cells. Therefore, a simple efficiency graph model is used for the simulation.

In this selected model, hydrogen consumption and efficiency of the fuel cell stack are both formulated as the function of power demand; these datasets are obtained through testing, which have been provided in ADVISOR (FC_ANL50H2) [27]. Scaling models are adopted to facilitate parameters’ sizing, as shown in Figure 2.

**(a)**

**(b)**

##### 2.3. Battery Pack

Benefiting from the high power density and high reliability, lithium-ion batteries have become the most commonly used energy storage system for FCHEV. The purpose of adding a battery pack into the propulsion system is to improve the drivability and efficiency performance of the vehicle.

For the model of battery pack, the circuit model is used for simulation. In this model, the battery is equivalent to a circuit with an ideal voltage source and an internal resistance. When the battery pack is operating, its current and power are formulated as follows:where is the open circuit voltage, the current, the internal resistance, and the power of battery pack. These parameters can be obtained from equations (5) to (7).where is the conversion efficiency for charge/discharge, the series number of cells, and the parallel number of cells. The state of charge (SOC) of battery at time step is defined aswhere is the capacity of battery pack.

and are the functions of SOC, as shown in Figure 3.

**(a)**

**(b)**

It is worth mentioning that the battery pack is composed of strings connected in parallel, with each string containing the same number of cells in series. More detailed descriptions about configuration and modelling of battery pack are available in [23].

##### 2.4. Drive Motor

For an FCHEV, the drive motor’s function mainly includes the following two parts: (a) converting the electrical energy from bus to mechanical energy to drive the vehicle and (b) recycling the vehicle’s braking power to charge the battery. A permanent magnet synchronous motor (PMSM) is adopted for its compact and high efficiency. The model of drive motor is based on its efficiency map, as shown in Figure 4.

#### 3. Long-Term Stochastic Drive Cycle

The drive cycle refers to the driving characteristics of a certain type of vehicle (passenger car, bus, and cargo vehicle, etc.) in a specific condition (urban, suburban, highway, etc.); it is generally a set of points representing speed versus time. There have been many systems of typical driving cycles around the world.

Although the representativeness of cycles will be considered in the design process of typical drive cycle, the optimization results obtained under few drive cycles cannot satisfy the randomness of real driving conditions [20]. Therefore, a long-term stochastic drive cycle is built to overcome this drawback. The constructed drive cycle is based on the typical drive cycles of heavy commercial vehicles, and cluster analysis and random process are adopted. The flowchart of this procedure is shown in Figure 5.

In this work, eight typical drive cycles for commercial vehicles are introduced as the DC library, which are numbered and shown in Figure 6. Twelve characteristic parameters are used for clustering, and the corresponding descriptions of these parameters are listed in Table 2. The characteristic parameters for each drive cycle are calculated and listed in Table 3. These parameters will be the raw data of clustering.

Hierarchical clustering and Euclidean distance are adopted to obtain the similar degree of the different drive cycles. The clustering results are shown in Figure 7. The library is divided into 5 classes with the benchmark of similarity to 1: the first category includes drive cycles 1, 4, 6, and 8, and drive cycles 2, 7, 3, and 5 are classified as the second to fifth categories, respectively. The drive cycles 8 (CHTC_HT), 2 (WVUCITY), 7 (NYCOMP), 3 (WVUINTER), and 5 (UDDSSHCV) are selected as the representatives of each category. These five drive cycles will be regarded as the originals in the long-term stochastic drive cycle.

The next step for this work is to build the combined drive cycle, which is generated from the original drive cycles. After that, the combined drive cycle will be divided into microtrips and microidles. A microtrip is a segment of the drive cycle where the velocity is not equal to zero, and a microidle is the segment where the velocity is equal to zero. The relationship between microtrips and microidles is depicted in Figure 8. Then, the divided segments are added to the stochastic drive cycle randomly. In this process, the microtrips and microidles are scaled with a stochastic procedure. The scaling processes for microtrips are for their velocity, and for microidles, it is for their length of time, the scaling factors are both defined as a stochastic number from 0.8 to 1.2. To follow the objective laws for drive cycles, the microtrip and microidle should be added in turn, which means after adding a microtrip, a microidle must be added as its follow-up. In addition, the long-term stochastic drive cycle can more effectively reflect the randomness of the drive cycle, the total length of stochastic drive cycle is set as 10 times that in the combined drive cycle in this work.

Another noteworthy point for this process is that after the scaling of a microtrip, three indexes should be detected: the maximum speed, maximum acceleration, and maximum deceleration. If any index is higher than the combined drive cycle, the scaling process should be repeated until satisfied.

After completing the construction of long-term stochastic drive cycle, evaluation indexes are introduced to evaluate the rationality. Figure 9 illustrates the distribution of two drive cycles; Table 4 depicts the comparison of characteristics of combined drive cycle and long-term stochastic drive cycle. The above charts show that the long-term stochastic drive can accurately reflect the characteristics of original cycles.

**(a)**

**(b)**

**(c)**

**(d)**

#### 4. DP-Based EMS

In general, the EMS for FCHEV has a great impact on efficiency. In this work, an optimal EMS based on dynamic programming algorithm is adopted to assess the economic potential of sizing.

Dynamic programming is a numerical method for solving multistage decision-making problems. It can provide the optimal solution for problems of any complexity level in the limits of computational capabilities [7].

For a controlled multistage decision-making problem, the state transfer function can be described aswhere is the state variable at time and is the control variable. For a control policy , the cost in time step is defined aswhere is the instantaneous cost function. The optimal cost is

The optimal control policy .

For the DP-based EMS for FCHEV, the SOC of battery is selected as the state variable, and the output power of fuel cell stack is adopted as the control variable. The state transfer function is defined as follows:where is the change rate of state in time , which can be obtained with equations (7) to (11).

For the structure of non-plug-in FCHEV, a constraint of terminal SOC is required to sustain the energy storage system. In principle, the constraint of terminal SOC can be taken into account in two different ways, namely, as a soft constraint or a hard constraint. In this study, a hard constraint is adopted by requiring that the energy storage stored at the end of the mission equal to the value at the start of the mission, as shown in equation (13). A more detailed explanation will be found in [28–30]. Other parameters under constraints are listed in equation (14). In particular, the and are set to 0.4 and 0.8, respectively.

#### 5. Structure of Parameters’ Sizing

The objective of sizing is to obtain the near-optimal sizes of maximum power of fuel cell stack and the number of battery cell . The sizes of propulsion should satisfy the following objectives: (a) meet the drivability and (b) minimize the operating cost of the vehicle during its life cycle.

For drivability, passenger cars typically have higher velocity and acceleration targets, while for commercial vehicles, equipped with large load mass, the vehicle’s climbing performance always receives more attention. Therefore, the sizing method is tested with three subobjects: (a) The peak power of fuel cell stack alone must be able to sustain the maximum velocity; (b) the peak power with fuel cell stack and battery pack together must meets the maximum climbing slope index at a given velocity; (c) the power with fuel cell stack and battery pack together should fulfil the power demand of long-term stochastic drive cycle. The first subobject is used to obtain the low limiting value of , and the others are used to determine . It should be noted that the for the subobject (c), the verification results are related to the selected EMS, DP-based strategy is employed for this work, which has been introduced in Section 4.

This sizing methodology is based on the generated long-term stochastic drive cycle, a parameters’ scanning model is added for sizing, and the dynamic programming algorithm is adopted for the EMS. The flowchart of this proposed methodology is shown in Figure 10.

In this sizing methodology, the range of fuel cell stack’s peak power is divided into independent nodes; the range of battery cell’s number is divided into independent notes; thses two groups of nodes constitute a crossover model with crossover points, and for each crossover point, there is an FCHEV model corresponding to it. The optimal cost for each size is obtained by DP-based EMS. In addition, the low limiting value of should meet the power demand of the vehicle at a maximum velocity (equation (15)), and the peak-powers provided by the battery and fuel cell should meet the requirements of the vehicle’s climbing performance (equation (16)). The peak power of the battery pack can be expressed as the product of its capacity and the maximum discharge rate equation (17):where is the power demand for maximum velocity, is the power demand for climbing, is the average efficiency of motor, is the charge/discharge efficiency of battery, and and are the efficiency of DC-AC inverter and DC-DC converter, respectively.

#### 6. Results and Discussion

This simulation is carried out with MATLAB; the interval is set to 10 kW and 25 for the fuel cell stack’s peak power and the battery cells’ number, respectively.

Under the current technical, lithium-ion batteries can maintain a discharge time of 5 s at the peak discharge rate of 30 C [31]. However, high rate discharge will seriously affect the life of the battery. In this paper, the maximum discharge rate of the battery is set as 20 C. In addition, the capacity of the battery cell used in the simulation is 3 Ah. Therefore, on the basis of equations (18) to (20), the minimum constraints of and will be estimated, as shown in Figure 11.

Figure 12 shows the hydrogen consumption of vehicles with different sizes under long-term stochastic drive cycle. For a fair comparison, all the hydrogen costs are computed with the SOC deviation. Two conclusions can be obtained according to this figure:(1)Starting from kW, the hydrogen consumption decreases with the increasing of the peak power of fuel cell stack, but it increases when the peak power exceeds 120 kW(2)For a determined , the hydrogen consumption decreases with the increasing of the number of battery cells, but the slope of decline gradually decreases

To understand the causes of the changing of hydrogen consumption, we consider three relevant performance indicators: the global propulsion efficiency, the braking-recovery-energy to negative-load-power ratio, and the global fuel cell stack efficiency.

The global propulsion efficiency is formulated as follows:where is the positive vehicle power, is the length of stochastic drive cycle, is the hydrogen consumption, and is the simulation step.

The braking-recovery-energy to negative-load-power ratio is defined as follows:where is the negative vehicle power and is the braking energy recovered by battery pack, which is formulated as follows:

The global fuel cell stack efficiency is defined as follows:

To make the figure clearer, four types of fuel cell stack peak power are selected for this work: The rest results of other types are visible in Figure 13.

**(a)**

**(b)**

**(c)**

Figure 14 shows the effect of and on defined relevant performance indicators. Figure 14(a) shows the maximum global propulsion efficiency is approximately 55.56%, which is reached when kW and . Increasing or decreasing the type of fuel cell will worsen this indicator. On the other hand, increasing the number of battery cells can improve the global efficiency, but its growth trend gradually slows down and eventually to level. Figure 14(b) shows the size of fuel cell stack has less affected in braking-recovery-energy, while the number of battery cells has a great impact on this indicator: when , limited by the physical constraints (e.g., the maximum SOC or maximum voltage), the battery pack is not enough to recover all the braking energy. This means that friction brakes will dissipate power from fuel cell stack, which is not friendly to the global efficiency. Figure 14(c) shows that the global fuel cell stack efficiency performs best when kW, and the increase in will benefit this indicator. Similarly, the growth trend gradually slows down with the increase of .

**(a)**

**(b)**

**(c)**

The results in Figure 10 clarifies that the reduction of hydrogen can be achieved with an appropriate , and the increase of will also benefits the fuel economy. It is beneficial to reduce the running cost, but the sharp increase in vehicle’s price led by the increase of fuel cell and battery is unacceptable for both manufactures and users. Vehicle’s price should also be considered as a further indicator.

To balance the cost of hydrogen consumption and component costs, a cost function which is defined as the cost per 100 km of the propulsion system is proposed, as shown in equation (18).

In this function, the vehicle’s cost is divided into fuel cell stack cost , battery cost , and hydrogen cost . The fuel cell stack cost and the battery cost are cycle-normalized cost, as expressed in equations (19) and (20) [28, 29]. In addition, the cost of hydrogen is the value of hydrogen consumed when the vehicle’s total driving distance reaches 100 km in the constructed long-term stochastic drive cycle, as shown in equation (21).where is the length of cycle, is the distance under a stochastic drive cycle, is the average travelled distance of the vehicle per year, is the vehicle lifetime, and is the yearly interest rate. Here, it is assumed that there are no components replaced during three-year operation, and the corresponding travelled distance is fifty thousand kilometres per year. These parameters above are given in Table 5. is the fuel cell stack price per kilowatt, is the battery price per kilowatt hour, and is the hydrogen price per gram. Three price conditions on the years 2020, 2025, and 2030 [32] are chosen for this study, as shown in Table 6.

Figure 15 shows the price impact on vehicle’s total cost. Figure 14(a) shows the vehicle total cost with the price in 2020. In this case, the unit prices of fuel cells and battery cell are higher; the optimal cost is 25.07 US dollars per 100 km with the corresponding size to = 70 kW and = 696. The cost previews of 2025 and 2030 are shown in Figure 14(b) and 14(c) respectively. On the preview of 2025, benefiting from the production increases and technology improves, the unit price of fuel cell and battery will decrease rapidly, and the optimal cost will drop to 21.46 US dollars per 100 km with the corresponding size to = 80 kW and = 640. In 2030, with the further decrease of component prices, the optimal cost will drop to 20.41 US dollars per 100 km. The optimal sizes and costs for different price conditions are summarized in Table 7.

**(a)**

**(b)**

**(c)**

Figure 16 gives the comparative results of optimal costs under three price conditions. Based on the price trends in future, the optimal cost gradually decreases from 2020 to 2030. Compared to 2020, the decreases have reached 14.4% and 18.59% in 2025 and 2030, respectively. The main reason is the sharp decrease of fuel cell and battery price. At the same time, benefiting from the decrease of fuel cell and battery price, the proportion of lifetime cost decreases rapidly. For example, based on the price performance in 2020, the proportion of lifetime cost is 25.56%, and it decreases to 18.26% and 14.35% in 2025 and 2030, respectively. Correspondingly, the hydrogen costs are 74.44%, 81.74%, and 85.65% for the price performance in 2020, 2025, and 2030, respectively. Along with the decreasing of hydrogen and components price, the hydrogen cost will be the more important performance index for vehicle cost in the near future.

#### 7. Conclusion

This paper presents a sizing methodology for a heavy-duty fuel cell commercial vehicle on the basis of long-term stochastic drive cycle, parameters’ scanning model, and DP-based EMS. The causes of hydrogen consumption under different sizes are studied in detail, and a discussion of the vehicle cost in the next decade is also given, on the basis of propulsion system and hydrogen price performance.

The main findings are as follows:(1)The hydrogen consumption will be minimized when the peak power of fuel cell stack reaches 120 kW, and the global propulsion efficiency will be improving with the battery cell number increase.(2)When the number of battery cells is small, the global propulsion efficiency is affected by the braking-recovery-energy and the fuel cell stack efficiency, while as the number of battery cells continues to increase, the main cause for that is the fuel cell stack efficiency. Limited by the maximum efficiency of fuel cell stack, the effect of increasing the battery cells number on the global propulsion efficiency will gradually decrease.(3)The size of = 70 kW and = 696 will be a better choice for the vehicle total cost, under the consideration of the cost performance in 2020. The predictive analysis shows that the larger fuel cell stack will be suitable for the FCHEV with the decrease of fuel cell and battery price.

#### Data Availability

The datasets and codes of this paper for the simulation are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this study.

#### Authors’ Contributions

Weiguang Zheng carried out conceptualization, formal analysis, writing, review, and editing. Weiwei Xin was responsible for software and wrote the original draft. Enyong Xu carried out conceptualization and data curation. Shuilong He carried out conceptualization and formal analysis. Jirong Qin carried out supervision and project administration. Heng Wang carried out methodology and validation. All authors have read and agreed to the published version of the manuscript.

#### Acknowledgments

This study was supported by Innovation-Driven Development Special Fund Project of Guangxi (Guike AA18242033), Liuzhou Science Research and Planning Development Project (nos. 2020GAAA0403 and 2019AD10203), Liudong Science and Technology Project (20200108), and Innovation Project of GUET Graduate Education (2019YCXS008).