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

本文试图通过不同精度的差分格式在较高Re数下的计算结果的对照比较,以期了解在高Re数下不同格式的优劣,达到合理选择计算格式,使其达到既计算精确,又相对节省时间的目的。     

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

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

双圆柱绕流的大涡模拟

基于FLUENT软件平台上的三维大涡模拟方法,对定常来流流场中的单圆柱及不同排列形式的双圆柱进行了水动力计算,并对圆柱周围的瞬时流场、时均流场以及圆柱受力进行了分析.分析结果表明,当Re=6×104时,串列双圆柱中上游圆柱受力明显大于下游圆柱,下游圆柱Gd、Cl的变化频率满足fd=2fl,两圆柱的脱涡频率近似相等,涡街的相位相同;并列双圆柱的Cd基本相等,且变化趋势有一定相似性,脱涡频率均接近单个圆柱,涡街的相位相反;交错45.排列的两个圆柱Cd基本相等,脱涡频率接近,下游圆柱的涡街相位明显滞后于上游圆柱.     

用SST-DES和SST-URANS方法数值模拟亚临界雷诺数下三维圆柱绕流问题

基于剪切应力运输(shear stress transport,SST)的分离涡模拟(detached-eddy simulation,DES)在近壁面的流动区域采用SST模型,而在其他区域采用亚格子模型求解流场。本研究中分别应用SST-DES和SST-URANS方法对亚临界雷诺数Re=3900下的三维圆柱绕流进行了数值模拟,比较分析了两种数值方法下圆柱表面的时均压力分布、圆柱后方的漩涡脱落特征、圆柱下游流场的时均速度剖面分布等,并将数值模拟结果同前人物理试验比较,结果表明SST-DES方法在亚临界雷诺数下的三维圆柱绕流数值模拟问题上比SST-RANS方法具有更好的效果。     

基于LES方法圆柱绕流三维数值模拟

该文采用计算流体软件CFX5中Large-eddy simulation(LES)模型计算了均匀流场中三维圆柱绕流的水动力特性。使用有限体积法对三维N-S方程进行求解。数值模拟着重研究了高雷诺数时展向各截面的压力、阻力、升力及涡管特性。数值计算结果表明:展向各截面柱体受力关于中截面对称且小于二维情况,柱体周围流场呈现明显的三维特性。     

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

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

不同长细比圆柱绕流的大涡模拟

自从人们对层流圆柱绕流现象有了系列研究及清楚的认识后,逐渐把目光投向湍流的圆柱绕流问题,但高雷诺数下的圆柱绕流具有很强的三维特性,因此在数值计算时模型轴向长度的选取是大家关心的问题。该文基于开源软件OpenFOAM运用LES计算了4种不同长细比圆柱的静止圆柱绕流(包括一个二维算例),重点分析和对比了结果中的一些基础参数如St数、平均阻力系数等,与实验结果吻合良好。文中还进一步分析了时均顺流向速度U在流场中沿流向和横向的分布,发现长细比为π的圆柱算例中,计算域的轴向尺度已经能满足圆柱下游靠近圆柱处的流场中的三维结构的发展。     

低雷诺数条件下并列三圆柱绕流的尾流模式和水动力系数研究

该文应用浸入边界法对雷诺数Re(28)100和间距比T/D(28)1.1-5.0条件下的等间距并列三圆柱绕流进行了数值模拟。研究发现,尾流模式与无量纲圆柱间距比T/D密切相关,在所研究的间距比范围内共出现六种尾流模式,分别为:单体尾流模式(T/D(28)1.1-1.2)、偏斜尾流模式(T/D(28)1.3)、第一类Flip-Flopping(FF1)尾流模式(T/D(28)1.4-1.5)、第二类Flip-Flopping(FF2)尾流模式(T/D(28)1.7-2.3)、同步尾流模式(T/D(28)2.5-3.2)和调制尾流模式(T/D(28)3.3-5.0)。通过分析近尾流的涡量图、升阻力系数时程曲线和不同圆柱所受升力的相位差,详细研究了各个尾流模式的特点和联系。     

两层黏性流体中直立圆柱体绕流的三维数值模拟

研究了两层流体中直立贯底圆柱体的三维黏性绕流问题.以不可压缩Navier-Stokes方程为控制方程,应用VOF方法追踪两层流体的内界面,建立了该问题的数值模拟方法.成功地数值模拟了两层流体中圆柱体黏性绕流产生的三维尾涡特征,表明了流体的密度分层效应对直立贯底圆柱体的尾涡特性和阻力系数都是有影响的.在海洋立管涡激振动的研究中,考虑流体的密度分层效应是重要的.     

圆柱尾迹流场中横向振荡翼型绕流的数值模拟

该文数值模拟了静止圆柱尾迹中横向振荡翼型绕流的涡结构、力及频谱特性,揭示前方尾迹流场对振荡翼绕流特性的影响,以了解鱼类游动时自主利用前方流场有利干扰的流体力学机理。研究表明圆柱尾迹的影响可导致翼型前缘涡脱落模态发生变化;翼型尾迹的频谱特性更加复杂;在一定的参数下翼型阻力系数可以整体下降。     

均匀流中横向振荡圆柱绕流场的数值分析

针对工程中关心的振荡圆柱绕流问题,采用同位网格系统中的SIMPLE算法求解N-S方程,以及动网格技术模拟圆柱的运动边界。对均匀来流中的横向振荡圆柱绕流场进行数值模拟,计算了振荡圆柱升力系数、阻力系数随时间的变化曲线,以及不同相位时的流线矢量图和涡量等值线分布图。为验证模型的可靠性,对静水中的振荡圆柱流场进行模拟。与实验及数值计算结果吻合较好。     

Numerical simulation of turbulent flow around a forced moving circular cylinder on cut cells

Fixed and forced moving circular cylinders in turbulent flows are studied by using the Large Eddy Simulation （LES） and two-equation based Detached Eddy Simulation （DES） turbulence models. The Cartesian cut cell approach is adopted to track the body surface across a stationary background grid covering the whole computational domain. A cell-centered finite volume method of second-order accuracy in both time and space is developed to solve the flow field in fluid cells, which is also modified accordingly in cut cells and merged cells. In order to compare different turbulence models, the current flow past a fixed circular cylinder at a moderate Reynolds number, Re = 3 900, is tested first. The model is also applied to the simulation of a forced oscillating circular cylinder in the turbulent flow, and the influences of different oscillation amplitudes, frequencies and free stream velocities are discussed. The numerical results indicate that the present numerical model based on the Cartesian cut cell approach is capable of solving the turbulent flow around a body undergoing motions, which is a foundation for the possible future study on wake induced oscillation and vortex induced vibration.     

低Reynolds数横向排列双圆柱绕流的POD-Galerkin谱方法数值模拟

该文将恰当正交分解(POD)方法应用到低Reynolds数(Re)二维横向排列双圆柱绕流问题.采用Fluent数值模拟得到Re=90时的流场数据库,在此基础上采用snapshot的方法构造了POD基.基于此POD基,构造了求解二维不可压缩Navier-Stokes方程组的低维Galerkin谱方法.采用该低维数值模拟方法,再现了低Reynolds数横向排列双圆柱绕流的多种复杂流动模式,其结果与已有的高精度直接数值模拟的结果符合较好,而计算量却非常小,显示了POD方法的有效性和优越性.     

Vanishing of three-dimensionality in the wake behind a rotationally oscillating circular cylinder

The flow behind a three-dimensional rotationally oscillating circular cylinder was studied by a numerical method. The computations were performed at a Reynolds number of 260, which is at a level that the flow wake has developed into a three-dimensional state called Mode-B. The purpose of this paper is to examine the influence of various rotational amplitudes (0.1-0.7) on the wake instability of the flow, while the oscillation frequency is fixed to the value of that measured in the wake of a stationary cylinder. The results show that the rotation with sufficiently high amplitude brings the flow back to its nominal two-dimensional state. Moreover, it is found that the value of the time-averaged drag and the RMS value of the lift are larger than those of a stationary circular cylinder.     

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.