三维机翼大迎角低速绕流及其涡控制的数值模拟

通过数值求解三维非定常可压缩雷诺平均Navier-Stokes方程，研究了三维机翼大迎角低速绕流及其带表面周期性吹吸气的流动控制，并探讨了非定常质量引射作为外激发手段对机翼气动力特性的影响，计算表明，在机翼前缘附近的小区域内施加周期性吹吸气进行流动控制，这类外激发手段对三维机翼的大迎角分离流动的控制是可行的，通过选择适当的控制参数，尤其是合适的周期性激振频率，可以有效地改善气动力特性。     

绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

无限翼展后掠翼大迎角绕流和涡控制的数值模拟

通过数值求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,研究了无限翼展直后掠翼在不同后掠角下的大迎角粘性分离流,探讨了后掠角对流场结构和升阻力特性的影响。基于对流动特性的机理分析,文中进一步数值模拟了无限展长后掠惭带表面吹吸气的大迎角绕流,研究了以非定常质量引射作为外激发手段对后掠翼前缘涡形成、发展和脱泻、以及对提高机翼升力特性的影响。     

多孔介质空心圆柱绕流的摄动解

本文主要讨论了绕多孔介质空心圆柱的有势流动问题。运用摄动方法,给出了速度势的四阶近似,然后得到了圆柱表面的无量纲的压力系数和阻力系数。结果发现在势流的情况下多孔介质圆柱受到附加阻力作用,并且在圆柱的喷射或吸收流体的特征速度与来流特征速度的比值是一个小量的时候,阻力系数与特征速度比成简单的线性关系。     

周期性变形圆柱绕流特征的数值分析

应用有限体积法求解不可压缩二维N-S方程，对在均匀来流中圆柱-椭圆柱周期性变形运动进行了数值模拟。对于两种典型状态，即雷诺数为40和100，本文描述了圆柱变形运动所引发的非定常流动随时间变化的过程，揭示了变形过程中柱体附近流场的变化主要取决于由壁面压力梯度和壁面的切向运动加速度所产生的壁面涡量流，分析了不同变形频率和不同变形量对升力和阻力、旋涡形成及其运动的影响，并发现了频率锁定现象。     

低雷诺数翼型绕流的电磁控制

对低雷诺数不可压翼型绕流在Lorentz力(电磁力)作用下的流场变化情况进行了实验探讨.发现当施加的Lorentz力与翼型表面流向相同时,其表面流体的分离被抑制,从而使尾流涡消失,且翼型阻力减小而升力增加;当Lorentz力与流向相反时,则其受阻力变大而升力降低.     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

垂直交叉双圆柱绕流数值模拟及涡结构分析

用分块耦合方法数值模拟了垂直交叉双圆柱绕流,两圆柱间距为圆柱直径的5倍,所取雷诺数为200.在交叉结构的中心流场比较复杂,用λ2方法准确地描述了其中的三维涡结构,发现上游圆柱规则的卡门涡街被下游圆柱截断成复杂的缠绕结构.与串列双圆柱相比,垂直交叉双圆柱绕流上下圆柱的升阻力特性、Strouhal数等都呈现出许多不同的特点,下游圆柱升力振幅小于上游圆柱,阻力平均值大于上游圆柱.本文旨在模拟出这种特定间距下的流场涡结构以及圆柱升阻力特性.     

激光等离子体点源减阻技术中入射能量对气动阻力的影响

数值模拟了来流马赫数为5时的圆球绕流冷流场,以及来流与头部驻点附近注入激光能量而产生的等离子体相互耦合形成的流场的演化过程。计算结果表明:施加局部能量点源能够改变原有的弓形激波结构,使其变为斜激波,从而减小气动阻力;气动阻力随入射激光能量的增大而逐渐减小,当入射激光能量为1.1 J时,气动阻力减小的百分比高达40%。     

表面吸气方法控制分离的数值模拟

本文使用数值模拟方法研究了吸气效应对分离的控制及增升作用。数值模拟采用的出发方程为可压流的Ｎ－Ｓ方程，格式为Ｂｅａｍ－ｗａｒｍｉｎｇ格式的改进型。主要研究结果表明：（１）在一定的吸气压强Ｃ＿（ｐｓ）范围内，绕圆柱流动的卡门涡街现象基本消失，流动从典型的非定常流趋于定常，圆柱的升阻比明显提高；（２）对于圆柱，吸气压强系数Ｃ＿（ｐｓ）＝－６时，综合效果比较好；（３）吸气孔位置开在圆柱上表面逆时针方向６０°时，圆柱的升阻比较高；（４）吸气可有效地提高翼型的升力     

基于被动控制的D型钝体减阻大涡模拟研究

钝体减阻在航空航天等多个领域内具有重要的应用潜力,也是重要的基础研究热点问题。为了精确捕捉流动控制的流动细节,并发展效果优良的流动控制减阻方案,论文围绕抽象出的D型钝体,采用大涡模拟方法,开展了被动控制减阻高精度数值模拟研究。首先基于前期研究成果,对D型钝体尾迹区剪切层附近放置一个光滑小圆柱的被动减阻方法开展了数值模拟,发现总阻力减小17.7%,与试验结果吻合很好,同时数值预测的速度场分布也与试验结果吻合良好。在此基础上,进一步提出了采用齿槽型表面结构的小圆柱对D型钝体尾迹区进行扰动,并开展了数值验证,发现总阻力减小21.4%,优于前期的减阻方法。最后研究了在三种雷诺数工况下两种小圆柱扰动情况下的减阻效果,均表现出良好的减阻效果,两种小圆柱扰动下总阻力最大降幅分别为19.6%和23.1%,同时基于大涡模拟计算结果对减阻流动机理进行了探讨。上述研究结果表明,通过进一步优化流场结构,可以得到更优的流动减阻方案。     

采用上游扰动体被动控制减阻技术的实验研究

通过风洞测压和水洞流动显示实验,研究了不同直径的扰动体尾迹下圆柱的流动结构和气动特性。结果表明当扰动体和圆柱串列布置时,流动呈现出两种流动模式:空穴流动和尾流撞击模式。当圆柱处于空穴流动模式时,其阻力可以大幅度降低。当扰动体直径为圆柱直径1/2时,圆柱阻力最大降低98%。综合绕流结构和圆柱表面压力分布的变化,指出阻力减小的原因主要是扰动体对圆柱前体的保护作用以及对后体流动的干扰作用,其中后者占主导地位。随着扰动体直径的增加,扰动体对圆柱前体的减阻效果增加。     

趋于临界马赫数的圆柱跨声速绕流特性分析

通过深化认识趋于临界马赫数Ma_cr的圆柱跨声速绕流特性,明确新型飞行器增升减阻设计的空气动力学理论依据。采用大涡模拟方法数值研究了来流马赫数Ma_∞为0.75和0.85、雷诺数Re为2×10⁵的圆柱跨声速绕流。结果表明:当Ma_∞趋于临界马赫数(Ma_cr≈0.9)时,圆柱的阻力下降且升力系数振荡被抑制;通过力的分解,得知圆柱的阻力减小来自旋涡力的影响,而非可压缩性;圆柱的阻力减小与其背压上升有关;剪切层初始阶段的对流马赫数Ma_c随Ma_∞的增加而增大,而增长率相反,这使得剪切层更为稳定、柱体背压更高。此外,由于Ma_∞=0.85时边界层分离点处的激波和尾迹处的激波向下游推移,使得近尾迹处的湍流脉动减弱,进而导致柱体的表面压力振荡和升力系数振荡被抑制。     

非定常尾迹输运对动叶气膜冷却流场影响

对不同尾迹宽度下涡轮动叶的气膜冷却流场进行非定常数值研究,研究了尾迹输运对动叶气膜冷却流场的影响机理.计算结果表明:非定常尾迹的总压损失和速度损失而形成低速区,当非定常尾迹输运到冷却孔附近时,尾迹对气膜层的干扰占主导地位,靠近气膜孔的叶片区域冷却效率降低明显.非定常尾迹的影响使得冷却气流被带到压力面叶片下游更广阔的区域,导致气膜冷却沿叶展方向的覆盖范围更广,从而使得压力面气膜冷却效率在尾迹通过时更高一些.当尾迹宽度增大时,对气膜冷却效率的影响程度增大.尾迹对压力面的影响比吸力面更大一些.      

静叶时序对压气机叶片附面层流动影响的数值研究

采用数值方法对两级低速压气机中径处的非定常流场进行模拟,针对压气机第2排静叶两个典型周向位置对动静叶干扰下的叶片附面层流动进行研究.建立尾迹与附面层干扰分析模型,结合叶片壁面摩擦力和近壁面附面层湍动能,详细分析了尾迹和势流干扰下静叶时序改变对叶片附面层流动产生的影响.对第2排静叶附面层的研究结果表明:静叶时序改变了尾迹在其叶排中的输运特征,能够降低壁面摩擦力和近壁面湍动能及其非定常最大波动幅值,影响吸力面附面层内动叶尾迹后沉寂区的宽度.在非定常条件下,尾迹能够诱导静叶层流附面层在尾迹干扰的局部范围内转捩发展为湍流状态,同时高湍流度尾迹的干扰具有抑制逆压梯度下附面层分离的作用,并能够延长层流区的范围.     

均匀来流中旋转振荡圆柱绕流的数值研究

通过求解二维非定常不可压缩 Ｎ-Ｓ方程,研究了均匀来流中旋转振荡圆柱的流动特性,系统地改变流动的振荡频率、振幅,分析了圆柱的受力特性及尾迹结构,发现了与横向、流向振荡的圆柱相似,旋转振荡也是流动控制的一个有效方法,有其自身的流动规律及特点。     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

双极性等离子体激励器圆柱绕流控制实验研究

在低速风洞中利用多级双极性等离子体激励器控制圆柱绕流的流动分离。实验风速U∞=10m/s,基于圆柱直径的雷诺数Re=2.8×10^4,在实验中将两组三级双极性等离子体激励器布置在圆柱模型肩部,利用粒子图像测速技术测量圆柱的尾流场。实验结果表明,采用定常和非定常激励均能抑制圆柱尾迹区,等离子体激励强度是影响激励器对圆柱绕流控制能力的重要因素;非定常脉冲激励耗电少,对流动控制能力强,效率明显高于定常激励,脉冲激励频率影响等离子体激励器对流动的控制能力。在实验风速为10m/s时,脉冲激励频率与圆柱涡脱落频率一致,流动控制效果较好。     

垂直交叉双圆柱绕流强迫振动数值模拟

采用虚拟体法数值模拟了平面单圆柱绕流强迫振动,以及空间垂直交叉双圆柱绕流下游圆柱强迫振动流场。验证了单圆柱强迫振动中的锁定状态以及相位突变现象,从而证实了该数值方法模拟振动流场的可靠性。研究了在雷诺数Re=150、间距为5倍圆柱直径、下游圆柱按正弦曲线振动时,对尾流场的影响。下游圆柱两端尾流场在振动的作用下,涡的横向间距增大,而中心尾流场由于受上游圆柱尾流的影响而保持原先的状态。下游圆柱在锁定区的振动使尾流场变得稳定,表现在流向二次涡结构的减少,以及尾流场的速度场频谱趋向单一化。