4.3.1. Sensitivity Analysis of Deflection Angles
The Krueger flap deflection region is considered first. “Theta” is defined as the deflection angle of the Krueger flap. The straight line where the trailing edge of the upper surface of the Krueger flap is located is the axis of rotation. The clockwise direction shown in the upper right corner of
Figure 13 is positive. The initial Krueger angle Theta of the high lift configuration is 30°.
The deflection direction of the Krueger flap is then considered. The purpose of this paper is to improve the stall state of the aileron region in the high lift configuration at large angles of attack. The degree of airflow upwash on the leading edge of the wing will be exacerbated with the local adverse pressure gradient increases if the Krueger flap angle is reduced. If the Krueger flap angle is reduced, the local adverse pressure gradient will increase with the degree of airflow upwash on the leading edge of the wing being exacerbated. Eventually, the stall in the aileron region becomes more serious. Therefore, the configurations of Theta 40°, 45°, and 50° are considered sensitivity analysis configurations. However, the definition of Theta is not enough for sensitivity analysis. The Krueger flap area where Theta changes must be initialized, so the spanwise length of the Krueger flaps redesigned is considered. “Length” is defined as the length of the outboard Krueger flap. The outboard Krueger flap will be redesigned, but it is not linear for the flow control of the Krueger flap to the wing along the flow direction shown in
Figure 13. Due to the three-dimensional effect with the wing sweep angle, the spanwise flow in the slot is strong, which causes the streamlines through the slot to be deflected towards the wing tip instead of flowing directly towards the trailing edge of the wing in the direction of the freestream. Therefore, the length should be extended from cycle 1 to cycle 2 to ensure that the flow control measure from the Krueger flap is effective for the entire aileron zone.
The CFD numerical simulation results with Mach number 0.15 and Reynolds number 8.90 × 10
7 are shown in
Figure 14 and
Figure 15. The sensitivity of the BWB high lift configuration studied in this section to the variation in the Krueger flap angle Theta is described from four aspects: stall state, lift characteristics, moment characteristics, and pressure distribution in the aileron zone. When Theta is 40°, there is a slight spanwise flow at the trailing edge of the outboard aileron, which is improved compared to the original high lift configuration, but the spanwise flow in the inboard aileron is still serious with a tendency for separation. When Theta is 45°, the flow in the outboard aileron is already very good, and there is only a slight spanwise flow in the inboard aileron. When the Theta is 50°, there is no obvious spanwise flow between the ailerons, and the improvement is extremely obvious compared with the original high lift configuration. Compared with the original high lift configuration, the adverse pressure gradient of Theta 50° drops more greatly, which greatly inhibits the flow separation in the aileron region, but causes a larger loss of local lift according to the comparison of pressure coefficient difference between Theta 30° and Theta 50°. Therefore, with the increase in the Krueger flap angle Theta from 30° to 50°, the streamline of the aileron shown in
Figure 14 is gradually straightened. Additionally, the stall state is gradually improved 18° angle of attack. This is because the negative pressure peak in C
p of
Figure 14 gradually decreases, resulting in the weakening of flow separation with a gradual decrease in the local adverse pressure gradient. This is because by increasing Theta at the same Alpha, the direction of the upwash flow induced by the Krueger flap on the leading edge of the wing will be more tangential to the direction of the upper wing surface, which leads to a result of gentler airflow over the leading edge of the wing. The negative pressure peak is reduced with a gentle flow. The increase in the Krueger flap angle Theta is very effective in ensuring the stall state of the aileron only from the point of view of suppressing the stall at the trailing edge. It should also be noted that the stall of the outboard aileron is easier to be improved.
Figure 15a shows the sensitivity results of the lift characteristic. The “INI” in
Figure 15 represents the initial high lift configuration. It can be seen that the lift curves in different Theta have very little difference at the angle of attack from −4° to 12°, and they all increase linearly. When the angle of attack is greater than 12°, the lift coefficient curves of the high lift configuration for each Theta are significantly different. The stall angle of attack did not change (22°), and the stalled process slowed down, which corroborates mutually with the aileron streamline state in
Figure 14. However, the maximum lift coefficient C
Lmax dropped by about 0.2 from Theta 30° to Theta 50°, which corroborates mutually with the C
p in
Figure 14. In short, with the increase in Theta, the stall characteristics become better, but the lift loss is gradually unacceptable.
Figure 15b shows the sensitivity results of the pitch moment characteristic. At the Alpha from −4° to 10°, the longitudinal static stability characteristics are good since the pitch moment coefficients decrease with the increase in the Alpha at each Theta. From Theta 40° to Theta 50°, when the Alpha is greater than 10°, the pitch moment coefficient curve moves slightly upward as the Theta increases. Therefore, the change in the Theta in the local Krueger flap will slightly change the pitch moment value under the condition of a large angle of attack, while it has little effect on the angle of minimum pitch moment. When Alpha is greater than 10°, the longitudinal static stability of the configurations of sensitivity analysis rapidly deteriorates. The angles of minimum pitch moment are 12°and these are 6° less than the original high lift configuration, which will greatly limit the available angle of attack when the aircraft takes off and lands. The outboard wing is located behind the center of mass of the aircraft and provides mainly nose-down pitch moments. Changing Theta will reduce the lift of the outboard wing, thus causing the nose-down pitch moment value provided by the outboard wing to reduce. However, the nose-up moment value provided by the inboard wing and the airframe in front of the center of mass has not changed. Therefore, the total pitching moment of the aircraft increases. This is why the angle of the minimum pitch moment decreases, and it is further analyzed in the next section with the spanwise lift distribution. The angle of attack of civil aircraft is not allowed to exceed the angles of minimum pitch moment. Therefore, it is also unacceptable for the BWB that the angles of minimum pitch moment of the configurations of sensitivity analysis reduce 4°compared to the initial high lift configuration.
In short, through the above sensitivity analysis, in terms of stall characteristics, locally increasing the Krueger flap angle can reduce the adverse pressure gradient of the wing and effectively improve the stall characteristics of the aileron zone, but in terms of aerodynamic characteristics, which would cause a certain loss of lift characteristics and move the nonlinear segment of the pitching moment upward as a whole.
In addition, the angles of the minimum pitch moment of the sensitive configurations are relatively advanced. It is speculated that the length of the outboard Krueger flap length determines the position of the angles of minimum pitch moment. Therefore, the sensitivity analysis of the length will be investigated in the next section, focusing on the relationship between length and the angle of minimum pitch moment.
4.3.2. Sensitivity Analysis of Length of the Outboard Krueger Flap
This section conducts a sensitivity analysis of the length of the outboard Krueger flap. In the previous section, changing part of the Krueger flap angle resulted in a large change in the longitudinal static stability of the BWB civil aircraft. (The angles of minimum pitch moment are 12°, and these are 6° less than the original high lift configuration.) However, only changing the deflection angle of the Krueger flap did not change the angle of minimum pitch moment from Theta 40° to 50°as shown in
Figure 15b. It is worth noting that changing the Theta also inevitably changes the length of the outboard Krueger flap. The length determines the range of effect of the outboard Krueger flap after changing the outboard Krueger flap declination angle, thereby affecting the lift distribution of the entire wing. Therefore, it is speculated that the length of the outboard Krueger flap determines the angles of the minimum pitch moment.
The length of the outboard Krueger flap defined as “Length” in the previous section needs to be parameterized now. As shown in
Figure 16, the sensitivity analysis of the length was carried out by taking the lengths 30%, 45%, and 60%.
The CFD numerical simulation results with Mach number 0.15 and Reynolds number 8.90 × 10
7 are shown in
Figure 17 and
Figure 18. The sensitivity of the BWB high lift configuration studied in this section to the variation in the length of outboard Krueger flap length is described from three aspects: stall state, lift characteristics, and moment characteristics in the aileron zone. The focus is on the change in moment characteristics.
Figure 17 shows the surface limiting streamlines of the sensitive configurations, which is similar to the situation in the previous section. As the length increases, the zone of perfect streamline gradually expands from the outboard aileron to the inboard aileron, and the stall state gradually becomes better. The adverse pressure gradient of C
p drops obviously. However, unlike in the previous section, the surface limiting streamlines of the outboard aileron zone are always perfect in this section about the sensitive analysis of the length. When the length is 30%, the flow of the outboard aileron zone is perfect, but the separation of the inboard aileron zone is serious. When the length is 60%, there is no obvious strong spanwise flow on the inboard and outboard ailerons, and the surface flow is very good. Thus, the spanwise length of the stall zone is controlled by the length.
Figure 18a shows the sensitivity results of the lift characteristic. It can be seen that as the length increases gradually, the effect of the high-lift device gradually becomes worse, which is also similar to the situation in the previous section. Therefore, the length is negatively correlated with the high lift effect. The local stall state at the high angle of attack in the aileron zone of each sensitive configuration has been improved, and the overall stall process has become moderate.
Figure 18b is the result of the pitch moment of sensitivity about length. It can be seen from
Figure 18 that as the Alpha changes from −4° to 10°, the pitch moment coefficients of different lengths have the same trend as the Alpha. When the angle of attack is smaller than the minimum pitch moment angle, the static margin is mainly maintained at around 6.8%. (Aircraft of blended wing body layout have a weaker capacity for longitudinal flattening than conventional layout plane due to the shorter body. Therefore, its static margin is smaller than in a general layout plane.) Thus, the longitudinal static stability characteristics are good. However, when the Alpha is greater than 10°, the angle of the minimum pitch moment of each sensitive configuration decreases with the increase in length. The calculation results are consistent with the speculation in the previous section. As the length increases from 0% to 60%, the angle of minimum pitch moment decreases from 18° to 12°, so the length of the outboard Krueger flap has a great influence on the longitudinal static stability of the BWB. The larger the angle of the minimum pitch moment, the greater the angle of attack available to the aircraft for safe flight. The pitch moment characteristic caused by the length is further analyzed with the spanwise lift distribution.
Figure 19 is the spanwise lift distribution of Theta sensitivity in
Section 4.3.1 and length sensitivity in
Section 4.3.2. The x coordinate is dimensionless, “b” represents the spanwise length of the aircraft, “C
L_Local” represents the local lift coefficient and “C” represents the local chord. As can be seen from
Figure 19, all the curves can be divided into two parts. The part of the outboard wing where the ailerons are located has less local lift than that of the inboard wing. In the part of the inboard wing where the flaps are located, the local lift is large because the flaps are in working condition. That is why the local lift has a large jump at the junction of the aileron zone and the flap zone. The outboard wing is highlighted by the dashed line, as this is the main study area. The black line represents the final configuration. It should be noted that Theta = 50°and length = 60% are the same configuration. According to the criterion of the control variable method, the spanwise lift distribution of the initial configuration (INI: Theta = 30°, length = 0%) is hidden. The “C
L” and “C
M” in
Figure 19 represent lift and pitch moment characteristics. It can be seen that the sensitivity analysis of Theta and length have almost similar C
L change (the difference between the maximum lift coefficients of any two adjacent sensitive configurations, such as the configurations Theta = 40°and Theta = 45° or the length of the configuration = 30% and length = 45%, is about 0.05) with the respective increase in Theta and length, but the C
M change is distinct. This is because the change in the C
L is mainly related to the lift loss, while the change in C
M is not only related to the change in lift loss but also the difference in the pressure center location. The pressure center location is determined by the spanwise lift distribution. It can be seen from the spanwise lift distribution in
Figure 19 and
Figure 20 that when Theta increases, the lift loss is almost uniformly distributed over the entire outboard wing. However, when the length increases, the lift loss is concentrated in a specific zone where the length changes. Therefore, the results of the pressure center location indicate that the diverse ways of changing the spanwise lift distribution are the reason why there is an extreme variation in C
M between the length and Theta.
In short, changing the length changes the confrontation between the nose-up moment and the nose-down moment in
Figure 20. The angle of attack corresponding to the C
Lmax of the local lift and the angle of minimum pitch moment is also changed. Increasing the length makes the pressure center location change greatly. That is the reason why the pitching moment coefficient C
M is very sensitive to length.
4.3.3. The Design of Aileron Surface Flow State
According to the sensitivity analysis before, it can be found that the stall characteristic of the aileron zone is sensitive to Theta and length. The increase in Theta suppresses the stall in the aileron zone. Additionally, the increase in the length reduces the spanwise length of the flow separation in the aileron zone. However, the increase in Theta degrades the lift characteristic seriously, while the increase in the length makes the lift characteristic and pitching moment characteristic significantly worse. Therefore, it is necessary to carry out a subtle design of the Theta value and length value of the Krueger flap, which not only effectively improves the stall characteristics of the aileron region but also makes the cost of lift characteristic and pitching moment characteristic acceptable. If the stall state of the entire aileron zone is guaranteed to be good, the cost of lift characteristic and pitching moment characteristic will be very high, such as the configuration with Theta 50°and length 60% of outboard Krueger flap. Considering the constraints such as aileron flow state, lift characteristic, pitching moment characteristic, achievability of Krueger flap motion mechanism, and wing structure, the following design parameters are finally obtained: Theta = 41.5° and length = 35%. The configuration with the parameters Theta = 41.5° and length = 35% of the outboard Krueger flap is defined as Design Configuration 1.
The CFD numerical simulation results with Mach number 0.15 and Reynolds number 8.90 × 10
7 are shown in
Figure 21 and
Figure 22. The design configuration 1 studied in this section is described from three aspects: stall state, lift characteristic, and moment characteristic in the aileron zone.
At Alpha = 18° and Alpha = 22°, relatively strong spanwise flow is generated in the entire aileron zone of INI, shown in
Figure 21a,c. However, the surface limiting streamline of the outboard aileron zone of Design Configuration 1 shown in
Figure 21b,d is relatively straight, implying that there is no strong spanwise flow, and the stall state is good. The limiting streamline on the surface of the inboard aileron has a slight spanwise flow trend, and the stall characteristic of the aileron zone in Design Configuration 1 compared with the INI has been significantly improved.
At Alpha changes from −4° to 16°, the difference in the lift curves shown in
Figure 22a between INI and Design Configuration 1 is very small. When the Alpha is greater than 16°, the lift coefficient curves have a little difference between INI and Design Configuration 1. The C
Lmax is reduced by 0.06, and the stall process slows down with the stall angle of attack maintained as 22°for Design Configuration 1. From Alpha −4°to Alpha 16°shown in
Figure 22b, the change trends of the pitch moment coefficient with the angle of attack between INI and Design Configuration 1 are consistent, and the static stability characteristics are good. When the Alpha is greater than 16°, the pitching moment of Design Configuration 1 gradually increases, and the angle of the minimum pitching moment is 16°, which is 2° smaller than the INI. Therefore, Design Configuration 1 pays a high cost of the pitching moment characteristic to eliminate the spanwise flow in the aileron zone.
In summary, Design Configuration 1 can effectively improve the spanwise flow in the aileron zone, and the lift characteristics are acceptable, but the cost of the pitch moment characteristics is high, which requires a further fine design to increase the angle of minimum pitch moment for the longitudinal static stability.
Based on the above design results and the previous sensitivity analysis results, it is found that there is a contradiction between aerodynamic coefficients and stall characteristics. At the large angle of attack, the better the stall characteristic of the BWB aileron zone, the worse the lift coefficient and pitch moment coefficient are under the condition that the airfoil of this aircraft can not be changed. The contradictions need to be reconciled to solve the stall problem in the aileron zone at the expense of acceptable aerodynamic characteristics at the appropriate location. However, the above analysis shows that Design Configuration 1 obtained only by designing the parameters Theta and length cannot perfectly reconcile this contradiction (the angle of minimum pitching moment is reduced by 2°). Therefore, it is vital to find a new solution to this problem.
It is worth noting that the stall problem in the outboard aileron zone is easier to be improved than in the inboard aileron zone. When the stall state of the inboard aileron zone is relatively great, in fact, the outboard aileron zone has been improved too much. Therefore, there is redundancy in the stall improvement in the outboard aileron zone. The longitudinal static stability loss caused by Design Configuration 1 (the angle of minimum pitching moment reduced by 2°) can be corrected by appropriately reducing the flow control to the outboard wing. That has no effect on the improvement of the stall characteristic of the outboard aileron zone by Design Configuration 1. As shown in
Figure 23, it is found that the slot between the Krueger flap and the wing at the wing tip is narrower (zone of the white arrow) than the slot near the blend wing body. Flow over the upper surface of the wing is induced by the lower surface of the Krueger flap, forming an upwash flow through the slot. It is a great method to solve the stall problem by studying flow control [
28,
29,
30]. Through the CFD analysis of the slot width, within a certain range, the narrower the slot, the stronger the flow control effect of the Krueger flap. Therefore, a refined design of the width of the Krueger flap slot is considered. Width is defined as the width of the Krueger flap slot. Through the iterative design of the width, the final Krueger flap was operated as follows under the condition that the Krueger flap kinematic mechanism is guaranteed to be achievable. On the one hand, based on Design Configuration 1, the width is doubled to slightly reduce the flow control of the outboard wing. On the other hand, based on the previous step, the spanwise distribution of the width is further finely adjusted. Based on the investigable results of the arm of force in
Figure 19 and
Figure 20 and the relationship between stall state and aerodynamic characteristics, the width is further increased by 20% of its width. Finally, Design Configuration 2, shown in
Figure 23, is obtained by slightly weakening the stall characteristic of the aileron zone in the above design and compensating for the loss of aerodynamic characteristics of Design Configuration 1.
The CFD numerical simulation results with Mach number 0.15 and Reynolds number 8.90 × 10
7 are shown in
Figure 24 and
Figure 25. Design Configuration 2 investigated in this section is described from three aspects: stall state, lift characteristic, and pitch moment characteristic in the aileron zone.
The separation flow in the aileron zone of INI shown in
Figure 24 is strong at Alpha 18° and Alpha 22°, so their local stall characteristics are bad. Although there is a slight spanwise flow trend on the surface limiting streamline of the aileron zone of Design Configuration 2, the overall streamline is relatively straight, and the stall characteristic is good all the time. Therefore, increasing the width has no obvious bad effect on the flow control of the Krueger flap to the wing. In addition, the flow separation, red circle line in
Figure 24, occurs first at the blend wing body at the Alpha 22° when the surface flow on the upper aileron is good so that the aircraft can obtain an obvious stall signal (lift drop) before the aileron is in hidden danger of control in the Design Configuration 2 at a high angle of attack.
At Alpha from −4°to 16°, the difference of the lift curves shown in
Figure 25a between INI and Design Configuration 2 is very small, and they all increase linearly. When the Alpha is greater than 16°, the lift coefficient curves have a little difference between INI and Design Configuration 2. The C
Lmax is reduced by 0.03 for INI and increased by 0.03 for Design Configuration 1. The stall angle of attack did not change, it was 22°, and the stall process slowed down. The outboard Krueger flaps are less costly to suppress flow separation in the aileron zone.
From the Alpha −4°to 16°shown in
Figure 25b, the change trends of the pitch moment coefficient with the angle of attack between the INI and the Design Configuration 2 are consistent, and the static stability characteristics are good. When the Alpha is greater than 16°, the pitching moment of Design Configuration 1 gradually increases, and the angle of the minimum pitching moment is 16°. However, the pitching moment of Design Configuration 2 continues to decrease and maintain the great static stability characteristic. When the Alpha is greater than 18°, the pitching moment of both the INI and the Design Configuration 2 starts to increase, so the angle of minimum pitch moment of the Design Configuration 2 is 18°, which is the same as the INI and better than the Design Configuration 1. At the same time, the aileron streamlines of Design Configuration 2 shown in
Figure 26 after widening the Krueger slot are slightly more deflected in the spanwise direction than that of Design Configuration 1, which is the price of the compromise of the stall characteristic to the aerodynamic characteristics. The pitch moment increase in Design Configuration 2 is slightly faster than that of the INI when the Alpha is greater than 18°, but the overall longitudinal static stability characteristic is basically the same as INI. In short, the design result of increasing the width is as expected.
In summary, appropriately increasing the width achieves better compromise design results from Design Configuration 1 to Design Configuration 2. The Design Configuration 2 can also effectively improve the flow separation in the aileron, with good lift characteristic, good pitching moment characteristic, and the angle of minimum pitch moment that is the same as the INI is 18°. Design Configuration 2, based on better inheriting the great aerodynamic characteristics of the INI, paid a very small price for the lift and pitching moment and completed the improvement of the stall characteristics of the aileron zone. It should be noted that the shape parameters of the Krueger flap can also slightly affect the aerodynamic characteristics, but redesigning the shape of the Krueger flap needs to consider the constraints of the wing structure, which is much more expensive than redesigning the width.
Finally, for the design of the aileron surface flow state characteristics of the three-dimensional Krueger flap, three important adjustment parameters are summarized as shown in
Figure 27, which can provide a reference for aerodynamic design work:
- (a)
Theta, the angle of the outboard Krueger flap, affects flow separation mainly;
- (b)
Length, the length of the outboard Krueger flap, affects pitch moment tremendously;
- (c)
Width, the width of the outboard Krueger flap slot, affects pitch moment slightly.