隔板深度对激波聚焦起爆影响的数值模拟研究

考虑几何结构参数对激波聚焦触发爆轰波的复杂影响,对H2/Air预混气的环形射流激波聚焦起爆现象开展了数值模拟研究,详细分析了不同隔板深度条件下的激波聚焦过程、流场演化特征以及爆轰波参数变化规律。研究结果表明,凹腔内激波聚焦诱导的局部爆炸以及隔板前缘处射流形成"卷吸涡"是引起爆轰波触发的两个重要机制,而隔板深度是影响环形射流激波聚焦起爆性能的关键因素。随着隔板深度的增加,凹腔内激波聚焦的强度逐步增强,回传的能量损失有所减小,进而导致爆燃转爆轰的距离与时间显著缩短。此外,当隔板深度由1 mm逐渐增加至3 mm时,爆轰波自持传播稳定性呈现出先降低后升高的变化趋势,产生这一现象的主要原因是爆轰波强度与三波点运动的相互作用。     

一种准二维柱面爆炸波加载装置的设计与验证

基于竖直爆轰管和径向Hele-Shaw Cell,设计并搭建了一套准二维柱面爆炸波加载装置,可以实现对Hele-Shaw Cell内部材料界面的径向冲击加载.竖直爆轰管内部的预混气体在底部点燃后,形成向上传播的冲击波,冲击波冲破爆轰管开口与Hele-Shaw Cell底板开孔之间的隔膜后,被Hele-Shaw Cell...     

Focusing of an oscillating shock wave emitted by a toroidal bubble cloud

An analysis of pressure-field dynamics is performed for an axially symmetric problem of interaction between a shock wave and a “free” bubble system (toroidal cluster) giving rise to a steady oscillating shock wave. The results of a numerical study of near-axis wave structure are presented for a focusing shock wave emitted by a bubble cluster. It is shown that the wave reflected from the axis has irregular structure. The Mach disk developing on the axis has a core of finite thickness with a nonuniform radial pressure distribution. The evolution of the Mach-disk core is analyzed, and the maximum pressure in the core is computed as a function of the gas volume fraction in the cluster. The effect of geometric parameters of the toroidal bubble cloud on the cumulative effect is examined.     

内聚锥形激波的非均匀强化与结构演变特征

以高超声速内转式进气道流动中的激波汇聚问题为背景,考虑工程实际中的来流和壁面几何条件这两个关键因素,分别提出了以来流攻角为研究参数的非轴向来流内锥流动模型,和以长/短轴比为研究参数的椭圆入口内锥流动模型.采用激波风洞实验观测和数值模拟相结合的方法,揭示了两类流动中激波的非均匀汇聚特征.结果表明:由来流攻角引起的激波初始沿周向强度分布的不均匀性会在汇聚过程中被放大,迎风面和背风面的激波差异不断加剧;来流攻角越大,初始激波强度不均匀性越强,在汇聚过程中激波面越容易出现不连续的拐折,且出现拐折后激波的汇聚效应会被削弱.由椭圆入口形成的等强度激波在初始时周向的几何不均匀性使激波在汇聚过程中出现沿长/短轴方向的强度差异,激波沿长轴方向上的强度增加更迅速;椭圆长/短轴比越大,激波初始几何不均匀性越强,在汇聚过程中长/短轴两个方向激波强度差异凸显得越快,波面越容易出现不连续的拐折,进而削弱激波的汇聚.在偏离轴对称达到一定程度时,这两种条件下的激波汇聚都会出现汇聚中心处从Mach反射向规则反射的转变.      

激波聚焦反射的实验和数值研究

 采用高速摄影技术和数值模拟方法,对入射激波在两种不同形状的抛物形反射器表面聚焦和反射的过程进行了研究,得到了激波聚焦反射过程的波系结构的实验阴影照片和数值计算结果,两者符合得很好。对激波聚焦形成的气体动力学焦点的特性进行了分析,结果表明,入射激波在两种反射器反射后聚焦所形成的气体动力学焦点均是由三波点在轴心处的会聚所导致的,气体动力学焦点位于相应的反射器壁面的几何焦点附近。不同的反射器中,激波聚焦前后的波系结构也不同,对较浅的反射器,入射激波反射前在反射器壁面形成了弓形激波,反射之后需要相对较长的时间完成聚焦,形成聚焦反射激波之后弓形激波仍未相交;对较深的反射器,入射激波反射后在更短的时间内聚焦,聚焦时弓形激波已经相交,聚焦反射激波之后的流场波系结构更加复杂。     

An experimental study of the diffusion-controlled self-ignition of hydrogen in a channel

The pulsed outflow of hydrogen into channels of circular and rectangular cross sections with a surface area of 20 mm² was experimentally studied. It was revealed that the shock wave formed during the outflow of a pulsed jet is the reason why it ignites at the contact surface. The range of initial pressures of hydrogen at which it ignites was determined and the dependence of the distance from the diaphragm at which a flame arises at the contact surface on the pressure in the shock wave front for circular and rectangular cross section channels was obtained.     

Der Einfluß einer zylindrisch konvergierenden Stoßwelle auf ein Funkenplasma

The influence of a cylindrically converging shock wave on a spark plasma is investigated. The shock wave is produced by a spark discharge in the axis of a cylindrical vessel. After the reflection at the wall the shock wave returns to the spark channel. The results of time resolved spectroscopic studies of the plasma state are compared with calculations on the basis of a gas-dynamic model. The influences of possible deviations from thermodynamic equilibrium are critically discussed.     

Numerical Study of Air Flow Induced by Shock Impact on an Array of Perforated Plates

This study is focused on the propagation behavior and attenuation characteristics of a planar incident shock wave when propagating through an array of perforated plates. Based on a density-based coupled explicit algorithm, combined with a third-order MUSCL scheme and the Roe averaged flux difference splitting method, the Navier–Stokes equations and the realizable k-ε turbulence model equations describing the air flow are numerically solved. The evolution of the dynamic wave and ring vortex systems is effectively captured and analyzed. The influence of incident shock Mach number, perforated-plate porosity, and plate number on the propagation and attenuation of the shock wave was studied by using pressure- and entropy-based attenuation rates. The results indicate that the reflection, diffraction, transmission, and interference behaviors of the leading shock wave and the superimposed effects due to the trailing secondary shock wave are the main reasons that cause the intensity of the leading shock wave to experience a complex process consisting of attenuation, local enhancement, attenuation, enhancement, and attenuation. The reflected shock interactions with transmitted shock induced ring vortices and jets lead to the deformation and local intensification of the shock wave. The formation of nearly steady jets following the array of perforated plates is attributed to the generation of an oscillation chamber for the inside dynamic wave system between two perforated plates. The vorticity diffusion, merging and splitting of vortex cores dissipate the wave energy. Furthermore, the leading transmitted shock wave attenuates more significantly whereas the reflected shock wave from the first plate of the array attenuates less significantly as the shock Mach number increases. The increase in the porosity weakens the suppression effects on the leading shock wave while increases the attenuation rate of the reflected shock wave. The first perforated plate in the array plays a major role in the attenuation of the shock wave.     

Z箍缩动态黑腔形成过程MULTI程序一维数值模拟

Z箍缩动态黑腔能够高效地将Z箍缩丝阵等离子体动能转换为黑腔辐射能,为驱动惯性约束靶丸聚变提供高品质的X射线辐射场.利用一维双温多群辐射磁流体力学程序MULTI-IFE,研究了"聚龙一号"装置驱动电流条件下的Z箍缩动态黑腔形成基本物理过程.数值模拟结果表明,在动态黑腔形成过程中,辐射热波的传播速度比冲击波的传播速度更快,比冲击波更早到达泡沫中心,使中心区域的泡沫在冲击波到达前就已具有较高的辐射温度.对于"聚龙一号"装置动态黑腔实验0180发次采用的负载参数,辐射热波和冲击波在泡沫中的传播速度分别约为36.1 cm/μs和17.6 cm/μs,黑腔辐射温度在黑腔形成初期约80 eV,在冲击波到达泡沫中心前可达100 eV以上,丝阵等离子体外表面发射的X射线能量集中在1000 eV以下.本文给出了程序采用的计算模型、美国"土星"装置丝阵内爆计算结果和"聚龙一号"装置动态黑腔实验0180发次模拟结果.     

激波作用下SF6气泡界面演化和射流发展的数值模拟

 数值研究了平面激波冲击氮气环境中SF₆气泡界面的Richtmyer-Meshkov不稳定性,重点关注其中的激波聚焦及射流的产生和发展过程。在入射激波马赫数为1.23的情况下,给出了压力、密度、数值纹影和涡量等物理量的演化图像,定量分析了流场中压力最大值、密度最大值、射流速度、环量和斜压力矩随时间的变化关系。计算结果表明,平面激波冲击SF₆气泡过程有很强的聚能效应,在气泡内部靠近下游极点处发生激波近似理想聚焦和点爆炸现象,直接导致出现二次波系以及向下游运动的细长射流结构。相比入射激波,二次波系产生斜压力矩和涡量的能力要弱得多。     

非均匀炸药冲击起爆和起爆后的行为

 介绍并分析了Campbell 等人及其他作者研究非均匀炸药冲击起爆和起爆后行为所获得的实验结果,但不涉及其冲击起爆条件。足够强的冲击波进入非均匀炸药后,爆轰将瞬时（指不经过感应时间）且直接（指不经过其他过程,如爆燃）被引发;非均匀炸药起爆后,其中传播的自始至终是一个不断增长的爆轰波,直至发展为正常爆轰,整个过程都是爆轰的增长（新定义）过程。不存在由反应冲击波不断增长并转变为爆轰波的所谓向爆轰的增长。所谓向爆轰的增长,实际上是爆轰的增长（按新定义）的初期;Craig原定义的爆轰的增长,实际上是爆轰的增长（按新定义）的后期;而所谓反应冲击波,实际上是增长中的初期爆轰波。爆轰的增长（按新定义）是所有猛炸药的特性,炸药反应不充分并逐渐趋于充分是爆轰的增长的化学机制。     

Shock wave-inertial microbubble interaction: methodology, physical characterization, and bioeffect study

A method of generating in situ shock wave-inertial microbubble interaction by a modified electrohydraulic shock wave lithotripter is proposed and tested in vitro. An annular brass ellipsoidal reflector (thickness = 28 mm) that can be mounted on the aperture rim of a Dornier XL-1 lithotripter was designed and fabricated. This ring reflector shares the same foci with the XL-1 reflector, but is 15 mm short in major axis. Thus, a small portion of the spherical shock wave, generated by a spark discharge at the first focus (F1) of the reflector, is reflected and diffracted by the ring reflector, producing a weak shock wave approximately 8.5 microseconds in front of the lithotripter pulse. Based on the configuration of the ring reflector (different combinations of six identical segments), the peak negative pressure of the preceding weak shock wave at the second focus (F2) can be adjusted from -0.96 to -1.91 MPa, at an output voltage of 25 kV. The preceding shock wave induces inertial microbubbles, most of which expand to a maximum size of 100-200 microns, with a few expanding up to 400 microns before being collapsed in situ by the ensuing lithotripter pulse. Physical characterizations utilizing polyvinylidene difluoride (PVDF) membrane hydrophone, high-speed shadowgraph imaging, and passive cavitation detection have shown strong secondary shock wave emission immediately following the propagating lithotripter shock front, and microjet formation along the wave propagation direction. Using the modified reflector, injury to mouse lymphoid cells is significantly increased at high exposure (up to 50% with shock number > 100). With optimal pulse combination, the maximum efficiency of shock wave-induced membrane permeabilization can be enhanced substantially (up to 91%), achieved at a low exposure of 50 shocks. These results suggest that shock wave-inertial microbubble interaction may be used selectively to either enhance the efficiency of shock wave-mediated macromolecule delivery at low exposure or tissue destruction at high exposure.     

激波作用不同椭圆氦气柱过程中流动混合研究

在激波与气柱相互作用问题中,压力与密度间断不平行产生的斜压涡量会引起流动的不稳定性,从而促进物质间的混合.本文基于双通量模型,结合五阶加权基本无振荡(WENO)格式,求解多组分二维Navier-Stokes方程,分析激波作用面积相同结构不同的椭圆气柱所致的流动和混合.数值结果清晰地显示了激波诱导Richtmyer-Meshkov不稳定性引起的气柱界面变形和波系演化.同时定量地从界面运动、界面结构参数变化(长度和高度)、气柱体积压缩率、环量及混合率等角度分析激波诱导的流动混合机制,研究椭圆几何构型对氦气混合过程的影响.结果表明,界面及相关参数的演化与气柱初始形状密切相关.当激波沿椭圆长轴作用于气柱时,气柱前端出现空气射流结构,且射流不断增长并渗透到下游界面,致使气柱分离成两个独立涡团,离心率越大,射流发展越快;同时激波作用气柱后在界面处产生不规则反射现象.圆形气柱界面演化与这种作用情形类似.当激波沿椭圆短轴作用于气柱时,界面上游出现类平面结构,随后平面上下缘处产生涡旋,主导流动发展,激波在界面作用产生规则反射,离心率越大,这些现象越明显.界面高度、长度、体积压缩率也因此有所差异.对界面演化、环量和混合率的综合分析表明,激波沿长轴作用于气柱且离心率较大时,流动发展较快,不稳定性导致的流动越复杂,越有利于氦气与环境介质的混合.     

Color Ultrafast Photography Study of a High-Current Discharge Initiated by Exploding a Short Wire in the Atmosphere

In this paper, we describe the ultrafast slit-scanning photography of light phenomena arising at the stage of the electrical breakdown in the atmosphere and in a channel initiated by a wire explosion. The characteristics and features of the photographic registration of the SFR-1 camera using color and spectrozonal negative aerial films sensitive in the visible and near infrared regions of the optical spectrum are presented. Using photograms in color, we were able to visualize small-scale inhomogeneities, a current-carrying channel and its “shell,” emissions of matter, and shock waves in the plasma. This technique made it possible to record a secondary shock wave in the plasma on the discharge axis during the explosion of short (length ≈15 mm) copper and nichrome wires with a diameter of 100–150 μm at a stored energy of ~100 J. At the moment of cumulation in the plasma, a “hot spot” with a minimum size of ~0.5 mm was formed on the discharge axis. It may be the source of narrow-beam coherent radiation with X-ray photon energies of 10–30 keV. The cumulation mechanism is discussed based on the regularities of the z-pinch phenomenon in a discharge initiated by a wire explosion.     

Shock wave-inertial microbubble interaction: a theoretical study based on the Gilmore formulation for bubble dynamics

The Gilmore formulation for bubble dynamics coupled with zeroth-order gas diffusion were used to investigate theoretically the cavitation activity produced by a modified XL-1 lithotripter [J. Acoust. Soc. Am. 105, 1997-2009 (1999)]. The model calculation confirms many of the basic features in bubble dynamics observed experimentally, in particular the strong secondary shock wave emission generated by in situ lithotripter shock wave-inertial microbubble interaction. In addition, shock wave-inertial microbubble interaction produced by a Dornier HM-3, the most commonly used clinical lithotripter, was evaluated. It was shown that the forced collapse of inertial microbubbles with strong secondary shock wave emission could be produced consistently, provided that an appropriate preceding shock wave and interpulse delay were used. Further, it was demonstrated that truncation of the tensile stress of the lithotripter shock wave could significantly reduce the large expansion of the bubble following shock wave-inertial microbubble interaction, which may alleviate the risk for vascular injury during shock wave exposure.     

Visualization of shock-vortex interaction radiating acoustic waves

Unsteady compressible flow fields past a wedge and a cone, evolved by propagation and interaction of shock waves, slip lines, and vortices, are studied by shadowgraphs and holographic interferograms taken during the shock tube experiment. The supplementary numerical calculation also presented time-accurate solution of the shock wave physics which was essential to recognize the similarity and dissimilarity between the wedge and the conical flows. The decelerated shock detained by the vortex interacts with the small vortexlets along the slip layer, producing diverging acoustics: this phenomenon is more distinct in the case of wedge flow for a given shock Mach number. The decelerated shock penetrated through the vortex core constitutes a transmitted shock, which eventually merges with the diaphragm shock that bridges the vortex pair/vortex ring. This phenomenon became remarkably salient in the case of conical flow.     

Shock tube exit flow fields through thin membranes

The flow field characteristics that form for a shock wave propagating through a membrane-like termination at the exit of a shock tube are studied. The strength of the shock wave reflected back into the tube, as well as the strength of the shock wave transmitted, is examined. Six different materials are used, ranging from copper and aluminium foils to a variety of elastic and plastic sheets, in a few cases with different initial pressure differentials. High-speed shearing interferometry imaging is done of the external flow. Three principal characteristics are present in the transmitted flow: a diffracted shock wave, an expansion wave and a re-compression shock wave. It is found that the prominence of these features varies depending on the material type. For the later flow development and material rupture, there are a number of principal characteristics: small vortices, secondary shock waves, a vortex ring, oblique waves and a Mach disc.     

Phase change and stress wave in picosecond laser–material interaction with shock wave formation

When background gas is present in pulsed laser–material interaction, a shock wave down to the nanoscale will emerge. The background gas will affect the phase change and explosion in the target. This study is focused on the void dynamics and stress wave in a model material (argon crystal) under picosecond pulsed laser irradiation. Our results show that existence of ambient gas and the shock wave significantly suppresses the void formation and their lifetime. Void dynamics, including their growing rate, lifetime, and size under the influence of ambient gas are studied in detail. All the voids undergo an accelerating and decelerating process in the growth. The collapsing process is almost symmetrical to the growing process. Higher laser fluence is found to induce an obvious foamy structure. Stress wave formation and propagation, temperature contour, and target and gas atom number densities are studied to reveal the underlying physical processes. Although the interaction of the plume with ambient gas significantly suppresses the void formation and phase explosion, no obvious effect is found on the stress wave within the target. Very interestingly, secondary stress waves resulting from re-deposition of ablated atoms and void collapse are observed, although their magnitude is much smaller than the directly laser-induced stress wave.     

激光冲击强化“残余应力洞”的形成机制

利用ABAQUS有限元软件进行了单个圆形高斯光斑的激光冲击强化数值模拟,分析材料表面光斑中心区域形成的"残余应力洞"现象,并通过分析材料的动态力学响应特征揭示了"残余应力洞"的形成机制。结果表明:在冲击波加载时,光斑边界处会产生很强的剪切应力,形成向四周传播的表面稀疏波和向材料内部传播的剪切波。当稀疏波同时传播到光斑中心,发生相遇、汇聚,使材料产生急剧的上下位移过程,造成冲击波加载塑性变形后的二次塑性变形。二次塑性变形中形成了较大的剪切塑性应变,并降低了冲击波加载阶段产生的轴向和径向塑性应变,使残余压应力降低,从而形成"残余应力洞"。