改进的格网惯导系统无阻尼综合校正方法

综合校正技术可作为抑制格网惯导系统（INS）导航误差的有效手段。加速度计零偏所造成的水平姿态误差是导致综合校正中陀螺漂移估计精度受限的重要因素。针对这一问题，提出了一种改进的无阻尼综合校正方法。首先，推导了格网坐标系框架下估计加速度计零偏和姿态的目标函数；其次，介绍了无阻尼条件下综合校正的两个核心方程；最后设计了无阻尼两点校策略。综合校正前，在多普勒计程仪（DVL）提供的速度辅助下完成加速度计零偏的估计和补偿，以此消除加速度计零偏所造成的水平姿态误差对综合校正中陀螺漂移估计精度的影响。在两次间断的外部位置和航向辅助下通过所设计的综合校正策略完成对陀螺漂移的估计。校正策略中所涉及的水平姿态误差在DVL辅助下由参数估计方法估计得到。仿真及实验结果表明：与现有的研究相比，所设计的综合校正方案进一步减少了DVL的辅助时间，同时由于准确地估计和补偿了加速度计零偏，陀螺漂移的估计精度显著提高，该方案在抑制导航误差方面具备更明显的优势。     

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

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

基于非结构嵌套网格的涵道螺旋桨数值模拟

基于非结构嵌套网格方法,通过求解可压缩Euler方程对涵道螺旋桨系统的非定常流场进行了数值模拟,并研究了桨尖间隙、螺旋桨位置、涵道尾扩张角等参数的影响规律.将计算区域分为相互重叠的旋转区域和静止区域,旋转区域包含螺旋桨和旋转的桨毂,静止区域包含涵道和部分静止的桨毂,采用嵌套网格算法真实地模拟了螺旋桨和涵道之间的相对旋转运动,得到了细致的螺旋桨流场.计算值和实验测量值的对比结果相一致.参数影响研究表明,减小桨尖间隙、在一定范围内后移螺旋桨或增加尾扩张角,均可以提高涵道螺旋桨的气动性能.     

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

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

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

离子推力器加速栅寿命概率性分析

交换电荷离子对加速栅极的溅射腐蚀是离子推力器的关键失效模式之一,基于交换电荷离子对加速栅溅射腐蚀的物理机理,对离子推力器加速栅工作寿命进行了概率性建模。利用该模型对20cm Xe离子推力器加速栅寿命和其达到预期寿命的可靠度进行了评估。结果显示加速栅的寿命近似服从高斯分布,当推力器工作环境压力近似6.7×10-3Pa时,加速栅工作寿命达到3kh的可靠度为0.9352。     

An efficient large-scale mesh deformation method based on MPI/OpenMP hybrid parallel radial basis function interpolation

An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption. Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation, the present method employs different surface points sets and influence radius for each volume point movement, accompanied with efficient geometry searching strategy. The deformed surface points, also called Control Points (CPs), are stored in each processor. The displacement of spatial points is interpolated by using only 20–50 nearest control points, and the local influence radius is set to 5–20 times the maximum displacement of control points. To shorten the searching time for the nearest control point clouds, an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique. Besides, an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications. Three 3D cases, including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements, are used to test the present mesh deformation method. The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 45° and more than 80% parallel efficiency with 1024 MPI processors. In addition, the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.     

网格自适应策略在高升力构型计算中的应用

提出了一种基于多尺度网格误差分析的新型混合网格自适应局部加密策略。在初始网格基础上,通过多重网格计算,比较次细网格和细网格的计算结果,揭示出网格敏感区域,然后对这些网格敏感区域加密,计算自适应加密后的网格得到最终计算结果。相比于常规网格自适应策略,提出的新策略特别适合多种流动现象并存的高升力构型绕流流场计算中的敏感区域判断。NASA高升力全展构型的数值试验充分表明这种新策略的有效性。     

COMPUTING REACTIVE FLOW ON A COARSE GRID

INTRODUCTIONConsider the Euler equations supplementedby an additional reactive equation which consistsof a scalar balance law for the mass fraction ofunburntgasρρuρEρY t+ρuρu2 + pρu E + upρu Y x=κ000-ρΘ( T) Ywhereρ is the density,ρu the momentum,ρEthe total specific energy,ρY the mass fraction ofunburnt gas( 0≤ Y≤ 1 ) ,κ a large number( thereaction rate) and T temperature.Θ( T) =01 　if T≤ Tignif T>Tign,Tign is the ignition temperature,E=pγ- 1 + u22 + q0 Y,p is …     

ERROR ANALYSIS IN DIFFERENTIAL QUADRATURE METHOD

ＥＲＲＯＲＡＮＡＬＹＳＩＳＩＮＤＩＦＦＥＲＥＮＴＩＡＬＱＵＡＤＲＡＴＵＲＥＭＥＴＨＯＤ￥ＨｅＢｅｉｑｉｎｇ（ＣｏｌｌｅｇｅｏｆＳｃｉｅｎｃｅ，ＮＵＡＡ２９ＹｕｄａｏＳｔｒｅｅｔ，Ｎａｎｊｉｎｇ２１００１６，Ｐ．Ｒ．Ｃｈｉｎａ）ＷａｎｇＸｉｎｗｅｉ（Ｄ...