Instabilities in the oscillatory flow of a complex fluid.

The dynamics of a fluid in a vertical tube, subjected to an oscillatory pressure gradient, is studied experimentally for both a Newtonian and a viscoelastic shear-thinning fluid. Particle image velocimetry is used to determine the two-dimensional velocity fields in the vertical plane of the tube axis, in a range of driving amplitudes from 0.8 to 2.5 mm and of driving frequencies from 2.0 to 11.5 Hz. The Newtonian fluid exhibits a laminar flow regime, independent of the axial position, in the whole range of drivings. For the complex fluid, instead, the parallel shear flow regime exhibited at low amplitudes [Torralba, Phys. Rev. E 72, 016308 (2005)] becomes unstable at higher drivings against the formation of symmetric vortices, equally spaced along the tube. At even higher drivings the vortex structure itself becomes unstable, and complex nonsymmetric structures develop. Given that inertial effects remain negligible even at the hardest drivings (Re < 10(-1)), it is the complex rheology of the fluid that is responsible for the instabilities observed. The system studied represents an interesting example of the development of shear-induced instabilities in nonlinear complex fluids in purely parallel shear flow.

Experimental demonstration of the Rayleigh acoustic viscous boundary layer theory.

Amplitude and phase velocity measurements on the laminar oscillatory viscous boundary layer produced by acoustic waves are presented. The measurements were carried out in acoustic standing waves in air with frequencies of 68.5 and 114.5 Hz using laser Doppler anemometry and particle image velocimetry. The results obtained by these two techniques are in good agreement with the predictions made by the Rayleigh viscous boundary layer theory and confirm the existence of a local maximum of the velocity amplitude and its expected location.

2D PIV measurements in the near field of grid turbulence using stitched fields from multiple cameras

We present measurements of grid turbulence using 2D particle image velocimetry taken immediately downstream from the grid at a Reynolds number of Re M = 16500 where M is the rod spacing. A long field of view of 14M x 4M in the down- and cross-stream directions was achieved by stitching multiple cameras together. Two uniform biplanar grids were selected to have the same M and pressure drop but different rod diameter D and crosssection. A large data set (10 4 vector fields) was obtained to ensure good convergence of second-order statistics. Estimations of the dissipation rate ε of turbulent kinetic energy (TKE) were found to be sensitive to the number of meansquared velocity gradient terms included and not whether the turbulence was assumed to adhere to isotropy or axisymmetry. The resolution dependency of different turbulence statistics was assessed with a procedure that does not rely on the dissipation scale η. The streamwise evolution of the TKE components and ε was found to collapse across grids when the rod diameter was included in the normalisation. We argue that this should be the case between all regular grids when the other relevant dimensionless quantities are matched and the flow has become homogeneous across the stream. Two-point space correlation functions at x/M = 1 show evidence of complex wake interactions which exhibit a strong Reynolds number dependence. However, these changes in initial conditions disappear indicating rapid cross-stream homogenisation. On the other hand, isotropy was, as expected, not found to be established by x/M = 12 for any case studied. © Springer-Verlag 2012.

Statistics and dynamics of turbulence-flame alignment in premixed combustion

The geometric alignment of turbulent strain-rate structures with premixed flames greatly influences the results of the turbulence-flame interaction. Here, the statistics and dynamics of this alignment are experimentally investigated in turbulent premixed Bunsen flames using high-repetition-rate stereoscopic particle image velocimetry. In all cases, the statistics showed that the most extensive principal strain-rate associated with the turbulence preferentially aligned such that it was more perpendicular than parallel to the flame surface normal direction. The mean turbulence-flame alignment differed between the flames, with the stronger flames (higher laminar flame speed) exhibiting stronger preferential alignment. Furthermore, the preferential alignment was greatest on the reactant side of the mean flame brush. To understand these differences, individual structures of fluid-dynamic strain-rate were tracked through time in a Lagrangian manner (i.e., by following the fluid elements). It was found that the flame surface affected the orientation of the turbulence structures, with the majority of structures rotating as they approached the flame such that their most extensive principal strain-rate was perpendicular to the flame normal. The maximum change in turbulent structure orientation was found to decrease with the strength of the structure, increase with the strength of the flame, and exhibit similar trends when the structure strength and flame strength were represented by a Karlovitz number. The mean change in orientation decreased from the unburnt to burnt side of the flame brush and appears to be influenced by the overall flame shape. © 2011 The Combustion Institute.

A new apparatus for modelling excavations

Underground space is commonly exploited both to maximise the utility of costly land in urban development and to reduce the vertical load acting on the ground. Deep excavations are carried out to construct various types of underground infrastructure such as deep basements, subways and service tunnels. Although the soil response to excavation is known in principle, designers lack practical calculation methods for predicting both short- and long-term ground movements. As the understanding of how soil behaves around an excavation in both the short and long term is insufficient and usually empirical, the judgements used in design are also empirical and serious accidents are common. To gain a better understanding of the mechanisms involved in soil excavation, a new apparatus for the centrifuge model testing of deep excavations in soft clay has been developed. This apparatus simulates the field construction sequence of a multi-propped retaining wall during centrifuge flight. A comparison is given between the new technique and the previously used method of draining heavy fluid to simulate excavation in a centrifuge model. The new system has the benefit of giving the correct initial ground conditions before excavation and the proper earth pressure distribution on the retaining structures during excavation, whereas heavy fluid only gives an earth pressure coefficient of unity and is unable to capture any changes in the earth pressure coefficient of soil inside the zone of excavation, for example owing to wall movements. Settlements of the ground surface, changes in pore water pressure, variations in earth pressure, prop forces and bending moments in the retaining wall are all monitored during excavation. Furthermore, digital images taken of a cross-section during the test are analysed using particle image velocimetry to illustrate ground deformation and soil-structure interaction mechanisms. The significance of these observations is discussed.

The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows

A series of flames in a turbulent methane/air stratified swirl burner is presented. The degree of stratification and swirl are systematically varied to generate a matrix of experimental conditions, allowing their separate and combined effects to be investigated. Non-swirling flows are considered in the present paper, and the effects of swirl are considered in a companion paper (Part II). A mean equivalence ratio of φ=0.75 is used, with φ for the highest level of stratification spanning 0.375-1.125. The burner features a central bluff-body to aid flame stabilization, and the influence of the induced recirculation zone is also considered. The current work focuses on non-swirling flows where two-component particle image velocimetry (PIV) measurements are sufficient to characterize the main features of the flow field. Scalar data obtained from Rayleigh/Raman/CO laser induced fluorescence (CO-LIF) line measurements at 103μm resolution allow the behavior of key combustion species-CH 4, CO 2, CO, H 2, H 2O and O 2-to be probed within the instantaneous flame front. Simultaneous cross-planar OH-PLIF is used to determine the orientation of the instantaneous flame normal in the scalar measurement window, allowing gradients in temperature and progress variable to be angle corrected to their three dimensional values. The relationship between curvature and flame thickness is investigated using the OH-PLIF images, as well as the effect of stratification on curvature.The main findings are that the behavior of the key combustion species in temperature space is well captured on the mean by laminar flame calculations regardless of the level of stratification. H 2 and CO are significant exceptions, both appearing at elevated levels in the stratified flames. Values for surface density function and by extension thermal scalar dissipation rate are found to be substantially lower than laminar values, as the thickening of the flame due to turbulence dominates the effect of increased strain. These findings hold for both premixed and stratified flames. The current series of flames is proposed as an interesting if challenging set of test cases for existing and emerging turbulent flame models, and data are available on request. © 2012 The Combustion Institute.

Development of a new apparatus for modeling deep excavation related problems in geotechnical centrifuge

Underground space is commonly exploited both to maximise the utility of costly land in urban development and to reduce the vertical load acting on the ground. Deep excavations are carried out to construct various types of underground infrastructure such as deep basements, subways and service tunnels. Although the soil response to excavation is known in principle, designers lack practical calculation methods for predicting both short- and long-term ground movements. As the understanding of how soil behaves around an excavation in both the short and long term is insufficient and usually empirical, the judgements used in design are also empirical and serious accidents are common. To gain a better understanding of the mechanisms involved in soil excavation, a new apparatus for the centrifuge model testing of deep excavations in soft clay has been developed. This apparatus simulates the field construction sequence of a multi-propped retaining wall during centrifuge flight. A comparison is given between the new technique and the previously used method of draining heavy fluid to simulate excavation in a centrifuge model. The new system has the benefit of giving the correct initial ground conditions before excavation and the proper earth pressure distribution on the retaining structures during excavation, whereas heavy fluid only gives an earth pressure coefficient of unity and is unable to capture any changes in the earth pressure coefficient of soil inside the zone of excavation, for example owing to wall movements. Settlements of the ground surface, changes in pore water pressure, variations in earth pressure, prop forces and bending moments in the retaining wall are all monitored during excavation. Furthermore, digital images taken of a cross-section during the test are analysed using particle image velocimetry to illustrate ground deformation and soil–structure interaction mechanisms. The significance of these observations is discussed.

Large-eddy simulations and measurements of a small-scale high-speed coflowing jet

Measurements and predictions are made of a short-cowl coflowing jet with a bypass ratio of 8:1. The Reynolds number is 300,000, and the inlet Mach numbers are representative of aeroengine conditions. The low Reynolds number of the measurements makes the case well suited to the assessment of large-eddy-simulation-related strategies. The nozzle concentricity is carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both laser Doppler anemometry and particle image velocimetry. The simulations are completed on 6× 106, 12× 106, and 50 × 106 cell meshes. To overcome near-wall modeling problems, a hybrid large-eddy-simulation-Reynolds-averaged-Navier-Stokesrelated method is used. The near-wall Reynolds-averaged-Navier-Stokes layer is helpful in preventing nonphysical separation from the nozzle wall.Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Lift and the leading-edge vortex

Flapping wings often feature a leading-edge vortex (LEV) that is thought to enhance the lift generated by the wing. Here the lift on a wing featuring a leading-edge vortex is considered by performing experiments on a translating flat-plate aerofoil that is accelerated from rest in a water towing tank at a fixed angle of attack of 15°. The unsteady flow is investigated with dye flow visualization, particle image velocimetry (PIV) and force measurements. Leading-and trailing-edge vortex circulation and position are calculated directly from the velocity vectors obtained using PIV. In order to determine the most appropriate value of bound circulation, a two-dimensional potential flow model is employed and flow fields are calculated for a range of values of bound circulation. In this way, the value of bound circulation is selected to give the best fit between the experimental velocity field and the potential flow field. Early in the trajectory, the value of bound circulation calculated using this potential flow method is in accordance with Kelvin's circulation theorem, but differs from the values predicted by Wagner's growth of bound circulation and the Kutta condition. Later the Kutta condition is established but the bound circulation remains small; most of the circulation is contained instead in the LEVs. The growth of wake circulation can be approximated by Wagner's circulation curve. Superimposing the non-circulatory lift, approximated from the potential flow model, and Wagner's lift curve gives a first-order approximation of the measured lift. Lift is generated by inertial effects and the slow buildup of circulation, which is contained in shed vortices rather than bound circulation. © 2013 Cambridge University Press.

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.

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.

Comparison of acoustic velocity perturbation measurements using PIV vs. two-microphone technique

Understanding combustion instabilities requires accurate measurements of the acoustic velocity perturbation into injectors. This is often accomplished via the use of the two microphone technique, as this only requires two pressure transducers. However, measurements of the actual velocities emerging from the injectors are not often taken, leaving questions regarding the assumptions about the acoustic velocity. A comparison of velocity measured at downstream of the injector with that of two-microphone technique can show the accuracy and limitations of two-microphone technique. In this paper, velocity measurements are taken using both particle image velocimetry (PIV) and the two-microphone technique in a high pressure facility designed for aeroengine injector measurements. The flow is excited using an area modulation device installed on the choked end of the combustion chamber, with PIV measurements enabled by optical access downstream of the injector through a quartz tube and windows. Acoustic velocity perturbations at the injector are determined by considering the Fourier transformed pressure fluctuations for two microphones installed at a known distance upstream of the injector. PIV measurements are realized by seeding the air flow with micrometric water particles under 2.5 bar pressure at ambient temperature. Phase locked velocity fields are realized by synchronizing the acquisition of PIV images with the revolution of the acoustic modulator using the pressure signal measured at the face of injector. The mean velocity fluctuation is calculated as the difference between maximum and minimum velocities, normalized by the mean velocity of the unforced case. Those values are compared with the peak-to-peak velocity fluctuation amplitude calculated by the two-microphone technique. Although the ranges of velocity fluctuations for both techniques are similar, the variation of fluctuation with forcing frequencies diverges significantly with frequency. The differences can be attributed to several limitations associated with of both techniques, such as the quality of the signal, the signal/noise ratio, the accuracy of PIV measurements and the assumption of isentropic flow of the particle velocity from the plenum through the injector. We conclude that two-microphone methods can be used as a reference value for the velocity fluctuation in low order applications such as flame transfer functions, but not for drawing conclusions regarding the absolute velocity fluctuations in the injector. Copyright © 2013 by ASME.

Displacement of particles in microfluidics by laser-generated tandem bubbles.

The dynamic interaction between laser-generated tandem bubble and individual polystyrene particles of 2 and 10 μm in diameter is studied in a microfluidic channel (25 μm height) by high-speed imaging and particle image velocimetry. The asymmetric collapse of the tandem bubble produces a pair of microjets and associated long-lasting vortices that can propel a single particle to a maximum velocity of 1.4 m∕s in 30 μs after the bubble collapse with a resultant directional displacement up to 60 μm in 150 μs. This method may be useful for high-throughput cell sorting in microfluidic devices.

Comparison of Experimental PIV Data and CFD Simulations for Flow in a Diesel Particulate Filter Inlet Diffuser

Flow maldistribution of the exhaust gas entering a Diesel Particulate Filter (DPF) can cause uneven soot distribution during loading and excessive temperature gradients during the regeneration phase. Minimising the magnitude of this maldistribution is therefore an important consideration in the design of the inlet pipe and diffuser, particularly in situations where packaging constraints dictate bends in the inlet pipe close to the filter, or a sharp diffuser angle. This paper describes the use of Particle Image Velocimetry (PIV) to validate a Computational Fluid Dynamic (CFD) model of the flow within the inlet diffuser of a DPF so that CFD can be used with confidence as a tool to minimise this flow maldistribution. PIV is used to study the flow of gas into a DPF over a range of steady state flow conditions. The distribution of flow approaching the front face of the substrate was of particular interest to this study. Optically clear diffusing cones were designed and placed between pipe and substrate to allow PIV analysis to take place. Stereoscopic PIV was used to eliminate any error produced by the optical aberrations caused by looking through the curved wall of the inlet cone. In parallel to the experiments, numerical analysis was carried out using a CFD program with an incorporated DPF model. Boundary conditions for the CFD simulations were taken from the experimental data, allowing an experimental validation of the numerical results. The CFD model incorporated a DPF model, the cement layers seen in segmented filters and the intumescent matting that is commonly used to pack the filter into a metal casing. The mesh contained approximately 580,000 cells and used the realizable ?-e turbulence model. The CFD simulation predicted both pressure drop across the DPF and the velocity field within the cone and at the DPF face with reasonable accuracy, providing confidence in the use the CFD in future work to design new, more efficient cones.