三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

刚性方程的二阶Gear方法的注记

构造了一个求解刚性方程的二阶差分格式,它是隐式的,并采用简单迭代求解,对于稳定系统,迭代收敛性条件与步长无关,而对不稳定系统,迭代收敛性条件对步长却有所限制,利用这一限制性条件去克服解刚性方程的二阶Gear方法时出现的“危险性”。     

超声速开式空腔流动的数值模拟

研究了不同长深比、不同雷诺数条件下,超声速开式空腔流动特性的变化。算法采用空间上5阶的WENO格式,时间上具有TVD性质的3阶Runge-Kutta方法,湍流模型选用B－L模型。数值模拟获得了不同的时刻的注场图像,压力波、涡以及空腔结构之间的相互耦合作用,使得空腔内存在典型的周期性压力脉动,但不同长深比空腔的流动性存在较大差别。     

几种常用格式和湍流模型在高速流动中的适用性研究

针对高超声速飞行器中存在的强激波、激波/边界层干扰、分离、湍流等复杂流动现象,对比分析了当前计算流体力学中的主要空间离散格式及湍流模型,发现不同格式对强激波的分辨率基本相同,Roe和LDE(low diffusion E-CUSP(convective upwind split pressure))格式对摩擦因数和传热系数的模拟优于其他格式;S-A(Spalart-Allmaras)一方程湍流模型计算的摩擦因数比k-ωSST(shear stress transport)两方程湍流模型高10%左右,而后者预测的分离区约为前者的2倍,且分离点靠前.     

常规装药近地爆炸冲击波效应数值模拟

针对现有经验公式在计算爆炸地面冲击波超压值时存在的局限,采用NND2M差分格式对TNT柱形装药在小爆高状态下的爆炸冲击波效应进行了数值模拟研究,并与试验实测值及经验公式计算值进行了比较,分析了存在差异的原因,证明了本文采用的计算是合理可靠的.     

强迫耗散非线性发展方程与完全平方守恒格式

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散非线性发展方程的显式准安全平方守恒差分格式。理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的，值得推广应用。     

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率,对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序,该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想,验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合,说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动。为了提高计算规模及计算速度,作者对程序进行了并行化并针对微机机群系统进行了并行优化。实际计算表明本程序具有较高的并行效率。      

旋转激光陀螺惯导系统误差传播特性分析

设计旋转式惯导系统(INS, Inertial Navigation System)最重要的工作是理清旋转调制误差自动补偿机理、掌握补偿前后误差的传播特性,在此基础上设计合理的转位方案,并对其进行仿真和试验验证.结合单轴旋转式惯导系统工作原理和误差传播方程,解释了旋转调制误差自动补偿的机理,利用理论和仿真两种手段分析了单轴旋转惯导系统中惯性元件常值误差、随机误差、标度因数误差和安装误差的传播特性,得到了旋转调制对惯性器件误差的调制效果,验证了单轴正反转停方案的有效性和应用的合理性.研究结果为旋转式激光陀螺(LG, Laser Gyroscope)惯导系统的设计提供了理论参考和设计依据.     

一种双钟型喷管液氧/甲烷发动机系统方案

根据双钟型喷管高度补偿特点及技术研究现状,提出了一种双钟型喷管液氧甲烷发动机系统方案,进行了双钟型喷管基弧段及延伸段面积比优化,并与其他系统方案进行了性能对比分析。研究表明,对于地面起动的芯级发动机,采用双钟型喷管是提高发动机综合比冲性能以及运载器有效载荷的有效途径。