弹道导弹制导计算中扰动引力的快速赋值

主动段扰动引力是引起弹道导弹制导方法误差的主要因素。因此,要提高导弹的制导精度,就必须能够在弹上实时计算扰动引力。但现有方法在计算快速性和存储量之间无法得到有效协调。为此,把广义延拓逼近思想引入有限元逼近方法中,将插值单元周围节点的信息也包含到单元内一点扰动引力的计算当中,建立了一种新的数学模型。对所选发射空域,在发射坐标系中进行了直角坐标划分。计算结果表明,这种方法能够更加精确地逼近弹道导弹主动段的扰动引力,在600 km×250 km×6 km的主动段飞行区域内,只需要保存60个节点数据,就能使由逼近误差导致的落点偏差小于10 m,是一般有限元逼近方法精度的4倍以上。     

在非结构网格上求解非线性航空声学问题的高精度有限体积法

介绍了有限体积数值方法,该方法适用于在任意非结构网格上求解线性和非线性的航空声学问题.本方法基于角点-中心的多参数格式,可在笛卡尔网格上达到六阶精度,对于可能的不连续性采用了自适应耗散.通过一系列算例研究了该方法的特性,结果表明:在模拟谐振型管中的噪声抑制中,所提出的方法是很有效的.     

驱动方腔内涡流的数值模拟

本文介绍一种求解涡度流函数形式的定常不可压缩Navier-Stokes方程的半隐式指数型差分格式。涡度方程用上述格式求解,而流函数方程用多层网格法求解。这种组合求解方式具有很好的稳定性及较快的收敛速度。本文对Re=100,400,1000,3000,5000的驱动方腔及Re=400,1000的驱动长方腔进行了数值模拟,结果与现有的计算相吻合。     

逆向射流多斜孔气膜冷却的数值模拟

为了使航空燃气轮机燃烧室在高于燃烧室材料许用温度的前提下能够正常工作,需对燃烧室内壁采取冷却措施,并对变吹风比和变角度下的逆向射流多斜孔气膜冷却进行了数值模拟。在贴体坐标下划分网格,采用k-ε模型模拟湍流流动,求解气膜冷却过程的数学模型。分别改变逆向射流同主流的夹角,使其为20°、30°、45°和60,°改变吹风比的变化范围为0.5到1.5等多种条件下,对多斜孔气膜冷却效果进行了数值模拟。在所研究的范围内得到了符合规律的结论。研究结果表明,当吹风比不变的情况下,减小逆向射流角度,用有效温比表示的冷却效果相应提高。当逆向射流角度不变时,增加吹风比,冷却效果也相应提高。     

涡轴发动机全状态实时气动热力学数值模拟

应用“外推”法得到了某型涡轴发动机起动状态的部件特性,采用容积法建立了表达发动机热力过程特性的非线性方程组,从而,解决了通常的部件法建模所带来的需要迭代的问题,建立了涡轴发动机全状态的气动热力学实时模型,给出了某发动机工作全过程的数值模拟结果。     

大曲率弯道和叶栅内的湍流数值模拟

本文在贴体曲线坐标系下应用非交错网格技术，发展了适于求解复杂流动的压力修正算法，采用速度的协变物理分量为动量方程的求解变量，推导得出一个适于任何曲线坐标系下控制单元体界面流通量的插值模式，可以有效地抑制不合理的波动压力场，计算表明，算法稳定收敛，计算值与试验值吻合良好，同时能预示通道内复杂的流动现象。     

激波作用下粉尘卷扬,点火与燃烧的数值研究

根据本文给出的守恒方程和计算方法，对激波与沉积可燃粉尘的相互作用进行了数值研究，讨论了激波作用下的颗粒的上扬、点火与燃烧的基本特征以及波后燃烧粉尘云的内部结构。     

三维可压缩混合层中扰动演化的研究

为了利用被动控制增强可压缩混合层中混合效率,本文用数值模拟的方法研究了三维可压缩混合层中扰动向下游的演化.通过在入口处引入不同的不稳定T-S波,给出了扰动随空间的演化过程,及马赫数对扰动演化的影响.计算结果表明在马赫数较小时,扰动波增长很快,当马赫数增大时,扰动所形成的流向涡的强度会减弱.这一结论与流动稳定性分析结果是一致的.     

漩涡发生器对燃油喷嘴的雾化性能影响的数值研究

为了探讨漩涡发生器(简称V.G)对燃油喷嘴雾化性能的影响,利用Gambit分别建立了原型喷嘴、加装2个漩涡发生器和加装4个漩涡发生器喷嘴的模型,并利用FLUENT软件进行了数值模拟,得到了三种模型的出口流场分布图。计算结果表明:在相同燃油压力下,加装漩涡发生器可以加快射流的混合速度,其中加装4个漩涡发生器的喷嘴的加速掺混能力最高,加装2个的其次,但仍比"原型"的情况有明显的改善;当燃油压力改变时,掺混能力会随着燃油压力的提高而提高。     

钝锥大攻角超声速分离流场的数值模拟

最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。