高精度差分格式及多尺度流场特性的数值模拟

简要介绍了作者近年发展的高精度差分格式,给出了两类格式精度的分析方法,分析了粘性面差分逼近所引入的耗散量低亏效应。文中按群速度将格式分为快型,慢型和混合型三类。对每类格式分析了高频波的传播特性。     

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

二阶时间精度的CFD/CSD耦合算法研究

非线性气动弹性分析中,涉及到非线性流体动力学(CFD)和非线性结构动力学(CSD)的耦合计算问题。本文分析比较了目前几种耦合算法:全耦合、松耦合和紧耦合,并从耦合边界能量传递守恒上就松耦合及紧耦合方法进行了时间精度的分析,得出传统松耦合中即使流体和结构子系统达到高阶时间精度,耦合算法的时间精度仍仅为一阶,紧耦合方法虽然可以达到二阶时间精度,但没有明显提高计算效率。随后,本文基于松耦合流程改进了耦合格式,分析表明其具有二阶时间精度,通过AGARD445.6机翼颤振模型的算例验证,说明该方法可以在保证计算精度的基础上明显提高计算效率,并保持了传统松耦合方法模块化的优点,在模拟非线性气动弹性问题时具有很高的优越性。     

基于全欧拉方程的二维平行剪切层声波产生和辐射的数值模拟

本文基于全欧拉方程对二维平行剪切层声波的产生和辐射进行了计算气动声学数值模拟。空间离散采用了频散相关保持有限差分格式,时间积分采用了低频散低耗散的龙格库塔法。数值模拟结果表明,当声源强度较小时,基于全欧拉方程的计算结果与解析解及基于线化欧拉方程的计算结果均符合得很好,当声源强度较大时(约140dB),线化欧拉与全欧拉的计算结果有较大的偏差,这说明此时非线性的影响不可忽略,线化方程不再适用。      

应用等效系统法评价飞行品质

本文介绍等效系统方法在评价飞行品质中的应用,包括等效系统方法的原理、所采用的数学模型和进行等效系统拟配的数学方法。最后,介绍了具体拟配过程中应当考虑的若干问题。     

色高斯噪声背景下相干信源DOA估计

为在色高斯噪声背景下估计相干信源DOA,提出一种基于四阶累积量的相干信源DOA估计算法。对各个阵元接收数据与参考阵元接收数据的四阶累积量进行排列,构造Toeplitz矩阵,该矩阵的秩仅等于信号个数,而与信号间的相干性无关。通过特征值分解得到信号子空间和噪声子空间,从而实现相干信源DOA估计。仿真实验结果表明,在色高斯噪声背景下,此算法能有效地估计出相干信源的DOA,并且相对于空间平滑算法具有更好的估计性能。     

多级嵌套L型阵列及其测向算法

针对基于偶数(2q)阶累积量的测向算法中测向性能提高有限的问题,提出了一种基于多级嵌套 L 型阵列的2维测向算法。首先利用阵列的多级嵌套结构和2q阶累积量,形成具有更多自由度的虚拟均匀面阵;然后使用2维平滑方法,恢复其2q阶累积量矩阵的秩; 采用2维MUSIC算法,进行方位角和俯仰角的估计。与常规的2q MUSIC 算法相比,所提算法不仅具有更好的测向精度,而且由于虚拟均匀阵包含更多的虚拟阵元,因此能够估计更多信源的方位角。另外,针对该L型阵列的最优配置问题,推导了各级子阵阵元数的最优和次优分配表达式。仿真结果表明这些结论的正确性。     

GPS/INS 组合导航系统降阶滤波器设计

本文研究了某型号飞机上ＧＰＳ／平台惯性组合导航系统中降阶卡尔曼滤波的设计方法。文中提出了一种降阶滤波器结构，并讨论了模型误差的伪随机噪声补偿方法。采用蒙特——卡洛分析法对降阶滤波器的实际误差进行了数字仿真，结果表明：本文提出的这种降阶滤波器能保证组合系统的导航精度，是一种具有实际使用价值的滤波器     

Revisiting the space-time gradient method: A time-clocking perspective,high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternativ...     

基于序列凸优化的高超声速滑翔式再入轨迹快速优化

针对具有热流、动压、过载以及多个禁飞区约束的再入轨迹优化问题,提出采用序列凸优化方法快速求解。利用归一化时间作为自变量解决终端时间自由问题,并引入辅助控制变量以减少序列优化结果中的高频振荡,在此基础上,通过线性化、离散化和非凸约束的凸化处理,将非凸非线性优化问题转化为二阶锥规划(Second Order Conic Programming,SOCP)问题,然后采用凸优化求解算法快速求解。数值优化结果与对比验证表明该方法能快速高效求解多约束条件下的再入轨迹优化问题,且计算效率和性能均优于传统的非线性规划方法。