Non-spherical collapse of an air bubble subjected to a lithotripter pulse

In order to better understand the contribution of bubble collapse to stone comminution in shockwave lithotripsy, the shockinduced and Rayleigh collapse of a spherical air bubble is investigated using numerical simulations, and the free-field collapse of a cavitation bubble is studied experimentally. In shock-induced collapse near a wall, it is found that the presence of the bubble greatly amplifies the pressure recorded at the stone surface; the functional dependence of the wall pressure on the initial standoff distance and the amplitude are presented. In Rayleigh collapse near a solid surface, the proximity of the wall retards the flow and leads to a more prominent jet. Experiments show that re-entrant jets form in the collapse of cavitation bubbles excited by lithotripter shockwaves in a fashion comparable to previous studies of collapse near a solid surface

MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver

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

Development and application of unified algorithms for problems in computational science

A framework is presented for developing computationally unified numerical algorithms for solving nonlinear equations that arise in modeling various problems in mathematical physics. The concept of computational unification is an attempt to encompass efficient solution procedures for computing various nonlinear phenomena that may occur in a given problem. For example, in Computational Fluid Dynamics (CFD), a unified algorithm will be one that allows for solutions to subsonic (elliptic), transonic (mixed elliptic-hyperbolic), and supersonic (hyperbolic) flows for both steady and unsteady problems. The objectives are: development of superior unified algorithms emphasizing accuracy and efficiency aspects; development of codes based on selected algorithms leading to validation; application of mature codes to realistic problems; and extension/application of CFD-based algorithms to problems in other areas of mathematical physics. The ultimate objective is to achieve integration of multidisciplinary technologies to enhance synergism in the design process through computational simulation. Specific unified algorithms for a hierarchy of gas dynamics equations and their applications to two other areas: electromagnetic scattering, and laser-materials interaction accounting for melting

Numerical simulations of non-spherical bubble collapse

A high-order accurate shock- and interface-capturing scheme is used to simulate the collapse of a gas bubble in water. In order to better understand the damage caused by collapsing bubbles, the dynamics of the shock-induced and Rayleigh collapse of a bubble near a planar rigid surface and in a free field are analysed. Collapse times, bubble displacements, interfacial velocities and surface pressures are quantified as a function of the pressure ratio driving the collapse and of the initial bubble stand-off distance from the wall; these quantities are compared to the available theory and experiments and show good agreement with the data for both the bubble dynamics and the propagation of the shock emitted upon the collapse. Non-spherical collapse involves the formation of a re-entrant jet directed towards the wall or in the direction of propagation of the incoming shock. In shock-induced collapse, very high jet velocities can be achieved, and the finite time for shock propagation through the bubble may be non-negligible compared to the collapse time for the pressure ratios of interest. Several types of shock waves are generated during the collapse, including precursor and water-hammer shocks that arise from the re-entrant jet formation and its impact upon the distal side of the bubble, respectively. The water-hammer shock can generate very high pressures on the wall, far exceeding those from the incident shock. The potential damage to the neighbouring surface is quantified by measuring the wall pressure. The range of stand-off distances and the surface area for which amplification of the incident shock due to bubble collapse occurs is determined

Cavitation Induction by Projectile Impacting on a Water Jet

The present paper focuses on the simulation of the high-velocity impact of a projectile impacting on a water-jet, causing the onset, development and collapse of cavitation. The simulation of the fluid motion is carried out using an explicit, compressible, density-based solver developed by the authors using the OpenFOAM library. It employs a barotropic two-phase flow model that simulates the phase-change due to cavitation and considers the co-existence of non-condensable and immiscible air. The projectile is considered to be rigid while its motion through the computational domain is modelled through a direct-forcing Immersed Boundary Method. Model validation is performed against the experiments of Field et al. [Field, J., Camus, J. J., Tinguely, M., Obreschkow, D., Farhat, M., 2012. Cavitation in impacted drops and jets and the effect on erosion damage thresholds. Wear 290–291, 154–160. doi:10.1016/j.wear.2012.03.006. URL http://www.sciencedirect.com/science/article/pii/S0043164812000968 ], who visualised cavity formation and shock propagation in liquid impacts at high velocities. Simulations unveil the shock structures and capture the high-speed jetting forming at the impact location, in addition to the subsequent cavitation induction and vapour formation due to refraction waves. Moreover, model predictions provide quantitative information and a better insight on the flow physics that has not been identified from the reported experimental data, such as shock-wave propagation, vapour formation quantity and induced pressures. Furthermore, evidence of the Richtmyer-Meshkov instability developing on the liquid-air interface are predicted when sufficient dense grid resolution is utilised

Comparison of multiphase SPH and LBM approaches for the simulation of intermittent flows

Smoothed Particle Hydrodynamics (SPH) and Lattice Boltzmann Method (LBM) are increasingly popular and attractive methods that propose efficient multiphase formulations, each one with its own strengths and weaknesses. In this context, when it comes to study a given multi-fluid problem, it is helpful to rely on a quantitative comparison to decide which approach should be used and in which context. In particular, the simulation of intermittent two-phase flows in pipes such as slug flows is a complex problem involving moving and intersecting interfaces for which both SPH and LBM could be considered. It is a problem of interest in petroleum applications since the formation of slug flows that can occur in submarine pipelines connecting the wells to the production facility can cause undesired behaviors with hazardous consequences. In this work, we compare SPH and LBM multiphase formulations where surface tension effects are modeled respectively using the continuum surface force and the color gradient approaches on a collection of standard test cases, and on the simulation of intermittent flows in 2D. This paper aims to highlight the contributions and limitations of SPH and LBM when applied to these problems. First, we compare our implementations on static bubble problems with different density and viscosity ratios. Then, we focus on gravity driven simulations of slug flows in pipes for several Reynolds numbers. Finally, we conclude with simulations of slug flows with inlet/outlet boundary conditions. According to the results presented in this study, we confirm that the SPH approach is more robust and versatile whereas the LBM formulation is more accurate and faster

Draping of Cluster Magnetic Fields over Bullets and Bubbles -- Morphology and Dynamic Effects

High-resolution X-ray observations have revealed cavities and `cold fronts' with sharp edges in temperature, density, and metallicity within galaxy clusters. Their presence poses a puzzle since these features are not expected to be hydrodynamically stable, or to remain sharp in the presence of diffusion. However, a moving core or bubble in even a very weakly magnetized plasma necessarily sweeps up enough magnetic field to build up a dynamically important sheath around the object; the layer's strength is set by a competition between `plowing up' of field and field lines slipping around the core. We show that a two-dimensional approach to the problem is quite generally not possible. In three dimensions, we show with analytic arguments and in numerical experiments, that this magnetic layer modifies the dynamics of a plunging core, greatly modifies the effects of hydrodynamic instabilities on the core, modifies the geometry of stripped material, and even slows the fall of the core through magnetic tension. We derive an expression for the maximum magnetic field strength, the thickness of the layer, and the opening angle of the magnetic wake. The morphology of the magnetic draping layer implies the suppression of thermal conduction across the layer, thus conserving strong temperature gradients over the contact surface. The intermittent amplification of the magnetic field as well as the injection of MHD turbulence in the wake of the core is identified to be due to vorticity generation within the magnetic draping layer. These results have important consequences for understanding the physical properties and the complex gasdynamical processes of the intra-cluster medium, and apply quite generally to motions through other magnetized environments, e.g., the ISM.Comment: For version of this paper with interactive 3D graphics and full-resolution figures, see http://www.cita.utoronto.ca/~ljdursi/draping/ . 19p, 26 figures, emulateapj format. Version accepted by ApJ - new references, improved figure

An interface capturing method for liquid-gas flows at low-Mach number

Multiphase, compressible and viscous flows are of crucial importance in a wide range of scientific and engineering problems. Despite the large effort paid in the last decades to develop accurate and efficient numerical techniques to address this kind of problems, current models need to be further improved to address realistic applications. In this context, we propose a numerical approach to the simulation of multiphase, viscous flows where a compressible and an incompressible phase interact in the low-Mach number regime. In this frame, acoustics is neglected but large density variations of the compressible phase can be accounted for as well as heat transfer, convection and diffusion processes. The problem is addressed in a fully Eulerian framework exploiting a low-Mach number asymptotic expansion of the Navier-Stokes equations. A Volume of Fluid approach (VOF) is used to capture the liquid-gas interface, built on top of a massive parallel solver, second order accurate both in time and space. The second-order-pressure term is treated implicitly and the resulting pressure equation is solved with the eigenexpansion method employing a robust and novel formulation. We provide a detailed and complete description of the theoretical approach together with information about the numerical technique and implementation details. Results of benchmarking tests are provided for five different test cases