高超声速滑翔飞行器再入轨迹快速、高精度优化

针对高超声速滑翔飞行器再入轨迹优化问题，提出一种基于稀疏差分法和网格细化技术的快速、高精度求解方法。该方法应用局部配点法将再入轨迹优化问题转化为非线性规划(NLP)问题，从两方面提高轨迹优化的效率和精度。一方面，引入一种高效的稀疏差分法计算NLP的一阶偏导数，提高NLP的求解效率；另一方面，提出一种基于新型广义二分网格的网格细化算法调整离散节点的数量和分布，使得方法能够采用较少的节点数目取得较高的优化精度，从而减小NLP的规模和计算量。应用该方法求解了高超声速滑翔再入轨迹优化问题，仿真结果表明所述方法能够快速生成一条严格满足各种约束的最优三维再入轨迹。在此基础上，研究了滑翔飞行器的再入落点区范围，进一步检验了该方法的有效性。     

升力偏置对共轴刚性旋翼前飞气动特性的影响

为了分析升力偏置对共轴刚性旋翼前飞气动特性的影响，建立了基于雷诺平均Navier-Stokes方程的计算流体力学方法进行共轴旋翼流场求解，采用嵌套网格方法模拟桨叶运动，采用双时间方法进行时间推进。针对不同升力偏置状态，采用基于"差量法"的共轴旋翼高效配平策略进行操纵量配平。通过对Harrington-1旋翼性能的计算，验证了方法的有效性。对比计算了共轴刚性旋翼在不同前进比和升力偏置量下的气动性能和流场特征，结果表明：双旋翼操纵量在小前进比状态有明显差别，在大前进比状态基本一致；在相同拉力状态，随着升力偏置量的增大，共轴旋翼升阻比先升高后降低，其阻力却不断增大，不同前进比状态的最大升阻比对应的升力偏置量不同；双旋翼相遇时桨叶拉力出现脉冲式波动，由于流场被前行桨叶所主导，因此后行桨叶拉力波动幅值更大，且波动幅值随升力偏置量的增加而增大。     

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

应用TTM网格研究冲压增程弹丸进气道内外流场

基于TTM贴体网格生成法,研究了某冲压增程弹丸超音速进气道内外流场数值模拟中单区域网格生成问题,流场求解采用了二阶隐式TVD格式。NACA0012翼型跨音速流场算例结果表明,所采用的算法和编制的程序可靠,能够适用于具有亚、跨、超音速三种流态共存的复杂流场数值计算。对超音速进气道数值模拟结果显示,得到的流场复杂波系结构是合理的,来流马赫数和攻角变化显著影响着进气道性能。     

统一系统载波捕获方法研究

介绍了现代统一系统中载波频率快速捕获的几种方法,并提出了统一系统载波快速捕获的一种算法。最后,通过对几种不同方法的计算时间与识别精度的比较,验证了该算法的高效性。     

一种求解轨迹优化问题的改进多分辨率技术

针对轨迹优化多分辨率技术可能发生细化遗漏或者细化失败的缺陷，设计了一种改进多分辨率技术。该方法通过在每次网格细化迭代过程中增加两层高分辨率节点并且引入节点检测算法，确保在轨迹的非光滑区域能够高效、持续地细化网格。针对原始多分辨率技术的初始网格的节点数目存在二进限制的弊端，定义了一种新的广义二分网格，使得多分辨率技术的初始网格节点数目可以为任意奇数。采用多个轨迹优化算例验证了本文方法的有效性并且与其它方法进行了对比。此外，仿真结果还表明，网格充分细化对目标函数影响不大，但是能够更加准确地满足路径约束和终端状态约束。     

基于结构网格的操纵面嗡鸣数值研究

利用无限插值法生成结构化动网格,耦合求解欧拉方程与结构方程,分别模拟了二维带操纵面翼型以及三维NASP机翼周围跨音速流动。对二维带操纵面翼型嗡鸣发生时的流场进行了分析,发现操纵面上的激波后会产生巨大的分离区,且分离区的变化过程滞后于激波变化过程,认为缝隙的存在对嗡鸣发生有促进作用;对三维NASP机翼的研究发现操纵面质量比的增加可以延缓嗡鸣发展。     

超声速民机发动机短舱布局对声爆的影响

超声速民机较高的飞行速度导致发动机短舱与机翼机身存在强烈的气动干扰，且高速喷流会改变主流的激波系结构，进而影响声爆强度。以超声速民机为研究对象，对不同短舱布局开展数值模拟研究，揭示发动机短舱位置、数量对近场过压信号特性的影响规律。数值模拟运用有限体积法在直角网格上求解流体控制方程进行，通过伴随自适应网格加密生成高分辨率网格以精确捕捉近场过压信号。结果表明：发动机进气道唇口激波、尾喷口后缘激波及尾喷流与机体机翼引起的激波系存在强烈的相互干扰，一定程度上增加了近场过压幅值，从而增强了声爆强度。短舱沿弦向前移及置于机翼上方可以有效降低声爆强度，沿展向外移通过抑制尾部各激波的合并也可有效降低声爆强度；在相同总推力前提下，相比三发构型，双发构型能有效降低后体激波强度，但较大尺寸的短舱引起较强的进气道唇口激波。综合考虑喷流噪声和气动阻力因素，翼下双发布局对于新一代超声速民机并非最佳选项。     

霜状冰结冰数值模拟研究

提出了一种三维霜冰结冰的数值模拟方法。基于欧拉-拉格朗日法计算空气-水滴两相流流场,获得了表面水滴撞击特性以及撞击量,求解质量守恒和能量守恒方程获得网格单元内的结冰量,三维霜冰的结冰计算程序采用用户自定义（UDF）编写,FLUENT提供的动网格功能实现了结冰后的网格重构。计算了4°攻角,液态水含量为1 g/m3,水滴直径20μm,结冰时间6 min,温度分别为-19.4°和-28.3°时NACA0012翼型的结冰,并与相关文献结果进行了对比,计算结果与文献实验结果吻合较好,验证了方法的合理性和可靠性。分析了不同的粒径、液态水含量、来流速度等对撞击特性以及结冰的影响,随着粒径、来流速度的增加,水滴收集系数、撞击区和结冰厚度增大,随着液态水含量的增加,水滴收集系数和撞击区是不变的,但是结冰厚度增加。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。