Inﬂuence of Novel Airframe Technologies on the Feasibility of Fully-Electric Regional Aviation

: The feasibility of regional electric aviation to reduce environmental impact highly depends on technological advancements of energy storage techniques, available battery energy density, and high-power electric motor technologies. However, novel airframe technologies also strongly affect the feasibility of a regional electric aircraft. In this paper, the inﬂuence of novel technologies on the feasibility of regional electric aviation was investigated. Three game-changing technologies were applied to a novel all-electric regional aircraft: active ﬂow control, active load alleviation, and novel materials and structure concepts. Initial conceptual design and mission analysis of the aircraft was performed using the aircraft design framework SUAVE, and the sensitivity of the most important technologies on the aircraft characteristics and performance were studied. Obtained results were compared against a reference ATR-72 aircraft. Results showed that an all-electric aircraft with airframe technologies might be designed with the maximum take-off weight increase of 50% starting from the battery pack energy density of 700 Wh/kg. The overall emission level of an all-electric aircraft with novel technologies is reduced by 81% compared to the ATR-72. On the other hand, novel technologies do not contribute to the reduction in Direct Operating Costs (DOC) starting from 700 Wh/kg if compared to an all-electric aircraft without technologies. An increase in DOC ranges from 43% to 30% depending on the battery energy density which creates a signiﬁcant market obstacle for such type of airplanes. In addition, the aircraft shows high levels of energy consumption which concerns its energy efﬁciency. Finally, the sensitivity of DOC to novel technologies and sensitivities of aircraft characteristics to each technology were assessed.


Introduction
Significant climate changes and potential environmental impact due to increased transportation in the near future have motivated many industries to focus on reducing CO 2 and NOx emissions. As a major transportation method, the aviation industry also follows the trend to reduce the emission of new generations of aircraft. Improvements in airframe and engine technologies increase aircraft efficiency and reduce their emission. However, a potential increase in air transportation may still lead to an increase in overall CO 2 and NOx emissions. Under Flightpath 2050 [1], The European Commission has set a future challenge for the new generation of aircraft to reduce their total emission. Figure 1 shows three schematic trends: if no advancements in aircraft technologies are present, if currently feasible technology advancements are achieved, and if novel technologies are introduced. This challenge leads to developing alternative environmentally-friendly energy sources as one of the main solutions for environmental impact reduction. Aircraft electrification has become one of the most popular approaches to reduce aircraft emissions. Today, many companies work in various directions to make electric flights available: improve battery energy capacity, modify and develop new propulsion systems, and introduce new aircraft configurations more applicable for future electrification. At the moment, research is highly focused on hybrid-and all-electric propulsion systems and aircraft configurations to reduce the environmental impact produced by aircraft. Friedrich introduced a simulation technique for a single-seat hybrid-electric demonstrator, designed the aircraft, and performed a scaling analysis to determine a large fuel savings impact for the small-and mid-scale sector of aircraft [2]. Hamilton investigated the effect of hybrid-electric aircraft from the operational perspective to reach an optimal aircraft operation with battery energy density constraints [3]. A lot of work is also dedicated to the introduction of distributed electric propulsion (DEP) energy networks and aircraft concepts featuring DEP. Kim summarized major contributions towards the development of DEP aircraft and relative technologies [4]. Finger focused on sizing methodologies, and aircraft design of future general aviation aircraft that could be not only hybrid-electric but also fully electric [5][6][7]. De Vries also developed initial sizing methodologies for hybrid-electric aircraft for conceptual aircraft design but introduced effects of DEP at early conceptual design stages [8,9]. Pornet also introduced a sizing methodology for hybridelectric aircraft and compared performance between the conventional reference aircraft and its hybrid version [10]. Sgueglia introduced an MDO framework for hybrid-electric aircraft with coupled derivatives and performed various optimization analyses to demonstrate the capabilities of the framework [11]. The work of Hepperle [12] addressed the potentials and limitations of all-electric aircraft specifically and introduced aircraft modifications that could potentially improve the aircraft performance at given battery energy densities. The ultimate goal of aircraft electrification is to achieve a fully-electric flight at adequate aircraft weight characteristics to maximize the emission reduction at the aircraft level. However, the availability of all-electric aircraft is limited to the General Aviation sector due to low battery energy density compared to the Jet-A fuel and relatively low maximum power capabilities of modern electric motors. Multiple all-electric concepts have already been introduced or are being developed at the moment. From existing airplanes, Pipistrel has already certified a twin-seat all-electric aircraft that has a maximum speed of 185 km/h and a maximum range of 139 km [13]. Bye Aerospace is currently performing flight tests on a one-seat general aviation aircraft that will be certified under the FAR Part 23 category. The eFlyer has a maximum cruise speed of 250 km/h and endurance of three hours [14]. In 2016, Siemens tested an all-electric energy network by retrofitting the Extra 300 aerobatic aircraft and set two world records [15]. The Equation Aircraft has also demonstrated an all-electric twin-seat amphibious aircraft with the empennage-mounted engine that can fly up to 240 km/h and up to 200 km [16] Airbus has also developed several versions of an experimental all-electric E-fan aircraft to test the capabilities of an all-electric aircraft [17]. The E-Fan concept had multiple versions, from a technology demonstrator to a production variant. MagniX has introduced an all-electric propulsion system and first retrofitted De Havilland Beaver and then 208B Cessna Grand Caravan. The electrified Caravan became the largest all-electric airplane until today [18]. It must be noted that the list of existing all-electric aircraft companies is incomplete and more companies exist and are at different design stages.
To achieve all-electric regional aviation faster, not only the energy source and propulsion system technologies must be improved, but also advanced airframe technologies need to be considered. Such novel technologies may significantly increase the aircraft performance characteristics and overall efficiency that will enable earlier integration of allelectric commercial aircraft into the market before. Liu introduced initial estimations of the impact of novel technologies for a range of aircraft from the short-range to the long-range [19]. Results showed that novel airframe technologies may significantly improve aircraft energy efficiency.
Under the Excellence Cluster Se 2 A (Sustainable and Energy Efficient Aviation), three energy-efficient aircraft are to be designed to cover the majority of commercial aircraft operations. Figure 2 shows three sample new energy-efficient airplanes: a short-range propeller aircraft that shall be applicable to the ATR-72, the medium-range aircraft having similar mission requirements to the Airbus A320, and the long-range aircraft to cover ranges similar to the Boeing B777. The initial cluster requirement for the short-range aircraft was to investigate the feasibility and availability of an all-electric short-range aircraft if novel airframe and propulsion technologies are introduced. A few novel technologies have been considered for the design of future regional aircraft. Hybrid laminar flow control (HLFC), load alleviation, and advanced materials and structure concepts. In addition, high-temperature superconducting motors and high-energy capacity batteries have been considered for the present research. An aircraft with similar top-level requirements as ATR-72 is designed considering the mentioned technologies and a battery-based full electric propulsion system.
The present work is divided into multiple sections. Section 2 describes novel technologies implemented in the new short-range aircraft. Section 3 describes methodologies, models, assumptions, and software used to perform novel technology assessments to design the aircraft. Section 4 describes the aircraft's top-level requirements (TLRs). Section 5 describes the aircraft concept selection and its initial sizing. Section 6 presents multidisciplinary design optimization studies of the SE 2 A aircraft to refine the design and make important design decisions. Finally, Section 7 assesses the influence of novel technologies on aircraft availability and the strength of each technology impact.

Hybrid Laminar Flow Control
The generation and extension of laminar flow over the aircraft surface significantly affect the overall aircraft drag, reducing aircraft weights, fuel burn, and operating costs. Preliminary estimations of an aircraft that features extended laminar flow along the wing, empennage, and fuselage demonstrated significant reduction (up to 50%) in overall drag and proved the importance of laminar flow control (LFC) [20]. Figure 3 shows the effects of the new wing design featuring the LFC and compares overall drag to the reference aircraft.
However, natural laminar flow (NLF) is limited along the shape, so to extend it even more, hybrid laminar flow control (HLFC) is required: the laminar boundary layer is not only extended using NLF design but also using the active suction technology. In this technology, the air is sucked from the aircraft's outer surface to delay the transition of the boundary layer and enable substantially higher percent laminarization compared to conventional surfaces. The skin of each aircraft component is split into two segments: a porous sheet and an inner sheet that supports the outer sheet. The inner sheet has orifices that suck the air from the boundary layer and delay transition. In each chamber, an individual pressure is adjusted by the throttle orifices, so that the pressure difference between the outside and the chamber delivers the locally desired amount of mass flow through the surface. Figure 4 shows a schematic image of the skin layout and the HLFC system for a wing section. The applied technology in this project is based on [20,21], which describes numerical approaches with active laminar flow control and also describes current progress in this technology.

Load Alleviation
Load alleviation introduces various techniques to reduce the bending moment experienced by the aircraft and that have a passive or active nature. Reduction of maximum bending moment enables the design of lighter wings for lower limit load factors, which will improve aircraft fuel efficiency.
Passive load alleviation solutions consider nonlinear stiffness material design [22], viscoelastic damping design [23], new structural concepts [24,25], and local morphing structures [22]. Nonlinear stiffness materials may improve the load distribution on the wing under low load cases (0.5 g to 1.5 g) and improve performance efficiency at those load cases. Finally, both aeroelastic tailoring and local morphing structures aim to extend the aeroelastic design space. Morphing is considered into two scenarios: deliberate structural non-linearity that affects the wing deformation and reduces effective angle-of-attack of the wing section under the load and change in airfoil shape to achieve load reduction [22].
The wing active load alleviation uses different types of flow control over the wing to achieve a more favorable wing load distribution and reduce the wing bending moment. Previously, researchers have approached the design of active load alleviation systems in different ways. Rossow et al. [26] investigated the feasibility of an aircraft featuring load alleviation technologies; Sun et al. [27] looked at active load alleviation from the control system perspective; and Zing et al. [28] demonstrated the experimental results of the high aspect-ratio wing wind-tunnel test with active load alleviation while using piezoelectric control. The work of Fezans et al. focused on novel sensors and control systems technologies to enable rapid and robust active load alleviation [29][30][31][32]. Finally, active load alleviation could also be reached by fluidic or micro-mechanical flow actuators that change the load distribution along the wing and reduce the bending moment [33].

New Materials and Structure Concepts
Novel structural concepts and materials are being developed to improve the aircraft structure in terms of stiffness and weight. Bishara et al. [34] describe advanced structural design with the integration of active flow control, which is directly applicable to the presented research. Under the excellence cluster, the reduction in the airframe weight is assumed to be reached by application of CFRP thin ply laminates. In addition, new structural concepts must be applied to satisfy HLFC requirements [22].

High-Power Superconducting Motors
One of the major problems in hybrid-and all-electric commercial aircraft deals with limitations of maximum available electric motor power. The motor resistance increases rapidly with the increase of its power, which rapidly diminishes the engine efficiency. One solution to preserve the motor efficiency and even increase it is to use superconducting motors. Such motors can be synchronous, trapped flux, and fully superconducting. They may also operate at low and high temperatures. To maintain required temperatures inside the motor, a cryogenic cooling system is required. High-temperature superconducting (HTS) motors can achieve high power-to-weight ratios with high power densities without major weight losses compared to low-temperature superconducting motors. Ref. [35] provides a descriptive survey of HTS motors and their future trends. For the present research, a high-temperature, fully superconducting motor was considered.

High Energy Capacity Batteries
Currently, the gravimetric energy density of modern batteries generally does not exceed 300 Wh/kg [36]. Such energy densities may be suitable for ultra-light or limited general aviation all-electric aircraft. However, they are insufficient to enable not only all-electric aircraft of the size of ATR-42 and larger, but even the hybrid-electric version of such aircraft is infeasible at the moment since the energy capacities are around 50 times less compared to kerosene-based fuels. On the other hand, the research of high-capacity batteries makes progress. The latest laboratory research demonstrated battery energy capacities of 800 Wh/kg on the cell level and 500 Wh/kg on the pack level for military applications [36]. In addition, the National Academy of Sciences predicts battery pack energy density to achieve 400-600 Wh/kg at Technology Readiness Level (TRL) 6 in the next 20 years (roughly, by 2035) and be commercially available in 30 years [37]. If the prediction becomes real and the linear trend remains, then batteries may reach TRL 6 with energy densities of 550-850 Wh/kg and be commercially available by 2060. To investigate potential scenarios for the far future and examine battery energy impact on the aircraft feasibility, it is assumed that the pack energy densities may range between 500 and 1100 Wh/kg.

Initial Aircraft Sizing
The conceptual design was performed using various tools. OpenVSP [38] and CATIA were used for the aircraft geometric modeling, and SUAVE [39] was used for more defined aircraft sizing, performance, and mission analysis. To have more capabilities from the aircraft performance analysis standpoint, SUAVE was extended to have classical performance analysis methods used in [40][41][42][43].
The aircraft sizing within SUAVE is performed iteratively. First, initial geometric specifications based on TRLs and any available information are input into the SUAVE.
The information includes the initial wing planform, its section properties, flaps characteristics, empennage and fuselage geometries, the definition of the propulsion system, and all required components such as batteries, gearboxes, power management systems (PMAD), cables, etc. Finally, initial guesses of aircraft weight are input. Then, the first iteration of constraint analysis using equations provided by Gudmundsson [41] is performed. The constraint diagram block includes rapid estimation methods for gas-turbine, piston, and electric aircraft power lapse with altitude and has options to either estimate required aerodynamic properties for the constraint diagram or have them as fixed inputs based on historical trends. Based on the design point selection, the wing loading and power-to-weight ratio are used to run the SUAVE mission analyses to estimate the aircraft performance and its required weight. Then, obtained maximum take-off mass (MTOM) is compared to the initial guess and updated if the tolerance is not reached. In addition, parameters such as the minimum drag coefficient (C Dmin ), Oswald efficiency (e), and maximum lift coefficient (C Lmax ) for clean and flapped configurations are input into the constraint diagram again to update all constraint curves and run the loop again. When the tolerance is reached, the program moves to the aircraft performance block to obtain performance plots for the given aircraft. Particular care was taken to size the wing flaps to meet take-off and landing field length requirements. Methods of Torenbeek [42] and Roskam [43] were implemented within SUAVE to analyze various types of leadingand trailing-edge devices. The empennage sizing within SUAVE is based on the fixed tail volume ratio and then updated and corrected separately based on the desired CG envelope. Performance analyses within SUAVE included take-off, all engines operative (AEO), and one engine inoperative (OEI) climb, cruise, descent, and landing. Figure 5 shows the sizing process within SUAVE.

Energy Network
In the SUAVE analysis framework, a propeller electric energy network with HTS motors was implemented. Figure 6 shows the all-electric energy network layout. Propeller modeling was performed using the cubic spline method described by Gudmundsson [41]. Propeller thrust can be described with a cubic spline based on design operating conditions. The thrust is defined by: where A, B, C, D are the coefficients required to fit the curve accurately. To obtain those coefficients, a system of four equations is established. The first equation is the equation for the static thrust. The second equation is derived for the cruise speed and the desired propeller efficiency. The third equation is the derivative of the curve at the cruise speed, which is equal to zero. Finally, the last equation is derived for the desired maximum cruise speed. The final system of equations for the thrust as a function of cruise speed becomes: Solving Equation (2), a general thrust expression is found, which then is used to find the propulsive efficiency of the propeller. Then, knowing efficiencies of components of the energy network shown in Figure 6, the total power at the battery is found to be used during the SUAVE mission analysis.
The propeller was sized using a combination of propeller diameter trends provided by [44] and the sizing method described by Torenbeek [42] where the required propeller RPM depends on the tip Mach number was used to estimate the maximum propeller RPM.
The assumption of the battery gravimetric energy density remains one of the most influential assumptions for aircraft design. Since the availability of a full-electric regional aircraft is challenging to achieve, an optimistic scenario of the energy density at TRL 6 was assumed. This way, if the aircraft is feasible, the task will be to investigate more realistic bounds of minimum energy densities. Otherwise, the availability of the all-electric airplane is highly unlikely. Consequently, the average optimistic pack energy density (ω B ) of 700 Wh/kg was assumed. To ensure safe battery operation, the state-of-charge of 20% was assumed.
High-temperature, fully superconducting motors were used for the airplane since the required power frequently reaches the order of Megawatts. All superconducting components also have cooling systems to enable desired operating temperatures. The cooling system power loss was estimated using the survey provided in [45] and assuming the Carnot efficiency ( ) of 30%. Table 1 summarizes all assumptions considered for the gas turbine engine and an electric propulsion network.

Weights Estimation
Weights estimation analysis was performed using the FLOPS [47] method programmed within SUAVE. To improve the accuracy of high aspect ratio wing weight estimation, FLOPS was coupled with the physics-based wing weight estimation tool EMWET [48]. The EMWET analysis was performed for a metal wing structure and was corrected by a factor based on the composite wing technology level. Loads used for the EMWET weight estimation were obtained within SUAVE using coupling with AVL [49] Vortex-Lattice method software. In the case where the aircraft wing span reaches Class-C airport regulation limits, an additional weight penalty function due to the wing folding was introduced, so the aircraft can maneuver at the airport and maximize its flight efficiency. The estimation method is based on [50] and accounts for the weight of the insert, the folding mechanism, and the pin to fix the wing. The wing weight corrected by the weight of the folding mechanism is defined by m wherem f old w is the mass of the wing with the folding mechanism,m non f old w is the mass ratio without the folding mechanism,m ins is the mass ratio of the folding insert,m mech f old is the mass ratio of the folding mechanism inside the insert, andm pin f old is the mass ratio of the folding pin. The insert mass ratio was defined based on the curve-fitting of three separate formulations available from [50]: folding at 32%, 48%, and 64% of the span wherem PWS is the mass ratio of the section outboard of the folding pivot point defined bȳ where C is the chord length, b is the wing span, t is the wing thickness-to-chord ratio, and y is the spanwise location of the folding insert. Subscripts r, t, and y correspond to the wing root, tip, and the folding insert. Finally, folding mechanism and pin masses are defined bȳ Figure 7 shows the wing weight penalty as a function of the wing joint location along the semi-span. HTS motor weight estimation was based on relations defined by Stuckl [35] and account for both the motor and the cooling system weight. The motor power-to-weight ratio is defined by where P Mot is the motor power in kW and the power-to-weight ratio P/W is defined in kW/kg. The cryogenic cooling system weight is defined by m cooling = 27.5P loss · e −1.225(log 10 (P loss )) where P loss is the cooling system input power loss in kW and m cooling is the cooling system weight in kg. Propeller weight was estimated using a relation described in [51] where the propeller weight depends on the number of blades, the activity factor, the motor power, propeller maximum RPM, and its material. A composite propeller with six blades was considered for the present design.
The gearbox weight was estimated using relations provided by [52] and is defined by where P is the power at the gearbox, K is the gearbox technology factor equal to 26 for future gearboxes, RPM motor is the motor power assumed at 4800 RPM, the RPM motor , and the m gearbox is the gearbox mass in kg.

Novel Airframe Technologies Assumptions
At the conceptual design stage, novel technologies have been integrated into the design process in terms of coefficients and correction factors. Laminar flow control is estimated by modifying the flow transition location along the wing profile. The desired goal within the SE 2 A for the flow transition is set to 70% for the wing if no disturbances such as the propeller wake or the location of other components is present. Otherwise, specific transition assumptions are prescribed. Wing folding also affects the application of the HLFC. The folding location has four general locations. If the wingspan is less than 36 m (Class C airport limitation), then no folding is applied, and the wing features the HLFC along the complete span. If the wingspan is longer than 36 m but not longer than 110% of the Part C requirement, then the tip folding is applied, the wing has HLFC until the wing folding location, and the folding tips feature natural laminar flow (NLF). If the folding location occurs at any position which splits the aileron, then the folding position is automatically moved to the position where the flap ends and the aileron start, so the aileron design is significantly simplified. The portion of the wing after the folding joint does not feature HLFC and only has NLF. Finally, if the wingspan is such that the aileron can be located at the folding portion of the wing, the folding joint is fixed to 18 m, and the flap is limited to the outboard position of 18 m. HLFC is also applied before the folding joint. The transition location of the NLF portion of the wing is assumed at 45% [53]. The fuselage laminar flow shall be preserved until the wing-body fairing, which is specified based on the wing position, which is determined by the center-of-gravity (CG) envelope and the empennage characteristics. The active load alleviation effect is estimated by reducing the aircraft limit load factor from the required 2.5 [54] to 2.0. For the load alleviation, it is also assumed that regulations of CS-25 may allow lower limit loads if sufficient maturity of the technology is reached. Advanced structure configurations and material effects are assumed to reduce the empty weight by 19% compared to the metal structure.

Top-Level Requirements of the SE 2 A SR Aircraft
Top-level requirements have been set to match the reference ATR-72 aircraft. Table 2 shows the summary of the SE 2 A SR TLRs.  Figure 8 shows a mission profile considered for the design. The mission includes the main flight segment, where the airplane cruises at 7300 m, and the required reserve segment which consists of the divert segment and the 30 min hold. In addition, 5% battery energy contingency was assumed after the reserve mission.

Initial Aircraft Design Concept Selection
In this work, two concepts were considered for the initial aircraft design. The first concept was the top-wing conventional configuration with wing-mounted engines. The battery was split between the fuselage belly and the engine pylons: one segment was located before the landing gear, and the other one was located after the gear. The battery distribution was based on the CG envelope to satisfy all possible passenger loading cases. Battery modules feature a 'quick-swap mechanism, so a mechanic can switch batteries between the flights using a mechanism similar to a high loader. The second configuration is the low-wing configuration with aft-mounted engines. Such configuration, unlike the first configuration, has a fully clear wing which enables laminar flow control along the entire wingspan. Engines are mounted at the aft of the fuselage. The vertical location of the engines approached their maximum to minimize the adverse effect from the wing at high angles of attack. Batteries also feature the 'quick-swap mechanism. Figure 9 shows the described configurations schematically modeled in OpenVSP. Each configuration has its assumptions, features, and limitations. Assumptions are based on research goals of the SE 2 A research groups that focus on the research and development of each technology separately. The first concept has a limited range of laminar flow along the wingspan. The region which is subjected to the propeller wash is assumed to have a fully turbulent flow. The remaining wing portion has 70% laminar flow due to the application of the HLFC. To have a favorable CG envelope, the location of the wing was set to 43% of the fuselage length. The second concept has limited laminar flow control before the landing gear at the wing pressure side. The remaining pressure side and the suction side are 70% laminar. It was also possible to shift the wing further aft to 53% of the fuselage length so that extra 10% could be laminarized compared to the conventional configuration.
Comparison among the concepts was performed using the sizing methodology as was described in Section 2. Due to the presence of various novel technologies and, as a consequence, increased aircraft development, and operating costs, the concept selection was based on the magnitude of Direct Operating Costs (DOC). DOC was calculated using the method presented in [55] where the DOC is divided into energy, maintenance, capital, crew costs, and fees. 2121 flight cycles per year were assumed for both aircraft. For the battery pack, 2000 cycles were assumed. The electricity price was taken from [56] for the year 2050 and is equal to 0.118 EUR/kWh that year. The airframe price was assumed similar to the ATR-72 [55] and is equal to 1595 EUR/kg. The remaining parameters for all required cost components were taken from [55]. A major concern is related to the modification of DOC due to novel technologies. It is unknown how much the DOC will increase due to an early stage of all technologies presented in this work. However, the complexity related to the introduction of the porous skin, the suction system and its structure, which is subject to maintenance complexities, NLF front portion, integration of various load alleviation technologies, and the structure to support the wing laminarization and load alleviation will increase maintenance costs by an order of magnitude. As an optimistic assumption, a scaling factor of 2 was used for the maintenance DOC gain. After the introduction of the maintenance DOC gain, the total DOC formulation becomes DOC Total = DOC Energy + DOC Crew + k g DOC Ma + DOC Cap + DOC Fees (11) where DOC Energy are costs of energy, DOC Crew are crew costs, DOC Ma are maintenance costs, DOC Cap are capital costs, DOC Fees are costs of fees, and k g is the maintenance cost gain factor. Finally, all costs were calculated in 2020 EUR. Due to the presence of an all-electric propulsion system, a conventional assumption of the best design point at the lowest power-to-weight ratio may not hold, as was previously shown by Finger [5] where the best design points were sensitive to the energy type of networks and could show points of minimum take-off or fuel weight at wing loadings larger than the one with minimum power-to-weight ratio. To ensure that a good initial guess for each concept is selected, an aspect ratio and wing loading trade study were performed using the initial sizing procedure described in Section 3. The procedure has multiple steps. First, the initial aircraft sizing is performed to find the design point corresponding to the maximum possible wing loading. Assuming that the selection of an arbitrary design point will not change the constraint diagram significantly, a set of sample wing loadings ranging from 60% of the maximum wing loading to the maximum possible wing loading were selected. Figure 10 shows the sweep of possible design points along the constraint diagram boundary. The aircraft is then sized for each selected wing loading using SUAVE, and the DOC is estimated. This approach is repeated for a selected set of aspect ratios to determine the combination of the wing loading and aspect ratio that may give the lowest desired aircraft characteristics important for the decision making. A sample set of weights and DOC for a low-wing configuration as a function of the wing loading and the AR is shown in Figure 11. Multiple trends can be observed from the Figure 11. First of all, the battery weight generally reduces with the AR, while the operating empty weight (OEW) and maximum take-off mass (MTOM) trends do not follow such trends. Moreover, DOC and weights show discontinuities for some aspect ratios and particular wing loadings. Points of the curve discontinuity correspond to the folding penalty, which is introduced starting from a specific wing loading.  Table 3 shows the summary of the high-wing and the low-wing configurations that have the lowest DOC. Based on the given results, the high-wing configuration has larger weights due to more battery required to complete the mission with less laminar flow.
Consequently, the airframe with less possible drag is more favorable. Consequently, the low-wing configuration is more efficient from the performance standpoint.  Due to an unconventional configuration of the low-wing configuration, a concern related to the aircraft stability and control may occur. Coe [57] has investigated the stability and control characteristics of an aircraft with aft-mounted propeller engines. Results demonstrated a reduction of longitudinal stability for this configuration due to the propeller-pylon interaction. However, minor destabilization can be compensated by a minor increase of the horizontal tail. Lateral stability is not significantly affected by aftmounted propellers and may improve directional control due to the boundary layer control introduced by the propeller slipstream. Riley [58] also provided a detailed description of business/commuter turboprop aircraft configuration with aft engines. Although the configuration is a twin-pusher, general information about handling qualities during the approach can be extracted. Generally, the aircraft can have acceptable handling qualities and can be certified, but more care to achieve proper dynamic stability and desired handling characteristics may be required.

SE 2 A SR Design Refinement Using MDO
To refine the initial design, a set of multi-disciplinary design optimizations (MDO) were performed. The main goal of the refinement is to reduce the aircraft DOC further. Multiple cases were considered. The first MDO featured all prescribed assumptions. The second MDO assumed that no folding penalty is present to estimate the influence of the folding penalty of the aircraft. This case is used only as a comparative study to investigate the folding mechanism effect on the optimal design. Finally, a special case of battery weight minimization was calculated to investigate the reduction of the battery weight and increase of DOC with respect to the aircraft optimized for the DOC. This case also directly influences the aircraft CO 2 emission.
Production of electricity also involves emission CO 2 and other greenhouse gases, since modern electricity generation is not a product of renewable energy sources. A major portion of electricity today comes from fossil fuels, consequently, an all-electric aircraft cannot be considered as a zero-emission one. Since the goal of the study was to design and assess the aircraft from the environmental perspective as well as performance and costs, proper metrics of preliminary emission rates must be introduced. Scholz describes a rapid and useful procedure of comparing kerosene-based and battery-electric emission [59,60] to calculate the CO 2 emission including additional factors for each energy source. For aircraft operating with kerosine-based fuels, additional chemicals other than CO 2 contribute to the emission level. In addition, the altitude at which the emission happens plays an important role. To account for secondary effects, the equivalent CO 2 emission level per flight can be described by m CO2,eq = (EI CO2 f km + EI NOx f km CF midpoint,NOx + CF midpoint,AIC )R km (12) where f km is the fuel flow per km, EI is the emission index for each species, CF is the characterization factor, and R km is the range in km. The EI CO2 is equal to 3.16 kg CO 2 per kg fuel [60]. The EI NOx is found using the Boeing Method 2 [61] where the index is computed based on existing data from the ICAO Engine Emissions Databank [62] where SGTP is the sustained global temperature potential, summarized in Table 4 and s is the forcing factor which depends on the altitude and is shown in Figure 12 for each species. The formulation given in Equation (12) was used for each flight segment of the reference ATR-72 and then added to obtain the total mission equivalent emission level. Segments with variable altitude such as climb and descent used average values of forcing factors s and emission index EI NOx . If the value of the forcing factor was for the altitude lower than the available data, the minimum value represented in the Figure 12 was used.  Battery emission, on the other hand, depends on the amount of energy used to charge them. The emission then depends on the share of fossil fuel to non-fossil fuel sources participating in the electricity production and the primary efficiency factor-the factor which describes the ratio between the energy required to produce electricity and the final output energy amount. Including all factors described above, the equivalent emission per flight becomes where x f f is the fossil fuel share ratio equal to 0.39 according to the forecast of 2050 [63], k PEF is the primary efficiency factor assumed equal to 2 [59], η charge is the battery charging efficiency equal to 0.9 [59], and C bat and C Jet−A are battery and Jet-A fuel energy densities, respectively. From Equation (15), the CO 2 emission is directly proportional to the battery mass, so the optimization for the minimum battery weight corresponds to the minimum emission case. All problems have similar design variables and constraints, as shown in Table 5. Due to the application of new airframe technologies, requirements for the sea-level climb, take-off, and landing are more significant compared to cruise and OEI requirements, as shown in Figure 10. Such behavior occurs due to a substantial reduction of the cruise drag, which strongly diminishes the required cruise power. Consequently, optimization problem constraints can be limited to these three main performance metrics. Additional geometric constraints of 51.8 m wingspan (similar to the wingspan of the NASA N+3 strut-braced wing concept [64]) and the wingtip of no larger than 1.4 m were introduced to limit further increase in the wing slenderness. Table 5 describes the formulation of the optimization problem. There, C r and C t are the wing root and tip chords, respectively, t/c is the wing thickness, and η max is the maximum throttle during the flight. To perform the MDO, SUAVE was coupled with MATLAB, and the Genetic algorithm was used to find optimal solutions. To account for the constraints, the penalty function similar to the one defined in [65] was used. The penalty function is defined by where y is the design variable, y c is the design variable constraint, µ is the unit step function equal to zero for y ≤ y c , and γ = 3. With the introduction of the penalty function, the objective function becomes where N is the total number of design variables. Sixty species per generation were set to have sufficient population size without major accuracy losses. Take-off and landing constraints were calculated using physics-based time-dependent performance formulations presented in [41]. The climb power-to-weight ratio constraint was calculated using a similar formulation used for the constraint analysis. Figure 13 shows the aircraft planforms for each optimization problem and Table 6 summarizes important properties of each optimized configuration. Limiting constraints for each optimization were similar to the constraint diagram: the solution was trying to minimize the power-to-weight ratio which reached the constraint limit of either the climb or the take-off. In addition, the top-of-climb throttle constraint was an additional constraint which dominated the others.  As shown in Table 6, the optimized aircraft with all imposed assumptions reduced all weights, slightly reduced the aspect ratio, increased the wingspan, reduced the taper ratio, and reduced the aircraft power-to-weight ratio. Minor differences in the field performance occurred. The DOC for the optimized aircraft reduced by 124 EUR/flight. In addition, the difference in emission between the initial and DOC-optimized designs is equal to 84 kg per flight which is not substantial. If the folding penalty is removed, then the wing has more freedom to increase its span and aspect ratio due to the absence of the folding weight penalty. The optimized aspect ratio increased by 1.23, the wingspan increased by 2.6 m compared to the initial design. However, the power-to-weight ratio increased due to a reduced wing planform area. All weights were further reduced compared to the aircraft configuration with the folding mechanism. The DOC was reduced by 206 EUR/flight compared to the initial configuration, which indicates a relatively strong effect of the folding mechanism on the aircraft. The emission level reduced by 201 kg compared to the initial design which is more substantial. Comparing the DOC optimal solutions with and without the folding mechanism, it becomes evident that wing folding affects the final design and cannot be neglected. If the battery minimization strategy is considered, then the wing planform approaches design limits. Due to a stronger reduction of the battery weight compared to the increase in the wing weight, minor reductions of the MTOM and empty weights are observed. Finally, DOC increased by 235 EUR/flight compared to the initial design or by 359 EUR/flight compared to the optimized configuration. The battery weight reduced by 2404 kg and 1883 kg compared to the initial and DOC-optimized aircraft, respectively. The emission level, however, reduced by 385 kg compared to the initial configuration. In addition, the wing span of the aircraft increases to 51.7 m and the folding portion from each side of the wing becomes 7.85 m. Such a large folding portion may affect the operational convenience at the airport and will make the aircraft taller when the wings are folded. That may create additional difficulties related to the aircraft storage.
If DOC and emission levels are compared against the reference aircraft, the increase in DOC demonstrates a substantial increase compared to the reference ATR-72. The most cost-effective aircraft version has 41% increase in DOC which is a substantial cost growth from the market standpoint. The DOC for the battery-optimal configuration increases by 46%. On the other hand, if emission levels are compared, the most DOC-efficient version has a reduction in equivalent CO 2 emission of 83% compared to the ATR-72, and the battery-optimal configuration has 86% emission reduction. From the comparison, the emission level of the all-electric aircraft reduces substantially, but such reduction of the emission comes at a significant cost increase. Moreover, DOC shows more sensitivity with respect to the configuration compared to the emission level. To minimize costs as much as possible, the DOC-optimal configuration was selected for future analyses.
Observing the wingspan deviations from the Part C airport requirements of 36 m, aircraft configurations minimized for DOC did not show major increases in wingspans. Consequently, it is possible to design an aircraft that will have no folding mechanism and will not have a significant DOC penalty so that the wing design complexity may be reduced substantially. An additional MDO study with a fixed wingspan constraint was performed to investigate the sensitivity of DOC to the configuration with a fixed span of 36 m. Table 7 demonstrates a comparison between the initial design, the optimal design without the span constraint, the optimal aircraft with a constrained span, and a reference ATR-72 to compare optimal models to their reference. The aircraft with a fixed span has a higher weight due to a restricted planform and, as a consequence, more battery is required for the mission. The battery weight increased by 1600 kg compared to the optimal design without the span constraint. The power-to-weight ratio also increased to satisfy the take-off requirement. However, the wing area is reduced, so it is less expensive to increase the motor power than increase the wing area. The DOC increased by 94 EUR/flight, which is not a large increase. At the same time, the wing has no folding mechanism and is easier to design and maintain. However, the emission level increases by a relatively small margin compared to the initial DOC-optimal configuration. At this moment, the configuration without folding is selected as the final decision.
The geometric summary of the selected aircraft is shown in Table 8. Single-slotted Fowler flaps were used for this aircraft. The empennage was initially sized within SUAVE using a constant tail volume fraction similar to existing aircraft in the class using data from Raymer [40] and then refined using AVL to ensure sufficient stability and trim during critical flight cases for the complete CG envelope obtained separately. The vertical tail was sized based both on the one-engine inoperative (OEI) case and the lateral stability condition of C N β > 0.01573 [66]. The propeller was positioned vertically as far as possible to reduce the possibility of the wing wake impinging on the propeller blades. The angle between the wing trailing edge and the propeller's lowest blade tip is equal to 14 deg, which should be sufficient for the majority of operational cases. The propeller features six blades and has a diameter of 3.56 m.
The CG envelope was created using available information about the stick-fixed and stick-free neutral points, payload clouds, and multiple baggage arrangement cases. The CG-range is equal to 21.5% mean aerodynamic chord (MAC) starting at 12% and ending at 33.5% MAC. Figure 14 shows the SE 2 A SR payload-range diagram compared to the reference ATR-72 [67]. In addition, the obtained payload-range diagram was validated using a Breguet range equation for an all-electric aircraft defined by where m 0 is the maximum take-off mass, g is the gravitational acceleration, and L/D is the cruise lift-to-drag ratio. For the range analysis, multiple options were considered. The first option is the case when the battery mass remains constant, independent of changes in payload weight. This way, the battery volume is fixed and equal to the one used for the maximum passenger case. The second option is to have sufficient battery volume to replace missing passengers with an extra battery. Such an option is more problematic to achieve due to the increased battery size and additional growth of the already large belly fairing but is useful to analyze for the range sensitivity analysis. The comparison shows that although the SE 2 A SR satisfies the harmonic range requirement, which matches the ATR-72, the ferry range is significantly limited due to the battery energy density and weight limitations. The absence of the payload extends the aircraft range to 1415 km, which is substantially smaller than the ATR-72 with its ferry range of 3380 km. Introduction of more battery which replaces the payload increases the ferry range until 1870 km, which is 455 km longer than the fixed battery weight case. However, this ferry range still does not introduce major range improvements compared to the ATR-72. Comparison between the SUAVE payload-range and Breguet payload-range diagrams show a difference of both harmonic and ferry ranges of 100 km. Such tendency happens due to the higher fidelity of SUAVE which includes all mission segments and calculates incremental energy consumption and aerodynamic characteristics, unlike the Breguet equation which considers a constant lift-to-drag ratio. On the other hand, minor deviations are expected and are satisfactory for the selected analysis.
To summarize current design outcomes, an all-electric SE 2 A SR aircraft is capable of achieving the harmonic range similar to the reference ATR-72 and satisfies all prescribed TLRs. Moreover, its emission level is reduced by 81% compared to the reference aircraft. On the other hand, the ferry range of the all-electric aircraft is limited to no more than 1870 km compared to 3380 km for the ATR-72 and has a 42% increase in DOC with respect to the reference which will create substantial market application problems for this type of airplane. Figure 15 shows a rendered image of the optimized SE 2 A SR aircraft.

Sensitivity Analysis of the Aircraft Characteristics to Technology Advancements
After performing a conceptual design of a regional all-electric aircraft and analyzing obtained results, several questions still remain: What is the influence of novel technologies on the aircraft configuration? How much do the deviations of novel technology assumptions affect the airplane configuration? What technologies affect the aircraft geometry the most?
A sensitivity analysis of the SE 2 A SR aircraft to the battery energy density with the absence of all novel airframe technologies was performed to determine what potential energy density could make all-electric regional aviation possible and how the aircraft characteristics are affected. The aircraft without novel airframe technologies assumed 5% laminar flow on the wing and empennage, 5% laminar flow for the fuselage, a limit load factor of 2.5, and no weight reduction due to composite materials. Technologies related to the propulsion system such as HTS motors remained since their absence immediately leads to the aircraft's infeasibility due to the Megawatt-level power required per motor. Furthermore, the performance characteristics of the aircraft without novel technologies were compared against the case if all technologies met the desired requirements. The sizing procedure for all configurations was performed for a similar mission profile and used similar methods to those described in Section 3 to ensure the satisfaction of all TLRs. Similar constraints based on TLRs were also implemented. As for the initial sizing, wing loading and aspect ratio sensitivity analyses were performed to determine a sweep of all possible aircraft that satisfy TLRs. For all possible configurations, the ones with minimum DOC were selected for each selected battery energy density. Finally, all characteristics were compared against the reference ATR-72 aircraft. Figure 16 shows the effect of battery energy density on the aircraft gross weight.
Results demonstrate that the absence of novel technologies significantly limits the feasibility of all-electric regional aviation unless substantially higher battery cell energy densities are reached. For instance, if the aircraft with all available airframe technologies can reach the MTOM of 40,000 kg at the pack energy density of 600 Wh/kg, the aircraft without airframe technologies can reach similar weights at 900 Wh/kg. In addition, weights approach an asymptotic value at high energy densities, so the effect of airframe technologies becomes more significant for high energy densities. For instance, if the pack energy density reaches 1000 Wh/kg, then the aircraft MTOW with all technologies becomes slightly less than 30,000 kg while similar weight is achieved by the aircraft without airframe technologies only at 1500 Wh/kg. To compare the DOC values, the aircraft without technologies did not include an additional maintenance cost penalty introduced for the aircraft with all airframe technologies. Moreover, a sensitivity analysis of the maintenance cost gain factor was performed to investigate its effect on the overall DOC and predict potential future DOC behavior. If DOC is compared, two trends can be observed. First, both aircraft have substantially higher DOC compared to the reference ATR-72, and the difference remains large even for high energy densities. On the other hand, the application of technologies affects the DOC. Comparing aircraft with and without technologies, there exists a point after which the DOC increases due to the airframe technologies and their maintenance complications and is equal to 980 Wh/kg. After this energy density, an aircraft without airframe technologies becomes cheaper to operate. Such a trend is a consequence of the maintenance cost gain and redistribution of dominant cost drivers, as shown in Table 9. For low battery energy densities, the most dominant DOC component is the capital cost which depends on the airframe and propulsion system weights. Due to significantly larger weights, these costs contribute to total DOC such that benefits of maintenance costs disappear. Moreover, more battery energy is required for the aircraft without technologies and slightly higher fees also contribute to a larger total DOC. On the other hand, with the increase in energy density, aircraft weights reduce substantially, so capital costs reduce as well. The overall distribution of costs becomes more evenly distributed for the aircraft without technologies, so lower maintenance cost effects become more significant. On the other hand, maintenance costs for the aircraft with technologies start dominating. That relationship leads to higher overall DOC for the aircraft with technologies at high battery energy densities. The present aircraft was designed having 700 Wh/kg of pack energy density which is still more beneficial than the aircraft without technologies. However, the maintenance penalty factor of 2 may be too optimistic. If the maintenance costs increase more, then the equilibrium point between the two DOCs will shift towards lower energy densities, as shown in Figure 16. The only possibility of approaching the reference aircraft DOC exists if maintenance cost gain is equal to 1. The maintenance cost of novel airframe technologies shall be considered later to ensure the financial applicability of novel technologies further. The effect of technology deviations on the aircraft weights and, as a consequence, applicability of all-electric aircraft is shown in Figures 17-19. In addition, points for extreme technology cases are summarized in Tables 10-13. A similar sizing approach, as discussed in Section 3, was used for the sizing with technology deviations. The wing loading, however, was fixed to the one that occurred after the MDO of the aircraft. In addition, to simplify the sensitivity and account for significant complexities to design an appropriate suction system for the fuselage, the fuselage was assumed fully turbulent for all technology sensitivities.      Based on obtained results, if the weight sensitivity on the technology is compared for two other technologies where one technology has a constant value and the other technology is fixed at either their minimum or maximum, the laminar flow control shows the largest influence on the MTOM. Structural weight reduction plays the second most important role, and the load alleviation affects the aircraft MTOM the least compared to the other two technologies. For instance, for the airframe reduction of 0%, if the maximum possible boundary layer transition is possible, then the change in MTOW due to the load alleviation is equal to 2435 kg. However, if for the same airframe reduction factor, the load alleviation has its minimum value, then the change in MTOM is equal to 2835 kg due to the HLFC technology. If technologies are combined, then the effect of the MTOM is magnified.

Conclusions
The present manuscript investigated the influence of novel technologies on the feasibility of all-electric short-range aircraft with passenger capacity and performance characteristics similar to ATR-72. Three novel airframe technologies were considered: hybrid laminar flow control, active load alleviation, and advanced materials and structure concepts. In addition, superconducting electric motors and high energy density batteries were considered for the aircraft energy system.
Multiple tools were used and improved to assess the feasibility question. SUAVE was used to perform the initial sizing and sensitivity analyses and was also coupled with MATLAB to perform MDO studies. EMWET and FLOPS weight estimation methods and a weight penalty function for folding wings were integrated into SUAVE. The cubic spline method for a propeller was used to enable the aircraft sizing. AVL was used to assess the aircraft stability and control.
Two aircraft configurations were considered for the design: the conventional configuration with a high wing and wing-mounted propellers and the low wing with propellers mounted at the aft of the fuselage. The low-wing configuration was selected because it had lower DOC, lower weights, and is possible to have certifiable stability and control characteristics.
To refine the selected concept and determine the aircraft characteristics' sensitivity to the folding penalty and the objective function, a set of MDO studies were performed. The configuration without folding wings was selected due to a relatively small increase in DOC due to the fixed wingspan and reduced complexity of the wing design and maintenance due to the absence of the wing folding mechanism.
The importance of airframe technologies was assessed by comparing aircraft configurations with and without them. The presence of airframe technologies may reduce the required battery pack energy density by 300 Wh/kg compared to their absence. The difference in required energy density also increases with the increase of the battery energy density. However, the application of technologies becomes more expensive than their absence if the battery pack energy density exceeds 980 Wh/kg for a given maintenance cost gain factor of 2.
Comparison of the all-electric aircraft to the reference ATR-72 showed substantial reduction of CO 2 emission by 81% of the all-electric aircraft. On the other hand, DOC increase of 42% presents a significant limit to market success of such aircraft. Moreover, neither increase in battery energy density, nor reduction of maintenance costs change aircraft DOC significantly. The only possibility to have a comparable DOC level may appear if maintenance cost gain is absent and battery pack energy density exceeds 1000 Wh/kg. Moreover, the all-electric aircraft consumes a lot of energy to complete the mission which creates major uncertainties related to its energy-efficiency from the operational standpoint.
Finally, a sensitivity analysis of technology availability on the aircraft characteristics was performed. Simulations were performed for each technology combination from the most optimistic to its absence. Results demonstrated high sensitivity of aircraft weights to technology levels where the HLFC has the most influence, followed by the aircraft weight reduction and finishing with the load alleviation. The combined effect of technologies magnifies the differences in weight even further.
Future research steps will include a detailed conceptual design of the SE 2 A SR aircraft. First of all, the design with available models for the laminar flow required power and additional weight penalties due to the implementation of all technologies and additional hardware will be performed. A more sophisticated assessment of the aircraft stability and control must be performed to ensure adequate flying qualities. Finally, higher fidelity analysis for the energy network and modeling of the motors and propellers will be considered.