A segregated spectral finite element method for the 2D transient incompressible Navier-Stokes equations

Abstract In this paper, the spectral element approximation and the velocity–pressuredecoupling method implementing the SIMPLE (semi-implicit method for pressure linked equations) algorithm are first combined to form a high-order segregated scheme for the solution of the two-dimensional transient incompressible Navier–Stokes equations. In contrast to previous segregated finite element methods based on the SIMPLE algorithm, the pressure equation is derived from the continuity equation using the element matrices to ensure convergence. High-order element basis functions are adopted, which greatly reduces the number of nodal points used in the calculation. The validation test that has an analytical solution demonstrates the high accuracy and convergence rate of the method. The flow in a lid-driven cavity with different inclination angles and the flow over a backward-facing step are investigated to further illustrate the performance of the scheme. The computed results are in excellent agreement with the benchmark solutions. The almost periodic solution for the flow of Re=10000 in a lid-driven square cavity is also captured by the present scheme.

[1]  R. J. Schnipke,et al.  An equal-order velocity-pressure formulation that does not exhibit spurious pressure modes , 1986 .

[2]  C. G. du Toit,et al.  A segregated finite element solution for non-isothermal flow , 2000 .

[3]  W. Q. Tao,et al.  IMPLEMENTATION OF CLEAR ALGORITHM ON COLLOCATED GRID SYSTEM AND APPLICATION EXAMPLES , 2004 .

[4]  Jie Shen,et al.  On error estimates of the projection methods for the Navier-Stokes equations: Second-order schemes , 1996, Math. Comput..

[5]  Jie Shen,et al.  Parallel spectral-element direction splitting method for incompressible Navier-Stokes equations , 2014 .

[6]  D. Gartling.  A test problem for outflow boundary conditions—flow over a backward-facing step , 1990 .

[7]  Yih-Ferng Peng,et al.  Transition in a 2-D lid-driven cavity flow , 2003 .

[8]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[9]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

[10]  Chuanju Xu,et al.  Analysis of iterative methods for the viscous/inviscid coupled problem via a spectral element approximation , 2000 .

[11]  Ya-Ling He,et al.  An Efficient Segregated Algorithm for Incompressible Fluid Flow and Heat Transfer Problems—IDEAL (Inner Doubly Iterative Efficient Algorithm for Linked Equations) Part II: Application Examples , 2008 .

[12]  Spencer J. Sherwin,et al.  Direct numerical simulations of the flow around wings with spanwise waviness at a very low Reynolds number , 2017 .

[13]  B. Armaly,et al.  Experimental and theoretical investigation of backward-facing step flow , 1983, Journal of Fluid Mechanics.

[14]  Atipong Malatip,et al.  Combined Streamline Upwind Petrov Galerkin method and segregated finite element algorithm for conjugate heat transfer problems , 2006 .

[15]  Spencer J. Sherwin,et al.  An adaptable parallel algorithm for the direct numerical simulation of incompressible turbulent flows using a Fourier spectral/hp element method and MPI virtual topologies , 2016, Comput. Phys. Commun..

[16]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[17]  Li-Chieh Hsu,et al.  A study of flow patterns for staggered cylinders at low Reynolds number by spectral element method , 2017 .

[18]  A. Patera.  A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[19]  Brandon E. Merrill,et al.  A spectrally accurate method for overlapping grid solution of incompressible Navier-Stokes equations , 2016, J. Comput. Phys..

[20]  Subhashree Mohapatra,et al.  A Non-Conforming Least Squares Spectral Element Formulation for Oseen Equations with Applications to Navier-Stokes Equations , 2016 .

[21]  C. Bruneau,et al.  The 2D lid-driven cavity problem revisited , 2006 .

[22]  I. Demirdzic,et al.  Fluid flow and heat transfer test problems for non‐orthogonal grids: Bench‐mark solutions , 1992 .

[23]  Suhas V. Patankar,et al.  A Calculation Procedure for Two-Dimensional Elliptic Situations , 1981 .

[24]  Ali Beskok,et al.  An algebraic factorisation scheme for spectral element solution of incompressible flow and scalar transport , 2010 .

[25]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[26]  Yusheng Feng,et al.  Computational Algorithms for Solving Spectral/$hp$ Stabilized Incompressible Flow Problems , 2016 .

[27]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[28]  Zhiguo Qu,et al.  An Efficient Segregated Algorithm for Incompressible Fluid Flow and Heat Transfer Problems—IDEAL (Inner Doubly Iterative Efficient Algorithm for Linked Equations) Part I: Mathematical Formulation and Solution Procedure , 2008 .

[29]  E. Erturk,et al.  Discussions on driven cavity flow , 2004, cs/0411048.

[30]  Anupam Dewan,et al.  A Multigrid-Accelerated Code on Graded Cartesian Meshes for 2D Time-Dependent Incompressible Viscous Flows , 2010 .

[31]  E. Erturk,et al.  Numerical solutions of 2‐D steady incompressible driven cavity flow at high Reynolds numbers , 2004, ArXiv.

[32]  E. Barragy,et al.  STREAM FUNCTION-VORTICITY DRIVEN CAVITY SOLUTION USING p FINITE ELEMENTS , 1997 .