对数值模拟气相爆轰中二阶段化学反应模型再研究

采用改进的高精度时-空守恒元解元算法(the space-time conservation element and solution element method,CE/SE method)和考虑组分的二阶段化学反应模型(Sichel的二步模型)对气相爆轰问题的数值模拟进行了分析.分析发现采用Sichel的二步模型得到的数值结果虽然比早期二阶段化学反应模型(旧二步模型)更接近实验值,但是仍然不能得到爆轰过程准确气体动力学参数.为此通过修改组分的质量分数分布形式对Sichel的二步模型进行了改造,然后采用新的二步模型对平面爆轰波进行了数值模拟.数值结果表明采用新的二步模型计算得到气体动力学参数更接近于实验值和基元反应模型的计算值,在计算精度上有较大提高.     

卫星轨道积分器的变步长外推算法

通过对卫星轨道的变步长数值积分算法预测公式和校正公式的误差进行计算,将预测值与校正值进行组合外推,消除掉误差的主要部分,并基于此提出了新的误差估计和步长选择公式。仿真试验表明,外推算法对数值计算的效率有所改善。     

局部流动分离诱导结构振动问题的数值模拟

基于非结构网格技术和SST-DES模型,采用流体与结构紧耦合策略,对飞机机身突出物（如雷达罩）下游的蒙皮振动过程进行了数值模拟计算。在求解NS方程过程中,使用有限体积法,结合有限元法求解结构运动方程。在时间域内模拟了不同结构参数和飞机不同迎角下的飞机蒙皮的结构响应,建立了亚音速条件下壁面流固耦合分析方法。计算结果表明：机身突出物诱导出的分离涡的确可以激发下游壁面的振动,其壁面振动的幅度与结构本身参数密切相关,振动主要为发生在壁面的法向上,随着迎角增加,壁面的振动减弱。通过基本原理和算例可以看出,方法可以预测不同条件下壁面的动态响应过程,适用范围宽,可以用于壁面结构的设计。     

一种基于试验的FLD算法

为方便科研和生产实践，本文提出了一种新的基于试验的FLD算法，即用幂次式建立FLD数学模型，经过对试验数据的拟合，计算出FLD曲线，进行实时绘图和相关性统计分析，试验和理论分析均表明，通过计算所得的FLD可以较真实地反映板材的局部成形性能，有一定的实用价值。     

激波风洞高超声速摩阻直接测量技术研究

介绍了在中国空气动力研究与发展中心（CARDC）激波风洞中进行的摩阻测量技术研究情况。在测量研究中,设计了压电型摩阻天平,为了提高摩阻天平的校准和风洞试验测量结果精度,便于风洞试验和校准之间安装的变换,本项研究的摩阻天平采用一种新结构,也就是测量表面和摩阻天平本体可以分离的分体式结构,由此确保在不同使用场合下,摩阻天平的测量表面或者校准加载块可拆卸和更换。验证性试验是在CARDC0．6m激波风洞中进行的,流场名义马赫数分别为8和10,单位雷诺数分别为2．85×10^7／m和1．58×10^7／m,试验中测量了带压缩拐角的进气道模型表面三个测点的摩擦阻力,也测量了摩阻测点及其附近热流,测量结果表明：模型表面的摩阻和热流与雷诺比拟准则符合得较好。     

关节臂式三维测量对复杂铸件验收的应用

针对大型复杂铸件毛坯加工前验收现有钳工划线方法操作难度大、耗时费力、效率低等难题,提出了一种关节臂式三维测量新型便捷的测量方法。首先从最终产品数模中提取待测物体数模坐标系,选取基准点坐标,建立工件坐标系,然后利用关节臂对铸件实物表面进行轨迹测量实施,最终按测量需要,进行数据处理,将扫描的铸件实物三维形貌与最终产品数模对比,可直观得出各铸造面铸造余量及分配情况。通过航天产品遥感器相机典型零件可见三镜支撑架和桁架的测量,对比了与钳工划线方法的差异,解决了传统划线方法无法克服的技术难题和难点。经实物验证表明,该方法对空间复杂铸件,可实现快速、准确的加工前验收,能够提高铸件产品验收质量和效率。     

星敏感器光轴测量的锥面误差模型及其应用

星敏感器是卫星高精度姿态测量系统中的重要器件,根据其测量可以建立不同形式的星敏感器测量模型,准确分析测量模型中测量误差的特性是保证卫星姿态确定精度的重要条件。本文针对广泛采用的星敏感器光轴矢量测量方程,从几何角度出发,结合实际情况,提出了一种新颖的星敏感器光轴测量的锥面误差模型。在两个锥面参数分别服从一定的概率分布的条件下,对测量误差特性进行了深入分析,确定出新的测量误差协方差矩阵。上述研究能够在不增加算法计算量的前提下,从新的角度,更为直观地建立了星敏感器光轴矢量测量模型。最后在仿真实验中,将新的锥面误差测量模型应用于基于扩展卡尔曼滤波（EKF）的姿态确定方法中,结果表明了利用该测量模型进行姿态确定的有效性。     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

典型气动布局高超声速飞行的气动力数值评估

采用自主开发的数值软件SPACER对典型气动布局在Ma=6.0(Re=1×107)高超声速巡航时的气动力性能进行了数值评估。评估的气动布局包括类乘波体(仿X-51A)、翼身融合体(仿ISR)、传统升力体(仿X-33)和轴对称锥形体(仿Fasthawk导弹)。由于升力体范围较广,还对一种相对扁平的升力体模型进行了数值评估。评估结果表明:在考察条件下,类乘波体具有较大的升阻比和较小的阻力,是很有潜力的高超声速巡航气动布局;翼身融合体和传统升力体需一定的改进和进一步的研究,如扁平升力体的气动性能可得到大幅提升;轴对称锥形体以其总阻力小,也具有一定优势。