绕方柱流脉动压力场的大涡模拟

采用有限差分法，对雷诺数为22000的绕方柱流的脉动压力场进行了三维大涡模拟。采用时间分解算子法，将N－S方程分为对流分步、扩散分步和传播分步，再分别进行计算，在大涡模拟中，用两层模型处理Smagorinsky模型在近壁区模拟结果的发散问题。计算的压力场以及阻力系数、升力系数和Strouhal数均与试验结果较为吻合。     

临界区雷诺数下翼型流场电磁力控制的脱体涡模拟

利用Maxwell方程直接数值计算表面包覆电极与磁极翼型产生的三维Lorentz力分布,将其加入到流体控制方程的动量方程中,采用脱体涡模拟方法对临界区雷诺数下受Lorentz力作用下翼型绕流场进行数值模拟,研究分析了电磁力作用系数和攻角对翼型绕流场结构及其升阻力系数的影响机理和规律。结果表明,Lorentz力可以有效地改善翼型周围的流场结构,达到减阻、增升、消涡以及延缓和抑制其失速的目的,因此是翼型的一种有效流动控制手段。     

亚临界雷诺数下电磁力控制圆柱体绕流场特性的脱体涡模拟

该文利用Maxwell方程直接数值计算表面包覆电极与磁极圆柱体产生的电磁力分布,将其加入到动量方程中,采用脱体涡模拟(DES)方法,在亚临界雷诺数(Re=5×104)下,对电磁力作用下圆柱体在弱电解质中的绕流场结构及其升阻力特性进行了数值模拟与分析.结果显示:电磁力可以增大圆柱体近壁面的流体动能,延缓流动分离,减弱绕流场的三维效应;同时电磁力还可以显著减小圆柱体阻力,抑制其阻力与升力脉动幅值,而且随着电磁力作用参数的增加,圆柱体平均阻力及其升力系数均方值均呈指数衰减特征.研究表明,在亚临界雷诺数下,电磁力可以有效地改善圆柱体绕流场结构,达到减小其阻力并抑制其升阻力脉动幅值之目的,因此是圆柱型结构一种有效流动的控制手段.     

用解NS方程和FW-H积分的混合方法计算圆柱绕流噪声

该文所研究的混合方法通过两个步骤模拟流动的声学远场.第一步,在包含所有声源的近场区域内,通过求解NS方程(对于湍流通过大涡模拟)获得非定常流场参数.第二步,通过带可穿透控制面的FW-H积分方法从流体动力近场获得声学远场.FW-H方程的三种积分形式,即三维的Farassat公式、二维的Guo公式和Locard公式被应用,并相互比较.研究了积分控制面位置和展向复制长度对声学模拟结果的影响.把数值模拟的ReD=150时圆柱层流绕流的声学远场和DNS结果进行了比较;把数值模拟的ReD=3900时圆柱湍流绕流的声学远场和实验结果进行了比较.两项比较都能相符合.通过比较FW-H方程三维积分形式和二维积分形式的结果,表明二维积分形式的结果可用作进一步三维精算前的简便估算.     

槽道中方形障碍物绕流的大涡模拟

本文应用二阶全展开ETG有限元离散格式与大涡模拟相结合的方法对固定在槽道一侧的方形障碍物的绕流进行了数值模拟，计算了雷诺数为40000的湍流情况下的绕流流场，将模拟结果与实验结果及他人的数值结果进行了对比，符合较好，表明大涡模拟与所论的具有较高精度的数值算法相结合，适合于具有大尺度涡的绕流运动流场的分析。     

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

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

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

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

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

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

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

高雷诺数条件下,无论来流是均匀还是非均匀,多方柱绕流情况下方柱受力及尾流的相互干扰都是相当复杂的.为了研究方柱受力及下游尾流的相互干扰,基于大涡模拟紊流模型对后品字等边布置的三方柱绕流进行了数值模拟,数值模拟结果表明,上游两个方柱的阻力系数要明显小于下游方柱,而上游两方柱升力系数远大于下游方柱,且下游方柱的升力系数基本...     

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

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

有限深两层流体中水波与水平圆柱潜体的相互作用

研究了有限深两层流体中表面波模态和内波模态的水波与水平圆柱潜体的相互作用问题。在线性势流理论框架内,建立了两种模态入射波作用下,水平圆柱体绕射势的多极展开求解理论。基于所建立的求解模型,对作用在圆柱体上的波浪力进行了数值计算分析,并且与均匀流中情况进行了比较。结果表明了,在入射波的某个频率范围内,流体的分层效应对圆柱体波浪力的影响是不可忽略的。同时,对流域内的绕射速度矢量场进行了数值分析。     

THE RADIATION AND DIFFRACTION OF WATER WAVES BY A BOTTOM-MOUNTED CIRCULAR CYLINDER IN A TWO-LAYER FLUID

An analytical method was proposed to analyze the radiation and diffraction of water waves by a bottom-mounted circular cylinder in a two-layer fluid. Analytical expressions for added mass and damping coefficients, as well as the wave excitation forces of the circular cylinder were obtained by an eigenfunction expansion method. The hydrodynamic forces on the bottom-mounted circular cylinder in a two-layer fluid include not only the added mass and damping coefficients, but also the wave forces of the surface and internal-wave modes. This is different from the case of a homogenous fluid. Some examples were given, showing that density stratification can have a relative large effect on these hydrodynamic forces over a wide range of frequencies.     

均匀流绕旋转圆柱分离流动的数值模拟

本文利用离散涡模型与边界层理论相结合的方法，研究了高雷诺数下，均匀流绕旋转圆柱的分离流动以及有关动力学特性，并验证了Magnus效应，数值模拟过程中，边界层方程用Keller盒式法求解，离散涡采用涡量均匀分布的圆形涡团模型并以涡核半径随时间增长来模拟涡旋的粘性扩散效应。在计算绕旋转圆柱的流动中，无须人为地引进非对称扰动，到一定时候就能自动形成交替脱落的涡街，旋转圆柱周向速度与来流速度比值α的变化范围为0.05-0.3。计算所得升阻系数、分离点的位置以及压力分布与理论和实验结果相符。     

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

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

圆柱非定常绕流及涡致振动的数值计算

用基于一般曲线坐标系和交错网格的差分法求解原始变量二维不可压粘性流体的N-S方程，计算了雷诺数从100到1×10     

绕圆柱非定常周期性涡旋脱落的数值模拟

利用非定常流函数涡量方程数值模拟圆柱突然起动尾流涡旋的形成及周期性脱落过程。对求解的流函数的一阶导数即速度项采用四阶精度的Hermitian公式，而方程的对流项则采用四阶精度的差分格式，并利用ADI方法迭代求解差分方程组。当雷诺数Re不大于40时，圆柱尾流为附体的两个对称涡，为定常解。当Re大于40后流动为非定常及非对称的，圆柱尾流呈现周期性涡旋交替脱落而形成著名的Karman涡街。选择Re=100为例，在初始条件未加任何扰动情况下，成功地模拟了圆柱非定常涡旋形成与脱落的完整过程（无量纲时间算到t=250及以上）。所计算的阻力系数与实验结果及其它数值方法的计算结果一致。约在t=200形成严格的Karman涡街。对涡量方程ADI求解方法的稳定性进行了分析。对流项采用四阶精度差分格式，若应用于定常问题，将极大提高数值求解的精度，若应用于非定常问题的求解，将对求解精度有所改善，其中时间空间两阶混合偏导数的处理是关键，有待进一步的数值实验。     

NUMERICAL SIMULATION OF AN OSCILLATING FLOW PAST A CIRCULAR CYLINDER IN THE VICINITY OF A PLANE WALL

Oscillating flow around a circular cylinder in the vicinity of a plane wall was investigated by solving the two-dimensional incompressible Navier-Stokes equations with a finite element Galarkin residual method. The effect of the gap G/D between the cylinder surface and the wall on the flow behavior was studied. For the case of G/D 〈 0.25, the periodicity in the flow is attributed to both the outer shear layer instability and the oscillating frequency. As G/D 〉 0.25, vortex shedding occurs and the periodicity in the flow is mainly due to the competition of the oscillating frcqucncy and the vortex shedding frequency from an isolated stationary cylinder.     

三维涡环与自由表面的粘性相互作用

本文用ＶＯＦ方法数值模拟了三维涡环与自由表面的粘性相互作用。重点考察了Ｆｒｏｕｄｅ数、Ｗｅｂｅｒ数、初始扰动及密度分层对涡环与自由表面相互作用的影响。结果表明：由于粘性的存在，使得涡环与自由表面的相互作用与无粘情形相比有许多不同。粘性的作用主要体现在涡量的耗散和新的二次涡环的产生；对自由表面变形而言，Ｆｒｏｕｄｅ数是主要的控制参数，Ｆｒｏｕｄｅ数越大，自由表面变形越显著，而Ｗｅｂｅｒ数的影响作用则     

CALCULATION OF VISCOUS FLOW AROUND CIRCULAR CYLINDER WITH THREE-DIMENSIONAL NUMERICAL SIMULATION

1 .　ＩＮＴＲＯＤＵＣＴＩＯＮＴｈｅｆｌｏｗａｒｏｕｎｄａｃｉｒｃｕｌａｒｃｙｌｉｎｄｅｒｉｓｖａｌｕａｂｌｅｉｎｏｉｌｃｏｍｐａｎｉｅｓｏｒｏｆｆｓｈｏｒｅｉｎｄｕｓｔｒｙｔｈａｔｒｅｌｙｉｎｍｏｓｔｃａｓｅｓｏｎｓｕｂｓｅａｐｉｐｅｌｉｎｅｏｒｒｉｓｅｒｓｔｏｔｒａｎｓｐｏｒｔｇａｓａｎｄｏｉｌｆｒｏｍｏｃｅａｎｇｒｏｕｎｄｓｔｏｌａｎｄｏｒｓｕｒｆａｃｅｐｌａｔｆｏｒｍ .Ｗｅｓｈｏｕｌｄｅｎｓｕｒｅｔｈａｔｔｈｅｓｅｏｆｆｓｈｏｒｅｓｔｒｕｃｔｕｒｅｓａｒｅｒｅｓｉｓ ｔａｎｔａｎｄｓｔａｂ…