一种弹体入水冲击波阴影成像可视化系统及其试验研究

扩大成像视野对于开展充水容器中弹体入水冲击波传播及弥散方面的可视化研究具有重要的实际意义。阴影成像技术适用于大视野实验,且对流场冲击波和扰动的可视化研究具有简单性和通用性,其中直接阴影成像最为简单,但可靠点光源的缺乏是阻碍其发展应用的瓶颈。因此基于国产短弧氙灯管,自制了短弧氙灯点光源,根据阴影成像原理,设计出一种弹体入水冲击波阴影成像可视化系统,详细介绍了其组成和运行原理。利用该系统对高速弹体入水进行了试验研究,获得了弹体入水冲击波的阴影成像和冲击波信号的压力时程曲线,通过阴影成像和冲击波信号相结合分析了弹体入水冲击波的传播特性,并进行了理论验证。结果表明:该弹体入水冲击波阴影成像可视化系统具有可靠性和设计的合理性。弹体高速入水后,初始冲击波的强度最大,随着冲击波的传播,冲击波强度逐渐降低,水中冲击波的传播速度不断降低,球形冲击波的半径逐渐增大。     

高速弹丸诱导斜爆轰激波结构实验研究

斜爆轰推进系统在高超声速推进领域具有广阔的应用前景,其释热迅速、比冲高、燃烧室结构简单的优点吸引研究人员的持续关注.然而,斜爆轰的地面试验同时涉及到高速试验环境模拟、燃料与氧化剂混合、高温燃烧流场结构测量等技术难点,当前国内外系统的试验研究仍然十分有限,难以支撑斜爆轰发动机的研制.为了研究自持传播的斜爆轰激波结构与波面流动特性,基于爆轰驱动二级轻气炮开展了高速弹丸诱导斜爆轰实验研究,使用直径30 mm球头圆柱形弹丸发射进入充满氢/氧可燃混合气体的实验舱中以起爆斜爆轰波,并采用两种阴影技术对实验流动结构进行测量.实验中在不同速度、不同充气压力下观察到三种弹丸诱导激波结构,即激波诱导燃烧、弹丸起爆爆轰波和相对弹丸驻定的斜爆轰波,实验舱充气压力下降则会造成爆轰横波尺度增加与波面流动失稳.实验中,斜爆轰激波角与理论分析结果吻合较好,弹丸气动不稳定带来较大的弹丸攻角会对激波角测量带来一定偏差.通过对斜爆轰波波面法向传播速度的测量发现,随着远离弹丸,斜爆轰传播速度由弹丸飞行速度衰减至接近实验气体CJ速度,弹丸速度的降低会加速斜爆轰波传播速度的衰减.     

环形激波绕射, 反射和聚焦的数值模拟研究

应用频散可控耗散格式对环形激波在圆柱形激波管内绕射、反射和聚焦的问题进行了数值模拟研究，研究结果表明环形激波形成强烈聚焦的关键因素是环形激波在圆柱形管道中向对称轴运动时，绕射激波就不断加速而不作通常情况下的衰减；不同马赫数的环形激波绕射也产生不同马赫数及形状的准柱形激波，导致聚焦效果和位置的差异；另外，环形激波聚焦于一个点而圆柱形激波聚焦于一条线，两者有本质不同。     

弹体超高速侵彻石灰岩靶体地冲击的数值模拟研究

为探究超高速动能武器的对地破坏效应及其影响因素,采用数值模拟方法对弹体超高速侵彻的地冲击规律进行了研究。首先,基于石灰岩静动态力学性能实验数据对材料模型参数进行了标定,并对已有弹体大范围着速侵彻石灰岩靶体进行了模拟,验证了所采用材料模型和数值模拟方法的合理性。随后,开展了钨合金长杆弹超高速侵彻石灰岩靶体的数值模拟,细致分析了地冲击传播的现象和机理:弹体超高速侵彻靶体时,弹靶交界面处会产生瞬时高压,并以应力波的形式在靶体中传播,对靶体内部造成破坏,且当弹体初速度高于3.0 km/s时,地冲击显著增强。最后,进一步研究了不同弹靶参数对地冲击的影响,发现从相对深度来看,弹体参数（弹体长径比、密度）对地冲击规律影响不大;而靶体特征特别是孔隙率对地冲击传播具有较大影响。     

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

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

用环形激波聚焦实现爆轰波直接起爆的数值模拟

利用基元反应模型和有限体积法对环形激波在可燃气体中聚焦实现爆轰波直接起爆进行了数值模拟。研究结果表明，标准状态下的氢气-空气混合气体在马赫数为3.1以上的环形激波聚焦产生的高温高压区作用下会诱发可燃气体的直接起爆形成爆轰波，爆轰波与激波和接触间断相互作用产生了复杂的波系结构；爆轰波爆点位置在对称轴上并不是固定的点，而是随着初始激波马赫数的变化而发生移动；可燃气体初始温度和压力对起爆临界马赫数都有影响，但是初始温度的影响大得多。     

Experimental investigation of converging shocks in water with various confinement materials

Fluid-solid coupling typically plays a negligible role in confined converging shocks in gases because of the rigidity of the surrounding material and large acoustic impedance mismatch of wave propagation between it and the gas. However, this is not true for converging shocks in a liquid. In the latter case, the coupling can not be ignored and properties of the surrounding material have a direct influence on wave propagation. In shock focusing in water confined in a solid convergent geometry, the shock in the liquid transmits to the solid and both transverse and longitudinal waves propagate in the solid. Shock focusing in water for three types of confinement materials has been studied experimentally with schlieren and photoelasticity optical techniques. A projectile from a gas gun impacts a liquid contained in a solid convergent geometry. The impact produces a shock wave in water that develops even higher pressure when focused in the vicinity of the apex. Depending on the confining material, the shock speed in the water can be slower, faster, or in between wave speeds in the solid. For solid materials with higher wave speeds than the shock in water, regions in the water is put in tension and cavitation occurs. Materials with slower wave speeds will deform easily.     

Numerical simulation of flow fields induced by a supersonic projectile moving in tubes

Numerical study of the shocked flows generated by a supersonic projectile released from a launch tube into a big chamber has been performed in this paper. Based on fixed Cartesian grids, the two-dimensional axisymmetric Euler equations are solved by the fifth-order WENO scheme implemented with moving boundary conditions. Both the level set technique and ghost fluid method are used for capturing the moving interface of the projectile implicitly. The numerical results show that complex shock phenomena exist in the transient shock flow, resulting from shock-wave reflection, shock-wave focusing, shock-wave/projectile interaction and shock-wave/contact surface interactions. The relationships between the acceleration of the projectile and the transient shock flow are also discussed in detail.     

Numerical investigation of axisymmetric shock wave focusing over paraboloidal reflectors

The problem of a plane shock wave incident to a paraboloidal reflector is numerically investigated. The numerical solver used is developed by an improved, implicit, upwind total variation diminishing scheme in a finite-volume approach. The real-gas effect is taken into account if high temperature occurs. The solver is validated on four test problems. The complicated flow fields of axisymmetric shock wave focusing for different-depth reflectors at various incident shock Mach numbers are studied. An interesting result of a maximum pressure happening at the reflector center is found. This is due to the occurrence of an implosion phenomenon. A maximum temperature might occur at the reflector center or at other locations, depending on the incident shock Mach number and the reflector depth. Moreover, vortical flows induced by shock wave focusing and their formation mechanism are explored. It was found that the vortices near the reflector are caused by a ring-shaped shock/slipline interaction. Owing to the slipline on the symmetry axis, a jet flow is induced, resulting in the formation of vortices near the symmetry axis. Received 13 January 1998 / Accepted 10 November 1998     

Cavitation phenomena in extracorporeal microexplosion lithotripsy

An experimental investigation was made of cavitation phenomena induced by underwater shock wave focusing applied to the extracorporeal microexplosion lithotripsy (microexplosion ESWL). Firstly an underwater microexplosion generated by detonation of a 10 mg silver azide pellet was studied and secondly underwater shock focusing and its induced cavitation phenomena were investgated. Underwater shock wave was focused by using a semi-ellipsoidal reflector in which a shock wave generated at the first focal point of the reflector was reflected and focused at the second focal point. It is found that an explosion product gas bubble did not produce any distinct rebound shocks. Meantime cavitation appeared after shock focusing at the second focal point where expansion waves originated at the exit of the reflector were simultaneously collected. A shock/bubble interaction is found to contribute not only to urinary tract stone disintegration but also tissue damage. The cavitation effect associated with the microexplosion ESWL was weaker in comparison with a spark discharge ESWL. The microexplosion ESWL is an effective method which can minimize the number of shock exposures hence decreasing tissue damage by conducting precise positioning of urinary tract stones.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

自然超空泡与尾喷流相互作用的数值模拟

认识带尾喷流和自然超空泡的水下高速航行体流体动力特性并发展其预报与控制方法仍是水动力学领域极具挑战性的课题．本研究采用CFD方法对尾喷流和自然超空泡之间的相互作用进行了数值研究．针对发动机欠膨胀超音速喷流,采用现有实验结果验证了基于两方程湍流模型的二维轴对称流动数值模型的可靠性．尾喷流在空气和蒸汽环境中流动的数值计算结果表明,由于蒸汽环境中背压较低,欠膨胀尾喷流无法及时形成压缩波,使得蒸汽环境中尾喷流的过膨胀区和气相扩散区的体积比空气中大;尾喷流很难形成规则的激波格栅,波系结构相对简单．针对携尾喷流的平头航行体超空泡流状态的数值模拟结果表明,尾喷流注入超空泡后迅速充满航行体周围的空腔区域;尾喷流与超空泡尾迹区域形成的回射流相互作用最终导致超空泡断裂,断裂过程中伴随着燃气泡的下泄现象;受空泡壁面约束,尾喷流难以在狭窄的超空泡空腔内完全膨胀,尾喷流的激波波系结构有显著的变化：在喷嘴附近受到压缩,在远离喷嘴区域受到超空泡水汽掺混的破坏;空泡内压强基本维持在饱和蒸汽压附近,没有显著增加．     

环形激波绕射、反射和聚焦流场的间断有限元模拟

采用间断有限元方法对环形激波在圆柱形激波管内绕射、反射和聚焦流场进行了数值模拟。将二维守恒方程的间断有限元方法发展到轴对称Euler方程,并对环形激波绕后台阶流动进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉运动激波在圆柱形激波管内传播的复杂流场结构;在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力,在聚焦点附近不会产生振荡或抹平间断现象。     

球形弹丸正撞击薄板防护屏碎片云特性研究

以质量、动量和能量守恒为基础,结合平面激波理论和热力学理论,对球形铝弹丸正撞击薄板铝防护屏形成的碎片云特性进行了建模,模型计算结果与实验结果吻合较好。利用该模型对多种工况下碎片云特性进行了计算,结果表明:(1)碎片云质心速度和膨胀速度随撞击速度和弹丸直径的增加而增大,随防护屏厚度增大而减小;膨胀半角随撞击速度和防护屏厚度的增大而增大,随弹丸直径的增大而减小;(2)速度和膨胀半角变化曲线具有相似性;(3)对于给定的弹丸和防护屏材料,碎片云中受到激波加载部分的材料各相态质量分数只与弹丸撞击速度有关。这些规律与实验规律吻合。     

激波的传播与干扰

激波的传播特性既取决于激波的产生条件,也与所处的传播环境密切相关。驱动条件、几何边界、介质的物理化学属性等发生变化时,都会引起激波传播特性的改变,而激波的变化反过来又会对其波及的流场产生影响。尽管激波传播及其干扰现象广泛存在于自然界和人类科技活动之中,其复杂机理的认识、规律的描述乃至应用潜力的挖掘仍有漫长的路要走。本文根据气体中激波传播和干扰现象以及与之相关的理论描述特征,在对激波传播以及反射、折射等基本现象进行简要阐述的基础上,重点围绕目前的热点问题,包括激波/激波干扰、激波/边界层干扰、激波与湍流作用、激波的聚焦与点火以及激波作用下气体界面不稳定性等研究进行了介绍和讨论,旨在对近年来该领域的进展及获得的成果做一个概述和归纳,期望对将来的深入研究有一个鉴借意义。     

Shock wave propagation and interactions

激波的传播特性既取决于激波的产生条件,也与所处的传播环境密切相关.驱动条件、几何边界、介质的物理化学属性等发生变化时,都会引起激波传播特性的改变,而激波的变化反过来又会对其波及的流场产生影响.尽管激波传播及其干扰现象广泛存在于自然界和人类科技活动之中,其复杂机理的认识、规律的描述乃至应用潜力的挖掘仍有漫长的路要走.本文根据气体中激波传播和干扰现象以及与之相关的理论描述特征,在对激波传播以及反射、折射等基本现象进行简要阐述的基础上,重点围绕目前的热点问题,包括激波/激波干扰、激波/边界层干扰、激波与湍流作用、激波的聚焦与点火以及激波作用下气体界面不稳定性等研究进行了介绍和讨论,旨在对近年来该领域的进展及获得的成果做一个概述和归纳,期望对将来的深入研究有一个鉴借意义.     

激波聚焦诱导点火和爆轰的数值研究

以二维轴对称多组分Euler方程为基础，采用非正交结构化网格和改进的波传播算法，模拟了激波在抛物形反射壁面聚焦反射诱导点火和爆轰的过程，描述了其流场形态。讨论了预混气组成、入射激波强度及反射壁面形状对点火和爆轰的影响。结果表明，激波在抛物形反射壁面顶点处聚焦反射可形成局部高温高压区域，该区域在一定条件下可点燃预混气甚至形成爆轰，其中低稀释剂浓度的预混气、较大的入射激波Mach数和较深的反射壁面有利于可燃预混气形成爆轰。     

椭球罩作用下的水下爆炸冲击波反射聚焦模型

基于冲击波传播、非线性反射和聚焦理论,建立水下爆炸冲击波在椭球罩作用下的反射聚焦模型。讨论自由传播、壁面反射和定向聚焦阶段的冲击波特性和压力计算方法,利用波前和波法线的近似方程构建压力场的数值计算域,进而模拟聚焦过程,并与现有实验结果进行对比,结果表明:所建模型可为正聚焦压力提供满足工程精度的预测,并能描述水下冲击波及产生的拉伸波聚焦过程中的一些细节;椭球罩能有效地聚焦水下冲击波,在动力学焦点附近获得有效增益区,在近轴方向上明显削弱冲击波压力衰减;理想条件下的动力学焦点一般位于几何焦点之前,但实际的反射罩变形和背向位移将使其发生后迁,甚至能越过几何焦点。     

Penetration of Cylindrical Projectiles into Saturated Sandy Media

The time and depth of vertical one-dimensional projectile penetration into sandy media in the near shore region are derived. A precise definition for the physical properties and for the behavior of the sandy medium following the projectile impact are evaluated. Three separate time intervals following projectile impact are identified. During the first 3 ms of penetration, the deviatoric friction stress is shown to be negligible and the integrated Mie–Grüneisen equation of state (or, equivalently, the Hugoniot-adiabat) may be applied to compute the normal penetration resistance force from the sand pressure. In order to compute sand pressure as a function of the sand density D by the integrated Mie–Grüneisen equation of state, the Mie–Grüneisen dimensionless constants γ₀ and s and the dimensional speed of sound C ₀ in the sandy medium are required. In order to illustrate the one-dimensional shock wave propagation in both wet and dry sands, Hugoniot data for wet and dry silica sands are evaluated by a three degrees of freedom algorithm to compute these required constants. The numerical results demonstrate that the amplitude of the shock wave pressure in the wet silica sand (41% porosity) is approximately one-third of the shock wave pressure amplitudes in the dry silica sands (22% and 41% porosity). In addition, the shock wave pressure dampens quicker in the wet sand than in the dry sands.     

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.     

Impact on the surface of a compressible liquid

The evolution of a weak, nearly plane shock wave produced by the impact on the plane boundary of a compressible liquid is considered. At the initial moment the liquid is at rest and occupies the lower half-plane. Then the points of its boundary get instantly velocities directed into the liquid domain. This leads to the formation of a shock wave the intensity of which is non-uniform due to a non-uniform distribution of the impact velocities. Initially the shock wave is plane but then it bends due to the non-linear effects and can later be focused. To analyze the liquid flow, the method of matched asymptotic expansions is used. For finite times the flow and the evolution of the shock wave are described within the framework of the acoustic approximation. For large times the flow becomes non-linear, and the form of the shock front depends essentially on the characteristics of the liquid flow behind it. If the non-uniformity of the impact velocity distribution is slight then the focusing of the shock wave is shown not to occur. The influence of viscosity of the liquid on the structure of its motion is discussed.