Premium
Numerical simulation of turbulent free surface flow with two‐equation k –ε eddy‐viscosity models
Author(s) -
Ferreira V. G.,
Mangiavacchi N.,
Tomé M. F.,
Castelo A.,
Cuminato J. A.,
McKee S.
Publication year - 2004
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.641
Subject(s) - turbulence , reynolds number , mechanics , turbulence modeling , boundary layer , mathematics , jet (fluid) , large eddy simulation , flow separation , discretization , reynolds averaged navier–stokes equations , reynolds stress equation model , k epsilon turbulence model , physics , classical mechanics , geometry , k omega turbulence model , mathematical analysis
This paper presents a finite difference technique for solving incompressible turbulent free surface fluid flow problems. The closure of the time‐averaged Navier–Stokes equations is achieved by using the two‐equation eddy‐viscosity model: the high‐Reynolds k –ε (standard) model, with a time scale proposed by Durbin; and a low‐Reynolds number form of the standard k –ε model, similar to that proposed by Yang and Shih. In order to achieve an accurate discretization of the non‐linear terms, a second/third‐order upwinding technique is adopted. The computational method is validated by applying it to the flat plate boundary layer problem and to impinging jet flows. The method is then applied to a turbulent planar jet flow beneath and parallel to a free surface. Computations show that the high‐Reynolds k –ε model yields favourable predictions both of the zero‐pressure‐gradient turbulent boundary layer on a flat plate and jet impingement flows. However, the results using the low‐Reynolds number form of the k –ε model are somewhat unsatisfactory. Copyright © 2004 John Wiley & Sons, Ltd.