The development of turbine exit flow in a swan-necked inter-stage diffuser

This paper describes both the migration and dissipation of flow phenomena downstream of a transonic high-pressure turbine stage. The geometry of the HP stage exit duct considered is a swan-necked diffuser similar to those likely to be used in future engine designs. The paper contains results both from an experimental programme in a turbine test facility and from numerical predictions. Experimental data was acquired using three fast-response aerodynamic probes capable of measuring Mach number, whirl angle, pitch angle, total pressure and static pressure. The probes were used to make time-resolved area traverses at two axial locations downstream of the rotor trailing edge. A 3D time-unsteady viscous Navier-Stokes solver was used for the numerical predictions. The unsteady exit flow from a turbine stage is formed from rotordependent phenomena (such as the rotor wake, the rotor trailing edge recompression shock, the tip-leakage flow and the hub secondary flow) and vane-rotor interaction dependant phenomena. This paper describes the time-resolved behaviour and three-dimensional migration paths of both of these phenomena as they convect downstream. It is shown that the inlet flow to a downstream vane is dominated by two corotating vortices, the first caused by the rotor tip-leakage flow and the second by the rotor hub secondary flow. At the inlet plane of the downstream vane the wake is extremely weak and the radial pressure gradient is shown to have caused the majority of the high loss wake fluid to be located between the mid-height of the passage and the casing wall. The structure of the flow indicates that between a high pressure stage and a downstream vane simple two-dimensional blade row interaction does not occur. The results presented in this paper indicate that the presence of an upstream stage is likely to significantly alter the structure of the secondary flow within a downstream vane. The paper also shows that vane-rotor interaction within the upstream stage causes a 10° circumferential variation in the inlet flow angle of the 2nd stage vane.

Properties of tetrahedral amorphous carbon prepared by vacuum arc deposition

The structural, optical, electrical and physical properties of amorphous carbon deposited from the filtered plasma stream of a vacuum arc were investigated. The structure was determined by electron diffraction, neutron diffraction and energy loss spectroscopy and the tetrahedral coordination of the material was confirmed. The measurements gave a nearest neighbour distance of 1.53 Å, a bond angle of 110 and a coordination number of four. A model is proposed in which the compressive stress generated in the film by energetic ion impact produces pressure and temperature conditions lying well inside the region of the carbon phase diagram within which diamond is stable. The model is confirmed by measurements of stress and plasmon energy as a function of ion energy. The model also predicts the formation of sp2-rich materials on the surface owing to stress relaxation and this is confirmed by a study of the surface plasmon energy. Some nuclear magnetic resonance, infrared and optical properties are reported and the behaviour of diodes using tetrahedral amorphous carbon is discussed. © 1991.

Properties of grain boundaries in bulk, melt processed YB2Cu3O7-δ fabricated using bridge-shaped seeds

Single grain REBa2C3uO7 ((RE)BCO, where RE is a rare earth element or yttrium) bulk superconducting materials have significant potential for a variety of engineering applications due to their ability to trap high magnetic fields. However, it is well known that the presence of grain boundaries coupled with a high angle of misorientation (typically 5�) significantly reduces the critical current density, J c , in all forms of high temperature superconducting materials. It is of considerable fundamental and technological interest, therefore, to investigate the grain boundary properties of bulk, film and tape (RE)BCO. We report a successful multi-seeding technique for the fabrication of fully aligned, artificial (0��misalignment) grain boundaries within large grain YBCO bulk superconductors using bridge-shaped seeds. The microstructure and critical current densities of the grain boundaries produced by this technique have been studied in detail.

Tsunami wave runup on the clustered islands using boussinesq-type model

This paper describes the implementation of the Boussinesq-type model and extends its application to the tsunami wave runup on the clustered islands (multiple adjacent conical islands), in turn, an extensively validated two-dimensional Boussinesq-type model is employed to examine the interaction between a propagating solitary wave and multiple idealised conical islands, with particular emphasis on a combination effect of two adjustable parameters for spacing interval/diameter ratio between the adjacent conical islands, S/D, and the rotating angle of the structural configuration,θ on maximum soliton runup heights. An extensive parameter study concerning the combination effect of alteringθ and S/D on the maximum soliton runup with the multi-conical islands is subsequently carried out and the distributions of the maximum runup heights on each conical island are obtained and compared for the twin-island cases. The worst case study is performed for each case in respect of the enhancement in the maximum wave runup heights by the multi-conical islands. It is found that the nonlinear wave diffraction, reflection and refraction play a significant role in varying the maximum soliton runup heights on multiconical islands. The comparatively large maximum soliton runups are generally predicted for the merged and bottom mounted clusteredislands. Furthermore, the joints of the clustered-merged islands are demonstrated to suffer the most of the tsunami wave attack. The conical islands that position in the shadow regions behind the surrounding islands are found to withstand relatively less extreme wave impact. Although, these numerical investigations are considerable simplifications of the multi conical islands, they give a critical insight into certain important hydrodynamic characteristics of the interaction between an extreme wave event and a group of clustered conical islands, and thus providing a useful engineering guidance for extreme wave mitigation and coastal development. Copyright © 2012 by the International Society of Offshore and Polar Engineers (ISOPE).

Occurrence of particle debris field during focused Ga ion beam milling of glassy carbon

To explore the machining characteristics of glassy carbon by focused ion beam (FIB), particles induced by FIB milling on glassy carbon have been studied in the current work. Nano-sized particles in the range of tens of nanometers up to 400 nm can often be found around the area subject to FIB milling. Two ion beam scanning modes - slow single scan and fast repetitive scan - have been tested. Fewer particles are found in single patterns milled in fast repetitive scan mode. For a group of test patterns milled in a sequence, it was found that a greater number of particles were deposited around sites machined early in the sequence. In situ EDX analysis of the particles showed that they were composed of C and Ga. The formation of particles is related to the debris generated at the surrounding areas, the low melting point of gallium used as FIB ion source and the high contact angle of gallium on glassy carbon induces de-wetting of Ga and the subsequent formation of Ga particles. Ultrasonic cleaning can remove over 98% of visible particles. The surface roughness (Ra) of FIB milled areas after cleaning is less than 2 nm. © 2010.

An integrated model of rolling tyre vibrations on a smooth road

We describe new results on the vibrations of rolling tyres, aimed at noise prediction for tyres of given design on a smooth road surface. This new approach incorporates our existing models, of smooth road-tyre interaction and belt vibration but includes additional features that are required for real tyre patterns. To this end, the model allows variable tread block size and grooves along the belt circumference; the density and angle of these grooves may also vary laterally. The key innovation is to treat the tyre belt as a laterally stacked series of rings, each of which is equipped with a set of viscoelastic springs around its circumference. It is shown how to use this construction to mimic the details of actual tyre patterns and, in conjunction with existing models, predict belt vibrations. The construction is applied to develop a ring discretisation for a real tyre that shows strong lateral variations. It is shown that the vibration amplitude is concentrated on a set of parallel lines in frequency-wavenumber space and that the tread pattern dictates the occurrence and spacing of these lines. Linkage to a boundary element calculation then allows quantification of the influence of tread parameters on radiated noise. Keywords: Vibration, tread pattern, tyre noise. Copyright © (2011) by the Institute of Noise Control Engineering.

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.

A body force-based method for prediction of multiple-pure-tone noise: Validation

Emissions, fuel burn, and noise are the main drivers for innovative aircraft design. Embedded propulsion systems, such as for example used in hybrid-wing body aircraft, can offer fuel burn and noise reduction benefits but the impact of inlet flow distortion on the generation and propagation of turbomachinery noise has yet to be assessed. A novel approach is used to quantify the effects of non-uniform flow on the creation and propagation of multiple pure tone (MPT) noise. The ultimate goal is to conduct a parametric study of S-duct inlets to quantify the effects of inlet design parameters on the acoustic signature. The key challenge is that the effects of distortion transfer, noise source generation and propagation through the non-uniform flow field are inherently coupled such that a simultaneous computation of the aerodynamics and acoustics is required to capture the mechanisms at play. The technical approach is based on a body force description of the fan blade row that is able to capture the distortion transfer and the blade-to-blade flow variations that cause the MPT noise while reducing computational cost. A single, 3-D full-wheel CFD simulation, in which the Euler equations are solved to second-order spatial and temporal accuracy, simultaneously computes the MPT noise generation and its propagation in distorted inlet flow. A new method of producing the blade-to-blade variations in the body force field for MPT noise generation has been developed and validated. The numerical dissipation inherent to the solver is quantified and used to correct for non-physical attenuation in the far-field noise spectra. Source generation, acoustic propagation and acoustic energy transfer between modes is examined in detail. The new method is validated on NASA's Source Diagnostic Test fan and inlet, showing good agreement with experimental data for aerodynamic performance, acoustic source generation, and far-field noise spectra. The next steps involve the assessment of MPT noise in serpentine inlet ducts and the development of a reduced order formulation suitable for incorporation into NASA's ANOPP framework. © 2010 by Jeff Defoe, Alex Narkaj & Zoltan Spakovszky.

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.

Identification of the steering control behaviour of five test subjects following a randomly curving path in a driving simulator

The paper is concerned with the identification of theoretical preview steering controllers using data obtained from five test subjects in a fixed-base driving simulator. An understanding of human steering control behaviour is relevant to the design of autonomous and semi-autonomous vehicle controls. The driving task involved steering a linear vehicle along a randomly curving path. The theoretical steering controllers identified from the data were based on optimal linear preview control. A direct-identification method was used, and the steering controllers were identified so that the predicted steering angle matched as closely as possible the measured steering angle of the test subjects. It was found that identification of the driver's time delay and noise is necessary to avoid bias in identification of the controller parameters. Most subjects' steering behaviour was predicted well by a theoretical controller based on the lateral/yaw dynamics of the vehicle. There was some evidence that an inexperienced driver's steering action was better represented by a controller based on a simpler model of the vehicle dynamics, perhaps reflecting incomplete learning by the driver. Copyright © 2014 Inderscience Enterprises Ltd.

Initiation of immersed granular avalanches.

By means of coupled molecular dynamics-computational fluid dynamics simulations, we analyze the initiation of avalanches in a granular bed of spherical particles immersed in a viscous fluid and inclined above its angle of repose. In quantitative agreement with experiments, we find that the bed is unstable for a packing fraction below 0.59 but is stabilized above this packing fraction by negative excess pore pressure induced by the effect of dilatancy. From detailed numerical data, we explore the time evolution of shear strain, packing fraction, excess pore pressures, and granular microstructure in this creeplike pressure redistribution regime, and we show that they scale excellently with a characteristic time extracted from a model based on the balance of granular stresses in the presence of a negative excess pressure and its interplay with dilatancy. The cumulative shear strain at failure is found to be ≃ 0.2, in close agreement with the experiments, irrespective of the initial packing fraction and inclination angle. Remarkably, the avalanche is triggered when dilatancy vanishes instantly as a result of fluctuations while the average dilatancy is still positive (expanding bed) with a packing fraction that declines with the initial packing fraction. Another nontrivial feature of this creeplike regime is that, in contrast to dry granular materials, the internal friction angle of the bed at failure is independent of dilatancy but depends on the inclination angle, leading therefore to a nonlinear dependence of the excess pore pressure on the inclination angle. We show that this behavior may be described in terms of the contact network anisotropy, which increases with a nearly constant connectivity and levels off at a value (critical state) that increases with the inclination angle. These features suggest that the behavior of immersed granular materials is controlled not only directly by hydrodynamic forces acting on the particles but also by the influence of the fluid on the granular microstructure.

Laboratory-scale fluidized bed rig for high-temperature tube wastage studies

As part of a study of the wear of candidate heat exchanger tube materials for use in fluidized bed combustors, two similar laboratory-scale rigs have been built and characterized. Specimens of selected alloys are carried on counter-rotating rotors immersed in a fluidized bed, and are exposed to particle impact velocities of up to approximately 3 ms-1 at temperatures up to 1000°C. The performance of this design of apparatus has been investigated in detail. The effects of several experimental variables have been studied, including angle of particle impact, specimen speed, position of the rotor within the fluidized bed, duration of exposure, bed material particle size, degradation of the bed material, degree of fluidization of the bed, and size of specimen. In many cases the results obtained with steel specimens at elevated temperatures are similar to those observed with polymeric specimens at low temperatures.

Computer generated hologram with geometric occlusion using GPU-accelerated depth buffer rasterization for three-dimensional display

We present a method of rapidly producing computer-generated holograms that exhibit geometric occlusion in the reconstructed image. Conceptually, a bundle of rays is shot from every hologram sample into the object volume.We use z buffering to find the nearest intersecting object point for every ray and add its complex field contribution to the corresponding hologram sample. Each hologram sample belongs to an independent operation, allowing us to exploit the parallel computing capability of modern programmable graphics processing units (GPUs). Unlike algorithms that use points or planar segments as the basis for constructing the hologram, our algorithm's complexity is dependent on fixed system parameters, such as the number of ray-casting operations, and can therefore handle complicated models more efficiently. The finite number of hologram pixels is, in effect, a windowing function, and from analyzing the Wigner distribution function of windowed free-space transfer function we find an upper limit on the cone angle of the ray bundle. Experimentally, we found that an angular sampling distance of 0:01' for a 2:66' cone angle produces acceptable reconstruction quality. © 2009 Optical Society of America.

Micro-vortex generator flow control for supersonic engine inlets

To investigate the flow control potential of micro-vortex generators for supersonic mixed-compression inlets, a basic model experiment has been designed which combines a normal shock wave with a subsonic diffuser. The diffuser is formed by a simple expansion corner, with a divergence angle of 6 degrees. The diffuser entry Mach numbers were M=1.3 and M=1.5 and a number of shock locations relative to the corner position were tested. Flow control was applied in the form of counter-rotating micro-vanes with heights of approximately 20% of boundary layer thickness. Furthermore, corner fences where employed to reduce sidewall effects. It was found that micro-vortex generators were able to significantly reduce the extent of flow separation under all conditions, but could not eliminate it altogether. Corner fences also demonstrated potential for improving the flow in rectangular cross section channels and the combination of corner fences with micro-vortex generators was found to give the greatest benefits. At M=1.3 the combination of corner fences and micro-vanes placed close to the diffuser entry could prevent separation for a wide range of conditions. At the higher diffuser entry Mach number the benefits of flow control were less significant although a reduction of separation size and an improved pressure recovery was observed. It is thought that micro-vortex generators can have significant flow control potential if they are placed close to the expected separation onset and when the adverse pressure gradient is not too far above the incipient separation level. The significant beneficial effects of corner fences warrant a more comprehensive further investigation. It is thought that the control methods suggested here are capable of reducing the bleed requirement on an inlet, which could provide significant performance advantages.