Bubble collapse in compressible fluids using a spectral element marker particle method. Part 2. Viscoelastic fluids

This paper is concerned with the development of a high‐order numerical scheme for two‐phase viscoelastic flows. In the companion paper, herein referred to as Part 1, the scheme is applied to the modelling of two‐phase Newtonian flows. The particular problem of the collapse of a 2D bubble in the vicinity of a rigid boundary is considered. Attention is given to the construction of the most general form of the compressible Oldroyd B model that is consistent with the compressible Newtonian and upper‐convected Maxwell models in the appropriate limits. The governing equations are discretized using the spectral element method, and the two phases are modelled using a marker particle method. A comprehensive set of results is presented for the problem of bubble collapse near a rigid wall, and qualitative agreement is obtained with other numerical studies and experimental observations. Viscoelastic effects that are predicted include increased bubble oscillation with increasing Weissenberg number and considerable bubble deformation and cusping near the wall. Most importantly, it has been shown that viscoelasticity has the ability to prevent jet formation and therefore is likely to have a mitigating effect on cavitation damage. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of a single bubble by compressible two‐phase fluids

The present work deals with the numerical investigation of a collapsing bubble in a liquid–gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non‐conservative transport equations of the pressure law coefficients. The original first‐order discretization is extended to higher order applying second‐order ENO reconstruction to the primitive variables. The derivation of the non‐conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale‐based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2‐D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser‐induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd.     

A numerical scheme for Euler–Lagrange simulation of bubbly flows in complex systems

An Eulerian–Lagrangian approach is developed for the simulation of turbulent bubbly flows in complex systems. The liquid phase is treated as a continuum and the Navier–Stokes equations are solved in an unstructured grid, finite volume framework for turbulent flows. The dynamics of the disperse phase is modeled in a Lagrangian frame and includes models for the motion of each individual bubble, bubble size variations due to the local pressure changes, and interactions among the bubbles and with boundaries. The bubble growth/collapse is modeled by the Rayleigh–Plesset (RP) equation. Three modeling approaches are considered: (a) one‐way coupling, where the influence of the bubble on the fluid flow is neglected, (b) two‐way coupling, where the momentum‐exchange between the fluid and the bubbles is modeled, and (c) volumetric coupling, where the volumetric displacement of the fluid by the bubble motion and the momentum‐exchange are modeled. A novel adaptive time‐stepping scheme based on stability‐analysis of the non‐linear bubble dynamics equations is developed. The numerical approach is verified for various single bubble test cases to show second‐order accuracy. Interactions of multiple bubbles with vortical flows are simulated to study the effectiveness of the volumetric coupling approach in predicting the flow features observed experimentally. Finally, the numerical approach is used to perform a large‐eddy simulation in two configurations: (i) flow over a cavity to predict small‐scale cavitation and inception and (ii) a rising dense bubble plume in a stationary water column. The results show good predictive capability of the numerical algorithm in capturing complex flow features. Copyright © 2010 John Wiley & Sons, Ltd.     

Development of new finite volume schemes on unstructured triangular grid for simulating the gas–liquid two‐phase flow

This paper presents a coupled finite volume inner doubly iterative efficient algorithm for linked equations (IDEAL) with level set method to simulate the incompressible gas–liquid two‐phase flows with moving interfaces on unstructured triangular grid. The finite volume IDEAL method on a collocated grid is employed to solve the incompressible two‐phase Navier–Stokes equations, and the level set method is used to capture the moving interfaces. For the sake of mass conservation, an effective second‐order accurate finite volume scheme is developed to solve the level set equation on triangular grid, which can be implemented much easier than the classical high‐order level set solvers. In this scheme, the value of level set function on the boundary of control volume is approximated using a linear combination of a high‐order Larangian interpolation and a second‐order upwind interpolation. By the rotating slotted disk and stretching and shrinking of a circular fluid element benchmark cases, the mass conservation and accuracy of the new scheme is verified. Then the coupled method is applied to two‐phase flows, including a 2D bubble rising problem and a 2D dam breaking problem. The computational results agree well with those reported in literatures and experimental data. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical analysis of gas bubbles in close proximity to a movable or deformable body

The growth and collapse of gaseous bubbles near a movable or deformable body are investigated numerically using the boundary element method and fluid–solid coupling technique. The fluid is treated as inviscid, incompressible and the flow irrotational. The unsteady Bernoulli equation is applied on the bubble surface as one of the boundary conditions of the Laplace’s equation for the potential. Good agreements between the numerical and experimental results demonstrate the robustness and accuracy of the present method. The translation and rotation of the rigid body due to the bubble evolution are captured by solving the six-degrees-of-freedom equations of motion for the rigid body. The fluid–solid coupling is achieved by matching the normal component of the velocity and the pressure at the fluid–solid interface. Compared to a fixed rigid body, the expansion of the bubble is not affected too much but much faster collapsing velocities during the collapsing phase of bubble can be observed when considering the motion of the rigid body. The rigid body is pushed away as the bubble grows and moved toward the bubble as the bubble collapses. The motion of two bubbles near a movable cylinder is also simulated. The large rotation of the cylinder and obvious deformation and distortion for the bubble in close proximity to a curved wall are observed in our codes. Finally, the growth and collapse of bubble near a deformable ellipsoid shell are also simulated using the combination of boundary element method (BEM) and finite element method (FEM) techniques. The oscillations of the ellipsoid shell can be observed during the growth and collapse of bubble, which much differs from the results obtained by only considering effects of a rigidly movable body on the bubble evolution.     

An alternate approach for modelling of transient vaporous cavitation

Simulation of cavitating flow has been a thrust area of research for long period due to its practical and economic importance. The major hurdle in developing a numerical model for such flows is the difficulty in representing the quick phase changes, in general, and the alternate change of flow from single phase to two phase and back, in particular. In this case, instability due to sharp variation of flow characteristics also restricts the development of numerical models. The present study demonstrates the use of a relatively simple formulation for the analysis of flow characteristics in a quasi‐rigid pipeline under abrupt phase changes due to cavitation. A popular scheme—MacCormack scheme—was used for developing a numerical solution for this problem. It uses the conservative form of the governing equations, viz. conservation of mass and momentum, the transport equation and the constitutive relationship. The model can handle variable properties of the water–vapour mixture, which is highly compressible. A newly introduced pressure under‐relaxation method overcomes the numerical instability due to sharp variation of flow characteristics during phase change. The model could predict the instant of occurrence of vapour pressure, duration of persistence of vapour pressure and the rise of pressure due to vapour collapse to satisfactory levels with published data and experimental results. Copyright © 2009 John Wiley & Sons, Ltd.     

A numerical continuous model for the hydrodynamics of fluid particle systems

In order to understand the hydrodynamic interactions that can appear in a fluid particle motion, an original method based on the equations governing the motion of two immiscible fluids has been developed. These momentum equations are solved for both the fluid and solid phases. The solid phase is assumed to be a fluid phase with physical properties, such as its behaviour can be assimilated to that of pseudo‐rigid particles. The only unknowns are the velocity and the pressure defined in both phases. The unsteady two‐dimensional momentum equations are solved by using a staggered finite volume formulation and a projection method. The transport of each particle is solved by using a second‐order explicit scheme. The physical model and the numerical method are presented, and the method is validated through experimental measurements and numerical results concerning the flow around a circular cylinder. Good agreement is observed in most cases. The method is then applied to study the trajectory of one settling particle initially off‐centred between two parallel walls and the corresponding wake effects. Different particle trajectories related to particulate Reynolds numbers are presented and commented. A two‐body interaction problem is investigated too. This method allows the simulation of the transport of particles in a dilute suspension in reasonable time. One of the important features of this method is the computational cost that scales linearly with the number of particles. Copyright © 1999 John Wiley & Sons, Ltd.     

The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

A generalized Rusanov method for the Baer‐Nunziato equations with application to DDT processes in condensed porous explosives

The paper addresses a numerical approach for solving the Baer‐Nunziato equations describing compressible 2‐phase flows. We are developing a finite‐volume method where the numerical flux is approximated with the Godunov scheme based on the Riemann problem solution. The analytical solution to this problem is discussed, and approximate solvers are considered. The obtained theoretical results are applied to develop the discrete model that can be treated as an extension of the Rusanov numerical scheme to the Baer‐Nunziato equations. Numerical results are presented that concern the method verification and also application to the deflagration‐to‐detonation transition (DDT) in porous reactive materials.     

Unified one‐fluid formulation for incompressible flexible solids and multiphase flows: Application to hydrodynamics using the immersed structural potential method (ISPM)

In this paper, we present a two‐dimensional computational framework for the simulation of fluid‐structure interaction problems involving incompressible flexible solids and multiphase flows, further extending the application range of classical immersed computational approaches to the context of hydrodynamics. The proposed method aims to overcome shortcomings such as the restriction of having to deal with similar density ratios among different phases or the restriction to solve single‐phase flows. First, a variation of classical immersed techniques, pioneered with the immersed boundary method (IBM), is presented by rearranging the governing equations, which define the behaviour of the multiple physics involved. The formulation is compatible with the “one‐fluid” formulation for two‐phase flows and can deal with large density ratios with the help of an anisotropic Poisson solver. Second, immersed deformable structures and fluid phases are modelled in an identical manner except for the computation of the deviatoric stresses. The numerical technique followed in this paper builds upon the immersed structural potential method developed by the authors, by adding a level set–based method for the capturing of the fluid‐fluid interfaces and an interface Lagrangian‐based meshless technique for the tracking of the fluid‐structure interface. The spatial discretisation is based on the standard marker‐and‐cell method used in conjunction with a fractional step approach for the pressure/velocity decoupling, a second‐order time integrator, and a fixed‐point iterative scheme. The paper presents a wide d range of two‐dimensional applications involving multiphase flows interacting with immersed deformable solids, including benchmarking against both experimental and alternative numerical schemes.     

Sliding of fine particles on the slip surface of rising gas bubbles: Resistance of liquid shear flows

In this paper a model was developed to describe the shear flow resistance force and torque acting on a fine particle as it slides on the slip surface of a rising gas bubble. The shear flow close to the bubble surface was predicted using a Taylor series and the numerical data obtained from the Navier–Stokes equations as a function of the polar coordinates at the bubble surface, the bubble Reynolds number, and the gas hold-up. The particle size was considered to be sufficiently small relative to the bubble size that the bubble surface could be locally approximated to a planar interface. The Stokes equation for the disturbance shear flows was solved for the velocity components and pressure using series of bispherical coordinates and the boundary conditions at the no-slip particle surface and the slip bubble surface. The solutions for the disturbance flows were then used to calculate the flow resistance force and torque on the particle as a function of the separation distance between the bubble and particle surfaces. The resistance functions were determined by dividing the actual force and torque by the corresponding (Stokes) force and torque in the bulk phase. Finally, numerical and simplified analytical rational approximate solutions for force correction factors for sliding particles as a function of the (whole range of the) separation distance are presented, which are in good agreement with the exact numerical result and can be readily applied to more general modelling of the bubble–particle interactions.     

The influence of viscoelasticity on the collapse of cavitation bubbles near a rigid boundary

An understanding of the phenomena associated with cavitation is important in many areas of science and engineering. This paper is concerned with the influence of viscoelasticity on the dynamics of cavitation bubbles near rigid boundaries. Viscoelastic effects are modelled using a Maxwell constitutive equation, and a generalized Bernoulli equation is derived. The governing equations are solved using the boundary element method in which both the bubble surface and the potential are represented by cubic splines. The numerical scheme is validated through comparisons with results in the literature for the inviscid case. The introduction of viscoelasticity introduces some interesting bubble dynamics including the occurrence of oscillations during collapse. Most importantly, it is shown that viscoelasticity can serve to suppress the formation of a liquid jet. The subsequent reduced pressures compared with the inviscid case suggest that viscoelasticity has a mitigating effect on cavitation damage.     

On high‐resolution schemes for solving unsteady compressible two‐phase dilute viscous flows

A high‐resolution numerical scheme based on the MUSCL–Hancock approach is developed to solve unsteady compressible two‐phase dilute viscous flow. Numerical considerations for the development of the scheme are provided. Several solvers for the Godunov fluxes are tested and the results lead to the choice of an exact Riemann solver adapted for both gaseous and dispersed phases. The accuracy of the scheme is proven step by step through specific test cases. These simulations are for one‐phase viscous flows over a flat plate in subsonic and supersonic regimes, unsteady flows in a low‐pressure shock tube, two‐phase dilute viscous flows over a flat plate and, finally, two‐phase unsteady viscous flows in a shock tube. The results are compared with well‐established analytical and numerical solutions and very good agreement is achieved. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical simulation of bubble and droplet deformation by a level set approach with surface tension in three dimensions

In this paper we present a three‐dimensional Navier–Stokes solver for incompressible two‐phase flow problems with surface tension and apply the proposed scheme to the simulation of bubble and droplet deformation. One of the main concerns of this study is the impact of surface tension and its discretization on the overall convergence behavior and conservation properties. Our approach employs a standard finite difference/finite volume discretization on uniform Cartesian staggered grids and uses Chorin's projection approach. The free surface between the two fluid phases is tracked with a level set (LS) technique. Here, the interface conditions are implicitly incorporated into the momentum equations by the continuum surface force method. Surface tension is evaluated using a smoothed delta function and a third‐order interpolation. The problem of mass conservation for the two phases is treated by a reinitialization of the LS function employing a regularized signum function and a global fixed point iteration. All convective terms are discretized by a WENO scheme of fifth order. Altogether, our approach exhibits a second‐order convergence away from the free surface. The discretization of surface tension requires a smoothing scheme near the free surface, which leads to a first‐order convergence in the smoothing region. We discuss the details of the proposed numerical scheme and present the results of several numerical experiments concerning mass conservation, convergence of curvature, and the application of our solver to the simulation of two rising bubble problems, one with small and one with large jumps in material parameters, and the simulation of a droplet deformation due to a shear flow in three space dimensions. Furthermore, we compare our three‐dimensional results with those of quasi‐two‐dimensional and two‐dimensional simulations. This comparison clearly shows the need for full three‐dimensional simulations of droplet and bubble deformation to capture the correct physical behavior. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical study of flows of complex fluids between eccentric cylinders using transformation functions

In this paper, we investigate fluid flows between eccentric cylinders by means of two stream‐tube analyses. The first method considers a one‐to‐one global transformation function that allows the physical domain to be transformed into a mapped domain, used as computational domain, that involves concentric streamlines. The second approach uses local transformations and domain decomposition techniques to deal with mixed flow regimes. Both formulations are particularly adapted for handling time‐dependent constitutive equations, since particle‐tracking problems are avoided. Mass conservation is verified in both formulations and the relevant numerical procedure can be carried out using simple meshes built on the mapped streamlines. Fluids obeying anelastic and viscoelastic constitutive equations are considered in the calculations. The numerical results are consistent with those in the literature for the flow rates tested. Application of the method to the K‐BKZ memory‐integral constitutive equation highlights significant differences between the model predictions and those provided by more simple rheological models. Copyright © 2002 John Wiley & Sons, Ltd.     

Compressible simulations of bubble dynamics with central-upwind schemes

This paper discusses the implementation of an explicit density-based solver, that utilises the central-upwind schemes for the simulation of cavitating bubble dynamic flows. It is highlighted that, in conjunction with the Monotonic Upstream-Centered Scheme for Conservation Laws (MUSCL) scheme they are of second order in spatial accuracy; essentially they are high-order extensions of the Lax–Friedrichs method and are linked to the Harten Lax and van Leer (HLL) solver family. Basic comparison with the predicted wave pattern of the central-upwind schemes is performed with the exact solution of the Riemann problem, for an equation of state used in cavitating flows, showing excellent agreement. Next, the solver is used to predict a fundamental bubble dynamics case, the Rayleigh collapse, in which results are in accordance to theory. Then several different bubble configurations were tested. The methodology is able to handle the large pressure and density ratios appearing in cavitating flows, giving similar predictions in the evolution of the bubble shape, as the reference.     

On the modeling and simulation of a laser‐induced cavitation bubble

Single cavitation bubbles exhibit severe modeling and simulation difficulties. This is due to the small scales of time and space as well as due to the involvement of different phenomena in the dynamics of the bubble. For example, the compressibility, phase transition, and the existence of a noncondensable gas inside the bubble have strong effects on the dynamics of the bubble. Moreover, the collapse of the bubble involves the occurrence of critical conditions for the pressure and temperature. This adds extra difficulties to the choice of equations of state. Even though several models and simulations have been used to study the dynamics of the cavitation bubbles, many details are still not clearly accounted for. Here, we present a numerical investigation for the collapse and rebound of a laser‐induced cavitation bubble in liquid water. The compressibility of the liquid and vapor are involved. In addition, great focus is devoted to study the effects of phase transition and the existence of a noncondensable gas on the dynamics of the collapsing bubble. If the bubble contains vapor only, we use the six‐equation model for two‐phase flows that was modified in our previous work [A. Zein, M. Hantke, and G. Warnecke, J. Comput. Phys., 229(8):2964‐2998, 2010]. This model is an extension to the six‐equation model with a single velocity of Kapila et al. (Phys. Fluid, 13:3002‐3024, 2001) taking into account the heat and mass transfer. To study the effect of a noncondensable gas inside the bubble, we add a third phase to the original model. In this case, the phase transition is considered only at interfaces that separate the liquid and its vapor. The stiffened gas equations of state are used as closure relations. We use our own method to determine the parameters to obtain reasonable equations of state for a wide range of temperatures and make them suitable for the phase transition effects. We compare our results with experimental ones. Also our results confirm some expected physical phenomena. Copyright © 2013 John Wiley & Sons, Ltd.     

Development of a hybrid particle‐mesh method for two‐phase flow simulations

A hybrid particle‐mesh method was developed for efficient and accurate simulations of two‐phase flows. In this method, the main component of the flow is solved using the constrained interpolated profile/multi‐moment finite volumemethod; the two‐phase interface is rendered using the finite volume particle (FVP) method. The effect of surface tension is evaluated using the continuum surface force model. Numerical particles in the FVP method are distributed only on the surface of the liquid in simulating the interface between liquid and gas; these particles are used to determine the density of each mesh grid. An artificial term was also introduced to mitigate particle clustering in the direction of maximum compression and sparse discretization errors in the stretched direction. This enables accurate interface tracking without diminishing numerical efficiency. Two benchmark simulations are used to demonstrate the validity of the method developed and its numerical stability. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulating surface tension with smoothed particle hydrodynamics

A method for simulating two‐phase flows including surface tension is presented. The approach is based upon smoothed particle hydrodynamics (SPH). The fully Lagrangian nature of SPH maintains sharp fluid–fluid interfaces without employing high‐order advection schemes or explicit interface reconstruction. Several possible implementations of surface tension force are suggested and compared. The numerical stability of the method is investigated and optimal choices for numerical parameters are identified. Comparisons with a grid‐based volume of fluid method for two‐dimensional flows are excellent. The methods presented here apply to problems involving interfaces of arbitrary shape undergoing fragmentation and coalescence within a two‐phase system and readily extend to three‐dimensional problems. Boundary conditions at a solid surface, high viscosity and density ratios, and the simulation of free‐surface flows are not addressed. Copyright © 2000 John Wiley & Sons, Ltd.