Premium
Assessment of the finite volume method applied to the v 2 − f model
Author(s) -
Mortensen Mikael,
Reif Bjørn Anders Pettersson,
Wasberg Carl Erik
Publication year - 2009
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.2091
Subject(s) - discretization , lagrange polynomial , interpolation (computer graphics) , finite volume method , quadrature (astronomy) , mathematics , gaussian quadrature , order of accuracy , finite element method , turbulence , mathematical analysis , solver , numerical analysis , physics , numerical stability , mathematical optimization , boundary value problem , classical mechanics , mechanics , motion (physics) , nyström method , polynomial , optics , thermodynamics
The objective of this paper is to assess the accuracy of low‐order finite volume (FV) methods applied to the v 2 − f turbulence model of Durbin ( Theoret. Comput. Fluid Dyn. 1991; 3 :1–13) in the near vicinity of solid walls. We are not (like many others) concerned with the stability of solvers ‐ the topic at hand is simply whether the mathematical properties of the v 2 − f model can be captured by the given, widespread, numerical method. The v 2 − f model is integrated all the way up to solid walls, where steep gradients in turbulence parameters are observed. The full resolution of wall gradients imposes quite high demands on the numerical schemes and it is not evident that common (second order) FV codes can fully cope with such demands. The v 2 − f model is studied in a statistically one‐dimensional, fully developed channel flow where we compare FV schemes with a highly accurate spectral element reference implementation. For the FV method a higher‐order face interpolation scheme, using Lagrange interpolation polynomials up to arbitrary order, is described. It is concluded that a regular second‐order FV scheme cannot give an accurate representation of all model parameters, independent of mesh density. To match the spectral element solution an extended source treatment (we use three‐point Gauss–Lobatto quadrature), as well as a higher‐order discretization of diffusion is required. Furthermore, it is found that the location of the first internal node need to be well within y + =1. Copyright © 2009 John Wiley & Sons, Ltd.