This paper investigates the hydrodynamic near-field of a NACA 16-616 aerofoil over a range of angles-of-attack, encompassing the pre-stall, stall and post-stall flow regimes. In both the static pressure and the pressure fluctuation results, it is shown that each flow regime is easily distinguished, and it is further shown that each regime has different spectral behaviour and boundary layer characteristics. It is found that the NACA 16-616 aerofoil stalls by an abrupt leading-edge mechanism, characterised by a sudden change in the static pressure and unsteady surface pressure spectra between $16^\circ $ and $17^\circ $ angles-of-attack, but of more interest is that there is a secondary yet significant trailing-edge flow separation mechanism occurring upstream of the trailing-edge and moving further upstream as the angle-of-attack increases in the pre-stall regime. A comparison is made between the spectra and coherence of the unsteady surface pressure of the NACA 16-616 aerofoil and those of the classic NACA 0012 aerofoil and shows that such a secondary mechanism has a significant impact for large pre-stall angles-of-attack on the unsteady surface pressure. This will have a significant impact on the radiated far-field sound, distinguishing the NACA 16-616 aerofoil from aerofoils that do not have this secondary mechanism. The existence and extent of this secondary trailing-edge separation mechanism is further shown by the hot-wire anemometry boundary layer velocity results that indicate separation within the pre-stall regime.

This paper describes a new efficient method for the construction of an approximately balanced aerodynamic Reduced Order Model (ROM) via the frequency domain using Computational Fluid Dynamics data. Time domain ROM construction requires CFD data, which is obtained from the DLR TAU RANS or Euler Linearised Frequency Domain (LFD) solver. The ROMs produced with this approach, using a small number of frequency simulations, are shown to exhibit a strong ability to reconstruct the system response for inviscid flow about the NLR7301 aerofoil and the FFAST wing; and viscous flow about the NASA Common Research Model. The latter demonstrates that the reduced order model approach can reconstruct the full order frequency response of a viscous aircraft configuration with excellent accuracy using a strip wise approach. The time domain models are built using the frequency domain, but also give promising results when applied to reconstruct non-periodic motions. Results are compared to time domain simulations, showing good agreement even with small ROM sizes, but with a substantial reduction in calculation time. The main advantage of the current model order reduction approach is that the method does not require the formation and storage of large matrices, such as in POD approaches.

A spacetime formulation is presented to solve unsteady aerodynamic problems involving large deformation or topological change such as store separation, slat and flap deployment or spoiler deflection. This technique avoids complex CFD meshing methods, such as Chimera, by the use of a finite-volume approach both in space and time, and permits a locally varying real timestep. The use of a central-difference scheme in the time direction can yield non-physical transient solutions as a consequence of information travelling backwards in time. Therefore, an upwind formulation is provided and validated against one-dimensional and two-dimensional test cases. A hybrid formulation (central in space, upwind in time) is also given and unsteady cases are computed for a spoiler and spoiler/flap deployment, with all three formulations compared, demonstrating that the use of an upwind time stencil yields more representative physical solutions and improves the rate of convergence.

In this paper, the accuracy and practical capabilities of three different reduced-order models (ROMs) are explored: an enhanced implicit condensation and expansion (EnICE) model, a finite element beam model, and a finite volume beam model are compared for their capability to accurately predict the nonlinear structural response of geometrically nonlinear built-up wing structures. This work briefly outlines the different order reduction methods, highlighting the associated assumptions and computational effort. The ROMs are then used to calculate the wing deflection for different representative load cases and these results are compared with the global finite element model (GFEM) predictions when possible. Overall, the ROMs are found to be able to capture the nonlinear GFEM behaviour accurately, but differences are noticed at very large displacements and rotations due to local geometrical effects.

Elevator failure may have fatal consequences for fighter aircraft that are unstable due to their high manoeuvrability requirements. Many studies have been conducted in the literature using active and passive fault-tolerant control structures. However, these studies mostly include sophisticated controllers with high computational load that cannot work in real systems. Considering the multi-functionality and broad operational prospects of fighter aircraft, computational load is very important in terms of applicability. In this study, an integrated fault-tolerant control strategy with low computational load is proposed without sacrificing the ability to cope with failures. This control strategy switches between predetermined controllers in the case of failure. One of these controllers is designed to operate in a non-failure condition. This controller is a basic controller that requires very little computational effort. The other controller operates when an asymmetric elevator failure occurs. This controller is a robust fault-tolerant controller that can fly the aircraft safely in case of elevator failure. The switching is decided by a failure detection system. The proposed integrated fault-tolerant control system is verified by non-linear F-16 flight simulations. These simulations show that the proposed method can cope with failures but requires less computational load because it uses a conventional controller in the case of no failure.

Landing gear are exposed to cyclic loads from the ground manoeuvres that aircraft perform in-service. Variability is observed in the loading magnitude associated with ground manoeuvres, along with the per-flight variability in ground manoeuvre occurrence and sequencing. Whilst loading magnitude variability has been widely characterised, significant assumptions are required regarding manoeuvre occurrence and sequencing when constructing landing gear load spectra for fatigue design. These assumptions are required due to the limited availability of data concerning ground manoeuvre occurrence and sequencing relating to aircraft in-service and require validation to facilitate the design of more efficient components. ‘Big-Data’ approaches, employing Automatic Dependent Surveillance-Broadcast (ADS-B) transponder data, enable aircraft ground tracks to be identified. This paper presents a methodology to characterise the variability in ground manoeuvre occurrence and sequencing using ADS-B data sourced from Flightradar24® for a wide-body aircraft fleet. Using statistics generated for the fleet, it was identified that significant variability exists in the occurrence and sequencing of turning and braking manoeuvres. The statistics also validate existing assumptions, including that the proportional share of left and right turning manoeuvres is equal. Finally, this paper discusses the utility of ADS-B datasets for constructing landing gear load spectra and monitoring of landing gear in-service.

This article presents a development of a simple analytical aerodynamic model capable of describing the effect of leading-edge vortices (LEVs) on the lift of rotating samara wings. This analytical model is based on the adaptation of Polhamus’ method to develop a sectional two-dimensional lift function, which was implemented in a numerical blade element model (BEM) of a rotating samara blade. Furthermore, wind tunnel experiments were conducted to validate the numerical BEM and to assess the validity of the newly developed analytical lift function. The results showed good agreement between the numerical model and the experimental measurements of rotational speed and rate of descent of the samara wing. The results were also compared with numerical predictions using BEM but adopting different lift coefficient expressions available in literature. This research contributed towards efficient aerodynamic modelling of the lift generated by LEVs on rotating samara wings for performance prediction, which could potentially be used in the design of bio-inspired rotary micro-air vehicles.

This paper investigates the synergies and trade-offs between passive aeroelastic tailoring and adaptive aeroelastic deformation of a transport composite wing for fuel burn minimisation. This goal is achieved by optimising thickness and stiffness distributions of constitutive laminates, jig-twist shape and distributed control surface deflections through different segments of a nominal “cruise-climb” mission. Enhanced aerostructural efficiency is sought both passively and adaptively as a means of aerodynamic load redistribution, which, in turn, is used for manoeuvre load relief and minimum drag dissipation. Passive shape adaptation is obtained by embedding shear-extension and bend-twist couplings in the laminated wing skins. Adaptive camber changes are provided via full-span trailing-edge flaps. Optimised design solutions are found using a bi-level approach that integrates gradient-based and particle swarm optimisations in order to tailor structural properties at rib-bay level and retrieve blended stacking sequences. Performance benefits from the combination of passive aeroelastic tailoring with adaptive control devices are benchmarked in terms of fuel burn and a payload-range efficiency. It is shown that the aeroservoelastically tailored composite design allows for significant weight and fuel burn improvements when compared to a similar all-metallic wing. Additionally, the trailing-edge flap augmentation can extend the aircraft performance envelope and improve the overall cruise span efficiency to nearly optimal lift distributions.

A recent consideration in aircraft design is the use of folding wing-tips with the aim of enabling higher aspect ratio aircraft with less induced drag while also meeting airport gate limitations. This study investigates the effect of exploiting folding wing-tips in flight as a device to reduce both static and dynamic loads. A representative civil jet aircraft aeroelastic model was used to explore the effect of introducing a wing-tip device, connected to the wings with an elastic hinge, on the load behaviour. For the dynamic cases, vertical discrete gusts and continuous turbulence were considered. The effects of hinge orientation, stiffness, damping and wing-tip weight on the static and dynamic response were investigated. It was found that significant reductions in both the static and dynamic loads were possible. For the case considered, a 25% increase in span using folding wing-tips resulted in almost no increase in loads.

By means of a broad survey of the developments in aeroelastic science during the past ten years and of comparable developments in associated fields, in so far as they involve elastic deformation, the present paper attempts to show how the individual subjects concerned are developing strong interconnections. At one time the subjects were for the most part regarded as discrete entities; now they are developing into component parts of an integrated whole —the dynamics of a deformable aeroplane.

It is suggested that in the future it may be necessary to discontinue the present methods of treatment of the component subjects as separate, and to adopt a standard method of attack on the problem treated as a unified whole, especially when adequate computational aids become available. A modified approach, involving an arrangement of the subjects in a frequency spectrum, is also discussed. Some speculations on the elastic provisions which may be necessary in the future are included.

An implicit cell-vertex upwind method is presented for the solution of the two-dimensional Euler equations for supersonic flows. It is based on the two-point difference scheme introduced by S.F. Wornom, but contains many significant developments. Particular attention is given to the implementation of boundary conditions, to ensure that no “numerical” boundary conditions are required, and “image” cells are avoided, making the scheme ideal for use in conjunction with a multiblock approach. Examples are given for quasi-one-dimensional nozzle flows, and in two dimensions, oblique shock reflection, a shock interaction problem, and a supersonic intake problem. The accuracy and efficiency of the scheme is demonstrated, along with its ease of use.

The precise birth date of the aerospace composites industry cannot readily be identified; perhaps one should really talk about its rebirth as the first aircraft relied on natural composites such as wood. McMullen gives 1946 as the date that work on cellulose based composites for aircraft use was abandoned in favour of much more stable inorganic reinforcement fibres. This change in the direction of approach was crucial to further developments and can be thought of as marking the start of the aerospace composites industry that can be seen today. Whatever the exact date the industry is now about 50 years old so this golden jubilee edition seems an appropriate place to look at the constraints on the use of composite materials and at recent work at Bristol University aimed at reducing these constraints.

A review of recent developments in helicopter noise reduction is presented. The paper includes a discussion of the present state of understanding of helicopter noise radiation mechanisms and a review of actual noise reductions achieved by the manufacturers over the past decade. An analysis of progress achieved in meeting certification limits is given, together with an evaluation of the impact on design. Future trends in helicopter noise reduction based on current noise research programmes are discussed and suggestions advanced for possible additional areas of study.

I am very appreciative of the honour that has been accorded me in being invited to deliver the Second Lanchester Memorial Lecture: appreciative, but all too conscious of my own inadequacy. This is especially true, since in the first Lecture we heard a tribute to a very great man from another such: and the prospect of following von Karman to this lectern is one to deter perhaps all but the most foolhardy. Since von Karman’s address discussed very fully the contributions of Lanchester to aeronautics, little is left for subsequent lecturers on this subject, and I have therefore chosen to speak of a problem which has been my special interest and concern throughout my professional life, namely, aeroelasticity—a field in which, as in so many others, Lanchester was a pioneer.

The ‘inverse method’ proposed in 1936 by Duncan, Collar and Lyon for the calculation of flutter boundaries is re-stated in a form suited to modern methods of dynamical and numerical analysis and programmed for use on a home computer. The method is extended to provide information on sub- and supercritical eigenvalues. For a microcomputer with 32K RAM, flutter problems of orders up to 35 may be accurately solved.

This paper describes the application of the theory of pulse responses to a linearised Euler scheme on a moving grid. The impulse responses of such a linear system are memory functions or temporal representations of the manner in which, and the time over which a perturbation remains active in the response of the system. The exact response of the linear system to an arbitrary input can then be predicted via convolution. The linearised Euler scheme is derived from a full nonlinear Euler method based on Jameson’s cell-centred scheme which has been modified to make it time-accurate and has the necessary terms to account for grid motion. The discrete form of the full Euler equations are linearised about the corresponding nonlinear steady mean flow assuming the unsteadiness in the flow and grid are small. The resulting set of differential equations for the flow perturbations are solved at each time step using a dual-time scheme. These equations can be easily solved to obtain pulse responses, which capture the entire frequency content of the linear system. These pulse responses then provide a simple and efficient method for calculating the system response to many different inputs via convolution. A further use of the pulse responses is their potential to allow a reduced-order linear model to be constructed, which preserves some of the accuracyof the original system. Preliminary work carried out in this area is presented. Results of calculations for heave, pitch and ramp test cases are presented.

A study of the laminar boundary layer on a yawed infinite wing in a compressible fluid is made. A method is presented for the investigation of compressibility effects by direct solution of the linearised equations of motion, first considerably simplified by the extension of a transformation due to Illingworth and Stewartson, which assumes a heat-insulated surface and a fluid of unit Prandtl number and a coefficient of viscosity which is proportional to the absolute temperature. As an example of this procedure the boundary layer near a stagnation point is calculated. The simplifications necessary for an extension of the momentum method, which is available for a yawed infinite wing in an incompressible fluid, are discussed, and this method is recommended for a first evaluation of the effect of compressibility.

Corrugated structures offer a potential solution for morphing wing skin applications due to their anisotropic behaviour that allows chordwise camber and length changes. Aerofoils with corrugated skins in the aft 1/3 of the chordwise section have been studied experimentally and computationally using various corrugation shapes and forms (sinusoidal, trapezoidal and triangular) at different Reynolds numbers. The study showed that the aerodynamic performance is highly dependent on corrugation amplitude, wavelength, gradient (combination of amplitude and wavelength) and Reynolds number. Evidence is given highlighting that penalties for having a non-smooth surface in the aft 1/3 of the chordwise section of an aerofoil can be eliminated for the lift curve slope and minimised for the zero lift drag coefficient.

When designing flight control laws using linearisations of an aircraft model about different flight conditions, some form of scheduling of the resultant gains will often be required to implement the controller over wide operating regions. In practice, the controller gains are often scheduled against relatively slowly-varying system states such as altitude or velocity. However, it may also be desirable to schedule gains against rapidly-varying states such as angle-of-attack, thereby generating a cyclic dependence through hidden coupling terms. Previous published work at Bristol has developed a numerical method of accounting for this dependence when scheduling state feedback gains against coupled states. The resulting ‘dynamic gain schedule’ is shown to significantly improve the transient response of the aircraft model during rapid manoeuvring and to reduce the chances of control surface actuator position limit saturation. In this paper, the novel design process, using eigenstructure assignment, is applied to a mathematical second-order longitudinal aircraft model which represents an approximate BAe Hawk wind-tunnel model. The dynamic gain scheduled controller is shown to work extremely well in practice when applied to the closed-loop experimental rig. Despite the highly nonlinear characteristics of the model aerodynamics and tailplane actuation system, as well as unmodelled high turbulence levels, dynamic gain scheduling demonstrates stable closed loop control even in regions where the nonlinearities are such that conventional gain scheduling fails to produce a stable response.