Next Article in Journal
Study on Failure Criteria and the Numerical Simulation Method of a Coal-Based Carbon Foam under Multiaxial Loading
Previous Article in Journal
Crew Interaction in Extended Space Missions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parametric Investigation of Leading-Edge Slats on a Blended-Wing-Body UAV Using the Taguchi Method

by
Spyridon Antoniou
1,2,
Stavros Kapsalis
1,2,
Pericles Panagiotou
1,2,* and
Kyros Yakinthos
1,2
1
Laboratory of Fluid Mechanics and Turbomachinery, Department of Mechanical Engineering, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
2
UAV Integrated Research Center (UAV-iRC), Center for Interdisciplinary Research and Innovation (CIRI), Aristotle University of Thessaloniki, 57001 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Submission received: 21 July 2023 / Revised: 9 August 2023 / Accepted: 15 August 2023 / Published: 17 August 2023

Abstract

:
The current study investigated the effect of leading-edge slats on the longitudinal stability at high angles of attack of a Blended-Wing-Body (BWB) Unmanned Air Vehicle (UAV). Using a Design of Experiments (DOE) approach and, specifically, the Taguchi method, four leading-edge slat design parameters were investigated on three different levels. These parameters were the slat semi-span, the rotation of the slat element, the extension forward of the leading edge and the downward drop below the leading edge. An L9 orthogonal array (OA) was used to investigate the influence of these key design parameters using three performance criteria, namely the angle at which pitch break occurs, the corresponding speed and the distance between the Neutral point of each configuration and the Neutral point of the reference platform. The investigation was conducted by using high-fidelity Computational Fluid Dynamics (CFD) methods for each of the nine configurations defined by the L9 OA, over a range of angles of attack between −4 and 16 degrees. Based on these results, and using a Signal-to-Noise ratio (SNR) analysis, two combinations were eventually derived, one that optimized pitch break angle and speed and one that optimized longitudinal stability. Finally, the Pareto Analysis of Variance (ANOVA) technique was conducted to define the contribution of each of the six design parameters on the selected performance criteria. More specifically, the semi-span seemed to have the most significant effect on pitch break angle and speed, whereas the rotation of the slat element was the most important parameter with regard to static stability.

1. Introduction

The Blended-Wing-Body (BWB) was introduced by Liebeck [1] as a novel and efficient configuration for high-speed subsonic commercial aircraft. Following this initial suggestion, several researchers investigated this configuration both for commercial airliner applications [2,3,4], as well as for cargo transport [5] and Unmanned Air Vehicle (UAV) applications [6,7,8]. Especially concerning the latter, and even though UAVs operate at lower Reynolds numbers than jet transports, adopting the BWB configuration can potentially lead to an enhanced aerodynamic efficiency compared to a typical tube-wing aircraft. Another big advantage is its greater internal volume, which allows for the installation of larger and heavier payloads.
BWBs are, by definition, tailless configurations; therefore, their design calls for a unique set of layout choices. Sweeping the wing aftwards is the most common choice, as it is even visible to the untrained eye [4,8]. In contrast to theoretical aerodynamics, where sweep is used to alleviate the compressibility effects and reduce wave drag [9], wing sweep is used in this case regardless of the operating Mach number to enhance longitudinal and lateral stability, by moving the aerodynamic to the rear and by providing a natural dihedral effect, respectively [10,11]. Moreover, by placing the control surfaces at the outboard section of the swept wing, the layout designer can achieve adequate moment arm and, consequently, controllability.
This design choice may address many key issues of BWB design but is also linked with a critical disadvantage. That is, the wing develops a pitch-up tendency at high angles of attack that can lead to total loss of controllability and lift [10,11]. This pitch-up tendency is commonly defined as unstable pitch break and often in swept wings with high aspect ratio (AR) [12,13]. The two main mechanisms that cause this phenomenon are the spanwise flow, i.e., the flow developing from the root towards the tip of the wing, which enhances the flow separation towards the tip, and the local upwash that is produced due to the swept geometry, leading the local sections close to the wingtip to reach their maximum lift coefficients sooner and then stall. Both mechanisms cause the tip to stall first [14,15,16]. On the one hand, as the flow at the inboard sections remains attached, the lift production keeps increasing as the angle of attack increases. On the other hand, though, tip stalling causes the center of pressure of the wing to move forward, affecting the pitching moment of the aircraft on a progressively negative way and, ultimately, leading to unstable pitch break. The control of the aircraft is also affected in a negative way, as the control surfaces located towards the tip progressively lose their effect. The aforementioned issues pose a serious threat to BWB operation at the flight segments where the aircraft operates high angles of attack, like take-off and landing.
To counter this destructive phenomenon, several methods of active and passive flow control have been proposed. The wing fences [17,18] stop the spanwise flow development by physically obstructing its path, forcing the flow to go parallel with the fence. The vortex generators [19], the vortilons [20] and the dog-tooth configurations [21] deal with the spanwise flow by creating a vortex over the suction side of the wing, which has a similar effect with the wing fences, creating a pseudo-fence. While all the above can lead to positive results, they are means of passive flow control and, as such, affect the airflow during the whole flight. This can lead to a drag increase which can compromise the aerodynamic efficiency and the performance of the aircraft.
Another possible technique to counter the pitch break is the use of leading-edge slats [13,14,22]. Deploying slats does not only lead to the creation of vortices which work as pseudo-fences but by increasing the speed of the local airflow as well. This way, the airflow remains attached to the wing at a greater percentage of its span, tip stall is delayed and the occurrence of an unstable pitch break is either avoided or delayed. The advantage of leading-edge slats is that, as a means of active flow control technique, they can be extended on-demand during take-off and landing, to counter the risk for unstable pitch break at high angles of attack. For the remainder of the flight, they remain retracted and have little to no effect on aerodynamic performance. Their main drawback is the weight and complexity penalty that comes along with slat integration [23]; therefore, if a layout designer is to consider this technique for their trade studies, they must have a set of trends and guidelines to properly design and exploit the full potential of their advantages.
The majority of research on leading-edge slats is focused mostly on the simulation of two-dimensional multi-element airfoils [24,25], their application on wing turbines [26,27], or their effect on the noise [28]. When the use of slats on aircraft wings is investigated [22], researchers focus on lift enhancement and the pitch attitude is often neglected, as a conventional aircraft can deal with the unstable pitch break of the wing by using its tail surfaces.
In the present work, a detailed parametric investigation of leading-edge slat configurations was performed. The objectives can be briefly described as follows:
  • To provide trends for the sizing of leading-edge slats on a tailless, BWB UAV configuration;
  • To incorporate critical aspects of longitudinal stability in the investigation;
  • To propose an efficient CFD-based methodology that can allow future researchers to replicate the analysis on different platforms during the conceptual and preliminary design phases, where both the time and the resources are usually limited.
The effect of key slat design parameters was examined using computational modeling to identify its respective effect on the angle where the pitch break appears, the speed at which the UAV encounters pitch-break (defined as Vstall in the current study) and the distance between the Neutral point of each configuration and the Neutral point of the reference platform (ΔxNP). By using the Taguchi methodology, four design parameters were examined on three different levels. The examined parameters were the slat semi-span (bs/2), the rotation of the slat element (δs), the extension forward of the leading edge (Δx) and the downward drop below the leading edge (Δy). Additionally, using the Pareto Analysis of Variance (ANOVA), the percentage contribution of the design parameters was defined for each of the optimized performance criteria.

2. Baseline Platform

The baseline platform utilized in this study is a fixed-wing, BWB (Blended-Wing-Body) UAV configuration. It closely resembles the RX-3 prototype, which underwent a comprehensive conceptual and preliminary design study [8]. The development process of the RX-3 was supported by applied aerodynamics investigations and parametric studies carried out at the Laboratory of Fluid Mechanics and Turbomachinery (LFMT) within the School of Mechanical Engineering at the Aristotle University of Thessaloniki (AUTH) in Greece. The design and sizing of the UAV, along with the aerodynamic, stability, and performance analyses, were conducted using established design methods and dedicated in-house tools. Moreover, the entire development process adhered to relevant airworthiness regulations. For reference, Figure 1 displays the external layout of the baseline platform, which also serves as the reference configuration in the context of the EURRICA project layout. The sweep of the wing section of the UAV was selected at 33 degrees, while the AR equaled to 8. More details about the design procedure, the tools and the layout can be found in [8]. The operating specifications of the baseline platform resembled those of tactical UAVs, as the latter are described in [29], and they are presented with detail in [8,30]. Some basic geometric, mass and performance data that are relevant to this study are presented in Table 1.

3. Tools and Methods

3.1. Slat Design

According to the corresponding literature [11,16,31,32], the key slat design parameters are the slat chord (cs), the extension forward of the leading edge (Δx or α), the downward drop below the leading edge (Δy or b), the size of the gap at the outlet of the slat (c), the rotation of the element (δS), the airfoil of the slats and the semi-span of the slats (bs/2), as can be seen in Figure 2.
The size of the gap at the outlet of the slat was directly linked to the positioning of the slat element resulted from Δx, Δy and δS and, as such, did not serve as an independent variable in the present work. The two dimensional (2D) profile of the slat (the airfoil and the chord of the slat) was linked to the corresponding profile of the main wing. In this work, the latter derived from a given baseline platform; therefore, no further analysis was conducted on the 2D slat profile. The placement of the inboard frame of the slat on the wing was heavily affected by structural limitations, i.e., by the spanwise locations that feature a wing rib that has been properly sized to accommodate the slat structure. However, as the structural analysis of the slat is beyond the scope of this work, the root of the wing served as the inboard frame and was kept constant during the whole study.
As a result, the number of design parameters of the slat elements that were used in the current trade study was reduced from seven to four, allowing for a more focused study and leading to a major reduction in the total computational time. The orthogonal array of the Taguchi methodology was defined by the semi-span of the slats (A), the rotation of the element (B), the extension forward of the leading edge (C) and the downward drop (D).
Τo produce the different slat configurations for all the combinations of the selected geometric parameters, the corresponding CAD models were created using commercial 3D CAD software. The main body of the UAV was created with the use of an in-house parametric tool [8] and the slat element was then subtracted from the main wing as a distinct solid body. In the end, the slat configuration was moved and rotated appropriately to generate the desired configuration. In Figure 3, the top view of two slat configurations is indicatively presented, as a result of the CAD modeling. Configuration 2 is presented on the right with blue, while configuration 6 is on the left with red.

3.2. Taguchi Methodology

In contemporary aircraft design, the widespread utilization of advanced computational methods and experiments is prevalent. These techniques enable crucial optimization and trade studies to take place, but they come with limitations due to their considerable demands on computational time, resources and human effort. Consequently, the design space is often constrained, as highlighted in the examination of winglet sizing studies that are presented in [30]. The complexity and the cost of the trade studies is further increased when rapid, lower-fidelity methods (e.g., vortex lattice methods) cannot be trusted to model complex 3D configurations, as is the case of the slats and their effect on the pitch break.
When the whole design space is under investigation, for a full factorial generic optimization problem, the number N of the possible combinations that occur for L number of levels of the design parameters and m number of design parameters is equal to the expression shown in Equation (1).
N = L m
To avoid the expenses associated with a full factorial analysis and simultaneously explore the entire design space, Taguchi introduced a standardized Design of Experiments (DOE) method [33]. This widely adopted approach in the industry facilitates the development of high-quality products by maximizing their robustness and making them less susceptible to uncontrollable factors, all at a relatively lower cost. Taguchi’s method achieves this by employing a special set of Orthogonal Arrays (OAs) that define the minimum required “experiments” for the chosen design parameters (Table 2). These experiments ensure that the conclusions drawn from the reduced “experimental” space remain valid across the entire design space. ensuring that the conclusions derived from the reduced “experimental” space are valid over the entire design space. It is important to clarify that the use of the Taguchi method itself is not a novel concept, as it is a well-established DOE method extensively employed in various industries. However, the proposed overall approach, which combines Taguchi with ANOVA to determine the percentage contribution of the leading-edge slats design parameters to stability and low-speed performance, is truly unique and original. This specific combination has not been previously discussed in published aeronautical literature.
The initial step in creating an orthogonal array involves determining the total degrees of freedom, which establishes the minimum number of required “experiments” to investigate the chosen design parameters. For each design parameter (A, B, C, etc.), the degrees of freedom are calculated as the number of levels minus one (Equation (2)), where i stands for the respective parameter (A, B, C, etc.). Therefore, the overall degrees of freedom are obtained by summing up the degrees of freedom of the main effects of the design parameters, along with an additional degree of freedom associated with the overall mean [34].
d f i = n u m b e r   o f   l e v e l s   o f   p a r a m e t e r   i 1
The second step involves determining the suitable layout for the Orthogonal Array (OA), which is determined by the total degrees of freedom and the number of design parameters. The number of rows in the OA must be at least equal to the total degrees of freedom, while the number of columns is based on the factors and their levels. The design engineers select the specific design parameters, considering relevant aspects in the field, theory, and common practice [10,35]. Concerning the current study, the design parameters defined in the previous section are presented in Table 3. The values of each design parameter level are selected based on reference platform and manufacturability restrictions and general guides for slats design [11,16,31,32]. In conclusion, the combination of the design parameters and their levels (four parameters with three levels each) led to an L9 orthogonal array for the current trade studies (Table 4). Thus, nine different slats configurations were examined and their effect on the aerodynamic performance of the examined platform was extracted. It is essential to note that by employing the Taguchi method, the total number of cases to be examined was significantly reduced from 81 to just 9, in contrast to the full factorial approach. This substantial reduction in cases drastically decreased the effort required for CAD and CFD analysis.

3.3. CFD Methodology

The aerodynamic behavior of each configuration was investigated using high-fidelity, CFD computations. In the present study, the steady-state Reynolds Averaged Navier–Stokes (RANS) equations were solved, with the adoption of an appropriate turbulence model, necessary for their closure. The model of Spalart and Allmaras (SA) [36] was selected, the one-equation eddy viscosity model. One transport equation (Equation (3)) was used for the calculation of the viscosity variable ν ~ . The use of the model is widely spread for external aerodynamics flows and has been experimentally validated in low freestream Reynolds number applications.
t ρ ν ~ + x i ρ ν ~ u i = G ν + 1 σ ν ~ x j μ + ρ ν ~ ν ~ x j + C b 2 ρ ν ~ x j 2 Y ν + S ν ~
The CAD models were imported in BETA ANSA pre-processing software (v19.1.0, Root, Switzerland) to properly discretize the surrounding control volume and generate the corresponding computational mesh (Figure 4). As the UAV was symmetrical, only half of the model was employed in the simulations, resulting in approximately 28 million computational cells for all tested cases. This study employed a combination of structured and unstructured meshing techniques. An unstructured mesh philosophy was used to discretize the majority of the computational domain, as it provided greater flexibility for discretizing the flow domain around the complex 3D configurations, such as the region between the slat element and the main body (Figure 4b). On the contrary, a structured mesh was applied in close proximity to the UAV surface. This may locally cause an increase in the number of elements, but was used to better capture the flow and its sharp changes (steep gradients) inside the boundary layer (Figure 5b). More specifically, concerning the boundary layer discretization, the y+ parameter did not exceed a value of 5 at any point of the model surface and, in general, was kept below 1. The reason for this was to ensure that the viscous sublayer region was adequately modeled [37]. Apart from thickness, the size of surface elements ranged from 1 mm at areas such as the trailing edge, to 30 mm at regions like the middle part of the main body. Concerning the leading edge curvature resolution, it was discretized into a minimum of 25 elements to ensure accurate modeling, both at the slat and the wing leading edge. The mesh size was based on the grid dependency study presented in [30] for the baseline configuration. The results are indicatively shown in Figure 6 for the Lift and Drag coefficients, with the red marker dictating the selected grid for this study.
The solution was performed with the Ansys CFX software (Release 18.2, Academic Multiphysics Campus Solution, Canonsburg, PA, USA), using a coupled pressure-based solver, and second-order spatial discretization schemes for the continuity, x-, y- and z-momentum and viscosity variable transport equations. The examined configurations were modeled in a wide range of angles of attack, which spanned between −4 degrees and 16 degrees, so that the pitch break occurred. The atmospheric conditions and the flight speed that were used for this investigation corresponded to take-off and landing segments of a typical flight, where the use of slats was expected, i.e., sea level conditions. It must be noted at this point that, for research purposes, a follow-up analysis could also be conducted at different flight altitudes to compare flight performance. However, the baseline configuration of this study was not expected to deploy slats at another flight phase and, consequently, this follow-up analysis was not carried out in this paper.
The boundary and initial conditions are summed up in Table 5, also including the turbulence boundary conditions at the inlet [38]. It is important to address the computational modeling’s fidelity at this stage. While utilizing a more sophisticated approach, such as employing the Unsteady RANS and advanced turbulence models or a Reynolds Stress Model (RSM), could potentially produce more accurate results in predicting the impact of slats on aerodynamics and stability, the authors of the current work chose not to focus on providing such detailed outcomes. As stated in the introduction, the primary goal was to conduct a comparative research study that strikes a balance between adequate fidelity (using a well-established turbulence model for aeronautical applications) and computational efficiency. This approach allows the identification of key trends and ensures the overall method’s reproducibility during the conceptual and preliminary design phases of UAV configurations, where time and resources are limited.

3.4. SNR and ANOVA

The Signal-to-Noise Ratio (SNR) was suggested by Taguchi for the analysis of the results. The SNR generally represents the efficiency factor and a quality measure connected with the loss function. By maximizing the SNR, the losses related to the procedure were minimized. The “experimental” observations were converted to SNR, following three basic, different definitions depending on the desirable direction of the results (Equations (4)–(6)). The logarithmic transformation of the observations in SNR strengthens the prediction of the performance criteria improvement.
Smaller   the   better :   S N R = 10 log 10 1 n y 2
Nominal   the   better :   S N R = 10 log 10 1 n Y ¯ s y 2
Bigger   the   better :   S N R = 10 log 10 1 n 1 y 2
In this specific case, the term y refers to the estimated performance criteria for each configuration and n to the number of configurations under investigation. Furthermore, Y ¯ is the mean value and S Y 2 is the variance of each performance criterion (Equations (7) and (8)) [39].
Y ¯ = 1 k y k
S Y 2 = 1 k ( y Y ¯ ) 2 k
Concerning the performance criteria of the current study, for the optimization of the pitch-break angle, the ‘Bigger the better’ definition was used, whereas for the optimization of the Vstall and the ΔxNP, the ‘Smaller the better’ definition was applied. For each performance criterion, the combination of the design parameters with the highest SNR was the optimum one and corresponded to the best performance with the minimum variation.
Additionally, the study included an ANOVA [34,40] to analyze the findings and draw conclusions concerning the optimal combination for each performance criterion. ANOVA is a statistical technique primarily employed for analyzing “experimental” data and examining the influence of design parameters on the overall variation of the results. In this research, a simplified version of ANOVA known as Pareto ANOVA [40] was utilized. Based on the Pareto principle, this method provides a quick and straightforward way to analyze the results of parametric trade studies. Importantly, the Pareto ANOVA does not require an F-test, allowing for the determination of the statistical significance of each design parameter for every performance criterion.
In the current study, an ANOVA was conducted for the angle where the pitch break appeared, the speed at which the UAV had pitch-break (which is named as Vstall in the current study) and the distance between the Neutral point of each configuration and the Neutral point of the reference platform (ΔxNP), taking into concern the SNR analysis results, yielding the importance and the percentage contribution of slat semispan (bS/2, Δs, Δx and Δy).
Following the methodology presented in [40], let i stand for the number of parameters and j stand for the number of levels of each design parameter so that i and j range from 1 to 4 and 1 to 3, respectively. Hence, following the well-established ANOVA methodology, the contribution of each design parameter is calculated using Equation (9).
P e r c e n t a g e   c o n t r i b u t i o n = S S i S S · 100
The terms S S i and S S represent the sum of squares of each design parameter and the total sum of squares, respectively, due to variation about the overall mean. The exact values of S S i and S S are calculated using Equations (10) and (11).
S S i = 1 3 ( S N R i j S N R ¯ ) 2
S S = 1 27 ( S N R i S N R ¯ ) 2
Finally, S N R i j is the average SNR of the ith parameter at the jth level, and the overall mean SNR is calculated using Equation (12).
S N R ¯ = 1 27 S N R i 27
The above procedure is conducted for each one of the performance criteria, yielding two percentage contributions in the current study.

4. Results

The basic aerodynamic coefficients, such as the lift and drag forces and the pitching moment, were extracted from the CFD simulations for the total of the configurations of the nine configurations. The lift and the drag coefficients were calculated using the reference area shown in Table 1. Moreover, Table 1 also presents the mean aerodynamic chord that was used for the calculation of the pitching moment coefficient. The lift coefficient, the drag polars of these configurations and the trim diagrams are presented in Figure 7, Figure 8 and Figure 9, compared with the baseline geometry without slats. The center of gravity at take-off was used for this analysis, which was located at 1790 mm from the UAV apex. As the UAV was designed to be inherently stable during the whole duration of a typical flight, the same conclusions can be extracted for the landing configuration as well.
The use of slats led to a slight reduction in the slope of the lift coefficient in the linear region of the curve. This slope reduction became more evident as the rotation of the slat element increased, reaching as much as 16% for configuration 9. This can be attributed to the fact that, at small angles of attack, the slat suction side was placed almost perpendicularly to the airflow and did not contribute to lift production. As the angle of attack increased, the gap of the slat became more effective, as the UAV did not reach its aerodynamic stall.
The effect of the slats on the lift coefficient became evident at angles of attack higher than 8 degrees. As can be seen in Figure 7, for all the of the examined configurations, the maximum lift coefficient either remained the same or was slightly increased, as was the case with configuration 4. However, it would be wrong to consider these angles of attack as operational because, as will be discussed in the next paragraph, pitch break already occurred. The total drag of the UAV increased because of the slats extension (Figure 8), but the authors did not consider this to be a major issue for the performance of the vehicle, as the devices are expected to be used during take-off and landing flight segments.
Figure 9 primarily shows that, as expected, the use of slats had a major impact on the angle of attack where pitch break occurred. The pitch break remained unstable, the stall progression of some configurations was more gradual and “soft” [9], depending on the corresponding design parameters. For example, the pitch break on configuration 2 (Figure 10, top) was delayed for about 2 degrees and seemed to be smoother than the baseline geometry, as the pitching moment caused smaller deviations to appear in this range. In contrast, configuration 8 (Figure 10, bottom) was expected to have pitch break between 4 and 8 degrees of angle of attack, while the Neutral point of the UAV moved forwards, leading closer to neutral static stability. An important note at this point is the fact that even though all configurations affected the static stability of the UAV, the use of slats did not lead to instability at any point.
The SNR analysis was conducted to define the optimum combination of the design parameters for each performance criterion (Pitch break angle, Vstall and ΔxNP). From the SNR analysis, the SNR of each configuration for each performance criterion was calculated (Table 6). It was found that the mean SNR for Pitch break angle was 18.99 Db, for the Vstall, it was −30.02 dB and for ΔxNP, it was 26.21 dB.
For the optimization of the pitch-break angle, the ‘Bigger the better’ definition of SNR was used, while for the optimization of the Vstall and the ΔxNP, the ‘Smaller the better’ definition was selected. The results of the SNR analysis for each parameter individually are presented in Table 7, Table 8 and Table 9. Delta was equal to the difference between the maximum and the minimum SNR value of each parameter. The rank of each design parameter revealed which one affected the response characteristic more and it was defined by the fraction of each parameter Delta divided by the sums of all the parameters’ Deltas.
The effect of the design parameters on the performance criteria is shown in plot form in Figure 11, Figure 12 and Figure 13. The optimum level of each parameter for each response characteristic was the one with the maximum SNR. The pitch break angle and the Vstall had similar behavior as they both took their maximum value at the first level for A (bS/2), B (δS) and C (Δx), and at the third level for D (Δy). ΔxNP performance metric took its maximum value at the first level for A (bS/2), B (δS) and D (Δy), and at the second level for C (Δx).
The optimum design parameter levels combinations for each examined performance criterion are summed up in Table 10. The combination A1B1C1D3 optimized the pitch break angle, according to the “Bigger is Better” condition and the corresponding SNR analysis. In a similar way, the same combination (A1B1C1D3) was the optimum for the Vstall metric, while the A1B1C2D1 was the optimum combination for ΔxNP, by the “Smaller is Better” definition of the SNR. The smaller semi-span and rotation of the slat elements seemed to be the preferable choice, while the translation of the slat element forward from the leading-edge can be kept at the smaller percentages for the optimization of the stability at this range of angles of attack.
The Pareto ANOVA results are presented in Figure 14, where the effect of design parameters on each performance criterion is specifically shown. According to the ANOVA, the pitch-break angle was mostly affected by the semi-span of the slat element (A) and the Δy (D), while the other two parameters had a similar, smaller contribution. Similarly, the semi-span (A) was the main factor that affected the Vstall metric, with the second most effective parameter being the Δy (D). For this metric, the rotation of the slat (B) contributed more than the Δx (C). The ΔXNP was affected mostly by the δs (B), while the semi-span (A) and the Δy (D) contributed approximately the same at this metric. Finally, the Δx (C) seemed to have the smallest effect compared to the other parameters.

5. Conclusions

This study presented a parametric, CFD-aided investigation of leading-edge slats conducted on a BWB UAV baseline platform. The slats were used primarily to enhance the longitudinal stability of the UAV at high angles of attack and expand its operating range. Using the Taguchi method, four different design parameters were examined and their effect on pitch break, static stability and overall performance of the UAV at high angles of attack was investigated. The various configurations were analyzed using computational modeling and the corresponding aerodynamic and stability coefficients were calculated. The combinations that maximized the pitch break angle and minimized Vstall and ΔxNP, respectively, were defined as a result of SNR analysis. Finally, a Pareto ANOVA was implemented to calculate the percentage contributions of the design parameters on each performance criterion. The main conclusions of this study are listed as follows:
  • The use of an optimized slat configuration can be used to considerably expand the operating range of a tailless UAV configuration. Using non-optimized slat design values, however, can lead in a reduction in aerodynamic performance.
  • The use of slats affects the slope of the pitching moment curve, but it does not result in an unfavorable behavior (static instability).
  • The optimum combination of design parameters for the maximization of pitch break angle is A1B1C1D3 and the optimum level values are bs/2 = 45%, δs = 0°, Δx = 1% and Δy = 5%.
  • The optimum combination of design parameters for the minimization of Vstall is A1B1C1D3 and the optimum level values are bs/2 = 45%, δs = 0°, Δx = 1% and Δy = 5%.
  • The optimum combination of design parameters for the minimization of ΔxNP is A1B1C2D1 and the optimum level values are bs/2 = 45%, δs = 0°, Δx = 3% and Δy = 1%.
  • The parameter that affects the pitch break angle and the Vstall the most is the bs/2 with a contribution factor of 38.34% and 35.24%, respectively.
  • The parameter that affects the ΔxNP the most is the δs with a contribution factor of 42.76%.
  • Τhe use of high-fidelity CFD simulations for the trade studies of the geometry of the leading-edge slats combined with the DOE resulted in a considerably lower amount of computational time. The design space is limited from eighty-one to just nine configurations by the implementation of the Taguchi methodology, resulting in an 89% reduction in the total computational effort.
The next steps of the current research will be focused on the addition of more geometric parameters in the parametric investigation, like the airfoil of the slat element. The use of a full-span slat will also be considered. The authors would also like to investigate the optimal slat design for each flight segment, like take-off, climb, approach and landing. This could also involve investigations at different flight altitudes to compare the slat performance. Finally, the structural implementations of the use of leading-edge slats can also be investigated, with emphasis on the actuation and the weight penalty for each of the alternative designs.

Author Contributions

Conceptualization, S.A. and S.K.; methodology, P.P. and S.K.; software, S.A.; validation, S.K., P.P. and K.Y.; investigation, S.A.; resources, S.A.; data curation, S.A. and S.K.; writing—original draft preparation, S.A. and S.K.; writing—review and editing, P.P. and K.Y.; visualization, S.A. and S.K.; supervision, P.P.; project administration, P.P.; funding acquisition, P.P. and K.Y. All authors have read and agreed to the published version of the manuscript.

Funding

The work presented in this paper was funded by the Hellenic Foundation for Research and Innovation (HFRI), partially in the framework of the EURRICA research project under the “2nd Call for H.F.R.I. Research Projects to support Post-Doctoral Researchers” (Project Number: 00745) and partially in the framework of the 3rd Call for HFRI PhD Fellowships (Fellowship Number: 6179).

Data Availability Statement

Data available on request due to restrictions related to the details of the baseline platform. The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liebeck, R.H. Design of the Blended Wing Body Subsonic Transport. J. Aircr. 2004, 41, 10–25. [Google Scholar] [CrossRef]
  2. Qin, N.; Le Moinge, A. Spanwise Lift Distribution for Blended Wing Body Aircraft. J. Aircr. 2005, 42, 356–365. [Google Scholar] [CrossRef]
  3. Reist, T.A.; Zingg, D.W. Aerodynamic design of blended wing body and lifting-fuselage aircraft. In Proceedings of the 34th AIAA—Applied Aerodynamics Conference, Session: Innovative Aerodynamic Concepts Designs, Washington, DC, USA, 13–17 June 2016. AIAA 2016-3874. [Google Scholar]
  4. Faggiano, F.; Vos, R.; Baan, M.; van Dijk, R. Aerodynamic Design of a Flying V Aircraft. In Proceedings of the 17th AIAA Aviation Technology, Integration, and Operations Conference, Denver, CO, USA, 5–9 June 2017; American Institute of Aeronautics and Astronautics Inc. (AIAA): Reston, VA, USA, 2017. AIAA 2017-3589. [Google Scholar] [CrossRef]
  5. Wakayama, S. Blended-wing-body optimization problem setup. In Proceedings of the 8th Symposium on Multidisciplinary Analysis and Optimization, Long Beach, CA, USA, 6–8 September 2000; p. 4740. [Google Scholar]
  6. Suewatanakul, S.; Porcarelli, A.; Olsson, A.; Grimler, H.; Chiche, A.; Mariani, R.; Lindbergh, G. Conceptual Design of a Hybrid Hydrogen Fuel Cell/BatteryBlended-Wing-Body Unmanned Aerial Vehicle—An Overview. Aerospace 2022, 9, 275. [Google Scholar] [CrossRef]
  7. Cummings, R.; Liersch, C.; Schütte, A. Multi-Disciplinary Design and Performance Assessment of Effective, Agile NATO Air Vehicles. Aerosp. Sci. Technol. 2020, 99, 105764. [Google Scholar] [CrossRef]
  8. Panagiotou, P.; Mitridis, D.; Dimopoulos, T.; Kapsalis, S.; Dimitriou, S.; Yakinthos, K. Aerodynamic design of a tactical Blended-Wing-Body UAV for the aerial delivery of cargo and lifesaving supplies. In Proceedings of the AIAA Scitech 2020 Forum, Orlando, FL, USA, 6–10 January 2020. [Google Scholar]
  9. Anderson, J.D. Fundamentals of Aerodynamics, 6th ed.; WCB/McGraw-Hill: Boston, MA, USA, 2011; ISBN 9780073398105. [Google Scholar]
  10. Raymer, D. Aircraft Design: A Conceptual Approach, Sixth Edition; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2018. [Google Scholar]
  11. Gudmundsson, S. General Aviation Aircraft Design, 2nd ed.; Butterworth-Heinemann: Oxford, UK, 2022. [Google Scholar]
  12. Shortal, J.A.; Maggin, B. Effect of Sweepback and Aspect Ratio on Longitudinal Stability Characteristics of Wings at Low Speeds; Defense Technical Information Center: Fort Belvoir, VA, USA, 1946; NACA TN-1093.
  13. Furlong, G.C.; Mchugh, J.G. A Summary and Analysis of the Low-Speed Longitudinal Characteristics of Swept Wings at High Reynolds Number; Defense Technical Information Center: Fort Belvoir, VA, USA, 1957; NACA TR-1339.
  14. Koven, W.; Graham, R.R. Wind-Tunnel Investigation of High-Lift and Stall-Control Devices on a 37 Degree Sweptback Wing of Aspect Ratio 6 at High Reynolds Numbers; Defense Technical Information Center: Fort Belvoir, VA, USA, 1948; NACA RM-L8D29.
  15. Haines, A.B. Some Notes on the Flow Patterns Observed over Various Swept-Back Wings at Low Mach Numbers (in the R.A.E. 10ftx7ft High Speed Tunnel; Cranfield University: Cranfield, UK, 1954; R&M No. 3192. [Google Scholar]
  16. Hoerner, S.F.; Borst, H.V. Fluid Dynamic Lift; Hoerner, L.A., Ed.; Hoerner Fluid Dynamics: Bricktown, NJ, USA, 1985. [Google Scholar]
  17. Dickson, J.K.; Sutton, F.B. The Effect of Wing Fences on the Longitudinal Characteristics at Mach Numbers Up to 0.92 of a Wing-Fuselage-Tail Combination Having a 400 Sweptback Wing with Naca 64a Thickness Distribution; National Advisory Committee for Aeronautics: Washington, DC, USA, 1955. [Google Scholar]
  18. Walker, M.; Bons, J. The Effect of Passive and Active Boundary-Layer Fences on Swept-Wing Performance at Low Reynolds Number. In Proceedings of the AIAA Scitech 2018 Forum, Kissimmee, FL, USA, 8–12 January 2018. [Google Scholar]
  19. Goharshadi, M.; Mirzaei, M. Delaying the Stall of a Low-Wing Aircraft Using A Novel Powerful Vortex Generator. Inventions 2022, 7, 95. [Google Scholar] [CrossRef]
  20. Phillips, E.; Wygnanski, I.J.; Menge, P.M.; Taubert, L. Passive and Active Leading Edge devices on a simple swept back wing. In AIAA Aviation 2019 Forum; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2019; p. 3393, NACA-RM-L52C25. [Google Scholar]
  21. Jaquet, B.M. Effects of Chord Discontinuities and Chordwise Fences on Low-speed Static Longitudinal Stability of an Airplane Model Having a 35 Degree Sweptback Wing; Defense Technical Information Center: Fort Belvoir, VA, USA, 1952; NACA-RM-L52C25.
  22. Alam, M.; Kara, K.; Alexander, A. Reduction of spanwise flow over a swept wing using an air curtain produced by rectangular slot. In AIAA AVIATION 2022 Forum; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2022; p. 3296. [Google Scholar]
  23. Rudolph, P.K. High-Lift Systems on Commercial Subsonic Airliners; NASA: Washington, DC, USA, 1996; NASA-CR-4746. [Google Scholar]
  24. Malik, M.; Lin, R. Transition prediction on the slat of a high-lift system. J. Aircr. 2004, 41, 1384–1392. [Google Scholar] [CrossRef]
  25. Renukumar, B.; Bramkamp, F.; Hesse, M.; Ballmann, J. Effect of flap and slat riggings on 2-D high-lift aerodynamics. J. Aircr. 2006, 43, 1259–1271. [Google Scholar]
  26. Steiner, J.; Viré, A.; Benetti, F.; Timmer, N.; Dwight, R. Parametric slat design study for thick-base airfoils at high Reynolds numbers. Wind Energ. Sci. 2020, 5, 1075–1095. [Google Scholar] [CrossRef]
  27. Wang, H.; Jiang, X.; Chao, Y.; Li, Q.; Li, M.; Zheng, W.; Chen, T. Effects of leading edge slat on flow separation and aerodynamic performance of wind turbine. Energy 2019, 182, 988–998. [Google Scholar] [CrossRef]
  28. Li, W.; Guo, Y.; Liu, W. On the mechanism of acoustic resonances from a leading-edge slat. Aerosp. Sci. Technol. 2021, 113, 106711. [Google Scholar] [CrossRef]
  29. Verstraete, D.; Palmer, J.L.; Hornung, M. Preliminary Sizing Correlations for Fixed-Wing Unmanned Aerial Vehicle Characteristics. J. Aircr. 2018, 55, 715–726. [Google Scholar] [CrossRef]
  30. Panagiotou, P.; Antoniou, S.; Yakinthos, K. Cant angle morphing winglets investigation for the enhancement of the aerodynamic, stability and performance characteristics of a tactical Blended-Wing-Body UAV. Aerosp. Sci. Technol. 2022, 123, 107467. [Google Scholar] [CrossRef]
  31. Olson, L.E.; McGowan, P.R.; Guest, C.J. Leading-Edge Slat Optimization for Maximum Airfoil Lift; NASA: Washington, DC, USA, 1979; No. A-7753. [Google Scholar]
  32. Weick, F.E.; Sanders, R. Wind-Tunnel Tests on Combinations of a Wing with Fixed Auxiliary Aerofoils Having Various Chords and Profiles; NASA: Washington, DC, USA, 1934; NACA-TR-472. [Google Scholar]
  33. Roy, R.K. Primer on the Taguchi Method, 2nd ed.; Society of Manufacturing Engineers: Deadborn, MI, USA, 2010; ISBN 9780872638648. [Google Scholar]
  34. Phadke, M.S. Quality Engineering Using Robus, 1st ed.; Prentice Hall: Englewood Cliffs, NJ, USA, 1989; ISBN 0-13-745167-9. [Google Scholar]
  35. Austin, R. Unmanned Aircraft Systems, UAVS Design, Development and Deployment; Wiley: Chichester, UK, 2010. [Google Scholar]
  36. Spalart, P.R.; Allmaras, S.R.; Reno, J. One-Equatlon Turbulence Model for Aerodynamic Flows. In Proceedings of the 30th Aerospace Sciences Meeting Exhibit, Reno, NV, USA, 6–9 January 1992; p. 23. [Google Scholar]
  37. Wilcox, D.C. Turbulence Modeling for CFD, 2nd ed.; DCW Industries: La Canada, CA, USA, 1998. [Google Scholar]
  38. Rumsey, C.L.; Spalart, P.R. Effective inflow conditions for turbulence models in aerodynamic calculations. AIAA J. 2009, 45, 2544–2553. [Google Scholar]
  39. Nelson, R.C. Flight Stability and Automatic Control; McGraw-Hill: New York, NY, USA, 1997; ISBN 0070462739. [Google Scholar]
  40. Vankanti, V.K.; Ganta, V. Optimization of Process Parameters in Drilling of GFRP Composite Using Taguchi Method. J. Mater. Res. Technol. 2014, 3, 35–41. [Google Scholar] [CrossRef]
Figure 1. External layout of the tactical BWB UAV baseline platform in millimeters (mm).
Figure 1. External layout of the tactical BWB UAV baseline platform in millimeters (mm).
Aerospace 10 00720 g001
Figure 2. Slat design geometric parameters definition [11].
Figure 2. Slat design geometric parameters definition [11].
Aerospace 10 00720 g002
Figure 3. Indicative comparison between two different configurations of the L9 array.
Figure 3. Indicative comparison between two different configurations of the L9 array.
Aerospace 10 00720 g003
Figure 4. Computational mesh (a) around the UAV and (b) slat detailed view.
Figure 4. Computational mesh (a) around the UAV and (b) slat detailed view.
Aerospace 10 00720 g004
Figure 5. Detailed view of the mesh (a) over the slat element and (b) the boundary layer.
Figure 5. Detailed view of the mesh (a) over the slat element and (b) the boundary layer.
Aerospace 10 00720 g005
Figure 6. Results of the grid-dependency study for the Lift (top) and Drag (bottom) coefficients.
Figure 6. Results of the grid-dependency study for the Lift (top) and Drag (bottom) coefficients.
Aerospace 10 00720 g006
Figure 7. Lift coefficient as a function of angle of attack for configurations 1–4 compared with baseline (top) and configurations 5–9 compared with baseline (bottom).
Figure 7. Lift coefficient as a function of angle of attack for configurations 1–4 compared with baseline (top) and configurations 5–9 compared with baseline (bottom).
Aerospace 10 00720 g007
Figure 8. Drag polars of the examined configurations.
Figure 8. Drag polars of the examined configurations.
Aerospace 10 00720 g008
Figure 9. Lift Coefficient vs. Pitching Moment Coefficient for all of the examined configurations.
Figure 9. Lift Coefficient vs. Pitching Moment Coefficient for all of the examined configurations.
Aerospace 10 00720 g009
Figure 10. Lift Coefficient vs. Pitching Moment Coefficient for configurations 2 (top) and 8 (bottom).
Figure 10. Lift Coefficient vs. Pitching Moment Coefficient for configurations 2 (top) and 8 (bottom).
Aerospace 10 00720 g010
Figure 11. Effect of the design parameters on the pitch break angle.
Figure 11. Effect of the design parameters on the pitch break angle.
Aerospace 10 00720 g011
Figure 12. Effect of the design parameters on Vstall.
Figure 12. Effect of the design parameters on Vstall.
Aerospace 10 00720 g012
Figure 13. Effect of the design parameters on ΔxNP.
Figure 13. Effect of the design parameters on ΔxNP.
Aerospace 10 00720 g013
Figure 14. Contribution of the design parameters on the performance criteria.
Figure 14. Contribution of the design parameters on the performance criteria.
Aerospace 10 00720 g014
Table 1. Baseline BWB UAV platform specifications.
Table 1. Baseline BWB UAV platform specifications.
MTOW260[kg]
Cruise speed160[km/h]
Take-off speed108[km/h]
Reference area (S)4.9[m2]
Quarter-chord sweep angle (Λc/4)33[deg]
Aspect Ratio (AR)8[-]
Mean aerodynamic chord1.005[m]
Table 2. Indicative comparison between a full factorial and a Taguchi experimental design.
Table 2. Indicative comparison between a full factorial and a Taguchi experimental design.
FactorsLevelsFull Factorial DesignTaguchi Design
224 (22)4
328 (23)4
4216 (24)8
72128 (27)8
15232.768 (215)16
4381 (34)9
Table 3. Design parameters and their levels under investigation.
Table 3. Design parameters and their levels under investigation.
ParametersLevels
123
A, Slat Semispan [% b/2]456075
B, δs [degrees]102030
C, Δx [% c ¯ ]135
D, Δy [% c ¯ ]024
Table 4. L9 Orthogonal Array.
Table 4. L9 Orthogonal Array.
ConfigurationABCD
1451020
2452042
3453064
4601044
5602060
6603022
7751062
8752024
9753040
Table 5. Boundary conditions for CFD modeling.
Table 5. Boundary conditions for CFD modeling.
ParameterTake-Off SetUnits
Freestream velocity100[km/h]
27.778[m/s]
Altitude0[m]
Ambient pressure101.3[kPa]
Ambient temperature15[°C]
Eddy Viscosity Ratio0.21[-]
Fractional Intensity0.01[-]
Table 6. Performance criteria and calculated SNR.
Table 6. Performance criteria and calculated SNR.
ConfigurationPitch Break AngleSNRVstallSNRΔxNPSNR
[°][dB][m/s][dB][m][dB]
11020.0029.73−29.460.0234.42
21020.0029.94−29.530.0329.90
31221.5829.33−29.350.0822.27
41221.5827.73−28.860.0428.40
5818.0632.66−30.280.0329.12
61221.5829.48−29.390.0921.11
7818.0632.17−30.150.0428.40
8818.0632.84−30.330.1020.35
9412.0443.67−32.800.0821.94
Table 7. Response table for pitch break angle.
Table 7. Response table for pitch break angle.
LevelA (bS/2)Β (δS)C (Δx)D (Δy)
120.5319.8819.8816.70
220.4118.7117.8719.88
316.0518.4019.2420.41
Deltamax4.471.482.013.71
Rank1432
Contribution ratio [%]38.3412.6817.2031.79
Table 8. Response table for Vstall.
Table 8. Response table for Vstall.
LevelA (bS/2)Β (δS)C (Δx)D (Δy)
1−29.45−29.49−29.73−30.85
2−29.51−30.04−30.40−29.69
3−31.09−30.51−29.93−29.51
Deltamax1.651.020.671.34
Rank1342
Contribution ratio [%]35.2421.8714.2928.61
Table 9. Response table for ΔxNP.
Table 9. Response table for ΔxNP.
LevelA (bS/2)Β (δS)C (Δx)D (Δy)
1−29.45−29.49−29.73−30.85
2−29.51−30.04−30.40−29.69
3−31.09−30.51−29.93−29.51
Deltamax1.651.020.671.34
Rank1342
Contribution ratio [%]35.2421.8714.2928.61
Table 10. Optimum design parameter combinations for the optimization of the performance criteria.
Table 10. Optimum design parameter combinations for the optimization of the performance criteria.
Design Parameters
A (Span)B (δs)C (Δx)D (Δy)
Pitch break angle
Optimum level1113
Optimum value45015
Contribution [%]38.3412.6817.2031.79
Vstall
Optimum level1113
Optimum value45015
Contribution [%]35.2421.8714.2928.61
ΔxNP
Optimum level1121
Optimum value45031
Contribution [%]26.2242.767.1823.84
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Antoniou, S.; Kapsalis, S.; Panagiotou, P.; Yakinthos, K. Parametric Investigation of Leading-Edge Slats on a Blended-Wing-Body UAV Using the Taguchi Method. Aerospace 2023, 10, 720. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace10080720

AMA Style

Antoniou S, Kapsalis S, Panagiotou P, Yakinthos K. Parametric Investigation of Leading-Edge Slats on a Blended-Wing-Body UAV Using the Taguchi Method. Aerospace. 2023; 10(8):720. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace10080720

Chicago/Turabian Style

Antoniou, Spyridon, Stavros Kapsalis, Pericles Panagiotou, and Kyros Yakinthos. 2023. "Parametric Investigation of Leading-Edge Slats on a Blended-Wing-Body UAV Using the Taguchi Method" Aerospace 10, no. 8: 720. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace10080720

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