Parallel finite element computation of unsteady incompressible flows

A parallel semi-explicit iterative finite element computational procedure for modelling unsteady incompressible fluid flows is presented. During the procedure, element flux vectors are calculated in parallel and then assembled into global flux vectors. Equilibrium iterations which introduce some ‘local implicitness’ are performed at each time step. The number of equilibrium iterations is governed by an implicitness parameter. The present technique retains the advantages of purely explicit schemes, namely (i) the parallel speed-up is equal to the number of parallel processors if the small communication overhead associated with purely explicit schemes is ignored and (ii) the computation time as well as the core memory required is linearly proportional to the number of elements. The incompressibility condition is imposed by using the artificial compressibility technique. A pressure-averaging technique which allows the use of equal-order interpolations for both velocity and pressure, this simplifying the formulation, is employed. Using a standard Galerkin approximation, three benchmark steady and unsteady problems are solved to demonstrate the accuracy of the procedure. In all calculations the Reynolds number is less than 500. At these Reynolds numbers it was found that the physical dissipation is sufficient to stabilize the convective term with no need for additional upwind-type dissipation. © 1998 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.     

A stabilized finite element method for the incompressible Navier–Stokes equations using a hierarchical basis

Stabilized finite element methods have been shown to yield robust, accurate numerical solutions to both the compressible and incompressible Navier–Stokes equations for laminar and turbulent flows. The present work focuses on the application of higher‐order, hierarchical basis functions to the incompressible Navier–Stokes equations using a stabilized finite element method. It is shown on a variety of problems that the most cost‐effective simulations (in terms of CPU time, memory, and disk storage) can be obtained using higher‐order basis functions when compared with the traditional linear basis. In addition, algorithms will be presented for the efficient implementation of these methods within the traditional finite element data structures. Copyright © 2001 John Wiley & Sons, Ltd.     

A least-squares finite element method for incompressible Navier-Stokes problems

A least-squares finite element method based on the velocity–pressure–vorticity formulation was proposed for solving steady incompressible Navier-Stokes problems. This method leads to a minimization problem rather than to the saddle point problem of the classic mixed method and can thus accommodate equal-order interpolations. The method has no parameter to tune. The associated algebraic system is symmetric and positive definite. In order to show the validity of the method for high-Reynolds-number problems, this paper provides numerical results for cavity flow at Reynolds number up to 10 000 and backward-facing step flow at Reynolds number up to 900.     

Two‐level finite element method with a stabilizing subgrid for the incompressible Navier–Stokes equations

We consider the Galerkin finite element method for the incompressible Navier–Stokes equations in two dimensions. The domain is discretized into a set of regular triangular elements and the finite‐dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual‐free bubble functions. To find the bubble part of the solution, a two‐level finite element method with a stabilizing subgrid of a single node is described, and its application to the Navier–Stokes equation is displayed. Numerical approximations employing the proposed algorithm are presented for three benchmark problems. The results show that the proper choice of the subgrid node is crucial in obtaining stable and accurate numerical approximations consistent with the physical configuration of the problem at a cheap computational cost. Copyright © 2008 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid: 3D validation

A three‐dimensional finite element method for incompressible multiphase flows with capillary interfaces is developed based on a (formally) second‐order projection scheme. The discretization is on a fixed (Eulerian) reference grid with an edge‐based local h‐refinement in the neighbourhood of the interfaces. The fluid phases are identified and advected using the level‐set function. The reference grid is then temporarily reconnected around the interface to maintain optimal interpolations accounting for the singularities of the primary variables. Using a time splitting procedure, the convection substep is integrated with an explicit scheme. The remaining generalized Stokes problem is solved by means of a pressure‐stabilized projection. This method is simple and efficient, as demonstrated by a wide range of difficult free‐surface validation problems, considered in the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

A finite element solution of the time-dependent incompressible Navier–Stokes equations using a modified velocity correction method

A finite element solution of the two-dimensional incompressible Navier–Stokes equations has been developed. The present method is a modified velocity correction approach. First an intermediate velocity is calculated, and then this is corrected by the pressure gradient which is the solution of a Poisson equation derived from the continuity equation. The novelty, in this paper, is that a second-order Runge–Kutta method for time integration has been used. Discretization in space is carried out by the Galerkin weighted residual method. The solution is in terms of primitive variables, which are approximated by polynomial basis functions defined on three-noded, isoparametric triangular elements. To demonstrate the present method, two examples are provided. Results from the first example, the driven cavity flow problem, are compared with previous works. Results from the second example, uniform flow past a cylinder, are compared with experimental data.     

A nodal‐based implementation of a stabilized finite element method for incompressible flow problems

The objective of this paper is twofold. First, a stabilized finite element method (FEM) for the incompressible Navier–Stokes is presented and several numerical experiments are conducted to check its performance. This method is capable of dealing with all the instabilities that the standard Galerkin method presents, namely the pressure instability, the instability arising in convection‐dominated situations and the less popular instabilities found when the Navier–Stokes equations have a dominant Coriolis force or when there is a dominant absorption term arising from the small permeability of the medium where the flow takes place. The second objective is to describe a nodal‐based implementation of the finite element formulation introduced. This implementation is based on an a priori calculation of the integrals appearing in the formulation and then the construction of the matrix and right‐hand side vector of the final algebraic system to be solved. After appropriate approximations, this matrix and this vector can be constructed directly for each nodal point, without the need to loop over the elements, thus making the calculations much faster. In order to be able to do this, all the variables have to be defined at the nodes of the finite element mesh, not on the elements. This is also so for the stabilization parameters of the formulation. However, doing this gives rise to questions regarding the consistency and the conservation properties of the final scheme, which are addressed in this paper. Copyright © 2000 John Wiley & Sons, Ltd.     

A Hermite finite element method for incompressible fluid flow

We describe some Hermite stream function and velocity finite elements and a divergence‐free finite element method for the computation of incompressible flow. Divergence‐free velocity bases defined on (but not limited to) rectangles are presented, which produce pointwise divergence‐free flow fields (∇· u _h≡0). The discrete velocity satisfies a flow equation that does not involve pressure. The pressure can be recovered as a function of the velocity if needed. The method is formulated in primitive variables and applied to the stationary lid‐driven cavity and backward‐facing step test problems. Copyright © 2009 John Wiley & Sons, Ltd.     

A 3D incompressible Navier–Stokes velocity–vorticity weak form finite element algorithm

The velocity–vorticity formulation is selected to develop a time‐accurate CFD finite element algorithm for the incompressible Navier–Stokes equations in three dimensions.The finite element implementation uses equal order trilinear finite elements on a non‐staggered hexahedral mesh. A second order vorticity kinematic boundary condition is derived for the no slip wall boundary condition which also enforces the incompressibility constraint. A biconjugate gradient stabilized (BiCGSTAB) sparse iterative solver is utilized to solve the fully coupled system of equations as a Newton algorithm. The solver yields an efficient parallel solution algorithm on distributed‐memory machines, such as the IBM SP2. Three dimensional laminar flow solutions for a square channel, a lid‐driven cavity, and a thermal cavity are established and compared with available benchmark solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

A high‐order accurate discontinuous Galerkin finite element method for laminar low Mach number flows

In this paper we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of laminar flow simulations at low Mach numbers using an implicit scheme. The algorithm is based on the flux preconditioning approach, which modifies only the dissipative terms of the numerical flux. This formulation is quite simple to implement in existing implicit DG codes, it overcomes the time‐stepping restrictions of explicit multistage algorithms, is consistent in time and thus applicable to unsteady flows. The performance of the method is demonstrated by solving the flow around a NACA0012 airfoil and on a flat plate, at different low Mach numbers using various degrees of polynomial approximations. Computations with and without flux preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers. The time accurate solution of unsteady flow is also demonstrated by solving the vortex shedding behind a circular cylinder at the Reynolds number of 100. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of free and forced convection incompressible flows using an adaptive parallel/vector finite element procedure

The numerical simulation of complex flows demands efficient algorithms and fast computer platforms. The use of adaptive techniques permits adjusting the discretisation according to the analysis requirements, but creates variable computational loads that are difficult to manage in a parallel/vector program. This paper describes the approach adopted to implement an adaptive finite element incompressible Navier–Stokes solver on the Cray J90 machine. Performance measurements for the simulation of free and forced convection incompressible flows indicate that the techniques employed result in a fast parallel/vector code. Copyright © 1999 John Wiley & Sons, Ltd.     

Gauge finite element method for incompressible flows

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J‐G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher‐order) finite elements. This method can achieve high‐order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright © 2000 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.     

Analysis of multifluid flows with large time steps using the particle finite element method

Multifluids are those fluids in which their physical properties (viscosity or density) vary internally and abruptly forming internal interfaces that introduce a large nonlinearity in the Navier–Stokes equations. For this reason, standard numerical methods require very small time steps in order to solve accurately the internal interface position. In a previous paper, the authors developed a particle‐based method (named particle finite element method (PFEM)) based on a Lagrangian formulation and FEM for solving the fluid mechanics equations for multifluids. PFEM was capable of achieving accurate results, but the limitation of small time steps was still present. In this work, a new strategy concerning the time integration for the analysis of multifluids is developed allowing time steps one order of magnitude larger than the previous method. The advantage of using a Lagrangian solution with PFEM is shown in several examples. All kind of heterogeneous fluids (with different densities or viscosities), multiphase flows with internal interfaces, breaking waves, and fluid separation may be easily solved with this methodology without the need of small time steps. Copyright © 2014 John Wiley & Sons, Ltd.     

Implementation of a stabilized finite element formulation for the incompressible Navier–Stokes equations based on a pressure gradient projection

We discuss in this paper some implementation aspects of a finite element formulation for the incompressible Navier–Stokes equations which allows the use of equal order velocity–pressure interpolations. The method consists in introducing the projection of the pressure gradient and adding the difference between the pressure Laplacian and the divergence of this new field to the incompressibility equation, both multiplied by suitable algorithmic parameters. The main purpose of this paper is to discuss how to deal with the new variable in the implementation of the algorithm. Obviously, it could be treated as one extra unknown, either explicitly or as a condensed variable. However, we take for granted that the only way for the algorithm to be efficient is to uncouple it from the velocity–pressure calculation in one way or another. Here we discuss some iterative schemes to perform this uncoupling of the pressure gradient projection (PGP) from the calculation of the velocity and the pressure, both for the stationary and the transient Navier–Stokes equations. In the first case, the strategies analyzed refer to the interaction of the linearization loop and the iterative segregation of the PGP, whereas in the second the main dilemma concerns the explicit or implicit treatment of the PGP. Copyright © 2001 John Wiley & Sons, Ltd.     

Transient solutions for three-dimensional lid-driven cavity flows by a least-squares finite element method

A time-accurate least-squares finite element method is used to simulate three-dimensional flows in a cubic cavity with a uniform moving top. The time- accurate solutions are obtained by the Crank-Nicolson method for time integration and Newton linearization for the convective terms with extensive linearization steps. A matrix-free algorithm of the Jacobi conjugate gradient method is used to solve the symmetric, positive definite linear system of equations. To show that the least-squares finite element method with the Jacobi conjugate gradient technique has promising potential to provide implicit, fully coupled and time-accurate solutions to large-scale three-dimensional fluid flows, we present results for three-dimensional lid-driven flows in a cubic cavity for Reynolds numbers up to 3200.     

A finite volume method for viscous incompressible flows using a consistent flux reconstruction scheme

An incompressible Navier–Stokes solver using curvilinear body‐fitted collocated grid has been developed to solve unconfined flow past arbitrary two‐dimensional body geometries. In this solver, the full Navier–Stokes equations have been solved numerically in the physical plane itself without using any transformation to the computational plane. For the proper coupling of pressure and velocity field on collocated grid, a new scheme, designated ‘consistent flux reconstruction’ (CFR) scheme, has been developed. In this scheme, the cell face centre velocities are obtained explicitly by solving the momentum equations at the centre of the cell faces. The velocities at the cell centres are also updated explicitly by solving the momentum equations at the cell centres. By resorting to such a fully explicit treatment considerable simplification has been achieved compared to earlier approaches. In the present investigation the solver has been applied to unconfined flow past a square cylinder at zero and non‐zero incidence at low and moderate Reynolds numbers and reasonably good agreement has been obtained with results available from literature. Copyright © 2006 John Wiley & Sons, Ltd.