Optimization of the Energy Storage of Series-Hybrid Propelled Aircraft by Means of Integer Differential Evolution

The possibility of realizing full electric or hybrid electric propulsion for aircraft has been considered due to the constant growth in the use of electric technologies in aircraft and the availability of high-power-density electrical machines and converters. In this paper, an optimized design approach is proposed with reference to the optimal trade-off between energy storage system sizing and the fuel mass of a series of hybrid aircraft. The problem is approached using an integer optimization algorithm based on differential evolution and by mixing both the flight mechanics and the electrical issues inherent to hybrid flights. This method has been validated by means of implementing numerical simulations and the results are reported and discussed in the paper.


Introduction
In the last two decades, increasing interest regarding the More Electric Aircraft (MEA) concept has been generated [1] with the aim of reducing fuel consumption, reducing emissions released into the atmosphere, increasing the redundancy and the reliability of onboard power systems, and reducing aircraft noise.
Recently, several research projects have led to the substantial electrification of onboard pneumatic components in order to increase the efficiency of flight systems-which has also introduced redundancies across the subsystems [2]-and the necessity to propose new power system architectures for the reduction of weight, increase of the reliability and to assure the correct performance of electrical drives [3]. Moreover, in order to increase the reliability and decrease the weight, DC-microgrids [4] seem to be the most suitable solution for future aircraft.
The current trend also envisages the use of electrical technology in aircraft propulsion [5][6][7][8] and both fully electric and hybrid propulsion aircraft prototypes have been presented. In the all-electric propulsion concept, the fuel and the thermal engine are completely substituted with energy storage systems (above all, fuel cells and battery) that are able to provide the entirety of the electric power required by the aircraft. The use of these new energy systems on board aircraft has been analyzed widely in the literature. In [9], an aircraft design that relies on fuel cells for propulsion and that is characterized by an unusual motor placement in order to maximize efficiency is presented, while in reference [10], the behavior of batteries that are used for the propulsion of small unmanned aircraft vehicles (UAV) at different temperatures is discussed. Additionally, Tariq [11] presents an overview of the battery systems for the MEA with an analysis of battery management systems for Li-ion batteries.
that the optimal value of the energy and power is achieved, but it is not guaranteed that the number of cells is optimal.
Starting from an assigned flight profile (cruise altitude, cruise speed, cruise range), the aim of this paper is to develop an optimization procedure which can be used during the pre-flight settings and is able to evaluate the optimal integration of the battery storage system and ICE to guarantee the correct operation during the execution of the flight profile.
The main novelty of the proposed optimization procedure regards the possibility to use it as a pre-mission energy management process. In fact, with the possibility to develop a modular battery pack and considering the mission profile, it is possible to calculate the optimal size of the battery pack. Through a new disposition of the module of the pack (each module can be constituted by a certain number of parallel and series cells), it is possible to minimize and reduce the useless weight of the aircraft.
The optimized design approach of the sizing of the energy storage system of a series-hybrid aircraft is detailed and validated in this paper. In the next section, the hybrid propulsion system of the aircraft is defined. In the third section, the dynamic equations of the aircraft profile are detailed and in Section 4 the battery system model is explained. The optimization sizing procedure is proposed in Section 5 and validated by simulation in Section 6 where the results are reported and discussed.

Hybrid Propulsion Scheme
A series-hybrid propulsion scheme is considered in the present work. The advantages of a series-hybrid propulsion are different, for example, the thermal engine can be controlled to let it operate at a fixed point during the recharging of the batteries, with the possibility to keep constant and maximize the efficiency. Furthermore, the series propulsion provides the possibility of having only electric propulsion of the propellers, avoiding the mechanical gear between the engine and the propeller, and increasing the efficiency (in particular when the propellers have a fixed pitch, and in this case, the variation of the speed within a certain range is necessary). Figure 1 shows a scheme of the series-hybrid architecture considered in the paper. Starting from an assigned flight profile (cruise altitude, cruise speed, cruise range), the aim of this paper is to develop an optimization procedure which can be used during the pre-flight settings and is able to evaluate the optimal integration of the battery storage system and ICE to guarantee the correct operation during the execution of the flight profile.
The main novelty of the proposed optimization procedure regards the possibility to use it as a pre-mission energy management process. In fact, with the possibility to develop a modular battery pack and considering the mission profile, it is possible to calculate the optimal size of the battery pack. Through a new disposition of the module of the pack (each module can be constituted by a certain number of parallel and series cells), it is possible to minimize and reduce the useless weight of the aircraft.
The optimized design approach of the sizing of the energy storage system of a series-hybrid aircraft is detailed and validated in this paper. In the next section, the hybrid propulsion system of the aircraft is defined. In the third section, the dynamic equations of the aircraft profile are detailed and in Section 4 the battery system model is explained. The optimization sizing procedure is proposed in Section 5 and validated by simulation in Section 6 where the results are reported and discussed.

Hybrid Propulsion Scheme
A series-hybrid propulsion scheme is considered in the present work. The advantages of a serieshybrid propulsion are different, for example, the thermal engine can be controlled to let it operate at a fixed point during the recharging of the batteries, with the possibility to keep constant and maximize the efficiency. Furthermore, the series propulsion provides the possibility of having only electric propulsion of the propellers, avoiding the mechanical gear between the engine and the propeller, and increasing the efficiency (in particular when the propellers have a fixed pitch, and in this case, the variation of the speed within a certain range is necessary). Figure 1 shows a scheme of the series-hybrid architecture considered in the paper. The architecture is based on a twin-propeller civil aircraft. Each propeller is controlled by a permanent magnet synchronous motor (to increase the reliability, it is also possible to use a motor with more than three phases), a cascade of a DC-DC converter, and a VSI links the battery modules to the motor. A single energy storage system and a single controlled rectifier for both the propellers are considered in the scheme. The reliability of the system can be increased by dividing the battery pack into two parts and using two different rectifiers ( Figure 2). In this way, the propulsion system is not completely compromised when a fault occurs in either the storage system or in one of the converters. The architecture is based on a twin-propeller civil aircraft. Each propeller is controlled by a permanent magnet synchronous motor (to increase the reliability, it is also possible to use a motor with more than three phases), a cascade of a DC-DC converter, and a VSI links the battery modules to the motor. A single energy storage system and a single controlled rectifier for both the propellers are considered in the scheme. The reliability of the system can be increased by dividing the battery pack into two parts and using two different rectifiers ( Figure 2). In this way, the propulsion system is not completely compromised when a fault occurs in either the storage system or in one of the converters.

Dynamics Equations of Aircraft Load Profile
To determine the dimension of the batteries and AVGAS, which must be used for the aircraft flight, it is necessary to model the aircraft trajectories, which can be divided into take-off, climb, cruise, and landing. The energy management strategy considered in this paper is based on the use of an all-electric propulsion during the take-off and climb phases and a mixed electric and series-hybrid propulsion during the cruise. The use of full electric propulsion during the take-off and the climb reduces the noise and pollution emissions of the aircraft, in particular when the aircraft is near the airport and the city; therefore, it is considered a so-called "green" integration of the airport in the city. The sizing of the battery and the determination of the fuel weight depend on the trajectories of the aircraft and their values' influence the performance of the aircraft. The modeling of the aircraft profile can be done using the equation of dynamics.

Take-Off Dynamics Equation
During the take-off, two different types of forces act upon the aircraft: the friction forces due to the weight of the aircraft and the contact between the tires and the ground and the aerodynamic drag. Considering Equations (A1-A3) in Appendix 1, it is possible to determine the length of the take-off using the following equation: During the take-off, it is assumed that the propulsion happens with a constant power equal to the maximum power of the propulsion system.

Climb Phase Dynamics Equation
Considering that the trust vector is parallel to the climb trajectory, the variation of the climb angle during the flight is negligible, and taking into account the dynamic Equation (A4), the duration of the climb [22] can be calculated by integrating the following equations: where P S is equal to:

Dynamics Equations of Aircraft Load Profile
To determine the dimension of the batteries and AVGAS, which must be used for the aircraft flight, it is necessary to model the aircraft trajectories, which can be divided into take-off, climb, cruise, and landing. The energy management strategy considered in this paper is based on the use of an all-electric propulsion during the take-off and climb phases and a mixed electric and series-hybrid propulsion during the cruise. The use of full electric propulsion during the take-off and the climb reduces the noise and pollution emissions of the aircraft, in particular when the aircraft is near the airport and the city; therefore, it is considered a so-called "green" integration of the airport in the city. The sizing of the battery and the determination of the fuel weight depend on the trajectories of the aircraft and their values' influence the performance of the aircraft. The modeling of the aircraft profile can be done using the equation of dynamics.

Take-Off Dynamics Equation
During the take-off, two different types of forces act upon the aircraft: the friction forces due to the weight of the aircraft and the contact between the tires and the ground and the aerodynamic drag. Considering Equations (A1)-(A3) in Appendix A, it is possible to determine the length of the take-off using the following equation: During the take-off, it is assumed that the propulsion happens with a constant power equal to the maximum power of the propulsion system.

Climb Phase Dynamics Equation
Considering that the trust vector is parallel to the climb trajectory, the variation of the climb angle during the flight is negligible, and taking into account the dynamic Equation (A4), the duration of the climb [22] can be calculated by integrating the following equations: where P S is equal to: and defines the specific excess power. The equations of Equation (A4) are referred to an ascendant climb phase and can also be considered for the descendent climb, with a variation of the signs in the equilibrium of the dynamic quantities.

Cruise Phase
The cruise phase is modeled supposing that the speed cruise and the altitude are kept constant; taking into account that during the cruise, the propulsion is electric but the energy is provided both by the batteries and the fuel. Starting from Equation (A5) and using some manipulations of the formula, it is possible to determine the Breguet range [27] equation, which depends on the initial and final value of the weight of the aircraft and also, on the fuel consumption: with η p equal to the product of thermal engine efficiency, propeller efficiency, electric generator and motor efficiency, battery conversion efficiency, and converter efficiency.

Modeling of the Electrical Storage System
The sizing of the energy storage system and the determination of the fuel mass necessary for the optimal thermal engine performance can be done by means of models of the electric components of the propulsion system. The electric motor and power converter can be modeled considering only the power density (kW/kg) and the efficiency, which depend on the load. The power densities of the electrical machines and converters are calculated considering the data of commercially available three-phase permanent magnet machines and converters with a high power density (all of the considered components are liquid cooled), the efficiency is modeled using a polynomial which depends on the required load power and is obtained from the datasheet of the commercial components.

Battery Modeling
The battery storage system is modeled according to [26]: where the meaning of the utilized symbols can be found in Figure 3. n s indicates the number of series cells, while n p is the number of parallel branches. The current of a single cell is I b,m ; the sum of the currents of each branch gives the total current of the battery pack. Without any current, the output voltage is V oc,sm (due to the internal voltage drop, this value decreases when the battery is not zero).
It is assumed that the battery temperature is kept constant (the variation of temperature could be due to the variation in the ambient temperature at different altitudes or to the losses generated during the operations). With this assumption, it is possible to avoid the dependence of battery parameters on the variable temperature, with a reduction of the complexity of the problem. Aerospace 2019, 6, x 6 of 20 It is assumed that the battery temperature is kept constant (the variation of temperature could be due to the variation in the ambient temperature at different altitudes or to the losses generated during the operations). With this assumption, it is possible to avoid the dependence of battery parameters on the variable temperature, with a reduction of the complexity of the problem.

Optimization Procedure
The optimization procedure is mainly based on the use of an integer differential evolution approach (IDEA) [21] with the aim to minimize (or maximize) certain parameters of the aircraft or the physical parameters of the aircraft trajectories. The optimization procedure is realized on more levels; the external level is inherent to the IDEA technique and acts on three variables: the number of parallel cells np and series cells ns, and the weight of fuel. At each iteration of the optimization, the procedure solves another three problems linked to the optimization of the trajectories. Figure 4 shows a flowchart of the procedure.

Pre-Processing to IDEA
Before initializing the procedure, a population of npop possible solutions is generated:

Generation of Mutations
The mutations (that is the new candidate solution) are generated according to the procedure represented schematically in Figure 4.

Optimization Procedure
The optimization procedure is mainly based on the use of an integer differential evolution approach (IDEA) [21] with the aim to minimize (or maximize) certain parameters of the aircraft or the physical parameters of the aircraft trajectories. The optimization procedure is realized on more levels; the external level is inherent to the IDEA technique and acts on three variables: the number of parallel cells n p and series cells n s , and the weight of fuel. At each iteration of the optimization, the procedure solves another three problems linked to the optimization of the trajectories. Figure 4 shows a flowchart of the procedure.

Pre-Processing to IDEA
Before initializing the procedure, a population of n pop possible solutions is generated:

Generation of Mutations
The mutations (that is the new candidate solution) are generated according to the procedure represented schematically in Figure 4.
The parameters ρ 1 , ρ 2 , and ρ 3 are random integer numbers comprised in the range (1,n pop ), α is a random positive real number below 1, and the function "int" rounds the number obtained in the square brackets to the nearest integer. If the obtained solution is positive, a candidate solution is found and processed inside the modeling of aircraft dynamics. The solution of Equations (A1) and (1) permits the calculation of the take-off time and the total energy, E to , required during this phase. It is imposed that the take-off starts at a level of 0 m and is completed at z to,f (about 10.6 m) [26]. The final condition of the battery pack is evaluated by The parameters ρ1, ρ2, and ρ3 are random integer numbers comprised in the range (1,npop), α is a random positive real number below 1, and the function "int" rounds the number obtained in the square brackets to the nearest integer. If the obtained solution is positive, a candidate solution is found and processed inside the modeling of aircraft dynamics.
The solution of Equations (A1) and (1) permits the calculation of the take-off time and the total energy, Eto, required during this phase. It is imposed that the take-off starts at a level of 0 m and is completed at zto,f (about 10.6 m) [26]. The final condition of the battery pack is evaluated by calculating the required electrical power. These values are used as initial values for the next step in the optimization problem.

Solutions of Optimization Problems of Climb Ascendant/Descendent
The optimization of the climb phase is more complex with respect to that of take-off. In fact, in this case, an objective optimization is carried out inside the IDEA. The goal of this optimization is the minimization of the total energy required for the climb phase. The simulated climb phase is assumed to possess a uniform acceleration. Considering the fact that the climb phase commences at zto,f and the final altitude is the cruise altitude, zcr, it is also considered that the final speed must be equal to the cruise speed (in this case, any acceleration during the cruise phase is avoided). In the hypothesis of

Solutions of Optimization Problems of Climb Ascendant/Descendent
The optimization of the climb phase is more complex with respect to that of take-off. In fact, in this case, an objective optimization is carried out inside the IDEA. The goal of this optimization is the minimization of the total energy required for the climb phase. The simulated climb phase is assumed to possess a uniform acceleration. Considering the fact that the climb phase commences at z to,f and the final altitude is the cruise altitude, z cr , it is also considered that the final speed must be equal to the cruise speed (in this case, any acceleration during the cruise phase is avoided). In the hypothesis of uniform acceleration, the altitude climb is divided in a certain number of divisions, in order to obtain that: where ∆v k is the variation of the velocity in each division. The used objective function is the total energy required in the climb: where the unknown variables are: x cl = ∆v 1 , . . . , ∆v n div , ∆z 1 , . . . , ∆z n div , T 1 , . . . , T n div , γ Then, the objective function becomes: In Equation (14), the sum of the energy calculated in each sub-interval considered for the climb phase is given. In particular, the numerator is presented as the product between the propeller trust, the velocity, and the time duration and this product is divided by the global efficiency of the electric propulsion chain. The variation in the altitude and aircraft speed must satisfy the following equations: The first two equations are related to the fact that the final altitude and the final aircraft speed must be equal to the desired cruise altitude and cruise speed; therefore, the optimal value obtained for the altitude and speed variation in each sub-interval must be related to these values. The third equation relates the climb angle to all the dynamic variables.
The inequalities are: The first two inequalities are used in order to limit the variation of the aircraft speed and altitude; the last two ensure that the maximum time for the climb phase and the climb angle are respected.
Starting from the final condition of the battery packs after take-off, the final state of charge and the final value of the open-circuit voltage V oc,sm are calculated. The same approach is used in the descendent phase of the climb.

Solutions of Dynamic Equations of the Cruise
The solution of Equation (A5) gives the operation characteristics of the cruise. In particular, the cruise range, x cr , is used as the objective function in the optimization procedure. According to the control strategy, during the cruise, all the power is obtained by the thermal engine coupled with the electric generator; therefore, it is possible to assume that the state of charge (SoC) of the battery is essentially constant and equal to the final values obtained in the ascendant climb phase, and will be modified only in the descendent climb phase. At each iteration, whether the power chosen for the ICE is sufficient to guarantee the required trust at the cruise speed is determined.
As previously mentioned, the landing phase is not considered in the simulations because a large part of the energy needed to decelerate the aircraft can be provided by the mechanical braking of the propeller.

Check of the Optimization Procedure Constraints and Choice of Objective Function
The output obtained by the solutions of mechanical flights and the electric state of the battery storage system are used to verify some constrains applied to the differential evolution optimization.
Specifically, two different problems are analyzed. In the first problem there is the minimization of difference between the cruise range and a target cruise value; the approach allows for a trade-off between the battery mass and fuel mass, which guarantees a certain cruise range. The adopted problem could be applied to a training course and in the application no fuel reserve is considered. The constrains are inherent to: the respect of the maximum take-off aircraft mass; -the final SoC of the battery; -the maximum duty-ratio between the output voltage battery pack and the required DC-link voltage; -the energy stored in the battery pack must be at least equal to the energy required in the take-off and in the climb phase; -the maximum value of the discharge current of the battery pack in the take-off and climb phases must be below the maximum discharge current of the battery cells used.
The second problem is related to the reduction in the battery costs: the minimization of the energy surplus contained in the battery (this surplus comes from the fact that the modules are defined as integer numbers) reduces the weight of the aircraft (or increases the amount of fuel with an increase in the cruise range) as well as the cost of battery installations.

Problem 2 min n p n s m s,c δ batt
The goal of Problem 1 is linked to the optimization of the fuel weight; it has to provide a trade-off between the battery mass and fuel mass which maximizes the cruise. Problem 2 is related to the reduction in the battery costs: the minimization of the energy surplus contained in the battery (this surplus comes from the fact that the modules are defined as integer numbers) reduces the weight of the aircraft (or increases the amount of fuel with an increase in the cruise range) as well as the cost of battery installations.

Crossover Operations and Check of Maximum Iterations Limit
As presented in [25], the crossover conserves the original characteristics of the population and also, introduces a new diversity, with an effect similar to genetic algorithms and other heuristics methods. The crossover works are as follows: if the generated solution satisfies all the constraints of the optimization problem and, therefore, becomes a possible new solution, the crossover introduces a new probabilistic constraint. That is, the replacement of the determined solution only occurs if a generated number is greater than a certain threshold value. If that happens, and the number of iterations of the system is below the maximum number of iterations, the procedure is finished. If the value of the obtained objective functions is greater than the value of an element of the starting populations, the new populations replace the old elements of the populations. Otherwise, the procedure is tried again. When the maximum number of iterations is reached, the minimum (or maximum, depending on the objective function) of the final populations is the optimal value.

Simulation Results
The performance of the proposed optimization procedure is verified by means of numerical simulations implemented in Matlab ©. The aircraft parameters are reported in Table 1; the parameters are the same as those of Tecnam P2006T [28], where it is supposed that the ICEs are substituted with an electric motor with the same power, and the empty aircraft mass is unchanged. The ICE considered for the hybrid flight has a maximum power less than the total power of the electric motor because it operates only in the cruise phase, when the required trust is usually below the maximum power of the propeller motors. The evaluation of the optimum value of rated power for the electric motor and of the maximum power of the ICE is not considered in this paper.
The simulation refers to the flight profile reported in Figure 5. In Table 2, the main parameters of the considered battery pack are reported. The first simulation reported demonstrates the convergence of the optimization problems. For each problem, other simulations using different cells [29] and two different voltages for the DC-link are carried out.

Simulation of Problem 1
The problem described by Equation (22) is simulated, imposing a maximum value of dmax = 6.5. An initial population of npop = 8 is generated, and the values α, ρ1, ρ2, ρ3 and the threshold for the crossover operations are obtained using the random generators of Matlab. The number of mutation

Simulation of Problem 1
The problem described by Equation (17) is simulated, imposing a maximum value of d max = 6.5. An initial population of n pop = 8 is generated, and the values α, ρ 1 , ρ 2 , ρ 3 and the threshold for the crossover operations are obtained using the random generators of Matlab. The number of mutation substitutions is set as 100. The cruise range is fixed as 180 km. The final optimized values obtained for 11 Ah batteries with a continuous discharge rate of 8 C and a DC-link voltage of 270 V, are reported in Table 3. Figures 6 and 7 show the trends of mutation substitutions of the battery module numbers (both in series and parallel) of the fuel mass and the cruise range. It is possible to note that the convergence of the method occurs at approximately 40 iterations. In Figure 8, the trends in SoC and the maximum current in the climb phase are reported.
In Table 4, the results obtained for different cells (with equal nominal voltage and parameters d 0 and d 1 ) and for an inverter DC-link of 270 V and 540 V are shown. of the method occurs at approximately 40 iterations. In Figure 8, the trends in SoC and the maximum current in the climb phase are reported. In Table 4, the results obtained for different cells (with equal nominal voltage and parameters d0 and d1) and for an inverter DC-link of 270 V and 540 V are shown.   current in the climb phase are reported. In Table 4, the results obtained for different cells (with equal nominal voltage and parameters d0 and d1) and for an inverter DC-link of 270 V and 540 V are shown.      It is evident that the cruise range is closer to the fixed value, but the final SoC is higher than that in the case of the DC-Link at 270 V. This happens because the number of cells is increased to respect the maximum duty cycle imposed for the DC-DC converter. It can be concluded that the excess stored energy can also be used for the reduction of pollutant emissions and to increase the cruise range.

Simulation of Problem 2
The problem described by Equation (18) is simulated, imposing a maximum value of d max = 6.5 and the number of iterations is 100. The parameters for the IDEA are calculated in the same manner as in Problem 1. The optimal solution obtained for 11 Ah batteries with a continuous discharge rate of 8 C and a DC-link voltage of 270 V is indicated in Table 5. Figures 9 and 10 show the trends of the mutation substitution of battery module numbers in series and in parallel, and the variation of the energy stored in excess and in the cruise range.
The results of simulations demonstrate that, in Problem 2, the type of optimization used enables the cruise range to be increased while reducing the energy stored in excess in the battery pack. This is due to the relationship between the fuel mass and the energy stored in the battery: if the fuel mass increases, more energy is needed for the take-off and climb phases, resulting in greater utilization of the battery and an increase in the cruise range as shown in Figure 11.
In Table 6, the results obtained for Problem 2 are reported and related to different types of cells and two voltage levels of the DC-link. With respect to the results of Problem 1, in this case it is possible to note that the use of a higher value for the DC-link improves the performance of the aircraft, especially for the maximization of the cruise range. Furthermore, the minimization of the residual energy stored in the batteries is also evident for the DC-link of 540 V.        As mentioned in Table 1, in the mission profile, the presence of at least one passenger must be guaranteed. From the calculation of the total mass required by the mission profile it is possible to obtain that for both; the best payloads are obtained with a DC-link of 270 V and with a battery of 16  As mentioned in Table 1, in the mission profile, the presence of at least one passenger must be guaranteed. From the calculation of the total mass required by the mission profile it is possible to obtain that for both; the best payloads are obtained with a DC-link of 270 V and with a battery of 16 Ah. Specifically, the maximum payload for Problem 1 is about 313 kg and is about 211 kg for Problem 2. This confirms the respect of all the constraints (electrical and mechanical) and highlights the possibility to effectively optimize the pre-flight settings of the aircraft.
Regarding the convenience of the hybrid technologies, it is possible to highlight the following results: with the series topologies adopted, the emissions during the take-off and climb phases are eliminated; these are the parts of the flight that occur close to built-up areas, such as cities; -to date, an economically beneficial way to replace a traditional piston engine aircraft with a hybrid propelled aircraft does not exist. In fact, excluding the cost of the electric powertrain and considering only the battery costs of $176/kWh [29], a cost of $30,000 regarding implementing a hybrid propelled aircraft for the first problem and a cost of $34,500 for the second one is obtained. Evidently, the introduction of a new battery type or the rapid increase of the battery energy density, will dramatically increase the viability of hybrid/all-electric technologies.

Verification of the Results Through the Real-Time Simulations of the Battery Pack
Real-time simulations of the performance of the battery pack during the take-off and climb phases are carried out with Simulink©. The battery cells are modeled according to [27] and using the parameters of the Kokam cells [30]. The simulations performed are related to the results obtained in Problem 1 and for the cases of an 11 Ah cell ( Table 4). The results for a DC-link voltage of 270 V are depicted in Figure 12, while the case of 540 V is shown in Figure 13.
As it is possible to note in both the cases, the final values of SoC are very close to the value obtained from the optimization. In this way, real-time simulations demonstrate the effectiveness of the proposed optimization procedure. introduction of a new battery type or the rapid increase of the battery energy density, will dramatically increase the viability of hybrid/all-electric technologies.

Verification of the Results Through the Real-Time Simulations of the Battery Pack
Real-time simulations of the performance of the battery pack during the take-off and climb phases are carried out with Simulink©. The battery cells are modeled according to [27] and using the parameters of the Kokam cells [30]. The simulations performed are related to the results obtained in Problem 1 and for the cases of an 11 Ah cell ( Table 4). The results for a DC-link voltage of 270 V are depicted in Figure 12, while the case of 540 V is shown in Figure 13.
As it is possible to note in both the cases, the final values of SoC are very close to the value obtained from the optimization. In this way, real-time simulations demonstrate the effectiveness of the proposed optimization procedure.

Conclusions
All electric aircraft and the hybrid electric aircraft are currently feasible technologies which can be applied to UAVs and passenger aircraft. A crucial problem is the correct integration and sizing of the storage systems utilized on board. This factor influences the performance of the aircraft in terms

Conclusions
All electric aircraft and the hybrid electric aircraft are currently feasible technologies which can be applied to UAVs and passenger aircraft. A crucial problem is the correct integration and sizing of the storage systems utilized on board. This factor influences the performance of the aircraft in terms of battery and fuel costs, flight distance range and flight time, and safety. An optimization procedure for a series-hybrid aircraft has been proposed in this paper based on an integer differential evolution approach, which taking into account the electrical and flight mechanical behavior of the aircraft, gives the optimal trade-off between the battery size and the mass of fuel. The series-hybrid architecture provides a take-off and climb phase (both ascendant and descendent) in the electric mode, with the goal of reducing the noise emission and pollution in the phase when the aircraft is nearest to the airport. The sizing optimization procedure has been applied to a real case and validated by implementing numerical simulations in Matlab, and the results have been verified by real-time scenario simulations. The results demonstrate the fast convergence of the method and the importance of the main electrical parameters (such as the DC-link voltage) in the application. The proposed procedure respected all the electrical and mechanical constraints and also determines the minimum payload which must be guaranteed in the aircraft. Future work will include the optimization of the power-train performance with the integration of the energy store optimization and the optimization of the drives control and the propulsion motor sizing procedure. Cruise range (m) x to Take-off range (m) z cr Cruise altitude (m) z to,f Take-off final altitude (m) γ Climb angle (rad) δ batt Energy density of battery module (J/kg) ∆v k Speed variation during climb interval (m/s) ∆z k Altitude variation during clim interval (m) ∆t k Duration of k-th climb phase sub-interval (m) η batt Battery efficiency η VSI VSI efficiency η dc DC-DC converter efficiency η mot Electric motor efficiency µ Friction coefficient Air density (kg/m 3 )

Take-Off Dynamics Equations
Considering that the vector which determines the trust force is parallel to the ground, the mathematical model of dynamics equations of the aircraft is, (A1) In this case, the variable W contains the constant weight of the battery and the weight of AVGAS, which is also kept constant because the take-off is realized only with the electric propulsion. The flight begins when the final velocity is reached: its value is 20% greater than the stall speed, defined as, Regarding the drag force, this value is calculated using the well-known formula: (A3)

Climb Phase Dynamics Equations
The dynamics equations of climb phase are Cruise Phase Dynamic Equations The dynamics equations of the cruise are