绕圆柱荷电气固两相流动的数值解法

提出了一种求解绕圆柱荷电气固两相流动的数值解法，即用时间推进法求定常解，用有限差分方法离散控制方程，用人工粘性法处理固相间断面，计算结果表明，此数值解法能够处理强固相间断面，并且可用于含强间断面的悬浮体两相淳动的数值计算。     

可压密相气固两相湍流流动的数值模拟方法

在ＤＩＰＳＡＲ算法的基础上,推导了适用于可压密相气固两相湍流流动数据会计的ＣＤＩＰＳＡＲ算法,利用提出的可压密相的两个ｋ－ε低湍流雷诺数模型,用ＣＤＩＰＳＡＲ算法,对竖直上升管中可压密相气固两相湍流流动进行了数值计算,同时考虑了密度变化和体积分数变化的影响,研究结果表明,随着颗粒体积分数的增加,颗料间及颗粒与壁面间的相互作用占了主导地位,使得单个颗粒脉动速度增大,但是气体及颗粒的速度差却随之减小,     

双拉格朗日模型模拟气固两相双圆柱绕流

以湍流剪切层中大尺度旋涡结构与颗粒群的相互作用及运动规律为研究对象,提出了一种基于离散涡方法与颗粒碰撞模型相耦合的双拉格朗日气固两相流模型.用次循环调整颗粒运动计算的时间步长,将两相计算耦合在一起,从而解决了离散涡方法求解单相流场与颗粒运动方程的时间相容性问题.首先用颗粒在兰金涡中的运动验证了所提的算法,随后将该方法应用于高雷诺数湍流旋涡结构中的颗粒水平输运过程研究,模拟了圆柱尾迹中颗粒的沉降过程及其被旋涡卷吸、扬起和团聚的过程,说明在有大尺度旋涡存在的流场中颗粒更容易被流体向下游榆运.当来流速度与重力方向一致时,颗粒分布特性受颗粒Stokes数(St)影响.当St小于1时,颗粒响应流场变化的时间短,随St增加颗粒对流场变化的响应变得迟钝,重力所起的作用更加明显;当St远远大于1时,颗粒的分布几乎不受旋涡的影响,而与自由沉降的分布十分接近.     

圆柱绕流的三维数值模拟

利用计算流体力学软件CFX－4,对粘性不可压缩流体的圆柱绕流进行了三维数值模拟,采用有限体积法和SIMPLE计算程式,利用不可压缩Navier-Stokes方程,模拟雷诺数在亚临界区内的绕流流动,并计算了流体的水动力特性。为克服数值模拟高雷诺数时的数值不稳定性,计算中采用了QUICK迎风格式,其对流项为三阶精度,其余项如扩散项等为二阶精度,圆柱两端边界采用周期性边界条件。计算结果表明,高雷诺数时圆柱周围的流动具有明显的三维特性,且沿柱长方向不同断面的升力和阻力系数并不相同。同时,对圆柱绕流进行了二维数值模拟,并与三维数值结果进行比较,发现三维模拟的升力和阻力系数均小于二维模拟。     

涡方法数值模拟高雷诺数圆柱绕流

圆柱绕流问题是二维物体粘性绕流的基本问题.二维N-S方程的各种差分求解格式往往难以用来实际计算高雷诺数流动问题。80年代中期以来,Stansby等发展了70年代初Chorin提出的随机点涡法,用网格涡方法计算点涡的对流速度,大大减小了计算量,笔者采用Stansby等发展了的随机点涡法,对雷诺数2000的圆柱绕流进行了数值模拟。     

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

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

圆柱绕流的离散涡数值模拟

首先用面元法模拟了圆柱绕流的势流解,得到了多个圆柱不同配置下的流线图.在此基础上按照一定分离规则引入点涡,就可以用离散涡方法很好地模拟多圆柱绕流.先后得到了流场的流型、涡量等值线图、速度矢量图和阻力升力系数.结果表明单圆柱和双圆柱后的流场有相当大的差别,后者的流型、速度矢量分布都发生了很大的变化,阻力和升力系数也显著增大.     

圆柱绕流远场涡对的数值模拟

基于离散涡方法求得非定常、不稳定流场,数值模拟了三种不同时刻高雷诺数下圆柱绕流结构的发展,从流谱图、等涡量线图和涡谱图可以清晰地看出从近场的初生卡门涡街,过渡到远场的二次涡街的过程,计算结果发现：远场离散涡有形成二个涡的涡对及三个涡的涡对的趋势,计算结果说明了流体运动中涡对结构的本质：由于来流是均匀的,没有加入任何拔动,当流体流过钝体时产生具有剧烈分离的不稳定流动,因此在远场形成的二次涡对及卡门涡     

数值模拟放置附属圆柱的主圆柱绕流

利用格子Boltzmann方法, 对主圆柱尾流区域内放置附属圆柱的绕流进行数值模拟. 结果表明: 放置单个附属圆柱时, 主圆柱所受阻力减小, 阻力
系数发生周期性改变; 放置两个附属圆柱的减阻效果更好.     

离心风机叶轮内气固两相流动的理论及数值解法

利用叶轮机械内Ｓ１、Ｓ２流面迭代方法对离心叶轮机械内气体流场进行数值分析，提出了用三维有限元数值插值技术和四阶Ｒｕｎｇｅ－Ｋｕｔａ法解高阶常微分方程组的新方法，并利用此方法求解离心叶轮机械内三维粒子的运动迹方程．以３００ＭＷ电站锅炉引风机为算例进行气固两相流数值分析，找出了不同颗粒粒径、进口位置对叶轮碰撞和磨损的影响，有助于以气动力学观点确定电站锅炉引风机防磨设计新方法．     

单方柱绕流的直接数值模拟研究

为了揭示柱体绕流的湍流流动机理,采用直接数值模拟方法对雷诺数为1200的单方柱绕流工况进行研究。首先通过与文献中斯特劳哈尔数、平均流速和表面压强系数等统计量进行对比,验证了数值方法的可靠性。其次采用本征正交分解方法系统地提取流场中的相干结构,结果表明：对于速度场的模态分解,第一阶模态代表平均速度场的特征,第二、三阶模态提取的是流场中的低频大尺度旋涡特征,第四、第五阶模态提取的是流场中的高频小尺度旋涡特征。最后基于联合概率密度函数分析速度梯度张量第二、第三不变量,发现方柱下游大致可分为两个流动阶段：发展阶段,流场以涡流层结构和耗散作用为主,涡流管结构逐渐生成;成熟阶段,流场中湍流结构伴随着高涡量拟能和高能量耗散率。     

脉动流条件下圆柱绕流特性研究

该文数值研究了脉动流条件下圆柱绕流场特性,对涡量分布、升阻力系数、升力系数频谱特性等进行了分析。结果表明,脉动流作用会使圆柱尾涡剪切层变薄,提高脉动频率使尾涡长度变短、脱离加快,提高无量纲脉动振幅使主导涡的脱落速度降低,使旋涡生成区域更靠近圆柱表面;脉动频率与无量纲脉动振幅的增大使升、阻力系数的振幅均增大,而且与升力系数相比,阻力系数的振幅更大,变化更快;升力系数频谱图存在多个主频,包括旋涡脱落频率和相位叠加频率,其中旋涡脱落频率的振幅随着脉动频率的增加而减小,随着无量纲脉动振幅的增加而增加。     

绕水平圆柱层流自然对流数值模拟

从Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和能量方程入手，利用有限差分方法摸拟绕水平圆柱定常层流自然对流流动，采用流函数涡量无量纲方程，非均匀网格剖分，二阶精度差分格式，非均匀圆柱表面温度分布，计算了瑞利数Ｒａ在１０４～１０６，普朗特数Ｐｒ在０．７～１００范围内的各种非均匀表面温度分布下的绕水平圆柱自然对流，特别是在圆往表面温度局部低于室温情况下，计算出了回流。全部结果与实际情况及已有的实验计算结果吻合。     

应用粒子图像速度场仪对泊肃叶流动及圆柱绕流的测量

为了考察粒子图像速度场仪（PIV）的测量精度及分辨复杂流动结构的能力，对泊肃叶流动和圆柱绕流两种典型流动进行了测量。将泊肃叶流动的速度廓线的测量值与分析解 进行了对比，发现在测试区域内90％的测量值与分析解吻合良好，具有较高的精度，而近壁面区域由于速度偏振，测量值低于分析解。在此基础上，利用PIV技术测量了Re为500时圆柱绕流流动的主流区、圆柱两侧的剪切层，以及圆柱后的回流区和旋涡区的分布，验证了该技术非常适合于研究复杂的流动结构。     

圆柱绕流的LBM模拟

Lattice Boltzmann Method(LBM)是一种近年来发展的一种数值方法。它具有并行效率高,边界处理简单的特点。本文采用一种能对曲线边界进行较好处理的方法,用LBM对Re=100圆柱绕流进行了计算,计算结果和经典结果一致。进一步,对柱群间复杂流场做了模拟,结果表明,此方法在处理复杂边界是有效的,并且具有较好的并行效率。     

用格子Boltzmann方法计算剪切流的方柱绕流问题

用格子Boltzmann方法计算来流为水平剪切流的方柱绕流问题, 得到了在不同速度梯度条件下方柱绕流的流线和等涡线图. 发现在圆柱尾部形成两排涡, 当来流速度梯度较大时, 两排涡有很大的不同. 计算结果表明, 用格子Boltzmann方法计算剪切流的方柱绕流问题是可行的, 计算结果与理论分析相符.     

圆柱绕流运动的GHM模拟

 基于颗粒流体动力学方法(GHM)对不同雷诺数(Re为80~1646)状态下的圆柱绕流问题进行模拟.结果显示,当Re=80时,在圆柱下游有一对较为对称的附着涡;随着雷诺数的增大,流动变得不稳定,圆柱后方的一对附着涡逐渐脱落消失,尾流逐渐变窄.通过圆柱绕流现象的算例研究,讨论了GHM的特点和需要进一步研究的问题.颗粒流体力学方法将流体视为一个离散的系统,将流场离散为弹性流体颗粒,采用赫兹碰撞理论研究颗粒之间的本构关系,处理颗粒所受到的碰撞力;采用帕斯卡原理和声速导数状态方程研究流场的密度差引起的压强差变化,处理颗粒受到的压差力.颗粒流体动力学方法不同于基于流体为连续体假设所建立的N-S方程,具有一定的原始创新意义.     

二维圆柱层流绕流及其控制数值模拟

采用有限体积法求解二维N-S方程,通过改变时间、空间计算参数,对雷诺数100的二维圆柱非定常流场进行了数值模拟,将所得结果与实验数据以及已有的计算结果进行了对比.并对在尾部多个位置加入不同长度阻隔板控制流动的圆柱流场进行了考察.结果发现:在合适的位置放置阻隔板,可以有效抑制尾迹中涡脱落的形成.     

圆柱面上微小后向台阶紊流绕流场的数值分析

利用分块联合求解的思维,在贴体正交曲线网格上,用ｋ－∈双方程模型和简化的ＳＯＬＡ方法对圆柱面上微小后向台阶的紊流绕流场进行了数值研究,探讨了台阶几何高度、安装位置、绕流雷诺数对流场的影响及圆柱上微上后向台阶扰流减阻的物理机制。