Aerodynamic Performance of a Gliding Swallowtail Butterfly Wing Model

In the present study, we perform a wind-tunnel experiment to investigate the aerodynamic performance of a gliding swallowtail-butterfly wing model having a low aspect ratio. The drag, lift and pitching moment are directly measured using a 6-axis force/torque sensor. The lift coefficient increases rapidly at attack angles less than 10° and then slowly at larger attack angles. The lift coefficient does not fall off rapidly even at quite high angles of attack, showing the characteristics of low-aspect-ratio wings. On the other hand, the drag coefficient increases more rapidly at higher angles of attack due to the increase in the effective area responsible for the drag. The maximum lift-to-drag ratio of the present modeled swallowtail butterfly wing is larger than those of wings of fruitfly and bumblebee, and even comparable to those of wings of birds such as the petrel and starling. From the measurement of pitching moment, we show that the modeled swallowtail butterfly wing has a longitudinal static stability. Flow visualization shows that the flow separated from the leading edge reattaches on the wing surface at α < 15°, forming a small separation bubble, and full separation occurs at α ≥ 15°. On the other hand, strong wing-tip vortices are observed in the wake at α ≥ 5° and they are an important source of the lift as well as the main reason for broad stall. Finally, in the absence of long hind-wing tails, the lift and longitudinal static stability are reduced, indicating that the hind-wing tails play an important role in enhancing the aerodynamic performance.     

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.     

Aerodynamic forces and flow fields of a two-dimensional hovering wing

This paper reports the results of an experimental investigation on a two-dimensional (2-D) wing undergoing symmetric simple harmonic flapping motion. The purpose of this investigation is to study how flapping frequency (or Reynolds number) and angular amplitude affect aerodynamic force generation and the associated flow field during flapping for Reynolds number (Re) ranging from 663 to 2652, and angular amplitudes (α _A) of 30°, 45° and 60°. Our results support the findings of earlier studies that fluid inertia and leading edge vortices play dominant roles in the generation of aerodynamic forces. More importantly, time-resolved force coefficients during flapping are found to be more sensitive to changes in α _A than in Re. In fact, a subtle change in α _A may lead to considerable changes in the lift and drag coefficients, and there appears to be an optimal mean lift coefficient around α _A = 45°, at least for the range of flow parameters considered here. This optimal condition coincides with the development a reverse Karman Vortex street in the wake, which has a higher jet stream than a vortex dipole at α _A = 30° and a neutral wake structure at α _A = 60°. Although Re has less effect on temporal force coefficients and the associated wake structures, increasing Re tends to equalize mean lift coefficients (and also mean drag coefficients) during downstroke and upstroke, thus suggesting an increasing symmetry in the mean force generation between these strokes. Although the current study deals with a 2-D hovering motion only, the unique force characteristics observed here, particularly their strong dependence on α _A, may also occur in a three-dimensional hovering motion, and flying insects may well have taken advantage of these characteristics to help them to stay aloft and maneuver. An erratum to this article can be found at     

An experimental study of a bio-inspired corrugated airfoil for micro air vehicle applications

An experimental study was conducted to investigate the aerodynamic characteristics of a bio-inspired corrugated airfoil compared with a smooth-surfaced airfoil and a flat plate at the chord Reynolds number of Re _C  = 58,000–125,000 to explore the potential applications of such bio-inspired corrugated airfoils for micro air vehicle designs. In addition to measuring the aerodynamic lift and drag forces acting on the tested airfoils, a digital particle image velocimetry system was used to conduct detailed flowfield measurements to quantify the transient behavior of vortex and turbulent flow structures around the airfoils. The measurement result revealed clearly that the corrugated airfoil has better performance over the smooth-surfaced airfoil and the flat plate in providing higher lift and preventing large-scale flow separation and airfoil stall at low Reynolds numbers (Re _C  < 100,000). While aerodynamic performance of the smooth-surfaced airfoil and the flat plate would vary considerably with the changing of the chord Reynolds numbers, the aerodynamic performance of the corrugated airfoil was found to be almost insensitive to the Reynolds numbers. The detailed flow field measurements were correlated with the aerodynamic force measurement data to elucidate underlying physics to improve our understanding about how and why the corrugation feature found in dragonfly wings holds aerodynamic advantages for low Reynolds number flight applications.     

Time resolved PIV analysis of flow over a NACA 0015 airfoil with Gurney flap

A NACA 0015 airfoil with and without a Gurney flap was studied in a wind tunnel with Re _c = 2.0 × 10⁵ in order to examine the evolving flow structure of the wake through time-resolved PIV and to correlate this structure with time-averaged measurements of the lift coefficient. The Gurney flap, a tab of small length (1–4% of the airfoil chord) that protrudes perpendicular to the chord at the trailing edge, yields a significant and relatively constant lift increment through the linear range of the C _L versus α curve. Two distinct vortex shedding modes were found to exist and interact in the wake downstream of flapped airfoils. The dominant mode resembles a Kàrmàn vortex street shedding behind an asymmetric bluff body. The second mode, which was caused by the intermittent shedding of fluid recirculating in the cavity upstream of the flap, becomes more coherent with increasing angle of attack. For a 4% Gurney flap at α = 8°, the first and second modes corresponded with Strouhal numbers based on flap height of 0.18 and 0.13. Comparison of flow around ‘filled’ and ‘open’ flap configurations suggested that the second shedding mode was responsible for a significant portion of the overall lift increment.     

An experimental comparison of two adaptation strategies in an adaptive-walls wind-tunnel

In the present work an experimental comparison is made between two adaptation strategies: the Judd's method and the Everhart's method. A NACA 0012 airfoil has been tested at Mach numbers up to 0.4: models with chords up to 200 mm have been tested in a 200 mm × 200 mm test section. The two strategies, though based on different theoretical approaches, show a fairly good agreement as far as c _p distribution on the model, lift and drag curves and residual interference are concerned and agree, in terms of lift curve slope and drag coefficient at zero lift, with the McCroskey correlation.A version of this paper was presented at the 2nd International Conference on Experimental Fluid Mechanics, Torino, Italy, July 4–8, 1994.This work has been supported by the Italian Ministry for University and Scientific and Technological Research (MURST 60%).     

Aerodynamic experimentation with ducted models as applied to hypersonic air-breathing vehicles

 A methodology of experimentation in high supersonic wind tunnels for studying aerodynamic characteristics of hypersonic flying vehicles powered by air-breathing engines is discussed. Investigations of such total aerodynamic forces as drag, lift and pitching moment at testing the models are implicit when the air flow through the model ducts is accomplished so that to provide the simulation of the external flow around the airplane and flow over the inlets, but the operating engines and, hence, the exhaust jets are not modeled. The methods used for testing such models are based on the measurement of duct stream parameters alongside with the balance measurement of aerodynamic forces acting on the models. In the tests, aerometric tools are used such as narrow metering nozzles (plugs), pitot and static pressure probes, stagnation temperature probes and pressure orifices in walls of the model duct. The aerometric data serve to determine the flow rate and momentum of the stream at the duct exit. The internal non-simulated forces of the model ducts are also determined using the conservation equations for energy, mass flow and momentum, and these forces are eliminated from the aerodynamic test results. The techniques of the said model testing have been well developed as applied to supersonic aircraft, however their application for hypersonic vehicles whose models are tested at high supersonic speeds, Mach number M _∞>4, implies some specific features. In the present paper, the results of experimental and theoretical study of these features are discussed. Some experimental data on aerodynamics of hypersonic aircraft models received in methodological tests are also presented. The tunnel experiments have been carried out in the Mach number range M _∞=2–6. Received: 25 July 1996/ Accepted: 14 December 1998     

Application of the NS-α Model to a Recirculating Flow

In this paper we investigate a subgrid model based on an anisotropic version of the NS-α model using a lid-driven cavity flow at a Reynolds number of 10,000. Previously the NS-α model has only been used numerically in the isotropic form. The subgrid model is developed from the Eulerian-averaged anisotropic equations (Holm, Physica D 133:215, 1999). It was found that when α ² was based on the mesh numerical oscillations developed which manifested themselves in the appearance of streamwise vortices and a ‘mixing out’ of the velocity profile. This is analogous to the Craik–Leibovich mechanism, with the difference being that the oscillations here are not physical but numerical. The problem could be traced back to the discontinuity in α ² encountered when α ² = 0 on the endwalls. A definition of α ² based on velocity gradients, rather than mesh spacing, is proposed and tested. Using this definition the results with the model show a significant improvement. The splitting of the downstream wall jet, rms and shear stress profiles are correctly captured a coarse mesh. The model is shown to predict both positive and negative energy transfer in the jet impingement region, in qualitative agreement with DNS results.     

Three-dimensional supersonic flow around a pointed body with upstream energy supply

The effect of a local source of energy on a three-dimensional supersonic flow and the aerodynamic characteristics of a pointed ogival body is numerically studied. The results obtained show that the position of the local source of energy upstream of the body on the axis or its deviation from the axis can affect significantly the aerodynamic characteristics of the body (drag, lift, and pitching moment) and the flight trajectory of the vehicle. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translanted from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 5, pp. 116–121, September–October, 1998.     

One Equation Model for Turbulent Channel Flow with Second Order Viscoelastic Corrections

A modified second order viscoelastic constitutive equation is used to derive a k–l type turbulence closure to qualitatively assess the effects of elastic stresses on fully-developed channel flow. Specifically, the second order correction to the Newtonian constitutive equation gives rise to a new term in the momentum equation involving the time-averaged elastic shear stress and in the turbulent kinetic energy transport equation quantifying the interaction between the fluctuating elastic stress and rate of strain tensors, denoted by P _w , for which a closure is developed and tested. This closure is based on arguments of isotropic turbulence and equilibrium in boundary layer flows and a priori P _w could be either positive or negative. When P _w is positive, it acts to reduce the production of turbulent kinetic energy and the turbulence model predictions qualitatively agree with direct numerical simulation (DNS) results obtained for more realistic viscoelastic fluid models with memory which exhibit drag reduction. In contrast, P _w  < 0 leads to a drag increase and numerical breakdown of the model occurs at very low values of the Deborah number, which signifies the ratio of elastic to viscous stresses. Limitations of the turbulence model primarily stem from the inadequacy of the k–l formulation rather than from the closure for P _w . An alternative closure for P _w , mimicking the viscoelastic stress work predicted by DNS using the Finitely Extensible Nonlinear Elastic-Peterlin fluid model, which is mostly characterized by P _w  > 0 but has also a small region of negative P _w in the buffer layer, was also successfully tested. This second model for P _w leads to predictions of drag reduction, in spite of the enhancement of turbulence production very close to the wall, but the equilibrium conditions in the inertial sub-layer were not strictly maintained.     

Irradiation-Induced Solute Clustering in a Low Nickel FeMnNi Ferritic Alloy

Understanding the radiation embrittlement of reactor pressure vessel (RPV) steels is required to be able to operate safely a nuclear power plant or to extend its lifetime. The mechanical properties degradation is partly due to the clustering of solute under irradiation. To gain knowledge about the clustering process, a Fe−1.1 Mn−0.7 Ni (at.%) alloy was irradiated in a test reactor at two fluxes of 0.15 and 9 ×10¹⁷ n _E > 1MeV .m^− 2.s^− 1 and at increasing doses from 0.18 to 1.3 ×10²⁴ n _E > 1MeV .m^− 2 at 300°C. Atom probe tomography (APT) experiments revealed that the irradiation promotes the formation in the α iron matrix of Mn/Mn and/or Ni/Ni pair correlations at low dose and Mn–Ni enriched clusters at high dose. These clusters dissolve partially after a thermal treatment at 400°C. Based on a comparison with thermodynamic calculations, we show that the solute clustering under irradiation can just result from an induced mechanism.     

Aerodynamics of intermittent bounds in flying birds

Flap-bounding is a common flight style in small birds in which flapping phases alternate with flexed-wing bounds. Body lift is predicted to be essential to making this flight style an aerodynamically attractive flight strategy. To elucidate the contributions of the body and tail to lift and drag during the flexed-wing bound phase, we used particle image velocimetry (PIV) and measured properties of the wake of zebra finch (Taeniopygia guttata, N = 5), flying at 6–10 m s⁻¹ in a variable speed wind tunnel as well as flow around taxidermically prepared specimens (N = 4) mounted on a sting instrumented with force transducers. For the specimens, we varied air velocity from 2 to 12 m s⁻¹ and body angle from −15° to 50°. The wake of bounding birds and mounted specimens consisted of a pair of counter-rotating vortices shed into the wake from the tail, with induced downwash in the sagittal plane and upwash in parasagittal planes lateral to the bird. This wake structure was present even when the tail was entirely removed. We observed good agreement between force measures derived from PIV and force transducers over the range of body angles typically used by zebra finch during forward flight. Body lift:drag (L:D) ratios averaged 1.4 in live birds and varied between 1 and 1.5 in specimens at body angles from 10° to 30°. Peak (L:D) ratio was the same in live birds and specimens (1.5) and was exhibited in specimens at body angles of 15° or 20°, consistent with the lower end of body angles utilized during bounds. Increasing flight velocity in live birds caused a decrease in C _L and C _D from maximum values of 1.19 and 0.95 during flight at 6 m s⁻¹ to minimum values of 0.70 and 0.54 during flight at 10 m s⁻¹. Consistent with delta-wing theory as applied to birds with a graduated-tail shape, trimming the tail to 0 and 50% of normal length reduced L:D ratios and extending tail length to 150% of normal increased L:D ratio. As downward induced velocity is present in the sagittal plane during upstroke of flapping flight, we hypothesize that body lift is produced during flapping phases. Future efforts to model the mechanics of intermittent flight should take into account that flap-bounding birds may support up to 20% of their weight even with their wings fully flexed.     

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.     

Separation flow control on a generic ground vehicle using steady microjet arrays

A model of a generic vehicle shape, the Ahmed body with a 25° slant, is equipped with an array of blowing steady microjets 6 mm downstream of the separation line between the roof and the slanted rear window. The goal of the present study is to evaluate the effectiveness of this actuation method in reducing the aerodynamic drag, by reducing or suppressing the 3D closed separation bubble located on the slanted surface. The efficiency of this control approach is quantified with the help of aerodynamic load measurements. The changes in the flow field when control is applied are examined using PIV and wall pressure measurements and skin friction visualisations. By activating the steady microjet array, the drag coefficient was reduced by 9–14% and the lift coefficient up to 42%, depending on the Reynolds number. The strong modification of the flow topology under progressive flow control is particularly studied.     

Experimental studies on drag reduction and rheology of mixed cationic surfactants with different alkyl chain lengths

Experimental studies of the effects of mixtures of cationic surfactants on their drag reduction and rheological behaviors are reported. Cationic alkyl trimethyl quaternary ammonium surfactants with alkyl chain lengths of C₁₂ and C₂₂ were mixed at different molar ratios (total surfactant concentrations were kept at 5 mM with 12.5 mM sodium salicylate (NaSal) as counterion). Drag reduction tests showed that by adding 10% (mol) of C₁₂, the effective drag reduction range expanded to 4–120 °C, compared with 80–130 °C with only the C₂₂ surfactant. Thus mixing cationic surfactants with different alkyl chain lengths is an effective way of tuning the drag reduction temperature range. Cryo-TEM micrographs revealed thread-like micellar networks for surfactant solutions in the drag reducing temperature range, while vesicles were the dominant microstructures at non-drag reducing temperatures. High extensional viscosity was the main rheological feature for all solutions except 50% C₁₂ (mol) solution, which also does not show strong viscoelasticity. It is not clear why this low extensional viscosity solution with relatively weak viscoelasticity is a good drag reducer. Received: 3 November 1999/Accepted: 5 January 2000     

Shear localisation with 2D viscous froth and its relation to the continuum model

Simulations of monodisperse and polydisperse (μ ₂(A) = 0.13±0.002) 2D foam samples undergoing simple shear are performed using the 2D viscous froth (VF) model. These simulations clearly demonstrate shear localisation. The dependence of localisation length on the product λV (shearing velocity V times the wall drag coefficient λ) is examined and is shown to agree qualitatively with published experimental data. A wide range of localisation lengths is found at low λV, an effect which is attributed to the existence of distinct yield and limit stresses. The general continuum model is extended to incorporate such an effect and its parameters are subsequently related to those of the VF model. A Herschel–Bulkley exponent of a = 0.3 is shown to accurately describe the observed behaviour. The localisation length is found to be independent of λV for monodisperse foam samples.     

Drag characteristics of a turbulent boundary layer over a flat plate with transverse square grooves

 The effect of transverse square grooves on the drag characteristics of a flat plate has been investigated using direct drag measurements in a tow tank. Two flat plates with 2.5 and 5 mm transverse square grooves, with groove spacing (s) to groove width (w) ratio, s/w=10, 20 and 40, were tested at plate length Reynolds numbers in the range 2–10×10⁶. There is an increase in drag over the “smooth” plate for all s/w configurations. Received: 31 December 1997/Accepted: 4 April 1998     

Temperature dependence of the interfacial tension of PS/PMMA, PS/PE, and PMMA/PE blends

The effect of temperature on the interfacial tension for PS/PMMA, PS/PE, and PMMA/PE was measured using the imbedded fiber retraction method. Interfacial tensions for PS/PMMA, PS/PE, and PMMA/PE were measured over temperature ranges of 160–250 °C, 140–220 °C, and 140–220 °C, respectively. The interfacial tension was found to follow a dependence of 3.6–0.013 T dyn/cm, 7.6–0.051 T dyn/cm and 11.8–0.017 T dyn/cm for PS/PMMA, PS/PE, and PMMA/PE, respectively. Comparison of the data with the mean field theory of Helfand and Sapse were made; however, a simple linear fit to the data described the temperature dependence in the experimental window as well as the predictions of the mean field theory. Received: 6 July 1999 Accepted: 23 March 2000     

An assessment of pressure drop and void fraction correlations with data from two-phase natural circulation loops

 Void fraction and pressure drop correlations play an important role in predicting the performance of natural circulation loops. Hence an assessment of the commonly used and often cited correlations for pressure drop and void fraction has been carried out with data from natural circulation loops. This assessment considered 33 void fraction correlations and 14 pressure drop correlations. The void fraction correlations were initially tested against the various limiting conditions. Only 14 correlations were found to satisfy at least two limiting conditions (i.e., at x = 0; α = 0 and at x = 1; α = 1) and were assessed against the data. This assessment showed that the Chexal et al. (1996) correlation is better than all the others considered. The assessment of pressure drop correlations were carried out with the Chexal et al. (1996) correlation for void fraction and Saha-Zuber model for the onset of subcooled boiling. This assessment showed that most correlations give predictions close to each other. Received on 17 January 2000     

Evaluation of drag correction factor for spheres settling in associative polymers

Drag correction factors are calculated for the creeping motion of spheres descending in various associative polymers of different concentration with various sphere-container ratios and Weissenberg numbers. The simple-shear rheology and linear viscoelasticity of these polymeric fluids have been previously presented and modeled with the BMP (Bautista–Manero–Puig) equation of state (Mendoza-Fuentes et al., Phys Fluids 21:033104, 2009). The drag on the sphere is initially kept nearly constant for small Weissenberg numbers, We < 0.1. As the Weissenberg number increases, We < 0.1, a reduction in drag is found. Experimental results show the presence of a critical Weissenberg number at which a drag reduction occurs. The reduction in the drag correction factor is associated to the onset of extension-thinning, which coincides with the formation of a negative wake. No increase in the drag correction factor was observed, due to the simultaneous opposing effects of extension-thickening and shear-thinning viscosity. The shape of the drag correction factor curve may be predicted considering the extensional properties of the solutions, as suggested elsewhere (Chen and Rothstein, J Non-Newton Fluid Mech 116:205–215, 2004).