Energy spectra and bubble velocity distributions in pseudo-turbulence: Numerical simulations vs. experiments

Direct numerical simulations (DNS) are performed to study the behavior of a swarm of rising air bubbles in water, employing the front tracking method, which allows to handle finite-size bubbles. The swarms consist of monodisperse deformable 4 mm bubbles with a gas fraction of 5% and 15%. This paper focuses on the comparison of the liquid energy spectra and bubble velocity probability density functions (PDFs) with experimental data obtained by phase-sensitive constant-temperature anemometry (CTA) and three-dimensional particle tracking velocimetry (PTV), respectively.     

A probabilistic method for field solutions in two-phase flows

Solid-particle motion and related transport phenomena in two-phase flow are fluctuating processes in space and time. A deterministic method can describe only partially the intrinsic physics of these processes. In this paper, the fluctuations of the flow parameters are modelled by considering the spatial correlations, and a probabilistic computational method for two-phase flow is presented. The probabilistic governing equations have been discretized in space using a finite volume method, and then solved by applying the Neumann expansion method. This last method is time efficient, and its convergence can be guaranteed even for large fluctuations. A liquid-solid particle mixture flow in a circular pipe is taken as an example. Computational results illustrate the merit of the probabilistic approach for the prediction of two-phase flow phenomena.     

A local adaptive grid procedure for incompressible flows with multigridding and equidistribution concepts

An adaptive grid solution procedure is developed for incompressible flow problems in which grid refinement based on an equidistribution law is performed in high-error-estimate regions that are flagged from a preliminary coarse grid solution. Solutions on the locally refined and equidistributed meshes are obtained using boundary conditions interpolated from the preliminary coarse grid solution, and solutions on both the refined and coarse grid regions are successively improved using a multigrid approach. For this purpose, suitable correction terms for the coarse grid equations are derived for all variables in the flagged regions. This procedure with Local Adaptation, Multigridding and Equidistribution (LAME) concepts is applied to various flow problems to demonstrate the accuracy improvements obtained using this method.     

An application of the wavelet analysis technique for the objective discrimination of two-phase flow patterns

This paper presents an application of the wavelet analysis technique for two-phase flow pattern identification by using the void fraction signals obtained from a multi-channel Impedance Void Meter (IVM) in a vertical-upward air–water flow. A new method for the objective discrimination of the two-phase flow pattern has been developed to provide information regarding the local energy of void fraction signals at a given scale on the joint time–frequency diagram. The void signals are processed with Continuous Wavelet Transform (CWT) to get the local wavelet energy coefficients map on the time–frequency diagram. The effective local wavelet energy and the effective scale are then calculated. Then the criteria for flow pattern identification are, finally, obtained. A series of void fraction measurements were conducted over a wide range of air–water vertical-upward flow condition to provide an extensive database to cover several types of flow patterns. The results show that the proposed method has a high precision for characterizing different flow regimes in two-phase flow, and is considerably more promising for the online recognition of two-phase flow patterns due to the short time of data processing.     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

A consistent incompressible SPH method for internal flows with fixed and moving boundaries

An improved incompressible smoothed particle hydrodynamics (ISPH) method is presented, which employs first‐order consistent discretization schemes both for the first‐order and second‐order spatial derivatives. A recently introduced wall boundary condition is implemented in the context of ISPH method, which does not rely on using dummy particles and, as a result, can be applied more efficiently and with less computational complexity. To assess the accuracy and computational efficiency of this improved ISPH method, a number of two‐dimensional incompressible laminar internal flow benchmark problems are solved and the results are compared with available analytical solutions and numerical data. It is shown that using smaller smoothing lengths, the proposed method can provide desirable accuracies with relatively less computational cost for two‐dimensional problems. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupled continuous and discontinuous finite element method for the incompressible flows

In this paper, we develop a coupled continuous Galerkin and discontinuous Galerkin finite element method based on a split scheme to solve the incompressible Navier–Stokes equations. In order to use the equal order interpolation functions for velocity and pressure, we decouple the original Navier–Stokes equations and obtain three distinct equations through the split method, which are nonlinear hyperbolic, elliptic, and Helmholtz equations, respectively. The hybrid method combines the merits of discontinuous Galerkin (DG) and finite element method (FEM). Therefore, DG is concerned to accomplish the spatial discretization of the nonlinear hyperbolic equation to avoid using the stabilization approaches that appeared in FEM. Moreover, FEM is utilized to deal with the Poisson and Helmholtz equations to reduce the computational cost compared with DG. As for the temporal discretization, a second‐order stiffly stable approach is employed. Several typical benchmarks, namely, the Poiseuille flow, the backward‐facing step flow, and the flow around the cylinder with a wide range of Reynolds numbers, are considered to demonstrate and validate the feasibility, accuracy, and efficiency of this coupled method. Copyright © 2016 John Wiley & Sons, Ltd.     

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines. Friction factors play an important role in the accurate calculation of pressure drop. Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall, gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions. In this paper by modification of a friction model available in the literature, an improved semiempirical model is proposed. The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved. Pressure gradient and velocity profiles are validated against experimental data. Using the improved model, the pressure gradient deviation from experiments diminishes by about 3%; the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.     

Numerical stabilities of loosely coupled methods for robust modeling of lightweight and flexible structures in incompressible and viscous flows

The growing interest to examine the hydroelastic dynamics and stabilities of lightweight and flexible materials requires robust and accurate fluid–structure interaction(FSI)models. Classically, partitioned fluid and structure solvers are easier to implement compared to monolithic methods;however, partitioned FSI models are vulnerable to numerical("virtual added mass") instabilities for cases when the solid to fluid density ratio is low and if the flow is incompressible.As a partitioned method, the loosely hybrid coupled(LHC)method, which was introduced and validated in Young et al.(Acta Mech. Sin. 28:1030–1041, 2012), has been successfully used to efficiently and stably model lightweight and flexible structures. The LHC method achieves its numerical stability by, in addition to the viscous fluid forces, embedding potential flow approximations of the fluid induced forces to transform the partitioned FSI model into a semi-implicit scheme. The objective of this work is to derive and validate the numerical stability boundary of the LHC. The results show that the stability boundary of the LHC is much wider than traditional loosely coupled methods for a variety of numerical integration schemes. The results also show that inclusion of an estimate of the fluid inertial forces is the most critical to ensure the numerical stability when solving for fluid–structure interaction problems involving cases with a solid to fluid-added mass ratio less than one.     

Simulation of the two-fluid model on incompressible flow with Fractional Step method for both resolved and unresolved scale interfaces

In the present paper, the Fractional Step method usually used in single fluid flow is here extended and applied for the two-fluid model resolution using the finite volume discretization. The use of a projection method resolution instead of the usual pressure-correction method for multi-fluid flow, successfully avoids iteration processes. On the other hand, the main weakness of the two fluid model used for simulations of free surface flows, which is the numerical diffusion of the interface, is also solved by means of the conservative Level Set method (interface sharpening) (Strubelj et al., 2009). Moreover, the use of the algorithm proposed has allowed presenting different free-surface cases with or without Level Set implementation even under coarse meshes under a wide range of density ratios. Thus, the numerical results presented, numerically verified, experimentally validated and converged under high density ratios, shows the capability and reliability of this resolution method for both mixed and unmixed flows.     

A height–flux method for simulating free surface flows and interfaces

A new technique for the numerical simulation of the free surface flows is developed. This technique is based on the finite element method with penalty formulation, and a flux method for surface advection. The advection part which is completely independent of the momentum solver is based on subdividing the fluid domain into small subvolumes along one of the co-ordinate axis. The subvolumes are then used to find the height function which will later describe the free surface. The free surface of the fluid in each subvolume is approximated by a line segment and its slope is calculated using the volume of the fluid in the two neighbouring subvolumes. Later, the unidirectional volume flux from one subvolume to its neighbouring one is calculated using the conservation laws, and the new surface line segments are reconstructed. This technique, referred to as the Height–Flux Method (HFM) is implemented to simulate the temporal instability of a capillary jet. The results of the numerical simulation well predict the experimental data. It is also shown that the HFM is computationally more efficient than the techniques which use a kinematic boundary condition for the surface advection.     

A multi-dimensional time-marching method of characteristics for viscous and heat-conducting flows

Summary A numerical scheme is presented which employs the characteristic surfaces in space-time for solving Navier-Stokes equations for compressible fluid flow. We consider the general case of a three-dimensional flow, a simplification of which yields the equations of the two-dimensional case. Emphasis is put on the method itself. We apply it to simulate a laminar hypersonic flow around a circular cylinder of a five-components gas mixture of nitrogen and oxygen with thermally perfect constituents and at chemical nonequilibrium. First, the partial differential equations are transformed into a standard form with directional derivatives, enabling to attain the compatibility conditions, including the viscosity terms. These conditions are discretized by approximating their integrals along the corresponding characteristic surfaces. The result is an explicit time-marching numerical scheme. Using a grid fitted between the shock and the cylinder, and starting from roughly estimated initial conditions, a steady solution is searched. A comparison is made with the solution obtained under the assumption of a perfect gas. Received 6 April 1999; accepted for publication 13 May 1999     

Advances in the simulation of multi‐fluid flows with the particle finite element method. Application to bubble dynamics

In this work, we extend the Particle Finite Element Method (PFEM) to multi‐fluid flow problems with the aim of exploiting the fact that Lagrangian methods are specially well suited for tracking interfaces. We develop a numerical scheme able to deal with large jumps in the physical properties, included surface tension, and able to accurately represent all types of discontinuities in the flow variables. The scheme is based on decoupling the velocity and pressure variables through a pressure segregation method that takes into account the interface conditions. The interface is defined to be aligned with the moving mesh, so that it remains sharp along time, and pressure degrees of freedom are duplicated at the interface nodes to represent the discontinuity of this variable due to surface tension and variable viscosity. Furthermore, the mesh is refined in the vicinity of the interface to improve the accuracy and the efficiency of the computations. We apply the resulting scheme to the benchmark problem of a two‐dimensional bubble rising in a liquid column presented in Hysing et al. (International Journal for Numerical Methods in Fluids 2009; 60 : 1259–1288), and propose two breakup and coalescence problems to assess the ability of a multi‐fluid code to model topology changes. Copyright © 2010 John Wiley & Sons, Ltd.     

A two‐step Taylor‐characteristic‐based Galerkin method for incompressible flows and its application to flow over triangular cylinder with different incidence angles

An alternative characteristic‐based scheme, the two‐step Taylor‐characteristic‐based Galerkin method is developed based on the introduction of multi‐step temporal Taylor series expansion up to second order along the characteristic of the momentum equation. Contrary to the classical characteristic‐based split (CBS) method, the current characteristic‐based method does not require splitting the momentum equation, and segregate the calculation of the pressure from that of the velocity by using the momentum–pressure Poisson equation method. Some benchmark problems are used to examine the effectiveness of the proposed algorithm and to compare with the original CBS method, and the results show that the proposed method has preferable accuracy with less numerical dissipation. We further applied the method to the numerical simulation of flow around equilateral triangular cylinder with different incidence angles in free stream. In this numerical investigation, the flow simulations are carried out in the low Reynolds number range. Instantaneous streamlines around the cylinder are used as a means to visualize the wake region behind, and they clearly show the flow pattern around the cylinder in time. The influence of incidence angle on flow characteristic parameters such as Strouhal number, Drag and Lift coefficients are discussed quantitatively. Copyright © 2009 John Wiley & Sons, Ltd.     

Interaction of a deformable solid with two-phase flows: An Eulerian-based numerical model for fluid-structure interaction using the level contour reconstruction method

We describe the formulation of a method for fluid-structure interaction involving the coupling of moving and/or flexible solid structures with multiphase flows in the framework of the Level Contour Reconstruction Method. We present an Eulerian-based numerical procedure for tracking the motion and interaction of a liquid-gas interface with a fluid-solid interface in the Lagrangian frame together with the evaluation of the fluid transport equations coupled to those for the solid transport, namely the left Cauchy-Green strain tensor field, in the Eulerian frame. To prevent excessive dissipation due to the convective nature of the solid transport equation, a simple incompressibility constraint for the strain field is enforced. A single grid structure is used for both the fluid and solid phases which allows for a simple and natural coupling of the fluid and solid dynamics. Several benchmark tests are performed to show the accuracy of the numerical method and which demonstrate accurate results compared to several of those in the existing literature. In particular we show that surface tension effects including contact line dynamics on the deforming solid phase can be properly simulated. The three-phase interaction of a droplet impacting on a flexible cantilever is investigated in detail. The simulations follow the detailed motion of the droplet impact (and subsequent deformation, breakup, and fall trajectory) along with the motion of the deformable solid cantilever due to its own weight as well as due to the force of the droplet impact.     

Pressure‐based unified solver for gas‐liquid two‐phase flows where compressible and incompressible flows coexist

We propose a pressure‐based unified solver for gas‐liquid two‐phase flows where compressible and incompressible flows coexist. Unlike the original thermo–Cubic Interpolated Propagation Combined Unified Procedure (CIP‐CUP) method proposed by Himeno et al (Transactions of the Japan Society of Mechanical Engineers, Series B, 2003), we split the advection term of the governing equations into a conservation part and into the rest. The splitting of advection term has two advantages. One is the high degree of freedom in choosing discretization schemes such as central‐difference schemes, upwind schemes, and Total Variation Diminishing (TVD) schemes. The other is the ease of implementation on unstructured grids. The advantages enable the analyses of various flows such as turbulent and supersonic ones in actual complicated boundaries. Therefore, the solver is useful for practical analyses. The solver was validated on the following test cases: subsonic single‐phase flows, incompressible single‐phase turbulent flows, and incompressible gas‐liquid two‐phase flows. With unstructured grids, we obtained the equivalent results as the ones with structured grids. After the validations, subsonic jet impinging on a water pool was calculated and compared with experimental results. It was confirmed that the calculated results were consistent with the experimental ones.