A multigrid method for the compressible Navier—Stokes equations coupled to the k—ε turbulence equations

Abstract

The steady compressible Navier—Stokes equations coupled to the k—ε turbulence equations are discretized within a vertex‐centered finite volume formulation. The convective fluxes are obtained by the polynomial flux‐difference splitting upwind method. The first order accurate part results directly from the splitting. The second order part is obtained by the flux‐extrapolation technique using the minmod limiter. The diffusive fluxes are discretized in the central way and are split into a normal and a tangential contribution. The first order accurate part of the convective fluxes together with the normal contribution of the diffusive fluxes form a positive system which allows solution by classical relaxation methods. The source terms in the low‐Reynolds k‐ε equations are grouped into positive and negative terms. The linearized negative source terms are added to the positive system to increase the diagonal dominance. The resulting positive system forms the left hand side of the equations. The remaining terms are put in the right hand side. A multigrid method based on successive relaxation, full weighting, bilinear interpolation and W‐cycle is used. The multigrid method itself acts on the left hand side of the equations. The right hand side is updated in a defect correction cycle.