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

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

淹没垂直圆柱周围水流结构及紊流特性试验研究

淹没垂直圆柱虽在各工程领域中均有广泛的实际应用,但还没有关于其绕流场内水流结构及紊流特性的系统研究。通过室内变坡水槽试验对淹没垂直圆柱周围的瞬时流速场进行了精细测量,分析了淹没率对垂直圆柱上下游竖直对称面内三维流速、紊动强度等水力参数垂线分布的影响程度,研究了向下水流、马蹄形旋涡和尾流旋涡等典型水流结构的冲刷力随淹没率的变化规律,揭示了淹没率影响垂直圆柱周围清水局部冲刷的机理。研究结果表明:增大淹没率能够减小垂直圆柱出露于水体中的有效长度,使得垂直圆柱对水体的阻流和扰动作用减弱,也使得垂直圆柱上游迎水面处向下水流、柱周马蹄形旋涡和下游尾流旋涡的冲刷能力被有效削弱,降低了垂直圆柱周围床面泥沙被旋涡体系卷起后顺水流失的概率。     

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

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

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

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

高雷诺数下双圆柱绕流的数值模拟

本文使用表面涡法研究高雷诺数下不同排列方式双圆柱绕流的流动状态，计算了双圆柱在并列、串列及级列的情况下的各种流动结构，涡街的变化及作用在圆柱上的受力情况。本文结果清楚地描述了双圆柱绕流复杂的流动状况，计算结果与实验显示的流动状况十分相似，斯特罗哈数和阻力系数与实验结果符合得很好。     

同轴套管间涡旋流动及强化传热数值模拟

相对旋转两同轴套管间的涡旋流动,能够带来二次流强化传热传质作用,在航空、水处理、生态保护、生物工程和膜分离等领域都具有广泛的应用价值。本文使用Fluent软件,对同轴套管间涡旋流动及传热特性进行数值模拟,考察了内管转速、内外管壁面温差等操作参数变化对同轴套管间流体传热性能的影响,分析了涡旋流动与传热效率之间的关联关系。模拟结果表明:内管转速增加在流场中形成泰勒涡,涡流扰动增大了高温壁面与流体间的热流密度,增强了流体传热效率。增大内外管壁面温差,也可加强流体传热性能,但其强化作用不及内管转速的强化作用显著。受流场中泰勒涡影响,流体速度、温度及热流密度沿轴向的分布都呈正弦状周期性波动,在相邻两涡交界面处,流体传热性能最好,在涡中心处的传热性能最差。     

Multi-mode of vortex-induced vibration of a flexible circular cylinder

The vortex-induced vibration of a flexible circular cylinder is investigated at a constant Reynolds number of 1 000.The finite-volume method on moving meshes is applied for the fluid flow, and the Euler-Bernoulli beam theory is used to model the dynamic response of a flexible cylinder.The relationship between the reduced velocity and the amplitude response agrees well with the experimental results.Moreover, five different vibrating modes appear in the simulation.From the comparisons of their vortex structures, the strength of the wake flow is related to the exciting vibrating mode and different vortex patterns arise for different vibrating modes.Only 2P pattern appears in the first vibrating mode while 2S-2P patterns occur in the other vibrating modes if monitoring at different sections along the length of the cylinder.The vibration of the flexible cylinder can also greatly alter the three-dimensionality in the wake, which needs further studies in our future work, especially in the transition region for the Reynolds number from 170 to 300.     

FLOW PAST TWO ROTATING CIRCULAR CYLINDERS IN A SIDE-BY-SIDE ARRANGEMENT

Measurements were performed using Particle Image Velocimetry(PIV) to analyze the modification of flow by the combined effects of the rotation and the Reynolds number on the flow past two rotating circular cylinders in a side-by-side-arrangement at a range of 425 ≤ Re≤ 1130,0 ≤α≤ 4(α is the rotational speed) at one gap spacing of T / d =1.11(T and d are the distance between the centers of two cylinders and the cylinder diameter,respectively).A new Immersed-Lattice Boltzmann Method(ILBM) scheme was used to st...     

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

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

双浮体自由垂向振荡响应特性

双圆柱体系统在波能装置和防波构件中得到应用，该系统是一个多自由度受迫阻尼振荡系统，阻力和激励力都来自于波浪力。本文采用线性波理论和特征函数展开皮配法，在平面入射波条件下，求解了双圆柱体垂荡运动辐射速度势和圆柱体固定不动时的绕射速度势，推导了附加质量、阻尼系数、波浪激励力表达式，列出了浮体垂荡运动微分方程。最后讨论了下浮体尺度及位置对系统响应特性的影响。数值计算表明无论下浮体尺度及位置如何变化，它只影响到浮体垂荡响应在低频端的值，下浮体的振幅值总是小于入射波的振幅值，而且它在一个较窄的频率范围内从较大的值很快变化接近于零。下浮体半径与上浮体半径之比的变化不仅影响到两个浮体的响应幅值，而且也影响到曲线的变化趋势；两浮体距离的变化对曲线的变化趋势影响不大，但对曲线的幅值有影响；下浮体厚度的变化只少许改变了响应曲线的幅值。     

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.     

明渠柱体上游马蹄涡的运动学特征研究

马蹄涡是水下柱体结构周围床面冲刷的主要动力。为掌握马蹄涡的运动学特征,引入了基于旋转强度的马蹄涡识别方法,提出了基于椭圆涡拟合的马蹄涡特征提取方法,其中,椭圆涡由Oseen涡和倾斜单向剪切流叠加而成。基于上述方法,通过开展明渠圆柱绕流PIV实验,对柱体上游对称面时均流场中马蹄涡的运动学特征进行了研究。结果表明,在柱体雷诺数大于104的紊流柱体绕流中,马蹄涡上游的流动分离点随柱体雷诺数的增大而缓慢向下游移动,但马蹄涡的位置保持不变,涡中心距柱体中轴约0.67D(D为柱体直径或宽度),距床面约0.06D;随着柱体雷诺数的增大,马蹄涡的半径不变而强度增加,使得马蹄涡作用于床面的剪切应力增大。研究结果建立了客观提取马蹄涡运动学特征的数学方法,得到了流动分离点、马蹄涡位置等运动学特征随水流条件的变化规律,揭示了明渠柱体绕流与其它流动中柱体绕流的差异,初步构建了马蹄涡的动力作用概念模型。     

RANS simulations of the U and V grooves effect in the subcritical flow over four rotated circular cylinders

This paper presents a CFD study about the effect of the V and U grooves in the flow over four cylinders in diamond shape configuration at subcritical flow conditions(Re= 41 000). The k- ? Realizable turbulence model was implemented to fully structured hexahedral grids with near-wall refinements. Results showed that the numerical model was able to reproduce the impinging flow pattern and the repulsive forces present in the lateral cylinders of the smooth cylinder array. As a consequence of the flow alignment induced by the grooves, a jet-flow is formed between the lateral cylinders, which could cause an important vortex induced vibration effect especially in the rear cylinder. The magnitudes of the shear stresses at the valleys and peaks for the V grooved cylinders were lower than those of the U grooved cylinders, but the separation points were delayed due the U grooves presence. It is discussed the presence of a blowing effect caused by counter-rotating eddies located near the grooves peaks that cause a decrease of the shear stresses in the valleys, and promote them at the peaks.     

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.     

OSCILLATION AMPLITUDE OF THE FREE END OF A FLEXIBLY MOUNTED RIGID CYLINDER IN TRIANGLE ARRANGEMENT

In this paper,the oscillation amplitude of the free end of a group of threecylinders are studied.The three cylinders are erected at the vertices of an equilateral trian-gle,and three arrangements are made in this study.One of the three cylinders is served asa working cylinder which is mounted by a flexible cantilever and two accelerometers areattached on the free end in two directions.The results of the three-cylinder arrangementare compared with those of two cylinders in tandem/side-by-side arrangement in the samecondition.The experimental results show that the rear cylinder of a group of three cylin-ders is affected more seriously than that of the corresponding two cylinders in the contrastarrangement.The main direction of oscillation is changed.The transversal oscillation isstronger than the longitudinal oscillation in the three-cylinder arrangement.On the con-trary,the transversal oscillation is much weaker than the longitudinal oscillation in two-cylinder arrangement.The largest transversal oscillat     

关于卡门涡街形状稳定性的一点分析

在圆柱绕流的卡门涡街中,单个涡对于来流有相对速度,以该速度观察,其周围的流线更类似无旋理论中有环量的圆柱绕流.因此在势流解的基础上,运用扰动理论分析了涡街形状的稳定性问题.结果表明,涡位置的纵横比有一变化范围,以卡门的经典值为下限,与涡相对速度,实际涡核的大小及其环量都有关系,但与不构成单调的关系.与实验数据和数值模拟的结果进行对比,表明论文的分析是合理的.     

Vortex shedding in the flow around two side-by-side circular cylinders of different diameters

In this paper,the 3-D turbulent flow around two side-by-side circular cylinders of different diameters,at sub-critical Reynolds number(Re=3 900),is numerically simulated by the large eddy simulation(LES).The spacing ratios(T/D)between the two cylinders are considered in four cases(T/D=1.2,1.5,1.8 and 2.7)to study the vortex shedding and turbulent properties in the flow field.The main results are focused on the drag and lift coefficients,the vortex shedding frequency,the coherent structure,and the scale properties.It is shown that when T/D is equal to 1.2,the vortex shedding of the main cylinder is strongly suppressed by the small cylinder,the drag and lift coefficients of the main cylinder are smaller than those in other three cases.While T/D is equal to 1.5,the vortex shedding of the main cylinder can be improved,the drag and lift coefficients of the main cylinder are larger than those in other three cases.The empirical mode decomposition(EMD)method is applied to decompose the velocity signals traced by the LES.It is shown that there is a linear relationship between the mean period and the mode in the semi-log coordinates.The vortex shedding period of the main cylinder is consistent with the period of the restructured coherent structures quantitatively.     

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

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.     

三方柱绕流的大涡模拟及频谱分析

高雷诺数条件下,无论来流是均匀还是非均匀,多方柱绕流情况下方柱受力及尾流的相互干扰都是相当复杂的。为了研究方柱受力及下游尾流的相互干扰,基于大涡模拟紊流模型对后品字等边布置的三方柱绕流进行了数值模拟,数值模拟结果表明,上游两个方柱的阻力系数要明显小于下游方柱,而上游两方柱升力系数远大于下游方柱,且下游方柱的升力系数基本上在零左右振荡,通过频谱分析发现了它们之间的区别。方柱的尾流存在着强烈的非对称相互干扰,且涡有明显的三维性。