ABSTRACT
MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock-bubble interaction, and bubble dynamics. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock-bubble, shock-droplet, and shock-water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas-liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock-bubble-vessel-wall and acoustic-bubble-net interactions are used to demonstrate the full capabilities of MFC.
ABSTRACT
We develop a shock- and interface-capturing numerical method that is suitable for the simulation of multicomponent flows governed by the compressible Navier-Stokes equations. The numerical method is high-order accurate in smooth regions of the flow, discretely conserves the mass of each component, as well as the total momentum and energy, and is oscillation-free, i.e. it does not introduce spurious oscillations at the locations of shockwaves and/or material interfaces. The method is of Godunov-type and utilizes a fifth-order, finite-volume, weighted essentially non-oscillatory (WENO) scheme for the spatial reconstruction and a Harten-Lax-van Leer contact (HLLC) approximate Riemann solver to upwind the fluxes. A third-order total variation diminishing (TVD) Runge-Kutta (RK) algorithm is employed to march the solution in time. The derivation is generalized to three dimensions and nonuniform Cartesian grids. A two-point, fourth-order, Gaussian quadrature rule is utilized to build the spatial averages of the reconstructed variables inside the cells, as well as at cell boundaries. The algorithm is therefore fourth-order accurate in space and third-order accurate in time in smooth regions of the flow. We corroborate the properties of our numerical method by considering several challenging one-, two- and three-dimensional test cases, the most complex of which is the asymmetric collapse of an air bubble submerged in a cylindrical water cavity that is embedded in 10% gelatin.
ABSTRACT
Shockwave lithotripsy repeatedly focuses shockwaves on kidney stones to induce their fracture, partially through cavitation erosion. A typical side effect of the procedure is hemorrhage, which is potentially the result of the growth and collapse of bubbles inside blood vessels. To identify the mechanisms by which shock-induced collapse could lead to the onset of injury, we study an idealized problem involving a preexisting bubble in a deformable vessel. We utilize a high-order accurate, shock- and interface-capturing, finite-volume scheme and simulate the three-dimensional shock-induced collapse of an air bubble immersed in a cylindrical water column which is embedded in a gelatin/water mixture. The mixture is a soft tissue simulant, 10% gelatin by weight, and is modeled by the stiffened gas equation of state. The bubble dynamics of this model configuration are characterized by the collapse of the bubble and its subsequent jetting in the direction of the propagation of the shockwave. The vessel wall, which is defined by the material interface between the water and gelatin/water mixture, is invaginated by the collapse and distended by the impact of the jet. The present results show that the highest measured pressures and deformations occur when the volumetric confinement of the bubble is strongest, the bubble is nearest the vessel wall and/or the angle of incidence of the shockwave reduces the distance between the jet tip and the nearest vessel surface. For a particular case considered, the 40 MPa shockwave utilized in this study to collapse the bubble generated a vessel wall pressure of almost 450 MPa and produced both an invagination and distention of nearly 50% of the initial vessel radius on a ðª(10) ns timescale. These results are indicative of the significant potential of shock-induced collapse to contribute to the injury of blood vessels in shockwave lithotripsy.
