A brief survey is made of the way the flow develops about wings with leading edges swept at about 50° as the stream Mach number rises from a subsonic to a supersonic value. The shock pattern which occurs may be complex. Three aspects of the wing flow are discussed in more detail, including the possible conditions for shock-induced separation of the surface boundary layer. The effect of the changing flow pattern on the overall wing lift, drag and pitching moment is commented upon.

In choosing the title, “The Development of All–Wing Aircraft,” as the subject of my lecture I run some risk of being accused of writing a company history rather than a paper of the broad scope ordinarily presented before this time–honoured institution. This is far from my intent, but being sincerely convinced that the all–wing aeroplane is a valuable step in the development of aeronautics, and desiring to contribute a maximum amount to the available data in the limited time at my disposal, my paper must be confined largely to experience gained by our company in its work on this subject.

The GOAHEAD (Generation of an Advanced Helicopter Experimental Aerodynamic Database for CFD code validation) consortium was created in the frame of an EU-project in order to create an experimental database for the validation of 3D-CFD and comprehensive aeromechanics methods for the prediction of unsteady viscous flows. This included the rotor dynamics for complete helicopter configurations, i.e. main rotor – fuselage – tail rotor configurations with emphasis on viscous phenomena like flow separation and transition from laminar to turbulent flow. The wind tunnel experiments have been performed during two weeks in the DNW-LLF on a Mach-scaled model of a modern transport helicopter consisting of the main rotor, the fuselage, control surfaces and the tail rotor. For the sake of controlled boundary conditions for later CFD validation, a closed test section has been used. The measurement comprised global forces of the main rotor and the fuselage, steady and unsteady pressures, transition positions, stream lines, position of flow separation, velocity profiles at the test section inlet, velocity fields in the model wake, vortex trajectories and elastic deformations of the main and tail rotor blades.

Most service failures in fatigue start on a free surface, as in many cases the most severe stress conditions exist there and are of a bi-axial stress nature, such as a drive shaft under combined bending and twisting. This particular case of multi-axial stress fatigue, because of its importance in design, has seen considerable investigation, although several aspects of the stress situation complicate the analysis of the test results with regard to the development of a general criterion of fatigue failure under multi-axial stress. The principal stress axes rotate with variations in the ratio of bending to twisting and the shaft is subjected to a stress gradient. Only the extreme portions of the shaft in the plane of bending are under the maximum bending stress, whereas the complete surface of the shaft is subjected to the maximum torsion stress.

The analysis of axisymmetrical thin shells has always been an important part of the theory of structures, and recently the problem of the cooling tower has aroused considerable interest for practical engineering reasons. Various solutions have been obtained, using either the finite difference or the finite element methods, with varying degrees of success. Amongst the latter, elements such as the TRIB 3c from the ASKA library and the axisymmetrical elements from the SABOR programme, with either a straight or a curved meridian line, have given reasonable displacements. But, by the nature of these elements, the stress distributions are discontinuous across element boundaries. Other existing elements, e.g. the arbitrarily curved triangular shell element SHEBA developed by Professor Argyris's team at the ISD in Stuttgart, are also suitable for the purpose and yield an excellent estimation of the stresses.

The lecture recalls great Manchester traditions in the mechanics of fluids and in related areas of mathematics which were brilliantly renewed from 1945 onwards through Sydney Goldstein’s work as Beyer Professor of Applied Mathematics. Goldstein brought to this important task not only his own personal distinction in those fields but also his great experience in leading interacting groups working in experimental and theoretical fluid dynamics with largely aeronautical objectives. The University has continued strongly to maintain these traditions.

The lecture celebrates also two fine pieces of research where gifted biologists have successfully used certain key mathematical ideas and methods in the mechanics of fluids to solve problems of importance in the field of animal locomotion. The first example is the work of C. J. Pennycuick, J. M. V. Rayner and G. R. Spedding in determining the character of the vortex wake behind a steadily moving bird in horizontal flapping flight. Typically, this takes the form of a “concertina wake” of fluctuating width but approximately constant circulation, exerting forces which are greater in the downstroke (which sheds a wide wake) and less in the upstroke (which sheds a narrower wake). Some notable advantages of the concertina wake are indicated.

Finally, work by Sir Eric Denton, Sir John Gray and Dr J. H. S. Blaxter on schooling behaviour in herrings and other clupeoid fishes is described. These fishes’ mechanisms for advantageously maintaining well coordinated motion of the school as a whole depend not only on eye vision but also on certain specialised organs admirably fitted to the sensing of hydrodynamic pressure fields (associated with the motions of other fishes) and of their lateral and longitudinal gradients.

In this work, stress-based and fracture mechanics criteria were developed to predict initiation and evolution, respectively, of intra- and inter-laminar cracking developed in composite laminates subjected to low velocity impact. The Soutis shear stress-strain semi-empirical model was used to describe the nonlinear shear behaviour of the composite. The damage model was implemented in the finite element (FE) code (Abaqus/Explicit) by a user-defined material subroutine (VUMAT). Delamination (or inter-laminar cracking) was modelled using interface cohesive elements and the splitting and transverse matrix cracks that appeared within individual plies were also simulated by inserting cohesive elements between neighbouring elements parallel to the fibre direction in each single layer. A good agreement was obtained when compared the numerically predicted results to experimentally obtained curves of impact force and absorbed energy versus time. A non-destructive technique (NDT), penetrant enhanced X-ray radiography, was used to observe the various damage mechanisms induced by impact. It has been shown that the proposed damage model can successfully capture the internal damage pattern and the extent to which it was developed in these carbon fibre/epoxy composite laminates.

The problem of designing a three-dimensional contraction arose recently in the course of the projected construction of a new low-speed (200 ft. /sec), closed-return wind-tunnel at the Northampton College of Advanced Technology. A three-dimensional contraction was required of contraction ratio 7·4:1 to provide smooth low turbulence conversion to the working section speed within an axial length of 126 inches. The inlet and outlet (working-section) dimensions are shown in Fig. 1 and the sections are similar. It was decided to design the contraction on the basis of an equivalent axi-symmetric contraction whose inlet and outlet circular areas are equal to those of the corresponding irregular octagons and then adapt the derived contraction wall profile to the actual sectional contours. The resulting equivalent inlet and outlet radii are then 55·4 in. and 20·4 in. respectively. The problem was thus reduced to the design of an axi-symmetric contraction of given area ratio within a specified axial length, with favourable velocity characteristics for low turbulence flow in the working section. This latter specification implies the absence of adverse pressure gradients on the walls near the working section and a small percentual difference between axial and wall velocities at the inlet to the working section.

A non-linear mathematical model for the study of the dynamics of an extensible cable subjected to aerodynamic forces generated by a uniform flow field is formulated. Solutions are found considering large displacement caused by suddenly applied loads (i.e., gusts, shock waves, turbulence) for a range of flow speeds and cable lengths. Transition from overdamped to oscillatory motion is observed when flow speed and cable length are increased and decreased respectively. The stability of the system is discussed.

Coupling an unsteady flow-solver with a structural model offers the opportunity to simulate aeroelastic behaviour of wings and rotor blades. The moving and deforming surfaces resulting from unsteady simulations require deforming meshes during the simulation, and it is common to use simple interpolation of surface displacements and velocities onto the initial undisturbed mesh. However, aeroelastic simulations can result in large displacements and deformations of solid surfaces, and simple interpolation of perturbations results in poor grid quality and possible grid crossover. A new interpolation technique is presented which is still simple in that it is driven solely by surface motion, but represents rotational effects near the solid surface, to maintain grid quality there. Furthermore, the scheme is fully analytic, so is very cheap computationally and results in grid speeds also being available analytically. Results, in terms of unsteady grid motion and flow solution, show the scheme to be effective and efficient.

In this study, development of a new fuel flow rate model for the climbing phase of flight was achieved using a genetic algorithm (GA) method. Two modelling approaches were performed using real flight data records (FDRs) from a medium-weight transport-category aircraft. The first model considered the dependency of fuel consumption only with respect to altitude, whereas the effects of both altitude and true airspeed (TAS) were included in the second model. The proposed models are improvements on existing models because the relationship between fuel flow rate, flight altitude, and TAS can be deduced using the derived formulations. Both modelling approaches were found to provide accurate results after performing an error analysis for fuel flow rate values. It was clear that incorporating the TAS effect into the second model enhanced the accuracy of the model, but the first model was also found to be appropriate for practical usage.

If surface fences are to be applied for measuring skin friction in three-dimensional boundary layers they must be calibrated for both magnitude and direction of the shear stress. Results of the calibration for fences of different height are given. Furthermore, a manufacturing process and a mounting procedure are described to obtain surface fences with identical calibration curves.

The paper provides a simplified picture of the processes occurring in liquid bi-propellant rocket motors. Droplet vaporisation, chemical reaction, and drag between gas and droplet are considered for a one-dimensional model. Quantitative theoretical results are presented for the variation of droplet size, droplet velocity, gas temperature and gas velocity within the chamber, particular attention being paid to the calculation of L*. The theory is applied to the German V2 rocket motor. Practical conclusions from, and extensions to, the theory are discussed.

In a gas turbine, ingestion of hot gas into the high-pressure turbine disc cavities could cause metal overheat. To prevent this, cool air is taken from the compressor and ejected through the cavities. However, this sealing flow also reduces the overall efficiency, and a compromise has to be found between the level of ingestion tolerated and the losses. Recent advances made in applying Computational Fluid Dynamics to such configurations are presented, with the aim of better understanding the physical phenomena and providing reliable design tools. First, results showing the pumping effect of the rotating disc are presented, including the influence of flow instabilities observed in both computational and experimental results. Second, the influence of the main annulus pressure asymmetries are analysed on a simplified representation of an available experiment, showing the combined influence of asymmetries generated by vanes and struts. Finally, a rim seal geometry representative of aero-engine design is studied in comparison to experiment, exhibiting the coupled influence of the cavity instabilities and annulus asymmetries.

A review of some methods of obtaining augmented lift coefficients in aerofoil sections precedes a description of lift by circulation control by blowing. The basic aerodynamic characteristics are considered as a prelude to a comparison of the merits of their application to helicopter rotors.

In calculating the performance of a helicopter rotor which uses high lift coefficients, the determination of the induced power has demanded an extension of the performance estimation curves normally used. A theoretical model has been developed to calculate the induced power factor which suggests a large increase in this parameter for rotor thrust coefficients above the current helicopter operating regime. Experimental results up to thrust coefficients of 0·1 which confirm the theoretical results are presented.

Some tentative ideas on the application of circulation control to produce a high speed helicopter are considered.

The problem of the vibrations induced in an aircraft structure excited by jet noise has been solved mathematically in general terms by Powell and others. The most promising methods for dealing with a practical structure are those based on the normal mode approach. In these methods the response of the structure in each of its normal modes is first determined and then the total response is obtained by a summation of individual responses. Use of these methods on practical structures is tedious because of the large number of modes which must be considered and also the difficulty of calculating theoretically the normal modes of a practical stiffened structure.

Adding a small aerodisk at the tip of the spike was proved to improve the spike’s performance in reducing the drag. If designed properly, the aerodisk spike can offer better performance in drag reduction in comparison with the pointed sharp-nosed spike. However, the impact of the aerodisk on the aeroheating reduction capabilities of the spike was not investigated in detail before. In the present work, a numerical investigation of the aerothermodynamic characteristics of spiked blunt bodies at Mach 6 flight conditions is conducted. It was found that, adding an aerodisk degrades the performance of the spike in reducing the aeroheating level. Nonetheless, the latter remains below that of the unspiked body.

In this paper, a mechanism of compressive fracture for elastic and elastic-plastic composite materials with interfacial adhesion defects is investigated. A classification of different approaches in modelling compressive response of layered materials is given. The analysis finds the upper and the lower bounds for the critical load. In order to achieve this, the problem of the internal fibre (layer) instability is considered within the scope of the exact statement based on the application of the model of a piecewise-homogeneous medium and the equations of the three-dimensional (3D) stability theory. The solution of the 3D problem is found for the most general case accounting for the bi-axiality of compressive loads. The characteristic determinants are derived for the first four fibre instability (microbuckling) modes, which are more commonly observed. Special attention is given to the calculation of critical loads for practical elastic and elastic-plastic layered materials.

Longitudinal vortices have tremendous practical utility for flow control (control by vortices) and in some cases have exhibited tremendous potential for causing harm if uncontrolled (i.e. control of vortices is required). Vortex control has thus far been carried out via multitudinous approaches in an empirical fashion, aided by the essentially inviscid nature of much of longitudinal vortex behaviour. Further refinement and several applications of vortex flow control require knowledge regarding the detailed flow physics of longitudinal vortices such as transition, transitional flow regimes, turbulence structure and modelling, and interaction with shock waves, other vortices and surfaces. This paper summarises vortex control applications and extant techniques for the control of longitudinal vortices produced by bodies, leading edges, tips, and intersections.