Description🔗

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.

References🔗

• 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

Usage🔗

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

Boundary file🔗

<patchName>
{
    type            patch;
    ...
}

Field file🔗

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

    // Inherited entries
    ...
}

where:

The inherited entries are elaborated in:

• epsilonWallFunctionFvPatchScalarField.H
• PatchFunction1.H

Further information🔗

Tutorial:

• $FOAM_TUTORIALS/verificationAndValidation/atmosphericModels/atmFlatTerrain

Source code:

• $FOAM_SRC/atmosphericModels/derivedFvPatchFields/wallFunctions/atmEpsilonWallFunction

API:

• Foam::atmEpsilonWallFunctionFvPatchScalarField

History:

• Introduced in version v2006