基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

柔性扑翼气动结构耦合特性数值研究

微型扑翼体积小、重量轻,其柔性变形对气动特性有显著的影响。通过求解雷诺平均N-S方程(ReynoldsAveraged Navier-Stokes,RANS)和结构动力学方程,对微型柔性扑翼飞行器的气动结构耦合特性进行了数值模拟研究。针对微型扑翼的大幅运动,发展了适用于扑翼的气动结构耦合数值计算方法,研究了微型扑翼的气动结构耦合特性。通过求解雷诺平均Navier-Stokes(RANS)方程得到微型扑翼的非定常气动特性;利用哈密顿原理(Hamilton Principle)推导了扑翼的结构动力学方程,采用结构有限元方法对该动力学方程进行离散并求解,得到扑翼的动态结构特性;采用松耦合方法进行迭代。计算结果与风洞实验结果相比吻合良好,验证了所发展方法的有效性。在此基础上研究了惯性力和关键运动参数对柔性扑翼气动及结构特性的影响规律,有助于比较详细、全面地了解微型扑翼的气动机理,为柔性扑翼的设计提供了参考依据。     

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

波浪形非均匀间隙封严结构影响涡轮性能的数值模拟

利用数值模拟的方法研究了波浪形非均匀间隙封严结构和均匀轴向间隙封严结构下轮缘封严气流对涡轮性能的影响.研究表明:燃气入侵与出流结构受到静盘、动盘及主流切向速度的影响,以低于动盘转速同向旋转,并改变了转子的进气条件,增强了压力面马蹄涡强度,因此对转子出口流场造成很大影响.封严气流与上游导叶尾迹的相互作用引起转子通道内熵增,造成涡轮效率的下降.与均匀轴向间隙封严结构相比,波浪形非均匀间隙封严结构使大的入侵与出流结构破碎为小的结构,对涡轮性能的负面影响减小,涡轮效率提高了0.9%.结果证明了波浪形非均匀间隙封严结构在具有较好的封严效果的同时提高了涡轮性能.      

平衡风速下十字形伞弹系统姿态摆动分析

以一定高度和速度飞行的母弹从其伞弹舱中抛撒伞弹系统群,各个伞弹系统的姿态摆动情况将影响其落点及落角,进而影响其作战效能。采用立式风洞试验和非定常数值仿真(CFD),分析平衡风速下十字形伞弹系统的姿态摆动情况。通过立式风洞试验,可确定柔性伞的外形、伞弹系统的阻力系数和摆动频率,并创新"拉拽式"试验模型的约束方法;采用基于三维N-S方程的CFX软件,进行伞弹系统非定常数值仿真,分析其摆动机理。结果表明:数值仿真得到的阻力系数、法向力系数、侧向力系数及其摆动频率,均与风洞试验结果相吻合,即数值仿真结果能够反映风洞试验中伞弹系统的摆动情况;伞弹系统的流场极不稳定,伞衣内部有一对方向相反、强度交替变化的旋涡,伞衣外部存在不稳定分离流动,二者相互关联,使得气动参数呈周期性波动,导致伞弹系统的姿态摆动。     

用于压气机流动计算的3种模型比较

为了更好地对压气机流动进行模拟,在课题组自行开发的结构化有限体积解算器上实现了用于压气机流场计算的混合平面法、谐波平衡法及相滞后法.以NASA Stage 35为例,对3种方法的计算结果进行了比较分析.结果表明:相滞后法的计算精度最高,混合平面法的计算精度最低;相滞后法与半环的双时间推进法结果相近,计算速度提高了20倍;相比混合平面法,谐波平衡法能准确地模拟动静叶间的非定常干涉及进出口参数变化;在谐波阶数达到5阶后,谐波平衡法计算得到的结果不随阶数变化,且与相滞后法的结果基本吻合;混合平面法的计算效率远高于另外两种方法,相滞后法与谐波平衡法在谐波阶数为5阶时的计算效率相当.      

跨声压气机内确定应力耦合问题研究

基于三维定常Denton程序发展了二维定常、二维非定常、三维非定常计算程序,以二维跨声压气机叶栅、三维跨声单级压气机近失速点为算例,进行确定应力耦合问题的研究.结果表明:在确定应力应用到时均方程的耦合处理中存在一些因素导致时均方程的解和非定常解的时均值有一定偏差,它们分别是空间差分离散引入的误差、黏性相关项引入的误差、掺混面处理引入的误差.在多排耦合问题中掺混面的处理方式至关重要,对耦合计算结果有较大的影响,简单的掺混面处理并不能获得理想的效果.      

轴向间距对涡轮时序效应影响的数值研究

时序效应可以改进多级叶轮机械的效率,为了探讨轴向间距对时序效应的影响,并应用到叶轮机械气动设计中,对1.5级静/动/静布局的涡轮叶片进行了3维非定常数值模拟。结果表明:保持轴向间距不改变,转子居中时涡轮效率整体最高,时均效率最大值为0.866,转子前移时次之,时均效率最大值为0.861,转子后移时效率最低,时均效率最大值仅为0.859,且转子相对位置发生改变并不影响效率-时序位置曲线;总轴向间距变大,时序效应变弱,效率提升幅度仅为0.7个百分点,效率-时序位置曲线随之改变。     

螺旋桨飞机滑流非定常数值模拟研究

采用非定常数值模拟方法,基于滑移网格技术,分析研究螺旋桨滑流对全机气动特性的干扰影响。用多块结构网格数值模拟软件,从搭接边界建立、搭接网格生成、流场插值传递3个方面对螺旋桨滑移搭接网格的数值模拟方法进行全面详细的介绍,并且阐述基于该技术的非定常数值模拟方法。在此基础上,开展螺旋桨飞机模型验证计算,分析研究螺旋桨滑流对全机流场和气动性能的影响,通过与风洞实验数据的对比分析,验证滑移搭接网格技术的可行性与准确性。研究结果表明,基于滑移搭接网格技术的螺旋桨滑流数值模拟,可以真实反映螺旋桨滑流流场特性,为螺旋桨类飞机的设计评估提供一种有效的数值模拟手段,在工程方面具有重要的实用价值。     

低马赫数下多凹腔燃烧室非稳态燃烧过程

作为稳定火焰的有效手段之一,凹腔构型在冲压发动机燃烧室研究中占有重要地位。在对以煤油为燃料的多凹腔燃烧室冷/热态流动特性分析的基础上,重点研究低进口马赫数条件下燃烧室点火起动初期非稳态过程。结果表明：上游凹腔内大涡结构有助于提高燃料的驻留时间,未燃混气被高速主流带入下游凹腔内继续反应,进一步提高燃烧效率;燃油喷射速度决定被卷吸进回流区的燃油质量分数的大小,进而影响燃烧效率高低;燃烧室点火起动初期出现了主流熄火、火焰逆流传播以及主流再着火等复杂现象,火焰逆流传播现象是在上游凹腔内燃料自燃与下游燃烧释热压缩来流两种机制共同作用下完成的。