An Optimal Sizing Design Approach of Hybrid Energy Sources for Various Electric Vehicles

In this paper, we present a discussion about green energy sources that have been widely utilized in electric vehicles (EVs). To achieve different requirements of various EVs, the correct sizing of energy sources is crucial so that the cost and output performance will be optimized. In this research, three energy sources, supercapacitors (SCs), lithium titanate oxide (LTO) batteries, and Nickel Manganese Cobalt (NCM) (or Li3) batteries, were considered for hybridization. An effective global search algorithm (GSA) was designed for optimal sizing of hybrid electric energy systems (HEESs). The GSA procedures were: (1) vehicle specification and performance requirements of energy sources, (2) determination of cost function and constraints, (3) GSA optimization with for-loops, (4) optimal results. Five examples of EVs, the electric sedan, long-distance electric bus, short-distance electric bus, electric forklift, and electric sports car, were analyzed for optimal hybrid energy combination under different criteria and specifications. The GSA effectively optimized the designs of energy sizing. The performance indices and vehicle requirements studied were the specific price, specific energy at a constant volume, specific energy at a constant mass, and specific power at a constant mass for three energy sources, SCs, LTO batteries, and Li batteries. The vehicle requirements including the maximum output power, vehicle acceleration, climbability, and maximum speed have been formulated as the design constraints. A numerical analysis of five types of EVs was analyzed for optimal sizing of the HEES and the optimal position of the DC/DC converter with the lowest cost function. The integrated system and control designs of the HESS using the GSA, more applications for green energy sources, and different types of EVs will be studied in the future.


Introduction
Due to the stringent transportation laws legislated, especially for fuel economy, electric vehicles (EVs) have attracted much more attention recently due to their characteristics of zero pollutants, low noise, high efficiency, and rapid acceleration compared to traditional engine vehicles [1]. The most popular EV type is the battery EV (BEV) with one type of battery, such as LiFePO4 batteries, lithium titanate oxide (LTO) batteries, nickel manganese cobalt (NMC) batteries, or nickel cobalt aluminium (NCA) batteries [2]. Each battery at the rear parts of dual energy sources, respectively. By properly controlling the control unit, the energy management is determined by two converters while the output voltages are stabilized at the downstream of the converters. The drawbacks of (1) are the extra cost and occupied volume of the converters. For passive type, the in/out energy flows of two energy sources are determined only by output voltages of the two sources without the energy management. Though the cost is lowest and the structure is simplest, the state-of-health (SOH) of both sources might decrease rapidly due to improperly drawn currents or lack of SOC control. For active-passive type 1, one DC/DC converter is at the output of the Li3 battery. As the output current (power) of the DC/DC converter is controlled while the output voltage is stabilized, the output power of the LTO is passively given. For active-passive type 2, one DC/DC is at the output of the LTO battery. As the output current (power) of the DC/DC is controlled while the output voltage is stabilized, the output power of the Li3 is passively given. Both active-passive types exhibit the same performance as that of the active type. However, the overall system efficiency is higher because the extra input/output power loss of the second converter is absent. Moreover, the system cost and occupied volume and weight are lower due to the single DC/DC converter. Though the cost is higher than that of the passive type, the SOH of dual energy sources can be extended to lower the cost. Of these four types, two active-passive types were chosen for system sizing optimization. The position of the single DC/DC converter was selected using the proposed GSA. Note that the system efficiency of hybrid energy sources is significantly higher than that of the single source, especially at the EoL period [23]. Hence, the single-source system was not considered in this research.
According to our literature review, the four configurations of HEESs are separated as shown in Figure 1. They are: (1) active, (2) passive, (3) active-passive type 1 and (4) active-passive type 2. For the active type (Figure 1a), two DC/DC converters are equipped at the rear parts of dual energy sources, respectively. By properly controlling the contro unit, the energy management is determined by two converters while the output voltages are stabilized at the downstream of the converters. The drawbacks of (1) are the extra cost and occupied volume of the converters. For passive type, the in/out energy flows of two energy sources are determined only by output voltages of the two sources without the energy management. Though the cost is lowest and the structure is simplest, the state-ofhealth (SOH) of both sources might decrease rapidly due to improperly drawn currents or lack of SOC control. For active-passive type 1, one DC/DC converter is at the output of the Li3 battery. As the output current (power) of the DC/DC converter is controlled while the output voltage is stabilized, the output power of the LTO is passively given. For active-passive type 2, one DC/DC is at the output of the LTO battery. As the output curren (power) of the DC/DC is controlled while the output voltage is stabilized, the output power of the Li3 is passively given. Both active-passive types exhibit the same perfor mance as that of the active type. However, the overall system efficiency is higher because the extra input/output power loss of the second converter is absent. Moreover, the system cost and occupied volume and weight are lower due to the single DC/DC converter Though the cost is higher than that of the passive type, the SOH of dual energy sources can be extended to lower the cost. Of these four types, two active-passive types were chosen for system sizing optimization. The position of the single DC/DC converter was selected using the proposed GSA. Note that the system efficiency of hybrid energy sources is significantly higher than that of the single source, especially at the EoL period [23] Hence, the single-source system was not considered in this research.

The Procedures of the GSA
In this section, we developed the GSA for HEESs of EVs. The procedure was separated into four steps: (1) vehicle specification and performance requirements of energy sources, (2) determination of cost function and constraints, (3) GSA optimization with forloops, (4) optimal results. The following is the procedure of the proposed GSA:

The Procedures of the GSA
In this section, we developed the GSA for HEESs of EVs. The procedure was separated into four steps: (1) vehicle specification and performance requirements of energy sources, (2) determination of cost function and constraints, (3) GSA optimization with for-loops, (4) optimal results. The following is the procedure of the proposed GSA: (1) Vehicle specification and performance requirements of energy sources; in this step, two phases are designed separately. The first phase is the performance indices of energy sources, while the second phase is the specification of requirements of the target vehicle. The performance indices of energy sources are: (a) specific energy (kWh/kg or kWh/L) (the energy stored in a constant mass or in a constant volume), as the higher the specific energy, the more electric energy is able to be released to drive the traction motor in an EV so that the traveling mileage is larger; (b) specific power (kW/kg or kW/L) (the maximized power that can be released in a constant mass or in a constant volume), as the higher the specific power, the more electric power is able to propel the traction motor so that the values of vehicle acceleration and maximum vehicle speed are higher; and (c) specific price (USD/kWh) (the cost for energy capacity of 1 kW-hr), as the lower the specific cost, at a constant stored energy, the lower the overall cost [24].
The second phase is the vehicle requirements and specification. The key indices are (a) traveling mileage; (b) maximal power; (c) climbability; and (d) acceleration. These indices are criteria for target cost function, which will be described later. For traveling mileage, it is correlated to the energy stored in the HEES. For a chosen specific driving cycle, the consumed energy for one cycle as well as the traveling distance can be calculated by the following two equations: where E, P, F, D, V, and t represent energy demand, power demand, traction force of the wheel, distance, vehicle velocity, and time step, respectively; the subscript i represents the driving time(s). Hence, the total distance decides the minimum requirement of energy storage. For maximal power, it is related to the tested initial speed, target speed, and the varied vehicle mass determined by the HEES sizing. The two following equations are listed as [3,4]: where P max , ∆E, ∆t, and M represent maximum power demanded, kinetic energy, demanded time(s), and mass, respectively; the subscripts veh, init., and goal. represent vehicle, initial, and goal, respectively. For climbability, it is correlated with the maximum power requirement as well. Modified from Equations (5)-(8), the calculated power is formulated as: where F, g, θ, µ, C d , ρ, A f , r, and τ represent force, gravity acceleration, road slope, rolling coefficient, aero drag coefficient, air density, frontal area, radius, and traction torque, respectively; the subscripts climb, roll, wind, and w represent climbing resistance, rolling resistance, wind resistance, and wheel, respectively. For acceleration, it can be determined by the profiles in Figure 2.
Appl. Sci. 2022, 12, x FOR PEER REVIEW (2) Determination of the cost function and constraints; to optimize the HEES iable of the cost function is set to be the total cost of the HEES. It can be listed as: where C is cost, and the subscripts , 3, and represent the HESS, Li3 and LTO battery set, respectively. Individual costs of dual energy sources are lis where , , , , , , , , and represent stored energy, size ra cific price (USD/kWh), power density (kw/kg), specific weight (kWh/kg), specific DC/DC converter (USD/kW), index for searching the DC/DC converter, index fo ing energy of HEES, and index for searching the size ratio, respectively.
Thus the cost function can be defined as follows: where , , and γ represent volume, weight, and penalty value, respectively.
(3) GSA optimization with for-loops; after deriving the cost function, perf (2) Determination of the cost function and constraints; to optimize the HEES, the variable of the cost function is set to be the total cost of the HEES. It can be listed as: where C is cost, and the subscripts HEES, li3, and LTO represent the HESS, Li3 battery, and LTO battery set, respectively. Individual costs of dual energy sources are listed as: (11) where E, SR, SP, PD, SW, SPDC, i, j, and k represent stored energy, size ratio, specific price (USD/kWh), power density (kw/kg), specific weight (kWh/kg), specific price of DC/DC converter (USD/kW), index for searching the DC/DC converter, index for searching energy of HEES, and index for searching the size ratio, respectively. Thus the cost function J can be defined as follows: subject to γ(i, j, k) = 0 if above equlation satisfied 10 6 if above equlation unsatisfied (17) where V, W, and γ represent volume, weight, and penalty value, respectively. (3) GSA optimization with for-loops; after deriving the cost function, performance requirements, and the performance indices of dual energy sources, the GSA was proposed, as shown in Figure 3. First, three for-loops were constructed. The most outer loop was the DC/DC position with two possible values: i = [0, 1]. The second loop was the required total energy of the HESS with bounded values from E HEES,min to E HEES,max . The increment of E HEES,inc separated for the index of E HEES was from j = 1 (E HEES,min ) to j = n (E HEES,max ).
The total number n was calculated by: n = [EHEES,max−EHEES,min] E HEES,inc + 1. The most inner loop was the size ratio with bounded values from 0 to 1. The increment of SR inc separated the index of SR from k = 1 (SR = 0) to k = m (SR = 1). The total number m was calculated from: m = SR inc + 1. Therefore, the possible combinations of DC/DC converter positions, total energies of the HEES, and size ratios were 2 × n × m cases. At each case with a determined DC/DC converter position, total HEES energy, and size ratio, the energy of the Li3 battery was E li3 = E HEES (j) × SR(k), while the energy of the LTO battery was Next, physical constraints were calculated based on Equations (13)-(16) including D max , P max , V max , and W max . The value of γ can be determined if any of the constraints are satisfied or unsatisfied. The three-dimensional database with indices i, j, and k was constructed.  (4) Optimal results; to determine the optimal solution at a specific DC/DC converter position and stored hybrid energy, the cost functions with various size ratios were compared. The cost function with a minimum value under a chosen size ratio was the best solution, while the size ratio was the optimal size ratio at a specific X (position of DC/DC converter). That is: Next, for the outer layer, using different DC/DC converter positions, the final optimal solution can be determined and chosen from 1 * . Figure 4 summarizes the flow chart of the GSA for optimal size ratio: (1) vehicle specification and performance requirements of energy sources, (2) determination of the cost function and constraints, (3) GSA optimization with for-loops, and (4) optimal results. Using the four steps, the optimal HEES will be found. (4) Optimal results; to determine the optimal solution at a specific DC/DC converter position and stored hybrid energy, the cost functions with various size ratios were compared. The cost function with a minimum value under a chosen size ratio was the best solution, while the size ratio was the optimal size ratio at a specific X (position of DC/DC converter). That is: Next, for the outer layer, using different DC/DC converter positions, the final optimal solution can be determined and chosen from J 1 * . Figure 4 summarizes the flow chart of the GSA for optimal size ratio: (1) vehicle specification and performance requirements of energy sources, (2) determination of the cost function and constraints, (3) GSA optimization with for-loops, and (4) optimal results. Using the four steps, the optimal HEES will be found.

Electric Sedan
With the rapid popularization of commercialized electric sedans, all international automakers have produced several types of EVs. However, the trade-off between high power (good acceleration and high maximum speed) and high energy (traveling mileage) makes a vehicle with a single energy source difficult to produce. To optimize the HEES, the cost function is listed as follows: where the optimized cost is calculated by the summation of the cost of a high power source (LTO), the cost of a high energy source (Li3), and the cost of a DC/DC converter. The constraints are as follows: where the constants 0.95, 1.22, and 50 represent battery efficiency, safety factor, and power demand at the constant velocity of 50 km/h, respectively. Because the battery efficiency was influenced by SOC, which varies for energy management but not for the sizing optimization, the efficiency was set as a constant average value of 95%. For the safety factor, 1.22, for the weight, cost, and volume, the references [25,26] were considered. The

Electric Sedan
With the rapid popularization of commercialized electric sedans, all international automakers have produced several types of EVs. However, the trade-off between high power (good acceleration and high maximum speed) and high energy (traveling mileage) makes a vehicle with a single energy source difficult to produce. To optimize the HEES, the cost function is listed as follows: where the optimized cost is calculated by the summation of the cost of a high power source (LTO), the cost of a high energy source (Li3), and the cost of a DC/DC converter. The constraints are as follows: where the constants 0.95, 1.22, and P 50 represent battery efficiency, safety factor, and power demand at the constant velocity of 50 km/h, respectively. Because the battery efficiency was influenced by SOC, which varies for energy management but not for the sizing optimization, the efficiency was set as a constant average value of 95%. For the safety factor, 1.22, for the weight, cost, and volume, the references [25,26] were considered. The subscripts HE and HP represent a high-energy and high-power battery, respectively. To satisfy the first constraint, P max , the maximal power of high energy or the high power were derived from the electric capacity divided by the specific weight and the specific price. The maximum cost of the HEES was a factor (1.22) multiplied by the summation of the total cost of the LTO batteries, Li3 batteries, and the DC/DC converter. The maximum weight of the HEES was an allowable factor (i.e., 1.22) multiplied by the summation of the total weight of the LTO batteries and the Li3 batteries. The minimum electric capacity was calculated by giving the required traveling distance at a constant speed of 50 km/h. With the boundary conditions of these four limitations, the optimal cost function concerning the vehicle performance requirement can be globally searched for.

Electric Forklift
Similar to Section 3.1, the cost function was the total cost of the HEES with a penalty factor. The major concern was the required maximum power while the forklift lifts up the goods. Hence, in Equation (22), the summation of the dual power sources multiplied by an efficiency factor (i.e., 0.95) needs to be higher than a preset maximum power value. 0.95 × (P HP + P HE ) > P max

Long-Distance Electric Bus
For the long-distance electric bus, the traveling distance was crucial. As Equation (20) shows, the cost function was the total cost of the HEES to be minimized. The major concern was the required maximum power. The summation of the stored energy needed to be higher than a preset value for a long traveling distance. Due to the limited space of the chassis to mount the HEES system, the summation of the occupied volume of dual sources multiplied by a safety factor (1.22) needed to be smaller than a preset space (V max ).

Short-Distance Electric Bus
For the short-distance electric bus, the traveling distance was a minor concern due to a fixed route for operation. Instead, the weight of the HEES, which influences the power consumption during frequent acceleration and deceleration, was crucial. Hence, the major concern was the HEES weight, which needed to be less than a preset value (W max ). Next, the summation of the HEES multiplied by an efficiency factor (i.e., 0.95) needed to be higher than a preset maximum power as well. Due to the limited chassis space to mount the HEES system, the summation of the occupied volume of dual sources multiplied by a safety factor (1.22) needed to be smaller than a preset space.

Electric Sports Car
For the electric sports car, the HEES weight was a key issue because of the stringent requirement of acceleration and deceleration. Hence, the major concern was the HEES weight, which needed to be less than a preset value (W max ). Next, the maximum vehicle speed was correlated to the maximum output power, so the summation of the HEES multiplied by an efficiency factor (i.e., 0.95) needed to be higher than a preset value. The limited space for the HEES was another concern, where the summation of the occupied volume of dual sources multiplied by a safety factor (1.22) must be smaller than a preset space.
From the constraints proposed for the five types of vehicles, the GSA is developed in Section 2 and will search for the optimal size ratio and the proper position of the DC/DC converter.

Simulation Results and Discussion
Firstly, the parameters and specifications of energy sources including the specific power, specific price, specific energy at constant volume, and specific energy at constant mass were searched and averaged for further use. Then, the performance with respect to size ratio for the cases of SC and Li3 battery and LTO and Li3 battery were simulated. Next, the sizing of five types of EVs was analyzed using the GSA.

Parameters and Specifications of Energy Sources
The performance indices of three energy storage devices (SC, LTO, and Li3 battery) needed to be searched in advance. Note that though the costs of materials of batteries were similar, the labor cost and management cost widely vary, which causes large differences in the final costs and influences the key indices. For the viewpoint of this purely academic research, we only retrieved the data from research papers rather than focusing on specific market areas or specific time periods. However, to avoid the large difference of a specific index, three data-cancellation rules were set to judge and delete the extreme values: (1) the maximum value is more than or equal to two times the minimum values; (2) the difference of average values before and after deleting the maximum or minimum value is more than 10%; (3) from the normal distribution, the Z score (standard score) is higher than 1.23.
In Table 1, for supercapacitors, the averaged value of specific power was 7.3 (kW/kg) [24][25][26], the value of specific price [24] was 10,000 (USD/kWh), the averaged value of specific energy at constant volume was 0.0085 (kWh/L) [27], and the averaged value of specific energy at constant mass was 0.00484 (kWh/kg) [24,27,28]. The specific price and specific power was relatively high. Note that not all data-cancellation rules were satisfied (e.g., for specific power, the maximum value was two times the minimum value and the average values before and after deleting 10 were 7.3 and 5.95 (an 18.50% difference); however, the Z-score was 1.069, which was below 1.10, so all data were used for calculation). 6.9 [26] Specific price (USD/kWh) 10,000 [24] 10,000 Power density (kWh/L) 0.0085 [27] 0.0085 Specific energy (kWh/kg) 0.005 [24] 0.00484 0.0065 [27] 0.00302 [28] In Table 2, for LTO batteries, the average value of specific power was 2.763 (kW/kg) [2,17,29,30]; the value of specific price was 578.66 (USD/kWh) [2,29,31]; the average value of specific energy at constant volume was 0.138 (kWh/L) [30][31][32]; and the average value of specific energy at constant mass was 0.0635 (kWh/kg) [2,17,[30][31][32][33]. Note that concerning the data-cancellation rules, for specific energy, the maximum value was nine times the minimum value, the average values before and after deleting 0.011 were 0.056 and 0.0635 (a 13.39% difference), and the Z-score was 1.62, which is higher than 1.10, so the data point 0.11 was deleted.

Optimization of Energy Sources
To fundamentally analyze four performance indices (mass, electric capacity, volume, and output power) with respect to size ratio and total price or electric capacity of the HEES, the GSA without the vehicle performance requirements (Equations (13)-(16)) was evaluated. There are two possible combinations: SC/Li3 and LTO/Li3 HEES. Figure 5 illustrates the performance of SC/Li3 with respect to total price and size ratio. In Figure 5a, at a fixed price, if the size ratio increases (more SC), the mass increases due to the high specific cost of the SC. In Figure 5b, at the fixed price, the electric capacity decreases as the SR increases due to the low energy density of the SC. In Figure 5c, at a fixed price, the volume increases due to the low power density at the constant volume. In Figure 5d, at a constant price, the output power increases as the SR increases due to the high power density of the SC. Table 3. The performance indices of the NCM (Li3) battery.

Optimization of Energy Sources
To fundamentally analyze four performance indices (mass, electric capacity, volume, and output power) with respect to size ratio and total price or electric capacity of the HEES, the GSA without the vehicle performance requirements (Equations (13)-(16)) was evaluated. There are two possible combinations: SC/Li3 and LTO/Li3 HEES. Figure 5 illustrates the performance of SC/Li3 with respect to total price and size ratio. In Figure 5a, at a fixed price, if the size ratio increases (more SC), the mass increases due to the high specific cost of the SC. In Figure 5b, at the fixed price, the electric capacity decreases as the SR increases due to the low energy density of the SC. In Figure 5c, at a fixed price, the volume increases due to the low power density at the constant volume. In Figure 5d, at a constant price, the output power increases as the SR increases due to the high power density of the SC. Another example is the LTO/Li3 HEES case. In Figure 6a, at a fixed price, if the size ratio increases (more LTO battery), the mass increases due to the lower specific energy at constant mass of the LTO battery referred to in Tables 2 and 3. In Figure 6b, at the fixed price, the electric capacity decreases as the SR increases due to the lower energy density of the LTO compared to Li3 batteries. In Figure 6c, at a fixed price, the occupied volume increases due to the lower specific energy at the constant volume. In Figure 6d, at a constant price, the output power increases as the SR increases because of the higher power density of the LTO batteries. Another example is the LTO/Li3 HEES case. In Figure 6a, at a fixed price, if the size ratio increases (more LTO battery), the mass increases due to the lower specific energy at constant mass of the LTO battery referred to in Tables 2 and 3. In Figure 6b, at the fixed price, the electric capacity decreases as the SR increases due to the lower energy density of the LTO compared to Li3 batteries. In Figure 6c, at a fixed price, the occupied volume increases due to the lower specific energy at the constant volume. In Figure 6d, at a constant price, the output power increases as the SR increases because of the higher power density of the LTO batteries. Appl. Sci. 2022, 12, x FOR PEER REVIEW 12 of 18 From another view point, to investigate the relationship of multidimensional performance indices, we chose the electric capacity as the y-axis variable. In Figure 7a, at a fixed capacity, if the size ratio increases (more SC modules), the total cost increases due to the extremely high specific cost of the SC referred to in Tables 1 and 3. In Figure 7b, at the fixed capacity, the occupied volume of the HEES increases as the SR increases due to the lower energy density of the SC compared to Li3 batteries. To remain at the same capacity, the volume increases. In Figure 7c, at a fixed capacity, the HEES mass increases due to the lower specific energy at the constant volume of the SC. In Figure 7d, at a constant capacity, the output power increases as the SR increases due to the higher power density of the SC module.  From another view point, to investigate the relationship of multidimensional performance indices, we chose the electric capacity as the y-axis variable. In Figure 7a, at a fixed capacity, if the size ratio increases (more SC modules), the total cost increases due to the extremely high specific cost of the SC referred to in Tables 1 and 3. In Figure 7b, at the fixed capacity, the occupied volume of the HEES increases as the SR increases due to the lower energy density of the SC compared to Li3 batteries. To remain at the same capacity, the volume increases. In Figure 7c, at a fixed capacity, the HEES mass increases due to the lower specific energy at the constant volume of the SC. In Figure 7d, at a constant capacity, the output power increases as the SR increases due to the higher power density of the SC module.
For the case of the LTO/Li3 HEES, in Figure 8a, at a fixed electric capacity, if the size ratio increases (more LTO batteries), the total cost slightly increases due to the higher specific cost of the LTO compared with Li3 batteries, referred to in Tables 2 and 3. In Figure 8b, at the fixed capacity, the occupied volume of the HEES increases as the SR increases due to the lower energy density of the LTO compared to Li3 batteries. To remain at the same capacity, the volume increases. In Figure 8c, at a fixed capacity, the HEES mass increases due to the lower specific energy at the constant volume of the LTO battery. In Figure 8d, at a constant capacity, the output power increases as the SR increases due to the higher power density of the LTO module. From the example above, the relationship of performance indices with respect to the combinations of LTO/Li3 HEES and SC/Li3 HEES was constructed. It can be applied for the vehicle designs in Section 4.2.2.

Case Study of Five Types of EVs
In this Section, the GSA developed in Section 2.2 and the vehicle performance requirements and constraints in Section 3 were simulated. Table 4 demonstrates the optimal sizing of the electric sedan. The vehicle performance requirements (constraints) were listed in Section 3.1. Table 4 shows that the case of LTO/Li3 combination was the best solution due to the lowest cost (33,111 USD) compared to that of the SC/Li3 case (33,842 USD). For the electric capacity, the optimal sizing was that the capacity of the LTO battery (as the high-power source) was 7.459 kWh, where the capacity of the Li3 battery (as the high-energy source) was 67.541 kWh. The size ratio was 9.9:90.1 where the total required capacity was 75 kWh. With the optimal sizing, the maximum output power of the LTO battery was 324.5 kW, while that of the Li3 battery was 139.8 kW.
The combined power exceeded the required power of 464 kW. The HEES mass (453 kg) was slightly higher than that of the SC/Li3 case (416 kg). Similarly, the occupied volume (222L) was slightly more than that of the SC/Li3 (212 L) case. The optimal position of the DC/DC converter was at the high-energy (Li3 battery) side. To prove that the GSA searched the optimal results in the electric sedan case, we manually and randomly selected cases that satisfied the physical constrains. For LTO/Li3 optimization, as the electric capacity was 7.55 and 67.45 kWh, respectively, the total price was calculated to be 33,151 USD. It increased by 0.12% compared to the GSA optimal result. The case with the highest price that satisfied the physical constraints was 33,181 USD, which was increased by 0.21%. In the case of SC/Li3 hybridization, the electric capacity was 0.33 and 74.64 kWh, respectively, where the price increased to 35,968 USD, which was increased by 6.3%. The case with the highest price that satisfied the physical constraints was 38,206 USD, which was increased by 12.89%. From another view point, to investigate the relationship of multidimensional performance indices, we chose the electric capacity as the y-axis variable. In Figure 7a, at a fixed capacity, if the size ratio increases (more SC modules), the total cost increases due to the extremely high specific cost of the SC referred to in Tables 1 and 3. In Figure 7b, at the fixed capacity, the occupied volume of the HEES increases as the SR increases due to the lower energy density of the SC compared to Li3 batteries. To remain at the same capacity, the volume increases. In Figure 7c, at a fixed capacity, the HEES mass increases due to the lower specific energy at the constant volume of the SC. In Figure 7d, at a constant capacity, the output power increases as the SR increases due to the higher power density of the SC module.  For the case of the LTO/Li3 HEES, in Figure 8a, at a fixed electric capacity, if the size ratio increases (more LTO batteries), the total cost slightly increases due to the higher specific cost of the LTO compared with Li3 batteries, referred to in Tables 2 and 3. In Figure  8b, at the fixed capacity, the occupied volume of the HEES increases as the SR increases due to the lower energy density of the LTO compared to Li3 batteries. To remain at the same capacity, the volume increases. In Figure 8c, at a fixed capacity, the HEES mass increases due to the lower specific energy at the constant volume of the LTO battery. In Figure 8d, at a constant capacity, the output power increases as the SR increases due to the higher power density of the LTO module. From the example above, the relationship of performance indices with respect to the combinations of LTO/Li3 HEES and SC/Li3 HEES was constructed. It can be applied for the vehicle designs in Section 4.2.2. In this Section, the GSA developed in Section 2.2 and the vehicle performance requirements and constraints in Section 3 were simulated. Table 4 demonstrates the optimal sizing of the electric sedan. The vehicle performance requirements (constraints) were listed in Section 3.1. Table 4 shows that the case of  Table 5 demonstrates the optimal sizing of the electric forklift [36]. The case of SC/Li3 combination was the best solution, as the case of LTO/Li3 was unable to satisfy the required maximum output power. The HEES cost was calculated to be 24,973 USD. For the electric capacity, the optimal sizing was that the capacity of the SC module (as the high-power source) was 1.357 kWh, where the capacity of the Li3 battery (as the high-energy source) was 28.643 kWh. The size ratio was 4.52:95.48, where the total required capacity was 30 kWh. With the optimal sizing, the maximum output power of the SC was 2046 kW, while that of the Li3 battery was 59.3 kW. The combined power exceeded the required power of 2000 kW. The HEES mass was 422 kg, while the occupied volume was 231 L. The position of the DC/DC converter was at the SC side. To prove that the GSA was useful in the electric forklift case, we manually selected cases that satisfied the physical constrains. For SC/Li3 hybridization, as the electric capacity was 1.5 and 28.5 kWh, respectively, the total price was 26,347 USD. It increased 5.5% compared to the GSA optimal result. In the case that the electric capacity was 2 (SC) and 28 (Li3) kWh, respectively, the price was 31,154 USD, which was increased by 24.75%. Table 6 demonstrates the optimal sizing of the long-distance electric bus. No matter LTO/Li3 or SC/Li3 cases, only using the Li3 battery was the best solution because of the long travelling distance. In this case, the Li3 battery cost was calculated to be 39,750 USD. For the electric capacity, the optimal sizing was that the capacity of the battery was 100 kWh. The HP:HE size ratio was 0:100, where the total required capacity was 100 kWh. The maximum output power of the battery was 442 kW. This power exceeded the required power of 300 kW. The HEES mass was 9212 kg, while the occupied volume was 200 L. There was no DC/DC converter for the HEES structure. Table 7 demonstrates the optimal sizing of the short-distance electric bus. The case of LTO/Li3 combination was the best solution, as the case of SC/Li3 was unable to satisfy the required electric capacity described in Section 3. The HEES cost was calculated to be 34,914 USD. For the electric capacity, the optimal sizing was that the capacity of the LTO module (as the high-power source) was 32 kWh, where the capacity of the Li3 battery (as the high-energy source) was 40 kWh. The size ratio was 44.4:55.6, where the total required capacity was 72 kWh. With the optimal sizing, the maximum output power of the LTO battery was 1392 kW, while that of the Li3 battery was 83 kW. The combined power exceeded the required power of 528 kW. The HEES mass was 13,702 kg, while the occupied volume was 331 L. The position of DC/DC converter was at the Li3 battery side. To prove that the GSA searched the optimal results in the short-distance electric bus case, we selected cases in which the size ratio of LTO:Li3 was 50:50 (both were 36 kWh) and the total price was 35,589 USD. It increased 1.93% compared to the GSA optimal result. The case with the highest price that satisfied the physical constraints was 36,264 USD, which was increased by 3.86%, where the electric capacities were 40 and 32 kWh, respectively. Table 8 demonstrates the optimal sizing of the electric sports car. The case of SC/Li3 combination was the best solution, as the case of LTO/Li3 was unable to satisfy the required maximum output power defined in Section 3. The HEES cost was calculated to be 45,315 USD. For the electric capacity, the optimal sizing was that the capacity of the SC module (as the high-power source) was 0.578 kWh, where the capacity of the Li3 battery (as the high-energy source) was 99.422 kWh. The size ratio was 0.58:99.42, where the total required capacity was 100 kWh. With the optimal sizing, the maximum output power of the SC was 871 kW, while that of the Li3 battery was 422 kW. The combined power exceeded the required power of 693 kW. The HEES mass was 595 kg, while the occupied volume was 298 L. The position of DC/DC converter was at the Li3 battery side.
Our simulation results showed that the designed GSA of the HEES will effectively search the optimal sizing and proper position of the DC/DC converter for various types of EVs.

Conclusions
The conclusions of this study are summarized as follows: (1) Performance indices and vehicle requirement: the average specific price, specific energy at a constant volume, specific energy at a constant mass, and specific power at a constant mass for three energy sources, SCs, LTO batteries, and Li3 batteries, were studied. The vehicle requirements including the maximum output power, vehicle acceleration, climbability, and maximum speed as the design constraints were formulated. (2) The GSA of the HEES: the global search method with four steps including parameter input, definition of cost function and constraints, global search, and optimal results was developed. (3) Numerical analysis of five types of EVs: by (1) and (2), five types of EVs including the electric sedan, electric forklift, long-distance electric bus, short-distance electric bus, and electric sports car were analyzed for optimal sizing of the HEES and the optimal position of the DC/DC converter with the lowest cost function. For the electric sedan, the optimal sizing was LTO+Li3 with a ratio of 9.9:90.1 and the DC/DC position was at Li3 side; for the electric forklift, the optimal sizing was SC+Li3 with a ratio of 4.5:95.5 and the DC/DC position was at SC side; for the long-distance electric bus, the optimal sizing was pure Li3 with no DC/DC converter; for the short-distance electric bus, the optimal sizing was LTO+Li3 with a ratio of 44.4:55.6 and the DC/DC position was at Li3 side; and for the electric sports car, the optimal sizing was SC+Li3 with a ratio of 0.58:99.42 and the DC/DC position was at Li3 side.
Utilizing these results, the integrated system and control designs of the HEES using the GSA, more applications for green energy sources, and different types of EVs will be studied in the future.