Large-eddy simulation of the flow past a circular cylinder at sub- to super-critical Reynolds numbers

Large-eddy simulation of turbulent flow past a circular cylinder at sub- to super-critical Reynolds numbers is performed using a high-fidelity orthogonal curvilinear grid solver. Verification studies investigate the effects of grid resolution, aspect ratio and convection scheme. Monotonic convergence is achieved in grid convergence studies. Validation studies use all available experimental benchmark data. Although the grids are relatively large and fine enough for sufficiently resolved turbulence near the cylinder, the grid uncertainties are large indicating the need for even finer grids. Large aspect ratio is required for sub-critical Reynolds number cases, whereas small aspect ratio is sufficient for critical and super-critical Reynolds number cases. All the experimental trends were predicted with reasonable accuracy, in consideration the large facility bias, age of most of the data, and differences between experimental and computational setup in particular free stream turbulence and roughness. The largest errors were for under prediction of turbulence separation.     

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

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

Large eddy simulations of a circular cylinder at Reynolds numbers surrounding the drag crisis

Large eddy simulations of the flow around a circular cylinder at high Reynolds numbers are reported. Five Reynolds numbers were chosen, such that the drag crisis was captured. A total of 18 cases were computed to investigate the effect of gridding strategy, turbulence modelling, numerical schemes and domain width on the results. It was found that unstructured grids provide better resolution of key flow features, when a ‘reasonable’ grid size is to be maintained.When using coarse grids for large eddy simulation, the effect of turbulence models and numerical schemes becomes more pronounced. The dynamic mixed Smagorinsky model was found to be superior to the Smagorinsky model, since the model coefficient is allowed to dynamically adjust based on the local flow and grid size. A blended upwind-central convection scheme was also found to provide the best accuracy, since a fully central scheme exhibits artificial wiggles, due to dispersion errors, which pollute the solution.Mean drag, fluctuating lift Strouhal number and base pressure are compared to experiments and empirical estimates for Reynolds numbers ranging from 6.31 × 10⁴ to 5.06 × 10⁵. In terms of the drag coefficient, the drag crisis is well captured by the present simulations, although the other integral quantities (rms lift and Strouhal number) show larger discrepancies. For the lowest Reynolds number, the drag is seen to be more sensitive to the domain width than the spanwise grid spacing, while at the higher Reynolds numbers the grid resolution plays a more important role, due to the larger extent of the turbulent boundary layer.     

低雷诺数下圆柱强迫振动特性光滑粒子流体动力学模拟

基于光滑粒子流体动力学(SPH)方法,开发了能够准确描述水流作用下圆柱强迫振动特性的数学模型。通过引入适合于无网格粒子法的开边界算法,来模拟出入流边界条件,建立了具有造流功能的SPH数值水槽。圆柱及计算域的上下边界均采用修正的动力边界条件进行模拟。借助于粒子位移矫正和压力修正算法,避免了圆柱周围流体粒子压力大幅震荡以及结构下游区域出现空腔等非物理性现象。使用典型的圆柱绕流数据,验证了所建SPH模型的计算性能,研究了固定圆柱在低雷诺数情况下的尾涡脱落模式和升阻力变化规律。明确了低雷诺数下强迫振动圆柱在频率锁定以及非锁定区间内的升力变化规律,量化了升力与外界激励频率之间的关系。     

Force and flow characteristics of a circular cylinder with uniform surface roughness at subcritical Reynolds numbers

The main purpose of this study is to establish a better understanding of the relationship between drag reduction and surface roughness. Experiments were conducted to measure the force and flow characteristics of a circular cylinder with different types of artificial surface roughness over the range 6 × 10³ < Re < 8 × 10⁴ (Re is based on the cylinder diameter D). The roughness cylinder was formed by covering the exterior surface of the cylinder with uniformly distributed (1) sandpaper, (2) netting, and (3) dimples. The roughness coefficient ranged from k/D = 0.0028 to 0.025 (k is the roughness height). A detailed quantitative measurement of the flow field around the cylinder using Particle Imaging Velocimetry (PIV) was carried out. The hydrodynamic force coefficients (drag and lift) of the rough cylinders are compared against those of a smooth cylinder measured under the same flow conditions. It is found that certain configuration of surface roughness significantly reduces the mean drag coefficient of the cylinder, particularly at large Reynolds numbers. In addition, the root-mean-square (r.m.s.) lift coefficient of the rough cylinders is considerably lower than that of a smooth cylinder.     

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

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

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%.     

The separation of flow past a circular cylinder on aβ-plane

Effects of the Ekman friction on the prograde (eastward) flows past a cylinder on a-plane are investigated when (=R ²/U, whereR is the cylinder radius andU the freestream speed)O(1) and(=2E _k ^1/2/R ₀·O(1) where is the non-dimensional beta parameter and the ratio of the square root of the Ekman numberE _k multiplied by 2 to the Rossby numberRo multiplied by the aspect ratio(=H/R, whereH is the fluid depth). Previous studies without the Ekman friction have shown that the-effect inhibits flow separation for pragrade flows through the asymptotic boundary condition by shifting the region of the adverse pressure gradient toward the rear stagnation point. It is found that the Ekman friction alleviates this-effect on the exterior flow. In the E_k ^1/4-boundary layer, on the other hand, Ekman friction suppresses the vorticity advection along the wall, which tends to make the boundary layer thickness thin and delay the flow separation. The Ekman friction thus affects flow separation in a complicated manner. Details of the boundary layer structures and the separation angles are described for 0.3< <4.0 and 0.1<<1.5.     

均匀流中直立圆柱体绕流三维数值模拟

研究直立贯底圆柱体的三维粘性绕流问题。以不可压缩Navier-Stokes方程为控制方程,采用有限体积法和SIMPLE算法,建立了数值模拟方法。考察在不同水平和垂直断面上,圆柱体绕流产生的尾涡和流动速度场的分布特性,成功地数值模拟了直立贯底圆柱体绕流场的三维特性。结果表明,在考虑重力影响的情况下,直立圆柱体周围的流动具有明显的三维特性,而且沿圆柱体轴向不同断面上的尾涡分布是不相同的。     

Flow induced vibrations of a sharp edge square cylinder in the wake of a circular cylinder

The response of an oscillating circular cylinder at the wake of an upstream fixed circular cylinder was classified by different researchers as galloping, wake induced galloping or wake induced vibration. Furthermore it is already known that a sharp edge square cylinder would undergo galloping if it is subjected to uniform flow. In this study the influence of the wake of a fixed circular cylinder on the response of a downstream square cylinder at different spacing ratios (S/D = 4, 8, 11) is experimentally investigated. The subject appears not to have received previous attention. The lateral displacements, lift forces and the pressure data from gauges mounted in the wake of the oscillating cylinder are recorded and analyzed. The single degree of freedom vibrating system has a low mass-damping parameter and the Reynolds number ranges from 7.7 × 10² to 3.7 × 10⁴.In contrast to that for two circular cylinders in tandem arrangement, the freely mounted downstream square cylinder displays a VIV type of response at all spacing ratios tested. There is no sign of galloping or wake induced galloping with the square cylinder. With increase at the spacing ratio the cross-flow oscillations decrease. It is shown that the vortices arriving from the upstream fixed circular cylinder play a major role on the shedding mechanism behind the downstream square cylinder and cause the square cylinder to shed vortices with frequencies above Strouhal frequency of the fixed square cylinder (St = 0.13). The VIV type of oscillations in the downstream square cylinder is most probably caused by the vortices newly generated behind the square cylinder.     

Large-eddy simulation of the turbulent near wake behind a circular cylinder: Reynolds number effect

The purpose of the present work is <!–<query>The highlights are in an incorrect format. Hence they have been deleted. Please refer the online instructions: http://www.elsevier.com/highlights and provide 3-5 bullet points.</query>–>to study the effect of the Reynolds number on the near-wake structure and separating shear layers behind a circular cylinder. Three-dimensional unsteady large-eddy simulation is carried out and two different subgrid scale models are applied in order to evaluate the turbulent wake reasonably. The Reynolds number based on the free-stream velocity and the cylinder diameter is ranging from Re = 5500–41,300, corresponding to the full development of the shear-layer instability in the intermediate subcritical flow regime. For a complete validation of this numerical study, hydrodynamic bulk coefficients are computed and compared to experimental measurements and numerical studies in the literature. Special focus is made on the variations of both the large-scale near-wake structure and the small-scale shear-layer instability with increasing Reynolds numbers. The present numerical study clearly shows the broadband nature of the shear-layer instability as well as the dependence of the shear-layer frequency especially on the high Reynolds numbers.     

Benchmark computations for flows around a stationary cylinder with high Reynolds numbers by RANS-overset grid approach

The aim of this paper is to evaluate the accuracy, stability and efficiency of the overset grid approach coupled with the RANS (Reynolds Averaged Navier-Stokes) model via the benchmark computations of flows around a stationary smooth circular cylinder. Two dimensional numerical results are presented within a wide range of Reynolds numbers (6.31 × 10⁴ ∼ 7.57 × 10⁵) including the critical flow regime. All the simulations are carried out using the RANS solver pimpleFoam provided by OpenFOAM, an open source CFD (Computational Fluid Dynamics) toolkit. Firstly, a grid convergence study is performed. The results of the time-averaged drag and lift force coefficients, root-mean square value of lift force coefficient and Strouhal number (St number) are then compared with the experimental data. The velocity, vorticity fields and pressure distribution are also given. One main conclusion is that the numerical solutions in regard to a fixed cylinderare not deteriorated due to the implementation of the overset grid. Furthermore, it can be an appealing approach to facilitate simulations of Vortex Induced Vibrations (VIV), which involves grid deformation. The present study is a good start to implement the overset grid to solve VIV problems in the future.     

横向受迫振荡圆柱低雷诺数绕流问题数值模拟

对雷诺数Re=200的圆柱在均匀来流中横向受迫振荡的问题进行数值模拟分析.通过改变运动圆柱的振荡频率和振幅,分析圆柱所受作用力及圆柱尾流泻涡结构的变化特性,振荡圆柱在非锁定状态和锁定状态下作用力曲线和圆柱尾流泻涡结构的变化情况,以及不同圆柱振荡频率和振幅情况下单个振荡周期内圆柱尾流泻涡结构模式之间的转换特征.     

Drag force on a circular cylinder in unseparated laminar high Reynolds number nonlinear random oscillatory flow

The stochastic properties of the drag force maxima on a circular cylinder subjected to nonlinear random waves are investigated. Unseparated laminar high Reynolds number flow is considered. A simplified approach based on second order Stokes waves is presented, including the sum-frequency effect only. It is demonstrated how a drag force formula valid for regular linear waves can be used to find the cumulative distribution function of individual drag force maxima for nonlinear irregular waves. Here the [Wang, 1968] drag force coefficient is used.     

Hydrodynamic forces on a partially submerged cylinder at high Reynolds number in a steady flow

The hydrodynamic forces on the stationary partially submerged cylinder are investigated through towing test with Reynolds number ranging from 5 × 10⁴ to 9 × 10⁵. Three test groups of partially submerged cylinders with submerged depths of 0.25 D, 0.50 D, and 0.75 D and one validation group of fully submerged cylinders are conducted. During the experiments, the hydrodynamic forces on the cylinders are measured using force sensors. The test results show a considerable difference in the hydrodynamic coefficients for the partially submerged cylinders versus the fully submerged cylinders. A significant mean downward lift force is first observed for the partially submerged cylinders in a steady flow. The maximum of the mean lift coefficients can reach 1.5. Two distinct features are observed due to the effects of overtopping: random distributions in the mean drag coefficients and a clear quadratic relationship between the mean lift coefficients and the Froude number appear in the non-overtopping region. However, the novel phenomenon of a good linear relationship with the Froude number for the mean hydrodynamic coefficients is clearly shown in the overtopping region. In addition, fluctuating hydrodynamic coefficients are further proposed and investigated. These results are helpful to have a better understanding of the problem and to improve related structural designs.     

Experimental investigation of Reynolds number effect on vortex induced vibration of rigid circular cylinder on elastic supports

Vortex induced vibration (VIV) of circular cylinders is widely accepted as depending weakly on Reynolds number (Re) and strongly on reduced velocity, mass ratio, and damping. Experiments are conducted in the regime right before transition from laminar to turbulent flow (high-Reynolds 2×l0⁴?4×l0⁴<Re<3.5×l0⁵?6×l0⁶) in the Low Turbulence Channel at the University of Michigan. Limited experimental data are available in that regime. This research was initiated in 2004 studying high-damping VIV to convert hydrokinetic energy from ocean/river currents to electricity using the VIVACE Converter.The following experimental observations are made: low and high Reynolds VIV significantly differ. The range of synchronization of the upper branch increases with increase in Reynolds number. Amplitude ratio (A/D) increases with Reynolds number within the upper branch. For high-Reynolds, A/D=1.9 is achieved routinely in spite of high damping. The lower branch disappears, overtaken by extended upper branch. High-Reynolds VIV enters into the 2P domain in the Williamson–Roshko (W–R) map beginning with the initial branch. Hysteresis is not observed in these experiments possibly because parameters remain in the 2P-domain. High Reynolds numbers have a stronger influence than mass ratio on A/D. High-damping high-Reynolds VIV amplitudes are double of those predicted by the modified Griffin plot by extrapolation.     

Numerical benchmark studies on drag and lift coefficients of a marine riser at high Reynolds numbers

Numerical benchmark studies on drag and lift coefficients of a marine riser have been organized by the 27th ITTC Ocean Engineering Committee. The purpose of the studies was to benchmark the capabilities of CFD methods through quantitative comparisons and validation studies against the model test results of a circular cylinder by MARIN. Studies were focused on the drag crisis phenomenon for the stationary smooth cylinder in the critical Reynolds number regime. Eight organizations have participated in the studies by using RANS, DES and LES methods. An overview of the model test results, test cases, submissions and comparison results are presented in this paper. Conclusions and recommendations are made for future studies.     

Forces on a circular cylinder near a plane wall

The drag and lift force are measured on circular cylinders fitted with end plates in a wind tunnel. The gap between the cylinder and the wall, G, the thickness of the turbulent boundary layer along the wall, δ, and the Reynolds number, Re, are varied in the following ranges: 0 < G/D < 2, 0.12 < δ/D < 0.97 and 4.8 × 10⁴ Re 3 × 10⁵. The lift and drag coefficients are presented in terms of a new variable G/δ.

It is found that the lift coefficient is governed by the gap to diameter ratio G/D while the drag coefficient is dominated by the ratio of gap to thickness of the boundary layer, G/δ.