Computational study of a pitching bio-inspired corrugated airfoil

The phenomenon of insect flight has been of scientific interest for many years and is more recently inspiring modern engineering devices such as Micro Aerial Vehicles (MAVs). Insect flight is characterized by unsteady fluid dynamics at low Reynolds numbers. The importance of viscous effects to the successful flapping flight of insects has been identified and with the current state of computing and Computational Fluid Dynamics (CFD) these effects can now be studied in detail. The present work attempts to simplify this complex phenomenon by considering symmetric oscillating rotational motion of a wing (pitching). What is of interest in this study is how the shape of a corrugated idealized insect wing affects the performance and flow characteristics around the pitching wing. Two dimensional CFD on an oscillating wing has been performed and reported. Measurements were taken to ensure the accuracy of the computational solution and the results validated against experimental PIV results. A range of frequencies and rotational amplitudes have been investigated. Lift and drag coefficients have been analyzed for all cases to quantify the effects of unsteady flow features on the performance of the oscillating wing. It was found that the wing shape used in this study resulted in the viscous features formed on the top of the wing exhibiting high sensitivity to the oscillating conditions and these influenced the performance of the wing. The flow features formed on the bottom of the wing remained similar throughout the cases tested. In the pitching regime this wing profile did not perform as well as published results for smooth airfoils in terms of thrust and propulsion efficiency. However this may be due to reduced frequency effects becoming important at our high pitching amplitude which need to be investigated further. There may be other oscillatory regimes that more accurately represent flapping flight in which the corrugated foil outperforms a smooth counterpart but these are yet to be investigated. Further research in this area may help answer the question as to how evolutionarily significant other benefits of a corrugated wing, such as being light and strong, are compared to its aerodynamic properties, the present results seem to favor the former.     

Dynamics of co-rotating vortices in a flow around a bio-inspired corrugated airfoil

Bio-inspired corrugated airfoils show favourable aerodynamic characteristics such as high coefficient of lift and delayed stall at low Reynolds numbers. Two-dimensional (2D) direct numerical simulation has been performed here on a corrugated airfoil at various angles of attack (0°, +5°, -5°) and Reynolds number of 280 to 6700. The objective is to analyse the pressure variation inside the corrugations and correlate it to the vortex movement across the corrugations and the overall aerodynamic characteristics of the corrugated airfoil. The flow characteristics have been examined based on the local Strouhal numbers in the corrugations of the airfoil. It is observed that the pressure variation in each corrugation is the result of vortex merging and separation in the corrugation which plays a major role in changing the flow characteristics. The Strouhal number of the flow is dictated by the most dominant local Strouhal number. The numerical results are further compared with experimental results obtained using particle image velocimetry, and the two set of results are found to match well. These results are significant because they elucidate the effect of corrugation, angle of attack, and Reynolds number on flow over a corrugated airfoil.     

Virtual flight simulation of a dual rotor micro air vehicle

In this paper, a new computational method is developed based on computational fluid dynamics (CFD) coupled with rigid body dynamics (RBD) and flight control law in an in-house programmed source code. The CFD solver is established based on momentum source method, preconditioning method, lower–upper symmetric Gauss–Seidel iteration method, and moving overset grid method. Two-equation shear–stress transport k ? ω turbulence model is employed to close the governing equations. Third-order Adams prediction-correction method is used to couple CFD and RBD in the inner iteration. The wing-rock motion of the delta wing is simulated to validate the capability of the computational method for virtual flight simulation. Finally, the developed computational method is employed to simulate the longitudinal virtual flight of a dual rotor micro air vehicle (MAV). Results show that the computational method can simulate the virtual flight of the dual rotor MAV.     

An experimental investigation of the flow structure over a corrugated waveform in a transpired air collector

An experimental investigation of the flow dynamics in a channel with a corrugated surface is presented. Particle image velocimetry was used to obtain two-dimensional velocity fields at three different locations along the channel length, over a range of Reynolds numbers. The results show a significant impact of the corrugation waveform on the mean and turbulent flow structure inside the channel. Strong bursting flow originating from the trough, sweeping flow from the bulk region and the vortex shedding off the crest were observed. Their interactions created a complex three-dimensional flow structure extended over almost the entire channel. The mean velocity profiles indicate a strong diffusion of shear. The profiles of various turbulent properties show the enhancement of turbulence in the vicinity of the waveform. It was found that the turbulence in the channel was almost entirely produced in this region above the corrugation trough. Significant momentum transfer from the corrugation wall by the turbulent velocity field was also observed. The mean and turbulent flow behaviour was found to be periodic with respect to the waveform over most of the channel length. The results show the presence of strong turbulence even at the Reynolds number that falls within the conventional laminar range.     

An extensive FPGA-based realization study about the Izhikevich neurons and their bio-inspired applications

This study is aimed at analysing damping and gyroscopic effects on the stability of parametrically excited continuous rotor systems, taking into account both external (non-rotating) and internal (rotating) damping distributions. As case-study giving rise to a set of coupled differential Mathieu–Hill equations with both damping and gyroscopic terms, a balanced shaft is considered, modelled as a spinning Timoshenko beam loaded by oscillating axial end thrust and twisting moment, with the possibility of carrying additional inertial elements like discs or flywheels. After discretization of the equations of motion into a set of coupled ordinary differential Mathieu–Hill equations, stability is studied via eigenproblem formulation, obtained by applying the harmonic balance method. The occurrence of simple and combination parametric resonances is analysed introducing the notion of characteristic circle on the complex plane and deriving analytical expressions for critical solutions, including combination parametric resonances, valid for a large class of rotors. A numerical algorithm is then developed for computing global stability thresholds in the presence of both damping and gyroscopic terms, also valid when closed-form expressions of critical solutions do not exist. The influence on stability of damping distributions and gyroscopic actions is then analysed with respect to frequency and amplitude of the external loads on stability charts in the form of Ince–Strutt diagrams.

     

Investigation on the planform and kinematic optimization of bio-inspired nano air vehicles for hovering applications

Wing shape and kinematics of flapping wing nano air vehicles are two important factors in their design process. These factors require an optimal design in terms of decreasing the needed aerodynamic power. Since, insects are regarded as the best natural flier in hovering flight, seven of their wings are considered in order to determine the best wing shape for hovering applications. Because of the difference in the original bio-inspired shape of these wings, two scenarios are studied, namely, considering the same wingspan and same wing surface. Using the quasi-steady approximation to model the aerodynamic loads and a basic gradient approach to optimize the kinematics of the wing, the optimum Euler angles, required aerodynamic power, and hence the best wing shape for each scenario are analytically determined. The results show that the wing shape and surface strongly impact the aerodynamic characteristics and performances of the chosen wing shapes. It is demonstrated that the twisted parasite wing shape is a good candidate to minimize the required aerodynamic power during hovering. The strategy used in this analysis can be used to evaluate the performance of any realistic wing shape design and could provide a guideline for selecting the best wing shape and kinematics for flapping wing nano air vehicles with hovering capabilities.     

An experimental study on boundary layer transition detection over a pitching supercritical airfoil using hot-film sensors

In the present work, experimental tests are conducted to study boundary layer transition over a supercritical airfoil undergoing pitch oscillations using hot-film sensors. Tests have been undertaken at an incompressible flow. Three reduced frequencies of oscillations and two mean angles of attack are studied and the influences of those parameters on transition location are discussed. Different algorithms are examined on the hot-film signals to detect the transition point. Results show the formation of a laminar separation bubble near the leading edge and at relatively higher angles of attack which leads to the transition of the boundary layer. However, at lower angles of attack, the amplification of the peaks in voltage signal indicate the emergence of the vortical structures within the boundary layer, introducing a different transition mechanism. Moreover, an increase in reduced frequency leads to a delay in transition onset, postponing it to a higher angle of attack, which widens the hysteresis between the upstroke and downstroke motions. Rising the reduced frequency yields in weakening or omission of vortical disturbances ensuing the removal of spikes in the signals. Of the other important results observed, is faster movement of the relaminarization point in the higher mean angle of attack. Finally, a time–frequency analysis of the hot-film signals is performed to investigate evolution of spectral features of the transition due to the pitching motion. An asymmetry is clearly observed in frequency pattern of the signals far from the bubble zone towards the trailing edge; this may reflect the difference between the transition and relaminarization physics. Also, various ranges of frequency were obtained for different transition mechanisms.     

An experimental study of the effects of pitch-pivot-point location on the propulsion performance of a pitching airfoil

An experimental investigation was conducted to characterize the evolution of the unsteady vortex structures in the wake of a pitching airfoil with the pitch-pivot-point moving from 0.16C to 0.52C (C is the chord length of the airfoil). The experimental study was conducted in a low-speed wind tunnel with a symmetric NACA0012 airfoil model in pitching motion under different pitching kinematics (i.e., reduced frequency k=3.8–13.2). A high-resolution particle image velocimetry (PIV) system was used to conduct detailed flow field measurements to quantify the characteristics of the wake flow and the resultant propulsion performance of the pitching airfoil. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged velocity distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about the behavior of the unsteady vortex structures. Both the vorticity–moment theorem and the integral momentum theorem were used to evaluate the effects of the pitch-pivot-point location on the propulsion performance of the pitching airfoil. It was found that the pitch-pivot-point would affect the evolution of the unsteady wake vortices and resultant propulsion performance of the pitching airfoil greatly. Moving the pitch-pivot-point of the pitching airfoil can be considered as adding a plunging motion to the original pitching motion. With the pitch-pivot-point moving forward (or backward), the added plunging motion would make the airfoil trailing edge moving in the same (or opposite) direction as of the original pitching motion, which resulted in the generated wake vortices and resultant thrust enhanced (or weakened) by the added plunging motion.     

Effects of camber angle on aerodynamic performance of flapping-wing micro air vehicle

This study investigates the unsteady aerodynamic characteristics of the cambered wings of a flapping-wing micro air vehicle (FW-MAV) in hover. A three-dimensional fluid–structure interaction solver is developed for a realistic modeling of large-deforming wing structure and geometry. Cross-validation is conducted against the experimental results obtained also in the present study to establish more accurate analyses of cambered wings. An investigation is carried out on the unsteady vortex structures around the wings caused by the passive twisting motion. A parametric study is then conducted to evaluate the aerodynamic performance with respect to the camber angle at three different flapping frequencies including normal operating conditions. The camber angles producing the largest thrust and highest propulsive efficiency are estimated at each flapping frequency, and their effects on aerodynamic performance are analyzed in terms of the stroke phase. The timing and magnitude of the passive twisting motion, which are dependent on the camber angle at the operating frequency, greatly affects the unsteady vortex structure. Consequently, the camber angle designed at the operating frequency plays a key role in enhancing the aerodynamic performance of FW-MAVs.     

Design,manufacturing, and flight testing of a fixed wing micro air vehicle with Zimmerman planform

An optimized and comprehensive method is used to design and manufacture a fixed wing micro air vehicle (MAV) with Zimmerman planform. The design process includes four stages which are the specification of the flight mission, determination of the best aspect ratio, identification of the optimum wing loading and thrust loading values, and estimation of the weight of the structural components of the MAV. To do this, various statistical and analytical methods are utilized. Based on an aerodynamic analysis, the results show that an optimum aspect ratio that maximizes the performance of the Zimmerman MAV for a well-defined cruise speed is determined. Considering six possible flights, a constraint analysis is performed and an optimum wing loading value is determined. It is shown that the computational method is beneficial to determine the exact masses for the structural components including the wing, fuselage, and vertical tail. Using the 3D panel method, the determination of the shape of the reflexed airfoil for the MAV is successfully done by minimizing the drag force and the angle of attack to use less powerful motor and avoid any stall effect, respectively. A stability analysis is then performed to check the safe flight of the designed vehicle. During test flight, the results show that the designed Zimmerman MAV satisfies the pre-defined specification. The final characteristics of the manufactured MAV are: wingspan of 44 cm, weight of 450 g, aspect ratio of 1.51, cruise speed of 20 m/s, and flight endurance of 20 min.     

The influence of airfoil kinematics on the formation of leading-edge vortices in bio-inspired flight

The formation process of leading-edge vortices has been investigated experimentally using Particle Image Velocimetry. Various airfoil kinematics have been tested, including asymmetric and peak-shifted plunging motions, and are evaluated for Re = 30,000 and a reduced frequency range of 0.2 ≤ k ≤ 0.33. By measuring the growth in the leading-edge vortex during the dynamic-stall process, the vortex pinch-off process is examined based on the concept of an optimal vortex formation time. The various kinematics are then evaluated with respect to their associated vortex strength, timing and convection into the wake.     

An experimental and numerical study into turbulent condensing steam jets in air

 Temperatures, velocities, and droplet sizes are measured in turbulent condensing steam jets produced by a facial sauna, for varying nozzle diameters and varying initial velocities (Re=3,600–9,200). The release of latent heat due to droplet condensation causes the temperature in the two-phase jet to be significantly higher than in a single-phase jet. At some distance from the nozzle, droplets reach a maximum size and start to evaporate again, which results in a change in sign of latent heat release. The distance of maximum size is determined from droplet size measurements. The experimental results are compared with semi-analytical expressions and with a fully coupled numerical model of the turbulent condensing steam jet. The increase in centreline temperature due to droplet condensation is successfully predicted. Received: 5 April 2000 / Accepted: 15 November 2000     

An experimental study on developing air-water two-phase flow along a large vertical pipe: effect of air injection method

The flow structure in a developing air-water two-phase flow was investigated experimentally along a large vertical pipe (inner diameter, D_h: 0.48 m, ratio of length of flow path L to D_h: about 4.2). Two air injection methods (porous sinter injection and nozzle injection) were adopted to realize an extremely different flow structure in the developing region. The flow rate condition in the test section was as follows: superficial air velocity: 0.02–0.87 m/s (at atmospheric pressure) and superficial water velocity: 0.01–0.2 0.01–0.2 m/s, which covers the range of bubbly to slug flow in a small-scale pipe (D_h about 0.05 m).

No air slugs occupying the flow path were recognized in this experiment regardless of the air injection methods even under the condition where slug flow is realized in the small-scale pipe. In the lower half of the test section, the axial distribution of sectional differential pressure and the radial distribution of local void fraction showed peculiar distributions depending on the air injection methods. However, in the upper half of the test section, the effects of the air injection methods are small in respect of the shapes of the differential pressure distribution and the phase distribution. The comparison of sectional void fraction near the top of the test section with Kataoka's correlation indicated that the distribution parameter of the drift-flux model should be modeled including the effect of D_h and the bubble size distribution is affected by the air injection methods. The bubble size distribution is considered to be affected also by L/D_h based on comparison of results with Hills' correlation.     

An experimental investigation on the surface water transport process over an airfoil by using a digital image projection technique

In the present study, an experimental investigation was conducted to characterize the transient behavior of the surface water film and rivulet flows driven by boundary layer airflows over a NACA0012 airfoil in order to elucidate underlying physics of the important micro-physical processes pertinent to aircraft icing phenomena. A digital image projection (DIP) technique was developed to quantitatively measure the film thickness distribution of the surface water film/rivulet flows over the airfoil at different test conditions. The time-resolved DIP measurements reveal that micro-sized water droplets carried by the oncoming airflow impinged onto the airfoil surface, mainly in the region near the airfoil leading edge. After impingement, the water droplets formed thin water film that runs back over the airfoil surface, driven by the boundary layer airflow. As the water film advanced downstream, the contact line was found to bugle locally and developed into isolated water rivulets further downstream. The front lobes of the rivulets quickly advanced along the airfoil and then shed from the airfoil trailing edge, resulting in isolated water transport channels over the airfoil surface. The water channels were responsible for transporting the water mass impinging at the airfoil leading edge. Additionally, the transition location of the surface water transport process from film flows to rivulet flows was found to occur further upstream with increasing velocity of the oncoming airflow. The thickness of the water film/rivulet flows was found to increase monotonically with the increasing distance away from the airfoil leading edge. The runback velocity of the water rivulets was found to increase rapidly with the increasing airflow velocity, while the rivulet width and the gap between the neighboring rivulets decreased as the airflow velocity increased.     

An experimental study on the effect of air bubble injection on the flow induced rotational hub

Modification of shear stress due to air bubbles injection in a rotary device was investigated experimentally. Air bubbles inject to the water flow crosses the neighbor of the hub which can rotate just by water flow shear stresses, in this device. Increasing air void fraction leads to decrease of shear stresses exerted on the hub surface until in high void fractions, the hub motion stopped as observed. Amount of skin friction decrease has been estimated by counting central hub rotations. Wall shear stress was decreased by bubble injection in all range of tested Reynolds number, changing from 50,378 to 71,238, and also by increasing air void fraction from zero to 3.06%. Skin friction reduction more than 85% was achieved in this study as maximum measured volume of air fraction injected to fluid flow while bubbles are distinct and they do not make a gas layer. Significant skin friction reduction obtained in this special case indicate that using small amount of bubble injection causes large amount of skin friction reduction in some rotary parts in the liquid phases like as water.     

Triaxial-in-a-plane soil stress gage for vehicle mobility applications

The approach used here consists of an axisymmetric gage with a tapered exterior. Raw outputs are normal stress in three directions in the plane perpendicular to the gage axis of symmetry. From these outputs, the time history of the complete state of stress in this plane can be determined. Of special interest is the plane of symmetry which is the vertical plane centered on one side of a wheeled or tracked vehicle proceeding in a straight line. The gage was placed from a berm on the side of the vehicle path, approximately 4 ft horizontally and at approximately 9 in. below the ground surface. Test vehicles were a 4 × 4 wheeled vehicle and a M113A2 armored personnel carrier. The measured stress results are largely consistent with expectations.     

Stabilization of flow boiling in a micro tube with air injection

Air injection as a stabilization method is evaluated for flow boiling in a micro tube. Pyrex glass tube coated by ITO film is employed as a test tube for flow visualization with water as a working fluid. Air bubble and liquid slug lengths are controlled by changing air and liquid mass velocities. Wall temperatures and inlet/outlet pressures show very large fluctuations during flow boiling without air injection. Severe reverse flow is also observed from flow visualization. On the other hand, wall temperature and inlet/outlet pressures as well as visualized flow patterns become very stable with air injection. In addition, much higher heat transfer coefficients are obtained for air injected cases. It is observed from the flow visualization that the flow becomes much stable and shows regular patterns.     

An experimental investigation of two-phase air–water flow through a horizontal circular micro-channel

This paper is a continuation of the authors’ previous work. Two-phase air–water flow experiments are performed in a horizontal circular micro-channel. The test section is made of a fused silica tube with an inner diameter of 0.15 mm and a length of 104 mm. The flow phenomena, which are liquid/unstable annular alternating flow (LUAAF), liquid/annular alternating flow (LAAF), and annular flow, are observed and recorded by a high-speed camera mounted together with a stereozoom microscope. A flow pattern map is presented in terms of the phase superficial velocities and is compared with those of other researchers obtained from different working fluids. Image analysis is performed to determine the void fraction, which increases non-linearly with increasing volumetric quality. It is revealed that the two-phase frictional multiplier data show a dependence on flow pattern rather than mass flux. Based on the present data, a new pressure drop correlation is proposed for practical applications. According to the present study, in general the data for the two-phase air–water flow characteristics are found to comply with those of working fluids other than air–water mixtures.     

The air cushion vehicle as a load-spreading transport device

The most effective fields of operation of the ACV are discussed, and it is shown that the vehicle effectively fills the role of a heavy transporter over weak terrain. A speculative discussion of the interaction between ACV skirts and various terrain is given in some detail. In conclusion short notes are given on some current experimental programmes, and on the work of some actual vehicles in the field.     

An experimental study of exit pressures for polymer melts

A slit die apparatus is used to measure exit pressures for five different polymer melts. Viscosity data obtained from the same apparatus agree well with values obtained from a cone-and-plate rheometer or a capillary rheometer. Except for a PVC sample where thermal degradation was found to occur, the exit pressures obtained by linear extrapolation of the measured pressure profiles are all positive, and increase with increasing shear stress. The values of the first normal stress difference calculated according to the exit pressure theory are of the right order of magnitude and in some cases correlate satisfactorily with values measured in a cone-and-plate rheometer. However, the high sensitivity of the exit pressure values to the method of extrapolation and the wild scatter of exit pressure data for some materials make it difficult to use the exit pressure method as a routine procedure for accurate determination of the first normal stress difference.