Large eddy simulation of cross-flow through a staggered tube bundle at subcritical Reynolds number

The Large Eddy Simulation (LES) technique is used to study the vortex shedding characteristics inside a staggered tube array consisting of six rows with intermediate spacings (S_L/D×S_T/D=1.6×3.6) at the subcritical Reynolds number of 8600 (based on the gap velocity). The filtered equations are discretised using the finite volume method in an unstructured, collocated grid arrangement with second-order accurate methods in space and time. The predictions of mean velocities and Reynolds stresses are in very good agreement with detailed LDA measurements performed in 17 stations along the depth of the array. The sizes of the recirculation zones behind the cylinders in the first and third row also compare favourably with available correlations. Two distinct and independent shedding frequencies are detected behind the first two rows, but the high-frequency component vanishes in the downstream rows. The corresponding Strouhal numbers agree well with measurements available in the literature for similar tube spacings. The lift coefficient as well as instantaneous flow patterns of the whole array are also examined.     

Flow-induced vibration of a circular cylinder subjected to wake interference at low Reynolds number

Two- and three-dimensional numerical simulations of the flow around two circular cylinders in tandem arrangements are performed. The upstream cylinder is fixed and the downstream cylinder is free to oscillate in the transverse direction, in response to the fluid loads. The Reynolds number is kept constant at 150 for the two-dimensional simulations and at 300 for the three-dimensional simulations, and the reduced velocity is varied by changing the structural stiffness. The in-line centre-to-centre distance is varied from 1.5 to 8.0 diameters, and the results are compared to that of a single isolated flexible cylinder with the same structural characteristics, m^?=2.0 and ζ=0.007. The calculations show that significant changes occur in the dynamic behaviour of the cylinders, when comparing the flow around the tandem arrangements to that around an isolated cylinder: for the tandem arrangements, the lock-in boundaries are wider, the maximum displacement amplitudes are greater and the amplitudes of vibration for high reduced velocities, outside the lock-in, are very significant. The main responsible for these changes appears to be the oscillatory flow in the gap between the cylinders.     

Numerical predictions of low Reynolds number flows over two tandem circular cylinders

Flows over two tandem cylinders were analysed using the newly developed collocated unstructured computational fluid dynamics (CUCFD) code, which is capable of handling complex geometries. A Reynolds number of 100, based on cylinder diameter, was used to ensure that the flow remained laminar. The validity of the code was tested through comparisons with benchmark solutions for flow in a lid‐friven cavity and flow around a single cylinder. For the tandem cylinder flow, also mesh convergence was demonstrated, to within a couple of percent for the RMS lift coefficient. The mean and fluctuating lift and drag coefficients were recorded for centre‐to‐centre cylinder spacings between 2 and 10 diameters. A critical cylinder spacing was found between 3.75 and 4 diameters. The fluctuating forces jumped appreciably at the critical spacing. It was found that there exists only one reattachment and one separation point on the downstream cylinder for spacings greater than the critical spacing. The mean and the fluctuating surface pressure distributions were compared as a function of the cylinder spacing. The mean and the fluctuating pressures were significantly different between the upstream and the downstream cylinders. These pressures also differed with the cylinder spacing. Copyright © 2004 John Wiley & Sons, Ltd.     

A three-dimensional simulation of a steady approach flow past a circular cylinder at low Reynolds number

The three-dimensional (3D) unsteady viscous wake of a circular cylinder exposed to a steady approach flow is calculated using a fractional-step finite-difference/spectral-element method. The calculated flow fields at Reynolds numbers of 100 (2D) and 200 (3D) are examined in detail. The flow field at Re = 100 is 2D as expected, while the flow field at Re = 200 has distinct 3D features, with spanwise wavelengths of about 3.75 cylinder diameters. The calculated results produce drag and lift coefficients and Strouhal numbers that agree extremely well with the experimental values. These 3D values at Re = 200 are in better agreement with experimental values than the results of a 2D calculation at Re = 200, which is expected. © 1998 John Wiley & Sons, Ltd.     

圆柱绕流流场结构的大涡模拟研究

为进一步揭示绕流现象的形成机理,本文分别对处于层流稳态区、尾流过渡区、剪切层转换区Re分别为26、200、1.4×105的三种典型流态下的单圆柱绕流进行了二维数值模拟研究。Re为26时应用层流模型直接求解N-S方程,而Re分别为200、1.4×105时使用大涡模拟的方法进行计算。数值模拟很好地再现了稳定的涡旋结构、周期性交替脱落的卡门涡街结构、不规则的涡旋结构,在此基础上分析了尾流结构的基本特征及其压强分布规律、平均的流场特性、积分参数(如升力系数、阻力系数、斯特劳哈尔等),并与有关研究成果进行了对比。研究发现,采用不同流动介质时流场特性有所差异,空气为介质时的计算结果更符合实验的成果,而水为介质时计算结果偏差较大,这主要是由尾流涡旋产生的不合理负压造成的。     

较高Re数圆柱尾流的控制

引入一个窄条作为控制件,在Re=3.0×10 3~2.0×10 4范围内对圆柱尾流进行控制实验。窄条长度与柱体长度相同,厚 度为柱体直径的 0.015~0.025倍,宽度为柱体直径的0.18倍. 窄条的两个长边 与柱中心轴平行, 而且三者共面. 控制参数为窄条位置, 可由间距(窄条到柱轴)比λ/(0.5D)和风向角β (窄 条面与来流的夹角)确定. 采用流动显示和热线测量方法,对控制和未控制尾流的流动状态, 平均速度分布和脉动速度情况,以及作用于柱体和控制件的总阻力进行了研究和比较. 研究结果证明, 当窄条位于柱体尾流中一定区域内时, 可有效抑制柱体两侧的旋涡脱落.有效控制后的尾流湍流度也相应减小. 在不同Re数下,找出了有效抑制旋涡脱落的窄条位置区域, 并用动量积分估计了作用于柱体和窄条上的总阻力与光圆柱阻力的比值及其随风向角的变 化. 对λ/(0.5D)=2.9情况,得到了减阻的风向角区域(β=0°~40°与180°附近)以及最大减阻率32%.以上事实表明,在近尾流局部区域施加小的干扰,可改变较高Re数圆柱尾流的整体性质.     

Vortex shedding and aerodynamic forces on a circular cylinder in linear shear flow at subcritical Reynolds number

Vortex shedding and aerodynamic forces on a circular cylinder in a linear shear flow with its axis normal to the plane of the velocity shear profile at subcritical Reynolds number are investigated experimentally. The shear parameter β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.27. The Strouhal number has no significant variation with the shear parameter. The time-mean base pressure increases and the fluctuating component of the base pressure decreases significantly with increasing shear parameter. Vortex shedding is suppressed by the velocity shear. Dislocation of the stagnation point takes place and this influences the pressure distribution around the cylinder together with the velocity shear. A mean lift force arises in the shear flow due to asymmetry of the pressure distribution, and it acts from the high velocity side to the low velocity side. In addition, the lift coefficient increases and the drag coefficient decreases with increasing shear parameter.     

Numerical investigation on vortex-induced vibration of an elastically mounted circular cylinder at low Reynolds number using the fictitious domain method

A direct numerical simulation of two-dimensional (2D) flow past an elastically mounted circular cylinder at low Reynolds number using the fictitious domain method had been undertaken. The cylinder motion was modelled by a two degree-of-freedom mass–spring–damper system. The computing code was verified against a benchmark problem in which flow past a stationary circular cylinder is simulated. Then, analyses of vortex-induced vibration (VIV) responses, drag and lift forces and the phase and vortex structures were carried out. Results show that the cylinder's non-dimensional cross-flow response amplitude reaches its summit of 0.572 in the ‘lock-in’ regime. The ‘2S’, instead of the ‘2P’, vortex shedding mode is dominated in the ‘lower’ branch for this 2D low-Re VIV. A secondary oscillation is observed in the lift force when ‘lock-in’ occurs. It is shown that this secondary component changes the phase, offset the energy input by the primary component and thus reduces the cylinder responses. Effects of the Skop–Griffin parameter on cylinder responses were also investigated.     

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.     

An experimental study of the dynamic characteristics in boundary layer in the flow past a two-dimensional circular cylinder at subcritical reynolds numbers

An experimental study of the dynamic characteristics of flow past a two-dimensional circular cylinder is described. The fluctuationsoof wall shear stress, surface-pressure and velocity of the flow are measured with hot-film, hot-wire and pressure transducer. The frequency feature of fluctuations of wall shear stress is given. The cross-correlation functions of these fluctuations at any two points are calculated. The experimental results reveal that there is an overall syncronous fluctuation, at the shedding frequency, in boundary layer in the flow past a two-dimensional circular cylinder at subcritical Reynolds number.     

Numerical investigation of fluid flow past circular cylinder with multiple control rods at low Reynolds number

Laminar flow past a circular cylinder with multiple small-diameter control rods is numerically investigated in this study. The effects of rod-to-cylinder spacing ratio, rod and cylinder diameter ratio, cylinder Reynolds number, number of control rods and angle of attack on the hydrodynamics of the main circular cylinder are investigated. Four different flow regimes are identified based on the mechanism of lift and drag reduction. The range of rod-to-cylinder spacing ratio where significant force suppression can be achieved is found to become narrower as the Reynolds number increases in the laminar regime, but is insensitive to the diameter ratio. The numerical results for the case with six identical small control rods at Re=200 show that the lift fluctuation on the main cylinder can be suppressed significantly for a large range of spacing ratio and various diameter ratios, while the drag reduction on the main cylinder is also achieved simultaneously. The six-control-rod arrangement has shown better performance in flow control than the arrangements with less control rods, especially in terms of force reduction at various angles of attack.     

Large-eddy simulation of turbulent flow and heat transfer in plane and rib-roughened channels

Large-eddy simulation results are presented and discussed for turbulent flow and heat transfer in a plane channel with and without transverse square ribs on one of the walls. They were obtained with the finite-difference code Harwell-FLOW3D, Release 2, by using the PISOC pressure-velocity coupling algorithm, central differencing in space, and Crank-Nicolson time stepping. A simple Smagorinsky model, with van Driest damping near the walls, was implemented to model subgrid scale effects. Periodic boundary conditions were imposed in the streamwise and spanwise directions. The Reynolds number based on hydraulic diameter (twice the channel height) ranged from 10 000 to 40 000. Results are compared with experimental data, k-? predictions, and previous large-eddy simulations.     

Adaptive mesh refinement method-based large eddy simulation for the flow over circular cylinder at ReD = 3900

With the development of computational power, large eddy simulation (LES) method is increasingly used in simulating complex flow. However, there still exist many factors affecting the LES quality and appropriate mesh resolution is among one of them. This work aims to develop an automatic procedure to refine the LES mesh by combining adaptive mesh refinement (AMR) and LES quality criteria. An LES refinement criterion is developed by estimating the proper grid length scale which meets the accuracy requirement of LES method. With this criterion, the baseline mesh is automatically refined with the AMR method. In this work, an efficient one-shot refinement strategy is also proposed to reduce the overall simulation time. Current AMR-based LES method is verified with the typical LES test case about the flow past circular cylinder at Re _D = 3900. Results show that the automatically refined mesh provides systematically better agreement with experimental results and with current method the balance between accuracy and computational expense for LES can be obtained.     

Numerical simulation of an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder

This article presents a numerical investigation of turbulent flow in an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder. The governing equations were discretized by the finite-volume method and SIMPLER method was applied to solve the equations on a staggered grid. The turbulent flow was numerically simulated using the standard k–ε, Abe–Kondoh–Nagano (AKN) and Shear Stress Transport (SST) turbulence models. The comparisons made between numerical results and experimental measurements showed that the SST model is superior to other models in the present calculation.Computations were performed for three different Reynolds numbers of 6000, 10 000 and 20 000 based on the cylinder diameter. To our knowledge, this study represents the first numerical investigation of the present flow configuration. The computational results were validated with the available experimental data of reattachment length, mean velocity distribution and wall static pressure coefficient in the turbulent blunt circular cylinder flows. Further, other characteristics of the flow, such as turbulent kinetic energy, pressure, streamlines, and the velocity vectors are discussed.The results show that the main characteristics of the turbulence flow in the separation region, such as reattachment length or velocity profiles, are nearly independent of the Reynolds number. The obtained results showed that a secondary separation bubble may appear in the main separation bubble near the leading edge. Furthermore, it was found that the turbulent kinetic energy has a large effect on the formation of the secondary bubble.     

Numerical study of the effects of particles on the near wake around a circular cylinder

In this work, we investigate the dynamics of the near wake in a turbulent flow going past a circular cylinder with/without particles at a moderate Reynolds number using a direct numerical simulation method. High-order finite-deference schemes are applied to solve for the bulk fluid properties, and a Lagrangian approach is adopted to track the individual particles. The single-phase flow is analysed and validated using previous experimental data. Two converged states, U- and V-shaped, are observed in the near wake, which are consistent with the experimental results. For the two-phase flow, the addition of smaller particles shortens the length of the recirculation region and causes a V-shaped profile to form behind the circular cylinder. Furthermore, the particles increase the drag force from the circular cylinder and suppress the vortex shedding frequency. An increase in the turbulent statistics in the very near wake and a decrease in the turbulent statistics further downstream are also observed.