Shock Wave Front Research Articles

Proton acceleration in plasma turbulence driven by high-energy lepton jets

Abstract The interaction of high energy lepton jets composed of electrons and positrons with background electron–proton plasma is investigated numerically based upon particle-in-cell simulation, focusing on the acceleration processes of background protons due to the development of electromagnetic turbulence. Such interaction may be found in the universe when energetic lepton jets propagate in the interstellar media. When such a jet is injected into the background plasma, the Weibel instability is excited quickly, which leads to the development of plasma turbulence into the nonlinear stage. The turbulent electric and magnetic fields accelerate plasma particles via the Fermi II type acceleration, where the maximum energy of both electrons and protons can be accelerated to much higher than that of the incident jet particles. Because of background plasma acceleration, a collisionless electrostatic shock wave is formed, where some pre-accelerated protons are further accelerated when passing through the shock wave front. Dependence of proton acceleration on the beam-plasma density ratio and beam energy is investigated. For a given background plasma density, the maximum proton energy generally increases both with the density and kinetic energy of the injected jet. Moreover, for a homogeneous background plasma, the proton acceleration via both turbulent fields and collisionless shocks is found to be significant. In the case of an inhomogeneous plasma, the proton acceleration in the plasma turbulence is dominant. Our studies illustrate a scenario where protons from background plasma can be accelerated successively by the turbulent fields and collisionless shocks.

Experimental study on attenuation effect of liquid viscosity on shockwaves of cavitation bubbles collapse

How to precisely control and efficiently utilize the physical processes such as high temperature, high pressure, and shockwaves during the collapse of cavitation bubbles is a focal concern in the field of cavitation applications. The viscosity change of the liquid will affect the bubble dynamics in turn, and further affect the precise control of intensity of cavitation field. This study used high-speed photography technology and schlieren optical path system to observe the spatiotemporal evolution of shockwaves in liquid with different viscosities. It was found that as the viscosity of the liquid increased, the wave front of the collapse shockwave of the cavitation bubble gradually thickened. Furthermore, a high-frequency pressure testing system was used to quantitatively analyze the influence of viscosity on the intensity of the shockwave. It was found that the pressure peak of the shockwave in different viscous liquid was proportional to Lb (L represented the distance between the center of bubble and the sensor measuring point), and the larger the viscosity was, the smaller the value of b was. Through in-depth analysis, it was found that as the viscosity of the liquid increased, the proportion of the shockwave energy of first bubble collapse to the maximal mechanical energy of bubble gradually decreased. The proportion of the mechanical energy of rebounding bubble to the maximal mechanical energy of bubble gradually increased. These new findings have an important theoretical significance for the efficient utilization of ultrasonic cavitation.

Damage behavior of the surrounding medium caused by different blast-induced dynamic and static loadings

The combination of the dynamic action of explosion stress wave and the static action of explosion gas expansion is what primarily drives rock-breaking blasting, and their different intensities lead to different blasting effects. Therefore, experiments were conducted to explore the damage behaviour of the surrounding medium caused by different blast-induced dynamic and static loadings using digital laser Schlieren and digital laser dynamic caustic systems. The results were as follows: In the air medium, when the explosives detonated, TATP (triacetone triperoxide) had the strongest quasi-static effect, and the highest velocity explosion shock wave and products compared with those of NHN (nickel hydrazine nitrate) and DDNP (diazodinitrophenol). Part of the TATP explosion product was ahead of the shock wave front, whereas for NHN and DDNP, their explosion products were behind the propagation of the shock wave front. The highest pressure peak of TATP was the result of the combined action of the shock wave and explosion product impacting the sensor. In the solid medium, NHN, the explosive with the strongest dynamic load, had the strongest stress wave and highest stress intensity factor peak of the main crack. The velocity reached the peak value at the early stage of propagation, and the dynamic action had a dominant effect on the cracking of the crack; the stronger the dynamic action, the higher the crack initiation energy. The stress intensity factor and propagation rate of the main crack of TATP with the strongest quasi-static effect decreased slower than those of the other explosives; therefore, the main crack propagation time and length were the longest. After cracking, the quasi-static effect dominated the process of crack propagation, where the stronger the quasi-static effect, the longer duration of the driving force of crack propagation. The research results provide a basis for customising explosives for different functions in practical engineering and realising the fine and efficient utilisation of explosion energy.

Understanding the unusual life of the Cartwheel galaxy using stellar populations

Collisional ring galaxies (RiGs) are the result of the impact between two galaxies, with one of them passing close to the centre of the other, piercing its gaseous and stellar disc. In this framework, the impact generates a shock wave front that propagates within the disc of the target galaxy soon after the encounter, producing a characteristic expanding ring-shaped structure. RiGs represent one of the most extreme environments in which we can study the physical properties of galaxies and the transformations they undergo during collisions. The paradigm RiG is the Cartwheel galaxy at z = 0.03. This galaxy has been the object of both theoretical and observational studies, but the details of the mechanisms that lead to its peculiar morphology of double rings with connecting spokes and to its physical properties are still far from clear. To shed light on the history of the Cartwheel galaxy, we performed a spatially resolved analysis as a function of galactocentric distance, exploiting spectroscopic data from VLT/MUSE observations combined with photometric data covering a large wavelength range, from the UV GALEX to the IR JWST/MIRI. Using full-index fitting of the stellar spectra, an analysis of the nebular emission, and joint full spectral and photometry fitting, we derived stellar ages, gas and stellar metallicities, and star formation histories (SFHs) in four spatially distinct regions of the galaxy. We find that, apart from the peculiar morphology, a large fraction of the Cartwheel galaxy is not affected by the recent impact from the companion bullet, and retains the characteristics of a typical spiral galaxy. On the contrary, the outer ring is strongly affected by the recent impact, and is completely dominated by stars formed not earlier than ∼400 Myr ago. Our picture suggests that the collision shock wave, while moving forward to the external region of the galaxy, drags the already formed stars, sweeping the inner areas outwards, as proposed by recent collision models. At the same time, the ages found in the external ring are older than the predicted timescale of the ring expansion after the collision.

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.

General backreaction force of cosmological bubble expansion

The gravitational-wave energy-density spectra from cosmological first-order phase transitions crucially depend on the terminal wall velocity of asymptotic bubble expansion when the driving force from the effective potential difference is gradually balanced by the backreaction force from the thermal plasma. Much attention has previously focused on the backreaction force acting on the bubble wall alone but overlooked the backreaction forces on the sound shell and shock-wave front, if any, which have been both numerically and analytically accomplished in our previous studies but only for a bag equation of state. In this paper, we will generalize the backreaction force on bubble expansion beyond the simple bag model. Published by the American Physical Society 2024

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.

Eulerian–Lagrangian multiscale numerical analysis of multimodal partial shedding dynamics

The objective of this paper is to investigate the multimodal partial shedding dynamics from a multiscale perspective of cloud cavitating flows under two distinct cavity shedding mechanisms, namely the re-entrant jet mechanism and the shock wave propagation mechanism. A two-way Eulerian–Lagrangian coupling algorithm is applied to capture the multiscale vapor topologies from microbubble to large-scale cavities. The large-scale cavity evolution is solved through large eddy simulations (LES) with the volume of fraction (VOF) method in Eulerian frame. The sub-grid microbubbles are tracked in Lagrangian frame based on the discrete bubble model (DBM) method. The predictions agree well with experimental observation of the periodical cavity evolution and microbubble dynamics under both the re-entrant jet mechanism and shock wave mechanism around a NACA66 hydrofoil. The numerical simulation provides detailed analysis of the cavitating turbulent flow on the microbubble behavior with emphasis on the spatial-temporal distribution characteristics of microbubbles. The results show that the number and mean size of microbubbles in the cavitation region increase gradually with the growth of attached sheet cavity, development of re-entrant jet and collapse of largescale cavity for both cavitation patterns. Meanwhile, microbubbles are mainly distributed on the largescale interfaces where have high value of vorticity and turbulent kinetic energy under the effect of re-entrant jet and vortex structures. And the probability density functions (PDFs) of microbubble exhibit gamma distributions with a dominant peak at approximately 50 μm for both shedding mechanisms. However, the shock wave formation and propagation process only occurs in the final stage of cavitating flow under shock wave mechanism causing the condensation of vapor and the decrease of the number and mean size of microbubbles. Moreover, the microbubbles are uniformly distributed along the streamwise and vertical directions behind shock wave front.

The shock wave evolution of the laser-induced breakdown under non-spherical symmetry using a 1064 nm nanosecond laser

In acoustics-related interdisciplinary areas, the shock wave of laser-induced breakdown has garnered significant attention. However, research on the propagation of shock waves in non-spherical symmetry is insufficient in both theoretical and experimental aspects. This paper aims to thoroughly study the evolution of underwater shock wave directivity by employing the method of shadowgraph. The shock wave front is determined by the dark fringes in the shadowgraph image and the normal propagation speed of the shock wave is calculated using Huygens principle. Subsequently, normal propagation speed is converted to pressure in different directions by employing the equation of state of the medium. It has been found that the spherical plasma produces an isotropic shock wave, whereas filamentary plasma generates a highly anisotropic one. To evaluate the anisotropic property of the shock wave, we introduce pressure directivity, which is defined as the pressure at any direction normalized by the maximum value. The temporal evolution of shock wave pressure directivity is obtained based on the shadowgraph images. In the case of filamentary plasma, there is a sudden transition of the pressure directivity in the axial from 40 ns to 165 ns, after which the pressure directivity is consistent with the hydrophone measurement. Based on the moving breakdown model of the plasma and the superposition principle, we propose a theoretical model to explain the experimental result of the pressure directivity. The outcome of our model exhibits considerable consistency with the experiment.

How Common can the Origin of Clusters Cha I, Cha II, ϵ Cha and η Cha be?

The trajectories of the groupings η Cha, ϵ Cha and Cha I, Cha II, constructed backward in time, have been studied. We concluded that the hypothesis about the joint formation of all four of these clusters from one molecular cloud cannot be completely excluded. However, 10–15 Myr ago, all these four clusters were located at approximately the same height above the plane of the Galaxy. Thus, the gas-dust clouds from which all these four clusters were formed were located on one broad front. It is possible that the appearance of the Cha I, Cha II, ϵ Cha and η Cha clusters may be associated with the impact on such a front of shock waves formed after supernova explosions in the Scorpius–Centaurus association. New estimates of the kinematic ages of the clusters Cha I and Cha II are obtained as 0.12 ± 0.19 and 0.05 ± 0.15, respectively. It is shown that the minimum size of the Cha I-north and Cha I-south clusters corresponded to the time of 0.55 ± 0.24 Myr and 0.04 ± 0.18 Myr ago, respectively, and approximately 1.5 Myr ago the distance between the trajectories of these two groupings was minimal.

Characterizing Lossy Dielectric Materials in Shock Physics by Millimeter-Wave Interferometry Using One-Dimensional Convolutional Neural Networks and Nonlinear Optimization

When a dielectric material undergoes mechanical impact, it generates a shock wave, causing changes in its refractive index. Recent demonstrations have proven that the modified refractive index can be determined remotely using a millimeter-wave interferometer. However, these demonstrations are based on the resolution of an inverse electromagnetic problem, which assumes that the losses in the material are negligible. This restrictive assumption is overcome in this article, in which a new approach is proposed to solve the inverse electromagnetic problem in lossy and shocked dielectric materials. Our methodology combines a one-dimensional convolutional neural network architecture, namely Inverse problem of Lossless or Lossy Shocked Wavefront Network (ILSW-Net), with a nonlinear optimization technique based on the Nelder–Mead algorithm to estimate losses within dielectric materials under a mechanical impact. Experimental results for both simulated and real signals show that our method can successfully predict the velocities and the refractive index while accurately estimating the shock wavefront.

Numerical Simulation of Shock Wave Propagation Law of Coal Dust Explosion in Complex Pipeline Networks.

In order to study the shock wave propagation law of coal dust explosion in complex structures, the diagonal pipe network model was established by using computational fluid dynamics software, and the numerical simulation of the coal dust explosion in the diagonal branch was carried out. The change law of the shock wave and its attenuation at the corner were analyzed. The results show that the pressure change process of coal dust explosion can be divided into three stages, namely, the stage where the pressure value at the ignition point is relatively high, the stage where the pressure value on both sides of the ignition source is maximum, and the stage where the pressure value at the shock wavefront is relatively maximum; the overpressure of the shock wave decays at the bifurcation point, and the larger the bifurcation angle, the greater the pressure attenuation amplitude. For the same bifurcation angle, the greater the initial pressure at the monitoring point, the greater its attenuation amplitude. Multiple regression fitting on the shock wave pressure values was performed before and after turning, and the function relationship between the shock wave pressure attenuation coefficient and the initial pressure, as well as different turning angles, was obtained. The fitting results can be used as a formula for calculating the ventilation network during the disaster process, which has theoretical significance for evaluating the consequences of coal dust explosion accidents.

Research on the blockage effect in a shock tube

The shock tube is a crucial experimental apparatus in the field of explosive impact. Limited research has been conducted on the visualization of the internal flow field in shock tubes under blockage conditions, and there is a lack of comprehensive understanding of the mechanisms behind blockage effects. This paper presents experimental and numerical investigations into the propagation of shock waves inside a shock tube and their interaction with obstructive elements. Blockage effect experiments were conducted on a shock tube test platform, where pressure data were obtained through a pressure monitoring system. Additionally, schlieren imaging technology was employed to visualize the evolution of shock wave fronts. A two-dimensional numerical model of the shock tube was established using Fluent software, and the computed results exhibited excellent agreement with the experimental data, confirming the accuracy of the numerical model. The mechanism of the blockage effect was revealed through analysis of pressure-time curves and images depicting the evolution of shock wave fronts. Furthermore, adjustments to parameters on the validated numerical model allowed for analysis and comparison of different blockage ratios. The results indicate that the disturbance to the flow field caused by the blockage effect is primarily due to the reflection of shock waves between the front surface of the obstructive element and the shock tube wall. This reflection leads to differences in pressure curves compared to free-field conditions, and such differences increase with higher blockage ratios.

Expansion kinetics of shock waves generated by laser ablation in copper across a transverse magnetic field

Spatio-temporal evolution of laser-induced shock wave (LISW), initiated by 6 ns laser ablation of Cu target in a transverse magnetic field (B) of diverse strength was investigated by optical beam deflection method. It was found that the average LISW velocity in the near field achieve a value 2.8 km/s which dropped to 2.1 km/s in the absence of B. The shock wave pressure near the target surface at the highest incident fluence reaches ∼ 50 MPa for applied B = 300 mT and nearly become halved without B. The decay of the near field shock front pressure estimated from the Sadowski shock-wave approximation pointed to a strong attenuation mechanism in the initial phase of the shock wave propagation. Interestingly, the shock wave front pressure swiftly increases with gradual increase of incident fluence that finally starts to saturate at high fluence regardless of applied B. The development of shock wave was further mapped with laser ablation of Cu target under transverse B. The ablation depth is found to enhance gradually with increase of B (0–300 mT) which could be attributed to the increase in plasma temperature under B.

Unsteady flow of droplet liquid in hydraulic systems of aircrafts and helicopters: models and analytical solutions

The subject of this study is the unsteady flow of liquid in pipelines, which are part of the design of airplanes and helicopters. This name means, first of all, the phenomenon of a sharp increase in pressure in the pipeline, which is known as a hydraulic shock. Although we have already learned to deal with this phenomenon in some parts of the systems, in many structural elements (flexible pipelines), inside which the working pressure reaches several hundred atmospheres, this phenomenon is still quite dangerous. As you know, the best way to deal with an unwanted phenomenon is through theoretical study. To date, there has been a huge amount of work in the direction of hydraulic shock research. This article does not fully cover these studies. It is limited to references to reviews and relevant works. Because the phenomenon of hydraulic shock has a significantly nonlinear character, analytical solutions of systems of equations corresponding to the simplest models were unknown until recently. This work presents, as an overview, already known analytical solutions describing the process of shock wave propagation. Most importantly, new achievements are given, both for the inviscid approximation and for considering internal viscous friction. It is shown that the internal friction within the considered model is negligible almost everywhere, except for the thin shock layer. The asymptotic is proportional to the tangent function and inversely proportional to the square root of the product of the Reynolds number and the dimensionless parameter characterizing the convection effect. Convection of the velocity field significantly affects the distribution of characteristics in hydraulic shock. If the self-similar solutions that were obtained earlier have a power-law character for the velocity distribution in the shock wave, then the simultaneous consideration in the model of convection and friction on the pipeline walls (according to the Weisbach-Darcy model) made it possible to obtain a distribution in the form of an exponential function that decays with increasing distance from the shock wave front. In addition, the work includes an original approach to solving a nonlinear system of differential equations that describes the propagation of a shock wave without considering the friction on the walls. Analytical solutions were obtained in the form of a function of pressure versus the velocity of shock wave propagation. Research methods. This work uses purely theoretical approaches based on the use of well-known fluid flow models, methods of analytical solution of differential equations and their systems, asymptotic methods, derivation of self-model equations, and finding their solutions. Conclusions. Analytical solutions of systems of differential equations were obtained, which describe models of hydraulic shock without considering viscous effects. A comparison of the obtained results with the results of other studies is given.

DETERMINATION OF THE AMOUNT OF LOAD ON PROTECTIVE STRUCTURES OF CIVIL DEFENSE AGAINST THE INFLUENCE OF A SHOCK WAVE ACCORDING TO DBN V.2.2-5:2023

The main types of explosions and their differences between each other are considered, as well as the difference in excess pressure in the air shock wave is shown. The definition of the calculated load on various load-bearing elements depending on the type of building, its location and structural inhomogeneities of building structures according to DBN B.2.2-5:20023 [1] is given. One of the extreme effects, the danger of which must be taken into account when designing protective structures of civil protection, is an external explosion, which can be caused by both man-made and military factors. The explosion generates a whole complex of loads and impacts. The pressure at the front of the air shock wave, as well as the duration of the compression and rarefaction phases, depend on the amount of explosives and the distance to the center of the explosion or cloud boundary. In general, a detonation explosion is characterized by the possibility of significant pressure on the front of the air shock wave and the short duration of the action (tenths and even hundredths of a second). The air shock wave, when interacting with the structure, affects its above-ground and underground parts. This process is complex: at the beginning of the interaction process, the air shock wave affects the external structures of the above-ground part of the building and the surface of the soil in front of the building, causing a compression wave in it, which, in turn, affects the front outer face of the underground part. Further, the air shock wave, passing through the structural inhomogeneities of building structures in the form of window and door openings, affects the internal load-bearing and non-load-bearing structures of the building, as well as the covering of the underground part. At the same time, compression waves in the soil affect the side and rear walls of the underground part.

Shunt capacitor array shock wave position sensor

The position of the shock wave front under an impact load can reveal many material properties. A design of a shunt capacitive shock wave position sensor has been investigated in this work, utilizing a series of electric probes connected to the shunt capacitor array to measure the moment when the shock wave front reaches a specific point on the sample, which can be used to determine the shock wave front's evolution law. Using the Simulink module in MATLAB software, a sensor circuit model was designed and simulated, and the influence of various conditions on the output signal's falling edge relaxation times were investigated. The results indicate that the relaxation times initially increase, then gradually flatten out, after that, when electrical probes are progressively triggered, the relaxation times will increase again. The relaxation time increase with the rising capacitance, self-inductance coefficient, and voltage drop, and it decrease with increasing coupling coefficient. In addition, it initially decreases and then increase as the initial output voltage increases. Comparing the output signals of the traditional series resistance shock wave position sensor under the same general parameters, it reveals that the voltage falling edge relaxation times for each stage of the sensor are in the order of 10−10 s magnitude, much smaller than the 10−8 s magnitude of the latter, indicating higher measurement accuracy for the sensor. In this work, it also provides a reliable theoretical basis for selecting circuit parameters of the shunt capacitor shock position sensor in experiments by simulation.

An immersed boundary method for the thermo–fluid–structure interaction in rarefied gas flows

An immersed boundary method for the thermo–fluid–structure interaction in rarefied gas flows is presented. In this method, the slip model is incorporated with the penalty feedback immersed boundary method to address the velocity and temperature jump conditions at the fluid–structure interface in rarefied gas flows within the slip-flow regime. In addition, the compressible flows governed by the Navier–Stokes equations are solved by using a high-order finite difference method; the elastic solid is solved by using the finite element method; the fluid and solid dynamics are solved independently, and the thermo–fluid–structure interaction is achieved by using a penalty feedback method in a partitioned way. To model the local rarefaction in the supersonic flow, an artificial viscosity is proposed by introducing the local Knudsen number to diffuse the sharp transition at the shock wave front. Several validations are conducted: the Poiseuille flow in a channel, the flow around a two-dimensional airfoil, a moving square cylinder in a channel, the flow around a sphere, and a moving sphere in quiescent flow. The numerical results from the present method show very good agreements with the previous published data obtained by other methods, confirming the good ability of the proposed method in handling the thermo–fluid–structure interaction in both weakly and highly compressible rarefied gas flows. Inspired by the micro/unmanned aerial vehicles in Martian exploration, the proposed method is applied to the aerodynamics of a flapping wing in rarefied gas flows in both two-dimensional and three-dimensional spaces to demonstrate the versatility of the proposed method for modeling flows involving large deformation and fluid–structure interaction.

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.

Study on cone angle of shockwave front in liquid composite protective structure

Study on cone angle of shockwave front in liquid composite protective structure