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.     

Extension of a high‐resolution scheme to 1D liquid–gas flow

This paper presents the extension of a high‐resolution conservative scheme to the one‐dimensional one‐pressure six‐equation two‐fluid flow model. Only mixtures of water and air have been considered in this study, both fluids have been characterized using simple equations of state, namely stiffened gas for the liquid phase and perfect gas for the gas phase. The resulting scheme is explicit and first‐order accurate in space and time. A second‐order version of the scheme has also been derived using the MUSCL strategy and slope limiters. Some numerical results show the good capabilities of this type of schemes in the solution of discontinuities in two‐fluid flow problems, all of them are based on water/air numerical benchmarks widely used in the two‐phase flow literature. Copyright © 2005 John Wiley & Sons, Ltd.     

Preconditioned Krylov subspace methods used in solving two‐dimensional transient two‐phase flows

This paper investigates the performance of preconditioned Krylov subspace methods used in a previously presented two‐fluid model developed for the simulation of separated and intermittent gas–liquid flows. The two‐fluid model has momentum and mass balances for each phase. The equations comprising this model are solved numerically by applying a two‐step semi‐implicit time integration procedure. A finite difference numerical scheme with a staggered mesh is used. Previously, the resulting linear algebraic equations were solved by a Gaussian band solver. In this study, these algebraic equations are also solved using the generalized minimum residual (GMRES) and the biconjugate gradient stabilized (Bi‐CGSTAB) Krylov subspace iterative methods preconditioned with incomplete LU factorization using the ILUT(p, τ) algorithm. The decrease in the computational time using the iterative solvers instead of the Gaussian band solver is shown to be considerable. Copyright © 1999 John Wiley & Sons, Ltd.     

Parallel finite element simulation of mooring forces on floating objects

The coupling between the equations governing the free‐surface flows, the six degrees of freedom non‐linear rigid body dynamics, the linear elasticity equations for mesh‐moving and the cables has resulted in a fluid‐structure interaction technology capable of simulating mooring forces on floating objects. The finite element solution strategy is based on a combination approach derived from fixed‐mesh and moving‐mesh techniques. Here, the free‐surface flow simulations are based on the Navier–Stokes equations written for two incompressible fluids where the impact of one fluid on the other one is extremely small. An interface function with two distinct values is used to locate the position of the free‐surface. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian–Eulerian domain. This allows us to handle the motion of the time dependent geometries. Forces and momentums exerted on the floating object by both water and hawsers are calculated and used to update the position of the floating object in time. In the mesh moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The non‐linear rigid body dynamics equations are coupled with the governing equations of fluid flow and are solved simultaneously to update the position of the floating object. The numerical examples includes a 3D simulation of water waves impacting on a moored floating box and a model boat and simulation of floating object under water constrained with a cable. Copyright © 2003 John Wiley & Sons, Ltd.     

Implementation of phase change thermodynamic probability for unsteady simulation of cavitating flows

The aim of this work is to investigate the non‐equilibrium effects of phase change in cavitating flows. For this purpose, the concept of phase change thermodynamic probability is used along with homogeneous model to simulate two‐phase cavitating flows. For simulation of unsteady behaviors of cavitation, which have practical applications, unsteady PISO algorithm based on the non‐conservative approach is utilized. For multi‐phase simulation, single‐fluid Navier–Stokes equations, along with the volume fraction transport equation, are employed. In this paper, phase change thermodynamics probabilities and cavitation model is briefly summarized. Thus, derivation of the cavitation model, starting from the basic thermodynamic equations to the mass and momentum conservation equations at a liquid–vapor two‐phase flow, is presented to explain the numerical model. Unsteady simulations of cavitation around a flat plate normal to flow direction are presented to clarify the accuracy of the model. The accuracy of the numerical results is good, and it is possible to apply this method to more complex geometries. Copyright © 2010 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.     

Founding editor Dr Philip M. Gresho

A numerical technique is presented for the approximation of vertical gradient of the non‐hydrostatic pressure arising in the Reynolds‐averaged Navier–Stokes equations for simulating non‐hydrostatic free‐surface flows. It is based on the Keller‐box method that take into account the effect of non‐hydrostatic pressure with a very small number of vertical grid points. As a result, the proposed technique is capable of simulating relatively short wave propagation, where both frequency dispersion and non‐linear effects play an important role, in an accurate and efficient manner. Numerical examples are provided to illustrate this; accurate wave characteristics are already achieved with only two layers. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of the unsteady behaviour of cavitating flows

A 2D numerical model is proposed to simulate unsteady cavitating flows. The Reynolds‐averaged Navier–Stokes equations are solved for the mixture of liquid and vapour, which is considered as a single fluid with variable density. The vapourization and condensation processes are controlled by a barotropic state law that relates the fluid density to the pressure variations. The numerical resolution is a pressure‐correction method derived from the SIMPLE algorithm, with a finite volume discretization. The standard scheme is slightly modified to take into account the cavitation phenomenon. That numerical model is used to calculate unsteady cavitating flows in two Venturi type sections. The choice of the turbulence model is discussed, and the standard RNG k–εmodel is found to lead to non‐physical stable cavities. A modified k–εmodel is proposed to improve the simulation. The influence of numerical and physical parameters is presented, and the numerical results are compared to previous experimental observations and measurements. The proposed model seems to describe the unsteady cavitation behaviour in 2D geometries well. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical simulation of 3D free surface flows,with multiple incompressible immiscible phases. Applications to impulse waves

A numerical method for the solution to the density‐dependent incompressible Navier–Stokes equations modeling the flow of N immiscible incompressible liquid phases with a free surface is proposed. It allows to model the flow of an arbitrary number of liquid phases together with an additional vacuum phase separated with a free surface. It is based on a volume‐of‐fluid approach involving N indicator functions (one per phase, identified by its density) that guarantees mass conservation within each phase. An additional indicator function for the whole liquid domain allows to treat boundary conditions at the interface between the liquid domain and a vacuum. The system of partial differential equations is solved by implicit operator splitting at each time step: first, transport equations are solved by a forward characteristics method on a fine Cartesian grid to predict the new location of each liquid phase; second, a generalized Stokes problem with a density‐dependent viscosity is solved with a FEM on a coarser mesh of the liquid domain. A novel algorithm ensuring the maximum principle and limiting the numerical diffusion for the transport of the N phases is validated on benchmark flows. Then, we focus on a novel application and compare the numerical and physical simulations of impulse waves, that is, waves generated at the free surface of a water basin initially at rest after the impact of a denser phase. A particularly useful application in hydraulic engineering is to predict the effects of a landslide‐generated impulse wave in a reservoir. Copyright © 2014 John Wiley & Sons, Ltd.     

NUMERICAL PREDICTION OF TWO-PHASE FLOW IN BUBBLE COLUMNS

A numerical model is described for the prediction of turbulent continuum equations for two-phase gas–liquid flows in bubble columns. The mathematical formulation is based on the solution of each phase. The two-phase model incorporates interfacial models of momentum transfer to account for the effects of virtual mass, lift, drag and pressure discontinuities at the gas–liquid interface. Turbulence is represented by means of a two-equation k–ϵ model modified to account for bubble-induced turbulence production. The numerical discretization is based on a staggered finite-volume approach, and the coupled equations are solved in a segregated manner using the IPSA method. The model is implemented generally in the multipurpose PHOENICS computer code, although the present appllications are restricted to two-dimensional flows. The model is applied to simulate two bubble column geometries and the predictions are compared with the measured circulation patterns and void fraction distributions.     

A direct numerical simulation of axisymmetric liquid–gas two‐phase laminar flows in a pipe

An accurate finite‐volume‐based numerical method for the simulation of an isothermal two‐phase flow, consisting of a liquid slug translating in a non‐reacting gas in a circular pipe is presented. This method is built on a sharp interface concept and developed on an Eulerian Cartesian fixed‐grid with a cut‐cell scheme and marker points to track the moving interface. The unsteady, axisymmetric Navier–Stokes equations in both liquid and gas phases are solved separately. The mass continuity and momentum flux conditions are explicitly matched at the true surface phase boundary to determine the interface shape and movement. A quadratic curve fitting algorithm with marker points is used to yield smooth and accurate information of the interface curvatures. It is uniquely demonstrated for the first time with the current method that conservation of mass is strictly enforced for continuous infusion of flow into the domain of computation. The method has been used to compute the velocity and pressure fields and the deformation of the liquid core. It is also shown that the current method is capable of producing accurate results for a wide range of Reynolds number, Re, Weber number, We, and large property jumps at the interface. Copyright © 2010 John Wiley & Sons, Ltd.     

A compressible single‐temperature conservative two‐phase model with phase transitions

A model for multidimensional compressible two‐phase flow with pressure and velocity relaxations based on the theory of thermodynamically compatible system is extended to study liquid–gas flows with cavitation. The model assumes for each phase its own pressure and velocity, while a common temperature is considered. The governing equations form a hyperbolic system in conservative form and are derived through the theory of a thermodynamically compatible system. The phase pressure‐equalizing process and the interfacial friction are taken into account in the balance laws for the volume fractions of one phase and for the relative velocity by adding two relaxation source terms, while the phase transition is introduced into the model considering in the balance equation for the mass of one phase the relaxation of the Gibbs free energies of the two phases. A modification of the central finite‐volume Kurganov–Noelle–Petrova method is adopted in this work to solve the homogeneous hyperbolic part, while the relaxation source terms are treated implicitly. In order to investigate the effect of the mass transfer in the solution, a 1D cavitation tube problem is presented. In addition, two 2D numerical simulations regarding cavitation problem are also studied: a cavitating Richtmyer–Meshkov instability and a laser‐induced cavitation problem. Copyright © 2014 John Wiley & Sons, Ltd.     

Slip effects on heat and mass transfer in MHD micropolar fluid flow over an inclined plate with thermal radiation and chemical reaction

The influence of partial slip, thermal radiation, chemical reaction and temperature‐dependent fluid properties on heat and mass transfer in hydro‐magnetic micropolar fluid flow over an inclined permeable plate with constant heat flux and non‐uniform heat source/sink is studied. The transverse magnetic field is assumed as a function of the distance from the origin. Also it is assumed that the fluid viscosity and the thermal conductivity vary as an inverse function and linear function of temperature, respectively. With the use of the similarity transformation, the governing system of non‐linear partial differential equations are transformed into non‐linear ordinary differential equations and are solved numerically using symbolic software MATHEMATICA 7.0 (Wolfram Research, Champaign, IL). The numerical values obtained for the velocity, microrotation, temperature, species concentration, skin friction coefficient and the Nusselt number are presented through graphs and tables for several sets of values of the parameters. The effects of various physical parameters on the flow and heat transfer characteristics are discussed.Copyright © 2011 John Wiley & Sons, Ltd.     

The numerical study of roll‐waves in inclined open channels and solitary wave run‐up

In this paper, we introduce a finite‐volume kinetic BGK scheme and its applications to the study of roll and solitary waves. The current scheme is based on the numerical solution of the gas‐kinetic Bhatnagar–Gross–Krook model in the flux evaluation across each cell interface. An intrinsic connection between the BGK model and time‐dependent, non‐linear, non‐homogeneous shallow‐water equations enables us to solve shallow‐water equations automatically with our kinetic scheme. The analytical solution, experimental measurements, and numerical calculations for problems associated with roll‐waves down an inclined open channel and solitary waves incident on a sloped beach are also presented. Copyright © 2005 John Wiley & Sons, Ltd.     

High‐order compact numerical schemes for non‐hydrostatic free surface flows

The development of a numerical scheme for non‐hydrostatic free surface flows is described with the objective of improving the resolution characteristics of existing solution methods. The model uses a high‐order compact finite difference method for spatial discretization on a collocated grid and the standard, explicit, single step, four‐stage, fourth‐order Runge–Kutta method for temporal discretization. The Cartesian coordinate system was used. The model requires the solution of two Poisson equations at each time‐step and tridiagonal matrices for each derivative at each of the four stages in a time‐step. Third‐ and fourth‐order accurate boundaries for the flow variables have been developed including the top non‐hydrostatic pressure boundary. The results demonstrate that numerical dissipation which has been a problem with many similar models that are second‐order accurate is practically eliminated. A high accuracy is obtained for the flow variables including the non‐hydrostatic pressure. The accuracy of the model has been tested in numerical experiments. In all cases where analytical solutions are available, both phase errors and amplitude errors are very small. Copyright © 2006 John Wiley & Sons, Ltd.     

A finite element thermally coupled flow formulation for phase‐change problems

A finite element, thermally coupled incompressible flow formulation considering phase‐change effects is presented. This formulation accounts for natural convection, temperature‐dependent material properties and isothermal and non‐isothermal phase‐change models. In this context, the full Navier–Stokes equations are solved using a generalized streamline operator (GSO) technique. The highly non‐linear phase‐change effects are treated with a temperature‐based algorithm, which provides stability and convergence of the numerical solution. The Boussinesq approximation is used in order to consider the temperature‐dependent density variation. Furthermore, the numerical solution of the coupled problem is approached with a staggered incremental‐iterative solution scheme, such that the convergence criteria are written in terms of the residual vectors. Finally, this formulation is used for the solutions of solidification and melting problems validating some numerical results with other existing solutions obtained with different methodologies. Copyright © 2000 John Wiley & Sons, Ltd.     

Two‐phase pressure drop caused by sudden flow area contraction/expansion in small circular pipes

Theoretical studies have been made to determine the pressure drops caused by abrupt flow area expansion/contraction in small circular pipes for two‐phase flow of air and water mixtures at room temperature and near atmospheric pressure. Two‐phase computational fluid dynamics (CFD) calculations, using Eulerian–Eulerian model (with the air phase being compressible for pipe contraction case) are employed to calculate the pressure drop across sudden expansion and contraction. The pressure drop is determined by extrapolating the computed pressure profiles upstream and downstream of the expansion/contraction. The larger and smaller tube diameters are 1.6 and 0.84 mm, respectively. Computations have been performed with single‐phase water and air, and two‐phase mixtures in a range of Reynolds number (considering all‐liquid flow) from 1000 to 12 000 and flow quality from 1.2 × 10^?3 to 1.6 × 10^?2. The numerical results are validated against experimental data from the literature and are found to be in good agreement. The expansion and contraction loss coefficients are found to be different for single‐phase flow of air and water, and they agreed reasonably well with the commonly used theoretical predictions. Based on the numerical results as well as experimental data, correlations are developed for two‐phase flow pressure drops caused by the flow area contraction as well as expansion. Copyright © 2010 John Wiley & Sons, Ltd.     

A high‐precision unstructured adaptive mesh technique for gas–liquid two‐phase flows

Adaptive mesh techniques are used widely in the numerical simulations of fluid flows, and the simulation results with high accuracies are obtained by appropriate mesh adaptations. However, gas–liquid two‐phase flows are still difficult to be simulated on adaptive meshes, especially on unstructured adaptive meshes, because the physical phenomena near gas–liquid interfaces are highly complicated and in general, not modeled appropriately on adaptive meshes. In this paper, a high‐precision unstructured adaptive mesh technique for gas–liquid two‐phase flows is developed and verified/validated. In the unstructured adaptive mesh technique, the PLIC algorithm is employed to simulate interfacial dynamic behaviors and, therefore, the reconstruction method for the interfaces in refined cells is developed, which satisfies the gas and liquid volume conservations and geometrical conservations of interfaces. In addition, the physics‐based consideration is performed on the momentum calculations near interfaces, and the calculation method with gas and liquid momentum conservations is developed. For verification, the slotted‐disk revolution problem is solved. As a result, the unstructured adaptive mesh technique succeeds in reproducing the slotted‐disk shape accurately and well maintaining the shape after one full‐revolution. The dam‐break problem is also simulated and the momentum conservative calculation method succeeds in providing physically appropriate results, which show good agreements with experimental data. Therefore, it is confirmed that the developed unstructured adaptive mesh technique is very efficient to simulate gas–liquid two‐phase flows accurately. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulation of sharp gas–liquid interface using VOF method and adaptive grid local refinement around the interface

The volume of fluid (VOF) method is used to perform two‐phase simulations (gas–liquid). The governing Navier–Stokes conservation equations of the flow field are numerically solved on two‐dimensional axisymmetric or three‐dimensional unstructured grids, using Cartesian velocity components, following the finite volume approximation and a pressure correction method. A new method of adaptive grid local refinement is developed in order to enhance the accuracy of the predictions, to capture the sharp gas–liquid interface and to speed up the calculations. Results are compared with experimental measurements in order to assess the efficiency of the method. Copyright © 2004 John Wiley & Sons, Ltd.     

On the efficiency of linearization schemes and coupled multigrid methods in the simulation of a 3D flow around a cylinder

This paper studies the efficiency of two ways to treat the non‐linear convective term in the time‐dependent incompressible Navier–Stokes equations and of two multigrid approaches for solving the arising linear algebraic saddle point problems. The Navier–Stokes equations are discretized by a second‐order implicit time stepping scheme and by inf–sup stable, higher order finite elements in space. The numerical studies are performed at a 3D flow around a cylinder. Copyright © 2005 John Wiley & Sons, Ltd.