Shock Waves最新文献

{"title":"A practical numerical simulation approach for explosions in large-scale complex urban environments","authors":"Z. L. Hu,&nbsp;L. L. Ma,&nbsp;H. Wu,&nbsp;Q. Fang","doi":"10.1007/s00193-025-01234-4","DOIUrl":"10.1007/s00193-025-01234-4","url":null,"abstract":"<div><p>In recent years, terrorist attacks and accidental explosions in urban environments have occurred frequently, causing severe damage, even collapse, of building structures, and have become a major concern of modern society. The need to design and evaluate the blast resistance of building structures is rising markedly. The utmost requirement is the determination of blast loads acting on building structures, i.e., the reflected overpressure of blast waves. To better keep the balance between computational efficiency and prediction accuracy of complex blast wave propagation and its interactions with buildings, a practical numerical simulation approach integrating multiple existing techniques including the multi-stage method, graded mesh, mapping, and un-refinement technique is proposed based on ANSYS/AUTODYN. Firstly, the propagation of blast waves is simplified into three stages, i.e., propagation in the free air from the explosion center to ground zero, propagation after the ground reflection, and interaction with building structures. These three stages are modeled by 1D uniform meshes and 2D/3D graded meshes with increasing mesh sizes. Then, the mapping technique, including mesh un-refinement, is adopted to transfer the predicted results at the previous stage into the next stage. The corresponding meshing strategy against the scaled distances <i>Z</i> (<span>(Z = R / root 3 of {W})</span>, where <i>R</i> is the distance between the detonation point and the target surface, <i>W</i> is the equivalent charge weight of TNT) for each stage is recommended through mesh sensitivity analyses. Finally, the proposed approach and mesh sizes are validated against four series of explosive tests for a single house, an intersection, and two city blocks by comparing with both the overpressures and impulses of blast waves. Additionally, two solvers, i.e., Euler FCT and Euler multi-material, are compared. The former solver is recommended due to its greater efficiency and accuracy. The present work could provide a helpful reference for the blast-resistant design and evaluation of urban building structures.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"301 - 327"},"PeriodicalIF":1.8,"publicationDate":"2025-06-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145169544","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Experimental and numerical studies of a near-field blast loading model for cylindrical charges","authors":"S.-L. Liang,&nbsp;J. Yu,&nbsp;L. Chen","doi":"10.1007/s00193-025-01229-1","DOIUrl":"10.1007/s00193-025-01229-1","url":null,"abstract":"<div><p>In the current practice of blast-resistant design, blast loads are determined by the Kingery and Bulmash charts in accordance with a database of free-air blasts of spherical charges and surface bursts of hemispherical charges initiated at the center. However, most charges are closer to cylinders in geometry. In addition, charge shapes and initiation configurations significantly affect blast loads under a near-field blast scenario. Therefore, it is imperative to develop a near-field blast loading model for cylindrical charges that can account for the effects of both charge shape and initiation configuration in blast-resistant design. Compared with incident blast loads, reflected blast loads are more relevant because the latter can be directly used for blast-resistant design. Accordingly, in this paper, experimental and numerical studies were performed to develop a near-field blast loading model for cylindrical charges in terms of the peak reflected overpressure and the maximum reflected impulse. Two series of tests were conducted with either one-end-initiated or both-end-initiated cylindrical charges to obtain reflected blast loads with different scaled distances. It was found that the spatial distribution of blast loads along the axial direction of the charges was extremely non-uniform. Then, high-efficiency numerical models were built using 2D to 3D mapping techniques. After being validated against experimental results, numerical models were employed to simulate the blast loads generated by cylindrical charges with different length-to-diameter ratios and initiation configurations (one-end, center, and both-end initiations) with scaled distances ranging from 0.2 to 1.0 m/kg<span>(^{mathrm {1/3}})</span>. To develop the blast loading model, the peak reflected overpressure and the maximum reflected impulse at the center of a rigid reflection surface were firstly determined by curve fitting as the benchmark blast loads, which were expressed as functions of scaled distance and length-to-diameter ratio, and then the benchmark blast loads were used to normalize the blast loads at different locations. Accordingly, the spatial distribution of blast loads can be described with the benchmark blast loads and a spatial load distribution function, in which the latter is determined by surface fitting of extensive numerical results. The results indicate that the blast loading model developed is able to predict the blast load with considerable accuracy.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"277 - 300"},"PeriodicalIF":1.8,"publicationDate":"2025-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145169802","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Soil-filled perimeter walls under blast","authors":"M. L. Ruiz-Ripoll,&nbsp;C. Roller,&nbsp;H. Dirlewanger,&nbsp;A. Stolz","doi":"10.1007/s00193-025-01223-7","DOIUrl":"10.1007/s00193-025-01223-7","url":null,"abstract":"<div><p>Soil-filled gabion systems can be used in many civil applications such as retaining walls against flooding and erosion or shoreline protection. In addition, the gabion systems provide good resistance in high dynamic loading scenarios such as blast events. These systems allow for a modular setup of easy-to-use perimeter walls with variable height and cross section, application as a gravity wall, and use of local filling material. The latter is the subject of the present paper. Depending on aggregate size and morphology, size distribution, and humidity, soil materials exhibit different material properties such as compaction parameters, cohesion, and the angle of friction among others. Each of these parameters directly affects the structure’s response under highly dynamic conditions. To understand the influence of varying soil parameters at varying loading conditions and thus to predict the structure’s behavior precisely, the authors investigated soil-filled perimeter walls experimentally and using hydrocode simulations. Since the soil’s properties primarily influence the wall’s behavior—at the resistance side—an extensive laboratory test campaign was required to characterize different soils. The experimental data serve for the derivation of dynamic material models and are complemented by numerical simulations. Furthermore, this paper describes the execution of near-field detonation and shock tube tests of soil-filled perimeter walls to analyze their load-bearing behavior under blast load. The experiments are evaluated with regard to the failure mechanism as well as the blast mitigation. Additionally, the blast mitigation effect is numerically investigated and the results are compared to the experiments.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 5","pages":"505 - 527"},"PeriodicalIF":1.8,"publicationDate":"2025-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00193-025-01223-7.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145493378","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A parametric analysis of streamwise vortices on a compression ramp at Mach 4","authors":"K. Dixit,&nbsp;R. R. Kumar,&nbsp;N. R. Vadlamani,&nbsp;N. Tsuboi","doi":"10.1007/s00193-025-01239-z","DOIUrl":"10.1007/s00193-025-01239-z","url":null,"abstract":"<div><p>Large-eddy simulations are conducted to investigate supersonic flow over a compression ramp at a free stream Mach number of 4.0 and a unit Reynolds number of <span>(4.56 times 10^{6})</span> per meter. Two ramp angles of <span>(15^circ )</span> and <span>(18^circ )</span> are considered along with three different ramp positions (P1, P2, and P3) from the plate leading edge, with the plate length increasing progressively from P1 to P3. Simulations reveal that with an increase in the plate length and ramp angle, the separation point shifts downstream, accompanied by an extended separation length. Furthermore, with an increase in the ramp angle and plate length, a higher Görtler number is observed upstream of the reattachment indicating a greater likelihood of Görtler instability. In particular, no streamwise vortices were observed for the 15P1 and 15P2 cases, while for the 18P3 case, increased instability resulted in the breakdown of streamwise vortices, driving the transition to turbulence. The wavelength of streamwise streaks decreased by approximately <span>(15%)</span> as the plate length increased by <span>(approx 80%)</span> from 18P1 to 18P3. Unsteady analysis revealed the role of spanwise secondary instabilities over these vortices, that trigger turbulent spots that propagate at a speed of <span>(approx 0.6 U_{infty })</span>. The peak value of the Stanton number is found to be <span>(approx )</span> 15–27% higher than the time and span-averaged value for the 15P3 and 18P3 cases, highlighting a strong effect of downwash due to streamwise vortices on the wall heating rate distribution. The unsteady data also reveal a negative correlation between the flow reattachment location and the Stanton number close to the reattachment point. An earlier reattachment is shown to increase the Stanton number and vice versa resulting in a <span>(approx 40%)</span> variation compared to the time-averaged value. The results from this study underscore the critical influence of plate length on the formation of streamwise vortices, with significant implications for wall heating rate distribution and flow transition dynamics.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"423 - 436"},"PeriodicalIF":1.8,"publicationDate":"2025-06-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00193-025-01239-z.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145165867","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Role of grain architecture in shock behavior and spalling behavior of Al metal-(text {Al}_{90}text {Sm}_{10}) metallic glass nanolaminates","authors":"S. Mishra,&nbsp;K. Vijay Reddy,&nbsp;S. Pal","doi":"10.1007/s00193-025-01231-7","DOIUrl":"10.1007/s00193-025-01231-7","url":null,"abstract":"<div><p>Molecular-dynamics-based simulations have been carried out for crystalline Al-<span>(text {Al}_{90}text {Sm}_{10})</span> metallic glass (MG) nanolaminates with different grain structures corresponding to varying values of shock intensities to analyze the structural evolution during shock-wave loading and spallation behavior of the nanolaminates. A transition from elastic–plastic behavior occurs in nanocrystalline NC-MG nanolaminates with increasing values of shock intensities when the shock traverses from the crystalline end to the MG end. On the other hand, an overdriven elastic front is observed for all values of shock intensities in columnar-grained CG-MG nanolaminates. When the shock-wave direction is reversed, a plastic wave dominates the shock profiles irrespective of the grain structures and shock intensity values. Adaptive common neighbor analysis (a-CNA) and dislocation analysis reveal that grain boundary-mediated plasticity is dominant in NC-MG nanolaminate specimens, while dislocation-mediated plasticity predominately governs the shock deformation behavior in CG-MG nanolaminates. The reflection of the rarefaction wave generated at the crystalline–amorphous interface aids in stacking fault generation in NC-MG nanolaminates but does not cause any structural changes in CG-MG nanolaminates. The spallation behavior of the nanolaminate specimens is significantly influenced by the grain structures and the presence of the free surfaces. The population of perfect icosahedral clusters <span>(langle 0,0,12,0rangle )</span> decreases during the passage of shock as determined using Voronoi cluster analysis.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"361 - 380"},"PeriodicalIF":1.8,"publicationDate":"2025-06-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145165824","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Effect of diffuse initial conditions on the Richtmyer–Meshkov instability","authors":"S. Pellone,&nbsp;T. Desjardins,&nbsp;K. Prestridge,&nbsp;J. Charonko","doi":"10.1007/s00193-025-01237-1","DOIUrl":"10.1007/s00193-025-01237-1","url":null,"abstract":"<div><p>In this study, we investigate the role of a diffuse interface on the Richtmyer–Meshkov (RM) instability by performing two-dimensional simulations of a single-mode perturbation (wavelength <span>(lambda )</span>) imposed on a diffuse interface (thickness <span>(delta )</span>) between air and <span>(hbox {SF}_6)</span>. By varying the ratio <span>(0.1le delta /lambda le 0.5)</span>, we examine the effect of the interfacial diffusion thickness on the baroclinic vorticity, perturbation growth, and fluid entrainment. The initial circulation is conserved with respect to <span>(delta )</span>, causing a reduction of the initial vorticity magnitude, thus resulting in a reduction of perturbation growth. In the linear regime, the diffusion layer delays perturbation growth, but in the nonlinear regime, the growth becomes insensitive to the initial diffusion thickness, as shown by our power-law scaling accounting for the redistribution of vorticity along the interface. The initial diffusion thickness increases the overall volume of the roll-up, but decreases its surface area. Introducing a new metric (the inter-fluid distance, <i>d</i>) reveals that initially thicker interfaces increase material separation and reduce strain rates within the roll-up structures, resulting in longer diffusion length scales. These structures undergo a gradual thinning over time, causing the inter-fluid distance to decrease to scales comparable to the strain-dominated diffusion length. Therefore, while the strain rate dominates the vortex-core evolution early on, the effect of diffusion may become important at later times, with this transition delayed for thicker initial interfaces.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"395 - 411"},"PeriodicalIF":1.8,"publicationDate":"2025-06-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145166434","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"In memoriam Dr. Robert (Rob) C. Ripley (1975–2024)","authors":"F. Zhang,&nbsp;L. Donahue,&nbsp;D. L. Frost","doi":"10.1007/s00193-025-01228-2","DOIUrl":"10.1007/s00193-025-01228-2","url":null,"abstract":"","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 5","pages":"453 - 457"},"PeriodicalIF":1.8,"publicationDate":"2025-06-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145493383","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical study of the influence of interfaces between explosives and steel plates","authors":"J. Liu,&nbsp;Z. Y. Sun,&nbsp;J. W. Yin,&nbsp;F. G. Zhang","doi":"10.1007/s00193-025-01230-8","DOIUrl":"10.1007/s00193-025-01230-8","url":null,"abstract":"<div><p>The interface between explosives and steel plates can vary in the clearance of a gap or the presence of a cushion, and the dimension of the interface region can also differ. These variations in the types of interface may affect the dynamic loading and energy absorption of steel plates driven by detonation. To investigate this issue, we conducted a numerical study on the influence of different interface types and thicknesses. Initially, we designed a simulation model of a detonation driving a steel plate, with one half featuring clearance between the explosive and the steel, and the other half filled with a cushion. We then carried out simulations to analyze the influence of varying clearance and cushion thickness on the dynamic loading and energy absorption of the steel plate. The results indicate that a small-thickness clearance between explosive and steel can increase the kinetic energy of the steel plate, and there may be an optimal clearance thickness to maximize the energy absorption of the steel plate. As the clearance thickness is increased, the first loading pressure in the steel decreases, and the spallation and recompression processes in the steel gradually transform into an approximately uniform loading process without fracture. On the other hand, filling the clearance with a cushion has a negative effect on the energy absorption of the steel plate, and the kinetic energy of the steel plate decreases nearly linearly with an increase of the cushion thickness. As the cushion thickness is increased, the first loading pressure in the steel decreases less, and the dynamic behaviors of spallation and recompression may occur. Lastly, we briefly discuss interfaces with uneven thickness, which should be strictly controlled to prevent the occurrence of unexpected phenomena.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"381 - 393"},"PeriodicalIF":1.8,"publicationDate":"2025-06-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145165428","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Acoustic measurements of laser-induced microshocks: time of arrival to yield estimations","authors":"E. R. Wainwright,&nbsp;S. W. Dean","doi":"10.1007/s00193-025-01227-3","DOIUrl":"10.1007/s00193-025-01227-3","url":null,"abstract":"<div><p>High-fidelity microphones can be used to characterize laser-generated microblast waves (“microshocks”) in tabletop experiments. This study probes both spherical and hemispherical microshocks, analogous to height-of-burst or surface-burst geometries, at distances of 1–15 cm and laser energies in the range of ~  300–630 mJ under face-on (<span>(0^{circ })</span>) or side-on (<span>(90^{circ })</span>) microphone orientations. We take a Kingery–Bulmash-style analysis approach and calculate the characteristic fitting parameters for time of arrival of the microshock. Blast waves from these laser energies cover scaled distances of ~  2–50 m/<span>(hbox {kg}^{mathrm {1/3}})</span>, roughly equivalent to the detonation of a few grams of TNT probed from several meters away. We compare the experimental results to BlastX simulations and tabulated data from a variety of sources. Under this experimental configuration, a 302-mJ laser pulse is equivalent to a TNT charge in the mass range 1–18 <span>(upmu )</span>g and the 628-mJ pulse is within the range 10–45 <span>(upmu )</span>g. This corresponds to a laser energy to shock coupling ratio when compared to 100% TNT equivalence of 1–24% and 7–29%, respectively. This work informs microblast scaling expectations for experiments using laser-induced shock waves as a microscale energetic characterization technique and provides connections between laboratory and free-field detonation testing.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"267 - 276"},"PeriodicalIF":1.8,"publicationDate":"2025-06-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145161732","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Viscoelastic materials evaluated for blast-resistant designs","authors":"M. Sutter,&nbsp;C. Thomas,&nbsp;A. D. Douglas,&nbsp;F. Dogan,&nbsp;C. E. Johnson","doi":"10.1007/s00193-025-01225-5","DOIUrl":"10.1007/s00193-025-01225-5","url":null,"abstract":"<div><p>Viscoelastic materials have extensive military applications due to their energy absorption capabilities, with the potential to reduce blast energy imposed on buildings, vehicles, and personnel. Based on current literature, limited information is available regarding the mitigation of blast energy related to these uses. The impact of thickness, nanoparticle addition, and layering variation was assessed in this study using commercially available viscoelastic materials in open-air blasts of Composition C4 to determine shock energy mitigation capabilities. Time-pressure waveforms were recorded to identify optimal changes in shock wave characteristics: reduced peak pressure, positive phase duration, and impulse, with increased rise times. Results were analyzed through trend and linear regression analysis to evaluate factors possibly influencing the behavior of the materials. Polyurethane-based materials reduced peak pressures by extending the positive phase duration, whereas silicone rubber maintained a similar duration with reduced peak pressures, suggesting differing energy dissipation mechanisms. Polyurethane was more effective due to its pressure reduction regardless of thickness, enabling thinner layers to be used to achieve similar results. Overall, thinner layers were more efficient, as diminished returns were evident by asymptotic points once reaching a 7-mm thickness. Incorporating graphene nanoplatelets increased energy transfer with peak pressure increases up to 16% in the polyurethane-based samples and impulse increases of 7.5% in the silicone rubber-based samples, making the baseline samples more effective. Layers alternating in material type reduced peak pressures up to 16% relative to baseline samples, with the most reduction occurring in the thicker layers. The alternate layering patterns proved pivotal in the results, those starting with silicone rubber being correlated to increases of 21% in positive phase duration and 6.5% in decay time.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 5","pages":"551 - 562"},"PeriodicalIF":1.8,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00193-025-01225-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145493319","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}