The aerodynamics of insect-like flapping are dominated by the production of a large, stable, and lift-enhancing leading-edge vortex (LEV) above the wing. In this paper the phenomenology behind the LEV is explored, the reasons for its stability are investigated, and the effects on the LEV of changing Reynolds number or angle-of-attack are studied. A predominantly-computational method has been used, validated against both existing and new experimental data. It is concluded that the LEV is stable over the entire range of Reynolds numbers investigated here and that changes in angle-of-attack do not affect the LEV’s stability. The primary motivation of the current work is to ascertain whether insect-like flapping can be successfully ‘scaled up’ to produce a flapping-wing micro air vehicle (FMAV) and the results presented here suggest that this should be the case.

The aircraft industry, as a whole, is striving to limit its impact on the environment. Improved engine design and operation may offer a reduction in emissions of a few percent. More efficient air traffic control (ATC) may offer a limited reduction in overall fuel burn. Improvements in aerodynamic design and materials available (e.g. on A350XWB, B787) might achieve a few percent increases in efficiencies. The use of alternative fuels is some way off. The ACARE objectives present a stiff challenge.

Our recent studies have shown that air-to-air-refuelling (AAR), well established in military circles, introduced to civil aircraft operations would provide fuel savings of the order of 30% – 40%. AAR will allow smaller (3,000nm range), more efficient (greener) aircraft, operating from shorter runways, to fulfil long-range route requirements. In addition, the ‘safety-net’ afforded by the availability of AAR will enable a host of hitherto borderline technologies to be accepted and utilised in future aircraft designs. Laminar flow will provide fuel savings and increased efficiency in its own right provided it is enabled within a civil AAR environment. Similarly, supersonic transport becomes an acceptable economic option.

The interesting article by Archer, Sapuppo and Betteridge leads to the problem of determining the optimum lift distribution along the wing during flapping motion. Adopting the assumptions made in Ref. 1, the problem may be reduced to one of minimising the induced drag for a specified, periodically varying, bending moment at the wing root. The wing-root bending moment should be of primary importance to a bird or, for that matter, to a mechanical flapping-machine, since it is at this point that most of the flapping energy is introduced.

We discuss the status, trends and long-term ambitions of Computational Fluid Dynamics (CFD) when applied away from the design point or concept, and therefore in the historically weak areas of CFD. This includes both undesirable flight conditions, such as stall, and undesirable products of the flow, such as noise. All pose the great challenge of turbulence, and accuracy is as dependent on the ideas behind the turbulence treatment as it is on computing power. A measured shift from Reynolds-averaged representations to large-eddy simulations will take place. Empiricism, both turbulence and engineering related, will recede only step by step over many years. Yet, CFD will make full use of every increase in computer power of this century. Increasingly, CFD will compete with flight tests, not just with wind tunnels, and will be validated by flight tests. Integration with other disciplines will allow us to predict crucial phenomena such as flutter, sonic fatigue, and pilot-induced oscillations. A gap will remain at any time between the phenomena amenable to a ‘grand challenge’ calculation and those amenable to a fully CFD-based design process, because certification involves thousands of conditions. The highest demand currently is in community and cabin noise for which industrial-accuracy methods are non-existent. On the other hand, gratifying progress has occurred in the area of stall and spin.

This is a discussion of the various types of flow which may occur on thin swept wings when the air stream is not attached to the wing surface everywhere. In particular, the behaviour of vortex sheets which arise from such mixed flows is examined.

This third technical note, on the background of the author's Main Lecture to the Society(1), illustrates the application of the technique given in ref. 2 on a simple example, and develops the basic ideas further, with particular reference to the effect of reinforcing flanges and initial loads, as may arise through differential thermal straining. A more ambitious example demonstrating the improvement achieved by the new triangle appears in ref. 1.

Jets issuing from rectangular slots cut in thin plates exhibit some unusual features, including unequal spreading rates in the spanwise and transverse directions, the appearance of velocity peaks near the “ends” of the jet and changing rates of centreline velocity decay in the downstream direction. This study examines the effects of aspect ratio on such flows. The flow field has been investigated using both total head tubes and hot wire anemometry. The results are presented in the form of three-dimensional plots of total pressure and contour maps of constant velocity, streamwise turbulence intensity and the Reynolds shear stresses. The decay of mean velocity and stream-wise turbulence intensity along the centreline are presented. The rates of spanwise spreading and the location of the velocity peaks at various downstream stations are discussed. If the aspect ratio is small enough, spanwise peaks in the mean velocity distribution are suppressed.

Comparative mixing of a Mach 2 elliptical free jet from a convergent-divergent elliptic nozzle with an aspect ratio of 2:1 in the presence of adverse and marginally favourable pressure gradients has been studied experimentally. It is found that the mixing of the elliptical jet is higher than that for the equivalent circular jet at all the levels of expansion. The decay of elliptic jet is significantly higher than the equivalent circular jet in all three zones of the jet field – the core, characteristic decay and fully developed regions. The reason for the faster decay of the elliptic jet is found to be the continuous variation in the size of the mixing-promoting vortices shed from the nozzle exit owing to its azimuthal asymmetry. The evolution of the jet and its axis-switching phenomenon has been studied using iso-pitot pressure contours taken at different axial locations in the plane normal to the jet axis. As expected, the elliptic jet spreads faster along the minor axis plane than the major axis plane, leading to axis-switching at all the levels of expansion studied. The axis-switching of the elliptic jet shifts upstream with increase in nozzle pressure ratio (NPR) from 4 to 5; from 5 to 7, it shifts downstream. But at marginally under-expanded condition of NPR 8, the axis-switching is found to shift slightly upstream. The occurrence of axis-switching in the elliptic jet indicates enhanced near-field mixing, compared to the equivalent circular jet. The shadowgraph pictures of the jet reveal that the waves prevailing in the elliptic jet are significantly weaker than those in the circular jet.

The design, construction, and validation of a new academic wind tunnel is described in detail. The wind tunnel is of a classical, blow-down type and generates a pressure-induced, turbulent separation bubble on a flat test surface by a combination of adverse and favorable pressure gradients. The Reynolds number, based on momentum thickness just upstream of separation, is Reθ ≃ 5,000 at a free-stream velocity of Uref = 25ms−1. The length of the separation bubble is estimated at 0°42 ± 0°02m by three different methods. Results of a numerical simulation demonstrate the absence of flow separation in the wind-tunnel contraction. This results in a turbulence level of about 0·05% in the test section. Oil-film visualisation experiments show that the flow near the wall is strongly three-dimensional in the recirculating region and that the topology of the limiting streamlines is consistent with experiments performed on configurations with fixed separation. Finally, spatial variations of the forward-flow fraction have been documented using a thermal-tuft probe and are shown to compare well with the results of the oil-film visualisation.

An improved variational optimization approach is established to optimize and analyse the propulsion efficiency of the high-altitude contra-rotating propellers for high-altitude airships based on the Vortex Lattice Lifting Line Method. The optimum radial circulation distribution, chord and pitch distribution are optimized under the maximum lift-to-drag ratio of aerofoils. To consider the effects of the actual Reynolds number and the Mach number of each aerofoil section, aerodynamics such as lift coefficient, drag coefficient and lift-to-ratio are obtained by interpolating a CFD database, which is established by numerical simulations under different Reynolds number, Mach number and angles-of-attack. The improved method is verified by validation cases on a high-altitude CRP using the three-dimensional steady Reynolds-averaged Navier-Stokes solver and moving reference frames technique. The optimization results of thrust, torque and efficiency for both the individual front/rear propeller and CRP are shown to agree reasonably well with the CFD results. Using the improved approach, the influence of blade numbers, diameter, rotation speeds, axial distance and torque ratio on the optimum efficiency of CRPs is illustrated in detail by conducting parametric studies.

Recently Messiter has proposed a first-order correction to simple Newtonian theory for the pressure distribution on the lower (compression) surfaces of lifting conical bodies. Although the basic theory holds for bodies with and without attached shock waves, solutions have so far only been obtained for bodies with detached shocks. In the present paper an approximate method of applying the theory to bodies with attached shocks is given. In spite of the approximations involved the calculated shock shapes and pressure distributions are in good agreement with some exact solutions for flat wings, except near the incidence for shock detachment. Like the detached shock case, the present solution can be applied to Nonweiler wings in certain off-design conditions. The combined results for the detached shock and for the attached shock enable the off-design behaviour of Nonweiler wings to be discussed in a systematic manner.

The purpose of a joint is to transfer load between the two parts being joined. As a result of this load transfer there will be a stress variation in the components in the joint region, as well as stresses in the joining medium (fasteners or adhesive). In fhjs work, the stress results of a two-dimensional finite element (FE) analysis are used to understand failure modes of a bolted joint in low-temperature cure carbon fibre reinforced plastic (CFRP) woven laminates loaded in tension, and to predict the bearing strength. Good agreement with experimentally observed damage modes and strengths is achieved with, in some cases, a difference of less than 10%.

Within the framework of a major European co-operative research project on rotorcraft aerodynamics and acoustics (Helishape), a parametric model rotor test was conducted in the open test section of the DNW using DLR's MWM test rig, and a highly instrumented model of a fully articulated ECF rotor equipped with blades of advanced design with two exchangeable blade tips. One set of blade tips (7A) was rectangular, while the other (7AD1) was a swept-back parabolic-anhedral shape. The objectives of this experimental research were to evaluate noise reduction techniques (conceptionally by variation of rotor speed, dedicated tip shapes and advanced aerofoils, as well as operationally, by identifying low noise, BVI-minimising, descent procedures) and to validate the partners’ aerodynamic and aeroacoustic codes. A comprehensive set of simultaneous acoustic and aerodynamic blade surface pressure data and blade dynamic and performance data was measured. In addition, valuable information on the tip-vortex geometry and blade-vortex miss distance was obtained by LLS flow visualisation. The experimental equipment, the test procedures, and the test matrix are briefly described. A survey of the main results is presented and the trends of the most important parameter variations for both rotors are discussed.

In this paper we study large angle rotational manoeuvres of spacecraft. The problem of attitude control is formulated and solved through potential functions, thus simplifying frequent reorientation manoeuvres. A combination of gas jet and FEEP thrusters is employed in order to obtain short settling times while meeting stringent pointing requirements. In this way, spacecraft reorientation is achieved in two integrated steps by means of a discrete/continuous control law. Autonomous avoidance of undesired space orientations is obtained and constraints due to the solar array pointing requirements are satisfied. A case study is discussed where the methodology is applied to a large space telescope.

This paper describes the development and validation of a high fidelity simulation model of the Bell 412 helicopter for handling qualities and flight control investigations. The base-line model features a rigid, articulated blade-element formulation of the main rotor, with flap and lag degrees of freedom. The Bell 412 HP engine/governor dynamics are represented by a second-order system. Other key features of the base-line model include a finite-state dynamic inflow model and lag damper dynamics. The base-line model gives excellent agreement with flight-test data over the speed range 15-120kt for on-axis responses. Prediction of off-axis responses is less accurate. Several model enhancement options were introduced to obtain an improved off-axis response. It is shown that the pitch/roll off-axis responses in transient manoeuvres can be improved significantly by including wake geometry distortion effects in the Peters-He finite-state dynamic inflow model.

Piezoelectric wafer active sensors (PWAS) are lightweight and inexpensive transducers that enable a large class of structural health monitoring (SHM) applications such as: (a) embedded guided wave ultrasonics, i.e., pitch-catch, pulse-echo, phased arrays; (b) high-frequency modal sensing, i.e., electro-mechanical impedance method; and (c) passive detection. The focus of this paper is on the challenges posed by using PWAS transducers in the composite laminate structures as different from the metallic structures on which this methodology was initially developed. After a brief introduction, the paper reviews the PWAS-based SHM principles. It follows with a discussion of guided wave propagation in composites and PWAS tuning effects. Then, the mechanical effect is discussed on the integration of piezoelectric wafer inside the laminate using a compression after impact. Experiments were performed on a glass fibre laminate, employing PWAS to measure the attenuation coefficient. Finally, the paper presents some experimental and multi-physics finite element method (MP-FEM) results on guided wave propagation in composite laminate specimens.

This paper presents a new flight trajectory optimisation method, based on genetic algorithms, where the selected optimisation criterion is the minimisation of the total cost. The candidate flight trajectories evaluated in the optimisation process are defined as flight plans with two components: a lateral flight plan (the set of geographic points that define the flight trajectory track segments) and a vertical flight plan (the set of data that define the altitude and speed profiles, as well as the points where the altitude and/or speed changes occur). The lateral components of the candidate flight plans are constructed by selecting a set of adjacent nodes from a routing grid. The routing grid nodes are generated based on the orthodromic route between the flight trajectory’s initial and final points, a selected maximum lateral deviation from the orthodromic route and a selected grid node step size along and across the orthodromic route. Two strategies are investigated to handle invalid flight plans (relative to the aircraft’s flight envelope) and to compute their flight performance parameters. A first strategy is to assign a large penalty total cost to invalid flight profiles. The second strategy is to adjust the invalid flight plan parameters (altitude and/or speed) to the nearest limit of the flight envelope, with priority being given to maintaining the planned altitude. The tests performed in this study show that the second strategy is computationally expensive (requiring more than twice the execution time relative to the first strategy) and yields less optimal solutions. The performance of the optimal profiles identified by the proposed optimisation method, using the two strategies regarding invalid flight profile performance evaluation, were compared with the performance data of a reference flight profile, using identical input data: initial aircraft weight, initial and final aircraft geographic positions, altitudes and speed, cost index, and atmospheric data. The initial and final aircraft geographic positions, and the reference flight profile data, were retrieved from the FlightAware web site. This data corresponds to a real flight performed with the aircraft model used in this study. Tests were performed for six Cost Index values. Given the randomness of the genetic algorithms, the convergence to a global optimal solution is not guaranteed (the solution may be non-optimal or a local optima). For a better evaluation of the performance of the proposed method, ten test runs were performed for each Cost Index value. The total cost reduction for the optimal flight plans obtained using the proposed method, relative to the reference flight plan, was between 0.822% and 3.042% for the cases when the invalid flight profiles were corrected, and between 1.598% and 3.97% for the cases where the invalid profiles were assigned a penalty total cost.

It is well known that the main purpose of the by-pass principle is to improve the propulsive efficiency of a simple jet engine by removing some of the energy left in the jet gases and using this to compress an extra quantity of air, known as the by-pass air, this air being ejected rearwards with the jet gases. In this way a greater mass of air is ejected rearwards at a lower jet velocity and thus a better propulsive efficiency is obtained. This is an extremely simplified view of the advantages of the by-pass engine, however, since an equivalent result of obtaining a lower jet velocity can be obtained by designing the jet engine for a lower combustion temperature. The by-pass principle is of advantage because it enables a higher propulsive efficiency to be obtained at the same time as employing a high combustion temperature and therefore a high basic cycle efficiency. If the component efficiencies of a gas turbine were 100 per cent, cycle efficiency would not depend upon combustion temperature at all, and there would thus be no advantage in principle in using the by-pass engine. In practice there would probably be some residual advantage left in that for a given thrust a lower engine weight could be obtained.

The inflation of a parachute encompasses problems of aerodynamics, dynamics, and elasticity. The most noticeable effects of the interaction of the events are the opening force and the filling time. Attempts to analyse the opening process and to predict the opening force date back to 1927. Since then a considerable number of methods have been published. Some of them are primarily analytical, others more empirical, and several deal with isolated problems. The earliest and most analytically oriented studies are by Scheube and O'Hara.

Computational aerodynamics, which complement more expensive empirical approaches, are critical for developing aerospace vehicles. During the past three decades, computational aerodynamics capability has improved remarkably, following advances in computer hardware and algorithm development. However, most of the fundamental computational capability realised in recent applications is derived from earlier advances, where specific gaps in solution procedures have been addressed only incrementally. The present article presents our view of the state of the art in computational aerodynamics and assessment of the issues that drive future aerodynamics and aerospace vehicle development. Requisite capabilities for perceived future needs are discussed, and associated grand challenge problems are presented.