Underwater explosion poses a significant threat to the structural integrity of ocean vehicles and platforms. Accurate prediction of the dynamic loads from an explosion and the resulting structural response is crucial to ensuring safety without overconservative design. When the distance between the explosive charge and the structure is relatively small (i.e., near-field explosion), the dynamics of the gaseous explosion product, i.e., the “bubble”, comes into play, rendering a multiphysics problem that features the interaction of the bubble, the surrounding liquid water, and the solid structure. The problem is highly nonlinear, as it involves shock waves, large deformation, yielding, contact, and possibly fracture. This paper investigates the two-way interaction between the cyclic expansion and collapse of an explosion bubble and the deformation of a thin-walled elastoplastic cylindrical shell in its vicinity. Intuitively, when a shock wave impinges on a thin cylindrical shell, the shell would collapse in the direction of shock propagation. However, some recent laboratory experiments have shown that under certain conditions the shell collapsed in a counter-intuitive mode in which the direction of collapse is perpendicular to that of shock propagation. In other words, the nearest point on the structural surface moved towards the explosion charge, despite being impacted by a compressive shock. This paper focuses on replicating this phenomenon through numerical simulation and elucidating the underlying mechanisms. A recently developed computational framework (“FIVER”) coupling a nonlinear finite element structural dynamics solver and a finite volume compressible fluid dynamics solver is used to complete this study. The solver utilizes an embedded boundary method to track the wetted surface of the structure (i.e. the fluid-structure interface), which is capable of handling large structural deformation and topological changes (e.g., fracture). The solver also adopts the level set method for tracking the bubble surface (i.e. the liquid-gas interface). The fluid-structure and liquid-gas interface conditions are enforced by constructing and solving one-dimensional multi-material Riemann problems, which naturally accommodates the propagation of shock waves across the interfaces. In this paper, mesh refinement study is made to examine the sensitivity of the results to various meshing parameters. The results show that the intermediate level of refinement is appropriate in terms of both the accuracy and the computation costs. Next, the deformation history of both the bubble and the structure are presented and analyzed to provide a detailed view of the counter-intuitive collapse mode mentioned above. We show that timewise, the structural collapse spans multiple cycles of bubble oscillation. Additional details about the time-histories of fluid pressure, structure displacement, and bubble size are presented to elucidate this dynamic bubble-structure interaction and the resulting structural failure.
more »
« less
Shock-induced bubble collapse near solid materials: effect of acoustic impedance
The fluid dynamics of a bubble collapsing near an elastic or viscoelastic material is coupled with the mechanical response of the material. We apply a multiphase fluid–solid coupled computational model to simulate the collapse of an air bubble in water induced by an ultrasound shock wave, near different types of materials including metals (e.g. aluminium), polymers (e.g. polyurea), minerals (e.g. gypsum), glass and foams. We characterize the two-way fluid–material interaction by examining the fluid pressure and velocity fields, the time history of bubble shape and volume and the maximum tensile and shear stresses produced in the material. We show that the ratio of the longitudinal acoustic impedance of the material compared to that of the ambient fluid, $Z/Z_0$ , plays a significant role. When $Z/Z_0<1$ , the material reflects the compressive front of the incident shock into a tensile wave. The reflected tensile wave impinges on the bubble and decelerates its collapse. As a result, the collapse produces a liquid jet, but not necessarily a shock wave. When $Z/Z_0>1$ , the reflected wave is compressive and accelerates the bubble's collapse, leading to the emission of a shock wave whose amplitude increases linearly with $\log (Z/Z_0)$ , and can be much higher than the amplitude of the incident shock. The reflection of this emitted shock wave impinges on the bubble during its rebound. It reduces the speed of the bubble's rebound and the velocity of the liquid jet. Furthermore, we show that, for a set of materials with $Z/Z_0\in [0.04, 10.8]$ , the effect of acoustic impedance on the bubble's collapse time and minimum volume can be captured using phenomenological models constructed based on the solution of Rayleigh–Plesset equation.
more »
« less
- NSF-PAR ID:
- 10202720
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 907
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Recent studies indicate that cavitation may play a vital role in laser lithotripsy. However, the underlying bubble dynamics and associated damage mechanisms are largely unknown. In this study, we use ultra-high-speed shadowgraph imaging, hydrophone measurements, three-dimensional passive cavitation mapping (3D-PCM), and phantom test to investigate the transient dynamics of vapor bubbles induced by a holmium:yttrium aluminum garnet laser and their correlation with solid damage. We vary the standoff distance ( SD) between the fiber tip and solid boundary under parallel fiber alignment and observe several distinctive features in bubble dynamics. First, long pulsed laser irradiation and solid boundary interaction create an elongated “pear-shaped” bubble that collapses asymmetrically and forms multiple jets in sequence. Second, unlike nanosecond laser-induced cavitation bubbles, jet impact on solid boundary generates negligible pressure transients and causes no direct damage. A non-circular toroidal bubble forms, particularly following the primary and secondary bubble collapses at SD = 1.0 and 3.0 mm, respectively. We observe three intensified bubble collapses with strong shock wave emissions: the intensified bubble collapse by shock wave, the ensuing reflected shock wave from the solid boundary, and self-intensified collapse of an inverted “triangle-shaped” or “horseshoe-shaped” bubble. Third, high-speed shadowgraph imaging and 3D-PCM confirm that the shock origins from the distinctive bubble collapse form either two discrete spots or a “smiling-face” shape. The spatial collapse pattern is consistent with the similar BegoStone surface damage, suggesting that the shockwave emissions during the intensified asymmetric collapse of the pear-shaped bubble are decisive for the solid damage.more » « less
-
more » « less
Accurate determination of high strain rate (>103 1/s) constitutive properties of soft materials remains a formidable challenge. Albeit recent advancements among experimental techniques, in particular inertial microcavitation rheometry (IMR), the intrinsic requirement to visualize the bubble cavitation dynamics has limited its application to nominally transparent materials. Here, in an effort to address this challenge and to expand the experimental capability of IMR to optically opaque materials, we investigated whether one could use the acoustic signature of the time interval between the bubble's maximum radius and first collapse time point, characterized as the bubble collapse time, to infer the viscoelastic material properties without being able to image the bubble directly in the tissue. By introducing a modified Rayleigh collapse time for soft materials, which is strongly dependent on the stiffness of the material at hand, we show that, in principle, one can obtain an order of magnitude or better estimate of the viscoelastic material properties of the soft material under investigation. Using a newly developed energy-based theoretical framework, we show that for materials stiffer than 10 kPa the bubble collapse time during a single bubble cavitation event can provide quantitative and meaningful information about the constitutive properties of the material at hand. For very soft materials (i.e., shear modulus less than 10 kPa), our theory shows that unless the collapse time measurement has very high precision and low uncertainties, the material property estimates based on the bubble collapse time only will not be accurate and require visual resolution of the full cavitation kinematics.
-
A new type of wave–mean flow interaction is identified and studied in which a small-amplitude, linear, dispersive modulated wave propagates through an evolving, nonlinear, large-scale fluid state such as an expansion (rarefaction) wave or a dispersive shock wave (undular bore). The Korteweg–de Vries (KdV) equation is considered as a prototypical example of dynamic wavepacket–mean flow interaction. Modulation equations are derived for the coupling between linear wave modulations and a nonlinear mean flow. These equations admit a particular class of solutions that describe the transmission or trapping of a linear wavepacket by an unsteady hydrodynamic state. Two adiabatic invariants of motion are identified that determine the transmission, trapping conditions and show that wavepackets incident upon smooth expansion waves or compressive, rapidly oscillating dispersive shock waves exhibit so-called hydrodynamic reciprocity recently described in Maiden et al. ( Phys. Rev. Lett. , vol. 120, 2018, 144101) in the context of hydrodynamic soliton tunnelling. The modulation theory results are in excellent agreement with direct numerical simulations of full KdV dynamics. The integrability of the KdV equation is not invoked so these results can be extended to other nonlinear dispersive fluid mechanic models.more » « less
-
more » « less
A laser pulse focused near the closed end of a glass capillary partially filled with water creates a vapor bubble and an associated pressure wave. The pressure wave travels through the liquid toward the meniscus where it is reflected, creating a fast, focused microjet. In this study, we selectively coat the hydrophilic glass capillaries with hydrophobic strips along the capillary. The result after filling the capillary is a static meniscus which has a curvature markedly different than an unmodified capillary. This tilting asymmetry in the static meniscus alters the trajectory of the ensuing jets. The hydrophobic strips also influence the advancing contact line and receding contact line as the vapor bubble expands and collapses. We present thirteen different permutations of this system which includes three geometries and four coating schemes. The combination of geometry and coatings influences the jet breakup, the resulting drop size distribution, the trajectory of the jet tip, and the consistency of jet characteristics across trials. The inclusion of hydrophobic strips promotes jetting in line with the channel axis, with the most effective arrangement dependent on channel size.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2020.810
