Nonlinear aeroelasticity of rotating and flapping wings - a review

The interaction between large deflections, rotation effects and unsteady aerodynamics makes the dynamic analysis of rotating and flapping wing a nonlinear aeroelastic problem. This problem is governed by nonlinear periodic partial differential equations whose solution is needed to calculate the response and loads acting on vehicles using rotary or flapping wings for lift generation. We look at three important problems in this paper. The first problem shows the effect of nonlinear phenomenon coming from piezoelectric actuators used for helicopter vibration control. The second problem looks at the propagation on material uncertainty on the nonlinear response, vibration and aeroelastic stability of a composite helicopter rotor. The third problem considers the use of piezoelectric actuators for generating large motions in a dragonfly inspired flapping wing. These problems provide interesting insights into nonlinear aeroelasticity and show the likelihood of surprising phenomenon which needs to be considered during the design of rotary and flapping wing vehicle

Optimum atomic spacing for AlAs etching in GaAs epitaxial lift-off technology

With respect to GaAs epitaxial lift-off technology, we report here the optimum atomic spacing (5-10 nm) needed to etch off the AlAs release layer that is sandwiched between two GaAs epitaxial layers. The AlAs etching rate in hydrofluoric acid based solutions was monitored as a function of release layer thickness. We found a sudden quenching in the etching rate, approximately 20 times that of the peak value, at lower dimensions (similar to2.5 nm) of the AlAs epitaxial layer. Since this cannot be explained on the basis of a previous theory (inverse square root of release layer thickness), we propose a diffusion-limited mechanism to explain this reaction process. With the diffusion constant being a mean-free-path-dependent parameter, a relation between the mean free path and the width of the channel is considered. This relation is in reasonable agreement with the experimental results and gives a good physical insight to the reaction kinetics.

KFVS (Kinetic Flux Vector Splitting) on moving grids for unsteady aerodynamics

Accurate numerical solutions to the problems in fluid-structure (aeroelasticity) interaction are becoming increasingly important in recent years. The methods based on FCD (Fixed Computational Domain) and ALE (Alternate Lagrangian Eulerian) to solve such problems suffer from numerical instability and loss of accuracy. They are not general and can not be extended to the flowsolvers on unstructured meshes. Also, global upwind schemes can not be used in ALE formulation thus leads to the development of flow solvers on moving grids. The KFVS method has been shown to be easily amenable on moving grids required in unsteady aerodynamics. The ability of KFMG (Kinetic Flux vector splitting on Moving Grid) Euler solver in capturing shocks, expansion waves with small and very large pressure ratios and contact discontinuities has been demonstrated.

Aerodynamics of delta wings with flaps at hypersonic speeds

Control surface effectiveness is an important parameter for any aeroplane. For a hypersonic aircraft, though the power required to operate the flaps is determined by low speed flying conditions, it is imperative to know the effect of flaps at hypersonic speeds. Hence, studies have been done on this topic by aerodynamicists for over 40 years. In spite of this, only a limited data is available in the literature on this subject. This paper discusses the experimental study of the effect of sweep on the aerodynamic characteristics of thin slab delta wings with flaps at hypersonic speeds. For the purpose of this investigation, a novel special thin six-component balance, which has a thickness of 4mm and can be housed inside wings with 8mm thickness, has been designed. The wings had a sweep of 76degrees, 70degrees and 65degrees, t/c of 0.053 and flaps with 12% of wing area and 12% of wing chord. Testing were done at Mach 8.2, Re number of 2.13 x 10(6) (based on chord), from alpha = -12degrees to 12degrees and flap angle of 20degrees, 30degrees and 40degrees. Separation lengths, measured from Schlieren pictures, clearly show that there is 'no appreciable' effect of sweep on them. Also, using a simple local flow field calculation, the separation has been identified to be transitional in nature. These features of separation reflect in the force data. Because of the small separation length, the flaps (inspite of their small size) were very effective in generating additional C-N, C-M and C-l, which increased with increase in flap angle. In general, the C-N, C-M and X-CP were unaffected by sweep for symmetric flap deflection at positive incidences and asymmetric flap case, For symmetric flap case at negative incidences, only C-N was not influenced by the sweep but C-M decreased and X-CP moved upstream as the sweep is decreased, The wing with lower sweep produces higher CA and lower (L/D)(max) for both symmetric and asymmetric flaps. The rolling moment and adverse yaw increased with decrease in sweep for asymmetric flap deflection. Newtonian theory is shown to be incapable of predicting the effect of sweep on C-l, C-n and on the incremental values of C-N, C-M and C-A. In conclusion, it can be said that a small flap is generally adequate for hypersonic aeroplanes provided they operate at altitudes where transitional and turbulent separation can be expected to occur. This would make the flaps effective and thus enable ample control authority.

Discrete Vortex Method-Based Model for Ground-Effect Studies

A discrete vortex method-based model has been proposed for two-dimensional/three-dimensional ground-effect prediction. The model merely requires two-dimensional sectional aerodynamics in free flight. This free-flight data can be obtained either from experiments or a high-fidelity computational fluid dynamics solver. The first step of this two-step model involves a constrained optimization procedure that modifies the vortex distribution on the camber line as obtained from a discrete vortex method to match the free-flight data from experiments/computational fluid dynamics. In the second step, the vortex distribution thus obtained is further modified to account for the presence of the ground plane within a discrete vortex method-based framework. Whereas the predictability of the lift appears as a natural extension, the drag predictability within a potential flow framework is achieved through the introduction of what are referred to as drag panels. The need for the use of the generalized Kutta-Joukowski theorem is emphasized. The extension of the model to three dimensions is by the way of using the numerical lifting-line theory that allows for wing sweep. The model is extensively validated for both two-dimensional and three-dimensional ground-effect studies. The work also demonstrates the ability of the model to predict lift and drag coefficients of a high-lift wing in ground effect to about 2 and 8% accuracy, respectively, as compared to the results obtained using a Reynolds-averaged Navier-Stokes solver involving grids with several million volumes. The model shows a lot of promise in design, particularly during the early phase.

CMC simulations of lifted turbulent jet flame in a vitiated coflow

Lifted turbulent jet diffusion flame is simulated using Conditional Moment Closure (CMC). Specifically, the burner configuration of Cabra et al. [R. Cabra, T. Myhrvold, J.Y. Chen. R.W. Dibble, A.N. Karpetis, R.S. Barlow, Proc. Combust. Inst. 29 (2002) 1881-1887] is chosen to investigate H-2/N-2 jet flame supported by a vitiated coflow of products of lean H-2/air combustion. A 2D, axisymmetric flow-model fully coupled with the scalar fields, is employed. A detailed chemical kinetic scheme is included, and first order CIVIC is applied. Simulations are carried out for different jet velocities and coflow temperatures (T-c) The predicted liftoff generally agrees with experimental data, as well as joint-PDF results. Profiles of mean scalar fluxes in the mixture fraction space, for T-c = 1025 and 1080 K reveal that (1) Inside the flame zone, the chemical term balances the molecular diffusion term, and hence the Structure is of a diffusion flamelet for both cases. (2) In the pre-flame zone, the structure depends on the coflow temperature: for the 1025 K case, the chemical term being small, the advective term balances the axial turbulent diffusion term. However, for the 1080 K case. the chemical term is large and balances the advective term, the axial turbulent diffusion term being small. It is concluded that, lift-off is controlled (a) by turbulent premixed flame propagation for low coflow temperature while (b) by autoignition for high coflow temperature. (C) 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Analysis of an internally mountable accelerometer balance system for use with non-isotropic models in shock tunnels

Knowledge of drag force is an important design parameter in aerodynamics. Measurement of aerodynamic forces at hypersonic speed is a challenge and usually ground test facilities like shock tunnels are used to carry out such tests. Accelerometer based force balances are commonly employed for measuring aerodynamic drag around bodies in hypersonic shock tunnels. In this study, we present an analysis of the effect of model material on the performance of an accelerometer balance used for measurement of drag in impulse facilities. From the experimental studies performed on models constructed out of Bakelite HYLEM and Aluminum, it is clear that the rigid body assumption does not hold good during the short testing duration available in shock tunnels. This is notwithstanding the fact that the rubber bush used for supporting the model allows unconstrained motion of the model during the short testing time available in the shock tunnel. The vibrations induced in the model on impact loading in the shock tunnel are damped out in metallic model, resulting in a smooth acceleration signal, while the signal become noisy and non-linear when we use non-isotropic materials like Bakelite HYLEM. This also implies that careful analysis and proper data reduction methodologies are necessary for measuring aerodynamic drag for non-metallic models in shock tunnels. The results from the drag measurements carried out using a 60 degrees half angle blunt cone is given in the present analysis.

Investigation of Missile-Shaped Body with Forward-Facing Cavity at Mach 8

A blunt-nosed hypersonic missile mounted with a forward-facing cavity is a good alternative to reduce the stagnation heating rates. The effects of a forward-racing cavity on heat transfer and aerodynamic coefficients are addressed in this paper. Tests were carried out in hypersonic shock tunnel HST2, at a hypersonic Mach number of 8 using a 41 deg apex-angle blunt cone. The aerodynamic forces on the test model with and without a forward-facing cavity at various angles of attack are measured by using an internally mountable accelerometer force balance system. Heat flux measurements have been carried out on the test model with and without a forward-facing cavity of the entire surface at zero degree angle of attack with platinum sensors. A numerical simulation was also carried out using the computational fluid dynamics code (CFX-Ansys 5.7). An important result of this study is that the smaller cavity diameter has the highest lift-to-drag ratio, whereas the medium cavity has the highest heat flux reduction. Theshock structure around the test model has also been visualized using the Schlieren flow visualization technique. The visualized shock structure and the measured aerodynamic forces on the missile-shaped body with cavity configurations agree well with the axisymmetric numerical simulations.

Force Measurements on a Reflex Cambered Delta Wing

AREFLEX spanwise cambered delta wing with a conical camber designed for M= 1.4, using the method of Ref. 1, was tested at the design Mach number as well as off-design Mach number M=0.15 and 2.3, respectively. The test results are compared with those of a plane wing and also with the available theoretical results at the design condition. At subsonic speed, the cambered wing has less lift at a given incidence and higher lift-to-drag ratio at a given lift than the plane wing, while at supersonic speeds, both of these quantities were less on the cambered wing. At supersonic speed, at the design incidence and Mach number, there is good agreement between results from theory and experiment. The center of pressure on the cambered wing is ahead of that on the plane wing at subsonic speed, while the reverse is true at supersonic speeds. Finally, it is found that over a useful range of lift the cambered wing is aerodynamically more efficient at subsonic speeds, and less so at supersonic speeds, than the plane wing.

A note on mass transfer in turbulent wall jets

THE PROCESS of mass transfer from saturated porous surfaces virtual origin ; exposed to turbulent air streams finds many practical applitransverse coordinate; cations. In many cases, the air stream will be in the form of a height of nozzle above flat plate--radial wall jet; wall jet over the porous surface. The aerodynamics of both plane and radial wall jets have been investigated in detail and a vast amount of literature is available on the subject [l-3].

Shock tunnel testing of a Mach 6 hypersonic waverider

A waverider is a lifting body configuration whose upper surface is parallel to the free stream, and the lower surface aerodynamically so designed, that the resulting shock at the design Mach number, is always attached with the leading edge of the vehicle. This prevents spillage from high pressure (lower) surface to the low pressure (upper) surface.In the present study a conical waverider has been designed, fabricated and tested at Mach 6 in the IISc hypersonic shock tunnel HST2. The measurements show that the waverider has a lift to drag ratio of 4.28 at the designed Mach number. Exhaustive FEM and CFD studies are also carried out to complement the force measurements in the tunnel.

Determination of optimal trajectory under design constraints for a satellite launch vehicle

An optimal pitch steering programme of a solid-fuel satellite launch vehicle to maximize either (1) the injection velocity at a given altitude, or (2) the size of circular orbit, for a given payload is presented. The two-dimensional model includes the rotation of atmosphere with the Earth, the vehicle's lift and drag, variation of thrust with time and altitude, inverse-square gravitational field, and the specified initial vertical take-off. The inequality constraints on the aerodynamic load, control force, and turning rates are also imposed. Using the properties of the central force motion the terminal constraint conditions at coast apogee are transferred to the penultimate stage burnout. Such a transformation converts a time-free problem into a time-fixed one, reduces the number of terminal constraints, improves accuracy, besides demanding less computer memory and time. The adjoint equations are developed in a compact matrix form. The problem is solved on an IBM 360/44 computer using a steepest ascent algorithm. An illustrative analysis of a typical launch vehicle establishes the speed of convergence, and accuracy and applicability of the algorithm.