skip to main content

Title: 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 more » 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. « less
; ; ;
Award ID(s):
1706003 1751487
Publication Date:
Journal Name:
Journal of Fluid Mechanics
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Shock waves from underwater and air explosions are significant threats to surface and underwater vehicles and structures. Recent studies on the mechanical and thermal properties of various phase-separated elastomers indicate the possibility of applying these materials as a coating to mitigate shock-induced structural failures. To demonstrate this approach and investigate its efficacy, this paper presents a fluid-structure coupled computational model capable of predicting the dynamic response of air-backed bilayer (i.e. elastomer coating – metal substrate) structures submerged in water to hydrostatic and underwater explosion loads. The model couples a three-dimensional multiphase finite volume computational fluid dynamics model with amore »nonlinear finite element computational solid dynamics model using the FIVER (FInite Volume method with Exact multi-material Riemann solvers) method. The kinematic boundary condition at the fluid-structure interface is enforced using an embedded boundary method that is capable of handling large structural deformation and topological changes. The dynamic interface condition is enforced by formulating and solving local, one-dimensional fluid-solid Riemann problems, which is well-suited for transferring shock and impulsive loads. The capability of this computational model is demonstrated through a numerical investigation of hydrostatic and shock-induced collapse of aluminum tubes with polyurea coating on its inner surface. The thickness of the structure is resolved explicitly by the finite element mesh. The nonlinear material behavior of polyurea is accounted for using a hyper-viscoelastic constitutive model featuring a modified Mooney-Rivlin equation and a stress relaxation function in the form of prony series. Three numerical experiments are conducted to simulate and compare the collapse of the structure in different loading conditions, including a constant pressure, a fluid environment initially in hydrostatic equilibrium, and a two-phase fluid flow created by a near-field underwater explosion.

    « less
  2. 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 motionmore »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.« less
  3. Observation of high-speed reactive flows using laser-absorption-based imaging techniques is of interest for its potential to quantitatively reveal both gas-dynamic and thermochemical processes. In the current study, an ultraviolet (UV) laser-absorption imaging method based on nitric oxide (NO) is demonstrated to capture transient flows in a shock tube. A tunable laser was used to generate a continuous-wave UV beam at 226.1019 nm to coincide with a strong NO absorption feature. The UV beam was expanded to a 20-mm diameter and routed through the shock tube to image the flow adjacent to the end wall. Time-resolved imaging was realized using amore »Lambert HiCATT high-speed UV intensifier coupled to a Phantom v2012 high-speed camera. Static absorbance measurements of 1.97% NO/Ar mixtures were first performed to validate the proposed imaging concept, showing good agreement with values predicted by a spectroscopic model. UV laser-absorption images of incident and reflected shock waves captured at 90 kHz temporal resolution are then reported. Translational temperature profiles across the incident and reflected shocks calculated from absorbance images show reasonable agreement with calculated values. After the passage of the reflected shock wave, the flow near the end wall was monitored to probe the development of the end-wall thermal boundary layer. Thermometry measurements across the thermal boundary layer show good agreement with analytical solutions. This study demonstrates the potential of UV laser-absorption imaging in high-speed flow fields, to be applied to more complex applications in the future.« less
  4. Erosion of solid surfaces due to cavitation has been studied for decades. However, it has been a long debate that which mechanism, namely shockwaves, microjets towards the surface, or both, during the cavitation bubble collapse is the primary factor responsible for that erosion. In this project we investigate the small-scale mechanisms of material erosion induced by the collapse of a single cavitation bubble close to a wall. More specifically, our experimental setup includes modification of the initial nucleus size, the maximum bubble radius, the stand-off distance to the wall, the material softness, and the initial flow temperature. We record themore »evolution of the bubble using high speed cameras as well as the local impacts on the materials. With the help of specifically designed cold-wires, we also measure the temperature in the liquid and in the bubble. Two different methods are used to generate the bubble: (i) an acoustic shockwave of variable intensity, (ii) a YAG laser, which may introduce a high temperature at the start. We also combine the two methods in which the laser initially creates a nucleus, then the shockwave triggers the expansion of the bubble. The objectives of the project are included in this paper, while some first results will be presented at the CAV2018 conference.« less
  5. Abstract. Plume-SPH provides the first particle-based simulation ofvolcanic plumes. Smoothed particle hydrodynamics (SPH) has several advantagesover currently used mesh-based methods in modeling of multiphase freeboundary flows like volcanic plumes. This tool will provide more accurateeruption source terms to users of volcanic ash transport anddispersion models (VATDs), greatly improving volcanic ash forecasts. The accuracy ofthese terms is crucial for forecasts from VATDs, and the 3-D SPH modelpresented here will provide better numerical accuracy. As an initial effortto exploit the feasibility and advantages of SPH in volcanic plume modeling,we adopt a relatively simple physics model (3-D dusty-gas dynamic modelassuming well-mixed eruption material,more »dynamic equilibrium and thermodynamicequilibrium between erupted material and air that entrained into the plume,and minimal effect of winds) targeted at capturing the salient features of avolcanic plume. The documented open-source code is easily obtained andextended to incorporate other models of physics of interest to the largecommunity of researchers investigating multiphase free boundary flows ofvolcanic or other origins.

    The Plume-SPH code ( also incorporates several newly developed techniques inSPH needed to address numerical challenges in simulating multiphasecompressible turbulent flow. The code should thus be also of general interestto the much larger community of researchers using and developing SPH-basedtools. In particular, the SPHε turbulence model is used to capturemixing at unresolved scales. Heat exchange due to turbulence is calculated bya Reynolds analogy, and a corrected SPH is used to handle tensile instabilityand deficiency of particle distribution near the boundaries. We alsodeveloped methodology to impose velocity inlet and pressure outlet boundaryconditions, both of which are scarce in traditional implementations of SPH.

    The core solver of our model is parallelized with the message passinginterface (MPI) obtaining good weak and strong scalability using novel techniquesfor data management using space-filling curves (SFCs), object creationtime-based indexing and hash-table-based storage schemes. These techniques areof interest to researchers engaged in developing particles in cell-typemethods. The code is first verified by 1-D shock tube tests, then bycomparing velocity and concentration distribution along the central axis andon the transverse cross with experimental results of JPUE (jet or plume thatis ejected from a nozzle into a uniform environment). Profiles of severalintegrated variables are compared with those calculated by existing 3-D plumemodels for an eruption with the same mass eruption rate (MER) estimated forthe Mt. Pinatubo eruption of 15 June 1991. Our results are consistent withexisting 3-D plume models. Analysis of the plume evolution processdemonstrates that this model is able to reproduce the physics of plumedevelopment.

    « less