Modification of static-wing tip vortex via a slender half-delta wing

The modification of the tip vortex generated by a rectangular NACA 0012 wing via a tip-mounted slender half-delta wing (HDW) was attempted experimentally at Re=2.81×10⁵. In addition to the increase in lift with increasing HDW deflection, compared to the baseline wing, the roll-up process of the tip vortex was also found to be significantly modified, as a result of the breakdown of the HDW vortex. The addition of the HDW also caused an increased total drag. Fortunately, the lift-induced drag was found to be reduced compared to its baseline counterpart for 0° and 5° HDW deflections. The change in the lift-induced drag also translates into a virtually unchanged profile drag, regardless of HDW deflection. In short, the largest lift-induced drag reduction achieved by the zero-deflection HDW resulted in an improved lift-to-drag ratio, at high angle-of-attack range, compared to the baseline wing.     

Computational modeling of vortex breakdown control on a delta wing

We present an effort to model the development and the control of the vortex breakdown phenomenon on a delta wing. The pair of the vortices formed on the suction side of a delta wing is the major contributor to the lift generation. As the angle of attack increases, these vortices become more robust, having high vorticity values. The critical point of a delta wing operation is the moment when these vortices, after a certain angle of attack, are detached from the wing surface and wing stall occurs. In order to delay or control the vortex breakdown mechanism, various techniques have been developed. In the present work, the technique based on the use of jet-flaps is numerically investigated with computational fluid dynamics by adopting two eddy-viscosity turbulence models. The computational results are compared with the experimental data of Shih and Ding (1996). It is shown that between the two turbulence models, the more advanced one, which adopts a non-linear constitutive expression for the Reynolds-stresses, is capable to capture the vortex breakdown location for a variety of jet exit angles. The performance assessment of the models is followed by the investigation of the effect of the jet-flap on the lift and drag coefficients.     

Investigation on vortex shedding of a swept-back wing

Effects of Reynolds number and angle of attack on the vortex shedding of a finite swept-back wing are experimentally studied. The cross-sectional profile of the wing is NACA 0012, and the sweep-back angle is 15° The time series signals detected by hot-wire in the wake region shows four distinct behaviors: laminar, subcritical, transitional, and supercritical. The derived Strouhal number curves are significantly different in these four behaviors. In addition, the statistical properties of turbulence, that is, the power spectrum density function, probability density function, correlation coefficient, Lagrangian integral time scales, and length scales are also presented in this paper.     

Near-field tip vortex behind a swept wing model

The near-field flow structure of a tip vortex behind a sweptback and tapered NACA 0015 wing was investigated and compared with a rectangular wing at the same lift force and Re=1.81×10⁵. The tangential velocity decreased with the downstream distance while increased with the airfoil incidence. The core radius was about 3% of the root chord c _r, regardless of the downstream distance and α for α<8°. The core axial velocity was always wake-like. The core Γ_c and total Γ_o circulation of the tip vortex remained nearly constant up to x/c _r=3.5 and had a Γ_c/Γ_o ratio of 0.63. The total circulation of the tip vortex accounted for only about 40% of the bound root circulation Γ_b. For a rectangular wing, the axial flow exhibited islands of wake- and jet-like velocity distributions with Γ_c/Γ_o=0.75 and Γ_o/Γ_b=0.70. For the sweptback and tapered wing tested, the inner region of the tip vortex flow exhibited a self-similar behavior for x/c _r≥1.0. The lift force computed from the spanwise circulation distributions agreed well with the force-balance data. A large difference in the lift-induced drag was, however, observed between the wake integral method and the inviscid lifting-line theory.     

Improving vortex models via optimal control theory

Low-order inviscid point vortex models have demonstrated success in capturing the qualitative behavior of aerodynamic forces resulting from unsteady lifting surface maneuvers. However, the quantitative agreement is often lacking for separated flows as a result of the ambiguity in the edge conditions in this fundamentally unsteady process. In this work, we develop a model reduction framework in which such models can be systematically improved with empirical results. We consider the low-order impulse matching vortex model in which, in its original form, Kutta conditions are applied at both edges to determine the strengths of single point vortices shed from each edge. Here, we relax the Kutta condition imposed at the plate׳s edges and instead seek the time rate of change of the vortex strengths that minimize the discrepancy between the model-predicted and high-fidelity simulation force histories, while the vortex positions adhere to the dynamics of the low-order model. A constrained minimization problem is constructed within an optimal control framework and solved by means of variational principles. The optimization approach is demonstrated on several unsteady wing maneuvers, including pitch-up and impulsive translation at a fixed angle of attack. Additionally, a stitching technique is introduced for extending the time interval over which the model is optimized.     

Effects of a trapped vortex cell on a thick wing airfoil

The effects of a trapped vortex cell (TVC) on the aerodynamic performance of a NACA0024 wing model were investigated experimentally at Re = 10⁶ and 6.67×10⁵6.67\times 10^{5}. The static pressure distributions around the model and the wake velocity profiles were measured to obtain lift and drag coefficients, for both the clean airfoil and the controlled configurations. Suction was applied in the cavity region to stabilize the trapped vortex. For comparison, a classical boundary layer suction configuration was also tested. The drag coefficient curve of the TVC-controlled airfoil showed sharp discontinuities and bifurcative behavior, generating two drag modes. A strong influence of the angle of attack, the suction rate and the Reynolds number on the drag coefficient was observed. With respect to the clean airfoil, the control led to a drag reduction only if the suction was high enough. Compared to the classical boundary layer suction configuration, the drag reduction was higher for the same amount of suction only in a specific range of incidence, i.e., α = −2° to α = 6° and only for the higher Reynolds number. For all the other conditions, the classical boundary layer suction configuration gave better drag performances. Moderate increments of lift were observed for the TVC-controlled airfoil at low incidence, while a 20% lift enhancement was observed in the stall region with respect to the baseline. However, the same lift increments were also observed for the classical boundary layer suction configuration. Pressure fluctuation measurements in the cavity region suggested a very complex interaction of several flow features. The two drag modes were characterized by typical unsteady phenomena observed in rectangular cavity flows, namely the shear layer mode and the wake mode.     

Large-eddy simulation of wing tip vortex in the near field

Large-eddy simulation with filtered-structure-function subgrid model and implicit large-eddy simulation (ILES without explicit subgrid model) using high-order accuracy and high resolution compact scheme have been performed on the tip vortex shedding from a rectangular half-wing with a NACA 0012 airfoil section and a rounded wing tip. The formation of the tip vortex and its initial development in the boundary layer and the near field wake are investigated and analysed in detail. The physics, why the tip vortex, which is originally turbulent in the boundary layer, is re-laminarised and becomes stable and laminar rapidly after shedding in the near field, is revealed by this simulation. The computation also shows the widely used second-order subgrid model is not consistent to six-order compact scheme and would degenerate the six-order LES results to second-order. Therefore, high-order schemes, grid refinement and six-order subgrid models are critical to LES approaches.     

Effects of a vectored trailing edge jet on delta wing vortex breakdown

Flow visualization was used to study the effects of a vectored trailing edge jet on the leading edge vortex breakdown of a 65° delta wing. The experimental results indicated that there is little effect of the jet on the leading edge vortex breakdown when the angle of the vectored jet is less than 10°. With the increase of the vectored angle ß, the effect of the jet on the flow becomes stronger, i.e., the jet delays the leading edge vortex breakdown in the direction of the vectored jet, and accelerates breakdown of the other leading edge vortex. Moreover, the effect of the jet control tends to be weaker with the angle of attack.     

Three-dimensional vortex formation on a heaving low-aspect-ratio wing: Computations and experiments

This paper addresses by means of high-resolution numerical simulations and experimental quantitative imaging the three-dimensional unsteady separation process induced by large-amplitude heaving oscillations of a low-aspect-ratio wing under low-Reynolds-number conditions. Computed results are found to be in good agreement with experimental flow visualizations and PIV measurements on selected cross-flow planes. The complex unsteady three-dimensional flow structure generated during dynamic stall of the low-aspect-ratio wing is elucidated. The process is characterized by the generation of a leading-edge vortex system which is pinned at the front corners of the plate and which exhibits intense transverse flow toward the wing centerline during its initial stages of development. This vortex detaches from the corners and evolves into an newly found arch-type structure. The legs of the arch vortex move along the surface toward the wing centerline and reconnect forming a ring-like structure which is shed as the next plunging cycle begins. Vortex breakdown, total collapse and reformation of the wing tip vortices are also observed at various stages of the heaving motion. At the relatively high value of reduced frequency considered, these basic flow elements of the complex three-dimensional dynamic stall process are found to persist over a range of Reynolds numbers.     

An experimental study of the unsteady vortex structures in the wake of a root-fixed flapping wing

An experimental study was conducted to characterize the evolution of the unsteady vortex structures in the wake of a root-fixed flapping wing with the wing size, stroke amplitude, and flapping frequency within the range of insect characteristics for the development of novel insect-sized nano-air-vehicles (NAVs). The experiments were conducted in a low-speed wing tunnel with a miniaturized piezoelectric wing (i.e., chord length, C = 12.7 mm) flapping at a frequency of 60 Hz (i.e., f = 60 Hz). The non-dimensional parameters of the flapping wing are chord Reynolds number of Re = 1,200, reduced frequency of k = 3.5, and non-dimensional flapping amplitude at wingtip h = A/C = 1.35. The corresponding Strouhal number (Str) is 0.33, which is well within the optimal range of 0.2 < Str < 0.4 used by flying insects and birds and swimming fishes for locomotion. A digital particle image velocimetry (PIV) system was used to achieve phased-locked and time-averaged flow field measurements to quantify the transient behavior of the wake vortices in relation to the positions of the flapping wing during the upstroke and down stroke flapping cycles. The characteristics of the wake vortex structures in the chordwise cross planes at different wingspan locations were compared quantitatively to elucidate underlying physics for a better understanding of the unsteady aerodynamics of flapping flight and to explore/optimize design paradigms for the development of novel insect-sized, flapping-wing-based NAVs.     

Reynolds number effects on leading edge vortex development on a waving wing

The waving wing experiment is a fully three-dimensional simplification of the flapping wing motion observed in nature. The spanwise velocity gradient and wing starting and stopping acceleration that exist on an insect-like flapping wing are generated by rotational motion of a finite span wing. The flow development around a waving wing at Reynolds number between 10,000 and 60,000 has been studied using flow visualization and high-speed PIV to capture the unsteady velocity field. Lift and drag forces have been measured over a range of angles of attack, and the lift curve shape was similar in all cases. A transient high-lift peak approximately 1.5 times the quasi-steady value occurred in the first chord length of travel, caused by the formation of a strong attached leading edge vortex. This vortex appears to develop and shed more quickly at lower Reynolds numbers. The circulation of the leading edge vortex has been measured and agrees well with force data.     

Effect of apex strake incidence-angle on the vortex development and interaction of a double-delta wing

The vortex flow characteristics of a sharp-edged delta wing with an apex strake was investigated through the visualization and particle image velocimetry (PIV) measurement of the wing-leeward flow region, and the wing-surface pressure measurement. The wing model was a flat-plate, and 65°-sweep cropped-delta wing with sharp leading edges. The apex strake was also a flat-plate wing with a cropped-delta shape of 65°/90° sweep, and it can change its incidence angle. The flow Reynolds number was 2.2 × 10⁵ for the flow visualization and 8.2 × 10⁵ for the PIV and wing-surface pressure measurements. The physics of the vortex flow in the wing-leeward flow region and the suction-pressure distribution on the wing upper-surface were interrelated and analyzed. The effect of a positive (negative) strake incidence-angle was the upward movement of the strake and wing vortices away from (downward movement of the strake and wing vortices toward) the wing-upper surface and the delayed (enhanced) coiling interaction between them. This change of vortex flow characteristics projected directly on the suction pressure distribution on the wing upper-surface.     

Non-uniform recovery of vortex breakdown over delta wing in response to blowing along vortex core

The transient characteristics of vortical structure over delta wing are studied experimentally when subject to single along-core blowing perturbation. Two half delta wing models with different sweep angle = 60° and = 75° are investigated in this study. For = 75°, the transient location of the onset of vortex breakdown moves upstream monotonously toward the unperturbed location. However, for =60°, there exists a chordwise region where the upstream propagation of the onset location of vortex breakdown is temporarily delayed. This delay causes the recovery process (upstream propagation) of the onset location of vortex breakdown to be non-uniform. In fact, this non-uniform recovery may take on different appearance such as a plateau or overshoots and will last for several times of the convective time scaleC/U . Furthermore, the location of chordwise region corresponding to this delay depends strongly upon the angle of attack (AOA). In additions, the non-uniform recovering characteristic of the onset of vortex breakdown may not be observed at high AOA if the blowing rate is too low. The mechanism governing the non-uniform recovering characteristic is clearly verified through the LDA measurement and the phase-locked flow visualization. Evidently the mutual interaction between the primary vortex structure and the secondary vortex is the key mechanism that leads to non-uniform recovering character over delta wing with sweep angle of 60°The authors are grateful for the support of the this investigation from National Science Foundation of the Republic of China under the grant no. NSC-83-0424-E-005-006.     

Flutter suppression of a plate-like wing via parametric excitation

The present work is motivated by the well known stabilizing effect of parametric excitation of some dynamical systems such as the inverted pendulum. The possibility of suppressing wing flutter via parametric excitation along the plane of highest rigidity in the neighborhood of combination resonance is explored. The nonlinear equations of motion in the presence of incompressible fluid flow are derived using Hamilton's principle and Theodorsen's theory for modeling aerodynamic forces. In the presence of air flow, the bending and torsion modes possess nearly the same frequency. Under parametric excitation and in the absence of air flow, each mode oscillates at its own natural frequency. In the neighborhood of combination resonance, the nonlinear response is determined using the multiple scales method at the critical flutter speed and at slightly higher airflow speed. The domains of attraction and bifurcation diagrams are obtained to reveal the conditions under which the parametric excitation can provide stabilizing effect. The basins of attraction for different values of excitation amplitude reveal the stabilizing effect that takes place above a critical excitation level. Below that level, the response experiences limit cycle oscillations, cascade of period doubling, and chaos. For flow speed slightly higher than the critical flutter speed, the response experiences a train of spikes, known as ‘firing,’ a term that is borrowed from neuroscience, followed by ‘refractory’ or recovery effect, up to an excitation level above which the wing is stabilized. The results of the multiple scales method are verified using numerical simulation of the original nonlinear differential equations.     

Numerical capture of wing tip vortex improved by mesh adaptation

This paper presents a mesh adaptation procedure linked to a finite volume solver, the goal of which is to increase the precision of the numerical simulation of a wing tip vortex flow. The adaptation scheme is applied to hexahedron meshes and hybrid meshes made up of tetrahedrons and prisms. To evaluate the ability of each type of element to capture the physics of a tip vortex, a specific test case is studied and results obtained numerically from this test case are compared with experimental results. The error estimator of the adaptation scheme is derived from a solution scalar variable. It is shown that the element anisotropy as well as the adaptation algorithms used have an impact on the precision of the solution. Adaptation of hexahedrons allows a better capture of the tip vortex far from the vortex root, even though the adaptation of those hexahedrons barely changes the number of nodes used to achieve a specified precision, contrary to the adaptation of hybrid meshes. Copyright © 2009 John Wiley & Sons, Ltd.