The effects of a second tunnel on the propagation of ground-borne vibration from an underground railway

Accurate predictions of ground-borne vibration levels in the vicinity of an underground railway are greatly sought in modern urban centers. Yet the complexity involved in simulating the underground environment means that it is necessary to make simplifying assumptions about this environment. One such commonly-made assumption is to model the railway as a single tunnel, despite many underground railway lines consisting of twin-bored tunnels. A unique model for two tunnels embedded in a homogeneous, elastic full space is developed. The vibration response of this two-tunnel system is calculated using the superposition of two displacement fields: one resulting from the forces acting on the invert of a single tunnel, and the other resulting from the interaction between the tunnels. By partitioning of the stresses into symmetric and anti-symmetric mode number components using Fourier decomposition, these two displacement fields can by calculated with minimal computational requirements. The significance of the interactions between twin-tunnels is quantified by calculating the insertion gains that result from the existence of a second tunnel. The insertion-gain results are shown to be localized and highly dependent on frequency, tunnel orientation and tunnel thickness. At some locations, the magnitude of these insertion gains is greater than 20dB. This demonstrates that a high degree of inaccuracy exists in any surface vibration-prediction model that includes only one of the two tunnels. © 2012 Springer.

Eigenmodes of lined flow ducts with rigid splices

This paper presents an analytic expression for the acoustic eigenmodes of a cylindrical lined duct with rigid axially running splices in the presence of flow. The cylindrical duct is considered to be uniformly lined except for two symmetrically positioned axially running rigid liner splices. An exact analytic expression for the acoustic pressure eigenmodes is given in terms of an azimuthal Fourier sum, with the Fourier coefficients given by a recurrence relation. Since this expression is derived using a Greens function method, the completeness of the expansion is guaranteed. A numerical procedure is described for solving this recurrence relation, which is found to converge exponentially with respect to number of Fourier terms used and is in practice quick to compute; this is then used to give several numerical examples for both uniform and sheared mean flow. An asymptotic expression is derived to directly calculate the pressure eigenmodes for thin splices. This asymptotic expression is shown to be quantitatively accurate for ducts with very thin splices of less than 1 % unlined area and qualitatively helpful for thicker splices of the order of 6 % unlined area. A thin splice is in some cases shown to increase the damping of certain acoustic modes. The influences of thin splices and thin boundary layers are compared and found to be of comparable magnitude for the parameters considered. Trapped modes at the splices are also identified and investigated. © 2011 Cambridge University Press.

Magnetically actuated artificial cilia: the effect of fluid inertia.

Natural cilia are hairlike microtubule-based structures that are able to move fluid on the micrometer scale using asymmetric motion. In this article, we follow a biomimetic approach to design artificial cilia lining the inner surfaces of microfluidic channels with the goal of propelling fluid. The artificial cilia consist of polymer films filled with superparamagnetic nanoparticles, which can mimic the motion of natural cilia when subjected to a rotating magnetic field. To obtain the magnetic field and associated magnetization local to the cilia, we solve the Maxwell equations, from which the magnetic body moments and forces can be deduced. To obtain the ciliary motion, we solve the dynamic equations of motion, which are then fully coupled to the Navier-Stokes equations that describe the fluid flow around the cilia, thus taking full account of fluid inertial forces. The dimensionless parameters that govern the deformation behavior of the cilia and the associated fluid flow are arrived at using the principle of virtual work. The physical response of the cilia and the fluid flow for different combinations of elastic, fluid viscous, and inertia forces are identified.

Influence of fluid viscosity on the response of buried structures in earthquakes

The use of high viscous pore fluid has been widely established to match the rate of excess pore pressure generation and subsequent dissipation in dynamic centrifuge tests. The appropriate viscosity is linked to the geometric and gravity scaling factors which corresponds to the use of pore fluid of 'N' cSt in a 'N'g centrifuge test. The use of either water (1 cSt) or pore fluid lower than 'N' cSt can influence the behaviour of soil liquefaction in a centrifuge test. In this paper, the floatation of a tunnel following soil liquefaction is investigated using pore fluids with two different viscosities. The results show that the uplift displacement of the tunnel is significantly affected by the pore fluid viscosity. © 2010 Taylor & Francis Group, London.

Surface tension-induced global instability of planar jets and wakes

The effect of surface tension on global stability of co-flow jets and wakes at a moderate Reynolds number is studied. The linear temporal two-dimensional global modes are computed without approximations. All but one of the flow cases under study are globally stable without surface tension. It is found that surface tension can cause the flow to be globally unstable if the inlet shear (or equivalently, the inlet velocity ratio) is strong enough. For even stronger surface tension, the flow is re-stabilized. As long as there is no change of the most unstable mode, increasing surface tension decreases the oscillation frequency. Short waves appear in the high-shear region close to the nozzle, and their wavelength increases with increasing surface tension. The critical shear (the weakest inlet shear at which a global instability is found) gives rise to antisymmetric disturbances for the wakes and symmetric disturbances for the jets. However, at stronger shear, the opposite symmetry can be the most unstable one, in particular for wakes at high surface tension. The results show strong effects of surface tension that should be possible to reproduce experimentally as well as numerically.

Characteristics of blood flow in microchannels and relevant impact on modelling blood behaviour in biochip separators

This paper reports a perspective investigation of computational modelling of blood fluid in microchannel devices as a preparation for future research on fluid-structure interaction (FSI) in biofluid mechanics. The investigation is carried out through two aspects, respectively on physical behaviours of blood flow in microchannels and appropriate methodology for modelling. The physics of blood flow is targeted to the challenges for describing blood flow in microchannels, including rheology of blood fluid, suspension features of red blood cells (RBCs), laminar hydrodynamic influence and effect of surface roughness. The analysis shows that due to the hyperelastic property of RBC and its comparable dimension with microchannels, blood fluid shows complex behaviours of two phase flow. The trajectory and migration of RBCs require accurate description of RBC deformation and interaction with plasma. Following on a discussion of modelling approaches, i.e. Eulerian method and Lagrangian method, the main stream modelling methods for multiphase flow are reviewed and their suitability to blood flow is analysed. It is concluded that the key issue for blood flow modelling is how to describe the suspended blood cells, modelled by Lagrangian method, and couple them with the based flow, modelled by Eulerian method. The multiphase flow methods are thereby classified based on the number of points required for describing a particle, as follows: (i) single-point particle methods, (ii) mutli-point particle methods, (iii) functional particle methods, and (iv) fluid particle methods. While single-point particle methods concentrate on particle dynamic movement, multipoint and functional particle methods can take into account particle mechanics and thus offer more detailed information for individual particles. Fluid particle methods provide good compromise between two phases, but require additional information for particle mechanics. For furthermore detailed description, we suggest to investigate the possibility using two domain coupling method, in which particles and base flow are modelled by two separated solvers. It is expected that this paper could clarify relevant issues in numerical modelling of blood flow in microchannels and induce some considerations for modelling blood flow using multiphase flow methods. © 2012 IEEE.

Geometry and interaction of structures in homogeneous isotropic turbulence

A strategy to extract turbulence structures from direct numerical simulation (DNS) data is described along with a systematic analysis of geometry and spatial distribution of the educed structures. A DNS dataset of decaying homogeneous isotropic turbulence at Reynolds number Reλ = 141 is considered. A bandpass filtering procedure is shown to be effective in extracting enstrophy and dissipation structures with their smallest scales matching the filter width, L. The geometry of these educed structures is characterized and classified through the use of two non-dimensional quantities, planarity' and filamentarity', obtained using the Minkowski functionals. The planarity increases gradually by a small amount as L is decreased, and its narrow variation suggests a nearly circular cross-section for the educed structures. The filamentarity increases significantly as L decreases demonstrating that the educed structures become progressively more tubular. An analysis of the preferential alignment between the filtered strain and vorticity fields reveals that vortical structures of a given scale L are most likely to align with the largest extensional strain at a scale 3-5 times larger than L. This is consistent with the classical energy cascade picture, in which vortices of a given scale are stretched by and absorb energy from structures of a somewhat larger scale. The spatial distribution of the educed structures shows that the enstrophy structures at the 5η scale (where η is the Kolmogorov scale) are more concentrated near the ones that are 3-5 times larger, which gives further support to the classical picture. Finally, it is shown by analysing the volume fraction of the educed enstrophy structures that there is a tendency for them to cluster around a larger structure or clusters of larger structures. Copyright © 2012 Cambridge University Press.

The rise heights of low-and high-Froude-number turbulent axisymmetric fountains

We present the results of an experimental investigation across a broad range of source Froude numbers, 0. 4 ≤ Fr 0 ≤ 45, into the dynamics, morphology and rise heights of Boussinesq turbulent axisymmetric fountains in quiescent uniform environments. Typically, these fountains are thought to rise to an initial height, z i, before settling back and fluctuating about a lesser (quasi-) steady height, z ss. Our measurements show that this is not always the case and the ratio of the fountain's initial rise height to steady rise height, λ = z i/z ss, varies widely, 0. 5 ≈ λ ≈ 2, across the range of Fr 0 investigated. As a result of near-ideal start-up conditions provided by the experimental set-up we were consistently able to form a vortex at the fountain's front. This enabled new insights into two features of the initial rise of turbulent fountains. Firstly, for 1. 0 ≈ Fr 0 ≈ 1. 7 the initial rise height is less than the steady rise height. Secondly, for Fr 0 ≈ 5. 5, the vortex formed at the fountain's front pinches off, separates from the main body and rises high above the fountain; there is thus a third rise height to consider, namely, the maximum vortex rise height, z v. From our observations we propose classifying turbulent axisymmetric fountains into five regimes (as opposed to the current three regimes) and present detailed descriptions of the flow in each. Finally, based on an analysis of the rise height fluctuations and the width of fountains in (quasi-) steady state we provide further insight into the physical cause of height fluctuations. © 2011 Cambridge University Press.

The night purging of a two-storey atrium building

We examine the fluid mechanics of night purging in a two-storey naturally ventilated atrium building. We develop a mathematical model of a simplified atrium building and focus on the rate at which warm air purges from each storey and the atrium by displacement ventilation into a still cool night environment of a constant temperature. To develop a first insight into how the geometry of the building influences the rate at which warm air purges from each storey via the atrium we neglect heat exchange with the fabric (so there is no thermal buffering) and furthermore assume that the warm air layers in each storey and the atrium are of uniform temperature. The plumes of warm air that rise from the storeys into the atrium, causing the atrium to fill with warm air, have a very strong influence on the night purge. Modelling these as axisymmetric turbulent plumes, we identify three forms of purging behaviour. Each purge is characterised by five key times identified in the progression of the night purge and physical rationale for these differing behaviours is given. An interface velocity deficit and volumetric purge deficit are introduced as measures of the efficiency of a night purge. © 2010 Elsevier Ltd.

Sensitivity analysis of a time-delayed thermo-acoustic system via an adjoint-based approach

We apply adjoint-based sensitivity analysis to a time-delayed thermo-acoustic system: a Rijke tube containing a hot wire. We calculate how the growth rate and frequency of small oscillations about a base state are affected either by a generic passive control element in the system (the structural sensitivity analysis) or by a generic change to its base state (the base-state sensitivity analysis). We illustrate the structural sensitivity by calculating the effect of a second hot wire with a small heat-release parameter. In a single calculation, this shows how the second hot wire changes the growth rate and frequency of the small oscillations, as a function of its position in the tube. We then examine the components of the structural sensitivity in order to determine the passive control mechanism that has the strongest influence on the growth rate. We find that a force applied to the acoustic momentum equation in the opposite direction to the instantaneous velocity is the most stabilizing feedback mechanism. We also find that its effect is maximized when it is placed at the downstream end of the tube. This feedback mechanism could be supplied, for example, by an adiabatic mesh. We illustrate the base-state sensitivity by calculating the effects of small variations in the damping factor, the heat-release time-delay coefficient, the heat-release parameter, and the hot-wire location. The successful application of sensitivity analysis to thermo-acoustics opens up new possibilities for the passive control of thermo-acoustic oscillations by providing gradient information that can be combined with constrained optimization algorithms in order to reduce linear growth rates. © Cambridge University Press 2013.