气体体积粘性对二维环形激波聚焦的影响

基于体积粘性系数ζ的分子运动论和连续介质理论,对二维环形激波聚焦（马赫数Ma=2.0）的体积粘性效应进行数值研究,结果表明：对于热完全气体,体积粘性使得激波汇聚中心点处的压力减小、温度增加、密度减小,聚焦点物理参数的改变量分别可达20%、10%、30%,体积粘性效应对环形激波聚焦的影响是不可忽略的;与转动模态相比,在振动模态下环形激波聚焦的体积粘性效应更为明显,因为激波聚焦点附近的体积粘性应力ζ▽·V与热力学压力p达到同一数量级,从而显著改变了流动参数.     

长试验时间激波风洞测力技术研究

中国科学院力学研究所复现飞行条件高超声速激波风洞JF12的落成突破了毫秒级试验时间的瓶颈,有效试验时间超过100ms.因此,对于JF12长试验时间激波风洞的测力试验,基于应变天平技术较为成熟、结构简单等优点,我们考虑采用传统的应变计天平.但是,激波风洞来流冲击所带来的惯性力干扰导致天平测力系统产生低频振动,传统内置应变天平的结构刚度很难保证信号有足够的处理周期,这大大限制了激波风洞测力模型的尺寸和重量.针对这个难题,基于JF12激波风洞的运行特点及对测力天平刚度特性的特殊要求,优化设计了应变天平的测力单元结构以适用于这种脉冲动态测力试验,相应加工制造了大刚度、低干扰、高灵敏度的系列脉冲型应变天平,结构形式包含了杆式和盒式,最大载荷(法向力)从1kN到30kN,以满足不同尺度飞行器的测力试验需求.同时,我们应用不同尺度的测力模型对研制的脉冲型天平在JF12激波风洞进行了一系列动态气动力测量试验,以进一步评估JF12系列脉冲型应变天平的结构特性和测力性能.      

爆轰波与激波对撞爆轰状态演化机理研究

对当量比氢氧混合气体中爆轰波与激波的正面对撞过程进行了二维数值研究. 采用了二阶精度NND差分格式与改进的两阶段化学反应模型,并以数值x-t纹影图以及烟迹图记录了对撞过程. 数值研究表明,透射爆轰波受到膨胀影响首先会衰减,甚至产生局部解耦现象;然后由于三波点的碰撞又能再次重新耦合. 在爆轰波对撞过程中,由于燃烧不均匀性而产生的弱横波对爆轰胞格的形成起着重要作用.      

一种基于AUSM思想的通量分裂方法

根据对流迎风分裂(AUSM)思想提出一种通量分裂方法,称为K-CUSP格式.它与传统H-CUSP和E-CUSP格式的最大差异在于总能量的分裂:K-CUSP格式将无粘守恒通量中所有的运动学量分裂到对流项,所有的热力学量分裂到压力项,即总能量被分裂成动能和静焓.对于压力项的数值通量,采用一种新的界面构造方法.数值测试表明:①K-CUSP格式继承了FVS格式的简单性和稳健性.在激波后不易出现压力过冲,在膨胀区域没有振荡,优于AUSM和WPS格式;②K-CUSP格式继承了FDS格式的分辨率.激波间断的分辨率和H-CUSP、E-CUSP格式基本相同,接触间断的分辨率高于FVS格式,低于Roe、AUSM和WPS格式.AUSM和WPS格式在计算运动接触间断时,速度存在很大振荡,而新格式不存在振荡.     

钝锥模型高焓流动的再压缩效应研究

在钝体高焓流动中,热化学非平衡效应使膨胀区壁面压力显著低于冻结流动,相关研究表明,再压缩效应可以消除上述差异. 为详细分析此问题,数值求解了Euler 方程,化学反应源项采用有限速率模型并考虑了5个组分,17 个基元反应. 结果表明,随着再压缩角的增大,热化学非平衡效应引起锥面压力先降后增非单调变化,这是由于再压缩流场中放热复合反应与吸热分解反应的强度变化决定的,其宏观表现为比热比的空间不均衡分布.      

散心柱面胞格爆轰演化数值研究

采用有限体积方法,在自适应非结构网格上求解二维含化学反应Euler方程,数值研究了柱面胞格爆轰波演化现象．化学反应计算采用单步可逆总包反应模型．数值结果演示了散心柱面胞格爆轰波演化过程中胞格结构的分裂现象,获得了与实验结果定性一致的结果．胞格结构的分裂演化在点火区近场和远场显示了不同的特点,其中爆轰波传播过程中波阵面当地曲率的变化是控制胞格分裂演化行为的关键因素．数值结果也显示胞格结构的分裂现象来自于爆轰波前锋结构中横波的自组织行为,即沿爆轰波波面传播的小扰动发展成为横波的过程,这种现象与胞格爆轰波的不稳定性密切相关．     

Numerical investigation on evolution of cylindrical cellular detonation

Cylindrical cellular detonation is numerically investigated by solving two- dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.     

频散可控格式的一种推广形式及其在爆轰波马赫反射中应用

频散可控格式是高精度捕捉激波的新格式.对原有频散可控格式(DCD)进行了推广,给出了适用于网格点排列不规则情况下的DCD格式,使其具有更广泛的适应性,在此基础上构造了在非结构三角形网格下的DCD格式.用一些典型算例对推广后的DCD格式进行检验,结果表明该格式是合理和可靠的,并能够较好地保留原有DCD格式的二阶精度.把推广后的DCD格式和具有8种组分和20个化学反应的基元反应模型相结合,对氢氧爆轰波在直管道中的传播问题及楔面上马赫反射问题进行计算,计算结果和实验结果比较表明,在非结构三角形网格下的DCD格式能够有效地捕捉爆轰波,在爆轰波阵面上不会产生振荡或是抹平间断现象.     

一种易于并行求解Euler方程组的分区技术

在具有复杂边界的计算区域内，求解偏微分方程组时，经常需要分区和并行计算，分区方法直接关系到数值计算的并行化程度，本文在应用时间算子分裂方法求解Euler方程组的过程中，提出了一种非常容易实现并行化计算的分区技术。     

氢氧爆轰波在变截面扩张管道中传播的数值模拟

采用二阶精度频散控制耗散格式(DCD)和8组分20个方程的基元反应模型,对轴对称变截面管道中氢氧爆轰波传播进行数值模拟。结果表明,爆轰波传播至突变截面扩张管道时,由于稀疏波的作用可能会使爆轰波局部熄爆甚至完全熄爆,对于某些敏感度高的反应气体爆轰波可以二次起爆。而在渐变截面扩张管道爆轰波相对不易熄爆。