平面激波加载下砂墙结构的冲击响应特性

砂墙结构在爆炸安全防护领域具有广泛应用，为了研究激波加载下砂墙结构的冲击响应特性，基于水平激波管实验装置，开展平面激波冲击砂墙结构系列实验，采用高速纹影摄像系统捕捉流场中激波波系的演化过程和砂墙结构的运动过程。入射激波马赫数为1.827～2.413，相应入射激波载荷强度为0.378～0.724 MPa。砂墙结构利用铁砂、矾土、石英砂3种实验用砂制备，所制备砂墙结构孔隙度分别为56.6%、69.3%、56.6%。高速纹影照片显示:平面激波冲击砂墙结构发生反射和透射，伴随入射激波和透射激波的传播，在百微秒内，砂墙未产生显著运动，表现出显著的类固体动力学响应特性。基于冲击理论，确定了铁砂墙、矾土砂墙、石英砂墙的线性冲击关系，冲击关系中线性常数λ值量级为100，根据凝聚介质实用状态方程推断:较低强度载荷冲击作用下，砂墙主要产生体积变形，而由冲击引起的热能效应则可以忽略。     

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

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

楔形炸药窗口材料中的斜激波特性

采用二维不定常流理论中的极曲线方法,对冲击波进入楔形受试炸药/LiF窗口界面的冲击波波后流场状态进行理论分析,同时采用流体动力学软件LS-DYNA对模型的波后流场状态进行数值模拟。结果表明,对于所研究的模型,一维平面正冲击波进入楔形受试炸药/LiF窗口界面时,由于斜激波的存在,波后粒子运动方向发生偏转,偏转角为3.3°。数值模拟结果表明,在相同模型条件下,冲击波波后粒子偏转角在2.77°~3.03°之间。二者的差异源于极曲线理论中未考虑稀疏波的影响。     

惰性介质对钝感炸药爆轰的约束效应

首先用改进冲击波极曲线理论分析惰性介质对钝感炸药爆轰的约束作用类型。改进冲击波极曲线基于爆轰ZND模型建立在前导冲击波上,并且未反应炸药采用JWL状态方程,惰性介质采用p（ρ,T）形式状态方程。理论考察声速小于炸药CJ爆速且压缩性不同的6种典型惰性介质约束情况。然后用带三项式Lee-Tarver化学反应率的二维Lagrange流体力学方法数值模拟考察约束相互作用。数值考察约束介质的影响因素有:压缩性、厚度、典型双层介质组合约束。从数值结果看出,由介质压缩性的不同给出的约束作用方式共7种:其中6种出现在介质声速小于炸药CJ爆速条件下,可运用冲击波极曲线理论;另外一种出现在介质声速大于炸药CJ爆速条件下,不能使用冲击波极曲线理论。同时,介质厚度、双层介质组合方式也能够影响爆轰前导冲击波阵面形状以及爆轰化学反应流动状态。     

柱形汇聚激波冲击球形重气体界面的实验研究

采用高速纹影法实验研究了柱形汇聚激波与球形重气体界面相互作用的 Richtmyer-Meshkov不稳定性问题. 激波管实验段基于激波动力学理论设计, 将马赫数为1.2 的平面激波转化为柱形汇聚激波, 气体界面由肥皂膜分隔六氟化硫(内)和空气(外)得到. 采用高速摄影机在单次实验中拍摄激波运动的全过程, 对柱形激波的形成进行了实验验证, 并进一步观测了汇聚激波与球形气体界面相互作用过程中的波系发展和气体界面变形以及反射激波同已变形界面二次作用的流场演化. 结果表明: 当柱形汇聚激波穿过气泡界面以后, 气泡左侧界面极点沿激波传播方向保持匀速运动, 气泡右侧界面发展成为射流结构, 气泡主体发展成为涡环结构; 在反射激波的二次作用下, 流场中无序运动显著增强并很快进入湍流混合阶段.     

反射激波作用重气柱的Richtmyer-Meshkov不稳定性的实验研究

采用高速摄影结合激光片光源技术,研究了反射激波冲击空气环境中重气体(SF₆)气柱的Richtmyer-Meshkov不稳定性。通过在横式激波管试验段采用可移动反射端壁获得不同反射距离,实现了反射激波在不同时刻二次冲击处于演化中后期的气柱界面,得到了不同的界面演化规律。反射距离较小时,斜压机制对气柱界面形态演化的影响显著,界面衍生出二次涡对结构;反射距离较大时,压力扰动机制的影响显著,界面在流向上被明显地压缩,没有形成明显的涡结构。由气柱界面形态的时间演化图像得到了界面位置和整体尺度随时间的变化,对反射激波作用后气柱界面的演化进行了量化分析。     

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

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

反射激波作用下三维凹气柱界面演化的数值研究

激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例. 单次激波与二维气柱相互作用已得到广泛关注, 但是反射激波再次冲击气柱 (尤其是三维气柱) 的研究较少, 相关演化规律和机理尚不清楚. 反射激波再次冲击演化中的气柱界面会产生新的斜压涡量, 影响涡量的输运和分布, 从而影响界面的演化. 本文采用自主开发的HOWD (high-order WENO and double-flux methods) 程序, 研究了马赫数为1.29的平面激波冲击N激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例.单次激波与二维气柱相互作用已得到广泛关注,但是反射激波再次冲击气柱(尤其是三维气柱)的研究较少,相关演化规律和机理尚不清楚.反射激波再次冲击演化中的气柱界面会产生新的斜压涡量,影响涡量的输运和分布,从而影响界面的演化.本文采用自主开发的HOWD (high-order WENO and double-flux methods)程序,研究了马赫数为1.29的平面激波冲击N_2气柱(气柱外为SF_6)的演化过程,并考察了反射激波对二维和三维凹气柱界面演化的影响规律.在数值模拟中,选取了不同的反射距离(定义为气柱和反射边界的距离),得到了二维和三维凹气柱在反射激波冲击前后的完整演化图像,提取了气柱上特征点位置随时间变化的定量数据,重点分析了不同演化阶段气柱几何特征及斜压涡量分布的变化趋势.研究表明,反射距离决定着反射激波作用气柱时的激波形状和气柱形态,从而影响斜压涡量的生成和分布,进而改变气柱的不稳定性演化过程.对于三维气柱,不同高度截面上的斜压涡量分布不同,从而诱导出复杂的三维演化结构.     

基于立方型气体状态方程的激波风洞准一维流动数值研究

采用基于立方型气体状态方程的准一维流动数值模拟方法研究了反射式高焓激波风洞的真实气体流动,重点关注了高压真实气体效应对风洞全场流动时空结构和驻室区气流参数的影响,并以理论分析揭示了高压真实气体效应对激波管内流动的作用机理。研究表明：对于以冷高压气体驱动的激波风洞,使用考虑分子体积和分子间作用力的真实气体状态方程能够更准确地描述气体的状态和风洞内的流动状况。高压真实气体效应主要在冷驱动气体中发生作用,其作用效果主要是使当地声速增大,从而使得入射稀疏波和反射稀疏波的传播速度加快;另一方面,高压气体效应在高温气体效应较显著的被驱动气体中作用微弱,且对激波管产生激波的强度和激波后的流动状态影响甚微。稀疏波的加快传播改变了激波管波系的相干时空关系。提前抵达的稀疏波可在一定情况下侵蚀激波风洞的有效试验时间。对于所测试的激波风洞构型,在150 MPa氢气驱动110 k Pa氮气的工况下,高压效应导致的有效试验时间缩短约38%。适当加长驱动段长度和采用高温气体驱动均可有效减弱高压真实气体效应的影响。     

柱形装药条件下锥形水中爆炸激波管内的冲击波特性

锥形水中爆炸激波管是进行水中爆炸实验的一种装置,该装置能够通过较小装药量在相同距离处实现自由场水中较大装药量爆炸的冲击波峰值。为了获得柱形装药条件下锥形水中爆炸激波管内的冲击波特性,本文通过数值计算的方式,对不同圆锥角和不同柱形装药质量下锥形激波管内的冲击波传播过程进行了模拟,通过对不同工况下激波管内冲击波特性进行分析,发现其初始冲击波的衰减规律符合自由场水中的指数衰减形式,并拟合得到了与自由场水中爆炸相容的冲击波峰值、比冲量和能流密度经验公式;发现其二次脉动压力周期与炸药质量呈反常规的变化规律,并引入等效静水压深度解释了这一现象;发现其二次脉动压力幅值与初始冲击波幅值之比比自由场水中更大,而二次脉动压力的比冲量与初始冲击波冲量之比与自由场水中相当。     

Generalized and exact solutions for oblique shock waves of real gases with application to real air

Calculation of the oblique shock wave of real gases is a difficult and time consuming problem because it involves numerical solution of a set of 10 equations, two of which (i.e., the equation of state and enthalpy function)—if available—are of a very complicated algebraic form. The present work presents a generalized method for calculating oblique shock waves of real gases, based on the Redlich-Kwong equation of state. Also described is an exact method applicable when the exact equation of state and enthalpy function of a real gas are available. Application of the generalized and the exact methods in the case of real air showed that the former is very accurate and at least twenty times faster than the latter. An additional contribution of the study is the derivation of real gas oblique shock wave equations, which are of the same algebraic form as the well known ideal gas normal shock wave relations.     

非线性弹性介质中冲击波斜反射的研究

本文讨论了各向同性非线性超弹性介质在平面小应变下的冲击波斜反射问题。给出了本构关系、简单波解和冲击波解,并作为例子求解了入射冲击波在自由面的斜反射问题。     

Approximate equation of state of condensed substances

The functions entering into the equation of state of Mie-Gruneisen solids are constructed approximately. In [1] an approximate equation of state was proposed for solids whose adiabatic curves obey a linear relationship between the velocity of the shock wave D and the mass velocity U; the slope of the shock adiabatic curve was equal to 1.5. In the present work, the dimensionless variables proposed in [2] are used to construct an equation of state which is free from the above-mentioned limitation. The equation of state found is used to make calculations of the shock compression of porous metals, in particular, of copper, with values of the porosity not differing significantly from unity.     

Droplet break-up through an oblique shock wave

The interaction of a two-phase flow with a wedge where a stationary shock wave is initially settled is studied in a two-dimensional configuration. Before the introduction of the dispersed phase, the flow around the wedge is a supersonic one phase flow such as an attached stationary shock wave is present. Then, the dispersed phase is introduced upstream the initial position of the stationary shock wave. The purpose of this study is to point out two-phase and droplets break-up effects on the oblique shock wave. The two-dimensional equations are solved by a TVD scheme where fluxes are computed by using Riemann solver for the gas phase equations and also for the dispersed phase equations wich is an original approach due to the authors (Saurel et al. 1994). In addition to drag forces and heat and mass transfers, the process of droplets fragmentation based on the particle oscillation is considered. Accepted April 28, 1995     

On flows behind shock waves reflected from a solid wall

The flow fields associated with the interaction of a normal shock wave with a plane wall kept at a constant temperature were studied based on kinetic theory which can describe appropriately the shock structure and its reflection process. With the use of a difference scheme, the time developments of the distributions of the fluid dynamic quantities (velocity, temperature, pressure and number density of the gas) were obtained numerically from the BGK model of the Boltzmann equation subject to the condition of diffusive-reflection at the wall for several cases of incident Mach number:M ₁=1.2, 1.5, 2.0, 3.0, 4.0, 5.0 and 6.0. The reflection process of the shocks is shown explicitly together with the resulting formation of the flow fields as time goes on. The nonzero uniform velocity toward the wall occurring between the viscous boundary layer and the reflected shock wave is found to be fairly large, the magnitude of which is of the order of several percent of the velocity induced behind the incident shock, decreasing as the incident Mach number increases. It is also seen that a region of positive velocity (away from the wall) within the viscous boundary layer manifests itself in the immediate vicinity of the wall, which is distinct for larger incident Mach numbers. Some of the calculated density profiles are compared with available experimental data and also with numerical results based on the Navier-Stokes equations. The agreement between the three results is fairly good except in the region close to the wall, where the difference in the conditions of these studies and the inappropriateness of the Navier-Stokes equations manifest themselves greatly in the gas behavior.This article was processed using Springer-Verlag TEX Shock Waves macro package 1990.     

Decay of a plane shock wave in a two-parameter medium with arbitrary equation of state

Special curves, called shock polars, are frequently used to determine the state of the gas behind an oblique shock wave from known parameters of the oncoming flow. For a perfect gas, these curves have been constructed and investigated in detail [1]. However, for the solution of problems associated with gas flow at high velocities and high temperatures it is necessary to use models of gases with complicated equations of state. It is therefore of interest to study the properties of oblique shocks in such media. In the present paper, a study is made of the form of the shock polars for two-parameter media with arbitrary equation of state, these satisfying the conditions of Cemplen's theorem. Some properties of oblique shocks in such media that are new compared with a perfect gas are established. On the basis of the obtained results, the existence of triple configurations in steady supersonic flows obtained by the decay of plane shock waves is considered. It is shown that D'yakov-unstable discontinuities decompose into an oblique shock and a centered rarefaction wave, while spontaneously radiating discontinuities decompose into two shocks or into a shock and a rarefaction wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 147–153, November–December, 1982.     

Supersonic low-concentration dusty-gas flow past a plane cylinder in the presence of an oblique shock wave interacting with the bow shock

The motion of an inertial dispersed admixture near a plane cylinder immersed in a steady-state hypersonic dusty flow in the presence of an oblique shock wave interacting with the bow shock is considered. It is assumed that the free-stream particle mass concentration is small and the particles do not affect the carrier flow. The III and IV shock wave interaction regimes are considered. The gas flow parameters in the shock layer are calculated from the numerical solution of the full Navier-Stokes equations for the perfect gas. A TVD second-order finite-difference scheme constructed on the basis of a finite volume method is used. For calculating the dispersed-phase parameters, including the concentration, the full Lagrangian method is used. On a wide range of variation of the particle inertia parameters, the patterns of the particle trajectories, velocity, concentration, and temperature in the shock layer are studied. The possibility of aerodynamic focusing of the particles behind the shock wave intersection point and the formation of narrow beams with a high particle concentration is revealed. These beams impinge on the cylinder surface and result in a sharp increase in the local heat fluxes. The maximal possible increase in the heat fluxes caused by the particles colliding with the cylinder surface is estimated for the flows with and without the incident oblique shock wave.     

On the energy conservation and critical velocities for the propagation of a ＂steady-shock＂ wave in a bar made of cellular material

The propagation of shock waves in a cellular bar is systematically studied in the framework of continuum solids by adopting two idealized material models, viz. the dynamic rigid, perfectly plastic, locking （D-R-PP-L） model and the dynamic rigid, linear hardening plastic, locking （D-R-LHP-L） model, both considering the effects of strain-rate on the material properties. The shock wave speed relevant to these two models is derived. Consider the case of a bar made of one of such material with initial length L 0 and initial velocity v i impinging onto a rigid target. The variations of the stress, strain, particle velocity, specific internal energy across the shock wave and the cease distance of shock wave are all determined analytically. In particular the ＂energy conservation condition＂ and the ＂kinematic existence condition＂ as proposed by Tan et al. （2005） is re-examined, showing that the ＂energy conservation condition＂ and the consequent ＂critical velocity＂, i.e. the shock can only be generated and sustained in R-PP-L bars when the impact velocity is above this critical velocity, is incorrect. Instead, with elastic deformation, strain-hardening and strain-rate sensitivity of the cellular materials being considered, it is appropriate to redefine a first and a second critical impact velocity for the existence and propagation of shock waves in cellular solids. Starting from the basic relations for shock wave propagating in D-R-LHP-L cellular materials, a new method for inversely determining the dynamic stress-strain curve for cellular materials is proposed. By using e.g. a combination of Taylor bar and Hopkinson pressure bar impact experimental technique, the dynamic stress-strain curve of aluminum foam could bedetermined. Finally, it is demonstrated that this new formulation of shock theory in this one-dimensional stress state can be generalized to shocks in a one-dimensional strain state, i.e. for the case of plate impact on cellular materials, by simply making proper replacements of the elastic and plastic constants.