A volume-tracking method for incompressible multifluid flows with large density variations

A numerical technique (FGVT) for solving the time-dependent incompressible Navier–Stokes equations in fluid flows with large density variations is presented for staggered grids. Mass conservation is based on a volume tracking method and incorporates a piecewise-linear interface reconstruction on a grid twice as fine as the velocity–pressure grid. It also uses a special flux-corrected transport algorithm for momentum advection, a multigrid algorithm for solving a pressure-correction equation and a surface tension algorithm that is robust and stable. In principle, the method conserves both mass and momentum exactly, and maintains extremely sharp fluid interfaces. Applications of the numerical method to prediction of two-dimensional bubble rise in an inclined channel and a bubble bursting through an interface are presented. © 1998 John Wiley & Sons, Ltd.     

Hybrid finite element/volume method for 3D incompressible flows with heat transfer

An implicit hybrid finite element (FE)/volume solver has been extended to incompressible flows coupled with the energy equation. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centred finite volume (FV) method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centres and the auxiliary variable at vertices. The Generalized Minimal Residual (GMRES) matrix-free strategy is adapted to solve the governing equations in both FE and FV methods. The presented 2D and 3D numerical examples show the robustness and accuracy of the numerical method.     

Adaptive hybrid (prismatic–tetrahedral) grids for incompressible flows

Hybrid grids consisting of prisms and tetrahedra are employed for the solution of the 3-D Navier–Stokes equations of incompressible flow. A pressure correction scheme is employed with a finite volume–finite element spatial discretization. The traditional staggered grid formulation has been substituted with a collocated mesh approach which uses fourth-order artificial dissipation. The hybrid grid is refined adaptively in local regions of appreciable flow variations. The scheme operations are performed on an edge-wise basis which unifies treatment of both types of grid elements. The adaptive method is employed for incompressible flows in both single and multiply-connected domains. © 1998 John Wiley & Sons, Ltd.     

A divergence-free semi-implicit finite volume scheme for ideal,viscous, and resistive magnetohydrodynamics

In this paper, we present a novel pressure-based semi-implicit finite volume solver for the equations of compressible ideal, viscous, and resistive magnetohydrodynamics (MHD). The new method is conservative for mass, momentum, and total energy, and in multiple space dimensions, it is constructed in such a way as to respect the divergence-free condition of the magnetic field exactly, also in the presence of resistive effects. This is possible via the use of multidimensional Riemann solvers on an appropriately staggered grid for the time evolution of the magnetic field and a double curl formulation of the resistive terms. The new semi-implicit method for the MHD equations proposed here discretizes the nonlinear convective terms as well as the time evolution of the magnetic field explicitly, whereas all terms related to the pressure in the momentum equation and the total energy equation are discretized implicitly, making again use of a properly staggered grid for pressure and velocity. Inserting the discrete momentum equation into the discrete energy equation then yields a mildly nonlinear symmetric and positive definite algebraic system for the pressure as the only unknown, which can be efficiently solved with the (nested) Newton method of Casulli et al. The pressure system becomes linear when the specific internal energy is a linear function of the pressure. The time step of the scheme is restricted by a CFL condition based only on the fluid velocity and the Alfvén wave speed and is not based on the speed of the magnetosonic waves. Being a semi-implicit pressure-based scheme, our new method is therefore particularly well suited for low Mach number flows and for the incompressible limit of the MHD equations, for which it is well known that explicit density-based Godunov-type finite volume solvers become increasingly inefficient and inaccurate because of the more and more stringent CFL condition and the wrong scaling of the numerical viscosity in the incompressible limit. We show a relevant MHD test problem in the low Mach number regime where the new semi-implicit algorithm is a factor of 50 faster than a traditional explicit finite volume method, which is a very significant gain in terms of computational efficiency. However, our numerical results confirm that our new method performs well also for classical MHD test cases with strong shocks. In this sense, our new scheme is a true all Mach number flow solver.     

Solution of the incompressible Navier–Stokes equations by the method of lines

The numerical method of lines (NUMOL) is a numerical technique used to solve efficiently partial differential equations. In this paper, the NUMOL is applied to the solution of the two‐dimensional unsteady Navier–Stokes equations for incompressible laminar flows in Cartesian coordinates. The Navier–Stokes equations are first discretized (in space) on a staggered grid as in the Marker and Cell scheme. The discretized Navier–Stokes equations form an index 2 system of differential algebraic equations, which are afterwards reduced to a system of ordinary differential equations (ODEs), using the discretized form of the continuity equation. The pressure field is computed solving a discrete pressure Poisson equation. Finally, the resulting ODEs are solved using the backward differentiation formulas. The proposed method is illustrated with Dirichlet boundary conditions through applications to the driven cavity flow and to the backward facing step flow. Copyright © 2015 John Wiley & Sons, Ltd.     

A PRESSURE-CORRECTION METHOD FOR THE SOLUTION OF INCOMPRESSIBLE VISCOUS FLOWS ON UNSTRUCTURED GRIDS

An unstructured grid, finite volume method is presented for the solution of two-dimensional viscous, incompressible flow. The method is based on the pressure-correction concept implemented on a semi-staggered grid. The computational procedure can handle cells of arbitrary shape, although solutions presented herein have been obtained only with meshes of triangular and quadrilateral cells. The discretization of the momentum equations is effected on dual cells surrounding the vertices of primary cells, while the pressure-correction equation applies to the primary-cell centroids and represents the conservation of mass across the primary cells. A special interpolation scheme s used to suppress pressure and velocity oscillations in cases where the semi-staggered arrangement does not ensure a sufficiently strong coupling between pressure and velocity to avoid such oscillations. Computational results presented for several viscous flows are shown to be in good agreement with analytical and experimental data reported in the open literature.     

The (9,5) HOC formulation for the transient Navier–Stokes equations in primitive variable

We recently proposed an improved (9,5) higher order compact (HOC) scheme for the unsteady two‐dimensional (2‐D) convection–diffusion equations. Because of using only five points at the current time level in the discretization procedure, the scheme was seen to be computationally more efficient than its predecessors. It was also seen to capture very accurately the solution of the unsteady 2‐D Navier–Stokes (N–S) equations for incompressible viscous flows in the stream function–vorticity (ψ – ω) formulation. In this paper, we extend the scope of the scheme for solving the unsteady incompressible N–S equations based on primitive variable formulation on a collocated grid. The parabolic momentum equations are solved for the velocity field by a time‐marching strategy and the pressure is obtained by discretizing the elliptic pressure Poisson equation by the steady‐state form of the (9,5) scheme with the Neumann boundary conditions. In particular, for pressure, we adopt a strategy on the collocated grid in conjunction with ideas borrowed from the staggered grid approach in finite volume. We first apply this extension to a problem having analytical solution and then to the famous lid‐driven square cavity problem. We also apply our formulation to the backward‐facing step problem to see how the method performs for external flow problems. The results are presented and are compared with established numerical results. This new approach is seen to produce excellent comparison in all the cases. Copyright © 2007 John Wiley & Sons, Ltd.     

A hybrid finite volume/finite element method for incompressible generalized Newtonian fluid flows on unstructured triangular meshes

This paper presents a hybrid finite volume/finite element method for the incompressible generalized Newtonian fluid flow （Power-Law model）. The collocated （i.e. non-staggered） arrangement of variables is used on the unstructured triangular grids, and a fractional step projection method is applied for the velocity-pressure coupling. The cell-centered finite volume method is employed to discretize the momentum equation and the vertex-based finite element for the pressure Poisson equation. The momentum interpolation method is used to suppress unphysical pressure wiggles. Numerical experiments demonstrate that the current hybrid scheme has second order accuracy in both space and time. Results on flows in the lid-driven cavity and between parallel walls for Newtonian and Power-Law models are also in good agreement with the published solutions.     

An implementation of the Spalart–Allmaras DES model in an implicit unstructured hybrid finite volume/element solver for incompressible turbulent flow

In this paper, we describe an implicit hybrid finite volume (FV)/element (FE) incompressible Navier–Stokes solver for turbulent flows based on the Spalart–Allmaras detached eddy simulation (SA‐DES). The hybrid FV/FE solver is based on the segregated pressure correction or projection method. The intermediate velocity field is first obtained by solving the original momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centers and the auxiliary variable at vertices, making the current solver a staggered‐mesh scheme. The SA‐DES turbulence equation is solved after the velocity and the pressure fields have been updated at the end of each time step. The same matrix‐free FV method as the one used for momentum equations is used to solve the turbulence equation. The turbulence equation provides the eddy viscosity, which is added to the molecular viscosity when solving the momentum equation. In our implementation, we focus on the accuracy, efficiency and robustness of the SA‐DES model in a hybrid flow solver. This paper will address important implementation issues for high‐Reynolds number flows where highly stretched elements are typically used. In addition, some aspects of implementing the SA‐DES model will be described to ensure the robustness of the turbulence model. Several numerical examples including a turbulent flow past a flat plate and a high‐Reynolds number flow around a high angle‐of‐attack NACA0015 airfoil will be presented to demonstrate the accuracy and efficiency of our current implementation. Copyright © 2008 John Wiley & Sons, Ltd.     

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered finite volume/finite element method for solution of the two dimensional (2D) incompressible Navier-Stokes equations on unstructured grids.An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling.The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by joining the centroid of cells sharing the common vertex.For the temporal integration of the momentum equations,an implicit second-order scheme is utilized to enhance the computational stability and eliminate the time step limit due to the diffusion term.The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite element method (FEM).The momentum interpolation is used to damp out the spurious pressure wiggles.The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both velocity and pressure.The classic test cases,the lid-driven cavity flow,the skew cavity flow and the backward-facing step flow,show that numerical results are in good agreement with the published benchmark solutions.     

A finite difference analysis of the extrudate swell problem

The extrudate swell phenomenon of a purely viscous fluid is analysed by solving simultaneously the Cauchy momentum equations along with the continuity equation by means of a finite difference method. The circular and planar jet flows of Newtonian and power-law fluids are simulated using a control volume finite difference method suggested by Patankar called SIMPLER (semi-implicit method for pressure-linked equations). This method uses the velocity components and pressure as the primitive variables and employs a staggered grid and control volume for each separate variable. The numerical results show good agreement with the analytical solution of the axisymmetric stick-slip problem and exhibit a Newtonian swelling ratio of 13.2% or 19.2% for a capillary or slit die respectively in accordance with previously reported experimental and numerical results. Shear thinning results in a decrease in swelling ratio, as does the introduction of gravity and surface tension.     

A numerical‐variational procedure for laminar flow in curved square ducts

A new numerical method is presented for the solution of the Navier–Stokes and continuity equations governing the internal incompressible flows. The method denoted as the CVP method consists in the numerical solution of these equations in conjunction with three additional variational equations for the continuity, the vorticity and the pressure field, using a non‐staggered grid. The method is used for the study of the characteristics of the laminar fully developed flows in curved square ducts. Numerical results are presented for the effects of the flow parameters like the curvature, the Dean number and the stream pressure gradient on the velocity distributions, the friction factor and the appearance of a pair of vortices in addition to those of the familiar secondary flow. The accuracy of the method is discussed and the results are compared with those obtained by us, using a variation of the velocity–pressure linked equation methods denoted as the PLEM method and the results obtained by other methods. Copyright © 2004 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.     

An overlapping control volume method for the Navier–Stokes equations on non‐staggered grids

An algorithm, based on the overlapping control volume (OCV) method, for the solution of the steady and unsteady two‐dimensional incompressible Navier–Stokes equations in complex geometry is presented. The primitive variable formulation is solved on a non‐staggered grid arrangement. The problem of pressure–velocity decoupling is circumvented by using momentum interpolation. The accuracy and effectiveness of the method is established by solving five steady state and one unsteady test problems. The numerical solutions obtained using the technique are in good agreement with the analytical and benchmark solutions available in the literature. On uniform grids, the method gives second‐order accuracy for both diffusion‐ and convection‐dominated flows. There is little loss of accuracy on grids that are moderately non‐orthogonal. Copyright © 1999 John Wiley & Sons, Ltd.     

Schwarz domain decomposition for the incompressible Navier–Stokes equations in general co‐ordinates

This paper describes a domain decomposition method for the incompressible Navier–Stokes equations in general co‐ordinates. Domain decomposition techniques are needed for solving flow problems in complicated geometries while retaining structured grids on each of the subdomains. This is the so‐called block‐structured approach. It enables the use of fast vectorized iterative methods on the subdomains. The Navier–Stokes equations are discretized on a staggered grid using finite volumes. The pressure‐correction technique is used to solve the momentum equations together with incompressibility conditions. Schwarz domain decomposition is used to solve the momentum and pressure equations on the composite domain. Convergence of domain decomposition is accelerated by a GMRES Krylov subspace method. Computations are presented for a variety of flows. Copyright © 2000 John Wiley & Sons, Ltd.     

Wet/dry areas method for interfacial (free surface) flow

In this paper, a new numerical method is developed for two‐dimensional interfacial (free surface) flows, based on the control volume method and conservative integral form of the Navier–Stokes equations with a standard staggered grid. The new method deploys two continuity equations, the continuity equation of the mass conservation for better convergence of the implicit scheme and the continuity equation of the volume conservation for the equation of pressure correction. The convection terms (the total momentum flux) on the surfaces of control volume are accurately calculated from the wet area exposed to the water, and the dry area exposed to the air. The numerical results produced by the new numerical method agree very well with the analytical solution, experimental images and experimentally measured velocity. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity field by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.     

Numerical method for calculation of the incompressible flow in general curvilinear co‐ordinates with double staggered grid

A solution methodology has been developed for incompressible flow in general curvilinear co‐ordinates. Two staggered grids are used to discretize the physical domain. The first grid is a MAC quadrilateral mesh with pressure arranged at the centre and the Cartesian velocity components located at the middle of the sides of the mesh. The second grid is so displaced that its corners correspond to the centre of the first grid. In the second grid the pressure is placed at the corner of the first grid. The discretized mass and momentum conservation equations are derived on a control volume. The two pressure grid functions are coupled explicitly through the boundary conditions and implicitly through the velocity of the field. The introduction of these two grid functions avoids an averaging of pressure and velocity components when calculating terms that are generated in general curvilinear co‐ordinates. The SIMPLE calculation procedure is extended to the present curvilinear co‐ordinates with double grids. Application of the methodology is illustrated by calculation of well‐known external and internal problems: viscous flow over a circular cylinder, with Reynolds numbers ranging from 10 to 40, and lid‐driven flow in a cavity with inclined walls are examined. The numerical results are in close agreement with experimental results and other numerical data. Copyright © 2003 John Wiley & Sons, Ltd.     

用RNG K-E模式数值模拟180°弯道内的湍流分离流动

将Ｙａｋｈｏｔ与Ｏｒｓｚａｇ新近提出的ＲＮＧＫ－ε湍流模式推广应用于１８０°强曲率弯道内的湍流分离流动的数值模拟，计算在任意曲线坐标下进行，并采用速度协变分量作为求解变量以保证计算的高度稳定性，控制方程的求解采用通常的控制容积法，文中给出了详细的数值计算结果，并与实验结果进行了比较，结果表明，ＲＮＧＫ－ε湍流模式能有效地模拟有强曲率影响的湍流分离流动，展示了这一模式在工程湍流计算中的前景     

A study of penalty elements for incompressible laminar flows

A finite element model is developed based on the penalty formulation to study incompressible laminar flows. The study includes a number of new quadrilateral and triangular elements for 2-dimensional flows and a number of new hexahedral and tetrahedral elements for 3-dimensional flows. All elements employ continuous velocity approximations and discontinuous pressure approximations respecting the LBB condition of numerical instability. An incremental Newton–Raphson method coupled with the Broyden method is used to solve the non-linear equations. Several numerical examples (colliding flow, cavity flow, etc.) are presented to assess the efficiency of elements.