Shock Wave Propagation Research Articles

APPLICATION OF CABARET AND WENO SCHEMES FOR SOLVING THE NONLINEAR TRANSPORT EQUATION IN THE MODELING OF SOUND WAVE PROPAGATION IN THE ATMOSPHERE

The most convenient model for describing the phenomenon of shock wave propagation in the atmosphere is the extended Burgers equation. This work investigates the influence of the numerical scheme on the results of solving the equation, which accounts for the nonlinear nature of shock wave propagation in the atmosphere. This equation is a key component of the extended Burgers equation and defines the transformation of the perturbed pressure profile during its propagation. Two numerical schemes were applied for the solution: CABARET and WENO, which are quasi-monotonic finite difference schemes that allow for solutions without significant numerical oscillations. An analysis of the applicability of these schemes for solving the considered problem was conducted.

Improvements of the damage responses and failure modes of an RC vehicle over pass bridge under contact explosion using FEM-SPH coupled computational frameworks

Winds, earthquakes, and blasts are unpredictable and most challenging dynamic loads. However, the blast is relatively new, and its protection over infrastructure is a vital concern. Conversely, bridges play a critical role in the transport system, and their loss may affect the Nation’s growth. However, the literature says that experimental and numerical simulation is an excellent option to predict the blast impacts. Most of the data from the experiments are collected by photos or videography at high temperatures. Therefore, the reliability of those data is a genuine concern, particularly for bridges. Conversely, if correctly explored, numerical simulation is an option to get into the blast field with high accuracy. A mesh-based method like Finite Element Method (FEM) is excellent for numerical simulation up to a specific temperature. Beyond that, it starts a mesh tangling. Element erosion and remeshing processes may overcome this issue if the mesh-free method simulates. Smoothed Particle Hydrodynamics (SPH) is a mesh-free particle method that can naturally handle massive deformation due to adaptive features. Finally, the current study has coupled two discretisation systems to avoid instability and increase run time. Therefore, FEM-SPH coupling has been adopted to explore the direct damage-reduction, complete damage-time response, damaged contours of the shock wave propagations, damaged contours of particle formation, the effective plastic strain (EPS) and different failure modes of the Vehicle Over Pass typed bridge without and with the protective layer of non-Explosive Reactive Armour (nERA). Finally, the improvement of the bridge is eye-catching as the nERA effectively reduces the deformation, damage, particle formation, EPS, and failure modes.

Numerical Analysis of the Possibility of High-Speed Screw Feeding of Charges for Solid-Fuel Shock Wave Generators

Numerical Analysis of the Possibility of High-Speed Screw Feeding of Charges for Solid-Fuel Shock Wave Generators

Isotropization by shock waves generation in anisotropic hydrodynamics

Anisotropic hydrodynamics (aHydro) has proven successful in modeling the evolution of quark-gluon matter created in heavy-ion collisions. The hydrodynamic description of quark-gluon plasma has also been widely used to study sound phenomena, such as shock waves. It has recently been shown that initial fluctuations in energy density and supersonic partons can generate fairly strong shock waves. However, significant anisotropic properties of the system due to the rapid longitudinal expansion of matter have not been taken into account in such studies. Moreover, the process of isotropization and its characteristic time-scales were not considered along with the question of the shock waves formation. Previous studies on shock discontinuous solutions in anisotropic hydrodynamics assumed constant anisotropy, leading to flow refraction towards the anisotropy axis and flow acceleration, characteristics of rarefaction waves, indicating limitations in this approach. This paper investigates discontinuous solutions for normal shock waves without flow refraction, introducing two compression parameters for longitudinal and transverse pressures. The resulting analytical solutions, as well as numerical computations, provide an isotropization mechanism of the system.

Collapse of microbubbles over an elastoplastic wall

The collapse of a vapour bubble over a material surface has been widely studied over the past few decades, but a comprehensive and quantitative analysis of the cavitation dynamics and its effects on solid materials at the mesoscale (nanometre up to micrometre), which would be of particular interest in applications exploiting cavitation power, is still lacking. In this work, we adopt a diffuse interface model to describe the microbubble dynamics, and a dynamic plasticity model for the solid. The former is particularly suited to studying the rich phenomenology characterising bubble collapse at the mesoscale, which comprises transitions to supercritical conditions, emission and propagation of shock waves, generation of liquid microjets and topological transitions, whereas the latter is used to characterise the permanent plastic deformation caused by the bubble collapse, and has been augmented to consider inertial effects, to assess whether or not an interaction between elastic and plastic waves may influence the resulting deformation. Results concerning the collapse of a microbubble at different liquid overpressures and initial standoff ratios are discussed, and the elastoplastic wave propagation in the solid, together with plastic deformation, is studied for different cases, depending on elastic and plastic material parameters.

Numerical Investigation of Network-Based Shock Wave Propagation of Designated Methane Explosion Source in Subsurface Mine Ventilation System Using 1D FDM Code

In coal mining operations, methane explosions constitute a severe safety risk, endangering miners’ lives and causing substantial economic losses, which, in turn, weaken the production efficiency and economic benefits of the mining industry and hinder the sustainable development of the industry. To address this challenge, this article explores the application of decoupling network-based methods in methane explosion simulation, aiming to optimize underground mine ventilation system design through scientific means and enhance safety protection for miners. We used the one-dimensional finite difference method (FDM) software Flowmaster to simulate the propagation process of shock waves from a gas explosion source in complex underground tunnel networks, covering a wide range of scenarios from laboratory-scale parallel network samples to full-scale experimental mine settings. During the simulation, we traced the pressure loss in the propagation of the shock wave in detail, taking into account the effects of pipeline friction, shock losses caused by bends and obstacles, T-joint branching connections, and cross-sectional changes. The results of these two case studies were presented, leading to the following insights: (1) geometric variations within airway networks exert a relatively minor influence on overpressure; (2) the positioning of the vent positively contributes to attenuation effects; (3) rarefaction waves propagate over greater distances than compression waves; and (4) oscillatory phenomena were detected in the conduits connecting to the surface. This research introduces a computationally efficient method for predicting methane explosions in complex underground ventilation networks, offering reasonable engineering accuracy. These research results provide valuable references for the safe design of underground mine ventilation systems, which can help to create a safer and more efficient mining environment and effectively protect the lives of miners.

Lie symmetry method to discuss the optimal system of Lie subalgebras and similarity solutions for shock propagation in a perfectly conducting dusty gas with magnetic field in rotating medium

In this paper, the similarity solutions for the shock wave propagation in a perfectly conducting dusty gas under the effect of axial or azimuthal magnetic field in the case of rotating medium by using Lie symmetry method have been discussed. The optimal system of one-dimensional Lie subalgebras related to the hyperbolic system of nonlinear partial differential equations had been constructed by using general invariant function and general adjoint transformation matrix. The optimal system is the collection of inequivalent optimal classes, which provides crucial information regarding different types of the similarity solutions. On the basis of the optimal system, we have discussed all possible similarity solutions for inequivalent optimal classes individually. Also, we have found that the similarity solutions exist in two cases for power law shock path and in one case with exponential law shock path. The shock wave decay is due to an increment in the physical parameter such as shock Cowling number or non-idealness parameter or adiabatic exponent or rotational parameter. For smaller value of the ratio of density of solid particles to initial species density of the gas i.e. [Formula: see text], the shock strength decreases with an increment in the mass fraction of solid particles in the mixture [Formula: see text], but the shock strength decreases as [Formula: see text] increases for [Formula: see text]. Also, the shock strength increases with an increment in [Formula: see text] at constant [Formula: see text]. It is observed that the shock wave is stronger in the presence of axial magnetic field as compared to azimuthal magnetic field.

Impact of annular nanosecond plasma actuators on drag reduction in transonic flow

During the last few decades, plasma actuators have emerged as promising devices for aerodynamic flow control. This study focuses on the use of nanosecond plasma actuators for such purposes. A thermal phenomenological model is employed to simulate the effects of these actuators. The propagation of shock waves and their interactions for two specific geometries of plasma actuators, linear and annular plasma synthetic jet actuators, are examined here. A comparative analysis of the performance of these two configurations is presented. Furthermore, the geometric characteristics and temperature model are analyzed to provide insights that can be applied to practical problems. The influence of the actuators on a projectile in the transonic flow is also investigated. The results of the present study show that actuators placed in the conical and cylindrical regions of the object do not contribute to drag reduction. Conversely, actuators positioned at the boat-tail and base of the object effectively reduce drag. This drag reduction is primarily attributed to thermal disturbances in the separation area. Additionally, it is observed that the effects of shock waves and their interaction with stationary waves around the projectile are negligible in terms of drag force reduction.

An improved efficient adaptive method for large-scale multi-explosives explosion simulations

Shock wave caused by a sudden release of high-energy, such as explosion and blast, usually affects a significant range of areas. The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource. This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size. In this work, a variable-domain-size adaptive mesh enlargement (vAME) method is developed based on the proposed adaptive mesh enlargement (AME) method for modeling multi-explosives explosion problems. The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions, rather than in all directions as in AME method. A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy. An estimate of speedup ratio is analyzed for further efficiency comparison. Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave, reflected wave, and Mach wave. Additionally, the vAME method is employed to validate the accuracy, as well as to investigate the performance of the fluid field and shock wave propagation, considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass. The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations. The vAME method yields a competitive efficiency, increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method, respectively. It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.

Locally linearized physics-informed neural networks for Riemann problems of hyperbolic conservation laws

In this work, we demonstrate that physics-informed neural networks (PINNs) tend to propagate predicted shock wave information bidirectionally in time, which does not align with the actual evolution direction of solutions to hyperbolic conservation laws. This mismatch results in instability and hinders the reduction of the loss of governing equations, as well as the initial condition loss by the deep neural network. In order to tackle this problem, we simplify the complexity of the problem by constructing equivalent linear transport equations in the region of shock wave generation. The speeds of these linearized waves are governed by the Rankine–Hugoniot relations of conservation laws. This approach is termed the Locally Linearized PINNs method. Specifically, an appropriate shock wave detector is initially designed to identify domains where shock waves occur. Near shock waves, the original nonlinear equations are transformed into their linearized forms, thereby modifying the residual terms of the partial differential equations. Additionally, an equilibrium factor is introduced in fluid compression regions to reduce prediction errors and stabilize the training of deep neural networks. Numerical examples illustrate that Locally Linearized PINNs effectively address the challenge of predicting global solutions in PINNs and significantly improve shock-capturing performance for hyperbolic conservation laws.

FSI modeling and simulation of blood viscosity impacts on cavitation in mechanical heart valves

Heart valve replacements are critical for patients with valve malfunctions, and the tri-leaflet mechanical heart valve (tMHV) is one of the most durable options available. The tMHV is used to replace malfunctioning heart valves, restoring normal blood flow with exceptional durability, often lasting up to 20 years without needing replacement. This durability makes tMHVs particularly suitable for patients under 60. However, mechanical issues like cavitation can undermine the valve's functionality, posing risks to both its longevity and the overall efficacy of cardiovascular treatments. While previous studies have investigated some aspects of cavitation in these valves, the combined effects of blood viscosity and fluid-structure interaction (FSI) on cavitation dynamics remain insufficiently explored.This work models and numerically simulate the influence of blood viscosity on cavitation within tMHVs, using FSI principles. A detailed geometric model of the tMHV was developed, incorporating experimental data and non-Newtonian fluid behaviour to accurately replicate blood flow. Simulations were conducted in ANSYS Fluent R1® using a transient solver to capture the dynamic FSI, with the Carreau-Yasuda model representing blood's shear-dependent viscosity. Cavitation was observed at pressures as low as 4.5 Pa—well below the 15-20 Pa range typically reported—indicating a higher vulnerability than previously recognized. The simulations further showed significant vapour bubble formation, with the maximum vapour volume fraction reaching 0.953. High-speed leakage flows, peaking at 11 m/s during valve closure, were also noted, considerably exceeding velocities observed in earlier studies.These findings demonstrate that cavitation occurs when blood pressure drops below vapour pressure, causing vapour bubbles to form. These bubbles generate shock waves that can damage the valve surfaces and surrounding tissue. Insights from this study will aid in the design of next-generation tMHVs with optimised flow dynamics, potentially reducing cavitation risks and enhancing patient outcomes by minimising valve-associated complications.

Effects of Connecting Structures in Double-Hulled Water-Filled Cylindrical Shells on Shock Wave Propagation and the Structural Response to Underwater Explosion

In conventional double-hulled submarines, the connecting structures that facilitate the linkage between the two hulls are crucial for load transmission. This paper aims to elucidate the effect of these connecting structures on resistance to shock waves generated by underwater explosions. Firstly, a self-developed numerical solver is built for the one-dimensional water-filled elastically connected double-layer plate model. The shock wave propagation characteristics, shock response of structure, water cavitation, and impact loads transmitted through the gap water and the connecting structures are analyzed quantitatively. The results reveal that the majority of the shock impulse is transmitted by the gap water if the equivalent stiffness of the connecting structures is much less than that of the gap water. Then, a three-dimensional model of the double-hulled, water-filled cylindrical shell is constructed in Abaqus/Explicit, utilizing the acoustic-structural coupling methodology. The analysis focuses on the influence of the thickness and density distribution of the connecting structures on the system’s shock response. The results indicate that a densely arranged connecting structure results in a wavy deformation of the outer hull and a notable reduction in both the impact response and strain energy of the inner hull. When the stiffness of the densely arranged connecting structure is comparatively low, the internal energy and plastic energy of the inner hull are decreased by 16.5% and 24.1%, respectively. The findings of this research are useful for assessing shock resistance and for the design of connecting structures within conventional double-hulled submarines.

Enhancing safety in hydrogen storage: Understanding the dynamic process in hydrogen-methane mixtures during the pressurized leakage

Using natural gas pipelines for hydrogen-doped transportation represents a cost-effective solution for large-scale and long-distance hydrogen transport. However, concerns have been raised regarding the potential for leakage under high-pressure conditions, which could result in the spontaneous combustion of the leaking hydrogen. It is, therefore, imperative to understand the gas dynamics involved in methane-hydrogen mixtures during high-pressure leaks to ensure the safety of hydrogen storage. In the present work, the release process and the spontaneous ignition of hydrogen/methane mixtures with a hydrogen blending ratio from 0.1 to 1.0 from pressurized storage containers (20 MPa) have been numerically studied upon validated models. The findings indicate that spontaneous ignition occurs in the mixtures with any hydrogen blending ratios under the releasing pressure of 20 MPa, but the methane component is not involved in the initial combustion if the hydrogen blending ratio is less than 0.5. As the hydrogen blending ratio increases, more combustible components are involved in the initial ignition, and the location of the initial ignition turns closer to the leakage port. Meanwhile, as the hydrogen blending ratio rises from 0.1 to 1.0, the shockwave propagation is accelerated with the moment the shockwave first reflects is advanced from 58 to 60 μs to 22–24 μs, and the maximum temperature in the early leakage will exceed 3000 K. Furthermore, a higher hydrogen blending ratio leads to a more rapid pressure rise in the external space (reaching peak pressure at 191 μs for pure hydrogen) and forming a Mach disk structure with higher Mach numbers (around 7). These findings provide critical insights for understanding the combustion behavior of hydrogen/methane mixtures in high-pressure pipelines and offer valuable guidance for developing safety protocols and infrastructure designs.

Connectivity between the solar photosphere and chromosphere in a vortical structure

Context. High-resolution solar observations have revealed the existence of small-scale vortices, as seen in chromospheric intensity maps and velocity diagnostics. Frequently, these vortices have been observed near magnetic flux concentrations, indicating a link between swirls and the evolution of the small-scale magnetic fields. Vortices have also been studied with magneto-hydrodynamic (MHD) numerical simulations of the solar atmosphere, revealing their complexity, dynamics, and magnetic nature. In particular, it has been proposed that a rotating magnetic field structure driven by a photospheric vortex flow at its footprint produces the chromospheric swirling plasma motion. Aims. We present a complete and comprehensive description of the time evolution of a small-scale magnetic flux concentration interacting with the intergranular vortex flow and affected by processes of intensification and weakening of its magnetic field. In addition, we study the chromospheric dynamics associated with the interaction, including the analysis of a chromospheric swirl and an impulsive chromospheric jet. Methods. We studied observations taken with the CRisp Imaging SpectroPolarimeter (CRISP) instrument and the CHROMospheric Imaging Spectrometer (CHROMIS) at the Swedish Solar Telescope (SST) in April 2019. The data were recorded at quiet-Sun disc centre, consisting of full Stokes maps in the Fe I line at 6173 Å and in the Ca II infrared triplet line at 8542 Å, as well as spectroscopic maps in the lines of Hα 6563 Å and Ca II K 3934 Å. Utilising the multi-wavelength data and performing height-dependent Stokes inversion, based on methods of local correlation tracking and wavelet analysis, we studied several atmospheric properties during the event lifetime. This approach allowed us to interpret the spatial and temporal connectivity between the photosphere and the chromosphere. Results. We identified the convective collapse process as the initial mechanism of magnetic field intensification, generating a re-bound flow moving upwards within the magnetic flux concentration. This disturbance eventually steepens into an acoustic shock wave that dissipates in the lower chromosphere, heating it locally. We observed prolonged magnetic field amplification when the vortex flow disappears during the propagation of the upward velocity disturbance. We propose that this type of magnetic field amplification could be attributed to changes in the local vorticity. Our analysis indicates the rotation of a magnetic structure that extends from the photosphere to the chromosphere, anchored to a photospheric magnetic flux concentration. It appears to be affected by a propagating shock wave and its subsequent dissipation process could be related to the release of the jet.

An analysis of magnetogasdynamic shock wave propagation in a rotational axisymmetric self-gravitating nonideal gas

An analysis of magnetogasdynamic shock wave propagation in a rotational axisymmetric self-gravitating nonideal gas

Shock wave propagation in a real gas with or without gravitational field in the presence of magnetic field and monochromatic radiation via group invariance method

PurposeThis article aims to find the similarity solutions for the one-dimensional motion of spherical symmetric shock wave in non-ideal gas influenced by the azimuthal magnetic field and monochromatic radiation in the presence or absence of gravitational field. This paper also aims to study the effects of physical parameters on the strength of shock wave, and on the flow variables in the flow-field region behind the shock front.Design/methodology/approachThe Roche model is used to describe the gravitational field effects due to a massive nucleus at the point of symmetry. To derive the similarity solutions, the Lie group symmetry method has been used. Also, the numerical solutions to the present problem are obtained by using Rung–Kutta method of the fourth order with the use of Mathematica software. The effects of variation in the parameter of non-idealness of the gas, the gravitation parameter, the strength of the ambient magnetic field and the adiabatic index of the gas on the shock wave, and on the flow variables is discussed. A comparative study between with and without gravitational field is also, made.FindingsFor different choices of the arbitrary constants that appeared in the solution of infinitesimal generators, we have obtained seven distinct cases of similarity solutions. In the absence of the gravitational field, the similarity solution exists to the power and exponential law shock paths, but in the presence of gravitational field, the similarity solution exists to the power law shock path case only. In the absence of gravitational field, the shock strength is enhanced in the exponential law shock path case in comparison to the power law shock path case. It is found that the shock wave decays with an increase in the value of the adiabatic exponent, the strength of magnetic field, non-idealness of the gas or gravitational parameter.Research limitations/implicationsThe consideration of medium under the influence of gravitational field due to a heavy nucleus at the center and presence of magnetic field decrease the shock strength. This result may be helpful in designing space vehicle and jet engine.Practical implicationsThe result of the present study may be used in the analysis of data from the measurements by space craft in the solar wind and in neighborhood of the Earth’s magnetosphere.Social implicationsThe obtained results may be used for mankind.Originality/valueThe study of spherical shock wave propagation influenced by monochromatic radiation and azimuthal magnetic field in a non-ideal gas with or without gravitational field has yet to be discussed by any authors by using the Lie group symmetry method. In this article, we have discussed all possible cases of similarity solutions using the Lie group symmetry method, which is not studied by anyone as known to us.

Development of crystal optics for X-ray multi-projection imaging for synchrotron and XFEL sources.

X-ray multi-projection imaging (XMPI) is an emerging experimental technique for the acquisition of rotation-free, time-resolved, volumetric information on stochastic processes. The technique is developed for high-brilliance light-source facilities, aiming to address known limitations of state-of-the-art imaging methods in the acquisition of 4D sample information, linked to their need for sample rotation. XMPI relies on a beam-splitting scheme, that illuminates a sample from multiple, angularly spaced viewpoints, and employs fast, indirect, X-ray imaging detectors for the collection of the data. This approach enables studies of previously inaccessible phenomena of industrial and societal relevance such as fractures in solids, propagation of shock waves, laser-based 3D printing, or even fast processes in the biological domain. In this work, we discuss in detail the beam-splitting scheme of XMPI. More specifically, we explore the relevant properties of X-ray splitter optics for their use in XMPI schemes, both at synchrotron insertion devices and XFEL facilities. Furthermore, we describe two distinct XMPI schemes, designed to faciliate large samples and complex sample environments. Finally, we present experimental proof of the feasibility of MHz-rate XMPI at the European XFEL. This detailed overview aims to state the challenges and the potential of XMPI and act as a stepping stone for future development of the technique.

Explosion risks: Variety of deflagration-to-detonation transition scenarios in smooth tubes

In the framework of comprehensive assessment of explosion risks on board of spacecrafts and on the facilities of launch places, the paper is focused on the detailed analysis of particular scenarios of deflagration-to-detonation transition taking place in smooth tubes filled with acetylene-oxygen mixtures of different compositions. By means of precise numerical simulation it is demonstrated that various scenarios of detonation onset can take place depending on the mixture composition and its initial thermodynamic state. It is demonstrated that independent on the particular scenario always the basic mechanism of detonation onset via the formation of strong enough shock wave takes place. In more reactive mixtures the strong shock originates from the self-sustained process of joint pressure build up and reaction intensification exactly at the flame front. In less reactive mixtures the transient flow behavior leads to the shock waves generation and interaction. As a result, a brand new reaction kernel could arise in the area of shock waves interaction. In number of cases, that leads to the coupling between the shock wave and the newborn reaction front and results in the strong shock formation and further detonation onset.

Effects of internal intervention in tubes on shock wave attenuation, formation and propagation of flame during high-pressure hydrogen release

Experiments on the impacts of the combination of fin structure and expansion chamber on shock wave attenuation, hydrogen ignition, and flame characteristics during high-pressure hydrogen release were carried out. The fin is a four-edged platform hollowed out on all sides. The inflated upper and lower tube walls are located 85, 95, and 105 mm away from the nickel burst disk. With the incident shock wave traveling through the fin, the reflected shock and the superposition of multiple reflected shock waves appear, and the reflected shock wave intensity can equal that of the release pressure. The decrease of the incident shock wave intensity behind the fin leads to the decrease of the non-uniformity of the mixed layer temperature, further causing the decrease of the flame propagation velocity. The presence of fins increases the probability of ignition in the fin region and the critical release pressure for ignition is significantly lower than that in the barrier free tube, mainly because the fin inside intensifies the hydrogen/air mixture, resulting in a lower minimum ignition energy required for self-ignition. Furthermore, the decrease of flame propagation velocity leads to the increase of hydrogen consumption inside the tube and the decrease of shock overpressure outside.

Failure analysis of concrete gravity dam using coupled acoustic–structural analysis

The failure mechanism of a gravity dam under blast loading is much more complex than failure under other dynamic loads such as earthquakes. It is therefore important to evaluate the safety performance of these structures under blast load. The structural dynamic behaviour of Koyna dam, India, subjected to underwater explosion (UndEx) was investigated. The Johnson–Holmquist 2 model was used to capture the non-linear response of concrete. The dam–reservoir interaction was modelled using the coupled acoustic–structural technique, which employs a total wave formulation to incorporate the propagation of spherical shock waves. Dynamic responses such as displacement, damage profiles and principal stresses under various detonation depths, weights and reservoir levels were examined. The findings indicated that areas facing the explosion experience the most damage, while local damage occurred at the dam’s heel as the explosion depth increased. Additionally, displacement decreased notably as the TNT depth increased for a consistent explosive amount. The damage distribution in the dam was found to depend on both TNT depth and weight. Damage accumulation in the dam was substantially reduced by lowering the reservoir height (five different levels were assessed). The results provide valuable insights into the Koyna gravity dam’s response to UndEx under varying conditions.