CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes （RANS）方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会（DPW Ⅵ）。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明：迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。     

TB17钛合金高温压缩变形行为

通过Gleeble 3800热模拟试验机对TB17钛合金在变形温度860～980℃、应变速率为0.001～1 s~(-1)、最大变形量为70%下高温变形行为进行研究。通过材料参数与真应变之间的关系,利用Arrhenious本构方程关系式和Z参数建立流变应力和变形温度、应变速率和真应变三者之间的本构关系,并对组织进行分析。结果表明:TB17钛合金在应变速率为0.001～0.01 s~(-1)、变形温度为890～980℃下更容易发生连续动态再结晶,而在应变速率为0.1～1 s~(-1)下主要发生不连续动态再结晶;误差分析结果显示计算值与实测值平均相对误差为6%,说明建立的本构关系模型具有较高的准确度。     

基于足部安装MIMU的行人导航系统设计

针对无卫星环境中行人导航定位需求,设计了一种基于足部安装MIMU的行人导航系统,重点研究了行人导航鞋传感器选型及样机设计流程;并在此基础上设计了基于零速修正和航向修正的行人导航算法,研究了零速区间内最优航向角的提取方法,实现了对瞬时磁干扰产生的磁航向角的有效判断。最后,利用设计的行人导航系统开展了三维空间航迹测试实验,实验结果表明,在不增加高度约束的条件下,该行人导航系统定位精度为0.33%,能够满足坡路环境中行人的定位需求。     

千分尺示值误差分析及修正方法

详细分析了千分尺示值误差产生的原因,着重分析了丝杆螺距不均匀对示值误差的影响,介绍了示值误差的三种类型及修正方法.     

空间推进方法求解抛物化Navier-Stokes方程及其验证

采用伪时间空间推进方法求解完全气体PNS方程。PNS方程采用有限体积方法离散,在推进方向上采用一阶精度迎风格式,在推进面的展向和法向采用三阶精度的MUSCL插值和AUSM类通量构造格式;在推进面上的伪时间推进采用二维LU-SGS迭代方法。通过算例证明,空间推进方法能得到正确的压力、摩阻和热流分布,且计算时间比时间迭代方法快一个数量级以上。     

NS方程计算中耦合转捩自动判断的阻力精确计算方法初探

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合了流动转捩的自动判断以提高现有求解器预测翼型阻力的准确性.由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用简化的eN-数据库转捩判断方法分析层流边界层的解得到转捩点的位置,这样随着流场的迭代求解求解器自动判断转捩点的位置.在对NLF0416翼型的气动性能计算中考虑流动转捩的因素后得到的翼型升阻力特性和实验吻合较好,验证了本文方法的正确性.     

Z-A模型的修正及在预测本构关系中的应用

Zerilli-Armstrong(Z-A)模型是基于热激活位错运动理论建立的具有物理理论基础的本构模型,但是模型建立的初始假设,降低了其预测精度.为此通过引入考虑应变率和温度耦合效应的应变强化指数对Z-A模型进行修正,拓展了其适用范围,并将修正后的Z-A模型应用于描述IC10合金在不同温度和应变率条件下的塑性流变行为.与修正前Z-A模型的预测结果及试验结果比较表明:修正后Z-A模型的预测精度有较大程度的提高,最大平均相对误差值从6.41%降低到3.65%.      

侧向风速对飞机尾流运动的影响

为了使机场雷达能实时探测和定位飞机尾流,对飞机飞行环境进行合理假设后提出了一个仿真尾流的保守被动模型,该模型能很好地描述尾流中水蒸气等物质的运动演化规律。综合考虑大气环境中风速对尾流的影响,利用该仿真模型对尾流进行了实时仿真,得到了不同侧向风速情况下,不同时刻尾流的状态分布等重要特性。通过分析侧向风速对飞机尾流运动过程的影响,证明了当侧向风速为1.0~3.0 m/s时是最危险的,所得结论可为飞机飞行过程中避免尾流的影响提供理论依据。     

低马赫数流动中的预处理Euler方程的收敛特性

对预处理后Euler方程在低马赫数流动数值模拟中的收敛特性进行了细致研究。首先采用量级分析法研究了Weiss-Smith预处理方法对低速情形下Euler方程收敛性的影响;然后针对低速情形下预处理后的Euler方程组中各方程的收敛性,提出了一个改进的预处理矩阵。数值模拟结果表明：Weiss-Smith预处理方法可以很好地改善Euler方程在低速流动时的收敛特性,并且改进后的预处理矩阵在连续方程的收敛特性方面效果明显。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。