稠密颗粒床中点火传火过程的两维两相流数值模拟

本文结合火炮工程背景,建立伴随化学反应的稠密颗粒群气固两相非定常两维轴对称流和点火器内的非定常一维两相流数学模型,并考虑点火器对主装药床的耦合作用。采用MacCormack两步差分格式,模拟点火传火过程中火焰传播及其物理量沿轴向和径向变化规律。部分计算结果与实验结果相吻合。发现点火开始阶段存在明显的径向效应,当火焰波传到固壁后,轴向流占主要优势。     

高含能颗粒床中瞬变点传火过程的三维两相流数值模拟

结合火炮工程背景，建立高含能稠密颗粒床中非定常三维两相流与点火器中非定常一维两相流的数理模型，并考虑主装药床与点火器的相互耦合作用。采用分数步法（ＦＳＭ）数值求解。发现在点传火阶段，膛内存在明显的三维效应，并伴随波动现象，火药颗粒向壁面聚集。部分计算结果与试验结果相一致。     

片状随行装药的两相流数值模拟

给出了片状随行装药结构的一维两相流数学物理模型,并对某５７ｍｍ随行装药火炮进行了数值模拟。研究了片状随行装药火炮的膛压、压力波、气固相速度以及空隙率等参量的分布情况。此外,还对比研究了计算结果与实验结果之间的差异。结果表明,该模型的计算骅音符得较好。     

大长径比点火管高密实火药床点传火过程两相流的数值模拟

设计了某中心点火管,完成了该点火管的点传火实验,针对该点火管长径比大、装填密度高的特点,建立了点火管内气固两相流动和燃烧过程的一维两相流模型,并进行了数值模拟。计算结果与实验结果良好符合,说明计算模型能够准确描述点火管内的实际物理化学过程,计算程序参数取值合理,该计算程序可为此类点火管各种结构尺寸及装填条件下的点传火性能分析及优化计算提供充分的理论依据和方法。并且,根据计算结果初步分析了该结构及装填条件下点火管的点传火性能,为下阶段工程优化设计提供参考     

大长径比、高密实火药床点火管点传火特性的影响因素

利用计算程序对影响大长径比、高密实火药床点火管点火性能的主要因素进行了对比模拟,讨论了不同结构尺寸和装填条件下的点传火性能,通过对计算结果的分析,总结了结构参数和装填条件对大长径比、高密实火药床点火管点传火性能的影响规律：小孔直径、首孔高度、单位长度小孔面积及装填密度显著影响点火压力及点传火性能;小孔直径、首孔高度及单位长度小孔面积是影响破膜后泄压速度的主要因素,而装填密度会影响药床的透气性及火焰传播的通畅性,均对点火安全性、瞬时性及一致性产生重要影响作用。

     

铝粉/空气二维黏性两相爆轰的数值模拟

为了深入研究爆轰波形成和传播的机理与特性,建立了管内铝粉/空气二维黏性两相爆轰过程的数学模型,采用守恒元与求解元方法进行数值计算,并对其物理参数的分布进行了分析。结果表明:管内燃烧转爆轰的初期,压力沿径向变化明显,与壁面碰撞有明显的反射波;燃烧转爆轰的中后期,压力沿轴向变化明显,但径向效应仍不能忽视,碰撞形成的反射波对最终稳定爆轰波的形成影响较大。研究结果同时表明:不仅铝粉颗粒初始半径对爆轰波的形成与传播有一定的影响;而且气体的黏性作用在研究爆轰管内近壁面处流场时不容忽视。研究结果有利于进一步揭示铝粉燃烧转爆轰的机理。     

液体发射药电热化学发射内弹道一维两相流模型及数值模拟

在理论分析的基础上 ,结合实验条件建立了液体发射药电热化学发射内弹道一维两相流数学模型 ,该模型包含电弧等离子体生成和运动、等离子体卷吸液体发射药生成小液滴以及膛内多相流体的相互作用。利用该模型对电热化学发射内弹道过程进行了数值模拟 ,并将模拟结果与实验数据进行了分析比较 ,结果表明模拟结果与实验数据基本吻合。     

两相正激波松弛流的数值研究

对气、固两相正激波松弛流动进行了数值研究,方程中考虑了压力梯度及非匀速运动引起的附加质量力。给出了完整的计算方法,对三种不同固相容积比、不同固相颗粒直径和不同马赫数分别进行了数值计算,所得计算结果给出了气、固两相的速度、温度、压力、密度、固相容积比及熵值沿流向分布情况。文章对气、固两相参数在松弛区中的变化规律进行了讨论。     

气液两相流泄漏扩散的数值模拟

采用Eulerian-Lagrangian方法数值模拟过热液体丙烷从喷嘴释放后射流的发展过程。其中,对丙烷蒸汽与空气的气相混合物采用欧拉法求解,利用拉格朗日法跟踪丙烷液滴的运动轨迹,并加入颗粒随机轨道模型来考虑液滴的湍流扩散效应,数值计算所得的物理参量与实验测量结果相吻合。在此基础上,本文将此数值模型和方法应用于开放空间、常温高压下储存液化气(丙烷)的容器发生小孔泄漏时,两相射流在大气中的扩散过程,计算结果符合实际情况。文中还讨论了影响两相流扩散过程的主要因素,为预测可能发生的火灾爆炸事故提供理论依据。     

节水水嘴起泡器气液两相流数值模拟

为了研究节水水嘴起泡器内部两相流的流动规律,采用Fluent软件对其内部流场进行数值模拟.根据起泡器内部流场的流动特性,采用欧拉两相流模型以及RNG(re—normalization group)κ-ε湍流模型,分析起泡器出口截面气液两相体积分数和速度的分布特点.结果表明:增大入口水流速度可以加快分散出口截面气液两相的分布,缩短流体流动的稳定时间;整流网具有分散流体,降低流速的作用;错开整流网相邻层之间的网格可以改善出口截面的液相分布;本模型中整流网采用三层网格达到较好的出水效果.     

点火管破孔的随机过程的一维两相流计算

从宏观的角度考虑了破孔过程中的随机因素,提出了破孔压差的概念,把各个传火孔的破孔压差看作是服从正态分布的独立随机变量,建立了该过程的随机模型,并用Monte Carlo方法进行了随机模拟。所得结论对点火系统的设计安全性具有一定的指导意义。     

悬浮RDX炸药和铝颗粒混合粉尘爆轰的数值模拟

采用两相流方法对炸药颗粒直径为20.0 m时与铝颗粒混合物的爆轰波的发展与传播过程及爆轰波参数进行了数值计算。结果表明,在炸药粉尘中加入铝颗粒,可以大大提高爆轰波参数。当铝颗粒直径为3.4 m时,尽管铝颗粒的直径较炸药颗粒直径小,但由于炸药颗粒的点火温度低,二者的点火时间相差不多。如果铝颗粒的直径为7.0 m,由于铝颗粒的点火滞后于炸药颗粒的点火,混合颗粒粉尘中可能形成双波阵面的爆轰波。     

基于详细化学反应动力学的激光点火模型

为研究含能材料在激光作用下的点火特性,从基本守恒方程组出发,基于详细化学反应动力学模型,建立了含能材料激光点火的气相模型.利用所建立的模型对RDX在激光辐照下的点火过程进行了数值模拟,得到了点火过程中的瞬时温度分布和组分分布,对点火特性进行了分析.计算得到的不同激光功率密度下的RDX点火延迟时间与文献结果一致.     

考虑两相流内弹道的自行火炮发射动力学计算

为了更加准确地计算火炮弹丸的起始扰动,将两相流内弹道应用于发射动力学研究。建立了完整的自行火炮系统发射动力学方程组,包括火炮系统的体动力学方程组、弹丸在膛内运动的动力学方程以及两相流内弹道方程。编制了计算程序,实现了对某自行火炮发射过程的数值模拟,在准确计算内弹道过程的同时,获得了火炮的动力响应、弹丸膛内运动和起始扰动。部分模拟结果与实验实测结果吻合较好。计算表明,采用两相流内弹道模型将提高发射动力学计算精度。     

Numerical simulation and homogenization of two-phase flow in heterogeneous porous media

A mathematically rigorous method of homogenization is presented and used to analyze the equivalent behavior of transient flow of two incompressible fluids through heterogeneous media. Asymptotic expansions and H-convergence lead to the definition of a global or effective model of an equivalent homogeneous reservoir. Numerical computations to obtain the homogenized coefficients of the entire reservoir have been carried out via a finite element method. Numerical experiments involving the simulation of incompressible two-phase flow have been performed for each heterogeneous medium and for the homogenized medium as well as for other averaging methods. The results of the simulations are compared in terms of the transient saturation contours, production curves, and pressure distributions. Results obtained from the simulations with the homogenization method presented show good agreement with the heterogeneous simulations.     

流动系统中爆轰波传播特性的数值模拟

基于带化学反应的二维Euler方程,采用氢气/空气的9组分19步基元反应简化模型,对充有当量 比的氢气/空气预混气体的矩形爆轰流场中爆轰波的传播过程进行了数值模拟,讨论了均匀来流对爆轰波传 播的影响。数值结果表明,在均匀来流的影响下,上游方向上的燃烧强度大于C-J爆速,下游方向上的燃烧强 度小于C-J爆速;上游方向传播的爆轰波的阵面压力大于下游方向传播的爆轰波的阵面压力。所以,经典的 C-J爆轰理论并不适用于流动系统中爆轰波传播特性的研究。     

Hyperbolicity and one-dimensional waves in compressible two-phase flow models

Hyperbolic models for compressible two-phase flows including a conservative symmetric hyperbolic model are reviewed. The basis for a theory of shock waves is developed within the framework of the latter. The analysis of small amplitude discontinuities allows us to conclude that in general there are two types of shocks corresponding to two sound waves. The problem of transition between a pure phase and a mixture (the phase vacuum problem) is analysed. It is proved that for some models the smooth centred wave solution can not provide such a transition. Within the framework of our conservative model there is the possibility of constructing discontinuous solutions which can resolve the phase vacuum problem.PACS: 47.55Kf, 47.40.-xE. Romenski: On leave from Sobolev Institute of Mathematics, Russian Academy of Sciences, Novosibirsk 630090, Russia