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

This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （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 join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational 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 el- ement 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 veloc- ity 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.

Combined adaptive meshing technique and characteristic-based split algorithm for viscous incompressible flow analysis

A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.

A Staggered Grid Method for Solving Incompressible Flow on Unstructured Meshes

A finite volume method based unstructured grid is presented to solve the two dimensional viscous and incompressible flow.The method is based on the pressure-correction concept and solved by using a semi-staggered grid technique.The computational procedure can handle cells of arbitrary shapes,although solutions presented in this paper were only involved with triangular and quadrilateral cells.The pressure or pressure-correction value was stored on the vertex of cells.The mass conservation equation was discretized on the dual cells surrounding the vertex of primary cells,while the velocity components and other scale variables were saved on the central of primary cells.Since the semi-staggered arrangement can’t guarantee a strong coupling relationship between pressure and velocity,thus a weak coupling relationship leads to the oscillations for pressure and velocity.In order to eliminate such an oscillation,a special interpolation scheme was used to construct the pressure-correction equation.Computational results of several viscous flow problems show good agreement with the analytical or numerical results in previous literature.This semi-staggered grid method can be applied to arbitrary shape elements,while it has the most efficiency for triangular cells.

A class of nonconforming quadrilateral finite elements for incompressible flow

This paper focuses on the low-order nonconforming rectangular and quadrilateral finite elements approximation of incompressible flow.Beyond the previous research works,we propose a general strategy to construct the basis functions.Under several specific constraints,the optimal error estimates are obtained,i.e.,the first order accuracy of the velocities in H1-norm and the pressure in L2-norm,as well as the second order accuracy of the velocities in L2-norm.Besides,we clarify the differences between rectangular and quadrilateral finite element approximation.In addition,we give several examples to verify the validity of our error estimates.

A control volume based finite element method for simulating incompressible two-phase flow in heterogeneous porous media and its application to reservoir engineering

Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.

FINITE ELEMENT ANALYSIS OF COMPLEX INCOMPRESSIBLE FLOWS

The complex flow has been simulated with finite element method for Navier-Stokes equations. Both two and three dimensional, laminar and turbulent problems were included. The velocity-pressure decoupling method was adopted and related algorithm was discussed. Effects of turbulence were accounted by large eddy simulation.

OPTIMAL MIXED H-P FINITE ELEMENT METHODS FOR STOKES AND NON-NEWTONIAN FLOW

Presents an h -- p finite element methods based upon a mixed variational formulation for the three-field Stokes equations and linearized Non-Newtonian flow. Computation of the algebraic system generated from Problem H[sub h]; Methodology; Results and discussion.

High-Order Fully Discrete Energy Diminishing Evolving Surface Finite Element Methods for a Class of Geometric Curvature Flows

This article concerns the construction of high-order energy-decaying numerical methods for gradient flows of evolving surfaces with curvature-dependent energy functionals.The semidiscrete evolving surface finite element method is derived based on the calculus of variation of the semidiscrete surface energy functional.This makes the semidiscrete problem naturally inherit the energy decay structure.With this property,the semidiscrete problem is furthermore formulated as a gradient flow system of ODEs.The averaged vector-field collocation method is used for time discretization of the ODEs to preserve energy decay at the fully discrete level while achieving high-order accuracy in time.Extensive numerical examples are provided to illustrate the accuracy and energy diminishing property of the proposed method,as well as the effectiveness of the method in capturing singularities in the evolution of closed surfaces.

Phase-field modeling of dendritic growth under forced flow based on adaptive finite element method

A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.

STABLE FOURTH-ORDER STREAM-FUNCTION METHODS FOR INCOMPRESSIBLE FLOWS WITH BOUNDARIES

Fourth-order stream-function methods are proposed for the time dependent, incom- pressible Navier-Stokes and Boussinesq equations. Wide difference stencils are used instead of compact ones and the boundary terms are handled by extrapolating the stream-function values inside the computational domain to grid points outside, up to fourth-order in the noslip condition. Formal error analysis is done for a simple model problem, showing that this extrapolation introduces numerical boundary layers at fifth-order in the stream-function. The fourth-order convergence in velocity of the proposed method for the full problem is shown numerically.

Two-Scale Picard Stabilized Finite Volume Method for the Incompressible Flow

In this paper,we consider a two-scale stabilized finite volume method for the two-dimensional stationary incompressible flow approximated by the lowest equalorder element pair P_(1)−P_(1)which do not satisfy the inf-sup condition.The two-scale method consist of solving a small non-linear system on the coarse mesh and then solving a linear Stokes equations on the fine mesh.Convergence of the optimal order in the H1-norm for velocity and the L^(2)-norm for pressure are obtained.The error analysis shows there is the same convergence rate between the two-scale stabilized finite volume solution and the usual stabilized finite volume solution on a fine mesh with relation h=O(H^(2)).Numerical experiments completely confirm theoretic results.Therefore,this method presented in this paper is of practical importance in scientific computation.

Pseudostress-BasedMixed Finite Element Methods for the Stokes Problem in Rn with Dirichlet Boundary Conditions.I:A Priori Error Analysis

We consider a non-standard mixed method for the Stokes problem in Rn,n∈{2,3},with Dirichlet boundary conditions,in which,after using the incompressibility condition to eliminate the pressure,the pseudostress tensor s and the velocity vector u become the only unknowns.Then,we apply the Babuˇska-Brezzi theory to prove the well-posedness of the corresponding continuous and discrete formulations.In particular,we show that Raviart-Thomas elements of order k≥0 for s and piecewise polynomials of degree k for u ensure unique solvability and stability of the associated Galerkin scheme.In addition,we introduce and analyze an augmented approach for our pseudostress-velocity formulation.The methodology employed is based on the introduction of the Galerkin least-squares type terms arising from the constitutive and equilibrium equations,and the Dirichlet boundary condition for the velocity,all of them multiplied by suitable stabilization parameters.We show that these parameters can be chosen so that the resulting augmented variational formulation is defined by a strongly coercive bilinear form,whence the associated Galerkin scheme becomes well posed for any choice of finite element subspaces.For instance,Raviart-Thomas elements of order k≥0 for s and continuous piecewise polynomials of degree k+1 for u become a feasible choice in this case.Finally,extensive numerical experiments illustrating the good performance of the methods and comparing them with other procedures available in the literature,are provided.

A Numerical Comparison of Outflow Boundary Conditions for Spectral Element Simulations of Incompressible Flows

Outflow boundary conditions(OBCs)are investigated for calculation of incompressible flows by spectral element methods.Several OBCs,including essentialtype,natural-type,periodic-type and advection-type,are compared by carrying out a series of numerical experiments.Especially,a simplified form of the so-called Orlanski’s OBCs is proposed in the context of spectral element methods,for which a new treatment technique is used.The purpose of this paper is to find stable low-reflective OBCs,suitable and flexible for use of spectral element methods in simulation of incompressible flows in complex geometries.The computation is firstly carried out for a 2D simulation of Poiseuille-B´enard channel flow with Re=10,Ri=150 and Pr=2/3.This flow serves as a useful example to demonstrate the applicability of the proposed OBCs because it exhibits a feature of vortex shedding propagating through the outflow boundary.Then a 3D flow around an obstacle is computed to show the efficiency in the case of more general geometries.Among the tested OBCs,the advection-type OBCs are proven to have better behavior as compared with the others.

Finite Element Numerical Simulation and PIV Measurement of Flow Field inside Metering-in Spool Valve

The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.

THE FINITE VOLUME PROJECTION METHOD WITH HYBRID UNSTRUCTURED TRIANGULAR COLLOCATED GRIDS FOR INCOMPRESSIBLE FLOWS

In this article a finite volume method is proposed to solve viscous incompressible Navier-Stokes equations in two-dimensional regions with corners and curved boundaries. A hybrid collocated-grid variable arrangement is adopted, in which the velocity and pressure are stored at the centroid and the circumcenters of the triangular control cell, respectively. The cell flux is defined at the mid-point of the cell face. Second-order implicit time integration schemes are used for convection and diffusion terms. The second-order upwind scheme is used for convection fluxes. The present method is validated by results of several viscous flows.

ANISOTROPIC MULTISTAGE FINITE ELEMENT METHOD FOR TWO DIMENSIONAL VISCOUS TRANSONIC FLOW IN TURBOMACHINERY

A new method based on the anisotropic tensor force finite element and Taylor-Galerkin finite element is presented in the present paper.Its application to two-dimensional viscous transonic flow in turbomachinery improves the conver- gence rate and stability of calculation,and the results obtained agree well with the experimental measurements.

SUPG finite element method based on penalty function for lid-driven cavity flow up to Re = 27500

A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.

Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from 100 to 1000. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.

Polygonal Finite Element for Two-Dimensional Lid-Driven Cavity Flow

This paper investigates a polygonal finite element(PFE)to solve a two-dimensional(2D)incompressible steady fluid problem in a cavity square.It is a well-known standard benchmark(i.e.,lid-driven cavity flow)-to evaluate the numerical methods in solving fluid problems controlled by the Navier-Stokes(N-S)equation system.The approximation solutions provided in this research are based on our developed equal-order mixed PFE,called Pe1Pe1.It is an exciting development based on constructing the mixed scheme method of two equal-order discretisation spaces for both fluid pressure and velocity fields of flows and our proposed stabilisation technique.In this research,to handle the nonlinear problem of N-S,the Picard iteration scheme is applied.Our proposed method’s performance and convergence are validated by several simulations coded by commercial software,i.e.,MATLAB.For this research,the benchmark is executed with variousReynolds numbers up to the maximum Re=1000.All results then numerously compared to available sources in the literature.

Parallel finite element computation of incompressible magnetohydrodynamics based on three iterations

Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are presented.These approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency component.Optimal error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are given.Some numerical examples are implemented to verify the algorithm.