新的解跨音速Euler方程的隐式杂交方法

 本文提出一个新的杂交格式,它由Jameson-Turkel LU分解法和Harten TVD（total variation diminishing）格式的修正数值通量耦合而成。LU分解用于Euler方程隐式部分的离散化而Harten的数值通量控制所谓的稳态解残差。精心构造LU算子,可使矩阵算子L和U很容易逐点推进求逆;以Harten的数值通量计算杂交格式右边的通量,这样做在a=1/2的条件下,可保证收敛的稳态解具有二阶精度并能高分辨率地捕获激波。数值实验表明,本文格式对于稳态计算是有生命力的并能保证较高的激波分辨率。     

超声速进气道流场和边界层计算

对设计马赫数为3.5的混压式进气道,在小迎角下作了流场和边界层计算。用内激波捕捉法计算无粘流场,用积分法计算层流和湍流边界层。结果和实验数据以及特征解一致,边界层位移厚度也合理。     

不可压Navier—Stokes方程组的SUPG有限元数值解

本文从定常不可压NaVier-Stokes方程组出发,构造了SUPG加权剩余公式。为保证数值解的精度,本文对速度取八节点插值,保留了摄动项中的二阶导数项。从用本文方法所做的算例来看,计算结果是令人满意的。     

叶片通道内二次涡的数值模拟

采用三维Ｅｕｌｅｒ方程数值模拟叶片通道内流向二次涡,及其对叶片出口气流角的影响。算例表明：采用这种方法可以计算出流向二次涡;该涡造成的出口气流偏转角与实验值比较,在远离壁面处相差较小。此方法可用于叶栅通道内流向二次涡的初步分析     

环流喷管的轴对称湍流流场计算

对火箭发动机环流喷管进行了二维数值仿真,用数值方法求解Navier－Stokes方程,数值格式为显格式时间推进格式,湍流模型选择可靠性较好的两层代数湍流模型。由于所研究的流场边界形状复杂,计算中采用了分区的方法。通过多个算例的对比,确定了使环形喷管推力效率较高的几何形状。画出了环流喷管内部的计算网格图、速度矢量图、等马赫线图、等压线图及喷管上壁面及轴线的压力分布图。     

基于全局敏度方程方法的飞机方案优化设计

针对飞机方案分析模块之间存在耦合关系,从而使方案分析需要迭代求解这一特点,采用全局敏度方程方法计算参数之间的敏度,以克服差分方法所带来的高费用和不精确问题。该方法本身具有递归和并行特点,使得其与数值优化方法结合起来特别适合于像飞机这类复杂的内部耦合系统的优化。给出了该方法与约束变尺度法结合,用于某假想的喷气教练机方案优化的简单算例以及用于某4座公务机方案优化的实例。     

提高火药燃气驱动冲击式涡轮效率的途径探讨

探讨了以火药燃气作为工质的运载火箭中所用的冲击式涡轮效率提高的途径。给出了具有小反力度冲击式涡轮效率的计算方法和火药燃气中含有微粒流时冲式涡轮效率的迭代计算方法。对四种情况下计算结果进行了对比分析，计算结果与试验结果符合良好。     

N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

高精度的对流差分格式

对于广义对流差分格式,本文应用两种方法(截断误差分析法及事后估计误差法)导出了使其误差最小的最优α值。和以前的方法不同,它们是从所研究的对流-扩散方程本身出发导得的,而无需事先知道方程的解析解。大量算例结果表明,截断误差分析法有很高的精度。事后估计误差法的误差,除个别点外,也都比一般格式的误差要小或者小得多。但是,对于复杂的问题,前者的应用可能会有困难,而后者却很易实现。