基于三维非结构网格的连续伴随优化方法

基于非结构网格,采用Euler方程、连续伴随方法以及自由形面变形技术(Free-form deformation,FFD)结合序列二次规划(Sequential quadratic programming,SQP)优化算法,开发出了三维飞行器的气动外形优化代码。应用Fortran 90语言,动态存储数组和指针,节省了内存。采用Roe形式的二阶迎风格式和LU-SGS隐式时间离散格式求解连续伴随方程,同时采用一种新的参数化坐标求解方法,使得FFD技术更加鲁棒。和有限差分方法求得的梯度进行比较,验证了连续伴随方法所求梯度的正确性,实现了ONERA M6机翼的减阻优化,证明本文开发的优化代码能够实现三维气动外形的减阻优化。     

DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

基于风洞试验的非定常气动力微分方程建模方法

为准确描述大攻角非定常运动中非定常气动力的气动特性,研究了以微分方程为基础的非定常气动建模方法。在某战斗机大振幅偏航、滚转单自由度及偏航-滚转耦合运动风洞试验的基础上,对偏航力矩、滚转力矩建立微分方程模型,分析研究不同频率、同一侧滑角情况下,气动力的迟滞时间;然后对高频运动下的气动力进行拟合,进而建立高频运动下的气动力微分方程模型。研究表明:模型可精确预测不同机动下非定常气动特性;建模方法具有较强的工程可行性。     

月面环境三维激光SLAM技术

同步建图与定位（SLAM）可实现月球车在未知复杂月面环境下的定位与导航,月球表面由陨坑、石头等起伏地形构成,缺乏树木、建筑物等地面存有的显著特征,大量特征不显著的点云数据会对月球车定位精度和实时性造成影响。本文提出了一种针对月面环境的显著特征点云提取方法以及基于曲面定位能力估计的增量式优化算法,通过Fisher信息矩阵计算曲面定位能力指标,获取机器人位姿估计的不确定性测量,利用增量式的SLAM方案进行优化,用于提高定位精度与实时性。通过在Gazebo （物理仿真平台）仿真场景下的测试,验证了算法性能。     

基于傅立叶平均法下的连续小推力动力学分析

通过将小推力展开为偏近点角的傅立叶级数,并对高斯摄动方程在一个轨道周期上的平均,将原方程的推力转化为仅由14个傅立叶系数表示的控制变量。仿真计算表明,平均化后的高斯方程使计算量与牛顿积分相比显著减少,且对小推力而言有足够的精度。对利用平均化后的高斯方程计算轨道根数时产生误差的原因进行了研究,并进一步分析小推力的范围和小推力近似表达式对上述误差的影响,为今后小推力下非开普勒轨道动力学分析提供了理论依据和参数。     

轨道要素奇异问题的改进四元数方法

运用四元数方法可以在一定范围内解决描绘飞行器轨道运动时轨道要素的奇异问题。但四元数固有的双值性使得在对运动方程进行积分时,其正负选取很困难。为了解决这一问题,采用了改进约束方程的四元数方法,并用该方法描述了拉格朗日行星摄动方程,然后研究了地球扁率摄动对地球同步卫星轨道的影响。仿真结果表明：当偏心率小于1时,四元数可以很好地解决轨道要素奇异性问题。与改进的春分点轨道要素相比,四元数的方法有着更加明确的物理意义和几何意义,用四元数表示的运动变量方程的计算更为简单,积分计算效率更高,而且其计算误差也能达到精度要求。     

Fault diagnosis and isolation of the component and sensor for aircraft engine

Aircraft engine component and sensor fault detection and isolation approach was proposed,which included fault type detection module and component-sensor simultaneous fault isolation module. The approach can not only distinguish among sensor fault, component fault and component-sensor simultaneous fault, but also isolate and locate sensor fault and the type of engine component fault when the engine component fault and the sensor faults occur simultaneously. The double-threshold mechanism has been proposed, in which the fault diagnostic threshold changed with the sensor type and the engine condition, and it greatly improved the accuracy and robustness of sensor fault diagnosis system. Simulation results show that the approach proposed can diagnose and isolate the sensor and engine component fault with improved accuracy. It effectively improves the fault diagnosis ability of aircraft engine.      

非结构网格下可压缩黏性流体的隐式算法

采用牛顿线化和Roe近似Jacobians矩阵隐式算法,实现基于非结构网格的三维定常黏性流动问题的求解.分别对平板、跨声涡轮叶片Mark-Ⅱ和CW-22这三种具有代表性的流动进行数值模拟,并详细分析了流动细节,模拟结果与理论以及实验值吻合良好,准确预测了流场内激波、二次流等重要的流动现象.验证了基于非结构网格的可压缩流体隐式计算方法对于黏性流动问题的计算是正确的.      

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

CE／SE并行计算研究及在Euler方程中的应用

时空守恒元和解元方法(CE/SE)是一种高精度的数值求解方法.采用改进的CE/SE方法和并行分区技术,针对非结构网格,发展了一套求解二维Euler方程的并行程序.对NACA0012翼型的亚跨声速流动和多段翼型复杂绕流进行了分区并行求解,多区并行计算的结果与实验的结果吻合较好,表明了本文研究方法的正确性和有效性.同时文中的分区方法能实现各处理器之间的负载均衡,有效地节省计算时间.