爆轰产物压力势能缓释方法研究

为扩大脉冲爆轰技术的应用领域 ,抑制其振动、噪声等不利因素 ,提高能源转换效率的可能性 ,本文通过实验和数值计算 ,研究了单次爆轰产物通过孔盘将压力势能缓慢释放到相邻腔体的现象。实验在长 2m、内径 40mm的爆轰管中进行 ,数值模拟采用自适应的有限体积方法和有限速率化学反应模型。实验和数值模拟分析的结果均表明 ,在爆轰管与相邻腔体之间放置孔径适当的隔离孔盘 ,可以使爆轰产物的压力势能逐渐地传递到腔体内 ,造成随时间逐渐上升的压力分布。     

三维气相爆轰动态并行计算程序设计与开发

在三维气相爆轰数值研究中,网格精度和计算域的规模导致网格数占有非常庞大的计算资源,进而给数值模拟带来了极大的挑战。本文针对这一难题,采用5阶WENO格式对带化学反应Euler方程组进行空间离散,基于MPI（MessagePassingInterface）并行模式开发了高精度动态并行代码,并对爆轰波在带有障碍物的三维方形管道中的传播过程进行计算。计算结果表明,高精度动态并行计算能够很好的模拟三维气相爆轰波在大尺寸管道中的传播,不仅提高了计算效率,而且提高了爆轰波阵面的分辨率。与高精度静态并行相比,高精度动态并行计算减少了界面数据通信时间,从而进一步提高了计算效率。因此,高精度动态并行程序为探究三维气相爆轰新的物理机制提供有效的手段。     

Numerical simulation of spinning detonation in square tube

A single spinning detonation wave propagating in a square tube is simulated three-dimensionally with the detailed chemical reaction mechanism for hydrogen/air mixture proposed by Petersen and Hanson. The spinning detonation is composed of a transverse detonation rotating around the wall normal to the tube axis, triple lines propagating partially out of phase, and a short pressure trail. The formation of an unburned gas pocket behind the detonation front was not observed in the present simulations because the rotating transverse detonation completely consumed the unburned gas. The calculated profiles of instantaneous OH mass fraction have a keystone shape behind the detonation front. The numerical results for the pitch and track angle on the tube wall agree well with the experimental results. This paper is based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems at Montreal, Canada, from July 31 to August 5, 2005.     

旋转爆轰的三维结构和侧向稀疏波的影响

基于带化学反应的三维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对圆环形燃 烧室内的旋转爆轰进行了数值模拟,讨论了旋转爆轰波的三维结构及侧向稀疏波对旋转爆轰波阵面的影响。 数值结果表明,爆轰波能够以旋转方式沿预混气层稳定传播。在侧向稀疏波作用下,爆轰波阵面发生变形。 与理想的C-J爆轰相比,爆轰波强度和爆轰参数都有所下降。     

有限谱ENO格式在爆轰波数值模拟中的应用

使用有限谱ENO(EssentiallyNon Oscillatory)格式 ,采用V .P .Korobeinikov二步化学反应模型 ,对稀释的化学当量的氢气和氧气混合物中非定常自维持爆轰波进行了二维数值模拟 ,研究了波的产生和演化机理 ,分析了爆轰波的三波结构和传播过程。计算得到的爆轰波参数和结构与以前的计算和实验结果一致。研究表明 :有限谱ENO格式可以成功地模拟非定常自维持爆轰波。     

爆轰波通过扩张喷管的双曝光全息实验和数值研究

结合实验和数值模拟方法,对以脉冲爆轰发动机为背景的爆轰波通过扩张喷管的流动进行了系列研究。实验采用双曝光全息干涉方法对爆轰波绕射流场进行测量,得到了比传统的纹影法更清晰和可定量化的照片。发展了基于非结构四边形网格自适应有限体积程序,结合基元化学反应模型对扩张喷管中爆轰化学反应流场进行了数值模拟,模拟结果与实验照片吻合较好。实验和数值模拟结果表明,爆轰波绕射具有许多和激波绕射不同的流场特征,其中包括二次起爆现象、化学反应面与前导激波相脱离而引起的复杂流场等,同时初始压力和扩张角度变化也对爆轰波绕射过程产生较大影响,初始压力越低,化学反应区和前导激波分离现象越明显,且前导激波的曲率越大。     

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.     

Numerical study on three-dimensional flow field of continuously rotating detonation in a toroidal chamber

Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonationshock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interesting properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.     

一类爆轰合成用乳化炸药的爆轰反应特征

结合热分解特性、爆轰产物状态及数值计算探讨所制备的含Fe、Mn(Zn)元素乳化炸药的爆轰反 应特征。对该系列乳化炸药及其爆轰产物分别进行DSC和DTG实验、XRD检测和TEM 表征,通过比较不 同类型炸药及其爆轰产物的热分解特征、爆轰产物的成分和形貌、晶型畸变度等研究该类炸药爆轰反应特征; 并通过数值计算研究和佐证了不同温度和压强状态下,爆轰固体产物可能存在状态。结果表明:硝酸铵有利 于炸药爆轰时形成相对比较均匀和平稳的爆轰反应结构,而这种结构有利于可重复性合成比较单一和均匀 的爆轰产物。高密度炸药的爆轰产物比低密度炸药的相对纯净。低密度炸药不完全爆轰,导致了爆轰产物杂 质较多;数值计算结果表明,爆轰固体产物分布相图可以辅助分析爆轰固体产物可能存在状态及Mn(Zn)的 爆轰合成机制。     

2021-08期封面

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

聚奥-9C装药的传爆管殉爆

为了研究冲击波作用下引信传爆装置的响应规律,进行了以主发炸药为RDX-8701、被发装置为聚奥-9C（JO-9C）装药的传爆管（含导爆药柱）的殉爆实验。通过观测残留传爆药、壳体和见证块变形,判断传爆管的爆炸程度,分析了殉爆过程中JO-9C爆轰波的成长历程及传播规律。采用AUTODYN软件建立了殉爆实验有限元模型,计算模型中主要考虑了主发炸药产生的爆炸冲击波对传爆管的冲击响应。基于流固耦合方法,通过调整距离模拟计算得到了传爆管的临界殉爆距离和殉爆安全距离。结果表明,传爆管上端的侧角受到爆炸冲击后产生的爆轰波沿斜下方传播,使传爆药柱完全爆轰,并引起导爆药柱发生殉爆;数值模拟结果显示,JO-9C装药的传爆管临界殉爆距离为5.7 mm,殉爆安全距离为8.8 mm。     

爆轰波在楔面上反射数值分析

应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波在楔面反射进行了数值模拟,计算中氢氧混合物的化学反应采用了8种组分20个反应方程式,在处理化学反应引起的刚性问题时采用了时间算子分裂的方法,模拟了爆轰波在楔面反射由马赫反射向规则反射转变的过程,得到了反射转变临界角,同时考虑了初始压力和组分的影响,并和实验及理论分析结果进行了比较,结果是令人满意的。     

气相爆轰胞格结构和马赫反射数值模拟

采用简化二阶段化学反应模型和改进的高精度时空守恒方法 (CE/SE)对二维可燃气体DDT过程和爆轰波在楔面上的马赫反射进行了数值模拟 ,计算和实验结果的比较令人满意。     

Investigation of the detonation regimes in gaseous mixtures with suspended starch particles

The existence of a secondary discontinuity at the rear of a detonation front shown in experiments by Peraldi and Veyssiere (1986) in stoichiometric hydrogen-oxygen mixtures with suspended 20-m starch particles has not been explained satisfactorily. Recently Veyssiere et al. (1997) analyzed these results using a one-dimensional (1-D) numerical model, and concluded that the heat release rate provided by the burning of starch particles in gaseous detonation products is too weak to support a double-front detonation (DFD), in contrast to the case of hybrid mixtures of hydrogen-air with suspended aluminium particles in which a double-front detonation structure was observed by Veyssiere (1986). A two-dimensional (2-D) numerical model was used in the present work to investigate abovementioned experimental results for hybrid mixtures with starch particles. The formation and propagation of the detonation has been examined in the geometry similar to the experimental tube of Peraldi and Veyssiere (1986), which has an area change after 2 m of propagation from the ignition point from a 69 mm dia. section to a 53 mm 53 mm square cross section corresponding to a 33% area contraction. It is shown that the detonation propagation regime in these experiments has a different nature from the double-front detonation observed in hybrid mixtures with aluminium particles. The detonation propagates as a pseudo-gas detonation (PGD) because starch particles release their heat downstream of the CJ plane giving rise to a non-stationary compression wave. The discontinuity wave at the rear of the detonation front is due to the interaction of the leading detonation front with the tube contraction, and is detected at the farthest pressure gauge location because the tube length is insufficient for the perturbation generated by the tube contraction to decay. Thus, numerical simulations explain experimental observations made by Peraldi and Veyssiere (1986). Received 5 July 1997 / Accepted 13 July 1998     

Set-valued solutions for non-ideal detonation

The existence and structure of a steady-state gaseous detonation propagating in a packed bed of solid inert particles are analyzed in the one-dimensional approximation by taking into consideration frictional and heat losses between the gas and the particles. A new formulation of the governing equations is introduced that eliminates the difficulties with numerical integration across the sonic singularity in the reactive Euler equations. With the new algorithm, we find that when the sonic point disappears from the flow, there exists a one-parameter family of solutions parameterized by either pressure or temperature at the end of the reaction zone. These solutions (termed “set-valued” here) correspond to a continuous spectrum of the eigenvalue problem that determines the detonation velocity as a function of a loss factor.     

氢氧混合气体爆轰波的真实化学反应模型数值模拟

采用高精度的ENO格式和基于基元化学反应的真实化学反应模型求解氢氧混合气体一维爆轰波的精细结构。采用直接起爆方法得到稳定传播的爆轰波 ,计算的爆轰波阵面参数和实验相当符合。对爆轰波反应区化学反应的研究表明 ,参与反应的不同组分具有不同类型的变化特征。网格尺寸影响的研究表明 ,计算结果的精度随着网格尺寸的增加而增加 ,并能保持较好的收敛性。移动网格研究结果表明 ,网格运动速度和爆轰速度接近时 ,两者的相互作用对计算结果产生一定影响。     

驻定斜爆轰波并行数值模拟

采用多组分化学反应Euler方程组对驻定在高速飞行弹丸上的斜爆轰波流场进行了数值模拟。计算中分别采用TVD格式和基元反应模型,并基于并行编程模型MPI（message passing interface）实现了非结构网格上的并行计算,对流项和化学反应项用时间分裂法进行处理。计算结果表明并行计算能有效地提高计算速度,扩展计算规模,为进一步研究超驱爆轰推进技术奠定基础。     

Numerical study of three-dimensional detonation structure transformations in a narrow square tube: from rectangular and diagonal modes into spinning modes

Three-dimensional (3-D) detonation structure transformations from rectangular and diagonal modes into spinning modes in a narrow square tube are investigated by high-resolution simulation. Numerical simulations are performed with a Riemann solver of the HLLC-type, new cell-based structured adaptive mesh refinement data structure, high-order, parallel adaptive mesh refinement reactive flow code. A simplified one-step kinetic reaction model is used to reveal the 3-D detonation structure. The four different types of initial disturbances applied in the ZND profiles lead to the structures of rectangular in phase, rectangular out of phase, rectangular partial out of phase and diagonal, respectively, during the initial stages of detonation propagation. Eventually, all these detonation structures evolve into the self-sustained spinning detonations. The asymmetric disturbance leads to a stable spinning detonation much faster than the rest. The important features in the formation of spinning detonation are revealed using a 3-D visualization, and a remarkable qualitative agreement with experimental and numerical results is obtained with respect to the transverse wave dynamics and detonation front structures. The transverse wave collisions produce the unburnt gas pockets and the energy to sustain the detonation front propagation and distortion. The periodic pressure oscillation of front plays a complex role as it shifts the reaction zone structure with an accompanying change in the driving energy of transition and the detonation parameters which result in the more distorted front and the unstable detonation. Eventually, the unstable distorted detonation evolves into a spinning detonation.     

组合状态方程在数值计算上的应用及分析

利用建立在爆轰产物实验等熵线数据上的组合状态方程以及常用的Jones-Wilkins-Lee (JWL)和Becker-Kistiakowsky-Wilson (BKW)爆轰产物状态方程进行了爆炸数值分析。状态方程数值计算的压力波形与实验波形的对比表明:组合状态方程和JWL状态方程计算的波形的到时、峰值、形状能够与实验吻合;BKW状态方程等熵式计算的数值波形峰值基本能够与实验吻合,但形状吻合不佳。     

正向爆轰驱动高焓激波风洞的数值模拟

对充满氢氧可燃气体、带扩容腔的正向爆轰驱动的激波风洞进行了数值模拟。计算采用了欧拉方程,频散可控耗散差分格式（DCD）和改进的二阶段化学反应模型。在扩容腔附近采用二维轴对称计算模型,而在驱动段和被驱动段的直管道部分则采用一维计算模型。本文分析了爆轰波在管道中的传播、反射和绕射过程。计算结果表明扩容腔的尺寸对爆轰波的传播、反射、汇聚等起着决定性的作用;带扩容腔的正向爆轰驱动的激波风洞能够得到平稳的持续时间较长的气流,提高了实验的精确度和可重复性。