A robust multi‐grid pressure‐based algorithm for multi‐fluid flow at all speeds

This paper reports on the implementation and testing, within a full non‐linear multi‐grid environment, of a new pressure‐based algorithm for the prediction of multi‐fluid flow at all speeds. The algorithm is part of the mass conservation‐based algorithms (MCBA) group in which the pressure correction equation is derived from overall mass conservation. The performance of the new method is assessed by solving a series of two‐dimensional two‐fluid flow test problems varying from turbulent low Mach number to supersonic flows, and from very low to high fluid density ratios. Solutions are generated for several grid sizes using the single grid (SG), the prolongation grid (PG), and the full non‐linear multi‐grid (FMG) methods. The main outcomes of this study are: (i) a clear demonstration of the ability of the FMG method to tackle the added non‐linearity of multi‐fluid flows, which is manifested through the performance jump observed when using the non‐linear multi‐grid approach as compared to the SG and PG methods; (ii) the extension of the FMG method to predict turbulent multi‐fluid flows at all speeds. The convergence history plots and CPU‐times presented indicate that the FMG method is far more efficient than the PG method and accelerates the convergence rate over the SG method, for the problems solved and the grids used, by a factor reaching a value as high as 15. Copyright © 2003 John Wiley & Sons, Ltd.     

A stress‐based least‐squares finite‐element model for incompressible Navier–Stokes equations

In this paper we present a stress‐based least‐squares finite‐element formulation for the solution of the Navier–Stokes equations governing flows of viscous incompressible fluids. Stress components are introduced as independent variables to make the system first order. Continuity equation becomes an algebraic equation and is eliminated from the system with suitable modifications. The h and p convergence are verified using the exact solution of Kovasznay flow. Steady flow past a large circular cylinder in a channel is solved to test mass conservation. Transient flow over a backward‐facing step problem is solved on several meshes. Results are compared with that obtained using vorticity‐based first‐order formulation for both benchmark problems. Copyright © 2007 John Wiley & Sons, Ltd.     

An hp‐adaptive unstructured finite volume solver for compressible flows

The idea of hp‐adaptation, which has originally been developed for compact schemes (such as finite element methods), suggests an adaptation scheme using a mixture of mesh refinement and order enrichment based on the smoothness of the solution to obtain an accurate solution efficiently. In this paper, we develop an hp‐adaptation framework for unstructured finite volume methods using residual‐based and adjoint‐based error indicators. For the residual‐based error indicator, we use a higher‐order discrete operator to estimate the truncation error, whereas this estimate is weighted by the solution of the discrete adjoint problem for an output of interest to form the adaptation indicator for adjoint‐based adaptations. We perform our adaptation by local subdivision of cells with nonconforming interfaces allowed and local reconstruction of higher‐order polynomials for solution approximations. We present our results for two‐dimensional compressible flow problems including subsonic inviscid, transonic inviscid, and subsonic laminar flow around the NACA 0012 airfoil and also turbulent flow over a flat plate. Our numerical results suggest the efficiency and accuracy advantages of adjoint‐based hp‐adaptations over uniform refinement and also over residual‐based adaptation for flows with and without singularities.     

A pressure‐based unstructured grid method for all‐speed flows

An all‐speed algorithm based on the SIMPLE pressure‐correction scheme and the ‘retarded‐density’ approach has been formulated and implemented within an unstructured grid, finite volume (FV) scheme for both incompressible and compressible flows, the latter involving interaction of shock waves. The collocated storage arrangement for all variables is adopted, and the checkerboard oscillations are eliminated by using a pressure‐weighted interpolation method, similar to that of Rhie and Chow [Numerical study of the turbulent flow past an airfoil with trailing edge separation. AIAA Journal 1983; 21 : 1525]. The solution accuracy is greatly enhanced when a higher‐order convection scheme combined with adaptive mesh refinement (AMR) are used. Copyright © 2000 John Wiley & Sons, Ltd.     

Three‐dimensional modeling of non‐hydrostatic free‐surface flows on unstructured grids

A three‐dimensional numerical model is presented for the simulation of unsteady non‐hydrostatic shallow water flows on unstructured grids using the finite volume method. The free surface variations are modeled by a characteristics‐based scheme, which simulates sub‐critical and super‐critical flows. Three‐dimensional velocity components are considered in a collocated arrangement with a σ‐coordinate system. A special treatment of the pressure term is developed to avoid the water surface oscillations. Convective and diffusive terms are approximated explicitly, and an implicit discretization is used for the pressure term to ensure exact mass conservation. The unstructured grid in the horizontal direction and the σ coordinate in the vertical direction facilitate the use of the model in complicated geometries. Solution of the non‐hydrostatic equations enables the model to simulate short‐period waves and vertically circulating flows. Copyright © 2015 John Wiley & Sons, Ltd.     

Incompressible multiphase flow and encapsulation simulations using the moment‐of‐fluid method

A moment‐of‐fluid method is presented for computing solutions to incompressible multiphase flows in which the number of materials can be greater than two. In this work, the multimaterial moment‐of‐fluid interface representation technique is applied to simulating surface tension effects at points where three materials meet. The advection terms are solved using a directionally split cell integrated semi‐Lagrangian algorithm, and the projection method is used to evaluate the pressure gradient force term. The underlying computational grid is a dynamic block‐structured adaptive grid. The new method is applied to multiphase problems illustrating contact‐line dynamics, triple junctions, and encapsulation in order to demonstrate its capabilities. Examples are given in two‐dimensional, three‐dimensional axisymmetric (R–Z), and three‐dimensional (X–Y–Z) coordinate systems. Copyright © 2015 John Wiley & Sons, Ltd.     

A hybrid pressure‐based solver for nonideal single‐phase fluid flows at all speeds

This paper describes the implementation of a numerical solver that is capable of simulating compressible flows of nonideal single‐phase fluids. The proposed method can be applied to arbitrary equations of state and is suitable for all Mach numbers. The pressure‐based solver uses the operator‐splitting technique and is based on the PISO/SIMPLE algorithm: the density, velocity, and temperature fields are predicted by solving the linearized versions of the balance equations using the convective fluxes from the previous iteration or time step. The overall mass continuity is ensured by solving the pressure equation derived from the continuity equation, the momentum equation, and the equation of state. Nonphysical oscillations of the numerical solution near discontinuities are damped using the Kurganov‐Tadmor/Kurganov‐Noelle‐Petrova (KT/KNP) scheme for convective fluxes. The solver was validated using different test cases, where analytical and/or numerical solutions are present or can be derived: (1) A convergent‐divergent nozzle with three different operating conditions; (2) the Riemann problem for the Peng‐Robinson equation of state; (3) the Riemann problem for the covolume equation of state; (4) the development of a laminar velocity profile in a circular pipe (also known as Poiseuille flow); (5) a laminar flow over a circular cylinder; (6) a subsonic flow over a backward‐facing step at low Reynolds numbers; (7) a transonic flow over the RAE 2822 airfoil; and (8) a supersonic flow around a blunt cylinder‐flare model. The spatial approximation order of the scheme is second order. The mesh convergence of the numerical solution was achieved for all cases. The accuracy order for highly compressible flows with discontinuities is close to first order and, for incompressible viscous flows, it is close to second order. The proposed solver is named rhoPimpleCentralFoam and is implemented in the open‐source CFD library OpenFOAM^®. For high speed flows, it shows a similar behavior as the KT/KNP schemes (implemented as rhoCentralFoam‐solver, Int. J. Numer. Meth. Fluids 2010), and for flows with small Mach numbers, it behaves like solvers that are based on the PISO/SIMPLE algorithm.     

Large eddy simulations of wall‐bounded flows using a simplified immersed boundary method and high‐order compact schemes

This paper presents a solution algorithm based on an immersed boundary (IB) method that can be easily implemented in high‐order codes for incompressible flows. The time integration is performed using a predictor‐corrector approach, and the projection method is used for pressure‐velocity coupling. Spatial discretization is based on compact difference schemes and is performed on half‐staggered meshes. A basic algorithm for body‐fitted meshes using the aforementioned solution method was developed by A. Tyliszczak (see article “A high‐order compact difference algorithm for half‐staggered grids for laminar and turbulent incompressible flows” in Journal of Computational Physics) and proved to be very accurate. In this paper, the formulated algorithm is adapted for use with the IB method in the framework of large eddy simulations. The IB method is implemented using its simplified variant without the interpolation (stepwise approach). The computations are performed for a laminar flow around a 2D cylinder, a turbulent flow in a channel with a wavy wall, and around a sphere. Comparisons with literature data confirm that the proposed method can be successfully applied for complex flow problems. The results are verified using the classical approach with body‐fitted meshes and show very good agreement both in laminar and turbulent regimes. The mean (velocity and turbulent kinetic energy profiles and drag coefficients) and time‐dependent (Strouhal number based on the drag coefficient) quantities are analyzed, and they agree well with reference solutions. Two subfilter models are compared, ie, the model of Vreman (see article “An eddy‐viscosity subgrid‐scale model for turbulent shear flow: algebraic theory and applications” in Physics and Fluids) and σ model (Nicoud et al, see article “Using singular values to build a subgrid‐scale model for large eddy simulations” in Physics and Fluids). The tests did not reveal evident advantages of any of these models, and from the point of view of solution accuracy, the quality of the computational meshes turned out to be much more important than the subfilter modeling.     

Effects of element order and interface reconstruction in FEM/volume‐of‐fluid incompressible flow simulation

In previous studies, the moment‐of‐fluid interface reconstruction method showed dramatic accuracy improvements in static and pure advection tests over existing methods, but this did not translate into an equivalent improvement in volume‐tracked multimaterial incompressible flow simulation using low‐order finite elements. In this work, the combined effects of the spatial discretization and interface reconstruction in flow simulation are examined. The mixed finite element pairs, Q₁Q₀ (with pressure stabilization) and Q₂P _{? 1} are compared. Material order‐dependent and material order‐independent first and second‐order accurate interface reconstruction methods are used. The Q₂P _{? 1} elements show significant improvements in computed flow solution accuracy for single material flows but show reduced convergence using element‐average piecewise constant density and viscosity in volume‐tracked simulations. In general, a refined Q₁Q₀ grid, with better material interface resolution, provided an accuracy similar to the Q₂P _{? 1} element grid with a comparable number of degrees of freedom. Moment‐of‐fluid shows more benefit from the higher‐order accurate flow simulation than the LVIRA, Youngs', and power diagram interface reconstruction methods, especially on unstructured grids, but does not recover the dramatic accuracy improvements it has shown in advection tests. Published 2012. This article is a US Government work and is in the public domain in the USA.     

A low‐dimensional spectral approach for transient axisymmetric free‐surface flow inside thin cavities of arbitrary shape

A low‐dimensional spectral method is used to solve the transient axisymmetric free surface flow inside thin cavities of arbitrary shape. The flow field is obtained on the basis of the lubrication equations, which are expanded in terms of orthonormal functions over the cavity gap. The formulation accounts for nonlinearities stemming from inertia and front location. The work is of close relevance to the filling stage during die casting, and injection molding, or the flow inside annular (extrusion) dies. Both flows under an imposed flow rate, and an imposed pressure at the cavity entrance are examined. The influence of inertia, aspect ratio, gravity, and wall geometry on the evolution of the front, flow rate, and pressure is assessed particularly in the early stage of flow, when a temporal behavior of the ‘boundary‐layer’ type develops. The multiple‐scale method is applied to obtain an approximate solution at small Reynolds number, Re. Comparison with the exact (numerical) solution indicates a wide range of validity for the multiple‐scale approach, including the moderately small Re range. Copyright © 2002 John Wiley & Sons, Ltd.     

A multiblock operator‐splitting algorithm for unsteady flows at all speeds in complex geometries

This paper describes a non‐iterative operator‐splitting algorithm for computing all‐speed flows in complex geometries. A pressure‐based algorithm is adopted as the base, in which pressure, instead of density, is a primary variable, thus allowing for a unified formulation for all Mach numbers. The focus is on adapting the method for (a) flows at all speeds, and (b) multiblock, non‐orthogonal, body‐fitted grids for very complex geometries. Key features of the formulation include special treatment of mass fluxes at control volume interfaces to avoid pressure–velocity decoupling for incompressible (low Mach number limit) flows and to provide robust pressure–velocity–density coupling for compressible (high‐speed) flows. The method is shown to be robust for all Mach number regimes for both steady and unsteady flows; it is found to be stable for CFL numbers of order ten, allowing large time steps to be taken for steady flows. Enhancements to the method which allow for stable solutions to be obtained on non‐orthogonal grids are also discussed. The method is found to be very reliable even in complex engineering applications such as unsteady rotor–stator interactions in turbulent, all‐speed turbomachinery flows. Copyright © 2004 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.     

Time‐dependent turbulent cavitating flow computations with interfacial transport and filter‐based models

Turbulent cavitating flow computations need to address both cavitation and turbulence modelling issues. A recently developed interfacial dynamics‐based cavitation model (IDCM) incorporates the interfacial transport into the computational modelling of cavitation dynamics. For time‐dependent flows, it is known that the engineering turbulence closure such as the original k–ε model often over‐predicts the eddy viscosity values reducing the unsteadiness. A recently proposed filter‐based modification has shown that it can effectively modulate the eddy viscosity, rendering better simulation capabilities for time‐dependent flow computations in term of the unsteady characteristics. In the present study, the IDCM along with the filter‐based k–ε turbulence model is adopted to simulate 2‐D cavitating flows over the Clark‐Y airfoil. The chord Reynolds number is Re=7.0 × 10⁵. Two angles‐of‐attack of 5 and 8° associated with several cavitation numbers covering different flow regimes are conducted. The simulation results are assessed with the experimental data including lift, drag and velocity profiles. The interplay between cavitation and turbulence models reveals substantial differences in time‐dependent flow results even though the time‐averaged characteristics are similar. Copyright © 2005 John Wiley & Sons, Ltd.     

Multi‐material incompressible flow simulation using the moment‐of‐fluid method

This paper compares the numerical performance of the moment‐of‐fluid (MOF) interface reconstruction technique with Youngs, LVIRA, power diagram (PD), and Swartz interface reconstruction techniques in the context of a volume‐of‐fluid (VOF) based finite element projection method for the numerical simulation of variable‐density incompressible viscous flows. In pure advection tests with multiple materials MOF shows dramatic improvements in accuracy compared with the other methods. In incompressible flows where density differences determine the flow evolution, all the methods perform similarly for two material flows on structured grids. On unstructured grids, the second‐order MOF, LVIRA, and Swartz methods perform similarly and show improvement over the first‐order Youngs' and PD methods. For flow simulations with more than two materials, MOF shows increased accuracy in interface positions on coarse meshes. In most cases, the convergence and accuracy of the computed flow solution was not strongly affected by interface reconstruction method. Published in 2009 by John Wiley & Sons, Ltd.     

Augmented mixed finite element methods for a vorticity‐based velocity–pressure–stress formulation of the Stokes problem in 2D

In this paper, we consider an augmented velocity–pressure–stress formulation of the 2D Stokes problem, in which the stress is defined in terms of the vorticity and the pressure, and then we introduce and analyze stable mixed finite element methods to solve the associated Galerkin scheme. In this way, we further extend similar procedures applied recently to linear elasticity and to other mixed formulations for incompressible fluid flows. Indeed, our approach is based on the introduction of the Galerkin least‐squares‐type terms arising from the corresponding constitutive and equilibrium equations, and from the Dirichlet boundary condition for the velocity, all of them multiplied by stabilization parameters. Then, we show that these parameters can be suitably chosen so that the resulting operator equation induces a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well posed for any choice of finite element subspaces. In particular, we can use continuous piecewise linear velocities, piecewise constant pressures, and rotated Raviart–Thomas elements for the stresses. Next, we derive reliable and efficient residual‐based a posteriori error estimators for the augmented mixed finite element schemes. In addition, several numerical experiments illustrating the performance of the augmented mixed finite element methods, confirming the properties of the a posteriori estimators, and showing the behavior of the associated adaptive algorithms are reported. The present work should be considered as a first step aiming finally to derive augmented mixed finite element methods for vorticity‐based formulations of the 3D Stokes problem. Copyright © 2010 John Wiley & Sons, Ltd.     

A new positive‐definite regularization of incompressible Navier–Stokes equations discretized with Q1/P0 finite element

A new regularization method is proposed for the Galerkin approximation of the incompressible Navier–Stokes equations with Q1/P0 element, by newly introducing a square‐type linear form into the variational divergence‐free constraint regularized with the global pressure jump (GPJ) method. The addition of the square‐type linear form is intended to eliminate the hydrostatic pressure mode appearing in confined flows, and to make the discretized matrix positive definite and then non‐singular without the pressure pegging trick. Effects of the free parameters for the regularization on the solutions are numerically examined with a 2‐D driven cavity flow problem. Furthermore, the convergences in the conjugate gradient iteration for the solution of the pressure Poisson equation are compared among the mixed method, the GPJ method and the present method for both leaky and non‐leaky 3‐D driven cavity flows. Finally, the non‐leaky 3‐D cavity flows at different Re numbers are solved to compare with the literature data and to demonstrate the accuracy of the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

Interfacial dynamics‐based modelling of turbulent cavitating flows,Part‐2: Time‐dependent computations

The interfacial dynamics‐based cavitation model, developed in Part‐1, is further employed for unsteady flow computations. The pressure‐based operator‐splitting algorithm (PISO) is extended to handle the time‐dependent cavitating flows with particular focus on the coupling of the cavitation and turbulence models, and the large density ratio associated with cavitation. Furthermore, the compressibility effect is important for unsteady cavitating flows because in a water–vapour mixture, depending on the composition, the speed of sound inside the cavity can vary by an order of magnitude. The implications of the issue of the speed of the sound are assessed with alternative modelling approaches. Depending on the geometric confinement of the nozzle, compressibility model and cavitation numbers, either auto‐oscillation or quasi‐steady behaviour is observed. The adverse pressure gradient in the closure region is stronger at the maximum cavity size. One can also observe that the mass transfer process contributes to the cavitation dynamics. Compared to the steady flow computations, the velocity and vapour volume fraction distributions within the cavity are noticeably improved with time‐dependent computations. Copyright © 2004 John Wiley & Sons, Ltd.     

Compressive advection and multi‐component methods for interface‐capturing

This paper develops methods for interface‐capturing in multiphase flows. The main novelties of these methods are as follows: (a) multi‐component modelling that embeds interface structures into the continuity equation; (b) a new family of triangle/tetrahedron finite elements, in particular, the P₁DG‐P₂(linear discontinuous between elements velocity and quadratic continuous pressure); (c) an interface‐capturing scheme based on compressive control volume advection methods and high‐order finite element interpolation methods; (d) a time stepping method that allows use of relatively large time step sizes; and (e) application of anisotropic mesh adaptivity to focus the numerical resolution around the interfaces and other areas of important dynamics. This modelling approach is applied to a series of pure advection problems with interfaces as well as to the simulation of the standard computational fluid dynamics benchmark test cases of a collapsing water column under gravitational forces (in two and three dimensions) and sloshing water in a tank. Two more test cases are undertaken in order to demonstrate the many‐material and compressibility modelling capabilities of the approach. Numerical simulations are performed on coarse unstructured meshes to demonstrate the potential of the methods described here to capture complex dynamics in multiphase flows. Copyright © 2015 John Wiley & Sons, Ltd.