Displacement and mixing ventilation driven by opposing wind and buoyancy

The effect of an opposing wind on the stratification and flow produced by a buoyant plume rising from a heat source on the floor of a ventilated enclosure is investigated. Ventilation openings located at high level on the windward side of the enclosure and at low level on the leeward side allow a wind-driven flow from high to low level, opposite to the buoyancy-driven flow. One of two stable steady flow regimes is established depending on a dimensionless parameter F that characterizes the relative magnitudes of the wind-driven and buoyancy-driven velocities within the enclosure, and on the time history of the flow. A third, unstable steady flow solution is identified. For small opposing winds (small F) a steady, two-layer stratification and displacement ventilation is established. Exterior fluid enters through the lower leeward openings and buoyant interior fluid leaves through the upper windward openings. As the wind speed increases, the opposing wind may cause a reversal in the flow direction. In this case, cool exterior fluid enters through the high windward openings and mixes the interior fluid, which exits through the leeward openings. There are now two possibilities. If the rate of heat input by the source exceeds the rate of heat loss through the leeward openings, the temperature of the interior increases and this flow reversal is only maintained temporarily. The buoyancy force increases with time, the flow reverts to its original direction, and steady two-layer displacement ventilation is re-established and maintained. In this regime, the increase in wind speed increases the depth and temperature of the warm upper layer, and reduces the ventilation flow rate. If, on the other hand, the heat loss exceeds the heat input, the interior cools and the buoyancy-driven flow decreases. The reversed flow is maintained, the stratification is destroyed and mixing ventilation occurs. Further increases in wind speed increase the ventilation rate and decrease the interior temperature. The transitions between the two ventilation flow patterns exhibit hysteresis. The change from displacement ventilation to mixing ventilation occurs at a higher F than the transition from mixing to displacement. Further, we find that the transition from mixing to displacement ventilation occurs at a fixed value of F, whereas the transition from displacement to mixing flow is dependent on the details of the time history of the flow and the geometry of the openings, and is not determined solely by the value of F. Theoretical models that predic t the steady stratification profiles and flow rates for the displacement and mixing ventilation, and the transitions between them, are presented and compared with measurements from laboratory experiments. The transition between these ventilation patterns completely changes the internal environment, and we discuss some of the implications for the natural ventilation of buildings. © 2004 Cambridge University Press.

Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

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.

Comparison of les to LDA and PIV measurements of a small scale high speed coflowing jet

Measurements and predictions are made of a short cowl co-flowing jet with a bypass ratio of 8:1. The Reynolds number for computations and measurements are matched at 300,000 and the Mach numbers representative of realistic jet conditions with core and co flow velocities of 240m/s and 216m/s respectively. The low Reynolds number of the measurements makes the case well suited to the assessment of large eddy resolving computational strategies. Also, the nozzle concentricity was carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both two dimensional coincident LDA and PIV systems. The computational simulations are completed on a modest 12×106 mesh. The simulation is now being run on a 50×106 mesh using hybrid RANSNLES (Numerical Large Eddy Simulation). Close to the nozzle wall a k-l RANS model is used. For an axisymmetric jet, comparison is made between simulations which use NLES, RANSNLES and also a simple imposed velocity profile where the nozzle is not modeled. The use of a near wall RANS model is shown to be beneficial. When compared with the measurements the NLES results are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Numerical computation of laminar flow in a heated rotating cavity with an axial throughflow of air

A heated rotating cavity with an axial throughflow of cooling air is used as a model for the flow in the cylindrical cavities between adjacent discs of a high-pressure gas-turbine compressor. In an engine the flow is expected to be turbulent, the limitations of this laminar study are fully realised but it is considered an essential step to understand the fundamental nature of the flow. The three-dimensional, time-dependent governing equations are solved using a code based on the finite volume technique and a multigrid algorithm. The computed flow structure shows that flow enters the cavity in one or more radial arms and then forms regions of cyclonic and anticyclonic circulation. This basic flow structure is consistent with existing experimental evidence obtained from flow visualization. The flow structure also undergoes cyclic changes with time. For example, a single radial arm, and pair of recirculation regions can commute to two radial arms and two pairs of recirculation regions and then revert back to one. The flow structure inside the cavity is found to be heavily influenced by the radial distribution of surface temperature imposed on the discs. As the radial location of the maximum disc temperature moves radially outward, this appears to increase the number of radial arms and pairs of recirculation regions (from one to three for the distributions considered here). If the peripheral shroud is also heated there appear to be many radial arms which exchange fluid with a strong cyclonic flow adjacent to the shroud. One surface temperature distribution is studied in detail and profiles of the relative tangential and radial velocities are presented. The disc heat transfer is also found to be influenced by the disc surface temperature distribution. It is also found that the computed Nusselt numbers are in reasonable accord over most of the disc surface with a correlation found from previous experimental measurements. © 1994, MCB UP Limited.

Quantification of the benefits of a streamlined tether for a high-altitude-tethered balloon

A high-altitude tethered balloon (HATB) reaching a height of 20 km has numerous applications including communications, meteorological monitoring, surveillance and, for the current paper, geoengineering. The HATB configuration consists of a 20 km tether rising up from the ground through the troposphere, where the wind-regime can be turbulent and include high-wind velocities due to the jet-stream, up into the more stable stratosphere where the tether would be attached to a spherical balloon. This paper evaluates wind-excited vibration of a HATB and will investigate the advantages of using a streamlined instead of a circular tether profile. Streamlining the tether reduces drag but introduces stability problems and complicates the numerical modelling. Consequently, prior to a thorough investigation of the stability issues due to fluid-structure interaction, a 3D nonlinear lumped-mass HATB model is used to quantify the benefit of a streamlined tether. The benefit is quantified by comparing the system specifications - such as balloon size and the tension in the tether - required to meet certain design requirements driven by the Stratospheric Particle Injection for Climate Engineering (SPICE) project. The SPICE project is investigating the feasibility of climate engineering using a HATB.

Flow field measurements of pulverized coal combustion using optical diagnostic techniques

This paper demonstrates the application of laser Doppler velocimetry (LDV) and particle image velocimetry (PIV) techniques to a particle-laden reacting flow of pulverized coal. A laboratory-scale open-type annular burner is utilized to generate velocity profiles of coal particles and micrometric alumina particles. Pair-wise two-component LDV measurements and high-speed stereo PIV measurements provide three-dimensional velocity components of the flow field. A detailed comparison of velocities for alumina and coal particle seeding revealed differences attributed to the wide size distribution of coal particles. In addition, the non-spherical shape and high flame luminosity associated with coal particle combustion introduces noise to the Mie scatter images. The comparison of mean and RMS velocities measured by LDV and PIV techniques showed that PIV measurements are affected by the wide size distribution of coal particles, whereas LDV measurements become biased toward the velocity of small particles, as signals from large particles are rejected. This small-particle bias is also reflected in the spectral characteristics for both techniques, which are in good agreement within the range of frequencies accessible. PIV measurements showed an expected lack of response of large coal particles to the turbulence fluctuations. The overall good agreement between LDV and PIV measurements demonstrates the applicability of the high-speed PIV technique to a particle-laden, high luminosity coal flame while highlighting some of its limitations. © 2013 Springer-Verlag Berlin Heidelberg.

Wave propagation velocity under a vertically vibrated surface foundation

The ultimate objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of soil modulus as a function of strain. In field experiments, an excitation is applied on the ground surface using large-scale shakers, and the response of the soil deposit is recorded through receivers embedded in the soil. The focus of this paper is on the simulation and observation of signals that would be recorded at the receiver locations under idealized conditions to provide guidelines on the interpretation of the field measurements. Discrete models are used to reproduce one-dimensional and three-dimensional geometries. When the first times of arrival are detected by receivers under the vertical impulse, they coincide with the arrival of the P wave; therefore related to the constrained modulus of the material. If one considers, on the other hand, phase differences between the motions at two receivers, the picture is far more complicated and one would obtain propagation velocities, function of frequency and measuring location, which do not correspond to either the constrained modulus or Young's modulus. It is necessary then to conduct more rigorous and complicated analyses in order to interpret the data. This paper discusses and illustrates these points. Copyright © 2008 John Wiley & Sons, Ltd.

Wave propagation under vertically excited surface foundation

Recently developed equipment allows measurement of the shear modulus of soil in situ as a function of level of strain. In these field experiments, the excitation is applied on the ground surface using large scale shakers, and the response of the soil deposit is recorded through embedded receivers. The focus of this paper is on the simulation of signals which would be recorded at the receiver locations in idealized conditions to provide guidelines on the interpretation of field measurements. Discrete and finite element methods are employed to model one dimensional and three dimensional geometries, respectively, under various lateral boundary conditions. When the first times of arrival are detected by receivers under the vertical impulse, they coincide with the arrival of the P wave, related to the constrained modulus of the material, regardless of lateral boundary conditions. If one considers, on the other hand, phase differences between the motions at two receivers the picture is far more complicated and one would obtain propagation velocities, function of frequency and depth, which do not correspond to either the constrained modulus or Young's modulus. It is thus necessary to apply some care when interpreting the data from field tests based on vertical steady state vibrations. The use of inverse analysis can be considered as a way of extracting the shear modulus of soil from the field test measurements. © 2008 ASCE.

The effect of shear strength on the ballistic response of laminated composite plates

The ballistic performance of clamped circular carbon fibre reinforced polymer (CFRP) and Ultra High Molecular Weight Polyethylene (UHMWPE) fibre composite plates of equal areal mass and 0/90 lay-up were measured and compared with that of monolithic 304 stainless steel plates. The effect of matrix shear strength upon the dynamic response was explored by testing: (i) CFRP plates with both a cured and uncured matrix and (ii) UHMWPE laminates with identical fibres but with two matrices of different shear strength. The response of these plates when subjected to mid-span, normal impact by a steel ball was measured via a dynamic high speed shadow moiré technique. Travelling hinges emanate from the impact location and travel towards the supports. The anisotropic nature of the composite plate results in the hinges travelling fastest along the fibre directions and this results in square-shaped moiré fringes in the 0/90 plates. Projectile penetration of the UHMWPE and the uncured CFRP plates occurs in a progressive manner, such that the number of failed plies increases with increasing velocity. The cured CFRP plate, of high matrix shear strength, fails by cone-crack formation at low velocities, and at higher velocities by a combination of cone-crack formation and communition of plies beneath the projectile. On an equal areal mass basis, the low shear strength UHMWPE plate has the highest ballistic limit followed by the high matrix shear strength UHMWPE plate, the uncured CFRP, the steel plate and finally the cured CFRP plate. We demonstrate that the high shear strength UHMWPE plate exhibits Cunniff-type ballistic limit scaling. However, the observed Cunniff velocity is significantly lower than that estimated from the laminate properties. The data presented here reveals that the Cunniff velocity is limited in its ability to characterise the ballistic performance of fibre composite plates as this velocity is independent of the shear properties of the composites: the ballistic limit of fibre composite plates increases with decreasing matrix shear strength for both CFRP and UHMWPE plates. © 2013 Elsevier Masson SAS. All rights reserved.

Autonomous rigid body attitude synchronization

This paper studies some extensions to the decentralized attitude synchronization of identical rigid bodies. Considering fully actuated Euler equations, the communication links between the rigid bodies are limited and the available information is restricted to relative orientations and angular velocities. In particular, no leader nor external reference dictates the swarm's behavior. The control laws are derived using two classical approaches of nonlinear control - tracking and energy shaping. This leads to a comparison of two corresponding methods which are currently considered for distributed synchronization - consensus and stabilization of mechanical systems with symmetries. © 2007 IEEE.

Electronic properties of GaAs, InAs and InP nanowires studied by terahertz spectroscopy.

We have performed a comparative study of ultrafast charge carrier dynamics in a range of III-V nanowires using optical pump-terahertz probe spectroscopy. This versatile technique allows measurement of important parameters for device applications, including carrier lifetimes, surface recombination velocities, carrier mobilities and donor doping levels. GaAs, InAs and InP nanowires of varying diameters were measured. For all samples, the electronic response was dominated by a pronounced surface plasmon mode. Of the three nanowire materials, InAs nanowires exhibited the highest electron mobilities of 6000 cm² V⁻¹ s⁻¹, which highlights their potential for high mobility applications, such as field effect transistors. InP nanowires exhibited the longest carrier lifetimes and the lowest surface recombination velocity of 170 cm s⁻¹. This very low surface recombination velocity makes InP nanowires suitable for applications where carrier lifetime is crucial, such as in photovoltaics. In contrast, the carrier lifetimes in GaAs nanowires were extremely short, of the order of picoseconds, due to the high surface recombination velocity, which was measured as 5.4 × 10⁵ cm s⁻¹. These findings will assist in the choice of nanowires for different applications, and identify the challenges in producing nanowires suitable for future electronic and optoelectronic devices.

Flow field measurements of a series of turbulent premixed and stratified methane/air flames

This paper presents flow field measurements for the turbulent stratified burner introduced in two previous publications in which high resolution scalar measurements were made by Sweeney et al. [1,2] for model validation. The flow fields of the series of premixed and stratified methane/air flames are investigated under turbulent, globally lean conditions (φg=0.75). Velocity data acquired with laser Doppler anemometry (LDA) and particle image velocimetry (PIV) are presented and discussed. Pairwise 2-component LDA measurements provide profiles of axial velocity, radial velocity, tangential velocity and corresponding fluctuating velocities. The LDA measurements of axial and tangential velocities enable the swirl number to be evaluated and the degree of swirl characterized. Power spectral density and autocorrelation functions derived from the LDA data acquired at 10kHz are optimized to calculate the integral time scales. Flow patterns are obtained using a 2-component PIV system operated at 7Hz. Velocity profiles and spatial correlations derived from the PIV and LDA measurements are shown to be in very good agreement, thus offering 3D mapping of the velocities. A strong correlation was observed between the shape of the recirculation zones above the central bluff body and the effects of heat release, stoichiometry and swirl. Detailed analyses of the LDA data further demonstrate that the flow behavior changes significantly with the levels of swirl and stratification, which combines the contributions of dilatation, recirculation and swirl. Key turbulence parameters are derived from the total velocity components, combining axial, radial and tangential velocities. © 2013 The Combustion Institute.

Effects of the biodiesel fuel physical properties on the swirl stabilised spray combustion characteristics

An increasin g interest in biofuel applications in modern engines requires a better understanding of biodiesel combustion behaviour. Many numerical studies have been carried out on unsteady combustion of biodiesel in situations similar to diesel engines, but very few studies have been done on the steady combustion of biodiesel in situations similar to a gas turbine combustor environment. The study of biodiesel spray combustion in gas turbine applications is of special interest due to the possible use of biodiesel in the power generation and aviation industries. In modelling spray combustion, an accurate representation of the physical properties of the fuel is a first important step, since spray formation is largely influenced by fuel properties such as viscosity, density, surface tension and vapour pressure. In the present work, a calculated biodiesel properties database based on the measured composition of Fatty Acid Methyl Esters (FAME) has been implemented in a multi-dimensional Computational Fluid Dynamics (CFD) spray simulation code. Simulations of non-reacting and reacting atmospheric-pressure sprays of both diesel and biodiesel have been carried out using a spray burner configuration for which experimental data is available. A pre-defined droplet size probability density function (pdf) has been implemented together with droplet dynamics based on phase Doppler anemometry (PDA) measurements in the near-nozzle region. The gas phase boundary condition for the reacting spray cases is similar to that of the experiment which employs a plain air-blast atomiser and a straight-vane axial swirler for flame stabilisation. A reaction mechanism for heptane has been used to represent the chemistry for both diesel and biodiesel. Simulated flame heights, spray characteristics and gas phase velocities have been found to compare well with the experimental results. In the reacting spray cases, biodiesel shows a smaller mean droplet size compared to that of diesel at a constant fuel mass flow rate. A lack of sensitivity towards different fuel properties has been observed based on the non-reacting spray simulations, which indicates a need for improved models of secondary breakup. By comparing the results of the non-reacting and reacting spray simulations, an improvement in the complexity of the physical modelling is achieved which is necessary in the understanding of the complex physical processes involved in spray combustion simulation. Copyright © 2012 SAE International.

Sand column impact onto a Kolsky pressure bar

A laboratory-based methodology to launch cylindrical sand slugs at high velocities is developed. The methodology generates well-characterised soil ejecta without the need for detonation of an explosive; this laboratory-based tool thereby allows for the experimental investigation of the soil-structure events. The experimental set-up comprises a launcher with a cylindrical cavity and a piston to push out the sand slug. The apparatus is used to launch both dry and water-saturated sand slugs. High speed photography is used to characterise the evolution of the sand slugs after launch. We find that the diameter of the slugs remains unchanged, and the sand particles possess only an axial component of velocity. However, the sand particles have a uniform spatial gradient of axial velocity and this results in lengthening of the slugs as they travel towards their target. Thus, the density of the sand slugs remains spatially homogenous but decreases with increasing time. The velocity gradient is typically higher in the dry sand slugs than that of the water-saturated slugs. The pressure exerted by the slugs on a rigid-stationary target is measured by impacting the slugs against a direct impact Kolsky bar. After an initial high transient pressure, the pressure reduces to a value of approximately ρv 2 where ρ is the density of the impacting sand slug and v is the particle velocity. This indicates that loading due to the sand is primarily inertial in nature. The momentum transmitted to the Kolsky bar was approximately equal to the incident momentum of the sand slugs, regardless of whether they are dry or water-saturated. © 2013 Elsevier Ltd. All rights reserved.