On the effects of turbulence modelling in design optimisation for high-lift devices

Aerodynamic shape optimisation is being increasingly utilised as a design tool in the aerospace industry. In order to provide accurate results, design optimisation methods rely on the accuracy of the underlying CFD methods applied to obtain aerodynamic forces for a given configuration. Previous studies of the authors have highlighted that the variation of the order of accuracy of the CFD solver with a fixed turbulence model affects the resulting optimised airfoil shape for a single element airfoil. The accuracy of the underlying CFD model is even more relevant in the context of high-lift configurations where an accurate prediction of flow is challenging due to the complex flow physics involving transition and flow separation phenomena. This paper explores the effect of the fidelity of CFD results for a range of turbulence models within the context of the computational design of aircraft configurations. The NLR7301 multi-element airfoil (main wing and flap) is selected as the baseline configuration, because of the wealth of experimental an computational results available for this configuration. An initial validation study is conducted in order to establish optimal mesh parameters. A bi-objective shape optimisation problem is then formulated, by trying to reveal the trade-off between lift and drag coefficients at high angles of attack. Optimisation of the airfoil shape is performed with Spalart-Allmaras, k - ω SST and k - ε realisable models. The results indicate that there is consistent and complementary impact to the optimum level achieved from all the three different turbulence models considered in the presented case study. Without identifying particular superiority of any of the turbu- lence models, we can say though that each of them expressed favourable influence towards different optimality routes. These observations lead to the exploration of new avenues for future research. © 2012 by the authors.

The acoustic analogy in an annular duct with swirling mean flow

This paper is concerned with modelling the effects of swirling flow on turbomachinery noise. We develop an acoustic analogy to predict sound generation in a swirling and sheared base flow in an annular duct, including the presence of moving solid surfaces to account for blade rows. In so doing we have extended a number of classical earlier results, including Ffowcs Williams & Hawkings' equation in a medium at rest with moving surfaces, and Lilley's equation for a sheared but non-swirling jet. By rearranging the Navier-Stokes equations we find a single equation, in the form of a sixth-order differential operator acting on the fluctuating pressure field on the left-hand side and a series of volume and surface source terms on the right-hand side; the form of these source terms depends strongly on the presence of swirl and radial shear. The integral form of this equation is then derived, using the Green's function tailored to the base flow in the (rigid) duct. As is often the case in duct acoustics, it is then convenient to move into temporal, axial and azimuthal Fourier space, where the Green's function is computed numerically. This formulation can then be applied to a number of turbomachinery noise sources. For definiteness here we consider the noise produced downstream when a steady distortion flow is incident on the fan from upstream, and compare our results with those obtained using a simplistic but commonly used Doppler correction method. We show that in all but the simplest case the full inclusion of swirl within an acoustic analogy, as described in this paper, is required. © 2013 Cambridge University Press.

A simplified mathematical approach for modelling stack ventilation in multi-compartment buildings

A simple mathematical model of stack ventilation flows in multi-compartment buildings is developed with a view to providing an intuitive understanding of the physical processes governing the movement of air and heat through naturally ventilated buildings. Rules of thumb for preliminary design can be ascertained from a qualitative examination of the governing equations of flow, which elucidate the relationships between 'core' variables - flow rates, air temperatures, heat inputs and building geometry. The model is applied to an example three-storey office building with an inlet plenum and atrium. An examination of the governing equations of flow is used to predict the behaviour of steady flows and to provide a number of preliminary design suggestions. It is shown that control of ventilation flows must be shared between all ventilation openings within the building in order to minimise the disparity in flow rates between storeys, and ensure adequate fresh air supply rates for all occupants. © 2013 Elsevier Ltd.

Modelling of submarine landslides and generated water waves

The slide of unstable sedimentary bodies and their hydraulic effects are studied by numerical means. A two-dimensional fluid mechanics model based on Navier-Stokes equations has been developed considering the sediments and water as a mixture. Viscoplastic and diffusion laws for the sediments have been introduced into the model. The numerical model is validated with an analytical solution for a Bingham flow. Laboratory experiments consisting in the slide of gravel mass have been carried out. The results of these experiments have shown the importance of the sediment rheology and the diffusion. The model parameters are adjusted by trial and error to match the observed “sandflow”.

Modeling study to compare the flow and heat transfer characteristics of low-power hydrogen, nitrogen and argon arc-heated thrusters

A modelling study is performed to compare the plasma °ow and heat transfer char- acteristics of low-power arc-heated thrusters (arcjets) for three di®erent propellants: hydrogen, nitrogen and argon. The all-speed SIMPLE algorithm is employed to solve the governing equa- tions, which take into account the e®ects of compressibility, Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. The temperature, veloc- ity and Mach number distributions calculated within the thruster nozzle obtained with di®erent propellant gases are compared for the same thruster structure, dimensions, inlet-gas stagnant pressure and arc currents. The temperature distributions in the solid region of the anode-nozzle wall are also given. It is found that the °ow and energy conversion processes in the thruster nozzle show many similar features for all three propellants. For example, the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appear- ing near the cathode tip; the °ow transition from the subsonic to supersonic regime occurs within the constrictor region; the highest axial velocity appears inside the nozzle; and most of the input propellant °ows towards the thruster exit through the cooler gas region near the anode-nozzle wall. However, since the properties of hydrogen, nitrogen and argon, especially their molecular weights, speci¯c enthalpies and thermal conductivities, are di®erent, there are appreciable di®er- ences in arcjet performance. For example, compared to the other two propellants, the hydrogen arcjet thruster shows a higher plasma temperature in the arc region, and higher axial velocity but lower temperature at the thruster exit. Correspondingly, the hydrogen arcjet thruster has the highest speci¯c impulse and arc voltage for the same inlet stagnant pressure and arc current. The predictions of the modelling are compared favourably with available experimental results.

Yield drag in a two-dimensional foam flow around a circular obstacle: Effect of liquid fraction

Raufaste, C., Dollet, B., Cox, S., Jiang, Y. and Graner, F. (2007). Yield drag in a two-dimensional foam flow around a circular obstacle: Effect of liquid fraction. European Physical Journal E, 23 (2), 217?228 Sponsorship: Y.J. is supported by US DOE under contract No. DE-AC52-06NA25396. S.C. is supported by EPSRC (EP/D071127/1)

The Oscillatory Squeeze flow rheometer: Comprehensive theory and a new experimental facility

Binding, David; Bell, D.; Walters, K., (2006) 'The Oscillatory Squeeze flow rheometer: Comprehensive theory and a new experimental facility', Rheologica Acta 46 pp.111-121 RAE2008

Applications of finite element computer modelling and thermal infrared remote sensing to the study of geothermal anomalies

Buried heat sources can be investigated by examining thermal infrared images and comparing these with the results of theoretical models which predict the thermal anomaly a given heat source may generate. Key factors influencing surface temperature include the geometry and temperature of the heat source, the surface meteorological environment, and the thermal conductivity and anisotropy of the rock. In general, a geothermal heat flux of greater than 2% of solar insolation is required to produce a detectable thermal anomaly in a thermal infrared image. A heat source of, for example, 2-300K greater than the average surface temperature must be a t depth shallower than 50m for the detection of the anomaly in a thermal infrared image, for typical terrestrial conditions. Atmospheric factors are of critical importance. While the mean atmospheric temperature has little significance, the convection is a dominant factor, and can act to swamp the thermal signature entirely. Given a steady state heat source that produces a detectable thermal anomaly, it is possible to loosely constrain the physical properties of the heat source and surrounding rock, using the surface thermal anomaly as a basis. The success of this technique is highly dependent on the degree to which the physical properties of the host rock are known. Important parameters include the surface thermal properties and thermal conductivity of the rock. Modelling of transient thermal situations was carried out, to assess the effect of time dependant thermal fluxes. One-dimensional finite element models can be readily and accurately applied to the investigation of diurnal heat flow, as with thermal inertia models. Diurnal thermal models of environments on Earth, the Moon and Mars were carried out using finite elements and found to be consistent with published measurements. The heat flow from an injection of hot lava into a near surface lava tube was considered. While this approach was useful for study, and long term monitoring in inhospitable areas, it was found to have little hazard warning utility, as the time taken for the thermal energy to propagate to the surface in dry rock (several months) in very long. The resolution of the thermal infrared imaging system is an important factor. Presently available satellite based systems such as Landsat (resolution of 120m) are inadequate for detailed study of geothermal anomalies. Airborne systems, such as TIMS (variable resolution of 3-6m) are much more useful for discriminating small buried heat sources. Planned improvements in the resolution of satellite based systems will broaden the potential for application of the techniques developed in this thesis. It is important to note, however, that adequate spatial resolution is a necessary but not sufficient condition for successful application of these techniques.

On the modelling of the thermal interactions between a spray curtain and an impinging cold gas cloud

A mixed Lagrangian-Eulerian model of a Water Curtain barrier is presented. The heat, mass and momentum processes are modelled in a Lagrangian framework for the dispersed phase and in an Eulerian framework for the carrier phase. The derivation of the coupling source terms is illustrated with reference to a given carrier phase cell. The turbulent character of the flow is treated with a single equation model, modified to directly account for the influence of the particles on the flow. The model is implemented in the form of a 2 D incompressible Navier Stokes solver, coupled to an adaptive Rung Kutta method for the Lagrangian sub-system. Simulations of a free standing full cone water spray show satisfactory agreement with experiment. Predictions of a Water Curtain barrier impacted by a cold gas cloud point to markedly different flow fields for the upward and downward configurations, which could influence the effectiveness of chemical absorption in the liquid phase.

Three-dimensional free surface modelling in an unstructured mesh environment for metal processing applications

In the casting of metals, tundish flow, welding, converters, and other metal processing applications, the behaviour of the fluid surface is important. In aluminium alloys, for example, oxides formed on the surface may be drawn into the body of the melt where they act as faults in the solidified product affecting cast quality. For this reason, accurate description of wave behaviour, air entrapment, and other effects need to be modelled, in the presence of heat transfer and possibly phase change. The authors have developed a single-phase algorithm for modelling this problem. The Scalar Equation Algorithm (SEA) (see Refs. 1 and 2), enables the transport of the property discontinuity representing the free surface through a fixed grid. An extension of this method to unstructured mesh codes is presented here, together with validation. The new method employs a TVD flux limiter in conjunction with a ray-tracing algorithm, to ensure a sharp bound interface. Applications of the method are in the filling and emptying of mould cavities, with heat transfer and phase change.

The numerical modelling of multi-physics phenomena: The shape casting process

In this paper a computer simulation tool capable of modelling multi-physics processes in complex geometry has been developed and applied to the casting process. The quest for high-quality complex casting components demanded by the aerospace and automobile industries, requires more precise numerical modelling techniques and one that need to be generic and modular in its approach to modelling multi-processes problems. For such a computer model to be successful in shape casting, the complete casting process needs to be addressed, the major events being:-• Filling of hot liquid metal into a cavity mould • Solidification and latent heat evolution of liquid metal • Convection currents generated in liquid metal by thermal gradients • Deformation of cast and stress development in solidified metal • Macroscopic porosity formation The above phenomena combines the analysis of fluid flow, heat transfer, change of phase and thermal stress development. None of these events can be treated in isolation as they inexorably interact with each other in a complex way. Also conditions such as design of running system, location of feeders and chills, moulding materials and types of boundary conditions can all affect on the final cast quality and must be appropriately represented in the model.

PHYSICA - A software environment for the modelling of multi-physics phenomena

A brief description of a software environment in FORTRAN77 for the modelling of multi-physics phenomena is given. The numerical approach is based on finite volume methods but extended to unstructured meshes (ie. FV-UM). A range of interacting solution procedures for turbulent fluid flow, heat transfer with solidification/melting and elasto-visco-plastic solid mechanics are implemented in the first version of PHYSICA, which will be released in source code form to the academic community in late 1995.

Computational modelling of multi-physics processes on high performance parallel computer systems

The demands of the process of engineering design, particularly for structural integrity, have exploited computational modelling techniques and software tools for decades. Frequently, the shape of structural components or assemblies is determined to optimise the flow distribution or heat transfer characteristics, and to ensure that the structural performance in service is adequate. From the perspective of computational modelling these activities are typically separated into: • fluid flow and the associated heat transfer analysis (possibly with chemical reactions), based upon Computational Fluid Dynamics (CFD) technology • structural analysis again possibly with heat transfer, based upon finite element analysis (FEA) techniques.

Multiphysics modelling of the metals casting process

Metals casting is a process governed by the interaction of a range of physical phenomena. Most computational models of this process address only what are conventionally regarded as the primary phenomena-heat conduction and solidification. However, to predict the formation of porosity (a factor of key importance in cast quality) requires the modelling of the interaction of the fluid flow, heat transfer, solidification and the development of stress-deformation in the solidified part of a component. In this paper, a model of the casting process is described which addresses all the main continuum phenomena involved in a coupled manner. The model is solved numerically using novel finite volume unstructured mesh techniques, and then applied to both the prediction of shape deformation (plus the subsequent formation of a gap at the metal-mould interface and its impact on the heat transfer behaviour) and porosity formation in solidifying metal components. Although the porosity prediction model is phenomenologically simplistic it is based on the interaction of the continuum phenomena and yields good comparisons with available experimental results. This work represents the first of the next generation of casting simulation tools to predict aspects of the structure of cast components.

Nanoscopic modelling of the adhesion, indentation and fracture characteristics of metallic systems via molecular dynamics simulations

Large-scale molecular dynamics simulations have been performed on canonical ensembles to model the adhesion and indentation characteristics of 3-D metallic nano-scale junctions in tip-substrate geometries, and the crack propagation in 2-D metallic lattices. It is shown that irreversible flows in nano-volumes of materials control the behaviour of the 3-D nano-contacts, and that local diffusional flow constitutes the atomistic mechanism underlying these plastic flows. These simulations show that the force of adhesion in metallic nano-contacts is reduced when adsorbate monolayers are present at the metal—metal junctions. Our results are in agreement with the conclusions of very accurate point-contact experiments carried out in this field. Our fracture simulations reveal that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures, the nucleation of dislocations is shown to cause a brittle-to-ductile transition. Limiting crack propagation velocities are computed for different strain rates and a dynamic instability is shown to control the crack movement beyond this limiting velocity, in line with the recent experimental results.