Comparison of flow patterns in the near wake of a circular cylinder and a square cylinder placed near a plane wall

The flow patterns in the near wake of a cylinder (either circular or square in shape, D=25 mm) placed in the proximity of a fully developed turbulent boundary layer (thickness δ=0.4D) are investigated experimentally using particle image velocimetry (PIV). The effects of changing the gap height (S) between the cylinder bottom and the wall surface, over the gap ratio range S/D=0.1–1.0, have been investigated. The results show that both the ensemble-averaged and instantaneous flow fields are strongly dependent on S/D. The flow patterns for the two types of cylinders share many similarities with respect to the change in S/D, such as the reduced recirculation length and increased velocity fluctuation in the near wake with increasing S/D, as well as the trend of suppression of vortex shedding at small S/D and onset of vortex shedding at large S/D. However, developments of the shear layers, in terms of wake width, flow curvature, etc., are considerably different for these two types of cylinders. In general, the wake development and momentum exchange for the square cylinder are slower those for the circular cylinder at the same gap ratio. Correspondingly, it is shown that the periodic vortex shedding is delayed and weakened in the case of square cylinder, as compared to that of the circular cylinder at the same S/D.     

Computational study on near wake interaction between undulation body and a D-section cylinder

We present a numerical study on the hydrodynamic performance of undulation NACA0012 foil in the near wake of D-section cylinder. Computations are conducted using unsteady incompressible Navier-Stokes equations with a moving adaptive mesh based on laminar flow. Investigations are focused on the effect of distance ratio between foil tip and centre of cylinder (L/D≤2.0) on the thrust/drag performance of foil and cylinder at various foil undulation frequency (St). We found that, foil thrust coefficient (C_t) increases considerably with the appearance of cylinder and an optimal distance exists at which C_t reaches maxima. The maximum increment is about eleven times that of its counterpart of single foil, which is obtained at St=0.23 and L/D=0.5. Our results for the cylinder drag coefficient (C_d) observed the existence of optimal parametric map, combined with various gap ratios and foil frequencies. With these parameters, insertion of an undulation foil can significantly lead to the drag reduction indicating that undulating foil could work efficiently as a passive vortex control device for cylinder drag reduction.     

Effect of peak perforation on flow past a conic cylinder at <Emphasis Type="Italic">Re</Emphasis>=100: drag,lift and Strouhal number

A flow past a circular-section cylinder with a perforated conic shroud, in which the perforation is located at the peak of the conic disturbance as the shroud installed on the cylinder and uniformly distributed with several circular holes, is numerically simulated at a Reynolds number of 100. Two factors in the perforation are taken into account, i.e. the attack angle relative to the direction of incoming flow and diameter of holes. The effect of such perforation on the drag, lift and vortex-shedding frequency is mainly investigated. Results have shown that variation of the attack angle has a little effect, especially on the drag and vortex-shedding frequency, except in certain cases due to the varied vortex-shedding patterns in the near wake. The increasing hole diameter still exhibits a little effect on the drag and frequency of vortex shedding, but really reduces the lift, in particular at larger wavelength, such as the lift reduction reaching almost 66%–68% after introducing the perforation.     

Effect of Peak Perforation on Flow Past A Conic Cylinder at Re=100: Drag, Lift and Strouhal Number

A flow past a circular-section cylinder with a perforated conic shroud, in which the perforation is located at the peak of the conic disturbance as the shroud installed on the cylinder and uniformly distributed with several circular holes, is numerically simulated at a Reynolds number of 100. Two factors in the perforation are taken into account, i.e. the attack angle relative to the direction of incoming flow and diameter of holes. The effect of such perforation on the drag, lift and vortex-shedding frequency is mainly investigated. Results have shown that variation of the attack angle has a little effect, especially on the drag and vortex-shedding frequency, except in certain cases due to the varied vortex-shedding patterns in the near wake. The increasing hole diameter still exhibits a little effect on the drag and frequency of vortex shedding, but really reduces the lift, in particular at larger wavelength, such as the lift reduction reaching almost 66%–68% after introducing the perforation.     

Three-dimensional flow around a square cylinder near a wall

In this study, two- and three-dimensional numerical simulations were performed to investigate the effect of the flow structure in the wake of a square cylinder placed near a plane wall by applying a fully implicit finite-difference method to the Navier-Stokes equations. The gap ratio between the cylinder and the wall, G/D, was varied from 0.2 to 4 for the Reynolds numbers of 175, 185 and 250. The role of the 3D structure on the lift and drag coefficients and Strouhal number was investigated. The results were compared with those of the 2D numerical simulations. The deviations of the 3D flow structure of the cylinder-wall pair from that of a single cylinder were also reported. At Re=250, B type secondary vortices were determined in the wake region. At Re=175 and 185, transition from A type vortex to fully periodic B type vortices was observed when the cylinder was brought closer to the wall.     

Numerical Simulation on Flow Past Two Side-by-Side Inclined Circular Cylinders at Low Reynolds Number

A series of three-dimensional numerical simulations is carried out to investigate the effect of inclined angle on flow behavior behind two side-by-side inclined cylinders at low Reynolds number Re=100 and small spacing ratio T/D=1.5 (T is the center-to-center distance between two side-by-side cylinders, D is the diameter of cylinder). The instantaneous and time-averaged flow fields, force coefficients and Strouhal numbers are analyzed. Special attention is focused on the axial flow characteristics with variation of the inclined angle. The results show that the inclined angle has a significant effect on the gap flow behaviors behind two inclined cylinders. The vortex shedding behind two cylinders is suppressed with the increase of the inclined angle as well as the flip-flop gap flow. Moreover, the mean drag coefficient, root-mean-square lift coefficient and Strouhal numbers decrease monotonously with the increase of the inclined angle, which follows the independent principle at small inclined angles.

     

不同雷诺数下倾斜圆柱绕流三维数值模拟研究

基于计算流体力学(CFD)开源代码OpenFOAM开展了不同雷诺数(Re=100、1500和3900)和倾斜角度(-60°≤α≤60°)工况下倾斜圆柱绕流流场的三维数值模拟,研究倾斜圆柱绕流的三维瞬时及时均尾流流场、流线拓扑、升阻力系数与旋涡脱落频率随雷诺数及倾斜角度变化的规律,探讨在顺流向及逆流向情况下独立性原则对倾斜圆柱绕流的适用性。研究结果表明:随着圆柱倾角的增大,倾斜圆柱尾流产生较为明显的轴向流,尾流旋涡脱落受到明显抑制,细碎旋涡逐渐消失,尾流宽度随之减小;随着雷诺数的增大,圆柱尾流涡管发生明显的变形,展向掺混使得大量细碎旋涡产生,呈现出明显的三维特性。在不同雷诺数下,阻力系数均值、升力系数均方根及无量纲涡脱频率在一定倾角范围内符合独立性原则。     

Numerical evaluation of passive control of VIV by small control rods

Flow past a circular cylinder with multiple small control rods is studied by numerical simulation for Re_D ranging from 1161.3 to 6387.1. The Reynolds-Averaged-Navier–Stokes (RANS) equations and shear stress transport (SST) k − ω turbulence model are used to calculate the vortex field, while a fourth-order Runge–Kutta method is employed for evaluating the structure dynamics of the cylinder group. Comparisons with experimental results demonstrate the validation of this method. This study is concerned with the vortex induced vibration (VIV) suppression efficacy of small control rods placed around a main cylinder. The effects of control rod number, diameter ratio, spacing ratio and Reynolds number on the hydrodynamics and vibration responses of the main cylinder are investigated. The reduced percents of in-line and cross-flow amplitudes and the increased percents of the whole cross-sectional area of cylinders and the drag coefficient are used to give a comprehensive evaluation. Results of simulation indicate that placing small rods with appropriate number at appropriate locations can achieve good suppression effectiveness at a wide range of Reynolds number. The numerical result for the case with nine control rods, diameter ratio of 0.15 and spacing ratio of 0.6 shows the best suppression effect among the cases investigated in this study.     

Flow past a circular cylinder between parallel walls at low Reynolds numbers

The flow about a circular cylinder placed centrally inside a channel is studied numerically with an unstructured collocated grid finite volume method based on the primitive variable formulation. The distance between the channel walls is allowed to vary to change the blockage ratio. Simulations are carried out over a range of Reynolds numbers that are consistent with the two-dimensional assumption. The study confirms that transition to vortex shedding regime is delayed when the channel walls are close to the cylinder because of the interaction between the vortices from the channel wall and cylinder wake. In the unsteady vortex shedding regime, the wake pattern is opposite to the classic Karman street in respect of the positions of the shed vortices. The cylinder drag coefficient and Strouhal number are considerably increased at smaller gaps while the root-mean-squared lift coefficient is significantly decreased. Several important flow parameters are correlated with the input parameters, namely Reynolds number and blockage ratio.     

五圆柱体结构群绕流特性及互扰效应研究

基于Fluent流体计算平台,运用大涡模拟方法对亚临界雷诺数Re=3900下“X”形排列五圆柱体结构群三维绕流特性进行研究,主要分析来流攻角α与间距比L/D两个关键参数对五圆柱体结构群的尾流区三维涡结构演化与流体力系数的影响,并揭示其内在流动互扰机理。研究表明:来流攻角和间距比的变化对五圆柱体结构群流动控制及互扰效应的影响显著。在小间距比工况下,观察到柱体群间隙区域内流体高速流动的现象,导致五圆柱体之间的互扰作用十分强烈。间隙流对中间圆柱体和下游圆柱体有较强的冲击作用,对其表面的流体力分布特性有显著的影响。另外,大间距比工况下,当α=0°与L/D≥5.0工况时,柱体群尾流效应强于其间隙流效应。当α=22.5°与L/D=7.0时,位于下游与中间处的圆柱体流体绕流特性存在较大差异。而当α=45°与L/D≥6.0时,位于上游与中间处的圆柱体尾流区均会产生正负交替的漩涡结构。     

Direct numerical simulation of flow over a slotted cylinder at low Reynolds number

Direct numerical simulation was conducted to investigate the flow past a slotted cylinder at low Reynolds number (Re) of 100. The slotting of cylinder affects the boundary layer separation, vortex formation position, recirculation region length and wake width, which are determined by the type of slit. The streamwise slit (SS1), T-shaped slit (SS3) and Y-shaped slit (SS4) act as passive jets, while the transverse slit (SS2) achieves an alternate self-organized boundary layer suction and blowing. The flow rate in slits fluctuates over time due to the alternate vortex shedding and fluctuating pressure distribution around the cylinder surface. One fluctuation cycle of flow rate is caused by a pair of vortices shedding for SS2, SS3 and SS4, while it is created by each vortex shedding for SS1. The wall shear stress and flow impact on the slit wall partly contribute to the hydrodynamic forces acting on the slotted cylinder. Taking into account the internal wall of slit, the transverse slit plays the best role in suppressing the fluid forces with drag reduction of 1.7% and lift reduction of 17%.     

Experiments with flexible shrouds to reduce the vortex-induced vibration of a cylinder with low mass and damping

Experiments employing a low-mass-damping cylinder have been conducted to determine the vortex-induced vibration (VIV) response of four suppressors of the flexible-shroud family. The VIV suppressors were inspired in the concept of the Ventilated Trousers (VT), a flexible shroud composed of a flexible net fitted with three-dimensional bobbins. Reynolds number varied between 5 × 10³ and 25 × 10³, while reduced velocity varied from 2 to 26. The VIV dynamic response showed that the VT suppressed the peak amplitude of vibration down to 40% of that of a bare cylinder. Other flexible shrouds also achieved suppression, but not as efficiently. Drag was reduced during the VIV synchronization range, but remained above the value for a bare static cylinder thereafter. Spectral analysis of displacement and lift revealed that, depending on the geometry and distribution of the bobbins, the flexible shroud can develop an unstable behavior, capturing energy from the wake and sustaining vibrations for higher reduced velocities. PIV measurements of the wake revealed that the entrainment flow through the mesh is necessary to extend the vortex-formation length of the wake; this mechanism only occurs for the VT mesh.     

Analysis of wake behind a rotating propeller using PIV technique in a cavitation tunnel

A two-frame particle image velocimetry (PIV) technique is used to investigate the wake characteristics behind a marine propeller with 4 blades at high Reynolds number. For each of 9 different blade phases from 0° to 80°, 150 instantaneous velocity fields are measured. They are ensemble averaged to study the spatial evolution of the propeller wake in the region ranging from the trailing edge to one propeller diameter (D) downstream location. The phase-averaged mean velocity shows that the trailing vorticity is related to radial velocity jump, and the viscous wake is affected by boundary layers developed on the blade surfaces and centrifugal force. Both Galilean decomposition method and vortex identification method using swirling strength calculation are very useful for the study of vortex behaviors in the propeller wake region. The slipstream contraction occurs in the near-wake region up to about X/D=0.53 downstream. Thereafter, unstable oscillation occurs because of the reduction of interaction between the tip vortex and the wake sheet behind the maximum contraction point.     

A finite difference solution of the shear flow over a circular cylinder

The incompressible viscous shear flow past a circular cylinder is analyzed by solving two-dimensional Navier-Stokes equations and pressure Poisson equation using a finite difference method. The shear flow is calculated for Reynolds numbers from 80 to 1000, and shear parameters up to 0.25. The numerical results indicate that the vortex shedding persists at the shear parameters up to 0.25 for the present Reynolds number range. The Strouhal number and the drag coefficient decrease as the shear parameter increases. There is a transverse force acting from the high velocity side toward the low velocity side in the shear flow.     

The virtual mass of a rectangular flat plate of finite aspect ratio

The skin friction of a two-dimensional planing flat plate is made up of two opposing components; a drag force from the flow aft of the stagnation line and an opposing thrust force from the jet flow. This paper is concerned only with the drag term and the wake velocity defect which it causes in the water behind the transom.It is concluded that the skin friction is less than would be expected from a flat plate at ambient static pressure (D_fo say) and is approximately equal to D_fo (1 − C_R), where C_R is the normal force coefficient based on wetted area. The wake velocity decrement due to this drag is found to be significant, particularly for surface piercing propellers.     

动波浪壁圆柱绕流三维数值模拟

利用计算流体力学软件Fluent开展了三维动波浪壁圆柱绕流的数值模拟,建立了三维运动波浪壁圆柱模型,通过C语言自编程序实现波浪壁面的运动控制,并保证壁面变形时网格的高质量。在来流速度u=0.125 m/s、雷诺数Re=12 500的情况下,开展了动波浪壁波动速度w=0、0.062 5、0.125、0.187 5 m/s四个工况的计算分析,并比较了不同波动速度对流场结构、升力、阻力特性的影响。结果表明:动波浪壁圆柱能有效抑制流动的分离,消除交替脱落的尾涡,从而消除周期振荡的升力;在消除卡门涡街的同时,圆柱后驻点处的涡量值随波动速度增加而增加,其原因在于波形移动加大了壁面流体的速度,从而减小了圆柱前后的压力差,减小了阻力;随着波动速度的增大,平均阻力系数呈明显下降趋势,当波动速度为来流速度的1.5倍时,平均阻力系数相对于光滑圆柱下降了53.76%。     

On the skin friction and resulting wake of a two-dimensional planing plate

The skin friction of a two-dimensional planing flat plate is made up of two opposing components; a drag force from the flow aft of the stagnation line and an opposing thrust force from the jet flow. This paper is concerned only with the drag term and the wake velocity defect which it causes in the water behind the transom.It is concluded that the skin friction is less than would be expected from a flat plate at ambient static pressure (D_fo say) and is approximately equal to D_fo (1 ? C_R), where C_R is the normal force coefficient based on wetted area. The wake velocity decrement due to this drag is found to be significant, particularly for surface piercing propellers.     

亚临界区雷诺数下圆柱绕流场电磁力控制数值研究

在亚临界区雷诺数下,采用脱体涡模拟方法对弱电解质中电磁力作用下圆柱绕流场及其升阻力特性进行了数值模拟与分析。结果表明,电磁力可以提高圆柱体边界层内的流体动能,延缓圆柱体近壁面流动分离,减弱绕流场中流向和展向大尺度涡的强度,减小圆柱体阻力及其升力脉动幅值;当电磁力作用参数大于某个临界值后,流动分离角消失,在圆柱体尾部产生射流现象,电磁力产生净推力作用,出现负阻力现象,而且升力脉动幅值显著减小且接近于零。     

不等直径串列双圆柱体绕流的数值模拟

利用Fluent软件模拟雷诺数为200时不同间距比G/D和直径比d/D情况下的圆柱绕流现象。根据模拟结果分析G/D和d/D对圆柱体的涡脱落形态,大、小柱体的升力,阻力系数和St数的影响,结果表明涡脱落形态随着G/D和d/D的变化呈现不同的形式,在G/D小于临界间距比时呈现单一涡脱落形态,在G/D大于临界间距比时,呈现双旋涡脱落形态;临界间距比随着d/D增大而增大。在临界间距比附近大、小柱体的升力和阻力系数值及St数变化较大,大、小柱体的St数具有相同的变化规律,St随着d/D的增大而减小。     

两层粘性流体中圆柱体受迫振荡数值模拟

研究两层粘性流体中无限长水平圆柱体的受迫振荡问题。在湍流模式下,采用VOF方法追踪两层流体的内界面,基于动网格技术模拟圆柱体的运动边界,对均匀流中横向振荡圆柱体的绕流场进行了数值模拟。计算受迫振荡圆柱体的升力系数、阻力系数随时间的演化曲线和圆柱体的尾涡分布,以及圆柱体的受迫振荡激发两层流体内界面的扰动,并与均匀流体的情况进行了比较分析。研究表明,流体的两层分层效应对受迫振荡圆柱体的升阻力系数和尾涡分布特性都有显著影响,在水下输油气管道涡激振动特性的工程评估中,应考虑流体的密度分层效应。