The purpose of this section is to provide a foundational understanding as to how the flow behaves under dynamic TECS and LECS kinematics in order to provide the context for the detailed measurements and analysis of other test cases in the subsequent sections. An intermediate deflection rate was chosen here as a representative case to investigate the underlying flow mechanics responsible for the transient response.
3.1. Surface Pressure and Flow Behavior during a Dynamic Trailing-Edge Control Surface Deflection
Figure 6 presents the
response of a dynamically deflected TECS from
to
at an intermediate deflection rate of
. Upon initiation of the control surface deflection (a), the
responds immediately, reaching a maximum value of 1.4 at approximately 80% of the control surface deflection (b). This indicates that the flow adheres closely to the newly acquired camber line, even at significant deflection angles, throughout the control surface transient. Once the deflection ends, the flow separates, and the transient
gradually relaxes to a steady-state value, as depicted in
Figure 6 (d–f). These findings align with PIV and smoke flow visualization results from previous studies by Mancini and Rennie [
22,
23], where dynamically actuated TECS exhibited the formation of a large Trailing-Edge Vortex (TEV) upon reaching the end of the deflection, followed by the shedding of several smaller TEVs.
Figure 7 shows the chord-wise pressure distribution at various stages of the control surface deflection measured using the techniques described in
Section 2.3. Before the motion starts, there is a typical
distribution where the top and bottom surface pressures are nearly identical, and the stagnation pressure is at unity; see
Figure 7a. As the deflection begins (see
Figure 7b), the
on the top surface decreases and increases on the bottom surface with the development of a favorable static pressure gradient on the top surface during the control surface deflection. Thus, the flow remains momentarily attached despite the high deflection angles; see
Figure 7c. This was also seen in the PIV and smoke flow visualization experiments by Mancini and Rennie [
23,
29]. The largest pressure difference between the top and bottom surface is created just past halfway of the deflection history. This correlates well with the largest
peak noticed in
Figure 6.
Figure 7b,c also show that the pressure decreases (
) on the upper surface, approaching the hinge location, where the attached flow experiences a sharp turn down the control surface (acceleration due to control surface geometry). Following the end of the control surface motion, the
distribution on the top surface gradually flattens just aft of the leading edge; see
Figure 7d–f. This suggests that the attached flow gradually separates aft of the hinge and relaxes to a steady-state condition.
Figure 8 shows the pressure distribution as the TECS deflection is reversed (i.e., returning from
to
). The
distribution on the top surface shows a region of constant pressure on the front of the airfoil, indicating a region of separation around the leading edge present at the deflected position; see
Figure 8a. As the deflection progresses, this region of separated flow reforms to a single suction peak on the leading edge; see
Figure 8d. Towards the end of the deflection, the suction peak reduces in magnitude, and the pressure on the top surface and bottom surface again becomes nearly identical; see
Figure 8f. The flow on the bottom surface conforms with the changing camber during the deflection and remains attached at the undeflected position. The transition from a separated flow to the attached flow was thus relatively smoother with fewer and smaller vortical structures [
23].
Theodorsen’s unsteady aerodynamic theory suggests that the transient
peak produced within a relatively short convective time during TECS deflection arrives from the deflection rate-dependent lift sources (virtual camber and rotation-induced normal acceleration) and added mass [
30]. Thus, the relationship between the magnitude of these
peaks to the deflection rate will be investigated in the latter sections.
3.2. Surface Pressure and Flow Behavior during a Dynamic Leading-Edge Control Surface Deflection
Figure 9 shows the time-varying
when deflecting an LECS from
to
at
. An increase in lift is noticed immediately after the LECS begins to deflect (
Figure 9), much like in the case of TECS. Snapshots from smoke flow visualization experiments, presented in
Figure 10, show that once the control surface begins to move, the flow conforms to the changing geometry and remains attached. Once the control surface reaches the full deflection, the shear layer rolls up on the upper surface of the airfoil to create an LEV, which grows in size as it convects across the airfoil. This presence of low-core pressure on the top surface can be related to the initial lift peak in
Figure 9 (a–c). Concurrently there are also vortical structures that begin to form around the hinge on the bottom surface, which creates a negative lift component. Note that a region of re-circulation was found on the airfoil at higher deflection angles during the time-averaged analysis of statically deflected cases, presented in the static analysis paper [
19]. This occurrence of bottom-surface suction can be correlated to the lift being generated in the opposite direction in
Figure 9 (d–f).
While the initial LEV on the top surface travels to approximately 30% of the chord, there is additional shedding of smaller vortices from the leading edge. The initial LEV grows in size as it convects towards the trailing edge, forming a larger vortical structure on the upper surface of the airfoil; see
Figure 10. In comparison to the vortical structures seen on the bottom surface, the occurrence of vortex mixing can be related to the
to increase in the positive direction; see
Figure 9 (f–h). This relatively large region of low pressure on the upper surface of the airfoil is only momentary, hence the departure of LEV from the airfoil at
, and the bottom-surface suction contributes to opposing lift generation (dip in the lift for the second time shown in
Figure 9 (h–j)). After
, the shedding on the LEVs on the upper surface does remain relatively consistent, and
gradually decreases and relaxes to a steady-state condition where the flow on the upper surface separates entirely.
after
thus remains settled to a steady-state
. Towards the end of the transient effect, the smaller LEVs that are shed from the leading edge on the top surface were found to be much further away from the airfoil surface. Though the quantitative contributions to the overall lift from these vortices are not evaluated in this research, it can be said with confidence that the vortices further away from the airfoil have relatively less impact on the overall
than the initial LEV. Hence, the formation and convection characteristics of these vortices play a significant role in lift production over the airfoil.
The initial increase in
, seen in
Figure 9, reaches a maximum
of 0.55. This occurs at
. During this stage, the shear layer on the upper surface of the airfoil rolls up to form the initial LEV. As the control surface has just reached full deflection, the region of re-circulation and thus the bottom surface suction does not appear to have occurred. Hence, if the control surface could be returned to the un-deflected position after deflection, this
peak could be potentially used without the deleterious consequences (though relatively smaller when compared to the case of TECS) of opposing lift arriving from the loss of initial LEV on the top- and the bottom-surface suction (which forms after the control surface reaches the end of the deflection).
The
distribution during LECS deflection at the same rate of
is shown in
Figure 11. At an undeflected position, the
distributions over the top and bottom surfaces of the airfoil are identical, as expected nominally; see
Figure 11a. As the deflection begins, the suction peak on the top surface gradually increases to the end of the deflection. The shear layer beings to roll up and starts forming an LEV, shown by the peak in
on the top surface in
Figure 11b. This correlates well with the formation of the LEV seen from
Figure 10. After the maximum
peak,
then reduced to approximately 0.25.
Figure 11d–f shows that this is due to a significant reduction in suction pressure. This is due to the LEVs convecting far away from the top surface of the airfoil. There is a region of re-circulation that starts forming around the hinge on the bottom surface, shown by the peak in
on the bottom surface of the hinge (
Figure 11d–f), indicated by small changes to the bottom pressure. However, the largest contribution to lift remains dependent on the upper surface. There is then a small positive trend of
seen between
= 2.7 and 3.2. This is again due to an increase in the suction pressure on the top surface. A suction peak forms at the hinge on the bottom surface, creating a favorable pressure gradient on the bottom surface. Hence, the flow remains attached on the bottom surface throughout most of the transient motion.
When the control surface is returned from its deflected position, i.e.,
to
, the flow transitions from being completely separated to attached, and it is observed that the transition consists of formations of relatively smaller vortical structures, much like in the case of TECS. As the deflection begins, the fully separated boundary layer transitions to discrete vortex shedding, after which the flow attaches from the leading edge to the trailing edge.
Figure 12 shows surface pressure data during and after the deflection. At
, the flow on the top surface is completely separated, shown by relatively flat
on
Figure 12a.
on the bottom surface has a negative suction
peak just aft of the leading edge and a suction peak at the hinge, creating a region of favorable pressure gradient on the control surface. Aft of the hinge, the pressure is returning to free-stream pressure. As the deflection initiates, the inflated peak on the bottom surface aft of the leading edge and the suction peak seen at the hinge both gradually reduce in size. The flow settles to a steady condition much quicker (by three convective times) in this case, relative to the previous deflection case from
to
.
From the investigation of the general flow mechanics associated with dynamic deflections of LECS and TECS, two key conclusions can be drawn. The first is that the lift responds immediately with control surface deflection irrespective of the flow condition, i.e., separated or attached. For both LECS and TECS, it was observed that responded immediately when deflecting from attached to separated flow case and vice versa. The second conclusion is that the formation of vortical structures plays a vital role in the overall duration of the lift response. It was observed that the duration of the unsteady effect (i.e., time for the flow to settle to steady-state) was relatively large for the case of LECS from deflecting from to . This was not the case during deflection from to or in the case of TECS. As the only significant difference between the three cases is the presence of LEVs, the duration of the unsteady effect can thus be linked to the formation of LEVs, which provides a relatively large circulatory lift contribution. Thus, it can be concluded that the response of LECS is a trade-off of positive lift produced on the airfoil by LEVs on the top surface and negative lift generated by the bottom-surface suction.
3.3. Effect of Varying Actuation Rates
It is now apparent that part of the transient lift arises from sources that are dependent on the actuation rate of the control surface. The effects of actuation rates on the lift responses of TECS and LECS are studied in this section for the deflecting and returning conditions, i.e., to (attached to separated) and to (separated to attached). The actuation rates selected for this study represent a range of flow conditions around the airfoil, varying from near-steady conditions () to unsteady flow conditions ().
Figure 13 (left) shows
against convective time of dynamic TECS deflections at five different actuation rates. The magnitude of the
peak was found to increase significantly with actuation rates. Nearly 1.5 times more lift is produced during this transient period than during similar static deflections, especially at
= 0.71 and 0.54. This is expected to be from deflection rate-dependent lift sources: added mass (arising from initially accelerating the fluid during the motion), virtual camber (variation in normal perturbation velocity), and rotationally induced normal acceleration (from each instantaneous angular deflection). The combination of these actuation rate-dependent lift sources being proportional to acceleration leads to the abrupt rise in the aerodynamic force during the control surface deflection. As expected, during lower actuation rates, the
peaks are less obvious and are closer to static values.
At halfway through the control surface deflection (where the angular velocity is maximum and acceleration is approaching zero), deflection rate-dependent terms are evidently responsible for a large portion of the lift history. The negative trend in control surface acceleration, as it slows down to conclude its motion, likewise correlates to a local lift trough. Around
, the
values from all TECS cases converge to steady-state values. These results correlated well with the findings from a similar analysis of lift response during a rapidly actuated trailing-edge control surface by [
31,
32].
Figure 14 shows the pressure distribution over the airfoil during various stages of the TECS deflection (top to bottom) and various deflection rates (left to right). For the case of
(left column), it is observed that a favorable pressure gradient exists on the top section during the deflection, which keeps the flow attached to the control surface during the deflection. As the control surface approaches the end of its deflection (subplot d), the area between the top and bottom surface reduces. Post-deflection, the flow on the top surface begins separation (shown by the flat pressure distribution on subplot g) and eventually relates to the steady-state (subplot j). At the higher actuation rate, the area between the top and bottom surface is greater, resolving to a larger transient
.
Figure 13 (right), presented at the start of this section, also shows five deflection cases at various actuation rates for the LECS configuration. The
peaks for all cases vary far less than with TECS. The difference in
peak between the fastest rate of
and the slowest rate
is approximately 0.15. However, when compared to the steady values, the
peaks are generally around three to four times greater. The PIV vorticity study of dynamically actuated LECS at the same rates (
Figure 15) suggests that dynamic actuation of LECS leads to the formation of LEVs on the top surface, which are known to improve the lift on the airfoil while they remain over the airfoil. It is thus evident that, with LECS, even at slower actuation rates, it is possible to achieve relatively large transient control forces. For practical applications, this could mean that once gust is detected, LECS is deflected at a realizable rate and then returned back to the undeflected position at an even more leisurely rate. Unlike the slower actuation rates, multiple peaks are observed during faster actuation (≥
); see
Figure 13 (right). As discussed previously, with
, the formation and convection of the LEVs are responsible for the fluctuation in the
response. With higher actuation rates, large dominant LEVs are formed, which have different convection characteristics (they convect relatively faster) than smaller LEVs. Hence, the
fluctuations for
and
show a relatively significant
fluctuation during the transient stage of the
response.
Figure 15 displays a vorticity flow field at various stages of the control surface deflection. In all cases, the leading-edge deflection activates the formation of an LEV. LEVs are seen to grow in size as they convect downstream with the flow. It was found that fast actuation rates produced relatively large dominating LEVs, whereas slow actuation rates were found to be largely separated and featured small LEV shedding (
Figure 10). During a fast actuation rate,
, a classical LEV is formed, dominating the flow, and it is hypothesized that the pitch component of a shear layer aids LEV formation. At slow deflection rates, flow generally separates, and the dominance of LEV is not seen. This presents an interesting question as to whether the development of upper surface flow and LEVs are also dependent on the local angle of attack at the leading edge. The study by [
31], where TECS was dynamically actuated from
to
, addressed this particular question and found the transient lift/pitching moment response to be similar, whether the fore element was at an incidence of
or
. Investigations of the effect of actuation rates on the LEV size were also previously carried out by [
33,
34,
35,
36,
37,
38]. In these investigations, it was found that increasing the pitch rate delayed the formation of LEVs on the upper surface and made the LEVs more compact and stronger, thus producing a much larger transient
.
Experiments were performed in the same manner as actuating oppositely.
Figure 16 shows the deflection of TECS (left) and LECS (right) returning from a deflected position, 40
, to an undeflected position, 0
. As the deflection rate increases, there is a progression toward a larger amplitude of transient lift increment; this holds true whether the initial flow condition is attached or separated. Recall that the initial condition is not attached flow, yet there is still an instantaneous response. However, the unsteady lift response here differs in magnitude to the transient peak when compared to the lift response going from the attached flow case. Much like the previous case, deflections from separated to attached flow cases are similarly driven by rate-dependent forces, which considerably surpass the contribution from quasi-steady circulatory lift and moment and result in large lift peaks during the motion transient.
Figure 17 shows the PIV vorticity fields for various LECS deflection cases.
Figure 18 summarizes the
response for TECS and LECS when deflecting from attached flow case to separated flow case (red) and vice versa (blue). The green line on the figure shows the difference in lift response between steady state and transition from separated to attached flow case. Classical aerodynamic theory using methods of linear superposition would predict a similar result between control surface deflections from
to
and
to
; however, we see here that this is not the case. It is hypothesized that the differences (green line in
Figure 18) arise from the fact that the case of control surface deflection from
to
begins its motion from a fully developed steady state, whereas the case of control surface deflection from
to
experiences a brief period of “inviscid effectiveness,” where it reaches a maximum steady circulatory lift contribution. This is evidenced by the large overshoot of the green curve at early times (
). When approaching the steady-state condition, the green line converges to the steady value at around
6 for both TECS and LECS; see
Figure 18. Thus, it can be concluded that when the end state is attached flow, relaxation to the steady state is relatively faster.
It is now apparent that the actuation rate of control surfaces or reduced frequency does have a vital influence on the lift response of the airfoil for TECS configuration. The lift response with LECS was found to be almost independent of the motion rate itself. The analyses of the lift responses of LECS and TECS provided in the preceding sections are for a case where the airfoil is aligned to the flow direction, i.e., airfoil angle of attack of zero. This was done so that the focus of the investigation remained close to the immediate unsteady effects arising from varying control surface deflection rates and flow speeds. For TECS, the end condition has zero lift if and only if the control surface comes to rest at zero airfoil angle of attack with respect to the corresponding zero-lift line. When fixed-wing drones are flying in turbulence, the effective angle of attack is constantly changing; thus, TECS must be returned back to a neutral position immediately. On the contrary, LECS produced of nearly zero with the control surface at a positive angle. Thus, LECS could be used more like a fine-tuning device for gust response. However, to validate this hypothesis, further investigations are required at various airfoil angles of attack and deflection angles.
3.4. Examination of and on Individual Airfoil Elements
A comprehensive understanding of the fundamental forces acting on the control surface and stationary element of the wing is required to capture these forces theoretically. Thus, in the following section, the distribution over the two segments of the airfoil is investigated at variations in control surface deflection rates at a fixed .
Figure 19 shows the breakdown of
for the fixed and moving elements at four different actuation rates. It is observed that slow to moderate control surface deflections do not have much of an effect on the front element. However, at higher deflection rates, unsteady peaks are also seen on the stationary element. This aligns well with the findings from Mancini’s research, which featured a 50% TECS [
23].
Figure 20 shows the pitching moment about the quarter chord for TECS (left) and LECS (right) at various deflection rates when deflecting from
to
. It is observed that (as expected) the trends of
are similar to
in variations in actuation rates in the case of TECS. This was also the conclusion drawn by a similar previous study [
23]; faster actuation rates correlated to larger transient
. As the
in the stationary element is relatively closer to
than the
of the control surface, it is expected to see a negative
for the steady state. During the motion transient, the relative magnitude of
between the two airfoil elements changes, with the control surface accounting for almost 50% of
, which creates negative spikes in
.
In the case of LECS, the initial lift spike can be attributed almost exclusively to the LECS at sufficiently high actuation rates (>0.54); see
Figure 21. At lower actuation rates, there exist significant contributions to the transient lift by the stationary rear element as well. Contrary to the case of TECS, the stationary element has a significant contribution to the steady-state lift as well. Thus, the resulting global lift peak can be attributed nearly equally to both the active LECS and the rear stationary element contributions. Furthermore, in the case of LECS, there is no clear relationship between
and deflection rates, much like seen in the case of
. Irrespective of the deflection rate, there is initially a positive
, which then reverses and eventually relaxes up to a positive
just above zero.
From the preceding sections, a sound understanding of the dynamics of LECS and TECS has been established through experimental investigations. To translate these into control algorithms for small fixed-wing drones with LECS, TECS, or a combination of both, the lift responses of LECS and TECS are compared against existing and modified airfoil theories in order to potentially arrive at simple low-order solutions that can capture the complex dynamics within reasonable accuracy.
3.5. Analytical Solutions
The existing aerodynamic solution by Theodorsen is presented in Equation (
1) [
30], for the case of TECS deflection. The coefficients
depend on airfoil and control surface geometry and
is Theodorsen’s function.
Theodorsen’s potential flow solution predicts that an airfoil with a deflecting control surface hinged on the trailing edge produces unsteady non-circulatory forces on the front element, despite it remaining stationary throughout the motion. This is due to the distribution of sources and sinks used in the model and placed on the wing to satisfy the no-through-flow boundary condition and thus affects the velocity field everywhere in the flow. Physically, the assumption of incompressibility mandates that any local disturbance to the flow (e.g., physically deflecting a control surface) causes a pressure disturbance everywhere in the flow instantaneously. Theodorsen’s model assumes attached, inviscid flow at all times. These two assumptions might be valid during the control surface motion, but not before and after the deflection where the flow is separated.
Figure 22 displays a comparison between the experimentally acquired
and Theodorsen’s solution. It is evident that Theodorsen’s solution (shown by the solid black line) does not accurately capture the experimentally acquired data (solid red line). The magnitude of the
peak is underestimated for the case of
and overestimated for
and
. Additionally, the comparison reveals a significant overprediction by Theodorsen’s solution in the relaxation period of the lift response to the steady state. To gain further insight into this disparity, the individual lift constituents from Theodorsen’s solution (shown by dashed lines) are also plotted against convective time in
Figure 22. It is evident that the lift force contributions from the steady component are significantly high, as it fails to capture flow separation at the end of the control surface deflection. Likewise, the exponential transition to steady state from Theodorsen’s function, C(k), is not accurately captured.
Figure 23 shows the comparison between experimentally acquired
and the potential flow model derived by Jaworski (Equation (
2)). It is observed that the Jaworski model predicts the magnitude and occurrence of the
peak reasonably well for the fastest actuation case but fails to do so for the two slower rates. The fluctuation of
peaks observed after the initial
peak is not captured by the potential flow solution. As discussed earlier, the formation of LEVs on the top surface and formation of the recirculation region on the bottom surface play a significant role in contributing to the unsteady transient lift response of a dynamically actuated LECS. Even well after the control surface comes to the end of the deflection, the LEV formed during the motion convects across the chord of the airfoil, aiding in the generation of lift. Although the unpacking of individual force contributions is not possible with solely surface pressure measurements, it is clear that rate-dependent forces dominate the overall lift production. Hence, this complex balance of vortical structures shed on the top and bottom surface cannot be simply represented with potential flow solutions.
A modified approach (Equations (3) and (4)), adopted from [
23], was also used to predict the lift response, shown in green in
Figure 24. Across the three different actuation rates shown, the lift response prediction is improved. The magnitude of the lift responses and relaxation to the steady state are in good agreement with the experimental data. The modified approach treats the two segments of the airfoil, TECS, and the stationary fore elements separately, assuming that there is no coupling effect between the two. This is not the case in reality, as we have previously observed that deflection of the TECS causes changes in the pressure field on the stationary segment of the airfoil, and there is lift generated by the stationary element. It has been shown that the dynamic deflection of the control surface imposes transient lift generation over the stationary element in the preceding sections. From the experimental results, the magnitude of the lift generated by TECS was found to be relatively lower compared to the lift generated by the stationary element. It is observed there that with Equation (
4), which neglects the lift produced by the stationary element, it shows good agreement with experimental results. Coincidentally, ignoring the stationary element and treating TECS essentially as a flat-plate airfoil pivoting about its leading edge made the prediction of the overall lift response “better”. This was an interesting finding and was also found by [
23]. Hence, before any conclusions could be drawn regarding validation of this model, the approach would need to be tested across a much larger test matrix featuring a range of control surface sizes.
For the modified model, the constants , and , responsible for relating the geometrical kinematics of the control surface deflections, were removed, as the total lift is now a summation of lift from a rapidly actuating flat-plate airfoil and a stationary element (not producing any lift at the airfoil angle of attack). Thus, the chord of the airfoil is essentially reduced, so the model correlated well with the experimental data.
Figure 25 shows a comparison of
from experimental data, Jaworski’s solution, and a modified model. It is clear that the modified approach cannot be applied for the case of LECS to predict transient
. The modified approach, where the control surface and the stationary element are treated as two separate elements, ignoring all coupling effects between them, shows significant disparity with the experimental results. This is unlike the case of TECS, where the modified approach improved the accuracy of estimating the lift response over Theodorsen’s solution.
As found in the preceding sections, some lift contributions arise from the stationary element. This is because, as the vortex convects across the chord of the airfoil, it is imposing a time-varying low-pressure core, which is generating lift. Hence, there is a lift contribution arising from both the body and the control surface. With the modified approach, there is no lift contribution arising from the stationary element, and all the lift is a function of LECS alone. It overpredicts, because there is nothing to attenuate the lift response. Potential flow does not account for separation, and there is nothing stopping the generation of the steady lift well past the separation point, and the flow relaxes to a steady value.
characteristics of a dynamically deflected LECS were found to be significantly dependent on the convections of vortical structures on the top or bottom surface of the airfoil. These phenomena cannot be captured by potential flow assumptions. The development and shedding of these vortices were found to be dependent on the control surface deflection rate. Thus, rate-dependent lift sources would need to be added into the theoretical models to accurately capture the lift response of dynamically actuated LECS. This can arise from further experiments specifically focused on the dynamic characterization of the development and shedding of vortical structures.