From bd5a1942805f7e0d0e01109b905fb4be4527ebfd Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Fri, 14 Dec 2018 11:20:32 +0000 Subject: [PATCH] atmBoundaryLayer: Added documentation --- .../derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C | 4 ++-- .../derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.H | 3 +++ 2 files changed, 5 insertions(+), 2 deletions(-) diff --git a/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C b/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C index 6d2777b020..0e4b23cf49 100644 --- a/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C +++ b/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C @@ -171,7 +171,7 @@ tmp atmBoundaryLayer::k(const vectorField& p) const if (offset_) { - const scalarField z = (zDir_ & p) - zGround_; + const scalarField z((zDir_ & p) - zGround_); tk.ref() = pos0(z)*tk() + neg(z)*kLower_; } @@ -188,7 +188,7 @@ tmp atmBoundaryLayer::epsilon(const vectorField& p) const if (offset_) { - const scalarField z = (zDir_ & p) - zGround_; + const scalarField z((zDir_ & p) - zGround_); tepsilon.ref() = pos0(z)*tepsilon() + neg(z)*epsilonLower_; } diff --git a/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.H b/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.H index 08e3d00e6c..d5d4c01d75 100644 --- a/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.H +++ b/src/atmosphericModels/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.H @@ -82,6 +82,9 @@ Usage Zref | Reference height [m] | yes | z0 | Surface roughness height [m] | yes | zGround | Minimum z-coordinate [m] | yes | + Ulower | Velocity below the BL | no | + kLower | Turbulence k below the BL | no | + epsilonLower | Turbulence epsilon below the BL | no | \endtable Example of the boundary condition specification: