圆柱绕流的三维数值模拟

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

圆柱绕流运动的GHM模拟

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

三维方管中柱体可压绕流的大涡模拟

该文利用大涡模拟,对可压缩流体绕经安装在方管底部的长方体和矩形体的绕流现象进行了三维数值模拟.模拟结果显示了流体绕经障碍物后各涡的形成和耗散过程. 分析了不同Re时柱体绕流流场的结构.发现长方体壁面上的各个涡在小雷诺数时基本稳定,在大雷诺数时涡容易脱落,流场呈现复杂的湍流结构.通过改变柱体障碍物的高度,得到了绕流流场在不同阻塞比下的结构变化,发现湍流现象随着阻塞比的增大而增强.     

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

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

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

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

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

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

高雷诺数下双圆柱绕流诱发振动的数值模拟

基于表面涡方法研究了不同间隙比及归一化阻尼参数下,串列和并列双圆柱的流体诱导振动问题,并计算了圆柱的振动响应、流体力、振动频率等.结果表明:圆柱为刚性时的计算模拟结果与相关文献的实验结果吻合良好;在文中的研究条件下,串列双弹性圆柱具有屏蔽流的特征,涡脱落频率比刚性圆柱的大,其值主要取决于归一化阻尼而非间隙比;流体从并列双弹性圆柱分离的涡是对称的,与并列双刚性圆柱相比其平均力系数略小,但由于流体诱导振动造成的脉动力系数较大,故流体诱导振动对并列圆柱的动态响应有重要的影响.研究结果为亚临界雷诺数范围内热交换器管束的复杂流体诱导振动等问题提供了一种可能的计算方法.     

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

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

圆柱绕流的LBM模拟

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

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

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

方柱绕流的大涡模拟

采用物理分步法对方柱绕流的流动分离现象进行了大涡模拟(LES).物理分步法将N-S方程剖分为对流、扩散及波动过程,对不同的物理过程分别采用适宜的方法进行求解.对流过程采用逆风——中心差分格式求解;扩散过程采用中心差分格式求解;压力方程采用SOR迭代法求解.计算的方柱四周的平均压力系数的分布与物理模型试验结果吻合良好;模拟的流动现象也得到了物理模型试验的证实.     

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

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

可压缩气体合成射流微喷优化设计

借助Fluent软件,基于标准k-ε湍流模型和PISO算法,采用动网格方法构造动边界条件,在可压缩模型下,对二维、黏性、非定常气体合成射流微喷的频率特性及影响射流性能的3个关键因素(驱动频率、振动膜峰值速度和喷口尺寸)进行数值模拟,为合成射流微喷的优化设计提供了依据.结果表明,在可压缩模型下,合成射流微喷的频率特性仿真结果与理论计算结果较吻合,验证了其在可压缩气体模型下仿真的必要性;同时,对任一给定的合成射流微喷,驱动频率和喷口宽度存在最优值.     

电磁力控制圆柱体三维绕流场的脱体涡模拟

利用Maxwell方程组直接数值计算表面包覆电极与磁极圆柱体产生的电磁力分布,将其加入到动量方程中,采用脱体涡模拟(DES)方法,在雷诺数Re=3 900时,对电磁力作用下圆柱体在弱电解质中的绕流场结构及其升阻力特性进行了数值模拟与分析.结果表明,电磁力作用可提高圆柱体边界层内的流体动能,抑制流动分离的产生,减弱圆柱绕流场的三维特性,在电磁力作用参数达到某个临界值后,在圆柱体后方产生射流现象;同时,随着电磁力作用参数的增大,圆柱体压差阻力及其总阻力减小,但摩擦阻力增大,而且电磁力的作用还可以显著减小升力脉动幅值.     

两相流方柱绕流的离散涡模拟

运用离散涡方法对流体的方柱绕流进行了数值模拟,对涡团与颗粒相互作用条件下的流固两相绕流进行了计算,并探讨了不同初始位置处颗粒运动对尾流的形成所起的作用.研究表明,尾流中聚积的颗粒主要来源于方柱两侧,来流中夹杂的颗粒多与方柱两侧涡团作用并在方柱两侧聚积.模拟方法将对航天领域中非线性剥蚀问题的解决具有促进作用.     

室内空气流动的直接数值模拟

应用直接数值模拟(DNS)技术,采用高分辨率的网格,计算了一个三维非定常的室内通风流场.对比激光多普勒仪器的测量数据,DNS显示了优异的预测能力,特别是准确地捕捉速度峰值的能力.其计算结果能够详细地描述出室内空气湍流流动的演变和细微的漩涡结构,为设计室内通风环境和研究室内污染物扩散传输提供了强有力的工具.     

LU-AUSM~ 混合格式在三维可压缩流动数值分析中的应用

将一维AUSM 格式扩展到三维任意曲线坐标 .为了获得高阶精度 ,采用了三阶MUSCL格式 .将AUSM 格式与LU SGS格式结合 ,对凸包通道跨音速无粘流动、平面叶栅跨音速和超音速无粘流动 ,以及喷管超音速粘性流动进行了数值模拟 .计算结果与文献计算结果和试验数据相符很好 .     

可压缩气体合成喷射流数值模拟

为了更深入地研究合成喷射流性能与激励器结构参数的关系,在考虑气体可压缩性质的条件下,对合成喷射流进行了理论分析,得出了合成喷射流性能与激励器结构参数的函数关系式,并采用动边界条件模拟实际压电膜的振动,对合成喷外部射流场和激励器腔体内部射流场进行了联合数值模拟.计算结果与已有的实验结果及理论分析结果吻合较好,验证了该数值模拟方法的可行性.     

车辆外部流场的数值模拟研究

将车辆外部流场视为三维定常可压缩气体的紊流流动来进行数值模拟,求出了车辆外部三维的压力、流速分布情况,获得了车辆外形与气动阻力系数等车辆空气动力学特性的定量关系;探讨了数值模拟的数学模型与计算方法和提高计算精度与加速收敛的措施,并用测试结果验证了计算结果的正确性