Pub Date : 2025-09-15DOI: 10.1007/s00193-025-01251-3
D. L. Frost
{"title":"Selected topics from the 26th International Symposium on Military Aspects of Blast and Shock, Wollongong, Australia, December 3–8, 2023","authors":"D. L. Frost","doi":"10.1007/s00193-025-01251-3","DOIUrl":"10.1007/s00193-025-01251-3","url":null,"abstract":"","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 5","pages":"451 - 451"},"PeriodicalIF":1.8,"publicationDate":"2025-09-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145493333","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}
Pub Date : 2025-09-03DOI: 10.1007/s00193-025-01246-0
V. Monnier, L. Vilasi, V. Rodriguez, P. Vidal, R. Zitoun
We identify six scenarios of detonation re-initiation downstream of an obstacle with a single-centered hole in square and round tubes. The tubes have the same cross-sectional area of (16,hbox {cm}^2), and the hole has the same shape as the tube cross section, but different open area ratios. The test mixtures investigated were 2 H2 + O2 + 2 Ar and 2 H2 + O2. We used soot foil recordings, high-speed schlieren, and chemiluminescence imaging to obtain longitudinal and frontal views of the diffraction phenomena. Depending on the initial pressure, one supercritical, four critical, and one subcritical scenarios were observed. The supercritical and critical scenarios were more likely to occur in the square tube than in the round tube, all other parameters being equal. Most of these transient, three-dimensional effects of the cross-sectional shape can be observed even at initial pressures for which there is no effect on steady propagation, e.g., without an obstacle or far from it. This raises the question of what dimensionality refers to in experiments in terms of global and local detonation dynamics.
{"title":"Effects of the channel cross-section shape on detonation re-initiation downstream of a single-hole obstacle","authors":"V. Monnier, L. Vilasi, V. Rodriguez, P. Vidal, R. Zitoun","doi":"10.1007/s00193-025-01246-0","DOIUrl":"10.1007/s00193-025-01246-0","url":null,"abstract":"<div><p>We identify six scenarios of detonation re-initiation downstream of an obstacle with a single-centered hole in square and round tubes. The tubes have the same cross-sectional area of <span>(16,hbox {cm}^2)</span>, and the hole has the same shape as the tube cross section, but different open area ratios. The test mixtures investigated were 2 H<sub>2</sub> + O<sub>2</sub> + 2 Ar and 2 H<sub>2</sub> + O<sub>2</sub>. We used soot foil recordings, high-speed schlieren, and chemiluminescence imaging to obtain longitudinal and frontal views of the diffraction phenomena. Depending on the initial pressure, one supercritical, four critical, and one subcritical scenarios were observed. The supercritical and critical scenarios were more likely to occur in the square tube than in the round tube, all other parameters being equal. Most of these transient, three-dimensional effects of the cross-sectional shape can be observed even at initial pressures for which there is no effect on steady propagation, e.g., without an obstacle or far from it. This raises the question of what dimensionality refers to in experiments in terms of global and local detonation dynamics.\u0000</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 6","pages":"609 - 620"},"PeriodicalIF":1.8,"publicationDate":"2025-09-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00193-025-01246-0.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145772239","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}
Pub Date : 2025-09-02DOI: 10.1007/s00193-025-01241-5
J. W. Denny, D. G. Farrimond, S. Gabriel, S. K. C. Yuen, S. Rigby, S. Al-Hajj, G. S. Langdon
In the field of blast protection engineering, it remains challenging to validate large, complex numerical models and the implications of modelling assumptions relating to how structures are represented (e.g., geometric fidelity) are not well understood. This paper presents experimental work addressing these two issues, in the context of the 2020 Beirut explosion, which remains an important case study for understanding urban blast effects. A series of reduced-scale (1:250) blast tests examined shielding effects caused by the Beirut grain siloes and investigated the influence of the siloes’ geometric fidelity on blast loading. Rigid obstacles were constructed at two geometric fidelities: “rectangular” (i.e., cuboid) and “accurate”, with closer resemblance to the siloes. Pressure gauges were mounted at multiple locations but at fixed blast scaled distances to examine blast–obstacle interaction behaviour. Additionally, Viper::Blast was used to perform computational fluid dynamics analyses of the tests. Experimental findings confirmed significant shielding (reduced pressure and specific impulse) locally behind the siloes ((Z)(<) 3 m/kg1/3), although models indicated that these effects ceased further afield ((Z)(>) 5 m/kg1/3). Overall, blast wave parameters did not exhibit significant differences between the rectangular and accurate representation of the siloes geometry, except for minor differences (10%) in peak overpressures in localised zones. Numerical models confirmed that these discrepancies were caused by differing blast wave scattering, diffraction, and superposition behaviour attributed to the siloes outer geometry. The results suggest that city-scale blast loading analyses can yield reliable results through idealising structures as simplified cuboidal obstacles. These findings will be of direct relevance to blast protection practitioners and researchers concerned with modelling urban blast scenarios.
在防爆工程领域,验证大型复杂的数值模型仍然具有挑战性,并且与结构如何表示(例如几何保真度)相关的建模假设的含义尚未得到很好的理解。本文介绍了在2020年贝鲁特爆炸背景下解决这两个问题的实验工作,这仍然是理解城市爆炸影响的重要案例研究。通过一系列缩小尺寸(1:250)的爆炸试验,研究了贝鲁特粮仓的屏蔽效应,并研究了粮仓几何保真度对爆炸载荷的影响。刚性障碍物以两种几何保真度建造:“矩形”(即长方体)和“精确”,与筒仓更相似。压力表安装在多个位置,但在固定的爆炸比例距离,以检查爆炸-障碍相互作用的行为。此外,使用Viper::Blast对试验进行计算流体动力学分析。实验结果证实了筒仓后面局部的显著屏蔽(降低压力和比冲)((Z)(<) 3 m/kg1/3),尽管模型表明这些影响在更远的地方停止((Z)(>) 5 m/kg1/3)。总体而言,爆炸波参数在筒仓几何形状的矩形和精确表示之间没有显着差异,除了微小的差异(10%) in peak overpressures in localised zones. Numerical models confirmed that these discrepancies were caused by differing blast wave scattering, diffraction, and superposition behaviour attributed to the siloes outer geometry. The results suggest that city-scale blast loading analyses can yield reliable results through idealising structures as simplified cuboidal obstacles. These findings will be of direct relevance to blast protection practitioners and researchers concerned with modelling urban blast scenarios.
{"title":"The influence of obstacle geometric fidelity on blast wave propagation: a reduced-scale case study examining the role of the grain siloes in the 2020 Beirut explosion","authors":"J. W. Denny, D. G. Farrimond, S. Gabriel, S. K. C. Yuen, S. Rigby, S. Al-Hajj, G. S. Langdon","doi":"10.1007/s00193-025-01241-5","DOIUrl":"10.1007/s00193-025-01241-5","url":null,"abstract":"<div><p>In the field of blast protection engineering, it remains challenging to validate large, complex numerical models and the implications of modelling assumptions relating to how structures are represented (e.g., geometric fidelity) are not well understood. This paper presents experimental work addressing these two issues, in the context of the 2020 Beirut explosion, which remains an important case study for understanding urban blast effects. A series of reduced-scale (1:250) blast tests examined shielding effects caused by the Beirut grain siloes and investigated the influence of the siloes’ geometric fidelity on blast loading. Rigid obstacles were constructed at two geometric fidelities: “rectangular” (i.e., cuboid) and “accurate”, with closer resemblance to the siloes. Pressure gauges were mounted at multiple locations but at fixed blast scaled distances to examine blast–obstacle interaction behaviour. Additionally, Viper::Blast was used to perform computational fluid dynamics analyses of the tests. Experimental findings confirmed significant shielding (reduced pressure and specific impulse) locally behind the siloes (<span>(Z)</span> <span>(<)</span> 3 m/kg<sup>1/3</sup>), although models indicated that these effects ceased further afield (<span>(Z)</span> <span>(>)</span> 5 m/kg<sup>1/3</sup>). Overall, blast wave parameters did not exhibit significant differences between the rectangular and accurate representation of the siloes geometry, except for minor differences (10%) in peak overpressures in localised zones. Numerical models confirmed that these discrepancies were caused by differing blast wave scattering, diffraction, and superposition behaviour attributed to the siloes outer geometry. The results suggest that city-scale blast loading analyses can yield reliable results through idealising structures as simplified cuboidal obstacles. These findings will be of direct relevance to blast protection practitioners and researchers concerned with modelling urban blast scenarios.\u0000</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 6","pages":"673 - 699"},"PeriodicalIF":1.8,"publicationDate":"2025-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00193-025-01241-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145772122","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}
Pub Date : 2025-08-12DOI: 10.1007/s00193-025-01242-4
Z. Yang, H. Qiu, Z. Feng, J. Gai, M. Zhao, J. Li, Y. Qi
To develop a detonation combustion chamber intake channel with a short axial length, minimal resistance to the incoming flow, and effective suppression of back-propagating pressure waves, this study investigated a novel intake structure consisting of a U-shaped channel. The study investigated the effects of obstacles with different geometric parameters installed at different positions within the U-channel on the attenuation of pressure waves. To standardize the channel structure, the concept of a gas retention volume ratio is introduced and systematically studied as a key parameter. The findings reveal that the flow area ratio of the channel is the most significant factor influencing the attenuation of pressure waves, while variations in the gas retention volume ratio also affect the wave propagation process. Furthermore, the study reveals that installing obstacle structures in the downstream leg of the U-channel results in bidirectional anisotropy, characterized by different total pressure recovery coefficients for forward and reverse flows.
{"title":"Experimental study on the suppression of pressure waves utilizing a U-shaped channel","authors":"Z. Yang, H. Qiu, Z. Feng, J. Gai, M. Zhao, J. Li, Y. Qi","doi":"10.1007/s00193-025-01242-4","DOIUrl":"10.1007/s00193-025-01242-4","url":null,"abstract":"<div><p>To develop a detonation combustion chamber intake channel with a short axial length, minimal resistance to the incoming flow, and effective suppression of back-propagating pressure waves, this study investigated a novel intake structure consisting of a U-shaped channel. The study investigated the effects of obstacles with different geometric parameters installed at different positions within the U-channel on the attenuation of pressure waves. To standardize the channel structure, the concept of a gas retention volume ratio is introduced and systematically studied as a key parameter. The findings reveal that the flow area ratio of the channel is the most significant factor influencing the attenuation of pressure waves, while variations in the gas retention volume ratio also affect the wave propagation process. Furthermore, the study reveals that installing obstacle structures in the downstream leg of the U-channel results in bidirectional anisotropy, characterized by different total pressure recovery coefficients for forward and reverse flows.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 6","pages":"621 - 634"},"PeriodicalIF":1.8,"publicationDate":"2025-08-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145772240","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}
Pub Date : 2025-08-06DOI: 10.1007/s00193-025-01232-6
S. D. Ryan, N. J. McCormick, R. C. Ripley
It is important to understand the effects of explosions to ensure the safety of civilians and military personnel. Rapid City Planner (RCP) is a comprehensive software tool for predicting the effects of conventional and improvised explosive devices with geographic information system-based outcomes for munition safety, building/structure damage, protection of assets, and human vulnerability in real cities. Modelling cased explosives requires the consideration of casing fragmentation, which is modelled in RCP using three different methods having a different level of detail. The present study focuses on introducing and validating one of the methods, namely the fast primary fragmentation method. The fast fragmentation solver is used to simulate casing breakup of steel-cased cylindrical charges with various TNT- and RDX-based explosives. The predicted velocities were within less than one percent of theoretical Gurney velocities, and the fragment size distributions compared well with experimental data for each explosive. A TNT-filled artillery shell trial was used to validate the fast fragmentation solver in terms of polar distributions of initial fragment speed and number, as well as the spatial spread of fragment throw in a free-field environment. The RCP hydrocode solver is used to assess the equivalent bare charge model used by the fast fragmentation approach. The detailed and fast-modelling approaches produced similar peak overpressures, but underpredicted impulses, at standoff distances between 5 and 10 m (mid-to-far field range). The failure strain and fragment size distribution are shown to have little effect on the blast wave, whereas both would have a significant impact on subsequent fragment effects. Finally, a full-scale scenario was modelled in RCP to show blast and primary fragmentation effects in a coastal urban environment, including a demonstration of urban blast effects from blast pressure and primary fragment trajectories outcomes.
{"title":"Blast and fragmentation modelling in urban environments using Rapid City Planner","authors":"S. D. Ryan, N. J. McCormick, R. C. Ripley","doi":"10.1007/s00193-025-01232-6","DOIUrl":"10.1007/s00193-025-01232-6","url":null,"abstract":"<div><p>It is important to understand the effects of explosions to ensure the safety of civilians and military personnel. Rapid City Planner (RCP) is a comprehensive software tool for predicting the effects of conventional and improvised explosive devices with geographic information system-based outcomes for munition safety, building/structure damage, protection of assets, and human vulnerability in real cities. Modelling cased explosives requires the consideration of casing fragmentation, which is modelled in RCP using three different methods having a different level of detail. The present study focuses on introducing and validating one of the methods, namely the fast primary fragmentation method. The fast fragmentation solver is used to simulate casing breakup of steel-cased cylindrical charges with various TNT- and RDX-based explosives. The predicted velocities were within less than one percent of theoretical Gurney velocities, and the fragment size distributions compared well with experimental data for each explosive. A TNT-filled artillery shell trial was used to validate the fast fragmentation solver in terms of polar distributions of initial fragment speed and number, as well as the spatial spread of fragment throw in a free-field environment. The RCP hydrocode solver is used to assess the equivalent bare charge model used by the fast fragmentation approach. The detailed and fast-modelling approaches produced similar peak overpressures, but underpredicted impulses, at standoff distances between 5 and 10 m (mid-to-far field range). The failure strain and fragment size distribution are shown to have little effect on the blast wave, whereas both would have a significant impact on subsequent fragment effects. Finally, a full-scale scenario was modelled in RCP to show blast and primary fragmentation effects in a coastal urban environment, including a demonstration of urban blast effects from blast pressure and primary fragment trajectories outcomes.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 5","pages":"459 - 474"},"PeriodicalIF":1.8,"publicationDate":"2025-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145493382","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}
Pub Date : 2025-07-30DOI: 10.1007/s00193-025-01233-5
G. A. Shubin, V. N. Smirnov, P. A. Vlasov, V. S. Arutyunov
An experimental and modeling study of the autoignition of (hbox {CH}_{4})–(hbox {O}_{2})–(hbox {Ar}) mixtures with 0.25–4.0% (hbox {CH}_{textrm{4}}) and 2.0% (hbox {O}_{textrm{2}}) was performed at reflected shock wave conditions of (sim )2 atm and (sim )1600–2300 K. The process was monitored by recording the absorption time profiles of (hbox {CH}_{textrm{3}}) and OH radicals at 216.6 nm and 306.772 nm, respectively. The ignition delay time was determined in two ways: as the times it takes to reach the peak (hbox {CH}_{textrm{3}}) concentration or one-half of the maximum OH concentration. Kinetic simulations were carried out using a number of reaction mechanisms, and the predictions were compared to the measurements. An analysis of the sensitivity of the ignition delay time to the rate constants of various elementary stages was conducted, and the main reactions controlling the ignition process were identified. It was demonstrated that uncertainties in experimental conditions, such as the initial temperature and pressure rise rate, produce an effect comparable with that stemming from uncertainties in the rate constants of the key reactions.
0.25-4.0的(hbox {CH}_{4}) - (hbox {O}_{2}) - (hbox {Ar})混合气自燃试验与模型研究% (hbox {CH}_{textrm{4}}) and 2.0% (hbox {O}_{textrm{2}}) was performed at reflected shock wave conditions of (sim )2 atm and (sim )1600–2300 K. The process was monitored by recording the absorption time profiles of (hbox {CH}_{textrm{3}}) and OH radicals at 216.6 nm and 306.772 nm, respectively. The ignition delay time was determined in two ways: as the times it takes to reach the peak (hbox {CH}_{textrm{3}}) concentration or one-half of the maximum OH concentration. Kinetic simulations were carried out using a number of reaction mechanisms, and the predictions were compared to the measurements. An analysis of the sensitivity of the ignition delay time to the rate constants of various elementary stages was conducted, and the main reactions controlling the ignition process were identified. It was demonstrated that uncertainties in experimental conditions, such as the initial temperature and pressure rise rate, produce an effect comparable with that stemming from uncertainties in the rate constants of the key reactions.
{"title":"Methane ignition behind reflected shock waves as monitored via (hbox {CH}_{3}) and OH absorption","authors":"G. A. Shubin, V. N. Smirnov, P. A. Vlasov, V. S. Arutyunov","doi":"10.1007/s00193-025-01233-5","DOIUrl":"10.1007/s00193-025-01233-5","url":null,"abstract":"<div><p>An experimental and modeling study of the autoignition of <span>(hbox {CH}_{4})</span>–<span>(hbox {O}_{2})</span>–<span>(hbox {Ar})</span> mixtures with 0.25–4.0% <span>(hbox {CH}_{textrm{4}})</span> and 2.0% <span>(hbox {O}_{textrm{2}})</span> was performed at reflected shock wave conditions of <span>(sim )</span>2 atm and <span>(sim )</span>1600–2300 K. The process was monitored by recording the absorption time profiles of <span>(hbox {CH}_{textrm{3}})</span> and OH radicals at 216.6 nm and 306.772 nm, respectively. The ignition delay time was determined in two ways: as the times it takes to reach the peak <span>(hbox {CH}_{textrm{3}})</span> concentration or one-half of the maximum OH concentration. Kinetic simulations were carried out using a number of reaction mechanisms, and the predictions were compared to the measurements. An analysis of the sensitivity of the ignition delay time to the rate constants of various elementary stages was conducted, and the main reactions controlling the ignition process were identified. It was demonstrated that uncertainties in experimental conditions, such as the initial temperature and pressure rise rate, produce an effect comparable with that stemming from uncertainties in the rate constants of the key reactions.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 6","pages":"595 - 607"},"PeriodicalIF":1.8,"publicationDate":"2025-07-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145772238","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}
Pub Date : 2025-07-17DOI: 10.1007/s00193-025-01238-0
J. L. Marsh, W. J. Jackson, S. A. Bentil
The effects of temperature on cavitation behavior have been examined in non-biological engineering applications, e.g., water purification, propeller degradation. However, there is a lack of detailed results on how temperature may affect shock-induced cavitation modeling for biological systems. In particular, it is essential to establish if experiments of biological systems (and surrogate biological systems) conducted at room temperature are accurate representations of cavitation at body temperatures. Additionally, many existing works on biological cavitation utilize distilled water. While distilled water is purified, it is not guaranteed to be free from all ions and is still slightly conductive. Since the impact of ion concentration on cavitation behavior has also yet to be quantified, it is preferable to utilize deionized (DI) water for such experiments. As a result, the present study examines the effect of temperature on shock-induced cavitation using a novel shock tube model to visually record the cavitation events in deionized water. Results show a statistically significant relationship between temperature and cavitation bubble number. Although deionized water was used in this study, the results highlight the need to incorporate temperature into future simulations and experiments involving biological fluids in shock wave environments.
{"title":"Influence of temperature on shock-induced cavitation: consequences for blast-induced traumatic brain injury modeling","authors":"J. L. Marsh, W. J. Jackson, S. A. Bentil","doi":"10.1007/s00193-025-01238-0","DOIUrl":"10.1007/s00193-025-01238-0","url":null,"abstract":"<div><p>The effects of temperature on cavitation behavior have been examined in non-biological engineering applications, e.g., water purification, propeller degradation. However, there is a lack of detailed results on how temperature may affect shock-induced cavitation modeling for biological systems. In particular, it is essential to establish if experiments of biological systems (and surrogate biological systems) conducted at room temperature are accurate representations of cavitation at body temperatures. Additionally, many existing works on biological cavitation utilize distilled water. While distilled water is purified, it is not guaranteed to be free from all ions and is still slightly conductive. Since the impact of ion concentration on cavitation behavior has also yet to be quantified, it is preferable to utilize deionized (DI) water for such experiments. As a result, the present study examines the effect of temperature on shock-induced cavitation using a novel shock tube model to visually record the cavitation events in deionized water. Results show a statistically significant relationship between temperature and cavitation bubble number. Although deionized water was used in this study, the results highlight the need to incorporate temperature into future simulations and experiments involving biological fluids in shock wave environments.</p></div>","PeriodicalId":775,"journal":{"name":"Shock Waves","volume":"35 6","pages":"649 - 659"},"PeriodicalIF":1.8,"publicationDate":"2025-07-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145772241","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}
Pub Date : 2025-06-29DOI: 10.1007/s00193-025-01236-2
Y. Z. Liu, Z. Wang, P. Li
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.
{"title":"Design of a novel nozzle for rocket engines with full flow at sea level","authors":"Y. Z. Liu, Z. Wang, 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}
Pub Date : 2025-06-27DOI: 10.1007/s00193-025-01235-3
J. Klein, J. R. Klein, O. Samimi-Abianeh
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.
{"title":"One-dimensional model predictions for the detonation diffraction critical tube diameter","authors":"J. Klein, J. R. Klein, 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}
Pub Date : 2025-06-26DOI: 10.1007/s00193-025-01234-4
Z. L. Hu, L. L. Ma, H. Wu, Q. Fang
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 Z ((Z = R / root 3 of {W}), where R is the distance between the detonation point and the target surface, W 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.
近年来,城市环境中的恐怖袭击和意外爆炸事件频发,造成建筑结构严重损坏甚至倒塌,已成为现代社会关注的主要问题。建筑结构抗震设计与评价的需求日益增加。最大的要求是确定作用在建筑结构上的爆炸荷载,即爆炸波的反射超压。为了更好地平衡复杂冲击波传播及其与建筑物相互作用的计算效率和预测精度,基于ANSYS/AUTODYN,提出了一种综合多阶段法、梯度网格法、映射法、非细化技术等多种现有技术的实用数值模拟方法。首先,将爆炸波的传播过程简化为三个阶段,即从爆炸中心到地面零点的自由空气传播阶段、经过地面反射后的传播阶段和与建筑结构的相互作用阶段。这三个阶段分别采用一维均匀网格和二维/三维渐变网格模型,网格尺寸逐渐增大。然后,采用映射技术,包括网格不细化,将前一阶段的预测结果传递到下一阶段。通过网格敏感性分析,推荐每段对应距离Z ((Z = R / root 3 of {W}), R为爆点与目标表面之间的距离,W为TNT当量装药重量)的网格策略。最后,通过比较爆炸冲击波的超压和脉冲,在单个房屋、十字路口和两个城市街区的四组爆炸试验中验证了所提出的方法和网格尺寸。并对Euler FCT和Euler multi-material两种求解方法进行了比较。建议采用前一种求解方法,因为它具有更高的效率和准确性。本研究可为城市建筑结构的抗震设计和评价提供有益的参考。
{"title":"A practical numerical simulation approach for explosions in large-scale complex urban environments","authors":"Z. L. Hu, L. L. Ma, H. Wu, 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}