Development of a scalable finite element solution to the Navier–Stokes equations

A scalable numerical model to solve the unsteady incompressible Navier–Stokes equations is developed using the Galerkin finite element method. The coupled equations are decoupled by the fractional-step method and the systems of equations are inverted by the Krylov subspace iterations. The data structure makes use of a domain decomposition of which each processor stores the parameters in its subdomain, while the linear equations solvers and matrices constructions are parallelized by a data parallel approach. The accuracy of the model is tested by modeling laminar flow inside a two-dimensional square lid-driven cavity for Reynolds numbers at 1,000 as well as three-dimensional turbulent plane and wavy Couette flow and heat transfer at high Reynolds numbers. The parallel performance of the code is assessed by measuring the CPU time taken on an IBM SP2 supercomputer. The speed up factor and parallel efficiency show a satisfactory computational performance.

[1]  T. J. Hanratty,et al.  Interpretation of large-scale structures observed in a turbulent plane Couette flow , 1997 .

[2]  Nils Tillmark,et al.  An investigation of turbulent plane Couette flow at low Reynolds numbers , 1995, Journal of Fluid Mechanics.

[3]  R. Freund,et al.  QMR: a quasi-minimal residual method for non-Hermitian linear systems , 1991 .

[4]  A. Stroud,et al.  Gaussian quadrature formulas , 1966 .

[5]  Chun-Ho Liu,et al.  Turbulent Plane Couette Flow and Scalar Transport at Low Reynolds Number , 2003 .

[6]  A. Peters,et al.  PASTIS-3D — A Parallel Finite Element Projection Code for the Time- Dependent Incompressible Navier-Stokes Equations , 1994 .

[7]  H. J. Leutheusser,et al.  Plane-Couette flow between smooth and rough walls , 1991 .

[8]  Eugenio Oñate,et al.  A finite element formulation for incompressible flow problems using a generalized streamline operator , 1997 .

[9]  Hasan U. Akay,et al.  A domain decomposition based parallel solver for viscous incompressible flows , 1997, Parallel CFD.

[10]  John Kim,et al.  The structure of turbulence in a simulated plane Couette flow , 1991 .

[11]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[12]  A. Ecer,et al.  Parallel Computational Fluid Dynamics: Recent Developments and Advances Using Parallel Computers , 1998 .

[13]  P. Moin,et al.  Effects of the Computational Time Step on Numerical Solutions of Turbulent Flow , 1994 .

[14]  R. Mella,et al.  RAMSYS: A Parallel Code for the Aerodynamic Analysis of 3D Potential Flows around Rotorcraft Configurations , 1998 .

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

[16]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[17]  C. H. Liu,et al.  Development of a finite element solution for the unsteady Navier-Stokes equations using projection method and fractional-θ-scheme , 2001 .

[18]  Parallel computation of atmospheric pollutant dispersion under unstably stratified atmosphere , 1998 .