低Reynolds数下单排方柱绕流尾迹涡型的LBM模拟

针对工程中遇到的周期柱体绕流问题,采用格子Boltzmann方法(LBM)对低Reynolds数下单排周期排列方柱的定常绕流进行数值模拟,详细分析了柱体后面的尾迹涡型随着Reynolds数的变化情况,并通过对不同Reynolds数下尾迹涡型结构的分析,得到了二射流、三射流和四射流发生合并现象的临界Reynolds数,其中二射流出现合并的临界Reynolds数与文献的数值模拟结果一致,三射流和四射流出现合并的临界Reynolds数则偏大,但更接近实验结果。     

阻流比对平行壁面间圆柱绕流尾迹演化的影响

为了揭示阻流比对圆柱绕流尾迹演化特性的影响,采用多块O型网格,对低雷诺数下较宽阻流比范围的平行壁面间圆柱绕流进行数值模拟研究,讨论不同阻流比下壁面对圆柱绕流尾迹形态和尾迹转捩的影响,得到尾迹转捩临界雷诺数、双子涡和旋涡脱落的特征参数。研究结果表明:圆柱绕流尾迹转捩和尾迹形态受阻流比影响显著,尾迹转捩的第一、第二临界雷诺数随阻流比而增加,当阻流比超过临界值后,尾迹中出现2个附着于壁面的回流区;斯特劳哈尔数和阻力系数均值随阻流比而增加,而升力系数均方根在不同阻流比范围内则呈现出不同的变化规律。     

平面自由湍射流拟序结构的大涡模拟研究

对空间发展的平面不可压缩湍射流拟序结构的非定常演化过程进行了大涡模拟。采用标志物浓度等值线分布示踪平面射流气相流动 ,Reynolds数为 1130 0 ,模拟结果再现了平面射流中 Kelvin- Helm holtz不稳定性的发生以及展向大尺度涡的卷起、合并、破碎过程。捕获到了射流拟序结构剧烈相互作用的“偶极子”和“三极子”现象。在拟序结构的配对过程中 ,大尺度旋涡的尾迹混合形成了流向“发卡”型拟序涡结构。平面射流拟序结构的时空演化过程所呈现出的卡门涡街特征主要集中于初始段的后部和过渡段的前端。数值流场显示结果和染线法的实验结果吻合得很好     

不规则柱体绕流的数值研究

为研究不规则柱体绕流的流体力学特性,介绍了国内外柱体绕流研究的相关进展与研究方法,然后运用数值模拟软件STAR-CCM+建立了4个二维模型,并在雷诺数Re=8×10~5和3.6×10~6两种情况下进行了数值模拟。通过对比普通圆柱体绕流和改变截面形状后的柱体的绕流效果,发现不规则柱体对绕流产生了一定的影响,能够延缓尾涡的脱落,效果显著。     

超临界雷诺数下Spar平台大直径柱体绕流的三维数值模拟

基于Smagorinsky-Lilly模型,采用大涡模拟方法,研究了超临界雷诺数(Re)下大直径柱体的三维黏性绕流问题.在Re=6×107下,对某经典式Spar平台大直径柱体在均匀来流中的三维粘性绕流特性进行了数值模拟与分析.结果表明,尽管来流是均匀的,但在这种超临界Re下大直径柱体尾涡会产生显著的三维结构,使得沿柱体表面压力及阻力系数的分布特性发生明显变化.     

自动网格体系在柱体绕流大涡模拟中的适用性评估

为了评估基于snappyHexMesh方法生成的自动网格体系在二维柱体绕流大涡模拟中的适用性,比较了该自动网格体系与人工网格体系对于Re为3 900圆柱绕流和Re为22 000方柱绕流的数值模拟结果。通过设置合理的计算域以及数值格式,采用snappyHexMesh自动网格以及人工网格的算例都表现出良好的数值稳定性。将不同网格体系的数值模拟结果与物理试验结果进行对比,结果表明,采用snappyHexMesh网格可以提高数值求解效率;圆柱绕流对网格体系的变化比较敏感,不同密度的snappyHexMesh网格会显著影响圆柱气动力特征以及尾流区域的流场结果;snappyHexMesh网格体系可以准确预测方柱绕流,在方柱绕流大涡模拟中具有相较于圆柱绕流更好的适用性。     

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

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

用Lattice Boltzmann方法计算矩形柱的绕流问题

采用二维9速度Lattice Boltzmann模型, 数值模拟流体流过矩形截面柱体的绕流问题. 在边界处, 采用二阶精度的插值边界处理方法, 计算了流动的Strouhal数及柱体受到的升力和阻力系数, 给出了流场的流线和等涡线图. 使来流方向与矩形柱的长边方向平行, 计算结果表明, 改变矩形的长/宽比, 流场的Strouhal数随长/宽比呈线性变化.     

绕水翼超空化流场的数值模拟

基于Reynolds平均法,采用RNG k-ε模型及Mixture两相流模型,对绕水翼ys930的超空化流动进行了数值模拟,得到了绕水翼各个超空化阶段的超空化形态及各阶段的临界空化数与水翼攻角、来流雷诺数之间的关系,确定了临界空化数与来流雷诺数在不同攻角条件下的0.4次方经验关系式,给出了系数的取值区间.结果表明:超空化的发展经历了超空穴形成、汽液两相共存和完全发展3个阶段;各阶段的临界空化数随水翼攻角的增大而增大,基本呈线性关系;由于回射流的作用,在水翼前缘处形成一个逆时针旋涡,使逐渐长大的空泡脱离水翼吸力面并向下游脱落;在不同空化数下,空泡脱落的位置不同,导致空泡长度的变化.     

介质阻挡放电等离子体对圆柱绕流尾迹区流场影响实验研究

利用二维粒子图像测速系统研究了低速风洞实验中介质阻挡放电等离子体对圆柱绕流尾迹区流场的影响。对圆柱尾迹区时均流场和瞬时流场的分析表明,在圆柱表面布置的两对等离子体激励器可在圆柱尾迹区形成射流,进而改变尾迹区的速度分布和涡量分布。     

近壁圆柱绕流积沙线的形成机理研究

应用粒子图像测速(PIV)系统进行近壁圆柱绕流实验,探讨壁面积沙现象的机理.分别在两个动量损失厚度雷诺数1 092和796下,对两个不同直径的圆柱,取间隙比(圆柱与壁面间的距离与圆柱直径之比)0到1.5中的12个值进行实验,总结了不同情况下积沙线的形成特点.研究结果表明:积沙线形成的可能原因是由于圆柱对平板壁面边界层的影响,使得边界层流动在圆柱下游一到两倍圆柱直径的平板壁面处发生壁面分离,在分离区域流体速度减小,并且产生附着涡,流体运动挟带到此处的粒子速度减小,发生沉积而形成.     

壁面加热对微圆柱漩涡脱落的影响

建立二维数学模型,对比研究了低雷诺数下（0~100）壁面加热对直径分别为2 cm、20 μm及2 μm的微圆柱漩涡脱落的影响规律。深入探讨了直径、壁面与流体温差对圆柱绕流中漩涡生长和脱落规律、尾流区速度、温度场及涡街参数分布等的影响规律及其原因,并基于计算结果拟合出了低雷诺数下斯特劳哈尔数（Sr)与雷诺数（Re)的定量关系式。研究发现当圆柱直径降至2 μm时,与常规尺度圆柱相比微圆柱尾流区涡街的出现明显提前;而壁面加热减小了三者的差距,使三种不同直径圆柱背风区涡街出现的最低Re均降至40以下。微圆柱涡街的漩涡列距和间距之比h/l的值随着温差的增加而逐渐增大;相同温差下不同直径圆柱尾流区涡街h/l均随Re增加而减小,在相同Re下h/l值随圆柱直径的减小而明显增大。存在壁面加热时微圆柱绕流的Sr要高于无加热工况,两者之差最高可达20%左右。     

求解不可压缩流动的分步有限元格式

提出求解不可压缩 Navier- Stokes方程的分步有限元格式 ,该格式没有高阶微分项产生、程序编制简单 ,适用于非线性的多维复杂流动。应用该方法实际模拟了二维圆柱绕流的旋涡形成与脱落过程 ,得出了不同 Re情况下圆柱绕流的流速分布。计算得到的不同 Re下的旋涡脱落频率(Strouhal数 )与前人已有的经典解答符合良好     

明渠湍流中的主要相干结构模式

利用自主开发的高帧频明渠湍流粒子图像测速(PIV)系统测量了3种Reynolds数下的恒定均匀流时间序列流场,运用本征正交分解能够给出数据集的最优模态的特性,分析了明渠湍流的主要相干结构模式及其能量关系。湍动能主要集中在少数几个低阶模态上,第1阶模态含有约30%以上的湍动能。Q2(Q4)事件是明渠湍流中除平均流动外含能最多、最主要的大尺度结构,湍动能主要通过Q2(Q4)事件向其他结构传递。它并不完全是发夹涡群诱导的结构,相反从能量的角度来看,发夹涡群的产生与维持与它有关键关系。Q2(Q4)事件可能存在大范围独立维持的大尺度机制。同时本征正交分解(POD)技术成功提取了发夹涡群结构,流向范围达到约3.4h,垂向尺度从交叠层一直到水面附近,整个结构与壁面呈约10°夹角。发夹涡群载有约23%的平均湍动能,并主要集中在处于对数区的中间部分。其他含能较少的大尺度结构有水面涡旋、单双涡结构、独立涡包等。     

Spanwise flow and the attachment of the leading-edge vortex on insect wings

The flow structure that is largely responsible for the good performance of insect wings has recently been identified as a leading-edge vortex. But because such vortices become detached from a wing in two-dimensional flow, an unknown mechanism must keep them attached to (three-dimensional) flapping wings. The current explanation, analogous to a mechanism operating on delta-wing aircraft, is that spanwise flow through a spiral vortex drains energy from the vortex core. We have tested this hypothesis by systematically mapping the flow generated by a dynamically scaled model insect while simultaneously measuring the resulting aerodynamic forces. Here we report that, at the Reynolds numbers matching the flows relevant for most insects, flapping wings do not generate a spiral vortex akin to that produced by delta-wing aircraft. We also find that limiting spanwise flow with fences and edge baffles does not cause detachment of the leading-edge vortex. The data support an alternative hypothesis-that downward flow induced by tip vortices limits the growth of the leading-edge vortex.     

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

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

数值分析串列双圆柱绕流

为研究不可压缩流动中二维的串列双圆柱中不同间距对圆柱间相互作用和尾流特征的影响,选取间距比L/D(L为两圆柱中心间的距离,D为圆柱直径)为2、3、4共3个间距进行了数值分析.计算均在Re=200的非定常条件下进行.计算了圆柱的升阻力系数、尾涡脱落频率等描述绕流问题的主要参量.并与参考文献的数值结果比较,验证了用FLUENT分析的有效性.     

翼尖涡多阶段演化过程及其对气动力的影响

针对高雷诺数和中等马赫数下翼尖涡的近场演化问题,以NACA0012机翼为对象,采用大涡模拟方法,研究了三组不同的马赫数(0.3、0.45、0.6)和雷诺数(5×10~5、1×10~6、2×10~6)下翼尖涡中主涡和次级涡的演化特性以及其对机翼气动力的影响。研究发现,根据主涡和次级涡特性可将翼尖涡近场演化过程分为三阶段:在第一阶段中二者独立生长,主涡涡核处涡量先增后减,次级涡涡核处涡量和流向速度显著变化;在第二阶段中次级涡运动至机翼上表面与主涡相互作用融合,二者涡核处涡量变化分别趋于平缓并最终相同;在第三阶段中主涡与次级涡融合后的共转涡和新生成的二次融合涡离开机翼进入尾迹。马赫数影响主涡和次级涡涡核处流向涡量及"扭结"现象,但不影响主涡和次级涡涡核处流向速度和融合位置;雷诺数影响主涡和次级涡涡核处流向涡量、主涡涡核处无量纲流向速度以及"扭结"现象,但不影响次级涡涡核处流向速度和融合位置。在整个翼尖涡近场演化过程中,与第一阶段相比,第二阶段通过显著改变机翼上表面压力分布,诱导出强烈的下洗现象,主导影响着机翼的气动力,此外翼尖涡能抑制翼尖附近上表面流动分离,在一定程度上减轻其对气动力的不利影响。     

低雷诺数下椭圆翼型层流分离研究

通过采用计算流体力学方法对不同相对厚度的椭圆翼型在低雷诺数范围下进行了数值模拟，研究了椭圆翼型在低雷诺数下的层流分离现象及流场结构.结果表明:在低雷诺数下，薄椭圆翼型在小攻角时前缘出现层流分离泡是其具有高的升力系数及升阻比的原因，随厚度的增加，前缘层流分离泡逐渐消失，在后缘形成时均小泡.随着雷诺数升高，薄椭圆翼型时均分离泡都出现在前缘，但外形缩小，而在较低雷诺数下，薄椭圆翼型小攻角时没有发生转捩再附现象；同时层流分离泡的出现也对翼型后缘分离涡的尺寸和位置产生了重要影响.相对厚度和雷诺数通过影响椭圆翼型上表面层流分离泡的尺寸、位置以及后缘分离涡的形态结构，进而改变了气动特性.     

弹簧压并状态下悬架减振器绕流涡旋特性分析

为研究汽车行驶过程中减振器弹簧压并状态下翼子板内流场特性的变化,将该状态下的减振器简化为三维变截面圆柱模型,并建立变截面圆柱绕流三维流场模型,利用Transition SST四方程转捩模型模拟低、中、高3种车速对大、小圆柱绕流涡旋特性的影响.结果表明:绕流后尾涡的大小、形态、上升角均受圆柱直径、雷诺数及边界条件的影响,在变截面处验证“下洗”运动对N区边缘涡生长的直接作用及对L区涡旋分布的干扰作用;3种流速下适合绕流涡旋振动压电能量回收的最优夹角分别为±10°,±15°,±20°;在有界的高雷诺数流场下对变截面圆柱绕流涡旋重新分区,发现新的涡旋连接方式.