The atmEpsilonWallFunction is a boundary condition that provides a wall constraint on the turbulent kinetic energy dissipation rate (i.e. epsilon) and the turbulent kinetic energy production contribution (i.e. G) for atmospheric boundary layer modelling.


  • Theoretical expressions (tags:PGVB, RH):

    Parente, A., GorlΓ©, C., Van Beeck, J., & Benocci, C. (2011). Improved k–Ρ model and wall function formulation for the RANS simulation of ABL flows. J. of wind engineering and industrial aerodynamics, 99(4), 267-278. DOI:10.1016/j.jweia.2010.12.017

    Richards, P. J., & Hoxey, R. P. (1993). Appropriate boundary conditions for computational wind engineering models using the k-Ξ΅ turbulence model. In Computational Wind Engineering 1 (pp. 145-153). DOI:10.1016/B978-0-444-81688-7.50018-8


The condition requires entries in both the boundary and field files.

Boundary file

    type            patch;

Field file

    // Mandatory entries
    type            atmEpsilonWallFunction;
    z0              <PatchFunction1<scalar>>;

    // Inherited entries


Property Description Type Required Default
type Type name: atmEpsilonWallFunction word yes -
z0 Surface roughness length [m] PatchFunction1<scalar> yes -

The inherited entries are elaborated in:

  • epsilonWallFunctionFvPatchScalarField.H
  • PatchFunction1.H

Further information


Source code:



  • Introduced in version v2006