激波聚焦诱导点火和爆轰的数值研究

以二维轴对称多组分Euler方程为基础,采用非正交结构化网格和改进的波传播算法,模拟了激波在抛物形反射壁面聚焦反射诱导点火和爆轰的过程,描述了其流场形态。讨论了预混气组成、入射激波强度及反射壁面形状对点火和爆轰的影响。结果表明,激波在抛物形反射壁面顶点处聚焦反射可形成局部高温高压区域,该区域在一定条件下可点燃预混气甚至形成爆轰,其中低稀释剂浓度的预混气、较大的入射激波Mach数和较深的反射壁面有利于可燃预混气形成爆轰。     

激波与爆轰波对撞的数值模拟研究

用二阶精度NND差分格式和改进的二阶段化学反应模型模拟了爆轰波与激波的对撞过程,研究了不同强度入射激波对爆轰过渡区域的影响. 当对撞激波较弱时,透射爆轰波演变主要受流动膨胀作用的影响,可划分为对撞影响区、爆轰恢复区和稳定发展区3个阶段. 在爆轰恢复区和稳定发展区,前导激波压力经历一个过冲、然后向稳定爆轰过渡的过程,表现了爆轰波熄爆和再起爆的物理特征. 当对撞激波较强时,可燃混合气体的高热力学参数导致了更高的化学反应活化程度,形成了弱爆轰向稳定爆轰的直接转变.      

单次脉冲爆轰发动机工作过程的数值模拟

通过求解14组分和19个基元反应的CH4-O2-N2详细化学反应动力学机理的二维轴对称Navier-Stokes方程,对爆轰管内半球形高温火团引发的爆轰过程和爆轰波进入外流场后的全流场分布进行数值模拟。模拟结果显示了爆轰波在管内的成长、稳定传播、进入尾喷管后衰减为激波和进入外流场的全过程,以及爆轰管出口端附近区域的复杂涡与激波的相互作用。对轴线上的压力分布和封闭端的压力等进行了讨论,为脉冲爆轰发动机的开发研制提供参考信息。     

多管脉冲爆轰发动机管外复杂波系的研究

为了研究多管脉冲爆轰发动机外流场复杂波系变化特点,用CE/SE方法对7管脉冲爆 轰发动机内、外流场进行了计算. 7个爆轰管的排列方式是1个居中,另外6个绕其形成环 状排列,管内填充按化学当量比混合的氢气氧气混和物. 计算结果表明：CE/SE方法是一种 好的爆轰波模拟方法,能有效地捕捉激波等强间断;7管脉冲爆轰发动机外流场中存在着多 道相互作用相互影响的激波与膨胀波. 计算结果对多管脉冲爆轰发动机外流场复杂波系的特 性研究具有重要参考价值.     

多管脉冲爆轰发动机外三维流场的数值计算

为研究多管脉冲爆轰发动机黏性外流场变化特点,推导了三维黏性CE/SE方法,并对7管脉冲爆轰发动机的内、外流场进行了数值计算. 7个爆轰管的排列方式是1个居中、另外6个绕其形成环状排列. 管内流场采用一维模型进行处理; 管外流场采用三维模型进行计算, 得到多管脉冲爆轰发动机黏性外流场的变化规律. 计算结果表明,此多管脉冲爆轰发动机黏性外流场中存在旋涡及多道相互作用的激波与膨胀波, 该结果对多管脉冲爆轰发动机外流场特性研究具有一定的参考价值.      

爆轰传播研究的近代进展

综述了八九十年代国外爆轰传爆研究的主要方向和进展，包括直径效应、曲率效应、发散爆轰波的传播、爆轰冲击波动力学、爆轰传播的数值模拟和炸药爆轰性能的有关问题等。斩近分析理论和爆轰冲击波动力学的工作，对发散爆轰传播规律的认识正在不断深化，提出了一些有效的理论方法，从而为更精确的爆轰数值模拟技术，尤其是爆轰驱动等勺的精确计算开辟了新的道路。     

凝聚炸药非理想爆轰的数值模拟

本文对凝聚炸药非理想爆轰的数值模拟进行了系统介绍。对敏感、钝感炸药的理想与非理想爆轰的数值模拟方法、宏观反应速率方法及近年发展起来的ＤＳＤ（ＤｅｔｏｎａｔｉｏｎＳｈｏｃｋＤｙｎｔｉｍｉｃｓ）方法和波阵面追踪法,以及对这些方法的适应能力、存在的主要问题进行了分析与讨论。      

脉冲爆轰发动机的系统性能分析

采用简化的脉冲爆轰推进装置模型,利用热循环效率分析方法,推导了包含进气道总压恢复系数的热循环效率公式. 并在特定来流条件下,考察了一个爆轰循环中进气道总压恢复系数和燃烧室初始温度对热循环效率和比冲的影响. 研究发现,降低来流总压损失有助于提高热循环效率,而提高燃烧室初始温度能更有效地提高热循环效率. 据此,提出了多级重起爆脉冲爆轰发动机概念,利用在突扩截面上解耦的爆轰波的前导激波去再次压缩工质,进一步提高工质的热力学参数,从而提高脉冲爆轰装置的热循环效率. 推导了此种构型PDE的热循环效率计算公式,并对多级重起爆脉冲爆轰发动机进行了原理性论证. 研究结果表明,多级重起爆方法提高了燃烧室的爆前温度,从而有效地提高脉冲爆轰发动机热循环效率. 最后,关于出口工质的非完全膨胀的情况,做了定性的阐述,认为只有降低工质的出口压力,才能更有效增加工质的出口动能,从而提高热循环效率.      

脉冲爆轰发动机技术的原理和实验研究进展综述

脉冲爆轰发动机是一种新概念引擎,具有热循环效率高的显著优点。介绍了脉冲爆轰发动机的基本工作原理、优点和军用、民用上的广阔前景;综述了脉冲爆轰最新的实验研究进展;讨论了单循环爆轰的工作过程以及试验研究脉冲爆轰的主要难点。     

正庚烷云雾爆轰特性的研究

本文从理论和实验上研究了正庚烷的云雾爆轰特性。实验中采用红外技术、传感器技术、烟膜技术分别测量了正庚烷云雾引爆的延迟时间tDDT,爆速D,爆压p和爆轰波阵面结构。在理论上提出了含有化学反应项的液滴燃烧模型,进行了正庚烷云雾爆轰过程的数值模拟。计算结果和实验基本符合。     

带制退器膛口流场的三维数值模拟

利用三维任意坐标 系下Euler方程组和LU分解的数值方法模拟了不带膛口装置的某步枪膛口流场,计算结果与 实验结果比较吻合. 对某口径火炮带有制退器的膛口流场进行了数值模拟,其结果反映了制 退器内外流场的主要特征和冲击波波系之间的相互作用,在数值模拟的基础上分析 了3排侧 孔对制退器性能的影响. 本文提供了有关膛口流场的详细数据,将对制退器总体性能研究有 重要的参考意义.     

绕流直接数值模拟误差分析

对于无边界绕流问题的计算流体力学模拟通常是将物体置于“足够大”的槽道中,而通过不断改变槽道尺寸以及离散网格密度,后验对比方式来检查模拟误差。本文结合多种经典流场理论,提出一种简单的先验误差估计方法来确定槽道尺寸以及相应的网格分布。在此方法中,对于槽道尺寸的确定基于线性叠加原理（即在极小雷诺数下采用Stokes理论解叠加,而在其他雷诺数条件下采用势流理论解叠加）,来估计槽道尺寸对绕流结果的影响。而对网格尺寸与分布,则是使用多项式逼近中的基本误差分析工具,应用到速度边界层,远场势流,以及Rankine涡等简单流动,从而确定整个绕流问题中的离散误差。为了验证前面的理论分析结果,本文模拟了相当大雷诺数范围内的二维翼型以及三维圆球绕流,所得数值结果非常好地验证了理论分析。结果表明,对于Stokes流动问题,槽道尺寸需要大约100倍于物体特征尺寸来保证其结果与无边界绕流相差不超过1%;而在雷诺数超过大约100时,槽道尺寸只需10倍（二维绕流）或者5倍（三维绕流）于物体特征尺寸来达到同等精度。在此先验误差估计方法可应用于一般化的绕流问题。

     

前体涡发生器对轴对称高超声速进气道激波振荡流动的影响实验

激波振荡是高超声速进气道不起动过程中常见的流动现象,会显著降低进气道气流捕获与压缩效率、产生剧烈的非定常气动力载荷而危害飞行器安全. 从激波振荡的控制出发,实验研究了前体转捩带位置的涡发生器对轴对称高超声速进气道激波振荡流动的影响. 分别在起动和激波振荡两种进气道流态下,选择无、0.5 mm与1 mm高度涡发生器工况进行对比研究. 并采用高速纹影与壁面动态测压同步记录非定常流动特征. 结果表明,1 mm高度内的涡发生器对起动状态的进气道主流流场结构、壁面压强分布影响不显著. 但对于激波振荡流动,涡发生器会明显缩小外压缩面分离区运动范围,缩短振荡周期,提升振荡周期内壁面压强的时均值. 涡发生器的影响程度随其高度的增大而增强,其中振荡周期从无涡发生器的4 ms缩短到1 mm高度涡发生器的3.13 ms. 此外,0.5 mm高度涡发生器会使得进气道内部测点的压强振荡幅值整体下降,相比无涡发生器工况的下降幅度可达23%. 流场结构与壁面压强信号的分析表明,涡流发生器主要通过其产生的流向涡影响激波振荡流动,包含流向涡对下游边界层的扰动以及流向涡与分离区的相互干扰.      

气相爆轰波平衡胞格稳定性实验研究

本文详细介绍了研究气体爆轰的烟迹法以及采用所提出的方法得到的爆轰波通过柔性材料壁面前后平衡胞格形成、消失和重新恢复的全过程烟迹照片。实验按初压不同分五组进行。从烟迹照片中各时期的胞格尺寸、横波的变化出发,对上述过程中横波强度变化情况及柔性壁面对横波的作用和影响进行了分析,对平衡胞格恢复的条件进行了讨论。     

A DIRECT NUMERICAL SIMULATION OF CHANNEL FLOW TRANSITION UNDER AN EXTERNAL BODY FORCE

This research seeks to increase our understanding on the laminar-turbulent transition under an external body force. Direct numerical simulation by the spectral method with a weak formulation is used to solve the transient 3-D Navier-Stokes equations. Initial disturbances consist of the finite-amplitude 2-D Tollmien-Schlichting wave and two 3-D oblique waves. Competitions among different modes were computed during transition for different Richardson numbers.

It is found that the body force can modify the transition mechanism of flows between two vertical plates. The body force was found to hasten the formation of three-dimensional flow. A non-laminar flow induced by the body force may present when the background flow is still laminar.     

随机选取法计算化学反应流爆轰过程引起的参数涨落

利用ＲＣＭ（随机选取法）解反应流爆轰问题时,掌握影响流动参数随机涨落的主要因素是选择适当计算参数、保证得到正确结果的关键。各种流动参数的随机涨落相互关联,只需要搞清楚炸药未反应质量分数的随机涨落并做出分析,其它流动参数涨落的情况也就清楚了。给出了影响未反应质量分数涨落的因数,同时还给出了具体算例,检验所得结论的正确性。     

高压气体冲击流场的数值研究

本文针对高压气体流动产生强冲击波的特点,利用近年来公认较好的高精度、高分辨率总变差不增(TVD)格式对圆柱形高压区气体向大气膨胀过程中产生的冲击流场问题进行了数值模拟。数值试验表明:采用高分辨率的TVD差分格式能够很好地描述高压气体冲击流场的波系结构,且没有出非物理振荡,对激波的分辨率高。     

外加磁场对熔融半导体CZ法拉晶过程中流动场影响的数值计算

数值分析结果表明,外加磁场可以改变熔融半导体中的流型,几千高斯的磁场可以显著地减小熔体的流动,但对温度场影响不大。     

某水下航行体尾翼展开过程中流场的数值模拟

利用Pro／E软件对水下航行体进行了实体建模．使用Fluent 6．0软件,应用三维黏性定常不可压缩Navier-Stokes方程,k-ε两方程湍流模型,采用有限体积法对水下航行体尾翼展开过程中的流场和载荷进行了数值仿真计算,根据所得结果对水下航行体的展开规律进行了分析计算．对照实验结果,说明了该数值解法的正确性,同时表明了数值分析结果对解决工程问题具有实际指导意义,可应用于水下航行体折叠尾翼的设计和改进．