并列旋转双圆柱流动特性的数值模拟

该文基于k - ε湍流模型,采用Galerkin有限元法对并列旋转双圆柱的绕流特性进行了数值模拟,计算的雷诺数为 1550.为了考查两圆柱旋转和间距的相互作用,文中采用三种间距比分别是T/D = 1.2,1.6和3.0 (T 为两圆柱中心之间的距离,D为圆柱直径) 和一系列不同的旋转速度比 (|α| ≤ 2).计算显示,当 |α| = 0,即圆柱不转动时,对应三种间距有三种典型的流型,单钝体流型对应小间距、偏流对应中等间距和对称流对应大间距;当 |α| 达到临界值时,涡脱落得到了有效的抑制,流动趋于稳定,升力系数和阻力系数的脉动值趋于零;平均升力系数和阻力系数随着 |α| 的增大分别增大和减小.     

2-D eddy resolving simulations of flow past a circular array of cylindrical plant stems

In the present study, 2-D large eddy simulations(LES) are conducted for flow past a porous circular array with a solid volume fraction(SVF) of 8.8%, 15.4% and 21.5%. Such simulations are relevant to understanding flow in natural streams and channels containing patches of emerged vegetation. In the simulations discussed in the paper, the porous cylinder of diameter D contains a variable number of identical solid circular cylinders(rigid plant stems) of diameter d= 0.048 D. Most of the simulations are conducted at a Reynolds number of 2 100 based on the diameter D and the velocity of the steady uniform incoming flow. Though in all cases wake billows are shed in the regions where the separated shear layers(SSLs) forming on the sides of the porous cylinder interact, the effect of these wake billows on the mean drag is different. While in the high SVF case(21.5%), the total drag force oscillates quasi-regularly in time, similar to the canonical case of a large solid cylinder, in the cases with a lower SVF the shedding of the wake billows takes place sufficiently far from the cylinder such that the unsteady component of the total drag force is negligible. The mean amplitude of the oscillations of the drag force on the individual cylinders is the largest in a streamwise band centered around the center of the porous cylinder, where the wake to wake interactions are the strongest. In all cases the maximum drag force on the individual cylinders is the largest for the cylinders directly exposed to the flow, but this force is always smaller than the one induced on a small isolated cylinder and the average magnitude of the force on the cylinders directly exposed to the flow decreases monotonically with the increase in the SVF. Predictions of the global drag coefficients, Strouhal numbers associated with the wake vortex shedding and individual forces on the cylinders in the array from the present LES are in very good agreement with those of 2-D direct numerical simulations conducted on finer meshes, which suggests LES is a better option to numerically investigate flow in channels containing canopy patches, given that LES is computationally much less expensive than DNS at high Reynolds number. To prove this point, the paper also discusses results of 2-D LES conducted at a much higher Reynolds number, where the near-wake flow is strongly turbulent. For the higher Reynolds number cases, where the influence of the turbulence model is important, the effect of the sub-grid scale model and the predictive capabilities of the unsteady Reynolds averaged Navier-Stokes(RANS) approach to predict flow past porous cylinders are discussed.     

不等直径串列圆柱绕流大涡模拟

为研究背负式海底管线中增设的小直径附属管线对主管线的水动力影响,将大涡模拟中经典Smagorinsky亚格子模型与特征线算子分裂有限元法结合,并引入出口对流边界条件,完善了基于特征线算子分裂有限元的大涡模拟方法。通过自编程序数值模拟Re=1 000的单圆柱绕流,计算结果与相关文献吻合较好,验证了该算法计算圆柱绕流的有效性,并分析了Re=1 000时不同直径比、间距比情况下的串列双圆柱绕流,根据流场的不同涡脱落形态及两圆柱平均阻力系数、升力系数随直径比、间距比变化的规律得到了不同直径比条件下的临界间距范围。达到临界间距后,流场由单一涡脱落状态转变为双涡旋脱落状态。最后分析了两圆柱平均阻力系数及升力系数在临界间距后急剧增加的原因,为背负式海底管线的布局优化提供了理论依据。     

Numerical simulation of flow past twin near-wall circular cylinders in tandem arrangement at low Reynolds number

Fluid flow past twin circular cylinders in a tandem arrangement placed near a plane wall was investigated by means of numerical simulations. The two-dimensional Navier-Stokes equations were solved with a three-step finite element method at a relatively low Reynolds number of Re = 200 for various dimensionless ratios of 0.25 ≤ G/D ≤ 2.0 and 1.0 ≤ L/D ≤ 4.0, where D is the cylinder diameter, L is the center-to-center distance between the two cylinders, and G is the gap between the lowest surface of the twin cylinders and the plane wall. The influences of G/D and L/D on the hydrodynamic force coefficients, Strouhal numbers, and vortex shedding modes were examined. Three different vortex shedding modes of the near wake were identified according to the numerical results. It was found that the hydrodynamic force coefficients and vortex shedding modes are quite different with respect to various combinations of G/D and L/D. For very small values of G/D, the vortex shedding is completely suppressed, resulting in the root mean square (RMS) values of drag and lift coefficients of both cylinders and the Strouhal number for the downstream cylinder being almost zero. The mean drag coefficient of the upstream cylinder is larger than that of the downstream cylinder for the same combination of G/D and L/D. It is also observed that change in the vortex shedding modes leads to a significant increase in the RMS values of drag and lift coefficients.     

Numerical simulation of low-Reynolds number flows past two tandem cylinders of different diameters

The flow past two tandem circular cylinders of different diameters was simulated using the finite volume method. The diameter of the downstream main cylinder （D） was kept constant, and the diameter of the upstream control cylinder （d） varied from 0.1D to D. The studied Reynolds numbers based on the diameter of the downstream main cylinder were 100 and 150. The gap between the control cylinder and the main cylinder （G） ranged from 0.1D to 4D. It is concluded that the gap-to-diameter ratio （G/D） and the diameter ratio between the two cylinders （d/D） have important effects on the drag and lift coefficients, pressure distributions around the cylinders, vortex shedding frequencies from the two cylinders, and flow characteristics.     

亚临界雷诺数下波浪型圆柱绕流的数值模拟及减阻研究

采用大涡模拟数值方法对亚临界雷诺数Re=3000下波浪型圆柱的绕流现象进行研究分析。研究结果表明,波浪型圆柱的三维尾迹涡结构能得到很好的控制,它在轴向方向呈现周期性正负涡的交替分布特性。随着幅值的增大,波浪型圆柱表面的自由剪切层得以延展,使得旋涡的脱落发生在波浪型圆柱下游较远处,从而达到减阻的目的。波浪型圆柱与普通直圆柱相比,其幅值对平均阻力系数及脉动升力系数的减少及尾迹控制有着更重要的影响,最大减阻可达16%。     

三维圆柱绕流数值模拟湍流方法的选择

研究湍流模拟方法对三维圆柱绕流数值模拟精度的影响,分别采用雷诺平均法（RANS）中的κ-ω模型、SST模型以及大涡模拟法（LES）对亚临界区内雷诺数Re=3 900时的三维圆柱绕流进行数值计算,分析了圆柱体表面的受力、圆柱后流场时均速度特性与瞬时涡量分布情况。结果表明,当流体流过圆柱表面时,圆柱表面出现的与流速方向相反的负压力差区域使流体从圆柱表面分离,引起了不稳定的周期性交替脱落的湍流涡泄,从而在圆柱表面产生周期性波动的升力,同时在圆柱后近流场区域形成回流区。研究还发现,LES法对圆柱体的受力以及流场时均速度特性的模拟效果要优于κ-ω模型与SST模型;相较于前人利用浸入边界法得到的模拟结果,LES法的模拟精度也有了较大提升;同样,通过对瞬时流场涡量等值线图的分析,并与已有的模拟结果进行对比,发现LES法不但可以从整体上表现出漩涡的周期性脱落,而且对流场中不同位置的、复杂的小尺度湍流涡泄也描述得非常细致,得到的自由分离剪切层长度与湍流涡泄的卷曲度更符合湍流涡泄的特征。所以,在亚临界区,LES法对湍流的模拟效果相对较好。     

Large-eddy simulation of the flow past both finite and infinite circular cylinders at Re =3900

The flow past a finite circular cylinder with a height-to-diameter ratio of 1.5 and an infinite circular cylinder of the same diameter at a Reynolds number Re= 3 900 is investigated using the large eddy simulation(LES). The objective of the present study is to explore the differences of the flow mechanisms between the finite and infinite circular cylinders. It is shown that the free end of the finite circular cylinders affects the wake region significantly. The mean drag coefficient and the fluctuating lift coefficient of the finite circular cylinder are smaller than those of the infinite circular cylinder. The three-dimensional separation and the separated shear layer instability of the finite circular cylinder can obviously be observed. The existence of an arch vortex in the average flow downstream of the free end is demonstrated.     

螺旋侧板截面形状对Spar平台涡激运动的影响研究

处于深水中的Spar平台两侧不断产生周期性漩涡脱落,极易诱发涡激运动,缩短平台疲劳寿命。为探究不同截面形状螺旋侧板对漂浮式风力机Spar平台涡激运动的影响,采用大涡模拟(LES)方法,对亚临界雷诺数流动状态下附加不同截面形状螺旋侧板的平台进行研究,并与光滑平台进行对比,分析其水动力系数及尾涡发放特点。研究表明:螺旋侧板在有效抑制平台升力系数的同时,还会增加阻力系数,圆形截面侧板对平台升力系数的降幅最大,阻力系数增幅最小;螺旋侧板的存在使漩涡脱落频率整体减小,从单一主频发放变为主频和次频发放;圆形截面侧板的漩涡脱落点距离平台最远,发放最为迟缓;在研究范围内,发现圆形截面侧板抑涡效果最佳。     

双圆柱在不同排列方式下绕流流态模拟

双圆柱绕流伴随着流动分离、旋涡生成与脱落、旋涡间相互干扰等复杂问题,其流动形态和流动特征受圆柱相对位置影响。使用FLUENT流体软件,选取间距比1.75,2.5和4,在二维层流模型下,模拟了双圆柱串列、30°夹角错置、60°夹角错置和双圆柱并列绕流,分析了双柱绕流流态、旋涡脱落形态、升力、阻力系数随圆柱相对位置改变而变化的规律,并对比已有的试验成果和模拟成果,为桥梁建设和圆柱绕流理论研究提供了基础数据。     

Flow characteristics of the two tandem wavy cylinders and drag reduction phenomenon

This paper presents an extensive numerical study of 3-D laminar flow around two wavy cylinders in the tandem arrangement for spacing ratios (L/Dm ) ranging from 1.5 to 5.5 at a low Reynolds number of 1...     

闸墩绕流影响因素及其作用规律分析

为探究某水利工程渡槽附近不平稳水面的形成原因并为制定消除措施提供参考,使用COMSOL Multiphysics软件对高雷诺数( Re =5×106)条件下,不同宽度的顺直流道、收缩渐变段和弯曲流道下游的圆柱绕流和不同形态的闸墩绕流进行了数值模拟研究,分析了流道布局和闸墩形状尺寸等因素对绕流流场和结构受力特性的影响。圆柱绕流分析结果表明:在计算参数范围内,随着顺直流道宽度的减小,圆柱的升阻力系数均逐渐增大;流道渐变段的截面收缩对圆柱绕流产生的影响不大;流道弯曲导致的来流偏流系数越大,圆柱受力的非对称性越明显,升阻力系数均值越大。闸墩绕流计算结果表明:无墩尾闸墩的绕流流场和受力特性与圆柱绕流类似,随着墩身长度增加,闸墩的平均阻力系数减小;带墩尾的闸墩呈流线型,其升阻力系数与尺寸相当的圆柱和无墩尾闸墩相比均显著减小;给闸墩设置合适的墩尾是解决闸墩绕流脱涡的有效措施。     

LARGE-EDDY SIMULATION OF FLOW AROUND CYLINDER ARRAYS AT A SUBCRITICAL REYNOLDS NUMBER

The complex three-dimensional turbulent flows around a cylinder array with four cylinders in an in-line square configuration at a subcritical Reynolds number of 1.5 × 10^4 with the spacing ratio at L/D = 1.5 and 3.5 were investigated using the Large Eddy Simulation （LES）. The full field vorticity and velocity distributions as well as turbulent quantities were calculated in detail and the near wake structures were presented. The results show that the bi-stable flow nature was observed at L/D = 1.5 and distinct vortex shedding of the upstream cylinders occurred at L/D = 3.5 at Re = 1.5 × 10^4. The techniques of Laser Doppler Anemometry （LDA） and Digital Particle Image Velocimetry （DPIV） are also employed to validate the present LES method. The results show that the numerical predictions are in excellent agreement with the experimental measurements. Therefore, the full field instantaneous and mean quantities of the flow field, velocity field and vorticity field can be extracted from the LES results for further study of the complex flow characteristics.     

Numerical investigation of the concentric annulus flow around a cylindrical body with contrasted effecting factors

This paper studies the wall-bounded flow around a cylindrical at a high Reynolds numbers body in a determined computational domain, with simulations of the 3-D, turbulent concentric annulus flow in a straight pipe. Numerical results show that a reversing zone, appearing as a tongue zone with nested velocities higher than the surrounding area, exists behind the cylindrical body. The annulus space is a region of high velocity and low pressure. The zero velocity, of combined the X- velocity and the Y- velocity, exists in the cross sections and no vortex shedding is formed behind the attaching cylinders. Among all investigated effecting factors, the diameters of the attaching and the main cylinders affect the wake feature behind the cylindrical body while the main cylinder length does not affect the distribution tendency of the flow field. The diameters of the main cylinder and the pipe affect the pressure values and the distribution tendencies on the main cylinder surface. Obviously, the increase of the pipe diameter reduces the drag coefficient of the cylindrical body and the increase of the diameter of the main cylinder increases the drag coefficient greatly. The numerical investigation of the concentric annulus flow provides foundations for further improvements of the intricate flow studies.     

高雷诺数下圆柱绕流与大振幅比受迫振动的数值模拟

利用CFD方法,研究较高雷诺数下圆柱绕流及圆柱大振幅比横向受迫振动特性。在Navier-Stokes方程和k??湍流模型的基础上,对圆柱绕流分别进行了二维与三维的数值模拟,得到了圆柱绕流的升力系数和漩涡发放频率数等结果。研究结果表明,三维模拟的升力系数均小于二维模拟结果。利用Fluent的动网格技术进行大振幅比圆柱横向受迫振动的研究,得到圆柱在不同振动频率时的升力系数及漩涡发放规律和圆柱在特定振幅及频率下的"锁定"区间。     

特征线算子分裂有限元的圆柱绕流大涡模拟

为研究圆柱绕流流场特性,将大涡模拟与特征线算子分裂有限元相结合,建立了大涡模拟特征线算子分裂有限元模型,对单圆柱和串列双圆柱绕流问题进行了数值模拟,所得结果与现有研究结果吻合较好。模拟结果表明:对单圆柱绕流,随着雷诺数的增大,圆柱近尾流区上下交替的涡旋逐渐靠近通过圆柱几何中心的水平线,且涡脱落位置逐渐靠近圆柱。对Re=1 000的串列双圆柱绕流,临界间距在圆柱直径的2.25~2.5倍之间;当两圆柱间距小于临界间距时,上游圆柱后方无明显涡旋脱落,间隙处压力较稳定;大于临界间距时,有涡旋脱落,上游圆柱尾流区上下表面交替出现强负压区。     

利用分离盘控制隔水管涡激振动的数值模拟

涡激振动是导致深海隔水管疲劳失效的重要因素.为有效控制隔水管涡激振动,提出了利用分离盘控制涡激振动的方法.对加装0.25～2.0倍隔水管直径的分离盘后的隔水管流动进行了二维数值模拟.研究了亚临界状态下各种模型旋涡脱落的流场结构、脱落模态、隔水管升阻力系数和旋涡泄放频率的变化特征.结果表明:分离盘加装前与经典的实验和计算结果吻合很好;不同尺寸的分离盘对隔水管涡激振动控制效果具有明显差异;当分离盘的长度为1.0至1.5倍的隔水管直径时,平均阻力系数可以减小20%,尾迹涡街的频率也有了较大的减小,达到控制涡激振动的最优效果.     

高雷诺数下串列双圆柱绕流的影响因素研究

针对串列双圆柱绕流问题进行了数值研究,利用Fluent软件建立二维数值模型,探究了不同柱体间距、不同柱体表面粗糙度对串列双圆柱阻力系数、升力系数、St数的影响。采用SST k-ω湍流模型,对雷诺数为1×10~6,柱体间距L/D分别为1.75、1.90、2.20、2.50,表面粗糙度Rh在0~0.5%D范围内的串列双圆柱进行数值模拟。结果表明:柱体间距L/D在1.75~2.50范围内,光滑柱体表面下,St数、后柱平均阻力系数会随着柱体间距的增加而增大;同时,柱体间距对升力系数振幅也存在极大影响。改变柱体表面粗糙度,会使柱体间距为1.75的串列双圆柱前柱平均阻力系数在Rh=0.2%D时降低7%,而各柱间距下的后柱平均阻力系数会随粗糙度的增加而减小。研究成果能够为河道桥梁建设和圆柱绕流理论研究提供基础数据。     

隔水管附属管控制流动的风洞和水洞实验研究

通过对附属轴向管的隔水管表面压力分布和流体动力特性进行的风洞和水洞实验测试结果的对比,验证了利用低成本风洞实验代替高成本水洞实验来考察隔水管流体动力特性及漩涡特征的可行性。与经典圆柱实验的对比验证了本实验两套测试系统的可靠性。对附加管数N 4和10两种布局的附属管进行流动控制实验的结果表明,4根和10根附属管分别作为少管数和多管数情况的代表,对于流体力的影响差异较大,但均能取得一定程度的降低升阻力和抑制漩涡脱落的效果。多管数对于来流方向的敏感性要低很多。10根管能使Cd和Cl分别获得相对于单隔水管约36%和80%的降低。附属管作为一种消除流向影响以降低流体力和控制漩涡脱落的手段,为抑制隔水管涡激振动提供了参考价值。     

3-D numerical investigation of the wall-bounded concentric annulus flow around a cylindrical body with a special array of cylinders

Concentric annulus flow around a combinational cylindrical body with a special array of cylinders at five high Reynolds numbers is investigated numerically using Fluent 6.3.26 in this paper. The numerical results show a good agreement with the experimental data in regard to the axial velocity of the flow. This study focuses on the flow structure and the hydrodynamic characteristics based on the velocity distribution, the pressure distribution, streamlines and vectors under 1-D, 2-D and 3-D conditions. Meanwhile, some global parameters including the pressure coefficient, the drag coefficient and the lift coefficient are analyzed. Numerical results show that the high velocity region and the reverse wake zone with low velocity exist in some spaces due to the disturbance of the cylindrical body. Negative pressures appear in some regions. Neither a wide area vortex nor the vortex shedding appears in the wall-bounded domain. The drag along the axial direction is the main force acting on the cylindrical body in the pipe domain. The annulus flow around the cylindrical body is analyzed to reveal the hydrodynamic characteristics of the complex turbulent concentric annulus flow field due to the multi-effects in the pipeline.