用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

后掠角对后掠机翼边界层稳定性及转捩的影响

后掠机翼边界层流动稳定性及转捩对翼型的设计及优化有着重要的参考价值,而机翼后掠角是引起后掠机翼边界层横流失稳的关键参数之一.以NACA0012翼型为研究对象,通过求解三维可压缩Navier-Stokes方程计算了展向无限长后掠机翼的基本流场;通过求解Orr-Sommerfeld方程得到了扰动Tollmien-Schishting波演化的中性曲线及幅值曲线,研究了后掠角对后掠机翼边界层流动稳定性的影响;最后采用eN方法进行了转捩预测.研究发现,随后掠角的增大,横流强度和扰动幅值放大指数n均先增加后减小,且后掠角在40°~50°之间横流强度达到最大值.当后掠角在50°左右时,用转捩预测eN方法计算的幅值增长指数N值最大,导致转捩发生所需的初始扰动幅值最小,转捩最易发生.      

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

磁尾静电离子束流-密度漂移不稳定性的数值研究

本文研究了磁尾等离子体片边界层宽频带静电噪声的产生机制.模型等离子体由暖的背景电子、暖的向地球方向的离子束流和较冷的向尾方向的离子束流组成.结果表明,静电离子束流-密度漂移不稳定性可以在比较宽的频率范围内激发,在低频区大传播角方向上增长率最大,在高频区小传播角方向上增长率也比较大。最大增长率的方向取决于离子束流和密度漂移的速度比值.这些结果与磁尾观测到的宽频带静电噪声特征符合一致.      

NFFD控制点分布对气动外形优化的影响

基于非均匀有理B样条的自由型面变形(NFFD)技术具有对变形对象普适性和控制点影响区域局部性的特点,广泛应用于气动外形优化。本文通过扩展控制体和合理布置外侧控制点,实现了NFFD技术在参数化曲面并改变曲面形状时,同步变形控制体内的表面网格和空间网格,并从理论上保证了控制体边界内外的网格协调。基于离散伴随方法求取梯度,分别采用拟牛顿(QN)法和序列二次规划(SQP)优化方法,通过从初始翼型NACA0012到标准翼型EH1590的反设计,研究了NFFD控制点数量和分布对设计结果的影响。在某飞翼标模单点全机升阻比优化应用中,改进控制点分布后获得了更高的升阻比,收敛速度显著提高。     

不同厚度翼型动态失速涡运动数值研究

在低马赫数下, 对三种不同相对厚度的NACA系列基本厚度翼型在俯仰振荡运动中的动态失速现象进行了数值研究.数值模拟时采用双时间法和LU-SGS隐式解法相结合的模式求解了非惯性坐标系下的非定常纳维-司托克斯(NS)方程组, 空间离散采用Roe格式,并结合刚性动网格生成技术.采用全湍流计算,通过引入BL(Baldwin-Lomax)模型计入湍流影响.计算出的气动力迟滞曲线与实验结果变化趋势符合较好,表明了数值方法的有效性,同时通过比较和分析不同厚度翼型在轻失速和深失速下的流场结构,发现绕翼型的失速旋涡产生和发展规律是明显不同的.      

防冰表面的对流换热计算分析

利用动量和热边界层积分分析的方法对采用电热防冰的二维机翼对流换热系数进行了数值计算,比较了等温突捩以及变温光滑两种不同的边界层模型.层流向湍流转捩的起始点位置和过渡区域的长度由经验公式判断,并通过内外热流耦合迭代求解得到了防冰表面的平衡温度,与文献试验值进行了比较.所得结果表明:加入过渡区模型可较好地预测对流换热系数及防冰表面温度.同时,结果还表明环境压力变化对换热系数的的影响较大.      

翼型波浪水面巡航地面效应数值模拟

数值求解非定常不可压缩流动的雷诺平均N-S(Navier-Stokes)方程和标准k-ε湍流模型,用VOR(Volume Of Fluid)方法模拟波浪水面,模拟了NACA2410翼型在波浪水面上方飞行的流场.研究了余弦波浪水面的生成方法,选取合适的计算网格和时间步长,避免了在生成规则余弦波浪水面时由于数值耗散使波面形状衰减.比较了固壁波浪与水面波浪计算结果的差异,发现固壁波浪的结果更接近余弦曲线分布.研究了波浪等级对翼型气动性能的影响,分析计算结果发现:在规则的余弦水面波浪上方飞行,翼型气动力呈现周期性,给出了一个周期内气动力的变化过程,以及波长和波幅对气动力平均值和波动幅度的影响规律.      

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

低噪声风力机翼型设计方法及实验分析

为了研究风力机翼型的噪声特性,基于翼型泛函集成理论与翼型噪声计算模型,建立了低噪声翼型优化设计数学模型,提出在设计攻角情况下升阻比与噪声比值最大为目标函数,对优化后的新翼型CQU-DTU-B18翼型与NACA-64-618翼型在相同的风洞实验及风速条件下进行了噪声对比分析.研究表明,理论噪声计算模型虽然与实验数据有一定的偏差,但是翼型的升压级随频率的变化趋势是一致的,表明了翼型噪声计算模型的准确性;相比NACA-64-618翼型,CQU-DTU-B18翼型具有更低的噪声特性,从而验证了该设计方法的可行性.对于如何设计低噪声翼型及怎样降低翼型噪声具有重要的指导作用.      

沟槽面湍流减阻研究综述

对近20年来沟槽面湍流边界层特性、湍流拟序结构、湍流减阻及其机理的研究进展进行了综述.内容涉及沟槽面平板、旋成体、机翼等在压、跨、超音速流动情况下的实验研究成果;压力梯度、攻角、侧滑角等的影响;湍流猝发特性、紊动特性、近壁区带条结构的特征及减阻机理等方面的工作.为更有效地减少表面摩阻,必须深入开展对沟槽面湍流边界层特性、湍流拟序结构及湍流减阻机理等方面的研究.      

低雷诺数下翼型前缘流动分离机制的研究

采用高精度有限差分格式,对来流雷诺数为1.0×104,攻角为3°的二维翼型流动进行了直接数值模拟,研究了低雷诺数下翼型前缘流动的分离机制,描述了分离涡系的相互作用规律.计算结果表明:前缘椭圆弧靠近叶身位置存在吸力峰,流动在吸力峰内强逆压梯度的作用下发生分离;翼型上表面形成了包含驻留涡、脱落涡和二次涡的涡系结构,其尺度随时间不断变化,具有强烈的非定常性;表面压力分布曲线可以较好的描述翼型边界层流动.      

76°/40°双三角翼前缘涡破裂及其控制实验研究

流动显示结果表明,喷流能有效地推迟双三角翼前缘涡的破裂,且随着攻角的增大,前缘涡破裂位置逐渐推后,喷流极大地改善了大攻角情况下前缘涡的非对称破裂特性,能有效地克服可能出现的机翼的"摇滚"现象.另外,后缘喷流可以减弱乃至消除前缘涡混掺现象的发生,进而有利于飞行器的操纵.      

后掠翼身干扰区流动特性及改善措施研究

利用流动显示及表面压力测量方法研究了后掠翼身干扰区的流动特性,并研究了用小边条等措施改善干扰区的流动特性的效果.结果表明,随着不同机翼后掠角、不同迎角及不同Re数对干扰区流动特性的影响,流态可以从一涡系变成多涡系,由定常变成非定常,而且在一定的Re数以后涡系会紊流化;翼身干扰区上游的的逆压梯度是导致边界层分离的物理原因,利用面积很小的边条可以降低干扰区局部的逆压梯度,可以导致干扰区的旋涡很弱,甚至不出现,这是很有实际意义的.     

76°/40°双三角翼前缘涡破裂及其控制实验研究

流动显示结果表明,喷流能有效地推迟双三角翼前缘涡的破裂,且随着攻角的增大,前缘涡破裂位置逐渐推后,喷流极大地改善了大攻角情况下前缘涡的非对称破裂特性,能有效地克服可能出现的机翼的"摇滚"现象.另外,后缘喷流可以减弱乃至消除前缘涡混掺现象的发生,进而有利于飞行器的操纵.     

一种热噪声标准源等效输出噪声温度的计算方法

资用噪声功率即等效输出噪声温度是热噪声标准源的一个关键技术指标,其是匹配负载终端和有耗传输线噪声贡献的合成。为减小等效输出噪声温度的计算复杂度,采用了多项式拟合后分段叠代计算的方法,并分析了有限温度测量点的情况下,不同分段数对计算结果的影响。     

伞衣透气性对翼伞气动特性的影响

为了提高翼伞的飞行性能,需要研究伞衣织物透气性对翼伞气动特性的影响。使用不可压雷诺时均Navier-Stokes(RANS)方程模拟伞衣外部流场,建立了包含附加动量源项的多孔介质域控制方程模拟伞衣,对2种透气性材料模型和无透气性影响传统模型的气动特性和流场分布进行了二维和三维定常数值模拟。数值结果表明,求解多孔介质域控制方程可以得到较准确的伞衣透气速度,伞衣表面的湍流度急剧增加;使用较大透气量材料制作伞衣时,升力系数大幅下降,阻力系数大幅上升,同时会造成内腔泄压影响翼伞的外形保持;使用微透气量材料制作伞衣时,升力系数在小迎角时小于不透气模型,在大迎角时大于不透气模型,较小的透气速度能在大迎角时延缓边界层分离。     

亚临界雷诺数圆柱绕流远场气动噪声实验研究

针对亚临界雷诺数下圆柱绕流的气动噪声问题,在北京航空航天大学D5气动声学风洞中进行了圆柱绕流远场气动噪声的实验研究。在D5气动声学风洞实验段不同来流条件下,测量消声室中距离圆柱不同位置的自由场传声器的远场声压信号,经过快速傅里叶变换(FFT)获得亚临界雷诺数下圆柱绕流气动噪声随来流速度、接收点距离的变化规律。实验结果表明:圆柱绕流的远场气动噪声在垂直气流方向最大,随着偏离垂直方向角度的增大噪声逐渐减小,属于典型的偶极子噪声源辐射特性;远场气动噪声近似与来流速度的6次方成正比,与接收点到圆柱中心距离的2次方成反比,这说明偶极子类型的噪声源是圆柱绕流的主要噪声源。      

合成射流控制NACA0015翼型大攻角流动分离

为了研究合成射流激励器处于NACA0015翼型回流区时对其分离流动的控制,采用商用计算流体力学软件Fluent 6.1求解Reynolds平均Navier-Stokes方程,通过对翼型气动力特性、脱落漩涡结构以及射流孔口附近流动结构的分析,揭示了合成射流处于分离区时对边界层控制的机理.结果表明:当合成射流孔口处于回流区时仍可有效推迟翼面边界层分离点,缩小回流区范围,从而有效提高翼型的升力.当射流方向垂直于壁面,无量纲频率以及吹气速度比都等于1时,翼型平均升力系数提高40%左右.