2021-08期封面

铝粉反应模型是对悬浮铝粉尘气-固两相爆轰进行数值模拟研究的关键。通过考虑铝粉燃烧产物氧化铝（Al₂O₃）在高温下的分解吸热反应,改进了铝粉的扩散燃烧模型。将该模型嵌入到三维的气-固两相爆轰数值计算程序中,分别对铝粉/空气混合物以及铝粉/氧气混合物的爆轰进行了数值模拟,计算得到的稳定爆轰波速度与实验结果、文献值均吻合较好,误差小于5.5%,表明改进的铝粉反应模型适用于不同氧化气体氛围中铝粉尘爆轰的模拟计算。此外,对两相爆轰参数及爆轰流场的物理量分布进行分析,获得了铝粉反应模型对爆轰波结构的影响规律。     

惰性颗粒抑爆过程的数值模拟

对高温火球诱导的爆炸及惰性颗粒的抑爆过程进行了数值研究。由于扰动特征时间 ,两相弛豫特征时间和化学平衡特征时间的差异 ,这种带激波的两相化学反应流守恒方程具有刚性。利用处理包括组元守恒方程的全耦合TVD格式和可处理方程源项的Lax Wendroff Rubin格式分别求解气相和颗粒相 ,并根据各特征时间的比值用分步法处理方程的刚性。计算结果反映了均温燃烧火球在可燃环境中形成加速火焰和激波 ,进而成长为爆轰波的过程 ,以及惰性颗粒对爆轰波的抑制过程 ,反映了激波、化学反应和惰性颗粒之间的相互作用和耦合。计算结果表明 ,仅当颗粒浓度大于某值时才能有效抑制爆轰 ,否则无论多大的颗粒散布区域 ,皆不能使爆轰波完全抑制 ,爆轰波或是在其中发展为两相爆轰波 ,或是穿越该区域后重新成长为气相爆轰。     

爆轰波平掠惰性气体界面及其解耦现象的数值研究

采用基元化学反应模型和迎风TVD格式,数值研究了爆轰波平掠惰性气体界面时的物理现象及其作用机制,并用点隐方法克服化学反应源项引起的计算刚性。数值结果显示,当爆轰波平掠过惰性气体界面时,形成了爆轰波、界面、透射激波以及稀疏波相互作用的现象。在高N2比例稀释的可燃混合气体情况下,当爆轰波平掠过特定惰性气体界面时,它与惰性气体界面相互作用产生的稀疏波可以导致爆轰波的解耦。     

气相爆轰波正反射激波加速研究

实验采用压力传感器测量了指定点压力时间曲线。数值模拟基于二维反应欧拉方程和基元反应模型,采用二阶附加半隐的龙格-库塔法和5阶WENO格式分别离散时间和空间导数项,获得了指定点数值压力时间曲线。理论分析基于爆轰理论和激波动力学,分析了气相爆轰波反射过程所涉及的复杂波系演变并获得了反射激波速度。结果表明:本文数值模拟和理论计算定性上重复并解释了实验现象。气相爆轰波在右壁面反射后,右行稀疏波加速反射激波。其加速原因是:尽管激波波前声速减小,但激波马赫数增大,波前气流速度减小。在低初压下,可能还由于爆轰波后未反应或部分反应气体的作用,导致反射激波加速幅度比高初压下大。     

气体-燃料液滴两相系统爆轰的数值模拟

用两相流体力学模型对气体 燃料液滴系统进行了研究。数值模拟了点火后两相系统爆轰波的发展过程，得到爆轰波的结构和参数。数值模拟结果表明气体 燃料液滴系统爆轰波有较宽的反应区，因而两相爆轰波的曲率对爆速的影响效应十分明显。进行了燃料液滴尺寸对爆轰波的结构和参数的影响的数值模拟。除了很小的液滴外，燃料液滴在爆轰波前导激波面和ＣＪ面间不能完全气化。随着液滴尺寸的增加，燃料液滴在爆轰波前导激波面和ＣＪ面间释放出的能量随之减少，爆轰参数也随之下降。     

Modeling spherical explosions with aluminized energetic materials

This paper deals with the numerical solution and validation of a reactive flow model dedicated to the study of spherical explosions with an aluminized energetic material. Situations related to air blast as well as underwater explosions are examined. Such situations involve multiscale phenomena associated with the detonation reaction zone, the aluminium reaction zone, the shock propagation distance and the bubble oscillation period. A detonation tracking method is developed in order to avoid the detonation structure computation. An ALE formulation is combined to the detonation tracking method in order to solve the material interface between detonation products and the environment as well as shock propagation. The model and the algorithm are then validated over a wide range of spherical explosions involving several types of explosives, both in air and liquid water environment. Large-scale experiments have been done in order to determine the blast wave effects with explosive compositions of variable aluminium content. In all situations the agreement between computed and experimental results is very good.     

Intensification of a detonation wave in the zone of secondary reactions in a two-phase medium

A method is proposed for the numerical calculation of one-dimensional nonsteady-state flows of a mixture of a gas with particles, based on the separation of a system of differential equations for a two-phase medium into two subsystems. The problem is solved concerning the propagation of a plane detonation wave in a mixture of a detonating gas with particles, behind the front of which secondary chemical reactions are taking place between the vapors of the particle material and the detonation products. The velocity profiles of the gas and of the thermodynamic functions behind the detonation wave front are determined, and also the dependence of the detonation velocity on the distance to the point of initiation. The conditions for intensification of the detonation wave are obtained in the zone of secondary reactions.Leningrad. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 92–96, September–October, 1972.     

Modelling detonation waves in condensed energetic materials: multiphase CJ conditions and multidimensional computations

A hyperbolic multiphase flow model with a single pressure and a single velocity but several temperatures is proposed to deal with the detonation dynamics of condensed energetic materials. Temperature non-equilibrium effects are mandatory in order to deal with wave propagation (shocks, detonations) in heterogeneous mixtures. The model is obtained as the asymptotic limit of a total non-equilibrium multiphase flow model in the limit of stiff mechanical relaxation only (Kapila et al. in Phys Fluids 13:3002–3024, 2001). Special attention is given to mass transfer modelling, that is obtained on the basis of entropy production analysis in each phase and in the system (Saurel et al. in J Fluid Mech 607:313–350, 2008). With the help of the shock relations given in Saurel et al. (Shock Waves 16:209–232, 2007) the model is closed and provides a generalized ZND formulation for condensed energetic materials. In particular, generalized CJ conditions are obtained. They are based on a balance between the chemical reaction energy release and internal heat exchanges among phases. Moreover, the sound speed that appears at sonic surface corresponds to the one of Wood (A textbook of sound, G. Bell and Sons LTD, London, 1930) that presents a non-monotonic behaviour versus volume fraction. Therefore, non-conventional reaction zone structure is observed. When heat exchanges are absent, the conventional ZND model with conventional CJ conditions is recovered. When heat exchanges are involved interesting features are observed. The flow behaviour presents similarities with non ideal detonations (Wood and Kirkwood in J Chem Phys 22:1920–1924, 1950) and pathological detonations (Von Neuman in Theory of detonation waves, 1942; Guenoche et al. in AIAA Prog Astron Aeronaut 75: 387–407, 1981). It also present non-conventional behaviour with detonation velocity eventually greater than the CJ one. Multidimensional resolution of the corresponding model is then addressed. This poses serious difficulties related to the presence of material interfaces and shock propagation in multiphase mixtures. The first issue is solved by an extension of the method derived in Saurel et al. (J Comput Phys 228(5):1678–1712, 2009) in the presence of heat and mass transfers. The second issue poses the difficult mathematical question of numerical approximation of non-conservative systems in the presence of shocks associated to the physical question of energy partition among phases for a multiphase shock. A novel approach is used, based on extra evolution equations used to retain the information of the material initial state. This method insures convergence in the post-shock state. Thanks to these various theoretical and numerical ingredients, one-dimensional and multidimensional unsteady detonation waves computations are done, eventually in the presence of material interfaces. Convergence of the numerical hyperbolic solver against ZND multiphase solution is reached. Material interfaces, shocks, detonations are solved with a unified formulation where the same equations are solved everywhere with the same numerical scheme.     

固体炸药爆轰与惰性介质相互作用的一种扩散界面模型

提出一种保持热力学一致性的扩散界面模型，用来数值模拟固体炸药爆轰与惰性介质的相互作用问题。基于混合网格内各组分物质间可以达到力学平衡状态而不能达到热学平衡状态的假设，由混合网格能量守恒以及压力相等条件，推导出每种组分物质的体积分数演化方程。由此获得的扩散界面模型包括组分物质的质量守恒方程、混合物质的动量及总能量守恒方程，同时包括组分物质的体积分数演化方程和混合物质的压力演化方程。该扩散界面模型的主要特点是考虑了化学反应以及热学非平衡的影响。提出的扩散界面模型在物质界面附近不会出现物理量的非物理振荡现象、适用于任意表达形式的物质状态方程以及任意数目的惰性介质。     

可燃介质中飞行圆球诱导斜爆轰的流场结构

基于带化学反应的二维轴对称Euler方程,利用带有Superbee限制函数的波传播算法,对氢/氧/ 氮预混气中飞行圆球诱导斜爆轰进行了数值模拟。结果表明,在达姆科勒数Da略大于临界情形时,圆球诱导 的驻定爆轰是一个由强过驱斜爆轰、弱过驱斜爆轰、反应激波和惰性激波组合而成的复合结构。波后的圆球 绕流流场内存在2个亚音速区和1个超音速区。在圆球背风表面,还形成了第2道斜激波。     

Distinctive features of galloping detonation in a supersonic combustible-mixture flow under an inert gas layer

The problem of detonation initiation in a supersonic flow of a stoichiometric propane-air mixture occupying partially or completely the cross-section of a plane channel is considered. The initiation in the flow is produced by a step or a wall completely cutting off the flow. The study is conducted within the framework of one-stage combustion kinetics. A numerical method based on the Godunov scheme is employed. The critical conditions for detonation formation are determined in terms of the oncoming flow velocity. A previously unknown mechanism of detonation propagation is found; it is related with the presence of the combustible mixture in the wall layer under an inert gas layer. It is due to the formation of a complicated wave structure of the flow characterized by the penetration of a shock wave formed in the inert gas layer into a combustible mixture layer ahead of the detonation wave with the result that the latter layer is heated and ignited. The process as a whole is periodic in nature, as distinct from the conventional cellular detonation in a homogeneous fluid. Many problems arise in connection with the use of detonation in engines and other power plants. The most important among them are detonation excitation and stabilization in combustion chambers. The detonation initiation within a layer under conditions of unbounded space and a fluid at rest was experimentally investigated in [1]. In the case of a combustion chamber bounded in the transverse direction, some new effects accompanying the detonation might be expected.     

基于最小二乘法的反应流拉格朗日分析在冲击起爆中的研究

现有反应流拉格朗日分析方法在已知粒子速度情况下,求解方法仍有不足。本文针对这一情况,将反解法与自洽检验法相结合,提出了基于最小二乘法的反应流拉格朗日反解法,该方法的理论精度能够实现应力沿路径线的 M阶导数为零（M为迹线数）,并且该方法能够满足自洽检验法。为了验证该方法的有效性,利用该方法对一组小隔板冲击起爆试验数据进行了处理,对比了本文方法、试验处理以及传统的反解法处理结果,表明该方法不仅可以适当减小迹线上的偶然误差,还能够使得迹线函数更好地反映各物理量沿迹线的变化性态。     

Numerical prediction of droplet dynamics in turbulent flow,using the level set method

In the present article, the droplet dynamics in turbulent flow is numerically predicted. The modelling is based on an interfacial marker-level set (IMLS) method, coupled with the Reynolds-averaged Navier–Stokes (RANS) equations to predict the dynamics of turbulent two-phase flow. The governing equations for time-dependent, two-dimensional and incompressible two-phase flow are described in both phases and solved separately using a control volume approach on structured cell-centred collocated grids. The topological changes of the interface are predicted by applying the level set approach. The kinematic and dynamic conditions on the interface separating the two phases are satisfied. The numerical method proposed is validated against a well-known computational fluid dynamics problem. Further, the deformation and breakup of a single droplet either suddenly moved in air or exposed to turbulent stream are numerically investigated. In general, the developed numerical method demonstrates remarkable capability in predicting the characteristics of complex turbulent two-phase flows.     

Numerical predictions of oblique detonation stability boundaries

Oblique detonation stability was studied by numerically integrating the two-dimensional, one-step reactive Euler equations in a generalized, curvilinear coordinate system. The integration was accomplished using the Roe scheme combined with fractional stepping; nonlinear flux limiting was used to prevent unphysical solution oscillations near discontinuities. The method was verified on one- and two-dimensional flows with exact solutions, and its ability to correctly predict one-dimensional detonation stability boundaries was demonstrated. The behavior of straight oblique detonations attached to curved walls was then considered. Using the exact, steady oblique detonation solution as an initial condition, the numerical simulation predicted both steady and unsteady oblique detonation solutions when a detonation parameter known as the normal overdrive, , was varied. For a standard test case with a specific heat ratio of , a dimensionless activation energy of , and dimensionless heat release of , an oblique detonation with a constant dimensionless component of velocity tangent to the lead shock, , underwent transition to unstable behavior at . This is slightly higher than the threshold of predicted by one-dimensional theory; thus, the two-dimensionality renders the flow slightly more susceptible to instability. Received 4 August 1996 / Accepted 5 March 1996     

钝感炸药的超压爆轰与冲击起爆过程数值模拟

采用Hybrid反应率结合修正的JWL方程,研究了LX-17、超细TATB等钝感炸药的冲击起爆（SDT）过程,并计算了爆轰波的对碰现象。结果表明,该方法计算钝感炸药的冲击起爆过程与实验数据符合较好;计算爆轰波对碰区的峰值压力提高了10％。     

颗粒材料中致密波结构研究

采用一维两相流模型与相应颗粒构形应力函数,研究了致密波的形成及其结构．用简化两相流模型系统地讨论致密波对有关因素的依赖关系．分析指出：小于基体材料音速的致密波仅能在非理想颗粒材料中存在,从波前到波后,所有状态物理量光滑过渡．大于基体材料音速的致密波,波头可能存在间断．应力函数与致密粘性确定后,致密波速度决定致密波结构、宽度、终态压实度．采用一维两相流模型模拟了活塞驱动颗粒床形成致密波这一动态过程．用线方法（MOL）对该方程组求数值解．计算表明,经过短暂的非稳态过程,颗粒床中形成一稳态致密波．分析了活塞速度与初始孔隙率对致密波结构的影响,并对简化两相流模型与两相流模型的计算结果进行了对比．     

Shock-induced two-phase flows in an aligned baffle system filled with suspended particles

This paper describes the shock propagation through a dilute gas-particle suspension in an aligned baffle system. Numerical solution to two-phase flows induced by a planar shock wave is given based on the two-continuum model with interphase coupling. The governing equations are numerically solved by using high-resolution schemes. The computational results show the shock reflection and diffraction patterns, and the shock-induced flow fields in the 4-baffle system filled with the dusty gas.     

弯管内爆轰波传播的流场显示和数值模拟

采用激光纹影系统拍摄了爆轰波在不同位置的流场照片. 用二阶附加半隐的龙格- 库塔法和五阶WENO格式 分别离散欧拉方程时间和空间导数项，用基元反应来描述爆轰化学反应过程，获得了压力、 温度、典型组元质量分数分布及数值胞格结构和爆轰波平均速度. 结果表明：受壁面稀疏波 和压缩波影响，爆轰波阵面发生畸变. 但由于弯管曲率半径较大，未出现爆轰波熄灭. 靠近 凹壁面的激波强度大于凸壁面侧，且凹壁面侧的反应区宽度较凸壁面侧要窄. 弯管出口处的 三波点数目较入口处减少，爆轰波衰减. 在出口直段，受扰动的爆轰波可恢复为自持爆轰波. 爆轰波流场、胞格结构、平均爆轰波速度的计算和实验结果定性一致.     

FLAT-PLATE BOUNDARY-LAYER FLOWS INDUCED BY DUSTY SHOCK WAVE

ＦＬＡＴ－ＰＬＡＴＥＢＯＵＮＤＡＲＹ－ＬＡＹＥＲＦＬＯＷＳＩＮＤＵＣＥＤＢＹＤＵＳＴＹＳＨＯＣＫＷＡＶＥ（王柏懿）（陶锋）ＦＬＡＴ－ＰＬＡＴＥＢＯＵＮＤＡＲＹ－ＬＡＹＥＲＦＬＯＷＳＩＮＤＵＣＥＤＢＹＤＵＳＴＹＳＨＯＣＫＷＡＶＥ￥ＷａｎｇＢｏｙｉ；Ｔａ...