An Incompressible Flow Computation by a Hierarchical Iterative and a Modified Residual Method

계층적 반복과 수정 잔여치법에 의한 비압축성 유동 계산

  • Published : 2004.09.01

Abstract

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

Keywords

References

  1. Dendy, J.E., 'Two methods of Galerkin type Optimum L-2 Rates of Convergence for First order Hyperbolics,' SIAM J. Numer. Anal., Vol. 11, (1974), p.637-653
  2. WahIbun, L.B., 'A Dissipative Galerkin method Applied to some quasi-linear Hyperbolic equations,' R.A.I.R.O., Vol. 29, (1975), p.1037-1045
  3. Hughes, T.J.R. and Brooks, A.N.. 'A Multi-dimensional upwind scheme with no crosswind diffusion,' T.J.R. Hughes(Ed.), Finite Element Methods for Convection Dominated Flows, Vol. 34, (1979), p.19-35
  4. Brooks, A. and Hughes, T.J.R., 'Streamline upwind/Petrov-Galerkin formulation for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations,' Comput. Methods Appl. Mech. Engrg., Vol. 32, (1982), p.199-259
  5. Ilinca, F., Hetu, J.F. and Pelletier, D.. 'On stabilized finite element formulations for incompressible advective-diffusive transport and fluid flow problems,' Comput. Methods Appl. Mech Engrg., Vol. 188, (2000), p.235-255
  6. Mitchell, J.A. and Reddy, J.N., 'A hierarchical iterative procedure for the analysis of composite laminates,' Comput. Methods Appl. Mech. Engrg., Vol. 181, (2000), p.237-260
  7. Whiting, C. H. and Jansen, K.E., 'A stabilized finite element method for the incompressible Navier-Stokes equations using a hierarchical basis,' Int. J. Numer. Meth. Fluids, Vol. 35, (2001), p.93-116
  8. Kovasznay, L.I.G.. 'Laminar flow behind a two-dimensional grid,' Proceedings of the Cambridge PhiIosophicaI Society, Vol. 44 (1948), p.58-62
  9. Barrett, R., Berry, M., Chan, T., Demmel, J., Donato, J., Dongarra, J., Eijkhout, V., Pozo, R., Romine, C. and van der Vorst, H. Templates for the Solution of Linear Systems; Building Blocks for Iterative Methods, SIAM, Philadelphia (1994)
  10. Logan, D.L., A First Course in the Finiye Element Method, 2nd Ed., PWS-KENT Publishing Co., Boston (1992)
  11. Ghia, U., Ghia, K.N., and Shin, C. T.. 'High-Re Solutions for Incompressible FIow Using the Navier-Stokes Equations and a Multigrid Method,' J. of Comp. Physics, Vol. 48, (1982), p.387-411
  12. De Vahl Davis, G., 'Natural Convection of Air in a Square Cavity: A Bench Mark Numerical Solution,' Int. J. for Nemerical Methods in Fluids, Vol. 3, (1983), p.249-264