点火能对丙烷-空气预混气体爆炸过程及管壁动态响应的影响

在长12 m的无缝不锈钢直管中,通过改变初始点火能量,探究了点火能对封闭管道内丙烷-空气混合气体爆炸传播特性和激波对管壁动态加载的影响。结果表明,初始点火能对预混气体爆炸火焰传播规律以及管壁的动态响应有显著影响:点火能越大,爆炸越剧烈,爆炸压力峰值压力和管壁最大应变就越大,且压力波和管壁应变的发展一致。火焰在传播过程中受到管道末端反射波的作用会发生短暂熄灭和复燃;管壁承受冲击波加载,应变信号主要分布在0~781.25 Hz,管壁最大应变率大于10-3 s-1,实验工况下管壁应变属动态响应。     

矩形凹槽管道中激波传播的数值研究

基于二维Euler方程,结合五阶加权基本无振荡(weighted essentially nonoscillatory,WENO)格式以及自适应网格加密(adaptive mesh refinement,AMR)技术对入射激波在矩形凹槽管道内传播过程进行了数值模拟。数值结果清晰地显示了入射激波传播过程中与多个矩形凹槽作用以及在凹槽内变化的整个过程,且与已有的实验结果吻合较好。另外,结果还揭示了入射激波与单个凹槽作用时,会发生绕射产生膨胀波,还会发生碰撞从而诱导反射激波。膨胀波会导致入射激波压力降低,而反射激波则导致其升高,但膨胀波的影响占主导作用,因而入射激波波阵面强度出现振荡下降。     

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

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

可燃气体中激波与障碍物作用在下游形成爆轰波的数值研究

对平面激波和单个矩形障碍物作用的过程进行了数值模拟,研究了反射产生的上行爆轰波在下游可燃气体中形成爆轰波的过程。数值结果表明,下游爆轰波形成过程主要有2种模式:爆轰波直接绕射和绕射波在上壁面反射,这和已有的实验结果是一致的。通过研究下游爆轰波的形成过程受入射激波马赫数、混合气体的压力和管道尺度的影响,分析了上游爆轰波向下游传播的波动力学过程,讨论了2种形成过程的作用规律和控制因素,阐明了下游爆轰波的形成规律。     

柔性障碍物对甲烷空气爆炸波激励作用的实验研究

为研究柔性障碍物对甲烷空气爆炸波的激励效应,采用双向拉伸聚丙烯（biaxially oriented polypropylene, BOPP）薄膜作为柔性障碍物将管道内甲烷空气预混气体与空气隔开,对比障碍物前后火焰、激波变化,分析膜状柔性障碍物激励效应的机理。实验结果表明:这种具有一定承压能力的柔性障碍物对甲烷爆炸波产生的激励效应不可忽视,在膜片破裂前产生多次激波反射过程,可诱导湍流火焰形成,促使膜前爆炸压力提高,膜片破裂后,火焰在伴流作用下传播速度突增,并加速逐渐逼近前驱冲击波,致使膜后爆炸压力大幅提高;激励效应可使膜片前后最大爆炸压力相差5倍,火焰速度相差7倍;另外在膜片位置2.5 m后增设一道膜片,可增强这种激励效应,而增加膜片的实质是使激波火焰相互作用的次数增加。     

管长对带弯头管道内可燃预混气体爆炸规律及薄壁动态响应的影响研究

基于气云燃爆实验平台,开展了不同长度弯曲管道内丙烷-空气预混气体燃爆实验测试与分析,探索管道长度对弯曲管道内可燃气体燃爆特性及爆炸冲击波对管道薄壁加载效应的影响。实验分别采用光电传感器、压力传感器和应变传感器对管道内爆炸参数和管道薄壁的应力应变规律进行了测试,结果表明,管道越长,弯曲段两侧壁面超压和管壁最大应变越大,压力时程曲线与管道薄壁应变时程曲线具有良好的一致性;弯管一定程度上激励了火焰加速传播,火焰传播速度在90°拐点处降至最低;管道长度对管道内火焰燃爆规律的影响是整体的,管道长度越长,管内压力和火焰速度越大,最大爆炸压力出现在管道末端,冲击波对管壁的加载属于动态加载。     

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程，研究了形 成爆轰波的点火特性. 对理想化学当量比氢气/空气混合气体，在初始压强20kPa的条件下， 马赫数2.6-2.8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的，第 2个点火区是由入射激波在抛物面上发生马赫反射引起的. 这种条件下流场中会出现爆燃转 爆轰，起爆点分别分布在管道壁面、抛物反射面和第2点火区附近. 起爆机理分别为激波管 道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相 互作用. 不同的点火和起爆过程导致了不同的流场波系结构，同时影响了爆轰波传播的波动 力学过程.     

气相爆轰波冲击气固界面的透反射特性

为了研究气相爆轰波冲击气固界面过程中透射波和反射波的相关特性,建立爆轰波冲击气固界面的一维理论模型,对不同初始压力条件下爆轰波到达气固界面后的界面两侧的压力和界面速度变化进行分析。利用时空守恒元求解元方法对气相爆轰波冲击气固界面过程进行数值模拟,分析气体部分反射波的压力分布和速度变化规律及透射入固体中应力波的波形和波速特征,并搭建气相爆轰波冲击活塞实验装置进行进一步验证。结果表明:气体爆轰波到达气固界面后,在固体中透射指数形式的弹性波,并在界面处向气体区反射一道激波。爆轰波后的稀疏波与反射激波相交,削弱反射激波,最终形成稳定激波回传。气固界面在稀疏波和反射稀疏波的作用下,压力和速度逐渐下降,最终也形成稳定状态。在不同混气初始压力情况下,爆轰波冲击过程中产生的最高压力和爆压的比值基本保持不变。理论模型对特征点相关物理量的计算值和实验数据符合的较好。     

激波与可运动颗粒群相互作用反射与透射机理实验研究

对水平圆柱形激波管内可压缩性气体与颗粒群的相互作用进行了实验研究与理论分析。利用由压力传感器、信号放大器、示波器和计算机组成的压力测量系统对激波与颗粒群作用时的动态压力进行了测量。发现激波管内发生的是一个复杂的过程,包括激波与颗粒群作用时伴随了激波和膨胀波的反射与透射现象、激波和膨胀波与接触面的干涉、以及激波从激波管端壁的反射等现象。当颗粒装载比α=1时,透射激波被直径为6mm的颗粒群衍射并且有膨胀波紧随其后,因此由透射激波引起的第一个压力峰急剧下降。颗粒装载比和颗粒直径影响透射激波的衍射和衰减。     

波纹管道阻火器内火焰传播的实验与数值模拟研究

对乙烯-空气预混火焰在波纹管道阻火器中的传播与淬熄过程进行了实验和数值模拟研究,实验结果显示:当乙烯接近当量浓度时,预混气体爆炸压力变化过程可分为4个阶段,等压燃烧阶段、缓慢上升阶段、快速上升阶段和压力振荡阶段;在爆炸过程中,由于反射压力波和火焰相互作用的影响,超压值出现多次振荡,压力振荡阶段一般可以持续数十毫秒;乙烯-空气火焰传播速度随管径增加、阻火单元波纹高度减小呈递增趋势,而且随着阻火单元厚度的增加,阻火器的阻火能力明显提高,可以更有效地使火焰淬熄。数值模拟结果显示:在管道封闭端点火后,火焰面呈半球形并以层流扩散的方式向四周传播;当火焰传播到管道壁面时,在管道壁面的约束作用下,火焰面发生变形,壁面附近的火焰逐渐超过了管道轴线附近的火焰,最后形成了“郁金香”状的火焰结构;当爆燃火焰经过阻火单元时,高温已燃气体被其吸收大量热量,同时在反应区产生的稀疏波作用下,气体温度逐渐降低、化学反应速率迅速减小,最终导致火焰被熄灭。通过模拟计算结果可以看出,在整个爆炸过程中,火焰传播速度与爆炸压力波动均较为明显。并提出了孔隙率和阻火单元厚度对火焰传播的影响机制。基于传热学理论模型,并结合实验数据,得出了爆燃火焰速度与爆炸压力之间的关系,为工业装置阻火器的设计和选型提供更为准确的参考依据。     

Interaction of a reflected shock from a concave wall with a flame distorted by an incident shock

Observations are presented from calculations where a laminar spherical CH₄/air flame was perturbed successively by incident and reflected shock waves reflected from a planar or concave wall. The two-dimensional axi-symmetric Navier–Stokes equations with detailed chemistry were used. The computational results were qualitatively validated with experiments which were performed in a standard shock tube arrangement. Under the influence of the incident shock wave, a Richtmyer–Meshkov instability is induced in the flame, and the distorted flame finally takes the form of two separated elliptical burning bubbles in the symmetric cross plane. Then, under subsequent interactions with the shock wave reflected from the planar or the concave wall, the flame takes a mushroom-like shape. Transverse waves produced by the shock reflection from the concave wall can compress the flame towards the axis, and the focusing shock generated on the concave wall will lead to a larger mushroom-like flame than that induced by the planar reflection.      

Pressure wave propagation in liquid-filled tubes of viscoelastic material

The propagation of small-amplitude waves in a thick-walled long viscoelastic tube of variable cross-section, filled with a viscous incompressible fluid, is considered with account for wave reflection at the tube end in application to arterial pulse wave propagation. A solution is obtained in the form of expansions in a small parameter. The effect of the coefficient of wave reflection at the tube end and the wall material parameters on the fluid volume flow-rate and the tube wall displacement is investigated. It is shown that the volume flow-rate phase spectrum characteristics depend only slightly on the wall properties and can be used in clinical diagnostics for finding the reflection coefficient from pressure and flow-rate records.     

生物激波管及动物耐受性的实验研究

在冲击伤研究中,必须针对两种试验方式分别对致伤机理、耐受性和机体瞬时动力学响应等进行研究。只有这样才能得到近真实的研究结果。文中介绍了能够模拟爆炸波波形因素和满足动物冲击伤实验研究物理参数要求的BST-Ⅰ型生物激波管。它是国内第一台用于生物实验的激波管。文中有应用该管开口和端板的两种实验方式,选用51只成年雄性杂种狗作耐受性实验观察。对实验结果作了分析。     

Experimental and numerical studies on interactions of a spherical flame with incident and reflected shocks

Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock tube arrangement, in which a high-speed shadowgraph imaging system is used to record evolutions of the flame. Numerical simulations are conducted by using second-order wave propagation algorithms, based on two-dimensional axisymmetric Navier-Stokes equations with detailed chemical reactions. Qualitative agreements are obtained between the experimental and numerical results. Under actions of incident shock waves, Richtmyer-Meshkov instability responsible for the flame deformation is induced in the flame, and the distoned flame takes a barrel shape. Then, under subsequent actions of the shock wave reflected from a planar wall, the flame takes an inclined non-symmetrical kidney shape in a symmetric cross section, which means a mushroom-like shape of the flame comes finally into being. The vorticity direction in the ring cap has been altered by the reflected shock＇s action, which makes the head of the mushroom-like flame extend quickly to the side wall.     

Shock wave attenuation by grids and orifice plates

The interaction of weak shock waves with porous barriers of different geometries and porosities is examined. Installing a barrier inside the shock tube test section will cause the development of the following wave pattern upon a head-on collision between the incident shock wave and the barrier: a reflected shock from the barrier and a transmitted shock propagating towards the shock tube end wall. Once the transmitted shock wave reaches the end wall it is reflected back towards the barrier. This is the beginning of multiple reflections between the barrier and the end wall. This full cycle of shock reflections/interactions resulting from the incident shock wave collision with the barrier can be studied in a single shock tube test. A one-dimensional (1D), inviscid flow model was proposed for simulating the flow resulting from the initial collision of the incident shock wave with the barrier. Fairly good agreement is found between experimental findings and simulations based on a 1D flow model. Based on obtained numerical and experimental findings an optimal design procedure for shock wave attenuator is suggested. The suggested attenuator may ensure the safety of the shelter’s ventilation systems.     

可压缩流向涡与反向运动激波相互作用的实验

对可压缩流向涡与反向运动激波相互作用的现象进行了实验研究．实验在９４ｍｍ×９４ｍｍ的方截面激波管中进行．在实验段上游安装了一个有限翼展平直机翼．当入射激波通过机翼后，波后２区气流在模型翼尖诱导出一条流向涡．入射激波在激波管端壁反射后，形成的反射激波在观察窗处和流向涡发生作用．实验中拍摄了激波与流向涡作用全过程的纹影照片，观察到了一些和定常激波与旋涡相互作用不同的现象，并与数值计算结果进行了初步比较     

Attenuation of shock waves propagating in polyurethane foams

Shock wave attenuation in polyurethane foams is investigated experimentally and numerically. This study is a part of research project regarding shock propagation in polyurethane foams with high-porosities = 0.951 ~ 0.977 and low densities of ρc = 27.6 ~55.8 kg/m³. Sixty Millimeter long cylindrical foams with various cell numbers and foam insertion condition were installed in a horizontal shock tube of 50 mm i.d. and 5.4 mm in length. Results of pressure measurements in air/foam combination are compared with CFD simulation solving the one-dimensional Euler equations. In the case of a foam B fixed on shock tube wall, pressures at the shock tube end wall increases relatively slowly comparing to non-fixed foam, free to move and a foam A fixed on shock tube wall. This implies that elastic inertia hardly contributes to pressure build up. Pressures behind a foam C fixed on shock tube wall decrease indicating that shock wave is degenerated into compression wave. Dimensionless impulse and attenuation factor decrease as the initial cell number increases. The momentum loss varies depending on cell structure and cell number.     

方管内汽油-空气混合气体密闭爆炸和泄爆特性研究

为研究汽油-空气混合气体密闭爆炸和泄爆特性，采用可视化方管进行了两种爆炸模式实验研究，并基于壁面自适应局部涡黏(wall-adapting local eddy-viscosity,WALE)模型和Zimont预混火焰模型进行了数值模拟研究。结果表明:（1）泄爆工况超压-时序曲线峰值数量多于密闭爆炸工况，且泄爆工况超压-时序曲线存在剧烈的类似简谐振动的振荡，而密闭爆炸工况的爆炸超压特征参数显著高于泄爆工况；（2）密闭爆炸工况最大火焰传播速度明显小于泄爆工况，但前者在火焰传播初期即达到最大值，而后者在火焰传播末期才达到最大值；（3）密闭爆炸工况出现郁金香形火焰，而泄爆工况出现蘑菇形火焰，郁金香火焰的形成与管道内火焰锋面、流场和流场动压三者之间耦合效应相关，蘑菇形火焰由外部流场湍流和斜压效应的共同作用引起。