A poly‐region boundary element method for incompressible viscous fluid flows

A boundary element method (BEM) for steady viscous fluid flow at high Reynolds numbers is presented. The new integral formulation with a poly-region approach involves the use of the convective kernel with slight compressibility that was previously employed by Grigoriev and Fafurin [1] for driven cavity flows with Reynolds numbers up to 1000. In order to avoid the overdeterminancy of the global set of equations when using eight-noded rectangular volume cells from that previous work, 12-noded hexagonal volume regions are introduced. As a result, the number of linearly independent integral equations for each node becomes equal to the degrees of freedom of the node. The numerical results for square-driven cavity flow having Reynolds numbers up to 5000 are compared to those obtained by Ghia et al. [2] and demonstrate a high level of accuracy even in resolving the secondary vortices at the corners of the cavity. Next, a comprehensive study is done for backward-facing step flows at Re=500 and 800 using the BEM, along with a standard Galerkin-based finite element methods (FEM). The numerical methods are in excellent agreement with the benchmark solution published by Gartling [3]. However, several additional aspects of the problem are also considered, including the effect of the inflow boundary location and the traction singularity at the step corner. Furthermore, a preliminary comparative study of the poly-region BEM versus the standard FEM indicates that the new method is more than competitive in terms of accuracy and efficiency. Copyright © 1999 John Wiley & Sons, Ltd.

[1]  A. Spence,et al.  Is the steady viscous incompressible two‐dimensional flow over a backward‐facing step at Re = 800 stable? , 1993 .

[2]  M. Bercovier,et al.  A finite element for the numerical solution of viscous incompressible flows , 1979 .

[3]  C. W. Oseen,et al.  Neuere Methoden und Ergebnisse in der Hydrodynamik , 1927 .

[4]  François Thomasset,et al.  Implementation of Finite Element Methods for Navier-Stokes Equations , 1981 .

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

[6]  N. Tosaka,et al.  Numerical solutions of steady incompressible viscous flow problems by the integral equation method , 1986 .

[7]  J. L. Sohn,et al.  Evaluation of FIDAP on some classical laminar and turbulent benchmarks. [FluId Dynamics Analysis Package , 1988 .

[8]  Jacques Periaux,et al.  Analysis of laminar flow over a backward facing step : a GAMM-workshop , 1984 .

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

[10]  R. Sani,et al.  Résumé and remarks on the open boundary condition minisymposium , 1994 .

[11]  Brian Launder,et al.  Numerical methods in laminar and turbulent flow , 1983 .

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

[13]  C. Brebbia,et al.  Boundary Element Techniques , 1984 .

[14]  L. G. Leal,et al.  Low-Reynolds-number flow past cylindrical bodies of arbitrary cross-sectional shape , 1986, Journal of Fluid Mechanics.

[15]  I. Duff,et al.  Direct Methods for Sparse Matrices , 1987 .

[16]  A. Fafurin,et al.  A boundary element method for steady viscous fluid flow using penalty function formulation , 1997 .

[17]  Gary F. Dargush,et al.  Boundary element method for steady incompressible thermoviscous flow , 1991 .

[18]  R. Tanner,et al.  Numerical solution of viscous flows using integral equation methods , 1983 .

[19]  A. Acrivos,et al.  Stokes flow past a particle of arbitrary shape: a numerical method of solution , 1975, Journal of Fluid Mechanics.

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

[21]  P. K. Banerjee,et al.  Development of an integrated BEM approach for hot fluid structure interaction. Annual report, November 1986-November 1987 , 1991 .

[22]  P. Hood,et al.  A numerical solution of the Navier-Stokes equations using the finite element technique , 1973 .

[23]  George Em Karniadakis,et al.  Onset of three-dimensionality, equilibria, and early transition in flow over a backward-facing step , 1991, Journal of Fluid Mechanics.

[24]  A BOUNDARY ELEMENT FORMULATION USING THE SIMPLE METHOD , 1993 .

[25]  M. Bush.  Modelling two-dimensional flow past arbitrary cylindrical bodies using boundary element formulations , 1983 .

[26]  C. Pozrikidis.  Boundary Integral and Singularity Methods for Linearized Viscous Flow: Index , 1992 .

[27]  P. K. Banerjee,et al.  Boundary element methods in engineering science , 1981 .