一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

槽道内涡波流场展向涡的分布特征

通过2DPIV测量了槽道内涡波流场的瞬态速度矢量场,利用雷诺分解和旋转强度的λci准则进行了展向涡的识别;并从识别的逆时针和顺时针展向涡的尺度属性、力学属性和运动属性方面分析了展向涡的分布特征。结果表明,流场中的顺时针展向涡的数量大于逆时针展向涡的数量;在法向位置0.05~0.45范围内,逆时针展向涡的旋转强度呈抛物线形状变化,顺时针展向涡整体呈下降趋势;逆时针展向涡的平均直径和椭圆平均长轴长度随着法向位置的增加而缓慢增加,而顺时针展向涡趋于下降;逆时针展向涡的椭圆长轴倾角最大值为67.28°,顺时针展向涡的等效椭圆的离心率呈下降趋势,似圆性较好;逆时针展向涡对全局平均的〈vωz〉分量在法向位置0.15处出现最大值61%;展向涡对流场雷诺应力的影响较小;逆时针展向涡的数量密度在0~0.08法向位置范围内呈急剧上升趋势,并达到最大值0.36,顺时针展向涡数量密度整体呈倒U形分布;不同旋向的展向涡占总展向涡数量密度的比值都大于0.3;展向涡的流向对流速度整体上小于流场的平均速度,且在槽道中心区域差值较大;展向涡的法向对流速度随着法向距离增加,从负值上升到正值,在中心主流两侧分别发生喷射和扫掠2种现象。     

旋翼前后行侧桨-涡干扰噪声辐射特性

基于Beddoes修正尾迹模型,推导了一个新的旋翼桨-涡干扰(Blade-vortex interaction,BVI)噪声源方向性预测和BVI噪声辐射特性分析的模型,为更好地研究桨-涡干扰噪声的参数影响,提供了一套计算方法.该方法由移动的桨尖涡和旋转的桨叶相遇来确定桨-涡干扰的位置,由相交点的桨叶和涡线相对位置得出BVI夹角、轨迹马赫数等重要参数的计算公式.应用所建立的方法,首先针对有相关结果可供对比的BVI源发生位置进行了算例计算,以验证方法的有效性.然后,着重对直升机前飞状态的旋翼BVI噪声的发生进行了预测,深入分析了在旋翼前行侧和后行侧BVI噪声源的方向性和BVI噪声的辐射规律,并研究了前进比等参数对BVI噪声传播方向性的影响.在此基础上,得出了一些新的结论.     

基于投影法的双后掠乘波体气动性能

双后掠布局能有效改善乘波体低速时的气动性能不足。为了获得双后掠乘波体,目前常采用的是定前缘型线的吻切锥乘波体设计方法,但该设计方法存在设计过程复杂,激波出口型线与理论不一致等问题。而采用直接投影获得双后掠乘波体的设计方法可以解决上述问题。为了系统研究基于投影法的双后掠乘波体的气动性能,使用CFD方法分析了采用该方法生成的双后掠乘波体在高超声速与低速时的气动性能。结果表明,该方法获得的乘波体在高超声速下的气动性能与定前缘型线的双后掠乘波体相当。且此方法仍保留了高超声速下"波效应"引起大攻角非线性升力、低速下"涡效应"有效提高升阻比等双后掠乘波体的优良气动特征,为基于投影法的双后掠乘波体的工程应用提供了指导。     

轴向旋流器单通道甲烷/空气回火特征大涡模拟研究

为了研究轴向旋流器叶片单通道间的回火特征,采用大涡模拟方法对旋流器叶片吸力面侧和压力面侧点火后火焰面的发展过程进行了比较研究。结果表明:纯流动下,吸力面侧有回流区,压力面侧无回流区;吸力面侧点火后,火焰面随吸力面侧的回流区向上游发展;压力面侧点火后,压力面侧不会回火,但是会引燃吸力面侧,进而回到吸力面侧回火模式;火焰面传递的关键位置为内壁面或叶片尾缘。     

近失速状态下压气机转子叶尖旋涡流动研究

在低速大尺寸单级压气机实验台上,利用SPIV技术测量了近失速状态下压气机转子尖区多个截面的三维瞬态速度场。基于瞬态场测量结果,详细阐述了近失速状态下压气机转子叶尖泄漏涡的演化过程和角区旋涡的形成过程。测量结果表明角区旋涡是一种总体意义上的旋涡,其涡核是由多个涡团组成的,形成角区旋涡的一个关键机制是压气机的旋转运动对源于近吸力面的正负涡量的涡团具有选择性。     

三角翼受迫滚转运动气动特性的数值模拟研究

基于动网格的技术,利用有限体积法对N—S方程进行数值离散,选用湍流模型,对65°大后掠三角翼的受迫滚转运动进行了数值模拟,在此基础上．对滚转运动的气动力特性和流场结构进行了分析。计算结果表明：滚转运动时．三角翼上表面涡分布的非对称性产生横侧方向的偏航力矩和滚转力矩,偏航力矩、滚转力矩系数滞回曲线呈现典型的“单环”形态,俯仰力矩、法向力滞回曲线呈现典型的“单8”形态。     

低雷诺数下跨声速转子流动失稳及周向槽处理机匣扩稳

采用数值方法模拟了低雷诺数条件下NASA Rotor 37跨声速压气机转子内部流场。结果表明,附面层径向涡是该压气机转子流动失稳的一个很重要的原因。根据该压气机转子低雷诺数条件下流动失稳的特点,研究了周向槽处理机匣结构对其性能的影响。结果表明,引入处理机匣后,附面层径向涡得到一定程度的抑制,由附面层径向涡所引发的叶顶阻塞区有所减小,提高了压气机转子的失速裕度。     

非定常集中涡破裂形态及其对扰动的影响

本文中利用三角翼前缘涡襟翼的强制振动来推迟涡破裂。流场显示实验结果表明,在有强制扰动时,集中涡显示出两种不同形式的破裂,可分成六种不同的破裂形态。非定常集中涡的一个破裂过程涉及几种破裂形态之间的转换。涡襟翼振动产生的非定常效应对涡破裂的影响与三角翼后掠角有密切的关系,随后掠角增大非定常效应的影响变小。就实验迎角范围而言(α<35°),对50°后掠角三角翼,涡襟翼振动可显著推迟涡破裂;但对70°后掠角三角翼,振动却能促进涡的破裂。涡襟翼振动改变了集中涡的性态,使集中涡趋向于和前缘平行。振动对涡破裂的作用部分地与这种效应有关。     

平行四边形截面的旁路式双喉道气动矢量喷管数值研究

将异型出口设计与气动推力矢量喷管结合,提出了平行四边形截面的旁路式双喉道气动矢量喷管(bypass dual throat nozzle,BDTN),并与基准矩形截面的BDTN进行流场结构及气动性能的对比研究。三维数值计算表明:平行四边形构型与矩形构型具有相同的气动性能变化规律;相较于基准矩形BDTN构型,由于壁面倾斜导致流场结构变化,平行四边形构型在矢量状态下的俯仰矢量角及推力系数有所降低,但对非矢量状态下的推力系数和流量系数影响不大;壁面倾角是性能变化的重要因素,相同落压比时壁面倾角越小,矢量角越小,壁面倾角不小于60°时喷管的稳定矢量角均可达10°以上,最大矢量角均可达15°以上;平行四边形出口增强了尾喷流与环境气流的掺混,出口射流中心线速度衰减大大加快,有利于提高红外隐身特性。