The solution of the incompressible Navier-Stokes equations in general two-and three-dimensional domains using a multigrid method is considered. Because a great variety of boundary-fitted grids may occur, robustness is at a premium. Therefore a new ILU smoother called CILU (collective ILU) is describ
A multigrid method for an invariant formulation of the incompressible navier-stokes equations in general Co-ordinates
β Scribed by Oosterlee, C. W. ;Wesseling, P.
- Publisher
- Wiley (John Wiley & Sons)
- Year
- 1992
- Tongue
- English
- Weight
- 621 KB
- Volume
- 8
- Category
- Article
- ISSN
- 0748-8025
No coin nor oath required. For personal study only.
β¦ Synopsis
The stationary incompressible Navier-Stokes equations are discretized with a finite volume method in curvilinear co-ordinates. The arbitrarily shaped domain is mapped onto a rectangular block, resulting in a boundary-fitted grid. In order to obtain accurate discretizations of the transformed equations, some requirements on geometric quantities should be met. The choice of velocity components is also of importance, Contravariant flux unknowns and pressure p are used as primary unknowns on a staggered grid arrangement.
The system of discretized equations is solved with a non-linear multigrid algorithm, into which a smoother, called Symmetric Coupled Gauss-Seidel. is implemented. Cell by cell, all unknowns in the grid cell are updated by solving four momentum equations and a continuity equation simultaneously. The solution algorithm shows satisfying average reduction factors for several domains.
π SIMILAR VOLUMES
This paper describes a domain decomposition method for the incompressible Navier -Stokes equations in general co-ordinates. Domain decomposition techniques are needed for solving flow problems in complicated geometries while retaining structured grids on each of the subdomains. This is the so-called
The article presents a fast pseudo-spectral Navier-Stokes solver for cylindrical geometries, which is shown to possess exponential rate of decay of the error. The formulation overcomes the issues related to the axis singularity, by employing in the radial direction a special set of collocation point