Energy Management Optimization of Series Hybrid Electric Bus Using an Ultra-Capacitor and Novel E ﬃ ciency Improvement Factors

: The existing series hybrid electric bus (SHEB) uses an ultra-capacitor (UC) to extend battery life, mitigate vehicle weight, and reduce cost. However, previous studies did not clearly identify the operation timing and load of the UC for e ﬃ ciency improvement in an SHEB. This paper proposes novel e ﬃ ciency improvement factors, with their application criteria for the ideal operation timing and load of the UC in an SHEB. The factors are the threshold of the required power of the motor (TRPM), slope of the power split ratio (SPSR), and y-axis intercept of the power split ratio (YPSR). The TRPM determines the duration of using just the battery. The SPSR or YPSR determine the most e ﬃ cient load ratio between the battery and UC. The criteria for using them are set using particle swarm optimization. Manhattan, Braunschweig, and Orange County driving cycles were used to reﬂect various road load conditions. The results showed that the proposed factors and their setting criteria guarantee a signiﬁcant reduction in the fuel consumption and more energy-e ﬃ cient SHEBs. Author Contributions: Conceptualization, S.L., M.K.; Methodology, S.L., M.K.; Data curation, K.L., G.H., M.K.; Formal analysis, J.K., G.H.; Investigation, K.-J.L., M.K.; Validation, G.H., K.L.; Visualization, K.-J.L., G.H., M.K.;


Introduction
Global issues, such as the exhaust gas generated by internal combustion engines (ICEs), have a significant impact on environmental pollution and the human body. Although developers of electric vehicles (EVs) are focusing on becoming emission-free, they are taking considerable time to address the vehicles' disadvantages, such as the construction of charging infrastructure, limited driving distance, and long charging/waiting time. To overcome these barriers, hybrid electric vehicles (HEVs) offer realistic and practical solutions and have been studied worldwide. HEVs use multiple energy sources and are divided into three categories according to their structures and characteristics: series, parallel, and series-parallel types. They usually use motors, generators, batteries, and ICEs to operate the automobile, meeting the energy requirements of each component. In particular, a series hybrid electric vehicle (SHEV) can operate the engine most efficiently, regardless of the traction motor status; additionally, high regenerative energy absorption caused by the large traction motor is available [1]. In the SHEV, electrical energy flows from the generator to the motor through the battery. The effect of the resistance in the battery is significant, as it determines the charging and discharging efficiency of the battery; specifically, a high battery resistance decreases the charging and discharging efficiencies,  Li-ion Battery Internal Resistance (ohm) 0.17 @ SOC 0.5 Ultra-capacitor (ohm) 0.022 To evaluate the energy efficiency of the vehicle, the longitudinal vehicle dynamics model was used, as shown in Figure 2. The tractive resistance consists of rolling, air, gradient, and driving acceleration resistances [29]. Rolling resistance is the force acting opposite to the running direction when the vehicle moves on the road. The friction force is generated by tire hysteresis and the condition of the road surface; it is calculated by:  To evaluate the energy efficiency of the vehicle, the longitudinal vehicle dynamics model was used, as shown in Figure 2. The tractive resistance consists of rolling, air, gradient, and driving acceleration resistances [29].  Li-ion Battery Internal Resistance (ohm) 0.17 @ SOC 0.5 Ultra-capacitor (ohm) 0.022 To evaluate the energy efficiency of the vehicle, the longitudinal vehicle dynamics model was used, as shown in Figure 2. The tractive resistance consists of rolling, air, gradient, and driving acceleration resistances [29]. Rolling resistance is the force acting opposite to the running direction when the vehicle moves on the road. The friction force is generated by tire hysteresis and the condition of the road surface; it is calculated by: cos Rolling resistance is the force acting opposite to the running direction when the vehicle moves on the road. The friction force is generated by tire hysteresis and the condition of the road surface; it is calculated by: where N is the normal force acting on the center of the vehicle, and µ is the rolling resistance coefficient. When a slope θ is applied to the system, the normal force direction changes, and N is replaced by N cos θ as: Air resistance is the frictional force generated by air. As the speed of the vehicle v increases, the air resistance increases and can be calculated by: where ρ is the air density, A f is the vehicle front area, C D is the air resistance coefficient, and v w is the wind speed. The third resistance is the gradient resistance caused by the slope of the road and is calculated by: Adding the accelerating resistance of the vehicle in the running direction to the previous resistances, the longitudinal dynamics equation of the total tractive resistance can be obtained as follows: where the effects of the rotational inertia of the motor, wheel, and transmission are neglected as their influences are negligible considering the mass of the vehicle [30]. As a result, the tractive resistance power of the vehicle can be calculated by multiplying the total tractive resistance with the vehicle speed, as follows: To overcome the tractive resistance of the vehicle, the output torque of the motor can be calculated by: where t f is the first-order lag and τ b is the torque determined by the motor torque limitations [31]. The speed of the motor, ω m , multiplied by the output motor torque, becomes the output power of the motor, and its relationship with (6) is expressed as follows: where η f is the efficiency of the final gear between the motor and wheel. Then, the vehicle speed v can be expressed as: v = rw m g f (9) where r is the radius of the wheel, and g f is the final gear ratio between the motor and wheel.
To drive a wheel, the energy originates from the engine, and the total fuel consumption of the engine during the simulation can be calculated by: (10) where f c is the instant fuel consumption, based on the engine map in Figure 3a, ω e is the engine speed, τ e is the engine torque, and t f is the total simulation time. The engine is mechanically connected to the generator via gears, and its operational status depends on the EMS. The output power of the generator is defined as follows: P g = η g (ω e g ge , τ e /g ge )η ge P e (11) where g ge is the gear ratio between the generator and engine, η g is the efficiency of the generator, η ge is the efficiency of the gear, and P e denotes the output power of the engine.
where ge g is the gear ratio between the generator and engine, g  is the efficiency of the generator, ge  is the efficiency of the gear, and e P denotes the output power of the engine.

ESS Modeling
The SHEB uses the ESS of the battery and UC connected in parallel. The battery usually demonstrates an advantage in energy density over the UC; conversely, the power density of the UC is more advanced than the battery [2,32]. UCs are well known for their high capacitance and semipermanent lifetime [33]. Unlike batteries, which are based on chemical reactions, UCs are charged by simple ion transfers between electrodes, making it feasible to charge and discharge rapidly, with high efficiencies. Moreover, over-charging and over-discharging do not influence the lifetime decrement, unlike batteries. Consequently, UCs are widely used as an ESS, especially as secondary batteries [33,34]. This makes them the core technology of the next generation of ESSs; however, the usage timing and output power of UCs are not clearly presented for improving vehicle efficiency; this study defines them with optimization to take full advantage of the UC. In this study, a resistor and capacitor (RC) model was used for the UC; the steady-state output voltage of the UC can be calculated by [35]: where u Q is the UC charge, C is the capacitance of the UC, u R is the internal resistance of the UC, and u I is the input current of the UC [2,13]. The amount of charge withdrawn from the UC can be calculated by: and the SOC of the UC is expressed by: where i Q , Q , and m Q are the initial charge, existing charge, and maximum charge of the UC, respectively.
The battery must maintain adequate SOC conditions, and the appropriate battery SOC is set to 40-60%, considering the charging/discharging efficiencies and battery energy shortage/overflow [29,36,37]. We adopted the internal resistance equivalent circuit to create the battery model shown in Figure 4, and the battery output power b P was calculated using (15).

ESS Modeling
The SHEB uses the ESS of the battery and UC connected in parallel. The battery usually demonstrates an advantage in energy density over the UC; conversely, the power density of the UC is more advanced than the battery [2,32]. UCs are well known for their high capacitance and semi-permanent lifetime [33]. Unlike batteries, which are based on chemical reactions, UCs are charged by simple ion transfers between electrodes, making it feasible to charge and discharge rapidly, with high efficiencies. Moreover, over-charging and over-discharging do not influence the lifetime decrement, unlike batteries. Consequently, UCs are widely used as an ESS, especially as secondary batteries [33,34]. This makes them the core technology of the next generation of ESSs; however, the usage timing and output power of UCs are not clearly presented for improving vehicle efficiency; this study defines them with optimization to take full advantage of the UC. In this study, a resistor and capacitor (RC) model was used for the UC; the steady-state output voltage of the UC can be calculated by [35]: where Q u is the UC charge, C is the capacitance of the UC, R u is the internal resistance of the UC, and I u is the input current of the UC [2,13]. The amount of charge withdrawn from the UC can be calculated by: and the SOC of the UC is expressed by: where Q i , Q, and Q m are the initial charge, existing charge, and maximum charge of the UC, respectively. The battery must maintain adequate SOC conditions, and the appropriate battery SOC is set to 40-60%, considering the charging/discharging efficiencies and battery energy shortage/overflow [29,36,37]. We adopted the internal resistance equivalent circuit to create the battery model shown in Figure 4, and the battery output power P b was calculated using (15).
where R b and V oc are the internal resistance of the battery and open circuit voltage of the battery, respectively, and the battery current I b can be obtained by:  Figure 4. Equivalent battery model.
The SOC of the battery, h SOC , can also be calculated using (14), where the charge values of the battery, rather than the UC, should be used. When ESS charging occurs, the charging energy is saved first to the UC and then to the battery to take full advantage of the UC, which is more efficient than the battery. In addition, when the SOC of both the battery and UC reaches the minimum allowance level, the UC is not used; using it during this period can cause power supply failure, considering its energy density. Accordingly, the engine/generator and battery are used in this case to satisfy the required power of the motor. Finally, hybrid electric vehicles have a linkage point for combining power sources. DC/DC converters are used for proper control of the ESS terminal voltages and the efficiencies are assumed to be 95% [25,29,38,39].

Energy Management Strategy
In this study, we adopted a practical and robust thermostat control strategy (TCS) as the EMS of the SHEB. The benefits of this strategy are: (1) the efficiency of the vehicle becomes very high as it uses the most efficient operating point of the engine; (2) estimating the energy consumption minimization via ESSs is feasible because we can use the vehicle in the EV mode; and (3) eliminating the effects of other interferences is possible because it uses a single-engine operating point. With other complex EMSs, the developed efficiency improvements can be deduced from other efficiency improvement interferences, such as engine/generator operating conditions.
With this EMS, evaluating the efficiency increment by the UC itself is feasible because the engine/generator is only used for the ESS charger when the SOC reaches the minimum allowance level [24]. However, considering the application principle of the efficiency improvement factors in the next section, the proposed concept is also applicable to other real-time EMSs [40][41][42][43]. The TCS algorithm is given from (17) to (19), where e S , h SOC , and l SOC , represent the engine/generator on/off switching, upper limitation of the battery SOC, and lower limitation of the battery SOC, respectively [28].
The output power of the generator in (11) is modified by TCS as:  The SOC of the battery, SOC b , can also be calculated using (14), where the charge values of the battery, rather than the UC, should be used. When ESS charging occurs, the charging energy is saved first to the UC and then to the battery to take full advantage of the UC, which is more efficient than the battery. In addition, when the SOC of both the battery and UC reaches the minimum allowance level, the UC is not used; using it during this period can cause power supply failure, considering its energy density. Accordingly, the engine/generator and battery are used in this case to satisfy the required power of the motor. Finally, hybrid electric vehicles have a linkage point for combining power sources. DC/DC converters are used for proper control of the ESS terminal voltages and the efficiencies are assumed to be 95% [25,29,38,39].

Energy Management Strategy
In this study, we adopted a practical and robust thermostat control strategy (TCS) as the EMS of the SHEB. The benefits of this strategy are: (1) the efficiency of the vehicle becomes very high as it uses the most efficient operating point of the engine; (2) estimating the energy consumption minimization via ESSs is feasible because we can use the vehicle in the EV mode; and (3) eliminating the effects of other interferences is possible because it uses a single-engine operating point. With other complex EMSs, the developed efficiency improvements can be deduced from other efficiency improvement interferences, such as engine/generator operating conditions.
With this EMS, evaluating the efficiency increment by the UC itself is feasible because the engine/generator is only used for the ESS charger when the SOC reaches the minimum allowance level [24]. However, considering the application principle of the efficiency improvement factors in the next section, the proposed concept is also applicable to other real-time EMSs [40][41][42][43]. The TCS algorithm is given from (17) to (19), where S e , SOC h , and SOC l , represent the engine/generator on/off switching, upper limitation of the battery SOC, and lower limitation of the battery SOC, respectively [28].
The output power of the generator in (11) is modified by TCS as: where τ e, OOP and ω e, OOP are the optimal operating points of the engine torque and speed, respectively. Therefore, the required power of the motor, P r (electric input power of the motor), supplied by the battery, UC, and generator is expressed as follows: where η m and P u denote the efficiency of the motor and output power of the UC, respectively. However, the relation between P b and P u , which is very significant for vehicle efficiency, is not yet clearly defined and is explained in the next section.

Efficiency Improvement Factors
To maximize the efficiency of the SHEB, the UC was adopted, and novel fuel efficiency improvement factors are suggested. The first factor determines whether to use just the battery or not. The second and third factors define the load sharing ratio when simultaneously using the battery and UC.

Required Power of the Motor
The required power of the motor, P r , indicates the input power of the traction motor, and it becomes P b in (19) when a battery alone supplies that power. A large P b indicates a high current in (16), considering the almost constant values of Voc and R i . The battery efficiency is [27] and a high current decreases the efficiency. Accordingly, the threshold of the required power of the motor (TRPM), over which the use of the battery is not allowed, needs to be defined; this is the first efficiency improvement factor. Below the threshold, the battery is very efficient with a small power requirement such that using just the battery as an ESS is reasonable, as shown on the left hill of Figure 5a.

Power Split Ratio between the ESSs
The large power requirement indicates that the traction motor consumes considerable energy, and the required current in the motor is high. Therefore, the efficiency of using just the battery significantly decreases; however, because the UC efficiency is much better than that of the battery, it is capable of significantly reducing energy loss [2,44]. As the participation rate of the UC increases, the charging/discharging loss in the process of passing through the battery is reduced, improving the efficiency of the SHEB. However, an excessive sharing ratio to the UC causes the power supply failure mentioned in Section 2, and the battery should only cover the huge power requirement of the motor, despite the poor efficiency of the battery-only mode. Accordingly, the criteria of the power-split ratio between the battery and UC should be defined as the second and third efficiency improvement factors. The power-split ratio is defined as a linear function to reflect its variation to the motor power requirement, and the slope and y-axis intercept are selected as the second and third efficiency improvement factors, respectively. They are named the slope of the power split ratio (SPSR) and y-axis intercept of the power split ratio (YPSR). Consequently, the power proportion of the UC, r split , becomes: Sustainability 2020, 12, 7354 8 of 19 In Figure 5b, the power split ratio of the UC and battery are r split and 1 − r split , respectively, when the required power of the motor is above the first improvement factor. This is also valid on the right side of Figure 5a. As soon as the required power of the motor passes over the first efficiency improvement factor, the power supply is conducted using the UC and battery simultaneously, where the proportions of r split and 1 − r split change as per the required power of the motor. Therefore, the following equations are used for power supply between the UC and battery:

Particle Swarm Optimization
As mentioned above, TRPM, SPSR, and YPSR require adjusting and hence we adopted the PSO algorithm in Figure 6a; the calculation burden and speed of the PSO algorithm are more advantageous than other global searching algorithms, such as the genetic algorithm (GA). The PSO algorithm is a modern heuristic optimization technique, inspired by the swarm intelligence of birds and ants. A group using swarm intelligence shares the best optimal solution information of both the individual and the group. They interact with each other and help individuals modify their search path such that global and local searches can be conducted [45][46][47]. Recent research has used the PSO algorithm in energy management and the cost optimization of eco-friendly vehicles owing to its powerful performance and ease of implementation [10,48,49].
To verify the effectiveness of the proposed efficiency improvement factors, the particles calculate the objective function for each factor (design variable) and build a personal database. By sharing these data with other particles, the PSO algorithm updates the global best (G best ) and personal best (P best ) solutions in the group. Based on the global and personal best solutions, the particle determines the next movement and each particle moves according to the principle of vector sum, as shown in Figure 6b. The particle equations of the velocity v t p and position x t p of (TRPM, SPSR, YPSR) at time t are expressed as follows: The inertia weight coefficient, ω t p , acceleration coefficients, c 1 and c 2 , and random numbers, r 1 and r 2 , are used in the process of determining the velocity of the particles. The inertia weight coefficient is multiplied by the velocity of the previous step, which determines the magnitude of the inertia velocity. At this time, the velocity and position of the particles are individually determined by the data of the previous step and coefficients. In addition, a linear decreasing weight (LDW) is used for the damping definition such that ω t p decreases with time, as follows [50,51]: where ω start and ω end indicate the initial inertia and final inertia weights, respectively, and iter max is the maximum number of iterations. The gradual decrease in damping ratio makes it possible to determine global and local searching. As a result, particles are searched by iterating the calculations; the PSO algorithm repeats this until the optimum is satisfied. The coefficients and parameters applied to the problem are summarized in Table 2 [52][53][54][55][56].   Figure 6. PSO algorithm: (a) algorithm flowchart; (b) movement of one particle.

Problem Formulation
To apply PSO algorithm, the objective function is defined by the total fuel consumption, f c total , in (10) for minimization and the efficiency improvement factors, TRPM, SPSR, and YPSR, are selected as the design variables. The problem formulation with the constraints is as follows: where TRPM is set to 250 kW or less considering the required power of the motor, which is determined with test cycles. SPSR is defined as positive to increase the participation rate of the UC, regarding the motor power requirement. YPSR is defined between 0 and 1, considering the proportion of the UC. Other variables are used in the powertrain and ESS model in Section 2.

Simulation Results and Discussion
The simulation was performed based on the Manhattan, Braunschweig, and Orange County cycles, as shown in Figure 7. These are realistic urban bus driving cycles; the Manhattan cycle shows a low average speed, the Braunschweig cycle shows frequent bus stops and high average speed, and the Orange County cycle shows relatively irregular and medium speed driving patterns [57].
where TRPM is set to 250 kW or less considering the required power of the motor, which is determined with test cycles. SPSR is defined as positive to increase the participation rate of the UC, regarding the motor power requirement. YPSR is defined between 0 and 1, considering the proportion of the UC. Other variables are used in the powertrain and ESS model in Section 2.

Simulation Results and Discussion
The simulation was performed based on the Manhattan, Braunschweig, and Orange County cycles, as shown in Figure 7. These are realistic urban bus driving cycles; the Manhattan cycle shows a low average speed, the Braunschweig cycle shows frequent bus stops and high average speed, and the Orange County cycle shows relatively irregular and medium speed driving patterns [57].

Optimization Results
The PSO algorithm was conducted and the optimum for each cycle was determined, as shown in Figure 8, where (a) shows the optimization trajectory of the design variables TRPM, SPSR, and YPSR according to the number of iterations, and (b) presents the target (object function) trajectory of the total fuel consumption, according to the number of iterations. The variations in the design variables were clear and we could identify the evident decrease in total fuel consumption as the iteration progressed such that the optimal variables and object functions were determined, as shown in Table 3. In Figure 8a, the values of SPSR and YPSR appear to be small compared to that of TRPM; however, this was caused by the scale of the x-axis (kW) and the y-axis (YPSR maximum of 1) in Figure 5b. according to the number of iterations, and (b) presents the target (object function) trajectory of the total fuel consumption, according to the number of iterations. The variations in the design variables were clear and we could identify the evident decrease in total fuel consumption as the iteration progressed such that the optimal variables and object functions were determined, as shown in Table 3. In Figure  8a, the values of SPSR and YPSR appear to be small compared to that of TRPM; however, this was caused by the scale of the x-axis (kW) and the y-axis (YPSR maximum of 1) in Figure 5b.

Power Supply of ESSs
The three bus driving cycles show different speed profiles, and the required power of the motor is also different for each cycle. These are shown in Figure 9. The values are the same for both the conventional and newly proposed SHEBs.

Power Supply of ESSs
The three bus driving cycles show different speed profiles, and the required power of the motor is also different for each cycle. These are shown in Figure 9. The values are the same for both the conventional and newly proposed SHEBs.   The required power of the motor is supplied by the battery in the conventional SHEB, without the UC; however, the battery-only and load sharing modes based on the efficiency improvement factors in Figure 5 exist in the newly developed SHEB. Figure 10a shows the enlarged view of the Manhattan cycle between 545 and 630 s and the required power of the motor was divided by the battery and UC supply between 545 and 570 s because the TRPM was more than 44.49 kW (the threshold defined by the first improvement factor) in Table 3. The power-split ratio between the UC and battery also followed the SPSR and YPSR in Table 3, i.e., 1.795 × 10 −3 and 0.2575 (the parameters defined by the second and third improvement factors), respectively. Moreover, the required power of the motor was covered by the battery alone between 610 and 621 s, as the required power of the motor was less than the TRPM of 44.49 kW. This was also true for the Braunschweig and Orange County cycles; Figure 10b,c show the load sharing and battery-only modes with an enlarged view of the other cycles. The required power of the motor is supplied by the battery in the conventional SHEB, without the UC; however, the battery-only and load sharing modes based on the efficiency improvement factors in Figure 5 exist in the newly developed SHEB. Figure 10a shows the enlarged view of the Manhattan cycle between 545 and 630 s and the required power of the motor was divided by the battery and UC supply between 545 and 570 s because the TRPM was more than 44.49 kW (the threshold defined by the first improvement factor) in Table 3. The power-split ratio between the UC and battery also followed the SPSR and YPSR in Table 3, i.e., 1.795 × 10 −3 and 0.2575 (the parameters defined by the second and third improvement factors), respectively. Moreover, the required power of the motor was covered by the battery alone between 610 and 621 s, as the required power of the motor was less than the TRPM of 44.49 kW. This was also true for the Braunschweig and Orange County cycles; Figures 10b and c show the load sharing and battery-only modes with an enlarged view of the other cycles.

Battery Efficiency Comparisons
Generally, the UC transfers electrical energy at a higher efficiency than the battery, so that the energy loss caused by the charging and discharging efficiency could be significantly improved. However, the UC energy density limitation prohibits the general use of the UC such that optimization should be conducted. The PSO algorithm was adopted and the battery and UC were simultaneously used to improve battery efficiency. Figure 11a shows the battery efficiency of the Manhattan cycle of Figure 10a when comparing it to the conventional SHEB, without the UC, using the enlarged view.
Observing the difference between 545 and 570 s, the increment of the battery efficiency is evident, with a maximum difference of 9.7% (95.2% vs. 85.5% at 555 s). By contrast, the efficiency of the UC is also considered in Figure 11b and was more than 98% during this period. These phenomena were equally applicable to the Braunschweig and Orange County cycles, and the corresponding results are shown in Figure 11c-f. Generally, the UC transfers electrical energy at a higher efficiency than the battery, so that the energy loss caused by the charging and discharging efficiency could be significantly improved. However, the UC energy density limitation prohibits the general use of the UC such that optimization should be conducted. The PSO algorithm was adopted and the battery and UC were simultaneously used to improve battery efficiency. Figure 11a shows the battery efficiency of the Manhattan cycle of Figure 10a when comparing it to the conventional SHEB, without the UC, using the enlarged view. Observing the difference between 545 and 570 s, the increment of the battery efficiency is evident, with a maximum difference of 9.7% (95.2% vs. 85.5% at 555 s). By contrast, the efficiency of the UC is also considered in Figure 11b and was more than 98% during this period. These phenomena were equally applicable to the Braunschweig and Orange County cycles, and the corresponding results are shown in Figure 11c Figure 12 shows the battery SOC of the conventional SHEB, using the battery alone, and that of the newly proposed SHEB, using the battery and UC. The engine is turned on later in the case of the proposed SHEB than in a conventional SHEB, as the battery shares some load with the UC when a lot of power is required by the motor. In this study, to equally match the total energy consumption of the various cases, the initial and final SOCs were set as equal; the battery SOC started and ended at 60%, and the UC SOC started and ended at 100%. To adjust the final SOC at the end of the driving cycles, we operated the engine to charge the battery and UC to their initial states after driving, as shown in Figure 12.  Figure 12 shows the battery SOC of the conventional SHEB, using the battery alone, and that of the newly proposed SHEB, using the battery and UC. The engine is turned on later in the case of the proposed SHEB than in a conventional SHEB, as the battery shares some load with the UC when a lot of power is required by the motor. In this study, to equally match the total energy consumption of the various cases, the initial and final SOCs were set as equal; the battery SOC started and ended at 60%, and the UC SOC started and ended at 100%. To adjust the final SOC at the end of the driving cycles, we operated the engine to charge the battery and UC to their initial states after driving, as shown in Figure 12.

Comparison of Fuel Efficiency
For the Manhattan, Braunschweig, and Orange County cycles, the total fuel consumptions of the conventional SHEB shown in Figure 13 were 1378 g, 3761 g, and 3679 g, respectively, while those of the newly proposed SHEB after PSO algorithm optimization were 1284 g, 3548 g, and 3499 g, respectively.
Consequently, the total fuel consumption decreased by 6.82% in the Manhattan cycle, 5.66% in the Braunschweig cycle, and 4.89% in the Orange County cycle, which verifies the effectiveness of the proposed efficiency improvement factors in the SHEB with the UC. The efficiency changes due to the PSO algorithm are also presented in Figure 13. In each figure, the second bars indicate the total fuel consumption with the random starting TRPM, SPSR, and YPSR in Figure 8a,c,e and the third bars denote the total fuel consumption after optimization. The improvements were 3.16%, 3.37%, and 2.59% for the Manhattan, Braunschweig, and Orange County driving cycles, respectively. As a result, a SHEB, with high efficiency using a battery and UC, could be realized and verified. The number of transit buses in a city is huge, and buses generally run with similar driving patterns. In addition, we can know the traffic information in advance using the bus traffic information system. Therefore, defining the corresponding TRPM, SPSR, and YPSR according to the driving information with the SHEB would be a great help for fuel economy enhancement and CO 2 reduction in the city.

Comparison of Fuel Efficiency
For the Manhattan, Braunschweig, and Orange County cycles, the total fuel consumptions of the conventional SHEB shown in Figure 13 were 1378 g, 3761 g, and 3679 g, respectively, while those of the newly proposed SHEB after PSO algorithm optimization were 1284 g, 3548 g, and 3499 g, respectively. Consequently, the total fuel consumption decreased by 6.82% in the Manhattan cycle, 5.66% in the Braunschweig cycle, and 4.89% in the Orange County cycle, which verifies the effectiveness of the proposed efficiency improvement factors in the SHEB with the UC. The efficiency changes due to the PSO algorithm are also presented in Figure 13. In each figure, the second bars indicate the total fuel consumption with the random starting TRPM, SPSR, and YPSR in Figure 8 and the third bars denote the total fuel consumption after optimization. The improvements were 3.16%, 3.37%, and 2.59% for the Manhattan, Braunschweig, and Orange County driving cycles, respectively. As a result, a SHEB, with high efficiency using a battery and UC, could be realized and verified. The number of transit buses in a city is huge, and buses generally run with similar driving patterns. In addition, we can know the traffic information in advance using the bus traffic information system. Therefore, defining the corresponding TRPM, SPSR, and YPSR according to the driving information with the SHEB would be a great help for fuel economy enhancement and CO2 reduction in the city.

Conclusions
Energy management optimization of the SHEB using a UC was studied. Novel efficiency improvement factors for defining the power supply between the battery and UC-TRPM, SPSR, and YPSR-were proposed, their efficiencies were considered, and the PSO algorithm was adopted to verify the effectiveness of the proposed efficiency improvement factors. From the results, it was confirmed that realizing a higher and sustainable SHEB without changing much of the structure is feasible, which would not have been the case for the existing SHEB. It was confirmed that the total fuel consumption of the proposed SHEB was reduced significantly, by 6.82% in the Manhattan cycle, 5.66% in the Braunschweig cycle, and 4.89% in the Orange County cycle, when compared to the results of the existing SHEB. All results were obtained using the MATLAB/Simulink and Amesim co-

Conclusions
Energy management optimization of the SHEB using a UC was studied. Novel efficiency improvement factors for defining the power supply between the battery and UC-TRPM, SPSR, and YPSR-were proposed, their efficiencies were considered, and the PSO algorithm was adopted to verify the effectiveness of the proposed efficiency improvement factors. From the results, it was confirmed that realizing a higher and sustainable SHEB without changing much of the structure is feasible, which would not have been the case for the existing SHEB. It was confirmed that the total fuel consumption of the proposed SHEB was reduced significantly, by 6.82% in the Manhattan cycle, 5.66% in the Braunschweig cycle, and 4.89% in the Orange County cycle, when compared to the results of the existing SHEB. All results were obtained using the MATLAB/Simulink and Amesim co-simulation, and the proposed improvement factors with the UC are also applicable to other real-time EMSs, for better performance. In the future, a study comparing optimization algorithms to find more efficiency improvement factors will be carried out.