飞机结冰中水滴撞击特性欧拉法的网格影响分析（英文）

过冷水滴撞击特性计算研究是飞机结冰预测与防除冰系统性能分析的基础,欧拉法是过冷水滴运动及撞击特性计算的常用方法之一,分析欧拉法在计算复杂表面水滴运动撞击时的准确性具有重要意义。以NACA0012翼型、冰风洞风道、S形进气道及三段翼为研究对象,利用欧拉法采用不同网格计算获得了水滴运动及局部水收集系数,结果表明：当水滴运动没有受到上游部件的导流或者遮挡时,欧拉法计算结果具有网格无关性;当水滴运动受到上游部件影响发生轨迹偏转时,欧拉法计算结果具有网格相关性,不同网格计算得到的液态水含量及局部水收集系数不一致,需要进一步考虑水滴运动受到上游效应影响时网格对欧拉法计算结果的影响。     

随机场TPDF在非结构网格平台的开发及检验

为实现对燃烧室内部燃烧过程的高精度数值模拟,增强航空发动机燃烧室研究手段,采用开源的非结构网格平台Saturne,开发可以耦合详细化学反应机理,高精度解析湍流燃烧过程的概率密度函数输运方程（TPDF）湍流燃烧模型。采用随机场方法求解TPDF方程,在原有程序基础上研发了TPDF程序模块及配合的加速算法等模拟单元,发展了针对燃烧室燃烧过程模拟的软件功能。采用射流火焰Flame D和旋流火焰TECFLAM对新软件进行了测试,结果表明：模拟结果与试验数据趋势一致,精度尚可。对某型燃烧室性能模拟验证中,计算的出口温度分布与试验结果一致,平均温度误差小于3.8%,新的软件可进一步应用于燃烧室性能的研究。     

高速进气道抛罩过程内流场演化特性研究

为探究进气道整流罩打开过程中内流场的演化规律,针对下颌式进气道及其配套整流罩,基于重叠网格法求解N-S方程,详细分析了抛罩过程中内流场演化过程与特征流场结构。研究表明,堵罩时,整流罩前端发展出V形激波,整流罩尾部形成后向台阶且产生一对涡矢量方向相反的流向涡,开罩过程中,整流罩前端的V形激波与尾部的后向台阶均会对内流场流动产生影响,捕获流量在中间某时刻达到峰值且大于稳定状态。对于开罩过程流动准定常且内流道无大分离区的流场,进气道在V形激波打在唇口时完成起动且整流罩在此后不影响内流场流动。假设的匀加速转动与重构出的真实转动吻合较好,对于更精确的运动过程模拟可采用三次函数的角速度变化律。     

多维对流通量重构方法的改进

为了研究多维对流通量重构格式中不同形式主方向对激波稳定性、压力场计算方面的影响,分别采用以笛卡尔方向、密度梯度及压力梯度作为主方向进行计算,对结果分析后发现,单纯依赖物理量梯度的主方向会造成压力等值线曲折的非物理解。为了消除这种非物理现象,采用一种基于压力函数的混合方式来确定主方向。单层及展向多层高超声速圆柱绕流问题模拟分析的结果表明:混合方式所确定的主方向不仅可以消除笛卡尔坐标系确定的主方向计算带来的红玉现象,同时还可以消除单纯的密度或压力梯度作为主方向所带来的压力场"zigzag"的非物理现象;通过与结构网格结果的比较表明:网格相关格式驻点处热流的相对误差为29.2%,而混合主方向确定方法所计算驻点热流相对误差只为4%,表明此改进方法计算得到的壁面热流是合理的。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

动基座近舰面流场数值模拟

航母尾迹流场对舰载机的着舰有较大影响,所以需要对其流场特点进行研究,分析不同状态下舰载机气动特性的变化。采用嵌套网格技术对航母处于垂荡状态下无人机的着舰进行了模拟。首先,利用SFS2舰船进行数值计算,验证了舰船流场的数值模拟方法。然后,对比了单独无人机定常与非定常计算结果,表明所建立的嵌套网格适用于无人机流场的模拟。接着,对航母单相流和两相流的流场结果进行了分析,结果显示甲板下方的流动对甲板上方流场没有大的影响。因此,忽略了水的影响只对航母在空气流场中的特性进行研究,结果表明航母尾迹非定常特性明显,静止航母下滑轨迹上的速度均处于周期性波动状态,且波动幅值随着远离航母而逐渐衰减;而在垂荡情况下航母尾迹变得更加紊乱,水平方向速度波动的周期性减弱,但垂向速度的波动幅值进一步增大。对于静止航母,无人机在不同时刻着舰气动特性的变化也存在差异;当航母处于垂荡状态时,无人机的升力和俯仰力矩在短时间内会有更大的波动。     

涡轮平面叶栅非轴对称端壁优化设计

开发了一套造型灵活直观、网格生成速度快的涡轮平面叶栅非轴对称端壁优化设计工具,该工具的核心技术是非均匀有理B样条(NURBS)曲面造型和网格变形.在此基础上以商业软件Isight为优化驱动器,以CFX为求解器,搭建了非轴对称端壁优化设计流程.以Pack B涡轮平面叶栅为例,对其进行了非轴对称端壁优化设计.优化后涡轮平面叶栅总压损失系数减小了12.96%.结果表明:涡轮平面叶栅端部的静压分布改变削弱了涡轮平面叶栅通道中马蹄涡、通道涡的强度,提高了涡轮平面叶栅的气动性能.      

深空测控通信中GMSK体制非相干解调算法研究

针对深空测控通信中GMSK体制非相干解调损失较大的难点,提出了一种改进的GMSK信号非相干维特比解调算法。通过分析相位状态网格图中相位转移规律,建立理论仿真模型。通过原理样机的研制和测试,实测数据表明：该算法具有解调损失低、实现复杂度适中的优点;相比于理论的最佳相干解调算法,在误码率1×10 -4 量级下,实测仅损失0.6 dB。目前该算法已应用于国内某深空测控通信系统GMSK体制基带设备中,并成功解调出欧空局Herschel–Planck卫星数据。