Direct numerical simulation of turbulent wake behind a surface-mounted square cylinder

In this research, direct numerical simulation has been performed to study the turbulent wake behind a wall-mounted square cylinder with aspect ratio 4 and Reynolds number 12 000 (based on the free-stream velocity and obstacle side length) in a developing boundary layer. Owing to the relatively high Reynolds number and high aspect ratio of the cylinder tested, the wake is wide spread behind the cylinder and exhibits complex and energetic vortex motions. The lateral and tip vortex shedding patterns at different frequencies, coherent structures downstream of the obstacle, the production rate and distribution of turbulent kinetic energy, and the instantaneous pressure distribution in the wake region have been thoroughly investigated. In order to validate the numerical results, the first- and second-order flow statistics obtained from the simulations have been carefully compared against available wind-tunnel measurement data.     

Surface pressure on a cubic building exerted by conical vortices

Experiments were performed to study surface pressure on a cubic building underlying conical vortices, which are known to cause severe structural damage and failure. The focus is on the effects of turbulence in the incident flow. Three turbulent boundary layers were created in a boundary layer wind tunnel. A wall-mounted cube, i.e. a cube situated on the horizontal ground floor surface of the wind-tunnel test section, was used as an experimental model. The cube was subjected to the incidence flow at 40°. Steady and unsteady pressure measurements were performed on the cube surface. The analysis suggests that conical vortices developed above the top surface of the wall-mounted cube. A larger mean suction was observed on the top cube surface in the less turbulent boundary layer. With an increase in turbulence in the incoming flow, the strong suction zones decreased in size. The fluctuating pressure coefficient profiles retained their shape when the turbulence in the upstream flow of the cube increased. The fluctuating pressure coefficient was observed to be larger in more turbulent flows. The pressure fluctuations were larger on the cube surface underlying outer boundaries of the conical vortex. The fluctuating pressure coefficient under the conical vortex was three to four times larger than in the weak suction zone on the central area of the top cube surface. Close to the leading cube corner, the pressure spectra were dominated by a single low frequency peak. As the conical vortex developed, this primary peak weakened and a secondary peak emerged at a higher reduced frequency. There is a general trend of shifting the pressure spectra towards higher reduced frequencies when the turbulence in the undisturbed incident flow increases.     

The turbulent wake of a circular cylinder rotating about the streamwise axis

This paper presents measurements in the turbulent wake of a circular cylinder rotating with its axis normal to the free-stream velocity; in other words, the axis of rotation was parallel to the streamwise direction. All three mean velocities and six Reynolds stresses were obtained at three positions downstream of the cylinder, with and without rotation of the free-stream. Most emphasis is given to the latter results because of the better flow quality. The ratio of the circumferential velocity of the cylinder to the free-stream velocity — the swirl number — had a maximum value of 0.6. Measurements for two combinations of the free-stream and angular velocities showed the velocity deficit in the wake to be a multi-valued function of the swirl number, implying that the rotation affected the separation of the cylinder's boundary layer in a complex manner. In the turbulent wake, the rotation did not significantly alter the magnitudes of the normal stresses, but caused large changes to the shape of the profiles of the axial and cross-stream normal stresses. Eventually, the primary (cross-stream) shear stress became almost entirely positive, but there was no corresponding change to the (cross-stream) gradient of the streamwise mean velocity. Despite these alterations to the turbulence, the rotationally-activated generation terms in the Reynolds transport equations never dominated the terms that are common to the wakes of rotating and non-rotating cylinders.This work was supported by the Australian Research Council. Most of the data acquisition software was written by Mr J. J. Smith.     

Turbulent flow around a wall-mounted cube: A direct numerical simulation

An immersed-boundary method was employed to perform a direct numerical simulation (DNS) of flow around a wall-mounted cube in a fully developed turbulent channel for a Reynolds number Re = 5610, based on the bulk velocity and the channel height. Instantaneous results of the DNS of a plain channel flow were used as a fully developed inflow condition for the main channel. The results confirm the unsteadiness of the considered flow caused by the unstable interaction of a horseshoe vortex formed in front of the cube and on both its sides with an arch-type vortex behind the cube. The time-averaged data of the turbulence mean-square intensities, Reynolds shear stresses, kinetic energy and dissipation rate are presented. The negative turbulence production is predicted in the region in front of the cube where the main horseshoe vortex originates.     

Lagrangian particle dispersion in turbulent flow over a wall mounted obstacle

Large-eddy simulations (LES) of particle-laden turbulent flows are presented in order to investigate the effects of particle response time on the dispersion patterns of a space developing flow with an obstruction, where solid particles are injected inside the wake of an obstacle [Vincont, J.Y., Simoens, S., Ayrault M., Wallace, J.M., 2000. Passive scalar dispersion in a turbulent boundary layer from a line source at the wall and downstream of an obstacle. J. Fluid Mech. 424, 127–167]. The numerical method is based on a fully explicit fractional step approach and finite-differences on Cartesian grids, using the immersed boundary method (IBM) to represent the existence of solid obstacles. Two different turbulence models have been tested, the classical Smagorinsky turbulence model and the filtered structure function model. The dispersed phase was modelled either by an Eulerian approach or a Lagrangian particle tracking scheme of solid particles with Stokes numbers in the range St = 0–25, assuming one-way coupling between the two phases. A very good agreement was observed between the Lagrangian and Eulerian approaches. The effect of particle size was found to significantly differentiate the dispersion pattern for the inhomogeneous flow over the obstacle. Although in homogeneous flows like particle-laden turbulent channels near-wall particle clustering increases monotonically with particle size, for the examined flow over an obstacle, preferential concentration effects were stronger only for an intermediate range of Stokes numbers.     

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.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

Comparison of flow structures in the downstream region of a cylinder and sphere

An experimental investigation of flow structures downstream of a circular cylinder and sphere immersed in a free-stream flow is performed for Re = 5000 and 10,000 using qualitative and quantitative flow visualization techniques. The obtained results are presented in terms of time-averaged velocity vectors, patterns of streamlines, vorticity, Reynolds stress correlations and turbulent kinetic energy distributions. Flow data reveal that the size of wake flow region, the location of singular and double points, the peak values of turbulence quantities, such as Reynolds stress correlations, vorticity fluctuations and turbulent kinetic energy vary as a function of models’ geometry and Reynolds Numbers. The concentration of small scale vortices is more dominant in the wake of the sphere than that of the cylinder. The maximum value of turbulent kinetic energy (TKE) occurs close to the saddle point for the cylinder case while two maximum values of TKE occur along shear layers for the sphere one because of the 3-D flow behavior.     

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.     

Momentum and heat transport in a finite-length cylinder wake

This paper reports an experimental study of turbulent momentum and heat transport in the wake of a wall-mounted finite-length square cylinder, with its length-to-width ratio L/d = 3–7. The cylinder was slightly heated so that heat produced could be considered as a passive scalar. A moveable three-wire probe (a combination of an X-wire and a cold wire) was used to measure velocity and temperature fluctuations at a Reynolds number of 7,300 based on d and the free-stream velocity. Measurements were performed at 10 and 20d downstream of the cylinder at various spanwise locations. Results indicate that L/d has a pronounced effect on Reynolds stresses, temperature variance and heat fluxes. The downwash flow from the free end of the cylinder acts to suppress spanwise vortices and, along with the upwash flow from the cylinder base, makes the finite-length cylinder wake highly three-dimensional. Reynolds stresses, especially the lateral normal stress, are significantly reduced as a result of suppressed spanwise vortices at a small L/d. The downwash flow acts to separate the two rows of spanwise vortices further apart from the wake centerline, resulting in a twin-peak distribution in temperature variance. While the downwash flow entrains high-speed fluid into the wake, responsible for a small deficit in the time-averaged streamwise velocity near the free end, it does not alter appreciably the distribution of time-averaged temperature. It has been found that the downwash flow gives rise to a counter-gradient transport of momentum about the central region of the wake near the free end of the cylinder, though such a counter-gradient transport does not occur for heat transport.     

减阻工况下壁面周期扰动对湍流边界层多尺度的影响

通过在平板壁面施加不同频率振幅的压电陶瓷振子周期性扰动,进行了湍流边界层主动控制减阻的实验研究.在压电陶瓷振子最大减阻工况下(80 V和160Hz),使用单丝边界层探针对压电振子自由端下游2mm处进行测量,得到不同法向位置流向速度信号的时间序列.通过对比施加控制前后的多尺度分析,发现压电振子产生的扰动只对近壁区产生影响,使得近壁区大尺度脉动降低,小尺度脉动强度增大,而对边界层的外区则基本没有影响.进一步对大尺度和小尺度的脉动信号进行条件平均,发现压电振子产生的扰动对小尺度脉动的影响在时间相位上并不均匀,小尺度脉动强度在大尺度脉动为正时比在大尺度脉动为负时具有更明显的增加.这表明壁面周期扰动主要通过使大尺度高速扫掠流体破碎为小尺度结构,来影响相应的高壁面摩擦事件,从而达到减阻效果.      

TR-PIV measurement of the wake behind a grooved cylinder at low Reynolds number

A comparative study of the wakes behind cylinders with grooved and smooth surfaces was performed with a view to understand the wake characteristics associated with the adult Saguaro cacti. A low-speed recirculation water channel was established for the experiment; the Reynolds number, based on the free-stream velocity and cylinder diameter (D), was kept at Re_D=1500. State-of-the-art time-resolved particle image velocimetry (TR-PIV) was employed to measure a total of 20 480 realizations of the wake field at a frame rate of 250 Hz, enabling a comprehensive view of the time- and phase-averaged wake pattern. In comparison to the wake behind the smooth cylinder, the length of the recirculation zone behind the grooved cylinder was extended by nearly 18.2%, yet the longitudinal velocity fluctuation intensity was considerably weakened. A global view of the peaked spectrum of the longitudinal velocity component revealed that the intermediate region for the grooved cylinder, which approximately corresponds to the transition region where the shear layer vortices interact, merge and shed before the formation of the Karman-like vortex street, was much wider than that for the smooth one. The unsteady events near St=0.3-0.4 were detected in the intermediate region behind the grooved cylinder, but no such events were found in the smooth cylinder system. Although the formation of the Karman-like vortex street was delayed by about 0.6D downstream for the grooved cylinder, no prominent difference in the vortex street region was found in the far wake for both cylinders. The Proper Orthogonal Decomposition (POD) method was used extensively to decompose the vector and swirling strength fields, which gave a close-up view of the vortices in the near wake. The first two POD modes of the swirling strength clarified the spatio-temporal characteristics of the shear layer vortices behind the grooved cylinder. The small-scale vortices superimposed on the shear layers behind the grooved cylinder were found to be generated and convected downstream in the same phase, which would significantly reduce the fluctuating force on the cylinder surface.     

Numerical simulation of the turbulent wake behind a normal flat plate

The three-dimensional wake flow behind a flat plate placed normal to the free stream has been investigated by means of direct numerical simulations. The Reynolds number Re based on the homogeneous inflow velocity and the uniform width d of the plate was 750. Coherent vortices were alternately shed from the sides of the plate with a frequency corresponding to a Strouhal number 0.168. The wake was distinctly turbulent downstream of the plate whereas the mean recirculation bubble extended 1.96d downstream. A steady 2D mean flow and the accompanying Reynolds stresses were obtained by averaging in time and along the span of the plate. These Reynolds-averaged statistics exhibited the same qualitative features as corresponding data from cylinder wakes.     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.     

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.     

Development of wake of a rectangular cylinder in a curved stream

In the wake of a rectangular cylinder measurements of mean velocity and some turbulent stresses are carried out in a straight duct and in a curved duct. The difference in turbulent quantities in the wake of the body, in the straight duct an in the curved duct is significant especially in the downstream side of the wake. The shear stresses are more sensitive to curvature than the normal stresses.     

Passive Scalar Transport in a Turbulent Cylinder Wake in the Presence of a Downstream Cylinder

This work aims to investigate how the presence of a downstream cylinder affects the passive scalar transport in a cylinder wake. The wake was generated by two tandem brass circular cylinders of the same diameter (d). The cylinder centre-to-centre spacing L/d was 1.3, 2.5 and 4.0, respectively, covering the three typical flow regimes of this flow. The upstream cylinder was slightly heated. Measurements were conducted at x/d= 10 and Re (≡ dU _∞/ν, where U _∞ is the free-stream velocity and ν is the kinematic viscosity of fluid) = 7000. A three-wire probe consisting of an X-wire and a cold wire was used to measure the velocity and temperature fluctuations, while an X-wire provided a phase reference. The phase-averaged velocity vectors and vorticity display single vortex street behind the downstream cylinder, irrespective of the flow regimes. However, the detailed flow structure exhibits strong dependence on L/d in terms of the Strouhal number, the vortex strength and its downstream evolution. This naturally affects passive scalar transport. The coherent and incoherent heat flux vectors show significant variation for different L/d.     

结构振动对湍流近尾迹的影响

研究了圆柱绕流中流体与结构的相互作用,侧重结构振动对湍流尾迹的影响,用激光测振仪测量圆柱在升力方向的位移;用热线和LDA（二维）测量湍流的近尾迹,通过变化自由流的速度和圆柱体直径（特征尺寸）来变化雷诺数,用两个振动特性不同的（一个相对刚性,一个相对弹）圆柱来产生尾迹,研究固体结构振动对湍流近尾迹的平均速度场和湍流场的影响,结果表明,结构自由振动对湍流近尾迹场影响明显,该影响随雷诺数的变化不明显。     

Automated topology classification method for instantaneous velocity fields

Topological concepts provide highly comprehensible representations of the main features of a flow with a limited number of elements. This paper presents an automated classification method of instantaneous velocity fields based on the analysis of their critical points distribution and feature flow fields. It uses the fact that topological changes of a velocity field are continuous in time to extract large scale periodic phenomena from insufficiently time-resolved datasets. This method is applied to two test-cases : an analytical flow field and PIV planes acquired downstream a wall-mounted cube.     

Separated flow structures around a cylindrical obstacle in a narrow channel

A detailed experimental study is performed on the separated flow structures around a low aspect-ratio circular cylinder (pin-fin) in a practical configuration of liquid cooling channel. Distinctive features of the present arrangement are the confinement of the cylinder at both ends, water flow at low Reynolds numbers (Re = 800, 1800, 2800), very high core flow turbulence and undeveloped boundary layers at the position of the obstacle. The horseshoe vortex system at the junctions between the cylinder and the confining walls and the near wake region behind the obstacle are deeply investigated by means of Particle Image Velocimetry (PIV). Upstream of the cylinder, the horseshoe vortex system turns out to be perturbed by vorticity bursts from the incoming boundary layers, leading to aperiodical vortex oscillations at Re = 800 or to break-away and secondary vorticity eruptions at the higher Reynolds numbers. The flow structures in the near wake show a complex three-dimensional behaviour associated with a peculiar mechanism of spanwise mass transport. High levels of free-stream turbulence trigger an early instabilization of the shear layers and strong Bloor–Gerrard vortices are observed even at Re = 800. Coalescence of these vortices and intense spanwise flow inhibit the alternate primary vortex shedding for time periods whose length and frequency increase as the Reynolds number is reduced. The inhibition of alternate vortex shedding for long time periods is finally related to the very large wake characteristic lengths and to the low velocity fluctuations observed especially at the lowest Reynolds number.