The establishment of the UK Applied Aerodynamics Consortium in 2004 brought together many of the UK’s leading research groups to tackle challenging aerodynamic problems on the national computing facility, HPCx. This paper provides a brief history of some early pioneers of numerical simulation and highlights some key contributions to development in parallel processing that laid the foundations for today’s researchers. The transition from vector to massively parallel processing is discussed from a UK viewpoint along with technological barriers that could have a significant impact on future systems. Solutions to these barriers are already being sought and the paper discussed some of the novel technologies that may be deployed in the future. In its short history, the consortium has made substantial progress and this is briefly discussed with several highlights that illustrate the scientific output. Although a number of challenges are identified, particularly with respect to developing a comprehensive visualisation capability, the consortium is well placed to build upon its initial success.

Hot gas ingestion (HGI) can be a problematic feature of short take-off vertical landing (STOVL) aircraft during the descent phase of landing, or while on the ground. The hot exhaust gases from the downwards pointing nozzles can be re-ingested into the engine intakes, causing power degradation or reduced engine surge margin. The flow-fields that characterise this phenomenon are complex, with supersonic impinging jets and cross-flows creating large ground vortices and fountain up-wash flows.

A flow solver has been developed to include a suitable linear mesh deformation technique for the descending aircraft configuration. The code has been applied to predict the occurrence of HGI, by simulating experimental results from a 1/15th scale model of a descending Harrier. This has enabled an understanding of the aerodynamic mechanisms that govern HGI, in terms of the near-field and far-field effects and their impact on the magnitude of temperatures at the engine intake.

This paper presents three sets of CFD results. First a validation exercise shows predicted results from the twin-jet with intake in crossflow test-case. This is an unsteady Reynolds averaged Navier Stokes (URANS) solution for a static geometry (there is no moving mesh). This allows comparison with experiment. Secondly, a full descent phase URANS Spalart-Allmaras (SA) turbulence model calculation is done on an 8·5m cell mesh for half the flow domain of the Harrier model and test-rig without dams/strakes. This shows how the HGI flow mechanisms affect the engine intake temperature profiles, for the case where there are no flow control methods on the underside of the aircraft. Thirdly, the full descent phase URANS SA turbulence model calculation is done on a 22·4m cell mesh for the full flow domain of the Harrier model and test-rig, with the dam/strake geometry included in the structured mesh region.

Aspects of parallel simulation of rotor flows are considered. These flows can be extremely expensive for a compressible finite-volume CFD code, and parallelisation can be essential. The award of HPCx time through the UK Applied Aerodynamics Consortium has allowed large rotor simulations to be performed and wake grid dependence to be investigated. However, there are several issues that need to be investigated when considering very large simulations, including the grid generation process, the parallel flow-solver, including an effective mesh motion approach, and visualisation options. Details of these are presented here, with particular emphasis on the flow-solver parallel performance. A detailed performance analysis of the unsteady flow-solver has been undertaken and the code optimised to improve parallel performance, and details of the parallel scaling performance are presented. The parallel scaling of the code is very good on all the HPC architectures tested here, and this has been recognised by an HPCx Gold Star Capability Incentive award. Results of simulation of a fourbladed lifting rotor in forward flight are also presented, for two mesh densities. It is shown that the solution computed on the serial limit on mesh size, around four million cells, exhibits excessive diffusion, and is of limited use in terms of detailed flow features. The results on a very fine mesh, 32 million cells, have shown a much better solution resolution, and it is also demonstrated that the λ2 vortex core visualisation option is extremely useful.

Numerical analysis of the flow in weapon bays modelled as open rectangular cavities of length-to-depth (L/D) ratio of 5 and width-to-depth (W/D) ratio of 1 with doors-on and doors-off is presented. Flow conditions correspond to Mach and Reynolds numbers (based on cavity length) of 0·85 and 6·783m respectively. Results from unsteady Reynolds-averaged Navier-Stokes (URANS), large-eddy simulation (LES) and detached-eddy simulation (DES) are compared with the simulation methods demonstrating the best prediction of this complex flow. It was found that URANS was not able to predict the change of flow characteristics between the doors-on and doors-off configurations. In addition, the energy content of the cavity flow modes was much better resolved with DES and LES. Further, the DES was found to be quite capable for this problem giving accurate results (within 3dB of) experiments and appears to be a promising alternative to LES for modelling massively separated flows.

An Euler optimisation for a BWB configuration with winglets incorporating an array of three-dimensional shock control bumps is carried out by employing an efficient adjoint-based optimisation methodology. A high fidelity multi-block grid with over two million grid points is generated to resolve the shape of the 3D shock control bumps, the winglet as well as the overall BWB shape, which are parameterised by over 650 design variables. In order to perform such a large aerodynamic optimisation problem feasibly, the optimisation tools such as the flow solver and the adjoint solver have to be parallelised with a good parallel efficiency. This paper reports the parallel implementation efforts on the adjoint solver; especially on the calculation of the sensitivity derivatives, which has to be looped over the total number of design variables. Results from the optimisation of the wing master sections, winglet aerofoil sections and the three dimensional bumps indicate a significant improvement regarding the aerodynamic performance against the baseline geometry for the given planform layout of the aircraft.

This paper presents a critical analysis of the aerodynamic loads created by the airframe of a conventional helicopter. The airframe is modelled and computed with an implicit, multi-block, multi-grid parallel Navier-Stokes solver. The flow solver has been optimised and run on up to 200 parallel processors. The cases reported include the effects of angle-of-attack (positive and negative), the effects of yaw (starboard and port) and side flow. Finally, the effects of the support strut in the wind tunnel experiments have been evaluated. Data are shown for the lift, drag and side force coefficients at flight Reynolds numbers (Re = 30m). A case of 30 degrees yaw at a flight Reynolds number is shown. We conclude that with the use of top-end computer resources it is possible to calculate the aerodynamic coefficients with a good degree of accuracy if the flow is mostly attached.

This paper describes the work done to achieve high parallel performance for an unstructured, unsteady turbomachinery computational fluid dynamics (CFD) code. The aim of the work described here is to be able to scale problems to the thousands of processors that current and future machine architectures will provide. The CFD code is in design use in industry and is also used as a research tool at a number of universities. High parallel scalability has been achieved for a range of turbomachinery test cases, from steady-state hexahedral mesh cases to fully unsteady unstructured mesh cases. This has been achieved by a combination of code modification and consideration of the parallel partitioning strategy and resulting load balancing. A sliding plane option is necessary to run fully unsteady multistage turbomachinery test cases and this has been implemented within the CFD code. Sample CFD calculations of a full turbine including parts of the internal air system are presented.

The flow-field beneath a jet-borne vertical landing aircraft is highly complex and unsteady. large-eddy simulation is a suitable tool to predict both the mean flow and unsteady fluctuations. This work aims to evaluate the suitability of LES by applying it to two multiple jet impingement problems: the first is a simple twin impinging jet in cross-flow, while the second includes a circular intake. The numerical method uses a compressible solver on a mixed element unstructured mesh. The smoothing terms in the spatial flux are kept small by the use of a monitor function sensitive to vorticity and divergence. The WALE subgrid scale model is utilised. The simpler jet impingement case shows good agreement with experiment for mean velocity and normal stresses. Analysis of time histories in the jet shear layer and near impingement gives a dominant frequency at a Strouhal number of 0·1, somewhat lower than normally observed in free jets. The jet impingement case with an intake also gives good agreement with experimental velocity measurements, although the expansion of the grid ahead of the jets does reduce the accuracy in this region. Turbulent eddies are observed entering the intake with significant swirl. This is in qualitative agreement with experimental visualisation. The results show that LES could be a suitable tool when applied to multiple jet impingement with realistic aircraft geometry.