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

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

Wave–current interaction with a vertical square cylinder

A three-dimensional numerical model is developed in this study to investigate the problem of wave–current–body interaction. The model solves the spatially averaged Navier–Stokes equations. Turbulence effects are modeled by a subgrid-scale (SGS) model using the concept of large eddy simulation (LES). The model is employed to study the wave–current interaction with a square cylinder that is mounted on the bottom and vertically pierces the free surface. The force analysis demonstrates that the presence of waves can reduce both the strength and frequency of vortex shedding induced by a uniform current due to the nonlinear wave–current interaction. The free surface elevation, strain rates of the mean flow, and eddy viscosity are found to closely correlate with the mechanism of vortex shedding. It is also shown that when the vortex shedding is neglected in the calculation such as by the potential flow approach, one may significantly underestimate the magnitude of in-line force. The energy spectral analysis reveals that there exist initiating, growing, and decaying regions for shedding vortices around the cylinder. In the vortex initiating region, both coherent and turbulent structures are nearly two-dimensional that become three-dimensional in the vortex growing region. The kinetic energy of both coherent and turbulent motions is dissipated in the vortex decaying region, within which the mean flow gradually returns back to two-dimensional.     

基于致动线模型的错列式两风机尾流场数值模拟

风电场中风机之间存在十分复杂尾流相互干扰现象,尾流相互干扰效应对风机的功率输出、叶片载荷等产生十分显著的影响。采用致动线模型以及计算流体力学方法,研究两风机之间的复杂尾流干扰效应。在保持两风机纵向间距一样的情况下,考虑两风机在不同横向间距下,数值模拟两风机的部分尾流相互干扰现象,分析两风机的气动功率输出特性、尾流速度变化特性、风轮平面附近轴向诱导因子分布特性,尾涡结构以及尾流干扰效应。数值模拟结果表明:在上下游风机在沿着流向方向距离保持不变的情况下,随着横向间距的变化,上下游风机的尾流存在十分复杂的尾流相互干扰效应,对下游风机的气动功率输出以及两风机风电场的整体流场产生了显著的影响。     

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

指数剪切流作用下海洋立管振动特性研究

细长海洋立管在复杂来流作用下的涡激振动（VIV）是海洋工程领域备受关注的热点问题。采用双向流—固耦合方法,对立管在指数剪切来流作用下的涡激振动实施了数值仿真研究。基于均方根振幅包络图的显性模态分析发现,与线性剪切流工况相比,当流速较小时,立管模型在指数剪切流作用下振幅最大值相对较小;流速较大时,指数剪切流的非线性分布致使立管的振动响应增强。通过对立管模型均方根振幅包络峰、谷处的振动响应频率分析发现,在波峰位置处振动频率单一且较为稳定,而在相邻波谷位置处多频共存现象显著,两者有显著差异。沿管的轴向波形主要表现为驻波主导和驻波—行波混合模式,行波一般间歇发生,流速越大发生频次愈高,其在横流向的传播方向通常由高流速区段传递至低流速区段。指数剪切流作用下,沿管体轴向各截面位置处的旋涡脱落模式差异显著,高流速段尾流区的旋涡发放具有较强的周期性,涡管完整且与管轴存在一定偏离角度,低流速段蜂窝状离散涡旋较多。     

张力腿平台绕流数值模拟分析

在一定来流条件下,张力腿平台(tension leg platform,简称TLP)的立柱后缘出现周期性的交替旋涡脱落,致使立柱受到垂直于来流方向的升力和平行于来流方向的阻力作用,导致TLP产生大幅度往复运动,显著增加平台结构和系泊系统的负载。目前,关于单柱、多柱结构绕流问题的研究较多,但对于TLP绕流特性的研究较少,机理尚存不明确的地方。为研究TLP的绕流力变化情况和流场特征,开展了数值模拟分析。利用计算流体动力学数值模拟软件,基于雷诺平均(RANS)法和分离涡模拟(DES)法对TLP绕流场进行仿真分析,揭示了TLP的绕流特性。结果表明,在3种来流角度和多个折合速度V_r下,TLP绕流的流体力系数存在明显差异,升力系数时域曲线呈现脉动性。TLP的上、下游立柱间存在明显的相互作用,影响了旋涡的形成与发展。TLP的旋涡脱落大多集中在平台固有频率附近,且在V_r=7,来流角度为0°时,升力系数频谱峰值最大,旋涡脱落集中。     

Experimental and numerical study on vortex-induced motions of a deep-draft semi-submersible

An experimental study on vortex-induced motions (VIM) of a deep-draft semi-submersible (DDS) was carried out in a towing tank, with the aim to investigate the VIM effects on the overall hydrodynamics of the structure. In order to study the fluid physics associated with VIM of the DDS, a comprehensive numerical simulation was conducted to examine the characteristics of vortex shedding processes and their interactions due to multiple cylindrical columns. The experimental measurements were obtained for horizontal plane motions including transverse, in-line and yaw motions as well as drag and lift forces on the structure. Spectral analysis was further carried out based on the recorded force time history. These data were subsequently used to validate the numerical model. Detailed numerical results on the vortex flow characteristics revealed that during the “lock-in”, the vortex shedding processes of the upstream columns enhance the vortex shedding processes of the downstream columns leading to the rapid increase of the magnitude of VIM. In addition to the experimental measurements, for the two uniform flow incidences (0° and 45°) investigated, comprehensive numerical data of the parametric study on the VIM characteristics at a wide range of current strength will also serve as quality benchmarks for future study and provide guidance for practical design.     

Vortex shedding suppression for flow over a circular cylinder near a plane boundary

In this study, the Navier-Stokes equations and the pressure Poisson equation for two-dimensional time-dependent viscous flows are solved with a finite difference method in a curvilinear coordinate system. With this numerical procedure, the vortex shedding flow past a circular cylinder near a wall is investigated. The flow is calculated for a broad range of gap ratios for different Reynolds numbers ranging from 80 to 1000. Based on the numerical solutions, the vortex shedding is observed using various methods, and the mechanism for the vortex shedding suppression at small gap ratios is analyzed. The critical gap ratio at which the vortex shedding is suppressed is identified at different Reynolds numbers.     

Effects of mass and damping ratios on VIV of a circular cylinder

In this study the basic characteristics of the dynamic response and vortex shedding from an elastically mounted circular cylinder in laminar flow is numerically investigated. The Reynolds number ranges from 80 to 160, a regime that is fully laminar. The governing equations of fluid flow are cast in terms of vorticity. The two-dimensional vorticity transport equation is solved using a vortex method. Effects of important parameters on the system response and vortex shedding are investigated; these include: mass ratio, damping ratio, Reynolds number and reduced velocity. The numerical results show that a decrease in either the mass ratio or damping ratio of the system can lead to an increase in both the oscillation amplitude and the reduced velocity range over which lock-in occurs. The results also suggest that the mass-damping parameter may characterize the system response adequately, although the effect of changing mass ratio appears to be a little more profound compared to damping ratio. Vorticity contour plots suggest that the vortex shedding occurs in the 2S mode, although a wake structure similar to the C(2S) mode appears at distances 15–20 diameters downstream in the lock-in region. The simulation results are in good agreement with previously published data.     

High aspect ratio (L/D) riser VIV prediction using wake oscillator model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (L/D) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures.     

Numerical investigation of local scour below a vibrating pipeline under steady currents

Local scour below a vibrating pipeline under steady current is investigated by a finite element numerical model. The flow, sediment transport and pipeline response are coupled in the numerical model. The numerical results of scour depths and pipeline vibration amplitudes are compared with measured data available in literature. Good agreement is obtained. It is found that pipeline vibrations cause increases of scour depth below the pipeline. The scour pit underneath a two-degree-of-freedom vibrating pipeline is deeper than that under a pipeline vibrating only in the transverse flow direction. The effects of water depth are also investigated. The present numerical result shows that water depth has weak effect on the scour depth. However it does affect the time scale of the scour. The shallower the water depth is, the less time it requires to reaches the equilibrium state of scour. It is found that the vibration forces vortices to be shed from the bottom side of the pipeline. Then vortex shedding around a vibrating pipeline is closer to the seabed than vortex shedding around a fixed pipeline. This contributes to the increase of the scour depth.     

Viscous Effects on Wave Forces on A Submerged Horizontal Circular Cylinder

Numerical simulations are carried out for wave action on a submerged horizontal circular cylinder by means of a viscous fluid model, and it is focused on the examination of the discrepancies between the viscous fluid results and the potential flow solutions. It is found that the lift force resulted from rotational flow on the circular cylinder is always in anti-phase with the inertia force and induces the discrepancies between the results. The influence factors on the magnitude of the lift force, especially the correlation between the stagnation-point position and the wave amplitude, and the effect of the vortex shedding are investigated by further examination on the flow fields around the cylinder. The viscous numerical calculations at different wave frequencies showed that the wave frequency has also significant influence on the wave forces. Under higher frequency and larger amplitude wave action, vortex shedding from the circular cylinder will appear and influence the wave forces on the cylinder substantially.     

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

利用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的增大而减小。     

波浪与起伏水平板防波堤相互作用数值模拟

利用自主研发的基于紧致插值曲线CIP(constrained interpolation profile)方法的数学模型,开展规则波与起伏水平板防波堤相互作用的数值模拟研究。模型在笛卡尔直角坐标下建立,以CIP方法为流场基本求解器,分步求解Navier-Stokes方程,利用高精度的流体体积类型的THINC/SW (tangent of hyperbola for interface capturing with slope weighting)方法重构自由液面,采用浸入边界IBM(immersed boundary method)方法处理波浪与起伏板防波堤的耦合作用问题,通过动量源项造波方法模拟波浪的产生。重点关注波浪的浅水变形和板两端涡旋脱落的非线性现象,分析不同潜深、波要素下的板周围流场分布、板的运动响应和波浪的反透射系数。结果表明:起伏水平板主要通过能量反射、板上浅水变形和板两端的涡脱落消能,能有效减小板后波高,具有作为防波堤的可行性。     

Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.     

不同节点布置下二维涡激振动的无单元法参数化分析

尝试建立1种基于无单元法的涡激振动数值模拟算法,并给出适用于涡激振动分析的前处理自动布点方法.基于无单元法理论,使用动最小二乘法构造形函数,利用无单元伽辽金法,采用速度和压力分离模式,及手工布点和自动布点2种方法所得出的节点布置形式,对流场控制方程进行空间离散,模拟二维涡激振动的流场形态,并对VIV相关参数进行分析,计算不同节点布置情况下的升力系数(Cl)、曳力系数(Cd)及斯特罗哈数(St),并与物理模型实验结果进行对比.计算结果表明,无单元伽辽金法应用于立管VIV分析是可行的,且文中采用2种布点方法均能较好的模拟流场中泻涡脱落的形态,计算结果与传统方法和物理模型实验结果吻合良好.证明文中2种布点方法都能用于固定圆柱的二维VIV分析,但自动布点法能够更好地适用于复杂问题的计算及圆柱体在流场中的VIV动力响应分析.     

Bifurcation of barotropic flow and homogenization of potential vorticity in a rotating laboratory model

Laboratory experiments with a rotating tank confirm the bifurcation character of a barotropic flow driven by an inflow and an outflow described by Sakai (1986). The model, a circular basin with a topographic β-effect, simulates a mid-latitude oceanic feature. At a low Rossby number, stationary Rossby waves are observed which are symmetrical with a line connecting the inlet and the outlet. As the Rossby number increases, a bifurcation occurs and two kinds of vortex flows are observed. In the vortex, potential vorticity is almost uniform. In addition to the two vortex flows, a jet-like inertial flow can also be observed. In general, thre results of these experiments agree well with those of a low-order model and a numerical model.     

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

利用计算流体力学软件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%。     

Effects of Degrees of Motion Freedom on Free-Fall of A Sphere in Fluid

Free-fall of a sphere in fluid is investigated at a Galileo number of 204 by direct numerical simulations (DNS). We mainly focus on the effects of different degrees-of-freedom (DOFs) of the sphere motion during free-fall. The characteristics of free-fall are compared with those of flow past a fixed sphere. Additional numerical tests are conducted with constraints placed on the translational or rotational DOFs of the sphere motion to analyze different DOFs of sphere motion. The transverse motion contributes significantly to the characteristics of free-fall; it results in the retardation of the vortex shedding, leading to the decrease of the Strouhal number. In addition, the transversal sphere motion exhibits the tendency to promote the sphere rotation. On the contrary, the effects of the sphere rotation and vertical oscillations during free-fall are negligible.     

Numerical Investigation on Vortex-Induced Rotations of A Triangular Cylinder Using An Immersed Boundary Method

A numerical study of vortex-induced rotations(VIRs) of an equivalent triangular cylinder, which is free to rotate in the azimuthal direction in a uniform flow, is presented. Based on an immersed boundary method, the numerical model is established, and is verified through the benchmark problem of flow past a freely rotating rectangular body.The computation is performed for a fixed reduced mass of m~*=2.0 and the structural stiffness and damping ratio are set to zero. The effects of Reynolds number(Re=25-180) on the characteristics of VIR are studied. It is found that the dynamic response of the triangular cylinder exhibits four distinct modes with increasing Re: a rest position,periodic rotational oscillation, random rotation and autorotation. For the rotational oscillation mode, the cylinder undergoes a periodic vibration around an equilibrium position with one side facing the incoming flow. Since the rotation effect, the outset of vortex shedding from cylinder shifts to a much lower Reynolds number. Further increase in Re leads to 2 P and P+S vortex shedding modes besides the typical 2 S pattern. Our simulation results also elucidate that the free rotation significantly changes the drag and lift forces. Inspired by these facts, the effect of free rotation on flow-induced vibration of a triangular cylinder in the in-line and transverse directions is investigated. The results show that when the translational vibration is coupled with rotation, the triangular cylinder presents a galloping response instead of vortex-induced vibration(VIV).