非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.      

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ＥＮＯ（ＥｓｅｎｔｉａｌｙＮｏｎｏｓｃｉｌａｔｏｒｙＳｃｈｅｍｅ）型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶ＴＶＤＲｕｎｇｅＫｕｔａ方法．最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算     

低温气动阀动态响应特性的非线性分析方法

为了在液体火箭发动机试验过程中缩短低温气动阀的响应时间,减少非线性求解阀门动态响应的计算量,提升试验过程中的安全性与可靠性,根据低温气动阀的工作原理,构建阀门的动态响应非线性传递模型,并在此基础上,通过引入变步长系数与定义遗忘因子,在四阶龙格-库塔法基础上,构建了一种动态递推龙格-库塔法的非线性动力学求解方法,并提出响应速度的改进方案.最后通过实例的验证与仿真,为提升低温气动阀响应速度提供了一种方法指导,为提高液体火箭发动机试验技术提供了支持.     

仿真飞行轨迹的设计及应用

由于全面真实的飞行轨迹数据很难获取,而捷联惯导系统性能测试又必须要有飞行轨迹数据和惯性器件的输出作为基础,为此,介绍了一种仿真飞行轨迹的设计方法,将飞机飞行过程分解为若干基本运动状态,并建立了相应的模型,可根据实际需要组合各种基本运动状态构成一条完整的飞行轨迹.据此设计了一条仿真飞行轨迹并对该轨迹进行了应用仿真.结果表明,此法不仅能够较好地模拟飞机的实际飞行状态,更能有效地测试捷联惯导算法的精确性和可靠性,应用效果令人满意     

基于非结构网格流场计算的网格重排序

提出基于非结构网格流场计算时的网格重排序,优化了网格信息在内存中的存贮位置,使得计算过程中调用数据的处理速度加快,有效提高了计算效率。以M6机翼及DLR-F4翼身组合体的跨音速无粘流场为算例,分别通过显式四步龙格-库塔推进格式和隐式LU-SGS推进格式来验证网格重排序后对计算效率的提高作用。原始网格均用Delaunay方法生长,并分别使用了两种重排序方法进行比较。计算结果表明:重排序后的计算效率可以提高到未排序时的3到5倍,很具实用价值。     

叶轮机气弹耦合方程求解的时间推进方法评估

为了提高叶轮机颤振预测时的计算效率,同时兼顾求解的稳定性,有必要选取最适合的时间推进方法。基于降阶结构动力学方程的双向流固耦合方法常用于压气机气弹性能分析,针对气动弹性控制方程中结构动力学方程考查了多种时间推进方法对计算结果的影响。通过单自由度弹簧系统的数值求解探讨了不同时间推进方法的特点;选取实际压气机算例NASA Rotor 67进行颤振流固耦合分析,进行大量数值试验找到了各时间推进方法需要的最大时间步长,以此为依据对比出不同时间推进方法的计算效率。结果表明：采用单自由度弹簧系统和压气机颤振耦合求解得出的结论基本一致,说明时间推进方法的性能更多地与方法本身的数学性质有关,但是流固耦合求解时流场和结构的信息交互会带来一定影响。以实际颤振算例叶片瞬态响应振幅的对数衰减率作为衡量精度的标准,经典龙格-库塔方法是颤振耦合计算达到相同精度耗时最少的方法,其次是4阶隐式Adams方法和Newmark方法。     

旋翼桨叶绕流Euler方程数值模拟

本研究应用格心格式有限体积法求解欧拉方程,模拟绕旋翼浆叶的流动。旋翼尾流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角。选择O—H型式生成弦向和展向网格。通过算例计算为进一步寻求有效的Euler或NS方程数值解积累了经验。     

弹射救生数值仿真及不利姿态下救生性能分析

 国内弹射座椅救生性能预测分析主要依赖地面弹射试验。由于国内相关试验设备的缺乏,不利姿态下弹射座椅的救生性能很难通过试验获得,因此国内不利姿态下弹射座椅救生性能的分析研究比较薄弱。采用数值仿真的方法对弹射座椅的救生性能进行预测和分析,针对弹射座椅弹射过程的4个阶段分别建立了数学模型,其中自由飞阶段采用四元数法代替欧拉速率方程,以解决传统的六自由度方程出现的奇异性问题。采用四阶龙格-库塔法对数学模型进行求解,计算结果与地面弹射试验结果吻合较好。在此基础上,对某新型弹射座椅的弹射姿态轨迹进行计算,分析研究了弹射座椅在不利姿态下的救生性能,得到不利姿态下弹射座椅救生性能不佳的主要原因是座椅高速稳定性差以及飞机向下的牵连速度和飞机滚转角的影响,并根据分析结果提出了改善弹射座椅不利姿态下救生性能的基本方法。