非结构动网格技术及其在超声速飞行器头罩分离模拟中的应用

针对激波诱导物体运动的二维动力学模型,比较了采用流动定常、结合非定常弹道方程的所谓"准定常"计算方法和非结构动网格计算方法,通过流场和运动特性分析,指出发展有相对运动的非定常流动模拟方法的重要性;另外采用修正弹簧近似和局部重构的非结构动网格技术,耦合求解Euler方程及6DOF弹道方程,数值模拟高度H=25km、飞行马赫数Ma=6.0、迎角α=0°状态下火箭抛罩的动力学过程.研究表明,稠密大气层内头罩受到的气动力比惯性力大几个数量级,导致头罩的运动特性完全由气动力控制.分离过程中,头罩和火箭头部之间的气动干扰在时间和空间上都表现出强烈的非线性特征,头罩姿态变化剧烈.本文还对分离过程中部分位置上的流场进行定常模拟,比较了定常和非定常的气动力,所得结论对于分析地面风洞可控轨迹系统(CTS)试验具有指导意义.     

雷诺数变化对翼型边界层发展及失速特性的影响

翼型失速及其边界层发展是飞行器设计中的基础科学问题,而雷诺数变化对其影响很大。针对后缘失速翼型,采用Menter k-ω SST模型及耦合扰动放大因子输运方程的转捩模型,进行雷诺数变化对层流-湍流转捩边界层特性和失速特性的影响分析。结果表明:雷诺数增大时,对于转捩边界层,当地涡量雷诺数增大,转捩前移且分离泡减小,流动能量耗散减小,翼型整体表面剪切效应增强,动能更充沛,流动自持能力增强,压力分布可以维持较长距离的梯度抵抗分离能力增强;因此雷诺数增大使翼型失速迎角提高、升力系数增加。     

阵列式表面电弧等离子体气动激励控制三角翼流动分离的实验

为探索多路阵列式微秒脉冲表面电弧放电(μs-SAD,Microsecond pulse surface arc discharge)对尖前缘小后掠三角翼流动分离的控制效果和作用机理,首先通过放电测试和纹影测试对多路阵列式μs-SAD的激励特性进行研究,揭示其对流场的作用原理,进一步将多路阵列式μs-SAD用于三角翼流动控制,开展了小后掠三角翼流动分离控制低速风洞实验,研究了来流速度、激励电压和激励频率等参数对控制效果的影响规律。结果表明:多路阵列式μs-SAD能够快速放热,单路瞬间放电能量可达68mJ,在流场局部可诱导产生冲击波;机翼前缘多路阵列式μs-SAD能有效改善三角翼大迎角气动特性,当来流速度为30m/s时,使最大升力系数提高27.2%,失速迎角推迟4°;来流速度增大到40m/s时,流动控制效果减弱,使最大升力系数提高15.5%;存在最佳激励频率使无量纲频率F+=1时,控制效果最好;激励电压存在阈值,其随来流速度的增加而增大,当激励电压超过阈值电压继续增大时,流动控制效果不再增强。     

基于目视间隔条件下的碰撞风险分析

为研究目视间隔条件下进近过程中前后两飞机的碰撞风险问题,通过Reich碰撞模型和碰撞风险模型,引入目视误差作为影响因素之一,建立基于目视间隔条件下的风险评估模型.从纵向、横向、垂向三个维度分析前后两飞机碰撞的概率,求得目视条件下进近两机碰撞的概率,以此作为风险分析的结果.计算结果表明,利用模型计算出的碰撞风险符合国际民航组织规定的航空器空中碰撞概率规定值.故而可证明,目视间隔和目视进近可作为一种安全的空中交通管制运行模式缩小安全间距,提高空域利用率,增加飞行流量.     

不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

基于充气前缘技术的旋翼翼型动态失速抑制

动态失速的发生会在直升机旋翼桨叶和桨毂上产生高的交变扭转振动载荷,并限制直升机高速重载状态下的使用包线。本文利用计算流体力学(CFD)方法对基于充气前缘(ILE)技术的SC1095旋翼翼型动态失速抑制进行研究,分析了ILE抑制动态失速的控制机理,获得了ILE结构布置和充放气方式对动态失速的影响规律。研究表明:ILE可以有效抑制动态失速的发生;ILE最大膨胀程度越大,其抑制动态失速的效果越好,但膨胀程度过大后抑制效果开始减弱;ILE在翼型上仰至最大迎角时恰好达到最大膨胀状态,其对动态失速的抑制效果最好;ILE保持最大膨胀状态的时间长短对抑制效果影响不大;在翼型上仰至不同迎角时开始对ILE充气会对动态失速抑制有较大影响;ILE整流段与翼型连接位置对动态失速抑制有很大影响,整流段越长,抑制效果越好。     

大涵道比风扇/增压级进气畸变数值模拟研究

为探究进口总压畸变下大涵道比风扇/增压级内部流场的主要特征,基于三维彻体力模型的思想,开发了一套能够实现风扇/增压级内外涵联算功能的三维数值计算程序。利用该程序模拟了某大涵道比风扇/增压级在周向总压畸变进气下的三维流场。计算结果表明：大涵道比风扇单转子不同叶高处的畸变传递特征差异较大,转子出口总压畸变强度由叶根到叶尖逐渐降低,在叶尖处衰减为最小值1.5%;在转子出口相应诱导出的总温畸变强度由叶根到叶尖逐渐升高,在叶尖处达到最大值1.4%;进口周向总压畸变导致风扇转子总压比下降0.5%,而风扇转子出口形成的总压总温复合畸变导致增压级总压比下降2%;总压畸变在增压级中呈逐级衰减趋势,而高温畸变区的周向范围在逐级增加。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构,层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟,研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型,基于结构化拼接网格,对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性,雷诺数对翼型气动特性具有显著影响。随雷诺数降低,时均分离泡外形增大、位置后移,平均阻力系数增大,特别是在较低雷诺数下,翼型升阻力系数随时间出现振荡现象。进一步研究表明,造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低,流动黏性增大,导致分离剪切层速度梯度减小,流动发生转捩及再附位置后移,直至翼型表面不再发生转捩和再附。     

风洞T型冲击三通管道流场特性数值模拟

在24 m跨声速风洞进气管路T型冲击三通连接处和隔板位置,每隔一定时间会出现裂纹,影响了风洞的安全稳定运行。为了解裂纹出现的原因,同时为改进设计提供依据,采用CFD(computational fluid dynamics)对风洞现有三通和优化方案进行了对比。控制方程为三维黏性不可压缩Navier-Stokes方程。结果表明：无隔板时,三通内的流动为最常见的类型,支管内存在3个分离区,在转向过程中形成第1分离区,即马蹄涡;随后是第2分离区,一对反向旋转的旋涡,即Dean涡,及三通顶部壁面形成第3分离区。流场沿y=0 mm和z=0 mm平面基本呈对称分布。有隔板时,流场的左右和上下结构均不对称;在隔板和外侧壁面间的角点形成范围较小的驻涡,在支管内形成不稳定的螺旋状分离涡,致使气流振荡,从而使得三通连接和隔板处管壁出现裂纹。根据上述流态设计了管路优化方案和整流装置,优化后能有效减小或消除分离;其中去掉隔板最简单易行,可以解决裂纹的问题。     

介质阻挡放电等离子体对翼型流动分离控制的实验研究

在低速开口风洞中进行了等离子体激励器对NACA0015翼型流动分离控制的实验研究。采用PIV技术,对翼型绕流流场进行了测量,显示了施加等离子体激励后流场的变化。通过五分量天平对升力和阻力的测量,研究了激励电压和激励频率对翼型流动分离控制的规律。研究表明,低风速下在翼型前缘施加等离子体激励,能够有效地控制翼型流动分离,在来流为20m／s时,最大升力系数增加11％,失速迎角增加6°;在给定的流动状态下,激励电压和激励频率存在一个阈值,不同迎角下该阈值不同,迎角越大,分离越严重,对激励强度的要求也越高。