Next Article in Journal
Observability Analysis and Improvement Approach for Cooperative Optical Orbit Determination
Previous Article in Journal
A Machine Learning Approach for Global Steering Control Moment Gyroscope Clusters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Design of a Distributed Ducted Fan System

School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(3), 165; https://doi.org/10.3390/aerospace9030165
Submission received: 28 February 2022 / Revised: 13 March 2022 / Accepted: 13 March 2022 / Published: 17 March 2022
(This article belongs to the Section Aeronautics)

Abstract

:
The distributed propulsion system applied to the vertical take-off and landing aircraft must maintain the high performance in both hover and cruise flight. The gap of the power unit has an adverse effect on the distributed ducted fan system, especially in cruise flight. Therefore, a new distributed ducted fan system was proposed, which eliminated the power gap in design, and adjusted the contraction and expansion of the wake through the deflectable induced wing arranged behind the ducted fan to ensure the high efficiency of the distributed ducted fan system in different flight phases. Then, a multi-objective design method of the distributed ducted fan system was proposed, and the feasibility of the design method was verified by designing the inlet and outlet of the duct and the induced wing. Design results show that the performance change of the distributed ducted fan system mainly came from the change of the inlet. By increasing the length and height of the inlet, the flow separation was alleviated and the duct thrust was increased in hove flight, but the cruise drag became larger. The increase of the inlet height made the operating point of the blade far away in hover and cruise flight, which increased the difficulty of the multi-objective design. Compared with the distributed ducted fan system composed of the traditional circular ducted fan, the hovering power load was reduced by 3.703%, but the cruise efficiency was increased by 17.372%, and the spanwise space was reduced by 20% in the final design.

1. Introduction

In recent years, the electric vertical take-off and landing (VTOL) aircraft has become a research hotspot in aviation due to its characteristics of low emission, low noise, and convenient take-off and landing, such as Lilium jet of German Lilium Company and XV-24A Lightning Strike aircraft of American Aurora Company (Torrance, CA, USA). Both aircraft adopt the distributed propulsion system.
The distributed propulsion system can use boundary layer ingestion (BLI) [1] and wake-filling [2] to improve the propulsion efficiency, increase lift by increasing dynamic pressure [3], reduce requirements for traditional control surfaces and increase the redundancy of the control system, and reduce noise signatures [4,5], making it a competitive solution in future aviation [6]. Studies have shown that N3-X aircraft with turboelectric distributed propulsion (TeDP) [7] to be capable of achieving a 60% reduction in total energy consumption over current aircraft [8]. Wick studied three integration distributed propulsion concepts: Lower TE (Trailing Edge), Upper TE, and Imbedded configurations by using the computational fluid dynamics (CFD) method based on a modified boundary condition, and the results indicated that the Lower TE configuration was the best performer [9]. Kerho studied the influence of fan thrust angle and fan thrust level of the overwing distributed propulsion system using an actuator disk model [10]. Zhang used the momentum source method (MSM) to study the lift enhancement effect of the distributed ducted fan system on a hybrid wing body UAV [11]. However, the above studies simplify the propulsion system, and do not consider the influence of the airframe on the distributed propulsion system.
Both Lilium jet and XV-24A Lightning Strike aircraft used electric ducted fans as the power of the distributed propulsion system, because the ducted fan can generate thrust and reduce the tip losses of blades [12]. The ducted fan in low Mach number has higher aerodynamic efficiency than the propeller with the same disc load. Stubblefield used the lift line model to design the blade of the ducted fan [13]. Guo established an efficient blade design method for the ducted fan by modifying the blade element momentum design method through CFD results [14]. Yu used the panel method and CFD method to optimize the blade of the ducted fan [15]. Ye simplified the blade by MSM, and optimized the NASA duct by using response surface model and neural network model, respectively [16]. Ji used the MSM to analyze and optimize the duct shape of the ducted fan [17]. However, there is still a lack of complete research on the design of the duct and fans of the distributed ducted fan system, especially the multi-objective design in hover and cruise flight.
In order to better integrate with the wing, the distributed ducted fan system generally adopted a rectangular duct [10,18,19]. Bento used the actuator disk model and the steady and unsteady CFD methods to study the influence of the duct shape on the performance of the ducted fan [20]. Sun studied the propulsion characteristics of the distributed ducted fan system with different duct shapes [21]. Their research both showed that the corner of the square duct is prone to a flow separation, resulting in a reduction in the efficiency of the ducted fan. Wang optimized the inlet of the rectangular duct in the distributed propulsion system based on the panel method in hover flight [22]. However, since the panel method cannot predict the flow separation, the accuracy of the optimized design needs to be further improved.
At present, most of the research on the distributed propulsion system focuses on configuration characteristic analysis, and the related design research of the ducted fan focuses on the independent propulsion system and single objective, while the multi-objective design research of the distributed ducted fan system is less. Therefore, this paper first proposed a new distributed ducted fan system for the VTOL aircraft. By establishing the parameterization method of the inlet and outlet of the rectangular duct, and using the efficient MSM to obtain the performance, a multi-objective design method of the distributed ducted fan system was proposed based on NSGA-II algorithm and Kriging model. Finally, the design method was verified by the design comparison.

2. Numerical Simulation Methods

Numerical simulation methods of the ducted fan were verified by using a 150 ducted fan shown in Figure 1 and Figure 2. The duct of the 150 ducted fan has an inlet length of l in = 40   mm , an inlet radius of R in = 90   mm , a length of the blade section of l h = 100   mm , a hub radius of R h = 30   mm , a radius of the inner wall of R = 75   mm , a maximum radius of the outer wall of R max = 96   mm , an outlet length of l out = 100   mm , and an outlet expansion angle of α e = 3 ° . The blade number of the rotor is 10 with a blade tip clearance of 0.5 mm, and the blade number of the stator is 6.
The performance of the 150 ducted fan was tested under no incoming flow state, as shown in Figure 3. The ATI delta sensor was used to measure the aerodynamic force. The “Eagle tree” was used to measure the rotating speed. The power was supplied by a regulated power supply, and the DC power was converted into AC power by an electric regulator. Each test point was tested stably for 5 s, and the results were obtained by time-averaged processing. Since the “HD_wall” shown in Figure 1 in the test state basically did not generate thrust, it was cancelled during the test to facilitate installation.
The multiple reference frame (MRF) method and MSM [23] were used in the numerical simulation of the ducted fan based on Navier–Stokes (N–S) equations and the k-ω SST turbulence model. Moreover, the implicit solver formulation, Roe-FDS flux type, and second-order upwind spatial discretization were used. The MRF method adopted the structured grid with 1.04 million background cells and 1.92 million blade cells, as shown in Figure 4. The grid with 1.13 million cells used in the MSM is shown in Figure 5, and the momentum source term was added in the middle part of the duct to be equivalent to the effect of the blade.
In the MSM, the momentum source term S m was calculated through Equation (1), where blade thrust T p is calculated by the corresponding MRF method, and v o l is the total volume of the grid to which the momentum source term was added.
S m = T p / v o l
The calculation was carried out on a computer using Ryzen 3900X processor of AMD. The result comparison between the two CFD methods and experiment is shown in Figure 6 and Table 1. The errors of the two CFD methods are within 10% near the design point. However, the time consumption of the MSM is only 1/9 of the MRF method. Therefore, the MSM can be used to design the ducted fan, and the MRF method can be used for more accurate characteristic analysis.

3. Characteristic Analysis of the Distributed Ducted Fan System

3.1. Distributed Ducted Fan System

Traditional circular ducted fans were arranged to form the distributed ducted fan system shown in Figure 7 with a power gap of d = 5   mm . The duct of the distributed ducted fan unit is shown in Figure 8. The duct has an inlet length of l in = 50   mm , an inlet radius of R in = 85   mm , an outlet length of l out = 100   mm , a length of the blade section of l h = 100   mm , a radius of the inner wall of R = 75   mm , a hub radius of R h = 30   mm , and a maximum radius of the outer wall of R max = 97.5   mm .
In order to take into account the hover and cruise performance of the distributed ducted fan system, the variable area nozzle (VAN) [24] was adopted as shown in Figure 8. The expanded outlet was used in hover flight to increase the flow velocity of the blade section, and the contracted outlet was adopted in cruise flight to decrease the flow velocity of the blade section, so that the operating point of the rotor and stator in different flight phases was close.
The rotor and stator of the ducted fan were designed in hover flight, using airfoils as shown in Figure 9. The airfoil camber is constant, and the airfoil thickness decreases linearly from the root to the tip. The geometric information of the rotor and stator is shown in Figure 10. The blade number of the rotor is ten, the blade number of the stator is six, and blade tip clearance is 0.5 mm.
It can be seen from Figure 7 that the distributed ducted fan has a high degree of repeatability along the span. In order to improve the efficiency of the calculation and design, a ducted fan unit is selected and its two sides are set as symmetrical boundaries, as shown in Figure 11. The MRF method was used to analyze the characteristics of the distributed ducted fan system.

3.2. Performance Analysis

The thrust demand of the distributed ducted fan unit in hover flight is T = 40   N , and the hover performance of the distributed ducted fan unit is shown in Table 2, where p l is the power load calculated according to Equation (2). The outlet expansion angle of the ducted fan is α e = 3 ° .
p l = T Ω Q r
In cruise flight of V 0 = 50   m / s , the thrust demand of the distributed ducted fan unit is T = 4   N , and the cruise performance of the distributed ducted fan unit can be improved by adjusting the outlet expansion angle. The efficiency η of the ducted fan is calculated through Equation (3).
η = T V 0 Ω Q r
It can be seen from Figure 12 that at Ω = 8500 r/min, the thrust and efficiency of the distributed ducted fan unit are improved with the decrease of outlet expansion angle, and the best outlet expansion angle is α e = 5 ° . The efficiency of the distributed ducted fan unit under cruise thrust demand is only 35.136%, as shown in Table 3.
Comparing the pressure distribution of the distributed ducted fan unit in two flight phases, as shown in Figure 13 and Figure 14, it can be seen that the mutual influence between the ducted fans reduces the low pressure near the inlet. The low pressure changes little in hover flight, and the flow between the adjacent ducts accelerates in cruise flight, causing the low pressure area on the outer wall of the duct to move backward. In addition, there is a flow separation at the tail of the duct, as shown in Figure 15, which seriously affects the cruise performance of the distributed ducted fan unit.

4. New Distributed Ducted Fan System and Related Parameterization Methods

4.1. New Distributed Ducted Fan System

Since the power gap is not conducive to the cruise performance of the distributed ducted fan system, the duct is considered to be rectangular, as shown in Figure 16, to eliminate the power gap.
In addition, in order to achieve simpler outlet control, the distributed ducted fan system of “ducted fan + induced wing” shown in Figure 17 was adopted. The deflectable induced wing is arranged behind the ducted fan. The symmetrical deflection of the induced wing can realize the expansion and contraction of the wake, and the same direction deflection of the induced wing can be used to generate vertical aerodynamic force.
Based on the above analysis, the inlet of the distributed ducted fan unit was designed in the form of “rectangle to circle” and the outlet was designed in the form of “circle to rectangle”, as shown in Figure 18. In this way, the ducted fan units of the distributed ducted fan can be seamlessly connected from left to right, and the outer wall of the duct can be shaped like the upper surface of the wing, as shown in Figure 19.

4.2. Parameterization Method of the Inlet and Outlet

4.2.1. Hyperelliptic Equation

The hyperelliptic equation can realize the transition from circle to rectangle, which was widely used in the generation and parameterization of the inlet [22,25]. The hyperelliptic equation is as follows, where a is the major semi-axis, b is the minor semi-axis, and m is the exponent of the hyperelliptic equation.
( x a ) m + ( y b ) m = 1
The hyperelliptic equation is an elliptic equation when m = 2 , and when m > 2 and m increases, the hyperellipse is gradually transformed into a rectangle.

4.2.2. Parameterization Method of the Inlet

The inlet of the distributed ducted fan unit adopted the form of “up-down and left-right symmetry”, so only the 1/4 shape of the inlet was parameterized, as shown in Figure 20. The inlet consists of an external inlet and an internal inlet. The internal inlet is a “rectangle to circle” transition surface, and the external inlet is used to connect the internal inlet and the outer wall of the duct. It also can be seen from Figure 20 that the inlet is controlled by three curves: a ( t ) , b ( t ) , and m ( t ) .
The following fourth-order Bezier curves were used to parameterize the three control curves, where ( p 0 , p 1 , , p 19 ) are the design variables. The design variables were unitized according to the geometric information, where h ii is the height of the internal inlet, h in is the height of the external inlet, and w in is the inlet width.
x ( t ) = l in [ p 0 ( 1 t ) 4 + 4 p 1 t ( 1 t ) 3 + 6 p 2 t 2 ( 1 t ) 2 + 4 p 3 t 3 ( 1 t ) + p 4 t 4 ]
y a ( t ) = R + ( w in 2 R ) [ p 5 ( 1 t ) 4 + 4 p 6 t ( 1 t ) 3 + 6 p 7 t 2 ( 1 t ) 2 + 4 p 8 t 3 ( 1 t ) + p 9 t 4 ]
y b ( t ) = R + R [ p 10 ( 1 t ) 4 + 4 p 11 t ( 1 t ) 3 + 6 p 12 t 2 ( 1 t ) 2 + 4 p 13 t 3 ( 1 t ) + p 14 t 4 ]
y m ( t ) = R + ( ( w in 2 ) 2 + ( h ii 2 ) 2 R )       [ p 15 ( 1 t ) 4 + 4 p 16 t ( 1 t ) 3 + 6 p 17 t 2 ( 1 t ) 2 + 4 p 18 t 3 ( 1 t ) + p 19 t 4 ]
Since the starting point of the Bezier curve is the first control point and the end point is the last control point, the following variables were obtained.
p 0 = p 5 = p 10 = p 15 = 0
p 4 = p 9 = p 19 = 1
p 14 = ( h in / 2 R ) / R
In order to ensure the horizontal connection between the inlet and the blade section and the vertical connection between the inlet and the outer wall of the duct, the following variables were obtained according to the derivative properties of the Bezier curve.
p 6 = p 11 = p 16 = 0
p 3 = 1
Therefore, there were eight design variables ( p 1 ,   p 2 ,   p 7 ,   p 8 ,   p 12 ,   p 13 ,   p 17 ,   p 18 ) for the duct inlet.
After determining the three control curves, the information of the major and minor semi-axis of the hyperellipse was obtained directly from the corresponding Bezier curve. The exponent of the hyperelliptic equation was obtained by solving the Equation (14).
y m ( t ) = ( y b ( t ) 0.5 m ) 2 + ( y a ( t ) 0.5 m ) 2
The length of the internal inlet is l ii , and t ii can be obtained from Equation (15).
l ii = x ( t ii )
When the internal inlet was generated, the variable range of x ( t ) and y b ( t ) was t [ 0 ,   t ii ] , and the variable range of y a ( t ) and y m ( t ) was t [ 0 ,   1 ] . When the external inlet was generated, the variable range of x ( t ) and y b ( t ) was t [ t ii ,   1 ] .

4.2.3. Parameterization Method of the Outlet

The outlet of the distributed ducted fan unit also adopted the form of “up-down and left-right symmetry”, but did not divide the internal and external parts. The outlet design ensured that the end of the upper and lower outlet was horizontal.
Similarly, the fourth-order Bezier curves such as Equations (5)–(8) were used to parameterize the outlet, except that y b ( t ) was unitized by Equation (16).
y b ( t ) = R + ( h out 2 R ) [ p 10 ( 1 t ) 4 + 4 p 11 t ( 1 t ) 3 + 6 p 12 t 2 ( 1 t ) 2 + 4 p 13 t 3 ( 1 t ) + p 14 t 4 ]
The following variables were obtained according to the characteristics of the outlet and Bezier curve.
p 0 = p 5 = p 6 = p 10 = p 11 = p 15 = p 16 = 0
p 4 = p 9 = p 13 = p 14 = p 19 = 1
So, the eight design variables of the duct outlet were ( p 1 ,   p 2 ,   p 3 ,   p 7 ,   p 8 ,   p 12 ,   p 17 ,   p 18 ) .
When the outlet was generated, the variable range of x ( t ) , y a ( t ) , y b ( t ) , and y m ( t ) was t [ 0 , 1 ] .

5. Multi-Objective Design of the Distributed Ducted Fan System

Multi-objective design was carried out on the ducted fan unit and the corresponding induced wing in the distributed ducted fan system. The specific design steps are shown in Figure 21. Firstly, the multi-objective optimization design of the inlet and outlet of the duct was carried out to obtain the duct geometry. Then, the blade of the distributed duct fan unit was designed in hover flight. Afterwards, the chord length and the deflection angle of the induced wing in cruise flight were designed according to the optimal cruise performance. Finally, the deflection angle of the induced wing was adjusted in the hove flight.

5.1. Multi-Objective Optimization Design of the Inlet and Outlet

5.1.1. Optimization Objectives and Methods

Based on NSGA-II algorithm and Kriging model, the inlet and outlet of the distributed ducted fan unit were optimized. The outer wall of the duct was generated according to the stretching of the NACA0010 airfoil.
Referring to the ducted fan in Section 3.1, the length of the inlet and outlet remained unchanged, the inlet with the width of w in = 160   mm and the height of h in = 140   mm to keep the inlet area basically consistent, and the outlet with the width of w out = 160   mm and the height of h out = 130   mm in the basic distributed ducted fan unit.
Then, the effect of increasing the length and height of the inlet was studied. The inlet length was changed to l in = 70   mm and l in = 100   mm , and the inlet height was changed to h in = 150   mm and h in = 160   mm . Other geometric parameters remained unchanged.
Since the inlet and outlet adopted the symmetrical form, the optimization was carried out on the 1/4 distributed ducted fan unit. The calculation model is shown in Figure 22, and the left, right, and bottom are set as symmetrical boundary conditions. The grid has 0.44 million cells.
The MSM was used to get the performance of the distributed ducted fan unit. The blade thrust in hover flight is T p = 18   N and in cruise flight is T p = 7   N . Optimization objectives are the maximum duct thrust T d in hover flight and minimum duct drag D d in cruise flight.
The initial sample of the Kriging model is 100, and the optimization ends when the sample number exceeds 400. In NSGA-II optimization, the population size is 50, and the maximum evolution number is 1000. The optimization design was carried out in three rounds. Starting from the second round, the Pareto set obtained in the previous round was added to the initial sample in this round to improve the accuracy of the optimization design.

5.1.2. Optimization Results and Analysis

The comparison of different optimization design results is shown in Figure 23 and Figure 24. The results show that the increase of the duct thrust in hover flight contradicts the decrease of the duct drag in cruise flight. The selected design results shown in Figure 23 and Figure 24 were analyzed. “50_80_70” is the result of the basic distributed ducted fan unit, and its three numbers represent the length, width, and height of the 1/4 inlet. The duct thrust is shown in Table 4 and Table 5.
It can be seen that the difference of the duct in different optimization results is mainly concentrated at the inlet, and the change of the outlet is small, as shown in Figure 25 and Figure 26. The duct thrust of the basic distributed ducted fan unit in hover flight is poor, and increasing the length or height of the inlet is beneficial to improve the duct thrust, but it also increases the cruise drag.
The flow field information is shown in Figure 27, Figure 28, Figure 29, Figure 30, Figure 31. In hover flight, there is a strong flow separation at the inlet of the basic ducted fan unit, and the increase in the length and height of the inlet can alleviate the flow separation and increase the duct thrust. However, the increase in the length and height of the inlet expands the windward area and wetted area of the duct, which decreases the duct thrust in cruise flight.
The difference of the duct thrust in hover flight between “100_80_70” and “50_80_80” is only 2.368%, but the cruise drag of “50_80_80” increases by 28.195%. In contrast, the increase in inlet height increases the inlet area and has a more significant effect on cruise. In addition, the increase of the inlet area accelerates the flow velocity and reduces the pressure in the blade section in cruise flight, as shown in Figure 31 and Figure 32, which causes the operating point of the rotor and stator in different flight phases far away, adding the difficulty of the multi-objective design.

5.2. Blade Design

After obtaining the duct geometry of the distributed ducted fan system, the rotor and stator of the ducted fan were designed in hover flight. In order to compare with the distributed ducted fan system shown in Section 3.1, the same rotor and stator were used without redesign. The ducted fan unit in the distributed ducted fan system is shown in Figure 32.
The hover performance of the ducted fan unit is shown in Table 6. Since the selected design takes into account the cruise performance and is not the best result in hover flight, its power load is reduced by 5.285% compared to the ducted fan unit in Section 3.1.

5.3. Design of the Induced Wing

The main function of the induced wing is to contract the wake of the distributed ducted fan in cruise flight to reduce the flow velocity of the blade section and the inflow angle of the blade element, so that the blade designed in hover flight still has a higher efficiency in cruise flight.

5.3.1. Induced Wing Design in Cruise Flight

The induced wing is 10 mm away from the outlet of the ducted fan and ±67 mm away from the central axis of the duct, as shown in Figure 33, and adopts the NACA0012 airfoil. Induced wings with chord lengths of 120 mm and 150 mm were used respectively, and the deflection angle (expansion is positive) of the induced wing was changed to evaluate the cruise performance of the distributed ducted fan system.
The results shown in Figure 34 show that the total thrust and cruise efficiency of the distributed ducted fan system increase with the decrease of the deflection angle of the induced wing. In contrast, the longer induced wing at the same deflection angle makes the wake more fully contracted, and the flow velocity in the blade section is also lower, as shown in Figure 35. The more fully contracted wake makes the rotor and stator closer to the design state, so the cruise performance is better. When the deflection angle continues to increase, the cruise efficiency of the distributed ducted fan system is reduced due to the flow separation occurring on the induced wing. The induced wing with a chord length of 150 mm and a deflection angle of −12° was selected in the cruise flight.
The cruise performance of the designed distributed ducted fan system is shown in Table 7. Compared with the distributed ducted fan system in Section 3.1, the cruise efficiency is improved by 11.888%.

5.3.2. Induced Wing Design in Hover Flight

The influence of the deflection angle of the induced wing on the hover performance of the distributed ducted fan system was analyzed. As shown in Table 8, the performance corresponding to 1° deflection angle is the best. Compared with the ducted fan unit without the induced wing in Table 6, the hovering power load is increased by 1.955%, which is due to the expansion of the wake, as shown in Figure 36. It shows that the induced wing can be used as a contracted outlet to improve the cruise efficiency of the distributed ducted fan system, and can also be used as an expanded outlet to further improve the hover performance of the distributed ducted fan system.
The above design results show that, compared with the distributed ducted fan system in Section 3.1, the unit width of the distributed ducted fan system designed in this paper is reduced from 200 mm to 160 mm, and the spanwise space is reduced by 20%. Although the hovering power load is reduced by 3.434%, the cruise efficiency is improved by 11.888%, which verifies the feasibility of the multi-objective design method of the distributed ducted fan system in this paper.

6. Improved Distributed Ducted Fan System

6.1. Improvement Method

In order to make the wake continue to contract without the flow separation on the induced wing, the height of the duct outlet can be reduced to make a contracted wake before the deflection of the induced wing. At this time, the induced wing with a smaller chord length can be used to reduce the friction drag and improve the cruise performance.
The height of the duct outlet was reduced to h out = 120   mm , and the duct inlet remained unchanged. The duct profile of the improved distributed ducted fan system is shown in Figure 37. The induced wing with 120 mm chord length was adopted, and the height of the induced wing from the central axis of the duct was changed to ±60 mm.

6.2. Improvement Results

The effect of different deflection angles of the induced wing in cruise flight is shown in Figure 38. The best cruise deflection angle is −12°. The thrust and efficiency of the distributed ducted fan system are better than those before the improvement, which verifies the effectiveness of the improvement method. In hover flight, the power load is the highest when the deflection angle of the induced wing is 2°, as shown in Table 9, indicating that the induced wing makes up for the loss caused by the reduction of the outlet area.
In the improved distributed ducted fan system, the chord length of the induced wing is 120 mm, and the deflection angle is 2° in hover flight and −12° in cruise flight, as shown in Figure 39. The induced wing can further expand the wake in hover flight and contract the wake with no flow separation in cruise flight.
Table 10 and Table 11 show the performance comparison of different distributed ducted fan systems under the thrust demand in hover and cruise flight. The improved design improves the cruise efficiency by 4.902% under the condition that the hovering power load remains basically unchanged. Compared with the distributed ducted fan system in Section 3.1, the hovering power load of the improved design is only reduced by 3.703%, but the cruise efficiency is further increased by 17.372%. The performance of the improved distributed ducted fan system is shown in Figure 40.

7. Conclusions

In this paper, a distributed ducted fan system for the VTOL aircraft was proposed, which eliminates the influence of the power gap in design, and the contraction and expansion of the wake were more simply realized by arranging the deflectable induced wing behind the ducted fan. Then, the parameterization method of the inlet and outlet of the distributed ducted fan unit was established, and the multi-objective design method of the distributed ducted fan system was proposed. Finally, the feasibility of the design method was verified by the design and improvement of the distributed ducted fan system. The conclusions are as follows.
  • The power gap of the distributed ducted fan weakens the low pressure near the inlet and has an adverse effect on the distributed ducted fan system, especially in cruise flight, due to the flow separation at the tail of the duct;
  • In the multi-objective design of the inlet and outlet of the distributed ducted fan system, the outlet obtained by the design is basically the same, and the performance change mainly comes from the change of the inlet. By increasing the length and height of the inlet, the flow separation at the inlet is alleviated, and the duct thrust in hover flight is increased, but it increases the cruise drag at the same time. Compared with the increase of the inlet length, the increase of the inlet height on the one hand brings greater cruise drag, and on the other hand increases the flow velocity in the blade section, which makes the multi-objective design of the distributed ducted fan system more difficult;
  • The performance of distributed ducted fan system can be improved by expanding the wake in hover flight and contracting the wake in cruise flight by the induced wing. The design results show that compared with the distributed ducted fan system composed of the traditional circular ducted fan, the hovering power load of the distributed ducted fan system of the final design is only reduced by 3.703%, but the cruise efficiency is increased by 17.372%, and the spanwise space is reduced by 20%.
Future research will establish a full geometric parameterization model of the ducted fan and induced wing for an integrated design, and further develop the multi-objective optimization design of the overwing distributed ducted fan system.

Author Contributions

Conceptualization, Z.Z. and J.G.; Methodology, J.G.; Software, J.G.; Validation, J.G.; Formal analysis, J.G.; Investigation, J.G. and Z.Z.; Resources, Z.Z.; Data curation, J.G.; Writing—original draft preparation, J.G.; Writing—review and editing, J.G. and Z.Z.; Visualization, J.G.; Supervision, Z.Z.; Project administration, Z.Z.; Funding acquisition, Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Civil Aircraft Specific Project of China (MJ-2015-F-009), Shaanxi Provincial Key Research and Development Program (2021ZDLGY09-08), and Natural Science Basis Research Plan in Shaanxi Province of China (2022JQ-060).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smith, A.; Roberts, H.E. The jet airplane utilizing boundary layer air for propulsion. J. Aeronaut. Sci. 1947, 14, 97–109. [Google Scholar] [CrossRef]
  2. Ko, A.; Schetz, J.A.; Mason, W.H. Assessment of the Potential Advantages of Distributed-Propulsion for Aircraft. ISABE 2003–1094. 2003. Available online: https://www.researchgate.net/publication/238191818_Assessment_of_the_Potential_Advantages_of_Distributed_Propulsion_for_Aircraft (accessed on 1 February 2022).
  3. Stoll, A.M.; Bevirt, J.B.; Moore, M.D.; Fredericks, W.J.; Borer, N.K. Drag Reduction through Distributed Electric Propulsion. In Proceedings of the 14th AIAA Aviation Technology, Integration, and Operations Conference, Atlanta, GA, USA, 16–20 June 2014. [Google Scholar] [CrossRef] [Green Version]
  4. Gohardani, A.S.; Doulgeris, G.; Singh, R. Challenges of future aircraft propulsion: A review of distributed propulsion technology and its potential application for the all electric commercial aircraft. Prog. Aerosp. Sci. 2011, 47, 369–391. [Google Scholar] [CrossRef]
  5. Kim, H.D.; Perry, A.T.; Ansell, P.J. A Review of Distributed Electric Propulsion Concepts for Air Vehicle Technology. In Proceedings of the 2018 AIAA/IEEE Electric Aircraft Technologies Symposium, Cincinnati, OH, USA, 9–11 July 2018. [Google Scholar] [CrossRef] [Green Version]
  6. Lundbladh, A. Distributed Propulsion and Turbofan Scale Effects. ISABE 2005-1122. 2005. Available online: https://publications.lib.chalmers.se/records/fulltext/10074/local_10074.pdf (accessed on 1 February 2022).
  7. Alrashed, M.; Nikolaidis, T.; Pilidis, P.; Jafari, S. Utilisation of turboelectric distribution propulsion in commercial aviation: A review on NASA’s TeDP concept. Chin. J. Aeronaut. 2021, 34, 48–65. [Google Scholar] [CrossRef]
  8. Felder, J.L.; Tong, M.T.; Chu, J. Sensitivity of Mission Energy Consumption to Turboelectric Distributed Propulsion Design Assumptions on the N3-X Hybrid Wing Body Aircraft. In Proceedings of the 48th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Atlanta, GA, USA, 30 July–1 August 2012. [Google Scholar] [CrossRef] [Green Version]
  9. Wick, A.T.; Hooker, J.R.; Hardin, C.J.; Zeune, C.H. Integrated Aerodynamic Benefits of Distributed Propulsion. In Proceedings of the 53rd AIAA Aerospace Sciences Meeting, Kissimmee, FL, USA, 5–9 January 2015. [Google Scholar] [CrossRef]
  10. Kerho, M.F. Aero-Propulsive Coupling of an Embedded, Distributed Propulsion System. In Proceedings of the 33rd AIAA Applied Aerodynamics Conference, Dallas, TX, USA, 22–26 June 2015. [Google Scholar] [CrossRef]
  11. Zhang, Y.; Zhou, Z.; Wang, K.; Fan, Z. Influences of distributed propulsion system parameters on aerodynamic characteristics of a BLI-BWB UAV. J. Northwest. Polytech. Univ. 2021, 39, 17–26. [Google Scholar] [CrossRef]
  12. Black, D.; Rohrbach, C. Shrouded Propellers—A Comprehensive Performance Study. In Proceedings of the 5th Annual Meeting and Technical Display, Philadelphia, PA, USA, 21—24 October 1968. [Google Scholar] [CrossRef]
  13. Stubblefield, J.M. Numerically-Based Ducted Propeller Design Using Vortex Lattice Lifting Line Theory. Master Thesis, Mechanical Engineering Department, Massachusetts Institute of Technology, Cambridge, MA, USA, 2008. [Google Scholar]
  14. Guo, J.; Zhou, Z.; Li, X. An efficient design method for the blade of ducted propeller. Acta Aeronaut. Astronaut. Sin. 2021; in press. Available online: https://hkxb.buaa.edu.cn/CN/10.7527/S1000-6893.2021.25253(accessed on 1 February 2022).
  15. Yu, L.; Drukenbrod, M.; Greve, M.; Abdel-maksoud, M. Research on automatic optimization of ducted propeller design based on CFD techniques. Chin. J. Hydrodyn. 2013, 28, 438–444. [Google Scholar] [CrossRef]
  16. Ye, K.; Ye, Z.; Qu, Z. Aerodynamic optimization method for duct design. J. Aerosp. Power 2013, 28, 1828–1835. [Google Scholar] [CrossRef]
  17. Ji, L.; Li, J.; Fang, Y. The aerodynamic characteristics and parameter optimization design of ducted fan. Sci. Technol. Eng. 2019, 19, 245–251. [Google Scholar]
  18. Perry, A.T.; Ansell, P.J.; Kerho, M.F. Aero-propulsive and propulsor cross-coupling effects on a distributed propulsion system. J. Aircr. 2018, 55, 2414–2426. [Google Scholar] [CrossRef]
  19. Yu, D.; Ansell, P.J.; Hristov, G. Aero-Propulsive Integration Effects of an Overwing Distributed Electric Propulsion System. In Proceedings of the AIAA Scitech 2021 Forum, Virtual Event, 11–15 and 19–21 January 2021. [Google Scholar] [CrossRef]
  20. Bento, H.F.; Vries, R.D.; Veldhuis, L.L. Aerodynamic Performance and Interaction Effects of Circular and Square Ducted Propellers. In Proceedings of the AIAA Scitech 2020 Forum, Orlando, FL, USA, 6–10 January 2020. [Google Scholar] [CrossRef] [Green Version]
  21. Sun, P.; Zhou, Z.; Guo, J. Numerical analysis for propulsion characteristics of ducted fans in different shapes. J. Aerosp. Power, 2021; in press. [Google Scholar] [CrossRef]
  22. Wang, H.; Wang, Y.; Zhou, F.; Liu, H. Optimization design of inlet for distributed ducted fan propulsion system based on panel method. J. Propuls. Technol. 2021, 42, 2465–2473. [Google Scholar]
  23. Rajagopalan, R.; Zhang, Z. Performance and Flow Field of a Ducted Propeller. In Proceedings of the 25th Joint Propulsion Conference, Monterey, CA, USA, 12–16 July 1989. [Google Scholar] [CrossRef]
  24. Hall, C.A.; Crichton, D. Engine Design Studies for a Silent Aircraft. J. Turbomach. 2007, 129, 1653–1662. [Google Scholar] [CrossRef]
  25. Liu, L.; Song, Y.; Chen, H.; Chen, F.; Cui, K. Design and optimization of half flush-mounted s-shaped inlet. J. Propuls. Technol. 2014, 35, 1303–1309. [Google Scholar] [CrossRef]
Figure 1. Duct profile of the 150 ducted fan.
Figure 1. Duct profile of the 150 ducted fan.
Aerospace 09 00165 g001
Figure 2. Geometry of the 150 ducted fan.
Figure 2. Geometry of the 150 ducted fan.
Aerospace 09 00165 g002
Figure 3. Picture of the 150 ducted fan during the test.
Figure 3. Picture of the 150 ducted fan during the test.
Aerospace 09 00165 g003
Figure 4. Schematic of the blade grid in the MRF method.
Figure 4. Schematic of the blade grid in the MRF method.
Aerospace 09 00165 g004
Figure 5. Schematic of the gird in the MSM.
Figure 5. Schematic of the gird in the MSM.
Aerospace 09 00165 g005
Figure 6. Comparison of the thrust of the ducted fan obtained by different methods.
Figure 6. Comparison of the thrust of the ducted fan obtained by different methods.
Aerospace 09 00165 g006
Figure 7. Schematic diagram of the distributed ducted fan system.
Figure 7. Schematic diagram of the distributed ducted fan system.
Aerospace 09 00165 g007
Figure 8. Schematic of the variable area nozzle.
Figure 8. Schematic of the variable area nozzle.
Aerospace 09 00165 g008
Figure 9. Airfoils of the blade of the ducted fan.
Figure 9. Airfoils of the blade of the ducted fan.
Aerospace 09 00165 g009
Figure 10. Distribution of the chord length and torsional angle of the rotor and stator: (a) Rotor; (b) Stator.
Figure 10. Distribution of the chord length and torsional angle of the rotor and stator: (a) Rotor; (b) Stator.
Aerospace 09 00165 g010
Figure 11. Schematic diagram of the calculation unit of the distributed ducted fan system.
Figure 11. Schematic diagram of the calculation unit of the distributed ducted fan system.
Aerospace 09 00165 g011
Figure 12. Variation of cruise performance of the distributed ducted fan unit with different outlet expansion angles: (a) Thrust; (b) Efficiency.
Figure 12. Variation of cruise performance of the distributed ducted fan unit with different outlet expansion angles: (a) Thrust; (b) Efficiency.
Aerospace 09 00165 g012
Figure 13. Comparison of surface pressure of the distributed ducted fan unit in different flight phases: (a) Hover; (b) Cruise.
Figure 13. Comparison of surface pressure of the distributed ducted fan unit in different flight phases: (a) Hover; (b) Cruise.
Aerospace 09 00165 g013
Figure 14. Comparison of the pressure distribution of the duct in different flight phases: (a) Hover; (b) Cruise.
Figure 14. Comparison of the pressure distribution of the duct in different flight phases: (a) Hover; (b) Cruise.
Aerospace 09 00165 g014
Figure 15. Comparison of the velocity image of the distributed ducted fan unit in cruise flight: (a) z = 0 m; (b) y = 0 m.
Figure 15. Comparison of the velocity image of the distributed ducted fan unit in cruise flight: (a) z = 0 m; (b) y = 0 m.
Aerospace 09 00165 g015
Figure 16. Distributed ducted fan system with the rectangular duct.
Figure 16. Distributed ducted fan system with the rectangular duct.
Aerospace 09 00165 g016
Figure 17. Configuration of the new distributed ducted fan system.
Figure 17. Configuration of the new distributed ducted fan system.
Aerospace 09 00165 g017
Figure 18. Inlet and outlet of the distributed ducted fan unit: (a) “Rectangle to circle” inlet; (b) “Circle to rectangle” outlet.
Figure 18. Inlet and outlet of the distributed ducted fan unit: (a) “Rectangle to circle” inlet; (b) “Circle to rectangle” outlet.
Aerospace 09 00165 g018
Figure 19. Schematic diagram of the distributed ducted fan unit.
Figure 19. Schematic diagram of the distributed ducted fan unit.
Aerospace 09 00165 g019
Figure 20. Schematic diagram of the inlet parameterization.
Figure 20. Schematic diagram of the inlet parameterization.
Aerospace 09 00165 g020
Figure 21. Multi-objective design process of the distributed ducted fan system.
Figure 21. Multi-objective design process of the distributed ducted fan system.
Aerospace 09 00165 g021
Figure 22. Schematic of the 1/4 calculation model of the multi-objective optimization.
Figure 22. Schematic of the 1/4 calculation model of the multi-objective optimization.
Aerospace 09 00165 g022
Figure 23. Pareto front of the inlet of different lengths.
Figure 23. Pareto front of the inlet of different lengths.
Aerospace 09 00165 g023
Figure 24. Pareto front of the inlet of different heights.
Figure 24. Pareto front of the inlet of different heights.
Aerospace 09 00165 g024
Figure 25. Duct profile of the inlet of different lengths.
Figure 25. Duct profile of the inlet of different lengths.
Aerospace 09 00165 g025
Figure 26. Duct profile of the inlet of different heights.
Figure 26. Duct profile of the inlet of different heights.
Aerospace 09 00165 g026
Figure 27. Comparison of the velocity image of different optimization results of different inlets in hover flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Figure 27. Comparison of the velocity image of different optimization results of different inlets in hover flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Aerospace 09 00165 g027
Figure 28. Comparison of the pressure image of different optimization results of different inlets in hover flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Figure 28. Comparison of the pressure image of different optimization results of different inlets in hover flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Aerospace 09 00165 g028
Figure 29. Comparison of the total pressure image at the entrance of the blade section of different inlets in hover flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Figure 29. Comparison of the total pressure image at the entrance of the blade section of different inlets in hover flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Aerospace 09 00165 g029
Figure 30. Comparison of the velocity image of different optimization results of different inlets in cruise flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Figure 30. Comparison of the velocity image of different optimization results of different inlets in cruise flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Aerospace 09 00165 g030
Figure 31. Comparison of the pressure image of different optimization results of different inlets in cruise flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Figure 31. Comparison of the pressure image of different optimization results of different inlets in cruise flight: (a) 50_80_70; (b) 100_80_70; (c) 50_80_80.
Aerospace 09 00165 g031
Figure 32. Schematic diagram of the ducted fan unit.
Figure 32. Schematic diagram of the ducted fan unit.
Aerospace 09 00165 g032
Figure 33. Schematic diagram of the relative position between the ducted fan and induced wing.
Figure 33. Schematic diagram of the relative position between the ducted fan and induced wing.
Aerospace 09 00165 g033
Figure 34. Comparison of the cruise performance of the distributed ducted fan system with different deflection angles of the induced wing at Ω = 8500 r/min: (a) Thrust; (b) Efficiency.
Figure 34. Comparison of the cruise performance of the distributed ducted fan system with different deflection angles of the induced wing at Ω = 8500 r/min: (a) Thrust; (b) Efficiency.
Aerospace 09 00165 g034
Figure 35. Comparison of sectional velocity image of the distributed ducted fan system with −12° deflection angle of the induced wing with different chord lengths: (a) 120 mm; (b) 150 mm.
Figure 35. Comparison of sectional velocity image of the distributed ducted fan system with −12° deflection angle of the induced wing with different chord lengths: (a) 120 mm; (b) 150 mm.
Aerospace 09 00165 g035
Figure 36. Comparison of the wake of the distributed ducted fan system with and without the induced wing in hover flight: (a) Without the induced wing; (b) With the induced wing.
Figure 36. Comparison of the wake of the distributed ducted fan system with and without the induced wing in hover flight: (a) Without the induced wing; (b) With the induced wing.
Aerospace 09 00165 g036
Figure 37. Comparison of duct profile of the distributed duct fan system before and after the improvement.
Figure 37. Comparison of duct profile of the distributed duct fan system before and after the improvement.
Aerospace 09 00165 g037
Figure 38. Comparison of the cruise performance of the distributed ducted fan system with different deflection angles of the induced wing before and after the improvement at Ω = 8500 r/min: (a) Thrust; (b) Efficiency.
Figure 38. Comparison of the cruise performance of the distributed ducted fan system with different deflection angles of the induced wing before and after the improvement at Ω = 8500 r/min: (a) Thrust; (b) Efficiency.
Aerospace 09 00165 g038
Figure 39. Sectional velocity image of the improved distributed ducted fan system in two flight phases: (a) Hover; (b) Cruise.
Figure 39. Sectional velocity image of the improved distributed ducted fan system in two flight phases: (a) Hover; (b) Cruise.
Aerospace 09 00165 g039
Figure 40. Hover and cruise performance of the improved distributed ducted fan system: (a) Hover; (b) Cruise.
Figure 40. Hover and cruise performance of the improved distributed ducted fan system: (a) Hover; (b) Cruise.
Aerospace 09 00165 g040
Table 1. Characteristic comparison of different methods at Ω = 10,500 r/min.
Table 1. Characteristic comparison of different methods at Ω = 10,500 r/min.
MethodTp (N)Td (N·m)T (N)ErrorTime (h)
EXP--32.116--
MRF15.32118.40433.7255.001%1.8
MSM15.32119.98435.3059.930%0.2
Table 2. Hover performance of the distributed ducted fan unit.
Table 2. Hover performance of the distributed ducted fan unit.
αe (deg)Ω (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)pl (kg·kw−1)
310,85040.36916.1951.0612.5241.0723.349
Table 3. Cruise performance of the distributed ducted fan unit.
Table 3. Cruise performance of the distributed ducted fan unit.
αe (deg)Ω (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)η
−586103.9988.7780.6310.9200.64235.136%
Table 4. Comparison of different optimization results of different inlets in hover flight.
Table 4. Comparison of different optimization results of different inlets in hover flight.
ResultTd (N)Df (N)
100_80_705.7020.258
70_80_705.6210.239
50_80_704.6830.243
50_80_755.4230.241
50_80_805.8370.240
Table 5. Comparison of different optimization results of different inlets in cruise flight.
Table 5. Comparison of different optimization results of different inlets in cruise flight.
ResultDd (N)Df (N)
100_80_700.5320.635
70_80_700.5040.597
50_80_700.4520.565
50_80_750.5860.590
50_80_800.6820.613
Table 6. Comparison of the hover performance of the ducted fan unit.
Table 6. Comparison of the hover performance of the ducted fan unit.
ConfigurationΩ (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)pl (kg·kw−1)
Initial system10,85040.36916.1951.0612.5241.0723.349
Design11,00040.19817.5181.1003.0011.0913.172
Table 7. Cruise performance of the designed distributed ducted fan system.
Table 7. Cruise performance of the designed distributed ducted fan system.
ConfigurationΩ (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)η
Initial system86103.9988.7780.6310.9200.64235.136%
Design81803.9948.6620.5931.1180.60139.313%
Table 8. Comparison of the hover performance of the distributed ducted fan system with different deflection angles of the induced wing.
Table 8. Comparison of the hover performance of the distributed ducted fan system with different deflection angles of the induced wing.
αe (deg)Ω (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)pl (kg·kw−1)
011,00040.06816.4081.0872.4251.0963.200
111,00040.34716.2881.0832.3921.0973.234
211,00039.62316.2011.0822.3431.0963.179
Table 9. Comparison of the hover performance of the distributed ducted fan system with different deflection angles of the induced wing after the improvement.
Table 9. Comparison of the hover performance of the distributed ducted fan system with different deflection angles of the induced wing after the improvement.
αe (deg)Ω (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)pl (kg·kw−1)
011,00040.42216.7751.0932.5821.0973.211
111,00040.44716.5891.0902.4971.0973.221
211,00040.42016.4821.0882.4321.0963.225
311,00040.26716.4191.0872.4101.0953.216
Table 10. Comparison of the hover performance of different distributed ducted fan systems.
Table 10. Comparison of the hover performance of different distributed ducted fan systems.
ConfigurationΩ (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)pl (kg·kw−1)
Initial system10,85040.36916.1951.0612.5241.0723.349
Design11,00040.34716.2881.0832.3921.0973.234
Improved design11,00040.42016.4821.0882.4321.0963.225
Table 11. Comparison of the cruise performance of different distributed ducted fan systems.
Table 11. Comparison of the cruise performance of different distributed ducted fan systems.
ConfigurationΩ (rpm)T (N)Tr (N)Qr (N·m)Ts (N)Qs (N·m)η
Initial system86103.9988.7780.6310.9200.64235.136%
Design81803.9948.6620.5931.1180.60139.313%
Improved design80503.9918.3990.5741.0850.58041.240%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, J.; Zhou, Z. Multi-Objective Design of a Distributed Ducted Fan System. Aerospace 2022, 9, 165. https://doi.org/10.3390/aerospace9030165

AMA Style

Guo J, Zhou Z. Multi-Objective Design of a Distributed Ducted Fan System. Aerospace. 2022; 9(3):165. https://doi.org/10.3390/aerospace9030165

Chicago/Turabian Style

Guo, Jiahao, and Zhou Zhou. 2022. "Multi-Objective Design of a Distributed Ducted Fan System" Aerospace 9, no. 3: 165. https://doi.org/10.3390/aerospace9030165

APA Style

Guo, J., & Zhou, Z. (2022). Multi-Objective Design of a Distributed Ducted Fan System. Aerospace, 9(3), 165. https://doi.org/10.3390/aerospace9030165

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop