Centrifuge modelling of submarine landslide flows

Landslides occur both onshore and offshore, however little attention has been given to offshore landslides (submarine landslides). The unique characteristics of submarine landslides include large mass movements and long travel distances at very gentle slopes. Submarine landslides have significant impacts and consequences on offshore and coastal facilities. This paper presents data from a series of centrifuge tests simulating submarine landslide flows on a very gentle slope. Experiments were conducted at different gravity levels to understand the scaling laws involved in simulating submarine landslide flows through centrifuge modelling. The slope was instrumented with miniature sensors for measurements of pore pressure beneath the flow. A series of digital cameras were used to capture the flow in flight. The results provide a better understanding of the scaling laws that needs to be adopted for centrifuge experiments involving submarine landslide flows and gives an insight into the flow mechanisms. © 2010 Taylor & Francis Group, London.

Low-order modelling of ducted flames with temporally varying equivalence ratio in realistic geometries

The interaction between unsteady heat release and acoustic pressure oscillations in gas turbines results in self-excited combustion oscillations which can potentially be strong enough to cause significant structural damage to the combustor. Correctly predicting the interaction of these processes, and anticipating the onset of these oscillations can be difficult. In recent years much research effort has focused on the response of premixed flames to velocity and equivalence ratio perturbations. In this paper, we develop a flame model based on the socalled G-Equation, which captures the kinematic evolution of the flame surfaces, under the assumptions of axisymmetry, and ignoring vorticity and compressibility. This builds on previous work by Dowling [1], Schuller et al. [2], Cho & Lieuwen [3], among many others, and extends the model to a realistic geometry, with two intersecting flame surfaces within a non-uniform velocity field. The inputs to the model are the free-stream velocity perturbations, and the associated equivalence ratio perturbations. The model also proposes a time-delay calculation wherein the time delay for the fuel convection varies both spatially and temporally. The flame response from this model was compared with experiments conducted by Balachandran [4, 5], and found to show promising agreement with experimental forced case. To address the primary industrial interest of predicting self-excited limit cycles, the model has then been linked with an acoustic network model to simulate the closed-loop interaction between the combustion and acoustic processes. This has been done both linearly and nonlinearly. The nonlinear analysis is achieved by applying a describing function analysis in the frequency domain to predict the limit cycle, and also through a time domain simulation. In the latter case, the acoustic field is assumed to remain linear, with the nonlinearity in the response of the combustion to flow and equivalence ratio perturbations. A transfer function from unsteady heat release to unsteady pressure is obtained from a linear acoustic network model, and the corresponding Green function is used to provide the input to the flame model as it evolves in the time domain. The predicted unstable frequency and limit cycle are in good agreement with experiment, demonstrating the potential of this approach to predict instabilities, and as a test bench for developing control strategies. Copyright © 2011 by ASME.

Transient flows in enclosures: A generalised approach for modelling the effects of geometry, heat gains and wind

Transient flows in a confined ventilated space induced by a buoyancy source of time-varying strength and an external wind are examined. The space considered has varying cross-sectional area with height. A generalised theoretical model is proposed to investigate the flow dynamics following the activation of an external wind and an internal source of buoyancy. To investigate the effect of geometry, we vary the angle of the wall inclination of a particular geometry in which a point source of constant buoyancy is activated in the absence of wind. Counter-intuitively the ventilation is worse and lower airflow rates are established for geometries of increasing cross-sectional areas with height. We investigate the effect of the source buoyancy strength by comparing two cases: (1) when the buoyancy input is constant and (2) when the buoyancy input gradually increases over time so that after a finite time the total buoyancy inputs for (1) and (2) are identical. The rate at which the source heat gains are introduced has a significant role on the flow behaviour as we find that, in case (2), a warmer layer and a more pronounced overshoot are obtained than in case (1). The effect of assisting and opposing wind on the transient ventilation of an enclosure of constant cross-sectional area with height and constant heat gains is examined. A Froude number Fr is used to define the relative strengths of the buoyancy-induced and wind-induced velocities and five different transient states and their associated critical Fr are identified. © 2010 Elsevier Ltd.

Heat source modelling and natural ventilation efficiency

We compare natural ventilation flows established by a range of heat source distributions at floor level. Both evenly distributed and highly localised line and point source distributions are considered. We demonstrate that modelling the ventilation flow driven by a uniformly distributed heat source is equivalent to the flow driven by a large number of localised sources. A model is developed for the transient flow development in a room with a uniform heat distribution and is compared with existing models for localised buoyancy inputs. For large vent areas the flow driven by localised heat sources reaches a steady state more rapidly than the uniformly distributed case. For small vent areas there is little difference in the transient development times. Our transient model is then extended to consider the time taken to flush a neutrally buoyant pollutant from a naturally ventilated room. Again comparisons are drawn between uniform and localised (point and line) heat source geometries. It is demonstrated that for large vent areas a uniform heat distribution provides the fastest flushing. However, for smaller vent areas, localised heat sources produce the fastest flushing. These results are used to suggest a definition for the term 'natural ventilation efficiency', and a model is developed to estimate this efficiency as a function of the room and heat source geometries. © 2006 Elsevier Ltd. All rights reserved.

Long range exhaustion-a mathematical model for the axisymmetric air flow of a local exhaust ventilation hood assisted by a turbulent radial jet

Abstract-Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art Banged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air Bow to be reduced and provide a means of identifying and assessing the various parameters that control the air Bow. The mathematical model is formulated in terms of the Stokes steam function, ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained Bow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions. | Mathematical modelling techniques are used to predict the axisymmetric air flow pattern developed by a state-of-the-art flanged exhaust hood which is reinforced by a turbulent radial jet flow. The high Reynolds number modelling techniques adopted allow the complexity of determining the hood's air flow to be reduced and provide a means of identifying and assessing the various parameters that control the air flow. The mathematical model is formulated in terms of the Stokes steam function, Ψ, and the governing equations of fluid motion are solved using finite-difference techniques. The injection flow of the exhaust hood is modelled as a turbulent radial jet and the entrained flow is assumed to be an inviscid potential flow. Comparisons made between contours of constant air speed and centre-line air speeds deduced from the model and all the available experimental data show good agreement over a wide range of typical operating conditions.

Effects of exhaust inlet size on the fluid flow patterns created by an Aaberg exhaust hood

The aim of this paper is to investigate the effects of the size of the exhaust inlet on the fluid flow patterns created by an Aaberg exhaust hood. Streamlines and lines of constant speed modelling those created by an Aaberg hood are used to examine the effect of the inlet size on the size and profile of the hood's effective suction area. © 1992.

Magnetic shielding properties of a superconducting hollow cylinder containing slits: Modelling and experiment

This paper deals with the magnetic properties of bulk high temperature superconducting cylinders used as magnetic shields. We investigate, both numerically and experimentally, the magnetic properties of a hollow cylinder with two axial slits which cut the cylinder in equal halves. Finite element method modelling has been used with a three-dimensional geometry to help us in understanding how the superconducting currents flow in such a cut cylinder and therefore how the magnetic shielding properties are affected, depending on the magnetic field orientation. Modelling results show that the slits block the shielding current flow and act as an entrance channel for the magnetic flux lines. The contribution of the slits to the total flux density that enters the cylinder is studied through the angle formed between the applied field and the internal field. The modelled data agree nicely with magnetic shielding properties measured on a bulk Bi-2212 hollow cylinder at 77K. The results demonstrate that the magnetic flux penetration in such a geometry can be modelled successfully using only two parameters of the superconductor (constant J c and n value), which were determined from magnetic measurements on the plain cylinder. © 2012 IOP Publishing Ltd.

Low-order modelling of the response of ducted flames in annular geometries

The adoption of lean premixed prevaporised combustion systems can reduce NOx emissions from gas turbines, but unfortunately also increases their susceptibility to thermoacoustic instabilities. Initially, acoustic waves can produce heat release fluctuations by a variety of mechanisms, often by perturbing the equivalence ratio. If correctly phased, heat release fluctuations can subsequently generate more acoustic waves, which at high amplitude can result in significant structural damage to the combustor. The prediction of this phenomenon is of great industrial interest. In previous work, we have coupled a physics based, kinematic model of the flame with a network model to provide the planar acoustic response necessary to close the feedback loop and predict the onset and amplitude of thermoacoustic instabilities in a lab-scale, axisymmetric single burner combustor. The advantage of a time domain approach is that the modal interaction, the influence of harmonics, and flame saturation can be investigated. This paper extends this approach to more realistic, annular geometries, where both planar and circumferential modes must be considered. In lean premixed prevaporised combustors, fluctuations in equivalence ratio have been shown to be a dominant cause of unsteady combustion. These can occur, for example, due to velocity perturbations in the premix ducts, which can lead to equivalence ratio fluctuations at the fuel injectors, which are subsequently convected downstream to the flame surfaces. Here, they can perturb the heat release by locally altering the flame speed, enthalpy of combustion, and, indirectly, the flame surface area. In many gas turbine designs, particularly aeroengines, the geometries are composed of a ring of premix ducts linking a plenum and an annular combustor. The most unstable modes are often circumferential modes. The network model is used to characterise the flow response of the geometry to heat fluctuations at an appropriate location, such as the fuel injectors. The heat release at each flame holder is determined in the time domain using the kinematic flame model derived, as a function of the flow perturbations in the premix duct. This approach is demonstrated for an annular ring of burners on a in a simple geometry. The approach is then extended to an industrial type gas turbine combustor, and used to predict the limit cycle amplitudes. Copyright © 2012 by ASME.

Multi-scale modelling of granular avalanches

Avalanches, debris flows, and landslides are geophysical hazards, which involve rapid mass movement of granular solids, water and air as a single-phase system. The dynamics of a granular flow involve at least three distinct scales: the micro-scale, meso-scale, and the macro-scale. This study aims to understand the ability of continuum models to capture the micro-mechanics of dry granular collapse. Material Point Method (MPM), a hybrid Lagrangian and Eulerian approach, with Mohr-Coulomb failure criterion is used to describe the continuum behaviour of granular column collapse, while the micromechanics is captured using Discrete Element Method (DEM) with tangential contact force model. The run-out profile predicted by the continuum simulations matches with DEM simulations for columns with small aspect ratios ('h/r' < 2), however MPM predicts larger run-out distances for columns with higher aspect ratios ('h/r' > 2). Energy evolution studies in DEM simulations reveal higher collisional dissipation in the initial free-fall regime for tall columns. The lack of a collisional energy dissipation mechanism in MPM simulations results in larger run-out distances. Micro-structural effects, such as shear band formations, were observed both in DEM and MPM simulations. A sliding flow regime is observed above the distinct passive zone at the core of the column. Velocity profiles obtained from both the scales are compared to understand the reason for a slow flow run-out mobilization in MPM simulations. © 2013 AIP Publishing LLC.

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 - o 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 AIAA.

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.