并排方柱绕流的大涡数值模拟及显示

采用大涡模拟的方法对雷诺数为7.5×105并排双方柱绕流流场进行了数值模拟,并对所采用的模拟方法进行了实验验证.利用有限体积法和SIMPLE算法计算程式,对双方柱三维物理模型求解不可压缩的N-S方程.计算得到了方柱绕流的速度场和涡量场以及不同时刻的速度分量分布情况.对方柱后速度场、涡旋形成、脱落以及波动性等问题进行了分析,结果符合物理学规律.     

人工压缩法计算突扩管道流

使用人工压缩法计算了突扩管道流.对计算所得不同雷诺数下的等压线图和流线图进行了对比分析,结果表明,在扩张截面的拐点处形成了低压区,扩张后在侧壁右侧形成了涡,且Re数越小涡越小,Re数越小涡心的位置离侧壁越近.同时在计算中也发现,Re数较大时在出口处将形成回流区,收敛速度也变慢.     

矩形管道内驻波声流的数值模拟

基于边界层理论,建立边界速度数学模型研究矩形管道中Rayleigh声流现象.此数学模型采用3阶的谱元方法求解,驻波声场对流体流动的影响采用壁面处的声边界速度来表达,同时引入雷诺数来分析非线性项和黏性项的重要性.数值结果表明:在2维和3维情况下,声边界速度模型均与近似解相符.声边界速度模型和近似解的差异来源于对非线性项的处理.与近似解相比,声边界速度模型的优势在于能考虑流体流动的非线性效应且仅要求矩形管道的特征尺寸的2倍小于波长.在2维情况下,回流区的涡心位于管道高度的1/4;而在3维情况下,回流区的涡心则靠近壁面.在壁面附近,非线性项的影响不能忽略;而在上下2个涡心的中间位置,非线性项比黏性项更加重要.     

扰流对大口径内锥流量计性能影响的仿真研究

利用数值仿真方法,对800 mm口径、β值为0.726 8的大口径内锥流量计进行了CFD数值模拟研究。分别在内锥流量计的上、下游直管段处安装半月型孔板,通过扰流件对速度场的影响,分析内锥流量计的性能。研究结果表明:大口径内锥流量计流出系数几乎与雷诺数无关;如果在内锥流量计的上游5D、下游1D或更远直管段处安装扰流件,基本上不影响流出系数;若同时在上游和下游安装扰流件,此时下游的扰流件对速度场起到整流器的作用。     

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

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

二维直角管道中高Reynolds流动的格子Boltzmann模拟

运用9 Bit的格子Boltzamnn方法对二维直角管道中高Reynolds流动进行数值模拟, 得到了在不同雷诺数下不可压缩流动涡的涡心的位置分布及其演化斑图.     

流道弯曲对圆柱绕流尾迹流场的影响

为揭示流道弯曲对圆柱绕流尾迹流场的影响,对低雷诺数下弯曲流道内二维圆柱绕流进行数值模拟,讨论管道宽度、弯曲角度和弯曲距离对流场的影响,得到一定管道宽度下多组合弯曲角度和弯曲距离的流场特征参数及漩涡脱落特征。研究结果表明:弯管内圆柱绕流流场是弯管和圆柱之间相互作用的结果,弯曲段的存在使尾迹流场呈不对称分布;流道宽度对流场特征参数有显著影响,流道越窄,流场越复杂,圆柱前后压差越大,涡长约增长1.4倍,分离角约增大3°,圆柱后壁面会出现附壁涡;压力系数、阻力系数平均值和升力系数均方根随弯曲角度增大明显增大,在θ为90°时达到最大值;随弯曲距离增大,弯曲段对圆柱后尾迹流场的影响越来越小。     

叶片前缘旋流蒸汽冷却流动和传热的数值研究

通过求解三维Reynolds-Averaged Navier-Stokes方程和标准k-ω湍流模型,数值研究了旋流蒸汽冷却的基本原理,分析了冷气雷诺数和来流温比对流动和传热特性的影响,旨在阐明旋流蒸汽冷却的原理,总结其流动传热的变化规律。在此基础上对无量纲换热系数Nu、雷诺数Re和来流温比φ进行数值拟合,得到旋流蒸汽冷却的传热关联式。研究表明:冷气在旋流腔内的高速转动引起强烈的径向对流运动,使得换热增强;增大雷诺数能够增大冷气的涡量,有效提高旋流腔的换热系数,同时降低阻力系数;增大来流温比使得冷气的涡量增大密度减小,旋流腔的换热系数略有减小,阻力系数显著降低;综合换热因子随着雷诺数的增大而增大,随着来流温比的增大而减小;拟合的传热关联式与数值计算结果吻合良好,可以准确地预测蒸汽旋流冷却的换热系数。     

用LBM方法模拟壁面驱动粘性不可压半圆形空腔流

基于格子Boltzmann方法(LBM)数值模拟壁面驱动的粘性不可压半圆形空腔流. 采用具有二阶精度的曲线边界处理方法, 得到了不同雷诺数下的流线图、 涡线图及速度分量沿半圆形中心线的分布. 在小雷诺数的条件下, 流动状态仅由一个涡组成; 随着雷诺数的增加, 出现一个二级涡, 涡的大小与雷诺数有关. 数值结果表明, 格子Boltzmann方法简单有效, 适合处理该问题.     

前混合磨料水射流除鳞喷嘴混合腔内部流场

基于计算流体动力学多相流混合物模型,应用FLUENT软件对前混合磨料水射流除鳞喷嘴高压水与磨料混合腔内部流场进行数值模拟.比较了不同进料方式对流场均匀性的影响.分析了两侧为高压水入口条件下,磨料入口直径、高压水入口直径、高压水入口位置和角度以及收缩段锥角对流场混合均匀性的影响,得到了影响混合腔内部流场混合均匀性的合理结构参数.数值计算结果表明:入口速度一定的条件下,磨料中进式喷嘴混合腔的混合均匀性优于磨料侧进式喷嘴混合腔.随磨料入口和高压水入口直径增加,混合腔出口的射流速度均增加,但随高压水入口直径增加导致出口磨料浓度呈先增后减的趋势,磨料入口与高压水入口合理的质量流量比值约为3∶4,两侧高压水入口位置对流场混合均匀性影响较小,高压水入口角度和收缩段锥角均为30°时流场性能更佳.     

封闭方腔自然对流的涡结构和传热特性

为研究封闭方腔自然对流的涡结构和传热特性,对普朗特数,Pr=0.71,方腔高长比A=H/L=1,瑞利数Ra=1.58×109的二维封闭方腔自然对流进行了直接数值模拟. 给出了水平边界层发展阶段和垂直边界层转捩阶段的流场结构以及Nu数分布. 结果表明,压力梯度对方腔的水平速度变化起着决定性的作用;高Ra数下的自然对流在逆压梯度作用下水平方向形成一系列的涡,这些涡使得水平边界层流动形成分层结构;壁面Nu数与速度梯度?v/?y有着密不可分的关系.      

CLT桨的尾流场及梢涡特性数值分析

为了探究尾流收缩叶梢负载(CLT)桨的尾流场及梢涡特性,基于商业计算流体力学软件Star-CCM+对CLT桨P1727的尾流场进行了数值模拟.计算采用DDES(延迟分离涡)方法在进速系数J=0.5的工况下对三套网格进行了不确定度分析,确认了数值方法的可行性,然后对螺旋桨的敞水性能及尾流场进行了计算.研究表明:尾流场的网格细化对螺旋桨的推力及扭矩只有微弱影响,对桨后的流场具有较大影响;尾流区域轴向速度分布分为加速流与自由流,进速系数越大,加速流与自由流的分界线向外扩散的趋势越弱;相邻梢涡之间会发生融合,进速系数越大,融合得越晚,梢涡强度越弱;CLT桨相较于常规桨多出一个端板涡,并且会与梢涡融合;梢涡在向下游发展的过程中,会有明显的收缩现象.     

孔板诱发管道振动的关键影响因素研究

基于计算流体力学（CFD）数值模拟与无量纲的分析方法,结合时均化的Realizable k-ε湍流模型与瞬态的LES湍流模型,进行了孔板诱发管道振动的性能分析与关键影响因素的无量纲分析。结果表明：孔板对流场的扰动主要体现在孔板下游,且扰动大小主要受无量纲参数雷诺数、孔径比的影响;随着雷诺数的增加,管道振动增大,但增幅在减小;管道振动与孔径比之间存在负相关关系。     

一种新型结构叶片对风力机气动性能的影响

针对风力机在旋转过程中产生的叶尖涡影响风力机本身以及下游风力机气动性能的问题,提出了一种控制叶尖涡的策略,以减小叶尖涡对风力机本身及下游风力机气动性能的影响.以PhaseⅥ叶片的1/8模型为原始模型,在叶尖处和轮毂处同时开洞,用管道将洞连接的模型称作新模型.采用数值模拟的方法对来流风速从6m/s到20 m/s的15个工...     

孔板诱发管道振动的关键影响因素

基于计算流体力学(CFD)数值模拟与无量纲的分析方法,结合时均化的Realizable k-ε湍流模型与瞬态的LES湍流模型,进行了孔板诱发管道振动的性能分析与关键影响因素的无量纲分析。结果表明:孔板对流场的扰动主要体现在孔板下游,且扰动大小主要受无量纲参数雷诺数、孔径比的影响;随着雷诺数的增加,管道振动增大,但增幅在减小;管道振动与孔径比之间存在负相关关系。     

应用粒子图像速度场仪测试圆柱绕流涡的结构

利用粒子图像速度场仪（PIV）的2帧互相关技术,在风洞中测量雷诺数为500的圆柱绕流.通过展示一个涡的产生、发展及脱落的连续瞬时速度矢量场、流线图、速度量级等势线图和涡量等势线图,分析圆柱后涡系发展及涡量输运规律.     

有体积热源的矩形空腔内层流自然对流的数值模拟

对具有内热源方腔的稳态层流耦合自然对流换热进行了三维的数值模拟,采用的模拟代码基于连续介质计算力学的开源库OpenFoam,解决了自然对流换热与固体传热的耦合问题.Ra数的变化从105到109.对外壁面为常温、方腔内充满含体积热源流体的自然对流计算结果表明,温度场、速度场与非耦合的工况有很大差异.     

二阶迎风有限体积法方腔流数值模拟

介绍了一种基于非结构网格的二阶迎风有限体积离散格式，在对流项的离散过程中，为了达到二阶精度，在界面上对物理量咖作Taylor展开，在处理展开项中的梯度时经入NND格式的优点，克服了中心差分格式不稳定的缺点，并对方腔流动进行了系统的数值模拟．计算网格采用三角形网格，节点数为12960，单元数为25600．压强修正基于SIMPLEC方法．给出了雷诺数Re达到5000时的定常流动结果，与以往方腔流计算的标准解非常吻合，但所用的网格数要少．计算结果表明该二阶迎风有限体积算法有很好的收敛性和稳定性，且松弛因子对计算结果影响很小．     

应用FPM方法模拟封闭方腔自然对流

为验证有限粒子方法（FPM）能否很好地模拟自然对流问题,采用FPM方法对封闭方腔自然对流问题进行了数值模拟. 对FPM方法进行详细描述,并利用FPM方法对拉格朗日型的N-S方程进行离散,基于离散后的方程对瑞利数Ra分别为104,105,106的封闭方腔自然对流问题进行数值模拟,给出了不同瑞利数条件下的速度和温度云图以及中心线上量纲一的速度和温度分布曲线. 结果表明,FPM方法能够获得较准确的温度分布规律,当瑞利数较低时也能获得较准确的速度分布规律,但随着瑞利数的增大,速度场的模拟结果精度和稳定性变差,因此为了获得更准确和光滑的速度分布场,需要对FPM方法做进一步的改进.      

射流纵向涡强化换热的数值模拟

为研究射流引起的纵向涡对流动和换热性能的影响,采用数值方法模拟了三维矩形通道内有射流从底面进入时的定常、不可压层流对流换热,得到了纵向涡影响下的速度场和温度场以及沿流动方向局部Nusselt数的分布.以场协同原理为指导,分析了射流纵向涡强化换热的原因,并进一步研究了射流角对纵向涡的换热强化效果的影响.结果表明 纵向涡改善了通道内速度场和温度场的协同关系,强化了对流换热; 射流垂直底面入射时,纵向涡的换热强化效果较好.