Shock Wave Velocity Research Articles

Experimental study on the evolution of shock waves in the near field of cylindrical charge explosion

Cylindrical charges are frequently used in blasting engineering and military operations, often detonated very close to target structures. Obtaining the characteristics of the flow field distribution in the near field of explosions through experimentation consistently poses challenges. These challenges constrain efforts to elucidate the evolutionary processes of explosive shock waves. Employing the explosive schlieren experimental system yields a clearer depiction of the near field flow than previous studies, thereby facilitating a more comprehensive understanding of the characteristics of cylindrical charge explosions. The overall flow field adopts a droplet shape, with the interface of the detonation products and the shock wave front forming continuous concave shapes. The detonation process in cylindrical charges can be conceptualized as the detonation wave propelling the detonation products at supersonic speeds. The formation of explosive shock waves can be elucidated through the dynamics of oblique shock waves generated by two-dimensional curved surfaces. This mechanism influences the relationship between the deflection angles of detonation products and the shock wave angles resulting from different detonation velocities. Theoretical calculations of shock wave positions align with experimental images. Significant variations in shock wave velocities at different positions within cylindrical charges are observed, identifying patterns of these variations. The research findings clarify the formation mechanisms of explosive shock waves in cylindrical charges and deepen the understanding of shock wave evolution in the near field.

Shock response of two epoxy resins at up to 330 GPa pressure

Experimental studies of the shock wave properties of two epoxy resins with the same composition but different curing temperatures (160 and 200 °C) at up to 330 GPa pressure have been carried out. Laser interferometry was used to record particle velocity profiles at up to 73 GPa pressure while measuring the shock wave velocity. The release sound velocity was experimentally determined in the 3–73 GPa pressure range. Cumulative explosive shock wave generators were used to study the shock Hugoniot of epoxy resins at pressures above 100 GPa. It was shown that the shock compressibility data of both samples are approximated by a single shock Hugoniot within the experimental error. A kink on Hugoniot recorded close to 25 GPa pressure indicates a chemical decomposition in epoxy resin. Above this kink, a change in the shock wave front structure was recorded. Hugoniots of epoxy resin and unidirectional carbon/epoxy composite were compared at up to 370 GPa pressure.

Numerical study of cavitation shock wave emission in the thin liquid layer by power ultrasonic vibratory machining

In the field of power ultrasonic vibration processing, the thin liquid layer nestled between the tool head and the material serves as a hotbed for cavitation shock wave emissions that significantly affect the material's surface. The precise manipulation of these emissions presents a formidable challenge, stemming from a historical deficit in the quantitative analysis of both the ultrasonic enhancement effect and the shock wave intensity within this niche environment. Our study addresses this gap by innovatively modifying the Gilmore-Akulichev equation, laying the groundwork for a sophisticated bubble dynamics model and a pioneering shock wave propagation model tailored to the thin liquid layer domain. Firstly, our study investigated the ultrasound enhancement effect under various parameters of thin liquid layers, revealing an amplification of ultrasound pressure in the thin liquid layer area by up to 7.47 times. The mathematical model was solved using the sixth-order Runge–Kutta method to examine shock wave velocity and pressure under different conditions. our study identified that geometric parameters of the tool head, thin liquid layer thickness, ultrasonic frequency, and initial bubble radius all significantly influenced shock wave emission. At an ultrasonic frequency of 60 kHz, the shock wave pressure at the measurement point exhibited a brief decrease from 182.6 to 179.5 MPa during an increase. Furthermore, rapid attenuation of the shock wave was found within the range of R0-3R0 from the bubble wall. This research model aims to enhance power ultrasonic vibration processing technology, and provide theoretical support for applications in related fields.

Experimental and numerical simulation of the attenuation effect of blast shock waves in tunnels at different altitudes

Traffic engineering such as tunnels in various altitudinal gradient zone are at risk of accidental explosion, which can damage personnel and equipment. Accurate prediction of the distribution pattern of explosive loads and shock wave propagation process in semi-enclosed structures at various altitude environment is key research focus in the fields of explosion shock and fluid dynamics. The effect of altitude on the propagation of shock waves in tunnels was investigated by conducting explosion test and numerical simulation. Based on the experimental and numerical simulation results, a prediction model for the attenuation of the peak overpressure of tunnel shock waves at different altitudes was established. The results showed that the peak overpressure decreased at the same measurement points in the tunnel entrance under the high altitude condition. In contrast, an increase in altitude accelerated the propagation speed of the shock wave in the tunnel. The average error between the peak shock wave overpressure obtained using the overpressure prediction formula and the measured test data was less than 15%, the average error between the the propagation velocity of shock waves predicted values and the test data is less than 10%. The method can effectively predict the overpressure attenuation of blast wave in tunnel at various altitudes.

Blast wave interaction during explosive detonation in a variable cross-sectional charge

Purpose. The research aims to assess the efficiency and performance of solid media destruction in directed blasting of a charge with variable cross-sectional shape. Methods. Numerical modelling of the blast wave interaction process is performed using the finite element method based on the Euler-Lagrange algorithm. The Johns-Wilkins-Lee equation of state is used to determine the pressure-volume dependences of medium destruction. Assessment of solid medium destruction mechanism during a directed blasting of a charge with variable cross-sectional shape is carried out based on polarization-optical method on models made of optically active material. Findings. Experimental studies of solid medium destruction by the action of directed blasting with a variable cross-section charge made it possible to determine the direction of blast wave propagation and its amplitude in stress wave, influencing the intensity of radial crack network formation in superposition areas, and directed perpendicularly to explosive cavity. At the same time, the average peak pressure in collision zone of two shock waves in centre of spherical cavity is approximately 1.48 and 1.84 times higher than that in weakly blast-loaded areas. It has been found that when two shock waves collide and superimpose on each other, the intensity of their impact increases. Moreover, the shock wave velocity in collision zone is higher than that of the radial shock wave. Originality. It has been determined that the maximum pressure values on the explosive cavity wall at the initiation points sharply increase and then gradually stabilize as the blast stress waves propagate and have an arbitrary distribution pattern. Three areas should be considered: not superimposed, weakly superimposed, and strongly superimposed. At each point of detonation, the pressure on explosive cavity wall will be minimal, while in the charge centre in the spherical insert zone, on the contrary, it will be maximal. In this case, the pressure in the central superposition area is about 2.84 times greater than at the initiation ends, and the nature of distribution changes according to a linear dependence. Practical implications. The performed research findings can serve as a basis for development of effective parameters of resource-saving methods for stripping hard rocks of complex structure in the conditions of ore mines.

Improving prediction accuracy of laser-induced shock wave velocity prediction using neural networks

The velocity of laser-induced shock waves affects the efficiency and efficacy of laser-based processes. The ability to accurately estimate shock wave velocity is critical for optimizing experimental combinations, creating laser-based systems, and assuring desired results. Traditional approaches to predict shock wave velocity involve empirical equations and analytical models based on simplified assumptions. However, these methods often lack accuracy and fail to capture the complex dynamics of laser-matter interactions. To overcome these limitations, we used a combination of an artificial neural network and a genetic algorithm to predict shock wave velocity. In this method, the neural network structure is dynamically designed. The optimization method does this by modifying the neural network's weights and figuring out the network's structure on our behalf. Based on the findings, our suggested technique worked very well; it surpassed other comparison methods by achieving the lowest average errors in terms of RMSE and MAE, which are 4.38 and 3.74, respectively. Moreover, the analysis has shown that our proposed method has a high level of reliability in predicting impulsive wave velocity using a neural network.

Analysis of shock wave propagation in the atmosphere through generated sound waves

The paper presents the results of studies of shock waves propagating in the atmosphere. These results were obtained by recording sound waves generated by a shock wave. It has been shown that it takes approximately two seconds for a shock wave to form. The frequency of the sound wave generated by the shock wave depends on the speed of propagation of the shock wave. It was found that the shock wave accelerates as it propagates upward. This phenomenon can be used as a method for determining the velocities of shock waves or supersonic moving bodies.

Electrical properties of different materials studied by sub-microsecond underwater electrical explosions of single wires

Results of an experimental research and one-dimensional hydrodynamical simulations of critically damped sub-microsecond timescale underwater electrical explosions of wires made of 12 different materials are presented. Using current and voltage waveforms, streak shadow images of the shocks generated in water and wire expansion obtained by one-dimensional hydrodynamic simulations, the maximal values of the energy density, energy density deposition rates, and specific action integrals were determined. It is shown that for all study materials, the deposited energy density significantly exceeds the energy density required for the solid–liquid phase transition but is substantially smaller to induce a full liquid–vapor phase transition of the wire. At the time when the maximal value of the deposited power is realized, the deposited energy densities were found to be larger than the atomization energy for all materials. Estimates of the plasma parameters show that the explosion of the wires can be characterized by a high resistance and lowly ionized weakly coupled plasma. Three groups of materials were distinguished by either decrease, plateau, or increase in the resistance after the maximum of the deposited power. It was confirmed that the observed maximum Planckian temperature for all wire material does not exceed 6000 K due to the “bath” effect and that there is a correlation between the wire radial expansion and the strong shock wave velocities.

Mechanical properties and impact behavior of frozen clay: Insights from static mechanical tests, fly-plate tests, and split-Hopkinson pressure bar analysis

The presence of frozen clay as a natural protective material makes it a crucial layer of defense against potential impacts in various engineering projects. Studying the strength and deformation characteristics of frozen clay is, therefore, particularly important. In this study, static mechanical tests, the Hopkinson impact test, and fly-plate tests were conducted on frozen clay to identify its mechanical properties under an impact load. The uniaxial compression strength, flexural strength, and elastic modulus displayed a linear increase with the change in temperature. The fractal dimension was used to describe the failure characteristics of frozen clay, yielding values ranging from 1.5691 to 1.8785. At the same temperature, the fractal dimension exhibited a strain rate effect as the strain rate increased. A light gas gun system was then used to conduct fly-plate tests on frozen clay at varying temperatures (−3, −20 °C, and ordinary temperature). The impact process was meticulously analyzed, considering factors such as shock wave velocity, particle velocity behind the shock wave, impact pressure, and volume strain. Moreover, our investigation plotted the D–u (volume strain–particle velocity) and P–μ shock (impact pressure–shock wave velocity) adiabatic curves. Notably, we observed that samples with a higher initial strength exhibited an increased resistance to compression under an identical initial density and moisture content, resulting in a discernible leftward shift of the P–μ curve. The results provide a theoretical basis and technical support for similar projects in the future.

Inverse design of the radiation temperature for indirect laser-driven equation-of-state measurement

The theoretical design for the time profile of radiation temperature plays an important role in indirect laser-driven equation-of-state measurement, which severely relies on a large number of radiation hydrodynamic simulations. In this work, we provide a concise data-driven method for optimizing the radiation temperature profile, which combines a time-varying Volterra model with an improvement achieved by data generation via radiation hydrodynamic simulations utilizing random perturbations in a skew normal distribution as inputs. We find that the time-varying Volterra model can be used to investigate the time-dependent relationship between the radiation temperature and the key physical quantities of interest, such as shock-wave velocity and ablation drive pressure. With this method, we realize the inverse designs of the radiation temperature profiles for planar dynamic shock and ramp compressions according to the desired shock-wave velocity and drive pressure, respectively, which shows the advantage of practical application in experiments.

Characterization of the Dynamic Behavior of Multinanobubble System under Shock Wave Influence.

Shockwave-induced changes in nanobubbles cause cavitation erosion and membrane damage but can also be applied to biocarrier transport. Currently, research focuses on single nanobubbles; however, in reality, nanobubbles usually appear as a multibubble system. Therefore, this study proposes a method based on cutting and replicating to construct a multibubble model. This method can be widely applied to molecular dynamics (MD) models and enhance the customization capabilities of MD models. The dynamic behavior of a multinanobubble system with different numbers and arrangements of nanobubbles is investigated with the MD method under the influence of shock waves in a liquid argon system. The study also explores the range of influence between nanobubbles. The results show that in the case of two nanobubbles, when the distance between the bubbles is constant, the smaller the angle between the direction of the shock wave and the line connecting the bubbles, the greater is the influence between nanobubbles, and the moment of collapse of the nanobubbles farther away from the shock wave is slower. When three nanobubbles are arranged with a right offset, after the first bubble collapses, the effect on the other two bubbles is similar to the changes in bubbles when the angle of arrangement is 30° or 60°. Under a different arrangement, the change of shock wave velocity on the nanobubble size only affects its collapse time and contraction collapse rate. When the shock wave with a radian of about 2.87 or greater than 2.87 touches the bubbles, the collapse of the second nanobubble will not be affected.

Laser-induced shock wave two-dimensional extraction technology based on a beam deflection method

This paper introduces a method for analyzing the spatiotemporal progression of laser-induced shock waves using the beam deflection technique. This method allows for the accurate measurement of the shock wave evolution and can replace high-speed cameras. The results demonstrate the detection signals at various distances and energies, as well as the extraction and reconstruction of the shock wave velocities and propagation trajectories. The characteristic velocities of the shock waves propagating in air from various metals and energetic materials were measured and compared with the results obtained from high-speed cameras. The study also predicts the macroscopic detonation velocity of energetic materials based on the characteristic velocity. Overall, this approach offers a reliable and cost-effective method for studying the shock waves and has potential applications in various fields.

Shock Response of Polymethyl Methacrylate (PMMA) Under Explosive Loading

AbstractThe dynamic response of polymethyl methacrylate (PMMA) is well understood for one-dimensional planar impact shocks, but limited research has been performed on the response of PMMA under spherical shock loading. In this work, the shock decay of an explosively-driven shock wave into PMMA was experimentally measured. PMMA cubes of various geometries were explosively loaded with an RP-80 detonator to produce the explosive shock wave. High-speed schlieren imaging was implemented to measure the explosively-driven shock wave velocity throughout the PMMA cubes. Photon Doppler velocimetry (PDV) was used to measure the particle velocity imparted by the shock wave at the surface of the cubes. The material shock response was studied at distances from 21.91 to 133.3 mm from the explosive source. The particle velocity history measured by PDV was compared to the wave profile visualized in the high-speed images. The shock wave pulse amplitude decreased with increased distance from the source. The conducted experiments extend the PMMA shock Hugoniot relating to the lower shock and particle velocity regime.

Investigating built environment and traffic flow impact on crash frequency in urban road networks

While numerous studies have examined the factors that influence crash occurrence, there remains a gap in understanding the intricate relationship between built environment, traffic flow, and crash occurrences across different spatial units. This study explores how built environment attributes, and dynamic traffic flow characteristics affect crash frequency by focusing on proposed traffic density-based zones (TDZs). Utilizing a comprehensive dataset from Greater Melbourne, Australia, this research emphasizes on the dynamic traffic flow variables and insights from the Macroscopic Fundamental Diagram model, considering parameters such as shockwave velocity and congestion index. The association between the potential influencing factors and crash frequency is examined using a random parameter negative binomial regression model. Results indicate that the data segmentation based on TDZs is instrumental in establishing a more refined crash model compared to traditional planning-based zones, as demonstrated by improved goodness-of-fit measures. Factors including density (e.g., employment density), network design (e.g., road density and highway density), land use diversity (e.g., job-housing balance and land use mixture), and public transit accessibility (e.g., bus route density) are significantly associated with crash occurrence. Furthermore, the unobserved heterogeneity effects of the shockwave velocity and congestion index on crashes are revealed. The study highlights the significance of incorporating dynamic traffic flow variables in understanding crash frequency variations across different spatial units. These findings can inform optimal real-time traffic monitoring, environmental design, and road safety management strategies to mitigate crash risks.

Ultrafast Process Characterization of Laser-Induced Damage in Fused Silica Using Pump-Probe Shadow Imaging Techniques.

This study delves into the intricate dynamics of laser-induced damage in fused silica using a time-resolved pump-probe (TRPP) shadowgraph. Three typical ultra-fast processes, laser-induced plasma evolution, shockwave propagation and material fracture splashing, were quantitatively investigated. The results indicate that the diameter of plasma is proportional to the pulse laser energy and increases linearly during the pulse laser duration with an expansion rate of approximately 6 km/s. The maximum shockwave velocity on the air side is 9 km/s, occurring at the end of the pulse duration, and then rapidly decreases due to air resistance, reaching approximately 1 km/s around a 300 ns delay. After hundreds of nanoseconds, there is a distinct particle splashing phenomenon, with the splashing particle speed distribution ranging from 0.15 km/s to 2.0 km/s. The particle sizes of the splashing particles range from 4 μm to 15 μm. Additionally, the smaller the delay, the faster the speed of the splashing particles. Overall, TRPP technology provides crucial insights into the temporal evolution of laser-induced damage in fused silica, contributing to a comprehensive understanding essential for optimizing the performance and safety of laser systems.

Shockwave velocimetry using wave-based image processing to measure anisotropic shock emission

Noninvasive optical measurements of the shockwave propagation velocity using multiple pulse illumination allow deducing the shockwave pressure amplitude through Hugoniot relations and an appropriate equation of state of the medium. This technique is particularly useful for spatially resolved measurements near the shockwave emission site. Due to diffraction, however, a shockwave front can significantly change its morphology, rendering precise velocity measurements non-trivial. As solution we propose a wave front evolution (WaFE) velocimetry technique, which applies Huygens principle. We take a shadowgraph of the wave front at subsequent times as initial condition for the acoustic Helmholtz equation and numerically propagate the fronts in time. From the instance of time, when two subsequently taken wave front shadows numerically interfere and form one sharp wave front, the local shock velocity is obtained and the local shock pressure amplitude measured. With artificial test images, it is shown that this technique has excellent sub-pixel accuracy, robustness to noise, and can work with low contrast images and even overlapping and interfering wave fronts. The software is made available freely and can be applied to general shock front velocity measurements. We apply WaFE to determine the anisotropic shockwave emission from an elongated laser-induced plasma in water from shadowgraphs of the shockwave front imaged four times onto the same camera frame using multiple pulse illumination at a repetition rate of 60 MHz. The direction dependence of attenuation of the shockwave pressure amplitude is measured at distances of 50–300 μm to the plasma.

The impact of low-velocity shock waves on the dynamic behaviour characteristics of nanobubbles.

Low-velocity shock wave-induced contraction and expansion of nanobubbles can be applied to biocarriers and microfluidic systems. Although experiments have been conducted to study the application effects, the dynamic behavior characteristics of nanobubbles remain unexplored. In this work, we utilize molecular dynamics (MD) simulations to investigate the dynamic behavior characteristics of nanobubbles influenced by low-velocity shock waves in a liquid argon system. The DBSCAN (Density-Based Spatial Clustering of Applications with Noise) machine learning method is used to calculate the equivalent radius of nanobubbles. Two statistical methods are then utilized to predict the time series changes in the equivalent radius of nanobubbles without rebound shock waves. The piston velocity is analyzed using the bisection method to obtain the critical impact states of the nanobubble. The results show that at the low velocity shock wave (piston velocity of 0.1 km s-1), the shock wave pressure is small, the non-vacuum nanobubbles contract and expand in a circular shape, and the gas particles inside the bubble are not dispersed. In contrast, the vacuum nanobubbles collapse directly. As the shock wave rebounds upon impact, it triggers periodic contraction and expansion of the nanobubbles. The predictions indicate that the equivalent radius will vary within a small range according to the pre-predicted values in the absence of the rebound shock wave. Nanobubbles are present in four critical impact states: dispersed gaps, multiple smaller bubbles, two split bubbles, and a concave bubble.

The role of κ‐deformed Kaniadakis distributed electrons on the dust ion‐acoustic waves in charge‐varying dusty plasma

AbstractThe one‐dimensional dynamics of nonlinear electrostatic dust‐ion acoustic waves (DIAWs), in a plasma containing Kaniadakis‐distributed electrons, fluid ions, and variable‐charge dust grains has been investigated. The κ−generalized electron charging current is derived based on the orbit‐motion limited approach. The reductive perturbation technique has been used to derive a KdV Burger equation of propagation. We have highlighted the dependence between the values of the dust particles charge Qd0 and the validity domain of the Kaniadakis parameter κ, where the Burger term C has positive. By integrating the KdV Burger equation numerically and assuming a laboratory dusty plasma, a transition from oscillatory to monotonic shock waves is observed with increasing values of κ. We have also found that the phase velocity of DIA shock waves decreases, and the dissipation can prevail over that of dispersion as the electrons move further away from thermal equilibrium (κ is increased). The coexistence of compressive and rarefactive of monotonic shock profile is observed. It is found that an increase in κ leads to an increase of the compressive monotonic shock wave and decrease the rarefactive monotonic shock wave.

Effect of dimensionless vent ratio on the flame-shock wave evolution dynamics of blended LPG/DME gas explosion venting

To reveal the flow field-shock wave-flame coupling evolution law of LPG/DME blended gas explosion venting under different dimensionless vent ratio(KV), the numerical model for multi-blended gas explosion venting was established and research of blended gas explosion venting at KV = 0.02–0.21 was carried out based on experimental verification. The results show that, after the vent opens, the turbulence intensity increases, the velocity gradient of the flow field increases, the maximum explosion wind velocity(Vmax) of the outfield varies linearly with KV, and the flame of “mushroom shape” is elongated, and “break flame” occurs under the action of the relative motion of “primary and secondary vortex group”. When KV＞0.02, the flame propagation velocity increases locally because of the thermal-mass diffusion instability, and with increasing KV, the average flame propagation velocity decreases from 258.29 m/s to 226.58 m/s. In the process of explosion venting, a typical shock wave fluctuating pressure waveform curve is formed in the outfield. Under the combined action of the infield pressure venting effect and the outfield turbulence intensity, the maximum explosion overpressure in the outfield (Pmax) reaches maximum 423.97 kPa at KV = 0.07. Meanwhile, the larger the KV, the greater the shock wave velocity due to the lead role of the infield explosion venting shock wave and outfield flame compression shock wave on the shock wave propagation velocity.

Исследование течения, создаваемого скользящим разрядом, с помощью методов компьютерного зрения

A quantitative study has been made of the flow with shock waves generated in air by a sliding surface discharge lasting less than one microsecond. The high-speed flow was visualized using the shadowgraph method, the process was recorded at a rate of 124 000 frames/s, the exposure time was 1 μs. The aim of this work is to study the dynamics of the two discontinuities: the cylindrical shock wave and the contact surface generated by the discharge. Each experiment allowed several hundred images to be taken of a short-lived gas-dynamic process lasting up to 1 ms. A YOLOv8 convolutional neural network was trained and used to determine the positions of the discontinuities. A data set of 984 markups was labeled. The model on the mAP50 metric achieved 0.887 and the mAP50-95 was 0.557. The model was used to automatically measure the vertical dimensions of the contact discontinuity. It expands at times up to 0.4 - 0.8 ms to a vertical size of 5 - 11 mm. The x-t plots and the velocities of the cylindrical shock waves were measured. It is shown that up to 1 ms after the discharge, the flow development is due to the blast wind motion behind the shock wave. It is shown that the use of computer vision can significantly speed up the analysis of high-speed flow visualizations and the extraction of quantitative information.