水下超声速过膨胀燃气射流的流场特性

为研究水下超声速过膨胀燃气射流的流场特性,在压力水筒中开展了大扩张比锥形喷管的固体火箭发动机水下点火实验,并基于雷诺时均Navier-Stokes（RANS）方法和流体体积（VOF）模型进行数值求解,分析了过膨胀燃气射流与水介质的相互作用过程。研究表明：超声速过膨胀燃气建立射流通道后,射流核心区长度随喷管落压比的减少而减少;射流核心区剧烈振荡,表现为高频的膨胀和收缩,振荡频率随喷管落压比的减小而增加,范围为100~200 Hz;射流边界不断振荡,并伴随波系结构变化,当过膨胀程度较大时,激波进入喷管使其发生流动分离现象,流动分离点周期性往复移动;分离区内压力脉动没有显著的特征频率,主要集中在100~600 Hz的宽频带,锥形喷管水下流动分离的简易判据为喷管出口压力不低于环境背压的0.44倍。     

推进剂药浆混合均匀性研究

用Ｈａａｋｅ流变仪研究了高固体含量丁羟推进剂的混合均匀性,提出了用复合粘度和复合模量作混合均匀性表征参数的评价方法。研究结果表明,高固体含量丁羟三组元配方的预混阶段需１５ｍｉｎ～２０ｍｉｎ快速混合,总预混时间以３０ｍｉｎ～３５ｍｉｎ为宜。     

非线性规划的SUMT混合法及其应用

本文提出了非线性规划的内点法与外点法相结合的一种新方法－ＳＵＭＴ混合法,它既兼备内,外点法的功能,又克服了各自单一方法中的缺点。另外,本文还给出无约束问题两种解法的比较及证明,并应用上述方法,利用计算机对精密齿轮按最小转动惯量进行最优设计,取得了令人满意的成果。     

时间发展可压化学反应混合层直接数值模拟

基于一个简化的单步不可逆燃烧动力学模型,对时间发展二维可压化学反应平面混合流场,实现了高精度迎风紧致有限差分直接数值模拟。研究了化学反应放热与流体混合运动间的耦合作用,分析了前者对可压混合层厚度、拟序结构及产物生成率等方面的影响。     

模拟高度复杂流动的自适应网格算法

提出一种基于对等弧分布律进行扭转修正的自适应网格算法。该方法将三维网格系统分裂为３个单方向自适应网格子系统,在计算中与主控方程的求解格式相耦合,交替地使用。应用所构造的算法结合ＭａｃＣｏｒｍａｃｋ显式格式对喷流／绕流干扰流场进行了雷诺平均Ｎ－Ｓ方程模拟,给出了清晰的剪切层与激波波系结构。计算结果表明,该自适应网格算法可较大程度地提高ＭａｃＣｏｒｍａｃｋ格式的分辨率。     

关于建立高阶差分格式的问题

为了能在不太密的网格上捕捉到流场的细致结构。通常采用高阶精度的差分 式进行数值模拟。为能抑制计算中和在激波附近产生的虚假波动,本文从物理构思出发,提出了建立高阶格式的两个基本原则,作者称之为抑制波动的原则和稳定性原则。     

不同流速比超声速混合层气动光学效应研究

针对不同流速比混合层流场气动光学效应问题,首先采用大涡模拟数值方法进行了数值仿真,其次用光线追迹方法进行了气动光学效应仿真,最后对混合层的气动光学效应进行了评价。结果表明,混合层流场涡结构与光程差极小值存在一一对应关系,并且其流速比低的混合层气动光学效应更加严重。同时,研究了不同光学参数对斯特列尔比的影响,得到了相应的规律。     

不同壁温及差分格式下超燃冲压发动机的仿真

针对圆截面超燃冲压发动机相关实验进行数值仿真,采用了不同的壁面温度条件和差分格式,并探究其对仿真结果的影响规律。结果表明采用常用的绝热壁面边界条件会导致激波边界层交互位置处的回流区变大,进而导致流道收缩,壁面压力升高。采用高精度格式不仅不能更为精确地计算回流区的大小,反而会使回流区面积变得更大。通过改变壁面边界条件,采用300 K的等温壁面,可以使激波边界层交互位置处的回流区变得极小,使仿真结果与实验结果十分吻合。计算结果表明,针对当前构型与实验条件,盲目提高数值格式的精度并不一定会得到更好的结果,相反可能会使仿真与实验更加偏离,而适当修改壁面温度条件,即使数值格式精度较低,仍然可以得到很好的仿真结果。最后,针对与实验结果吻合的仿真结果,分析了圆截面超燃冲压发动机的流场特性,重点研究了凹腔剪切层及其质量交换特性。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

无隔道超声速进气道附面层排除特性飞行试验研究

为了分析评估某型歼击机无隔道进气道附面层的排除特性,设计搭建鼓包表面附面层压力梯度测量试验系统,进行了不同飞行高度、马赫数和姿态角等工况下的飞行试验。通过对飞行试验数据的整理、计算和对比分析同型号的缩比模型风洞试验结果,研究了无隔道进气道鼓包表面附面层排除特性。研究结果表明：稳定平飞时,在亚音速范围内,随着飞行高度的增加,鼓包构型对附面层的排除效果增大,而在超音速范围内,变化规律相反;在接近马赫数1.8及以上飞行工况下,鼓包表面附面层的扫除能力有所减弱,附面层气流分离加速,进而会造成较大的进气压力损失和畸变。单纯迎角飞行有利于增强附面层的排除能力;而带侧滑角飞行时,附面层压力系数曲线的拐点沿鼓包中心线平行向“背风面”偏移,偏移量与侧滑角成正比,进气道鼓包表面“迎风面”附面层排除能力增大,而“背风面”受气流分离影响而减弱。