Shock-wave Processes Research Articles

Numerical Study of Shock Wave Processes in Materials under the Influence of Ultrashort Laser Pulses Using the Baer–Nunziato Model

Numerical Study of Shock Wave Processes in Materials under the Influence of Ultrashort Laser Pulses Using the Baer–Nunziato Model

Numerical modelling of large elasto-plastic multi-material deformations on Eulerian grids

Abstract The paper deals with an Eulerian model for heterogeneous multiphase (multi-material) elasto-plastic media based on the diffuse interface method in the one-dimensional uniaxial strain approximation. In this approximation, an equilibrium hypoelastic Wilkins bi-material model is derived. This is carried out on the basis of an analogy between the elasto-plastic Wilkins and hydrodynamic Euler models. For the obtained model, a Godunov-type numerical method is developed based on the approximate Riemann solver HLLC. The results of the proposed Eulerian diffuse interface model are compared with reference solutions on moving grids with explicit tracking interfaces. It is shown that the present diffuse interface numerical model describes with good accuracy strong shock-wave processes in heterogeneous multiphase elasto-plastic media.

The vanishing pressure limits of Riemann solutions for the Aw-Rascle hydrodynamic traffic flow model with the logarithmic equation of state

Two kinds of Riemann solutions for the Aw-Rascle hydrodynamic traffic flow model with the logarithmic equation of state are explicitly obtained by using the combination between 1-rarefaction or 1-shock wave along with 2-contact discontinuity. The formation of vacuum state and delta shock wave is identified and analyzed when the perturbation parameter in the pressure term drops to zero, where the intrinsic cavitation and concentration phenomena are surveyed and explored concretely. Additionally, several numerical results displaying the formation process of vacuum state and delta shock wave are also presented by taking three different perturbation parameters for comparison.

Meso–Macro Energy Exchange in Shock-Wave Processes and Dynamic Strength of AB2 Steel

Meso–Macro Energy Exchange in Shock-Wave Processes and Dynamic Strength of AB2 Steel

Numerical study of the interaction of a shock wave with a layer of particles using the Baer–Nunziato model

The work is devoted to the application of the Baer–Nunziato model to study shock wave processes in particle beds. The numerical algorithm is based on an HLLC-like solver and takes into account the processes of establishing equilibrium at the interface, taking into account the effect of particle compaction. The problem statement corresponds to the full-scale experiment of B.C. Fan et al. As part of the work, qualitative and quantitative agreement was obtained with the results of natural experiments and calculations of other authors, and an explanation was proposed for the experimentally observed phenomena.

EXPERIMENTAL AND THEORETICAL STUDY OF SPACE DEBRIS IMPACTS ON SHIELDED TARGETS

A comprehensive study is reported on highand hypervelocity impacts of a steel ball simulating a space debris particle on shielded targets. Experimental studies of high-velocity impact of a steel ball were performed in the velocity range up to 2500 m/s. The obtained data were used to verify the mathematical model and numerical algorithm. Numerical modeling of space debris impact on a shielded target was carried out in the impact velocity range of 1400-7000 m/s using the finite element method implemented in the original EFES software package. The proposed failure algorithm can describe the fragmentation of the material and the formation of new contact boundaries without computational mesh distortion. The specific features of shock wave processes and the destruction of the target and impactor at different impact velocities were investigated.

Experimental study on gas and coal dust explosive overpressure and flame dynamic characteristics in an engineering-level test roadway

The continuous development of coal science and technology has made gas and coal dust explosion disasters an important factor that restricts efficient and intelligent coal mining, which seriously threatens the safe production process of coal mines. To explore the gas and coal dust explosive overpressure and flame propagation characteristics in an actual roadway, the dynamic characteristics of gas and coal dust mixed explosion propagation and evolution laws of explosion flames were investigated using an integrated explosion test system and a high-speed image acquisition system in an engineering-level test roadway with a length of about 700 m and a cross-sectional area of 7.2 m2. Experimental results showed that the peak overpressure measured at each measuring point during the propagation process of explosion shock wave in the roadway did not rise or fall monotonously but fluctuated. The power of explosion shock wave was significantly strengthened by adding coal dust, while the flame propagation speed sharply increased in a certain zone, which generally showed a first increasing and then declining trend. In addition, the flame was blue white after the gas in the roadway was ignited, developed in an irregular shape, and ignited the surrounding combustible gas soon, which further ignited the coal dust under the combined action of pressure wave and flame front. In this case, the flame was deep yellow on the whole. The gas and coal dust explosion flame propagated along the longitudinal section above the roadway, and the flame propagated at an accelerated speed on the transverse section due to the disturbance of obstacles. The study results will provide an important theoretical basis for the R&D of technical active explosion suppression equipment in coal mines and the improvement in their installation technologies.

Numerical Simulation and Measurement of Deformation Wave Parameters by Sensors of Various Types.

The results of applications of various methods for measuring the parameters of high-speed loading using a strain gauge, a fiber Bragg grating located on a metal measuring rod and an interferometer monitoring the movement of the free boundary of the end of the rod are presented. Numerical simulation confirmed the adequacy of the description of the shock-wave process according to experimental data and showed that, with the thickness of the adhesive layer fixing the fiber Bragg grating and the strain gauge on a dimensional rod up to 100 µm, the deformation parameters of the sensors correspond to the parameters of the stress-strain state of the rod. Experimentally, a good correspondence of the results of measuring the magnitude of the relative deformation at a pulse duration of 10-100 µs using sensors of various types is shown, and an estimate of the limit values of the measured values of the deformation wave parameters is given.

The growth and decay of nonlinear wave in polytropic reacting gas with small solid dust particles

The aim of this study is to examine the evolutionary process of shock waves along the characteristic path under the effect of dust particles in a polytropic reacting gas. Testing the nature of the governing quasilinear of this system as a coordinate, we reconstruct the governing equations as Bernoulli’s differential equation and obtain the results. The results in the characteristic plane can display a non-linear behavior in the physical plane. It is shown how the reacting gas parameter influences the evolutionary process of the compressive and expansive waves, respectively, in reacting gas flows. Also, it is shown how the existence of dust in reacting gas influences the growth and decay of the compressive and expansive waves. The transport equation leading to the evolution of shock waves is determined and used to derive the condition for the shock formation. The comparative study of the effect of reacting gas parameters and dust particles on the flow patterning and distorting of shock waves for planar, cylindrically and spherically symmetric flows is also performed.

Numerical and Analytical Investigation of Shock Wave Processes in Elastoplastic Media

Numerical and Analytical Investigation of Shock Wave Processes in Elastoplastic Media

Numerical and Analytical Investigation of Shock Wave Processes in Elastoplastic Media

The Wilkins model for an elastoplastic medium is considered. A theoretical analysis of discontinuous solutions under the assumption of one-dimensional uniaxial strain is performed. In this approximation, the material equations for the deviator stress tensor components are integrated exactly, and only the conservative part of the governing equations remains, which makes it possible to derive a class of exact analytical solutions for the model. To solve the full nonconservative system of equations (without assuming the uniaxial strain), a Godunov-type numerical method is developed, which uses an approximate Riemann solver based on integrating the system of equations along a path in the phase space. A special choice of path is proposed that reduces the two-wave HLL approximation to the solution of a linear equations. Numerical and exact analytical solutions are compared for a number of problems with various regimes of shockwave processes.

Test and damage assessment of shallow buried RC tunnel under explosion

As a major element of the transportation network, tunnels are unavoidably threatened by accidental loads such as vehicle bombs and tank truck explosions. The goal of this research is to explore the dynamic characteristics and damage assessment of tunnel structures under contact blast loads. First, three scaled-down reinforced concrete tunnel models were made, and the explosion test and static loading test were carried out successively to evaluate the axial residual bearing capacity, axial displacement and failure mechanism of the tunnel. Secondly, the finite element model is built by utilizing LS-DYNA, and the reliability of the finite element method is confirmed by comparing the data of the explosion test with the static loading test. At the same time, the calculation method for damage coefficient and the classification criteria for damage grade based on axial residual bearing capacity are presented. Then, based on the finite element method, the propagation process of the explosion shock wave in the tunnel and the damage mechanism of the tunnel are investigated. Finally, seven explosion scenarios are developed, the damage degree of these seven tunnels under the blast load is quantitatively analyzed, and further anti-blast design ideas are put forth. The study in this article may give an intended reference for the damage assessment, anti-explosion design and strengthening work of reinforced concrete tunnels.

Numerical investigation and analysis of indoor gas explosion: A case study of “6·13” major gas explosion accident in Hubei Province, China

Taking the ' 6·13 ′ major gas explosion accident in Shiyan, Hubei Province, China as an example, three problems were studied in this work: (1)The determination of the volume of natural gas involved in the explosion; (2)The propagation process of shock wave inside the building and the damage evolution process of the accident-related building; (3)The overpressure and fragment injury to the person outside the building. Through the numerical simulation in ANSYS/LS-DYNA software, the volume of natural gas involved in the explosion is determined to be 10240 × 1400 × 400 cm (length × width × height) from three perspectives: the damage to the building, the distribution of overpressure inside the building, and the TNT equivalent of the explosion energy. The simulation results are in good line with the accident, which verifies the effectiveness of the scheme and the accuracy of the numerical model. Based on the reasonable filling scheme, the propagation process of shock waves inside the building, the damage evolution process of the building, and the injury ranges of overpressure and fragments outside the building are analyzed. It can be found that the propagation of shock waves in confined space is complex and variable. The explosion shock waves are first reflected and superimposed in the watercourse, resulting in pressure rise. At about 8ms, the shock waves rushed into the first-floor space of the building, and the maximum overpressure was about 0.56 MPa. At about 50 ms, the shock waves rushed into the second-floor space, and the maximum overpressure was about 0.139 MPa. The first and second-floor slabs and infilled walls were almost completely destroyed. The interior walls of the infilled walls are mainly collapsed, and the exterior walls are ejection around the building as the center. The peak displacement and peak velocity of the interior walls of each floor are about 15% of the exterior walls. The fragments which cause fragment injury mainly come from the retaining wall above the watercourse, the maximum velocity is about 89 m/s, and the maximum displacement is 8.9 m. The safety distance of fragment injury is about 8.8 m, while the safety distance of overpressure injury is about 4.6 m. The lethal distance of fragment injury is greater than that of overpressure injury. Compared with the distance between different damage levels of overpressure injury, the difference in fragment injury is small. Therefore, the safety assessment at the engineering level only needs to consider the safety distance of fragment injury. This study can provide suggestions for evaluating the damage of natural gas cloud explosions in confined spaces and is helpful for accident investigation and safety protection.

Damage characteristics of cabin in navigational state subjected to near-field underwater explosion

An arbitrary Lagrange-Euler (ALE) method is used to establish a numerical model for the fluid-structure interaction of near-field underwater explosion while cabins are under sailing conditions. First, the ALE method is used to calculate shock wave and bubble pulsation loads of free field underwater explosion. The results show that the values of shock wave peak, bubble fluctuation period, maximum radius of bubble are similar to empirical formula calculation, which verifies the accuracy of ALE method in calculating explosion loads. On this basis, from the aspects of local damage, longitudinal strength and impact environment of cabin structure, the dynamic response characteristics of cabin structure under the coupling of impact and bubble load at a certain speed and different explosion distance are analyzed and compared with that at static condition. The results show that compared with the static state, the local damage is weakened, the longitudinal strength is increased, the difference of longitudinal strength between navigation and static state is 7%–12%, the shock spectrum value of the compartment corresponding to the explosion source is reduced, and the shock spectrum value of the adjacent compartment is slightly increased; When the cabin is sailing, the increase of the explosion distance within the maximum radius of bubble pulsation increases the action area of bubble pulsation load, and the effect of bubble pulsation load is enhanced, but the coupling enhancement of explosion load is difficult to appear, so that the damage of cabin structure is reduced. The working condition of the cabin suffering from near-field underwater explosion damage under the navigational state analyzed in this paper is closer to the reality, and more truly reflects the process of underwater explosion shock wave and bubble load coupling applied to cabins, which lays a foundation for the fuze warhead coordination design and relevant parameter setting when the torpedo and other underwater weapons give full play to the damage efficiency.

Laser shock dynamic compaction of tungsten carbide reinforced-aluminum matrix (WCp/Al) composites

Laser impact can generate high-strain impact loads. On this basis, a laser impact dynamic compaction (DC) of the tungsten carbide particles (WCp)/Al composite powder with large differences in stiffness and plasticity was performed in this study. Then, the influences of laser energy and WCp content on relative density (actual density/theoretical density), mechanical properties, and microstructure of the pressed billet surface, damage surface, and cross section of the Al matrix composites (WCp/Al) were analyzed. Along with the analysis of the high-strain compaction process, the spread of shock wave process in the compacted billets and the densification behavior of the powder compacts were evaluated via 2D multiparticle finite element method (MPFEM). Results showed that relative density of the compacted billets tended to increase and then decrease with the rise in of WCp content (5%, 10%, 15%, and 20%). In the DC of the composite powder by the laser impact, the pore filling was mainly filled by plastic deformation of the relatively soft Al particles, and the hard WCp did not undergo plastic deformation due to its high strength. A composite powder compact with high relative density was successfully prepared with a WCp content of 10% at 1800 mJ of laser energy. The eventual relative density reached 97.72%, and the overall microhardness increased by approximately 23.25%. The 2D MPFEM simulations could sufficiently predict the dynamic mechanical response in the DC by the laser impact and depict the strain and stress variation patterns at the particle-scale level. Finally, the densification mechanism of the composite powder billets was investigated. The simulation results pertaining to the relationship between the final laser energy and relative density were in significant agreement with the experimental results.

Study on the impulse mechanism of optical films formed by laser plasma shock waves

Abstract In a high-power laser system, when the surface pressure of the optical film caused by laser plasma shock wave is greater than the adhesion per unit area of the film layer, the film will produce mechanical damage, and in serious cases, the whole system may not work. Therefore, studying the formation mechanism of optical film surface pressure and impulse caused by laser plasma shock wave is the key to ensure the normal operation of the high-power laser system. In this article, by studying the relaxation process of laser plasma shock wave on the surface pressure of optical film, and using the time accumulation effect of various pressures on the surface of the optical film, the calculation model of impulse on the optical film’s surface formed by laser plasma shock waves was established, and the variation rules of the impulse I st {I}_{\text{st}} and impulse coefficient j j on the unit area of single-layer Al2O3 and HfO2 optical films with different parameters were obtained. When the incident laser wavelength λ \lambda was 1,064 nm, the energy E E was 0.1 J 0.1\hspace{.1em}\text{J} , the pulse width t p {t}_{\text{p}} was 10 ns 10\hspace{.25em}\text{ns} , the focal length of the focusing lens f f was 350 mm 350\hspace{.25em}\text{mm} , the distance between the film surface and the focal plane of the focusing lens z 0 {z}_{0} was 5 mm 5\hspace{.5em}\text{mm} , and the film radius R R was 5 mm 5\hspace{.5em}\text{mm} , the calculation and simulation results show that the impulse I st {I}_{\text{st}} of the two films was 1 0 − 4 N s 1{0}^{-4}\hspace{.5em}\text{N}\hspace{.5em}\text{s} order of magnitude, the impulse coefficient j of the two films was 1 0 − 5 N s/J 1{0}^{-5}\hspace{.5em}\text{N}\hspace{.5em}\text{s/J} , the Al2O3 film with small atomic number will obtain larger I st {I}_{\text{st}} and j j , I st {I}_{\text{st}} and j j of the two films increase with the increase of E E and f f , and I st {I}_{\text{st}} and j j of the two films decrease with the increase of z 0 {z}_{0} and t p {t}_{\text{p}} . In the total impulse transfer time ( t 0 {t}_{0} ), I st {I}_{\text{st}} and j j both increase with the increase of R R .

Investigation of Formation Process and Intensity of Coal and Gas Outburst Shockwave

The shock wave of a coal and gas outburst is a high-pressure and high-speed impact airflow formed rapidly after the outburst. The propagation destroys the ventilation facilities and causes the destruction of the ventilation system. The theoretical research on the outburst shock wave is of great significance. In order to deeply understand the formation mechanism of the outburst shock wave, this paper draws on the shock wave theory to theoretically analyze the microscopic formation process of the outburst shock wave. The main difference between the formation process of a coal and gas outburst shock wave and the formation process of a general shock wave is that the outburst shock wave has a solid–gas flow zone in the high-pressure zone. The calculation formulas of pressure, density, temperature and other parameters before and after the outburst shock wave are derived. After the outburst shock wave passes through, the pressure, temperature and density of the roadway air will change suddenly. The relationship expression between outburst gas pressure and outburst shock wave intensity is derived, which can reflect the role of pulverized coal in the formation process of a shock wave. In order to facilitate the understanding and calculation, the concept of equivalent sound velocity of coal-gas flow is proposed, and the direct calculation of the impact strength of a coal and gas outburst is attempted. This paper is helpful to improve the understanding of the essence of a coal and gas outburst shock wave. It is also of great significance to outburst disaster relief.

Study on the mechanism of surface pressure of optical films formed by laser plasma shock wave

In a high-power laser system, when the surface pressure of the optical film caused by laser plasma shock wave is greater than the adhesion per unit area of the film layer, the film will produce mechanical damage, and in serious cases, the whole system may not work. Therefore, studying the formation mechanism of optical film surface pressure caused by laser plasma shock wave and calculating the pressure is the key to ensure the normal operation of high power laser system. In this paper, by studying the relaxation process of shock wave on optical film surface pressure, a theoretical calculation model of shock wave on optical film pressure is established, and the variation law of pressure with different parameters is obtained, which reveals the mechanism of forming the optical film surface pressure. The calculation and simulation results show that the maximum pressure is 108 N m−2 during the laser pulse, and the pressure decreases with the increase of laser pulse time after the pulse, and the total action time of laser plasma and shock wave on the film is in the order of microseconds. The pressure increases with the increase of incident laser energy, focal length of focusing lens and incident laser pulse width, which increases with the decrease of the distance between the film surface and the focal plane of the focusing lens. The pressure changes more obviously with the incident laser energy and the distance between the film surface and the focal plane of the focusing lens than with the focal length of the focusing lens and the incident laser pulse width.

Implementation of Flux Limiters in Simulation of External Aerodynamic Problems on Unstructured Meshes

The study is dedicated to the peculiarities of implementing the flux limiter of the flow quantity gradient when solving 3D aerodynamic problems using the system of Navier–Stokes equations on unstructured meshes. The paper describes discretisation of the system of Navier–Stokes equations on a finite-volume method and a mathematical model including Spalart–Allmaras turbulence model and the Advection Upstream Splitting Method (AUSM+) computational scheme for convective fluxes that use a second-order approximation scheme for reconstruction of the solution on a facet. A solution of problems with shock wave structures is considered, where, to prevent oscillations at discontinuous solutions, the order of accuracy is reduced due to the implementation of the limiter function of the gradient. In particular, the Venkatakrishnan limiter was chosen. The study analyses this limiter as it impacts the accuracy of the results and monotonicity of the solution. It is shown that, when the limiter is used in a classical formulation, when the operation threshold is based on the characteristic size of the cell of the mesh, it facilitates suppression of non-physical oscillations in the solution and the upgrade of its monotonicity. However, when computing on unstructured meshes, the Venkatakrishnan limiter in this setup can result in the occurrence of the areas of its accidental activation, and that influences the accuracy of the produced result. The Venkatakrishnan limiter is proposed for unstructured meshes, where the formulation of the operation threshold is proposed based on the gas dynamics parameters of the flow. The proposed option of the function is characterized by the absence of parasite regions of accidental activation and ensures its operation only in the region of high gradients. Monotonicity properties, as compared to the classical formulation, are preserved. Constants of operation thresholds are compared for both options using the example of numerical solution of the problem with shock wave processes on different meshes. Recommendations regarding optimum values of these quantities are provided. Problems with a supersonic flow in a channel with a wedge and transonic flow over NACA0012 airfoil were selected for the examination of the limiter functions applicability. The computation was carried out using unstructured meshes consisting of tetrahedrons, truncated hexahedrons, and polyhedrons. The region of accidental activation of the Venkatakrishnan limiter in a classical formulation, and the absence of such regions in case a modified option of the limiter function, is implemented. The analysis of the flow field around a NACA0012 indicates that the proposed improved implementation of the Venkatakrishnan limiter enables an increase in the accuracy of the solution.

Simulation of an explosion-proof bin made of relaxation material

The article discusses the numerical simulation of a blast wave interaction with a relaxation multiphase material used in design of the destructible protection bin with combined method of suppression of the blast damage effects. Numerical simulation of shock-wave systems and processes occurring during the detonation of a condensed high explosive (HE) in semi-open (open at the top) explosion-proof bins was obtained, including when one which uses multiphase gas-liquid material (at a first approximation – water-bubble mechanical foam) as relaxation substance that absorbs the explosion energy and thus reduces the damaging effect of the blast waves.