Shock Waves最新文献

{"title":"Design of a novel nozzle for rocket engines with full flow at sea level","authors":"Y. Z. Liu,&nbsp;Z. Wang,&nbsp;P. Li","doi":"10.1007/s00193-025-01236-2","DOIUrl":"10.1007/s00193-025-01236-2","url":null,"abstract":"<div><p>To prevent flow separation under overexpanded conditions in traditional large-area-ratio nozzles of rocket engines at sea level, the method of characteristics for wall pressure control is adopted. This method, which is based on thrust-optimized contours, can be implemented to redesign the latter half of a divergent contour to ensure that the wall pressure of the new contour is not less than the critical separation pressure of 0.03 MPa. The newly generated nozzle is named the full-flow nozzle. Then, the design method is verified by simulations, and the performance of full-flow nozzles is evaluated. The results show that the method of wall pressure control can achieve the intended purpose, and the newly generated contour ensures that the nozzle is not only in the full-flow state at sea level but also able to withstand combustion chamber or ambient pressure fluctuations. The combustion chamber pressure is 8.5 MPa, and the specific heat ratio of hot gas is 1.144. Compared with the thrust-optimized contour with an area ratio of 40, in which the flow tends to separate at sea level, the full-flow nozzle can increase the area ratio to 60. Thus, the vacuum specific impulse can be increased by approximately 5.24 s. Compared with the thrust-optimized contour nozzle with an area ratio of 60, the vacuum specific impulse of the full-flow nozzle with an equal area ratio is decreased by 1.57 s.\u0000</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"413 - 421"},"PeriodicalIF":1.8,"publicationDate":"2025-06-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145171473","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":"One-dimensional model predictions for the detonation diffraction critical tube diameter","authors":"J. Klein,&nbsp;J. R. Klein,&nbsp;O. Samimi-Abianeh","doi":"10.1007/s00193-025-01235-3","DOIUrl":"10.1007/s00193-025-01235-3","url":null,"abstract":"<div><p>Detonation diffraction leads to either successful transmission of the detonation or quenching wherein the propagation mechanism is attenuated. The transmission behavior is governed by competing effects of energy release, curvature, and unsteadiness. There is a potentially unique critical diameter that will determine the diffraction outcome for every combustible mixture composition at each set of initial conditions. The critical diffraction diameter has been correlated to several detonation parameters to date; however, these correlations all have limitations. Analytical or quasi-analytical solutions to the diffraction problem, specifically those able to predict the critical diameter, are scarce. The present work develops several critical diameter models by uniting previous work on diffraction phenomena and the critical initiation energy problem. Curvature, decay rate, and energy-based models are established, and their critical diameter predictions are compared against a wide range of experimental critical diameter data. While detonation diffraction is a complex multifaceted phenomenon, a curvature-based one-dimensional model in this work is shown to accurately reproduce empirical critical diameter behavior at relatively low computational cost.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"329 - 347"},"PeriodicalIF":1.8,"publicationDate":"2025-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145170341","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":"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":"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":"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":"Design and flow characterization of the Mach 5 shock tunnel at New Mexico State University","authors":"L. Saenz,&nbsp;J. Frankel,&nbsp;A. Gross,&nbsp;F. Shu","doi":"10.1007/s00193-025-01226-4","DOIUrl":"10.1007/s00193-025-01226-4","url":null,"abstract":"<div><p>This paper provides a comprehensive discussion on the flow characterization and design of the Mach 5 shock tunnel facility at the Hypersonic Research Center of New Mexico State University (NMSU). It reviews the operational principles of low-enthalpy shock tunnels as well as the measurement techniques employed for the facility characterization. The material and thickness of the secondary diaphragm are shown to significantly affect the stability of stagnation properties. Stagnation conditions are determined through an analysis of pressure–time history data measured in the driven tube. Schlieren flow visualizations over a 10<span>(^circ )</span> half-angle straight cone and a sphere are used to estimate the freestream Mach number. Additionally, femtosecond laser electronic excitation tagging (FLEET) velocimetry is conducted to measure instantaneous velocities in the freestream and turbulent boundary layer flows within the test section. The shock tunnel has a total temperature ranging between 610 and 630 K, with a freestream Mach number of 5.1. The steady test time, as indicated by pitot pressure measurements, ranges from 2 to 2.5 ms, while velocimetry and wall-static pressure data suggest that driver gas arrival in the test section occurs approximately 30 ms after flow stabilization. The facility was made available for use in undergraduate courses in Fall 2022.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 4","pages":"437 - 450"},"PeriodicalIF":1.8,"publicationDate":"2025-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145162707","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}