3D two phase flows numerical simulations by SL‐VOF method

This paper describes the development of a semi‐Lagrangian computational method for simulating complex 3D two phase flows. The Navier–Stokes equations are solved separately in both fluids using a robust pseudo‐compressibility method able to deal with high density ratio. The interface tracking is achieved by the segment Lagrangian volume of fluid (SL‐VOF) method. The 2D SL‐VOF method using the concepts of VOF, piecewise linear interface calculation (PLIC) and Lagrangian advection of the interface is herein extended to 3D flows. Three different test cases of SL‐VOF 3D are presented for validation and comparison either with 2D flows or with other numerical methods. A good agreement is observed in each case. Copyright © 2004 John Wiley & Sons, Ltd.     

Validation of the S‐CLSVOF method with the density‐scaled balanced continuum surface force model in multiphase systems coupled with thermocapillary flows

Three numerical methods, namely, volume of fluid (VOF), simple coupled volume of fluid with level set (S‐CLSVOF), and S‐CLSVOF with the density‐scaled balanced continuum surface force (CSF) model, have been incorporated into OpenFOAM source code and were validated for their accuracy for three cases: (i) an isothermal static case, (ii) isothermal dynamic cases, and (iii) non‐isothermal dynamic cases with thermocapillary flow including dynamic interface deformation. Results have shown that the S‐CLSVOF method gives accurate results in the test cases with mild computation conditions, and the S‐CLSVOF technique with the density‐scaled balanced CSF model leads to accurate results in the cases of large interface deformations and large density and viscosity ratios. These show that these high accuracy methods would be appropriate to obtain accurate predictions in multiphase flow systems with thermocapillary flows. Copyright © 2016 John Wiley & Sons, Ltd.     

A simple and efficient error analysis for multi‐step solution of the Navier–Stokes equations

A simple error analysis is used within the context of segregated finite element solution scheme to solve incompressible fluid flow. An error indicator is defined based on the difference between a numerical solution on an original mesh and an approximated solution on a related mesh. This error indicator is based on satisfying the steady‐state momentum equations. The advantages of this error indicator are, simplicity of implementation (post‐processing step), ability to show regions of high and/or low error, and as the indicator approaches zero the solution approaches convergence. Two examples are chosen for solution; first, the lid‐driven cavity problem, followed by the solution of flow over a backward facing step. The solutions are compared to previously published data for validation purposes. It is shown that this rather simple error estimate, when used as a re‐meshing guide, can be very effective in obtaining accurate numerical solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

A further work on multi‐phase two‐fluid approach for compressible multi‐phase flows

This paper is to continue our previous work Niu (Int. J. Numer. Meth. Fluids 2001; 36 :351–371) on solving a two‐fluid model for compressible liquid–gas flows using the AUSMDV scheme. We first propose a pressure–velocity‐based diffusion term originally derived from AUSMDV scheme Wada and Liou (SIAM J. Sci. Comput. 1997; 18 (3):633—657) to enhance its robustness. The scheme can be applied to gas and liquid fluids universally. We then employ the stratified flow model Chang and Liou (J. Comput. Physics 2007; 225 :240–873) for spatial discretization. By defining the fluids in different regions and introducing inter‐phasic force on cell boundary, the stratified flow model allows the conservation laws to be applied on each phase, and therefore, it is able to capture fluid discontinuities, such as the fluid interfaces and shock waves, accurately. Several benchmark tests are studied, including the Ransom's Faucet problem, 1D air–water shock tube problems, 2D shock‐water column and 2D shock‐bubble interaction problems. The results indicate that the incorporation of the new dissipation into AUSM⁺‐up scheme and the stratified flow model is simple, accurate and robust enough for the compressible multi‐phase flows. Copyright © 2008 John Wiley & Sons, Ltd.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

A locally DIV‐free projection scheme for incompressible flows based on non‐conforming finite elements

We present a projection scheme whose end‐of‐step velocity is locally pointwise divergence free, using a continuous ?₁ approximation for the velocity in the momentum equation, a first‐order Crouzeix–Raviart approximation at the projection step, and a ?₀ approximation for the pressure in both steps. The analysis of the scheme is done only for grids that guarantee the existence of a divergence free conforming ?₁ interpolant for the velocity. Optimal estimates for the velocity error in L²‐ and H¹‐norms are deduced. The numerical results demonstrate that these estimates should also hold on grids on which the continuous ?₁ approximation for the velocity locks. Since the end‐of‐step velocity is locally solenoidal, the scheme is recommendable for problems requiring good mass conservation. Copyright © 2005 John Wiley & Sons, Ltd.     

An unsteady single‐phase level set method for viscous free surface flows

The single‐phase level set method for unsteady viscous free surface flows is presented. In contrast to the standard level set method for incompressible flows, the single‐phase level set method is concerned with the solution of the flow field in the water (or the denser) phase only. Some of the advantages of such an approach are that the interface remains sharp, the computation is performed within a fluid with uniform properties and that only minor computations are needed in the air. The location of the interface is determined using a signed distance function, and appropriate interpolations at the fluid/fluid interface are used to enforce the jump conditions. A reinitialization procedure has been developed for non‐orthogonal grids with large aspect ratios. A convective extension is used to obtain the velocities at previous time steps for the grid points in air, which allows a good estimation of the total derivatives. The method was applied to three unsteady tests: a plane progressive wave, sloshing in a two‐dimensional tank, and the wave diffraction problem for a surface ship, and the results compared against analytical solutions or experimental data. The method can in principle be applied to any problem in which the standard level set method works, as long as the stress on the second phase can be specified (or neglected) and no bubbles appear in the flow during the computation. Copyright © 2006 John Wiley & Sons, Ltd.     

The capturing of free surfaces in incompressible multi‐fluid flows

By treating it as a contact discontinuity in the density field, a free surface between two immiscible fluids can be automatically ‘captured’ by the enforcement of conservation laws. A surface‐capturing method of this kind requires no special tracking or fitting treatment for the free surface, thereby offering the advantage of algorithm simplicity over the surface‐tracking or the surface‐fitting method. A surface‐capturing method based on a new multi‐fluid incompressible Navier–Stokes formulation is developed. It is applied to a variety of free‐surface flows, including the Rayleigh–Taylor instability problem, the ship waves around a Wigley hull and a model bubble‐rising problem to demonstrate the validity and versatility of the present method. Copyright © 2000 John Wiley & Sons, Ltd.     

Colocated schemes for the incompressible Navier–Stokes equations on non‐smooth grids for two‐dimensional problems

The accuracy of colocated finite volume schemes for the incompressible Navier–Stokes equations on non‐smooth curvilinear grids is investigated. A frequently used scheme is found to be quite inaccurate on non‐smooth grids. In an attempt to improve the accuracy on such grids, three other schemes are described and tested. Two of these are found to give satisfactory results. Copyright © 2000 John Wiley & Sons, Ltd.     

An unfitted interior penalty discontinuous Galerkin method for incompressible Navier–Stokes two‐phase flow

A discontinuous Galerkin method for the solution of the immiscible and incompressible two‐phase flow problem based on the nonsymmetric interior penalty method is presented. Therefore, the incompressible Navier–Stokes equation is solved for a domain decomposed into two subdomains with different values of viscosity and density as well as a singular surface tension force. On the basis of a piecewise linear approximation of the interface, meshes for both phases are cut out of a structured mesh. The discontinuous finite elements are defined on the resulting Cartesian cut‐cell mesh and may therefore approximate the discontinuities of the pressure and the velocity derivatives across the interface with high accuracy. As the mesh resolves the interface, regularization of the density and viscosity jumps across the interface is not required. This preserves the local conservation property of the velocity field even in the vicinity of the interface and constitutes a significant advantage compared with standard methods that require regularization of these discontinuities and cannot represent the jumps and kinks in pressure and velocity. A powerful subtessellation algorithm is incorporated to allow the usage of standard time integrators (such as Crank–Nicholson) on the time‐dependent mesh. The presented discretization is applicable to both the two‐dimensional and three‐dimensional cases. The performance of our approach is demonstrated by application to a two‐dimensional benchmark problem, allowing for a thorough comparison with other numerical methods. Copyright © 2012 John Wiley & Sons, Ltd.     

Fourth‐ and tenth‐order compact finite difference solutions of perturbed circular vortex flows

In this study, high‐order compact finite difference calculations are reported for 2D unsteady incompressible circular vortex flow in primitive variable formulation. The fourth‐order Runge–Kutta temporal discretization is used together with fourth‐ or tenth‐order compact spatial discretization. Dependent on the perturbation initially imposed, the solutions display a tripole, triangular or square vortex. The comparison of the predictions with the detailed spectral calculations of Kloosterziel and Carnevale (J. Fluid Mech. 1999; 388 :217–257) shows that the vorticity fields are very well captured. The spectral resolution of the present method was quantified from the decomposition of the vorticity distribution in its azimuthal components and compared with reported spectral results. Using identical grid resolution to the reference results yields negligible differences in the main features of the flow. The perturbation amplitude and its first harmonic are virtually identical to the reference results for both fourth‐ or tenth‐order spatial discretization, as theoretically expected but seldom a posteriori verified. The differences between the two spatial discretizations appear only for coarser grids, favouring the tenth‐order discretization. Copyright © 2005 John Wiley & Sons, Ltd.     

SIMPLE‐type preconditioners for cell‐centered,colocated finite volume discretization of incompressible Reynolds‐averaged Navier–Stokes equations

This paper contains a comparison of four SIMPLE‐type methods used as solver and as preconditioner for the iterative solution of the (Reynolds‐averaged) Navier–Stokes equations, discretized with a finite volume method for cell‐centered, colocated variables on unstructured grids. A matrix‐free implementation is presented, and special attention is given to the treatment of the stabilization matrix to maintain a compact stencil suitable for unstructured grids. We find SIMPLER preconditioning to be robust and efficient for academic test cases and industrial test cases. Compared with the classical SIMPLE solver, SIMPLER preconditioning reduces the number of nonlinear iterations by a factor 5–20 and the CPU time by a factor 2–5 depending on the case. The flow around a ship hull at Reynolds number 2E9, for example, on a grid with cell aspect ratio up to 1:1E6, can be computed in 3 instead of 15 h.Copyright © 2012 John Wiley & Sons, Ltd.     

A compact finite difference method on staggered grid for Navier–Stokes flows

Compact finite difference methods feature high‐order accuracy with smaller stencils and easier application of boundary conditions, and have been employed as an alternative to spectral methods in direct numerical simulation and large eddy simulation of turbulence. The underpinning idea of the method is to cancel lower‐order errors by treating spatial Taylor expansions implicitly. Recently, some attention has been paid to conservative compact finite volume methods on staggered grid, but there is a concern about the order of accuracy after replacing cell surface integrals by average values calculated at centres of cell surfaces. Here we introduce a high‐order compact finite difference method on staggered grid, without taking integration by parts. The method is implemented and assessed for an incompressible shear‐driven cavity flow at Re = 10³, a temporally periodic flow at Re = 10⁴, and a spatially periodic flow at Re = 10⁴. The results demonstrate the success of the method. Copyright © 2006 John Wiley & Sons, Ltd.     

A high‐order mass‐lumping procedure for B‐spline collocation method with application to incompressible flow simulations

This paper presents new developments of the staggered spline collocation method for cost‐effective solution to the incompressible Navier–Stokes equations. Maximal decoupling of the velocity and the pressure is obtained by using the fractional step method of Gresho and Chan, allowing the solution to sparse elliptic problems only. In order to preserve the high‐accuracy of the B‐spline method, this fractional step scheme is used in association with a sparse approximation to the inverse of the consistent mass matrix. Such an approximation is constructed from local spline interpolation method, and represents a high‐order generalization of the mass‐lumping technique of the finite‐element method. A numerical investigation of the accuracy and the computational efficiency of the resulting semi‐consistent spline collocation schemes is presented. These schemes generate a stable and accurate unsteady Navier–Stokes solver, as assessed by benchmark computations. Copyright © 2003 John Wiley & Sons, Ltd.