The application of classical POD and snapshot POD in a turbulent shear layer with periodic structures

The classical and snapshot proper-orthogonal-decomposition was applied to data taken in a one-stream mixing layer in a narrow channel. Due to this particular geometry the flow develops large periodic structures. POD-analysis of simultaneously measured velocity components in spanwise direction identify as largest mode not only their periodic fraction, but also higher Fourier modes of the two-dimensional fluctuation. The energy content of the plane motion reaches values of about 90%. The amplitude of small three-dimensional vortices embedded in higher POD modes is correlated with the phase of the large structures, which indicates their influence on the entire turbulent motion. Application of scalar snapshot POD on phase averaged data of the entire flow field allows separation into modes. The eigenvalues and eigenvectors show identical distribution for theu- andv-component. Comparison of streakline plots of the reconstructed velocity field from different numbers of modes with flow visualization exhibits that the largest physical structure is described by only the first two modes. This is also supported by calculation of the vorticity component inz-direction. The total energy content of the largest structure is approximately 60%.Supported by the Deutsche Forschungsgemeinschaft under Grant No. Fi 178/28.     

Near Field Round Jet Flow Downstream from an Abrupt Contraction Nozzle with Tube Extension

Round air jet development downstream from an abrupt contraction coupled to a uniform circular tube extension with length to diameter ratio L/D?=?1.2 and L/D?=?53.2 is studied experimentally. Smoke visualisation and systematic hot film velocity measurements are performed for low to moderate Reynolds numbers 1130?<?Re _b ?<?11320. Mean and turbulent velocity profiles are quantified at the tube exit and along the centerline from the tube exit down to 20 times the diameter D. Flow development is seen to be determined by the underlying jet structure at the tube exit which depends on Reynolds number, initial velocity statistics at the tube exit and the presence/absence of coherent structures. It is shown that the tube extension ratio L/D as well as the sharp edged abrupt contraction influence the initial jet structure at the tube exit. For both L/D ratios, the presence of the abrupt contraction results in transitional jet flow in the range 2000?<?Re _b ?<?4000 and in flow features associated with forced jets and high Reynolds numbers Re _b ?>?10⁴. The tube extension ratio L/D downstream from the abrupt contraction determines the shear layer roll up so that for L/D?=?1.2 flow visualisation suggests the occurrence of toroidal vortices for Re _b ?<?4000 whereas helical vortices are associated with the transitional regime for L/D?=?53.2. Found flow features are compared to features reported in literature for smooth contraction nozzles and long pipe flow.     

Steady flow past a torus with aspect ratio less than 5

Steady flow past a torus with an aspect ratio less than 5 and its axis aligned with the flow is studied numerically by solving the steady, axisymmetric Navier–Stokes equations. The wake structure behind tori exhibits diverse behaviours. The detached recirculating zone on the axis, the attached recirculating zone, and the detached recirculating zone behind the torus tube may appear individually or concurrently, depending on the aspect ratio and the Reynolds number. A wake structure map is summarized based on the observed flow behaviours. Six flow regimes with different wake behaviours are identified and the corresponding flow regime map is plotted, which include the no-recirculating-zone regime, the single-detached-recirculating-zone regime, the single-attached-recirculating-zone regime, the two-recirculating-zone regime I, the two-recirculating-zone regime II, and the three-recirculating-zone regime. Over the range of aspect ratio 1.9<AR<2.4, the detached wake initially increases but then decreases in size with Reynolds number, and eventually disappears at Reynolds numbers beyond a critical value (depending on the aspect ratio). The underlying mechanisms of the onset and disappearance of the recirculating zones are discussed in terms of vorticity accumulation and base bleed. The recirculating zone first occurs when the maximum vorticity on the surface of the torus exceeds about 5. The detached recirculating zone on the axis of the torus disappears once the flow rate through the hole of the torus is beyond a certain threshold. In addition, the present results suggest that different transition modes to non-axisymmetric flow for tori with different aspect ratios reported in the literature may result from the wake structures prior to the transition.     

The role of shape and relative submergence on the structure of wakes of low-aspect-ratio wall-mounted bodies

The effects of shape and relative submergence (the ratio of flow depth to obstacle height, d/H) are investigated on the wakes around four different low-aspect-ratio wall-mounted obstacles at Re _H  = 17,800: semi-ellipsoids with the major axes of the base ellipses aligned in the streamwise and transverse directions, and two cylinders with aspect ratios matching the ellipsoids (H/D = 0.89 and 0.67, where D is the maximum transverse dimension). Particle Image Velocimetry was used to interrogate the flow. Streamwise features observed in the mean wake include counter-rotating distributions of vorticity inducing downwash (tip structures), upwash (base structures), and horseshoe vortices. In particular, the relatively subtle change in geometry produced by the rotation of the ellipsoid from the streamwise to the transverse orientation results in a striking modification of the mean streamwise vorticity distribution in the wake. Tip structures are dominant in the former case, while base structures are dominant in the latter. A vortex skeleton model of the wake is proposed in which arch vortex structures, shed from the obstacle, are deformed by the competing mechanisms of Biot-Savart self-induction and the external shear flow. The selection of tip or base structures in the ellipsoid wakes is caused by tilting of the arch structures either upstream or downstream, respectively, which is governed by ellipsoid curvature. An inverse relationship was observed between the relative submergence and the strength of the base structures for the ellipsoids, with a dominant base structure observed for d/H = 1 in both cases. These results demonstrate a means by which to achieve significant modifications to flow structure and thereby also to transport mechanisms in the flow. Therefore, this work provides insight into the modeling and control of flow over wall-mounted bodies.     

Instabilities in viscoelastic flow through an axisymmetric sudden contraction

Kinematics and dynamics of the viscoelastic flow in an axisymmetric 4 : 1 sudden contraction geometry are studied for a highly elastic polyisobutylene (PIB) based polymer solution (referred to as PIB-Boger fluid). The critical conditions for the onset of the elastic instabilities and the dynamics of the resulting secondary flows are measured for various flow rates. The spatio-temporal characteristics of the flow are determined by instantaneous pressure measurements and streakline photography. The nonlinear dynamics of the global flow field both upstream and downstream of the contraction plane are systematically examined. New dynamic flow behavior and elastic instabilities downstream of the contraction plane are reported. It is shown that the instantaneous pressure measurements along with flow visualization can be used as an effective tool to characterize viscoelastic flows in complex geometries.     

Experimental investigation of the wave-induced flow around a surface-touching cylinder

The wave-induced flow around a circular cylinder near both a rigid wall and an erodible bed is experimentally investigated using Particle Tracking Velocimetry (PTV). The aim of this study is to gain quantitative information on the local mean flow, the vorticity dynamics and the evolution of the erodible bed. The flow is characterized in terms of the Keulegan–Carpenter (KC), Reynolds (Re) and Ursell (U_r) numbers. The effects of changing these parameters over the ranges 1<KC<31, 3×10³<Re<2.6×10⁴ and 1.5<U_r<152 are investigated. For KC<1.1 the flow does not separate. When KC increases, the flow becomes unstable and large-scale vortical structures develop. The dimensionless intensity (|Γ^⁎|) depends non-monotonically on KC, with a local maximum at KC=17, and the dimensionless area of the same macrovortex (A^⁎) follows a somewhat similar law. Although the dimensionless boundary layer thickness (δ^⁎) exhibits some discontinuities between KC regimes, it decreases with KC at x/D=0.5, as x/D=1 weakly depends on KC and can be regarded as constant (δ^⁎=0.7) and then, increases with KC when moving away from the cylinder. These findings are used to interpret the physics governing the flow around a cylinder touching a wall and are compared with available results from the literature (Sumer et al., 1991). The evolution of the scour mechanism occurring over an erodible sandy bed is also investigated. The validity of some empirical formulas in the literature is also tested on the basis of the available dataset. The empirical relationships of Cevik and Yuksel (1999) and Sumer and Fredsøe (1990) for the dimensionless scour depth (S/D) agree well with our results. The dimensionless scour width (W_s/D) is predicted well by Sumer and Fredsøe's (2002) empirical equation for KC<23, whereas Catano-Lopera and Garcia's (2007) formula is more accurate for higher values of KC.     

Predicting pressure drop in open-cell foams by adopting Forchheimer number

Interconnected struts arranged in 3-D foam structures pose a challenge in understanding fluid flow, which is significantly different from that in traditional porous media. Different flow regimes (Darcy, transition and weak inertia regimes) and thus, different flow laws in open-cell foams are used. The impact of characteristic lengths’ choices based on both, morphological and hydraulic parameters on flow law formulation has been studied. Ambiguities in definitions and measurements of several key parameters have been shown and limitations in the use of some parameters have been pointed out.An equivalent Reynolds number in the form of Forchheimer number (F_o) has been proposed to establish the friction factor relationship in order to avoid any morphological ambiguities. This number takes into account hydraulic characteristics of viscous and inertia regimes simultaneously. It has been observed that when F_o < 0.1, the flow through open-cell foams remains in the Darcy regime while the occurrence of weak inertia regime dominates when F_o > 1. Transition regime occurs in a narrow range of flow velocity when 0.1 < F_o < 1. The limits of transition for regime identification are found to be independent of foam morphologies. The form drag coefficient varies in relation with foam morphological parameters and is not a “universal” constant.Empirical correlations have been derived to predict hydraulic characteristics and friction factor data for different strut shapes and porosities. An excellent agreement has been obtained between predicted and numerical/experimental flow data.     

Experimental study on a non-dilute two-phase coflowing jet: Dynamics of particles in the near flow field

The dynamics of particles in multi-phase jets has been widely studied due to its importance for a broad range of practical applications. The present work describes an experimental investigation on an initially non-dilute two-phase jet, aimed at improving the understanding in this field. A two-color PDPA has been employed to measure simultaneously the velocity and size of particles. The measurements are post-processed to check the reliability of the results and to derive information on particle volume flux as an indication of their concentration. Acoustic forcing is applied in order to control coherent structures, which are responsible for mixing and transport phenomena, and also to get phase-locked measurements. Phase-averaged statistics enabled to freeze the jet structure, not visible in the time-averaged data. The results along the jet centerline confirm that drag forces and the spread angle of the jet initially control particle dispersion, very near the nozzle exit (x/D < 4). However, as the vortical structures evolve forming tongue-shaped structures, the total particle volume flux is augmented when these structures connect with the main stream (x/D > 5). This is due to an increase of the number of smaller size particles, even when a decrease of the number of larger size particle is observed. Further analysis at five cross-stream sections across two consecutive vortices confirm that small particles are convected around the coherent structure and then incorporated to the main stream, increasing the particle concentration at the jet core. On the other hand, the number of larger particles (as well as their contribution to axial volume flux) starts to decay in regions of high azymuthal vorticity. This behaviour is partly ascribed to the transversal lift force, associated to the large spatial gradients observed in these regions. Saffman and Magnus forces have been estimated to be comparable or even greater than radial drag forces. The results suggest that the Saffman force might accelerate particles in radial direction, inducing a high radial volumetric flow rate from high to low axial velocity regions.     

Direct numerical simulation of a turbulent Couette-Poiseuille flow,part 2: Large- and very-large-scale motions

A direct numerical simulation dataset of a fully developed turbulent Couette-Poiseuille flow is analyzed to investigate the spatial organization of streamwise velocity-fluctuating u-structures on large and very large scales. Instantaneous and statistical flow fields show that negative-u structures with a small scale on a stationary bottom wall grow throughout the centerline due to the continuous positive mean shear, and they penetrate to the opposite moving wall. The development of an initial vortical structure related to negative-u structures on the bottom wall into a large-scale hairpin vortex packet with new hairpin vortices, which are created upstream and close to the wall, is consistent with the auto-generation process in a Poiseuille flow (Zhou et al., J. Fluid Mech., vol. 387, 1999, pp. 353–396). Although the initial vortical structure associated with positive-u structures on the top wall also grows toward the bottom wall, the spatial development of the structure is less coherent with weak strength due to the reduced mean shear near the top wall, resulting in less turbulent energy on the top wall. The continuous growth of the structures from a wall to the opposite wall explains the enhanced wall-normal transport of the streamwise turbulent kinetic energy near the centerline. Finally, an inspection of the time-evolving instantaneous fields and conditional averaged flow fields for the streamwise growth of a very long structure near the centerline exhibits that a streamwise concatenation of adjacent large-scale u-structures creates a very-large-scale structure near the channel centerline.     

Coherent vortical and straining structures in the finite wall-mounted square cylinder wake

Topological aspects of the turbulent wake of a finite, surface-mounted, square-cross-section cylinder of h/d = 4 are addressed by decomposing the velocity field into a quasi-periodic coherent part and the unresolved incoherent fluctuations. The three-dimensional large scale structure is educed through a reconstruction of planar phase-averaged PIV measurements using the simultaneously sampled surface pressure difference on opposing sides of the obstacle as a phase reference. A topological model for the vortex structure is educed and mean streamwise wake vorticity is explained in terms of the connections between initially vertical structures shed alternately from either side of the obstacle, rather than previously proposed ‘tip’ vortex structures generated at the obstacle free-end. The coherent structure educed accounts for a significant portion of the fluctuating energy in the wake. The turbulent field is further analyzed by finding Lagrangian straining structures that form by induction of the coherent vorticity field, and these structures are related to the energy transfer from the base phase-averaged flow since they act to stretch incoherent vorticity fluctuations in their neighbourhood.     

Finite deformation and stability behaviour of cylindrical membranes subjected to symmetric line loads2

An ‘exact’ analysis of the complete non-linear load-deflection and stability behaviour of cylindrical membranes without end ‘shear walls’ subjected to longitudinal Symmetric line loads is presented. The analysis includes low as well as high profile structures. In order to determine the lateral stability behaviour, infinitesimal lateral displacements are superimposed on the symmetric finite deflection field.Results indicate that such structures may experience vertical and/or lateral instability depending on their initial geometry. Membranes with initial central angles, θ₀ ? 90° are stable, both vertically and laterally, for all load values. For 42.23° < θ₀ < 90° the structure is laterally stable but becomes vertically unstable at a certain ‘limit load’ W?_V while for θ₀< 42.23°, the structure becomes laterally unstable at a load value w?_l<W?_V.The analysis admits contact between vertical segments of the membrane to either side of the line load as well as between the membrane and the horizontal surface next to the supports.     

Two-phase flow laden with spherical particles in a microcapillary

Solid–liquid two-phase flow in a finite Reynolds number range (2 < Re < 12), transporting neutrally-buoyant microspheres with diameters of 6, 10, and 16 μm through a 260-μm microcapillary, is investigated. A standard microparticle-tracking velocimetry (μ-PTV) that consists of a double-pulsed Nd:YAG laser, an epi-fluorescent microscope, and a cooled-CCD camera is used to examine the flow. The solid particles are visualized in view of their spatial distributions. We observe a strong radial migration of the particles across the flow streamlines at substantially small Re. The degree of particle migration is presented in terms of probability density function. Some applications based on this radial migration phenomena are discussed in conjunction with particle separation/concentration in microfluidic devices, where the spatial distribution of particles is of great importance. In doing so, we propose a particle-trajectory function to empirically construct the spatial distribution of solid particles, which is well correlated with our experimental data. It is believed that this function provides a simple method for estimating the spatial distribution of particles undergoing radial migration in solid–liquid two-phase flows.     

Bubble migration in a turbulent boundary layer

The wall void peaking distribution observed in an upward turbulent bubbly boundary layer along a flat plate is generated by bubbles that move towards the plate, come into contact with the wall and then slide along it. This transverse ‘migration’ has been studied using flow visualization, high speed video and particle tracking techniques to measure the trajectories of mono-disperse air bubbles at very low void fractions. Investigations have been performed at four Reynolds numbers in the range 280 < Re_θ < 3000, covering both the laminar and turbulent regimes, with mono-disperse bubbles of mean equivalent diameter between 2 mm and 6 mm. Lagrangian statistics calculated from hundreds of trajectories show that the migration only occurs in the turbulent regime and for bubble diameters below some critical value: 3.5 mm < d_eqcrit < 4 mm. Above this size (We > 3), the interface deformation is such that bubbles do not remain at the wall, even when they are released at the surface. Also, bubble migration is shown to be non-systematic, to have a non-deterministic character in the sense that trajectories differ significantly, to increase with Reynolds number and to take place on a short time scale. A series of experiments with isolated bubbles demonstrates that these results are not influenced by bubble–bubble interactions and confirm that two-way coupling in the flow is limited. Flow visualizations show that the migration originates with the capture of bubbles inside the large turbulent structures of the boundary layer (‘bulges’). The bubbles begin to move towards the wall as they cross these structures, and the point at which they reach the wall is strongly correlated with the position of the deep ‘valleys’ which separate the turbulent ‘bulges’. The analysis of the mean Lagrangian trajectories of migrating bubbles confirms these observations. Firstly, the average time of migration calculated from these trajectories coincides with the mean transit time of the bubbles across the structures. Secondly, once the trajectories have been scaled by this transit time and the boundary layer thickness δ, they all have the same form in the region y/δ < 0.4, independent of the Reynolds number.     

Study of the Log-Layer Structure in Wall Turbulence Over a Very Large Range of Reynolds Number

Studies of the logarithmic layer structure in turbulent boundary layers are presented that span three orders of magnitude change in Reynolds number. The experiments considered used two separate laboratory scale facilities, as well as the atmospheric surface layer at the SLTEST facility in Utah. Several experimental techniques were used in order to probe the three-dimensional nature of the flow structures. The main focus is on two-point correlation statistics at a given z/δ, which are found to agree well over all Reynolds numbers when scaled with an outer length-scale. Large-scale coherence recently noted in the logarithmic region of laboratory-scale boundary layers is also found to be present in the atmospheric surface layer flow. Recent findings regarding the influence of these large scale motions on the near-wall region are also presented.     

Turbulence measurements in the bubbly flow region of hydraulic jumps

A hydraulic jump is characterized by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air–water flow measurements were performed in a large-size facility using two types of phase-detection intrusive probes: i.e. single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent integral time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of Chanson [H. Chanson, Bubbly flow structure in hydraulic jump, Eur. J. Mech. B/Fluids 26 (3) (2007) 367–384], providing the turbulent scales of the eddy structures advecting the air bubbles in the developing shear layer. The length scale L_xz is an integral air–water turbulence length scale which characterized the transverse size of the large vortical structures advecting the air bubbles. The experimental data showed that the dimensionless integral turbulent length scale L_xz/d₁ was closely related to the inflow depth: i.e. L_xz/d₁ = 0.2–0.8, with L_xz increasing towards the free-surface.     

On thermal convection between two coaxial square containers inclined by 45°

This paper is concerned with flow visualization, velocity, streakline, and heat transfer measurements on thermal convection between two coaxial square containers inclined by 45°. It is found that there is no critical Rayleigh number (Ra)_c for the onset of thermal convection in the test section. There appear two great circulations, which are symmetrical with respect to the diagonal vertical plane of the two coaxial square containers, and enhance the heat transfer significantly. Received on 23 October 1997     

Vortex dynamics and scalar transport in the wake of a flat-plate controlled by a vibrating trailing-edge flap

We investigate the onset and development of vortical flow disturbances introduced into the wake of a horizontally fixed flat-plate by means of the controlled motion of a trailing edge flap. The vibrating mechanics of the flap allows for the introduction of both impulsive and harmonic weak amplitude velocity disturbances which are propagated downstream into the wake flow of the flat-plate. Quantitative experimental and numerical predictions of both steady and unsteady wake flow velocity resulting from different flapping frequencies are made at low Reynolds numbers (Re < 10⁴). Frequency response tests of the wake confirmed the existence of two dominant frequencies where the wake flow organises with a particular arrangement of downstream moving vortex structures. Numerical predictions of steady (unforced) and forced wake velocity profiles and kinetic energy profiles are in good agreement with the experimental results. In order to understand practical implications of the dominant vortex structures in scalar transport, we have extended the numerical part of the study solving for the concentration equation of a passive scalar being injected in particular regions of the physical domain. A spatial correlation between the trajectory of vortex structures and the scalar concentration downstream the wake is observed. Moreover, the onset of tip vortex structures produced during the forcing cycle seems to be responsible of a local increase of scalar concentration near the span wise flap ends.     

Coherent Structures and their Frequency Signature in the Separated Shear Layer on the Sides of a Square Cylinder

The purpose of the present work is to study the various specific time scales of the turbulent separating flow around a square cylinder, in order to determine the Reynolds number effect on the separating shear layer, where occurs a transition to turbulence. Unsteady analysis based on large eddy simulation (LES) at intermediate Reynolds numbers and laser doppler velocimetry (LDV) measurements at high Reynolds numbers are carried out. The Reynolds number, based on the cylinder diameter D and the inflow velocity U _o , is ranging from Re?=?50 to Re?=?300,000. A special focus is performed on the coherent structures developing on the sides and in the wake of a square cylinder. For a large Reynolds number range above Re?≈?1,000, both signatures of Von Karman (VK) and Kelvin–Helmholtz (KH) type vortical structures are found on velocity time samples. The combination of their frequency signature is studied based on Fourier and wavelet analysis. In the present study, We observe the occurrence of KH pairings in the separating shear layer on the side of the cylinder, and confirm the intermittency nature of such a shear flow. These issues concerning the structure of the near wake shear layer which were addressed for the round cylinder case in a recent experimental publication (Rajagopalan and Antonia, Exp Fluids 38:393–402, 2005) are of interest in the present flow configuration as well.     

Effects of Permeable Cylinder on the Flow Structure in Deep Water

Flow behaviors around permeable cylinders were investigated using Particle Image Velocimetry technique in deep water. The height of deep water and free stream velocity were kept constant as h_w = 340 mm and U = 156 mm/s. To find out the effect of the permeable cylinders on the flow structure, eight different porosities (β = 0.4, 0.5, 0.6, 0.65, 0.7, 0.75, 0.8, and 0.85) were used. The results have indicated that the permeable cylinders are effective on the control of large-scale vortical structures downstream of the permeable cylinder. As the porosity increases, turbulent kinetic energy and Reynolds shear stress decrease. This means that the fluctuations in the wake region are significantly weakened by permeable cylinders. The permeable cylinders having the porosity higher than 0.6 do not pose an obstacle in the flow. Furthermore, for all diameter values of permeable cylinders, it can be concluded that the flow structures downstream of the permeable cylinder show similar trend with each other.     

Instabilities in viscoelastic flow past a square cavity

Elastic flow transitions in viscoelastic flow past a square cavity adjacent to a channel are reported. The critical conditions for the onset of flow transitions and the qualitative and quantitative characterization of the secondary flows generated by the instability have been examined using streakline photography and instantaneous pressure measurements. Cellular type of instabilities inside the cavity is observed for flow rates beyond a critical value. Small and large scale eddies are observed at high flow rates. The flow inside the cavity and in the channel upstream and downstream of the cavity becomes weakly time-dependent for high flow rates.