From 5550139aec88d3eb4e7e964a2fb40eed1db5cebb Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Fri, 27 Mar 2015 15:15:16 +0000
Subject: [PATCH] atmBoundaryLayer: Removed the hard-coded limit on z0 Resolves
bug-report http://www.openfoam.org/mantisbt/view.php?id=1638
---
.../derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C b/src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
index bc666d79cef..5cd044eb341 100644
--- a/src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
+++ b/src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
@@ -71,7 +71,7 @@ atmBoundaryLayer::atmBoundaryLayer(const vectorField& p, const dictionary& dict)
flowDir_ /= mag(flowDir_);
zDir_ /= mag(zDir_);
- Ustar_ = kappa_*Uref_/(log((Zref_ + z0_)/max(z0_, scalar(0.001))));
+ Ustar_ = kappa_*Uref_/(log((Zref_ + z0_)/z0_));
}
@@ -134,7 +134,7 @@ tmp<vectorField> atmBoundaryLayer::U(const vectorField& p) const
scalarField Un
(
(Ustar_/kappa_)
- *log(((zDir_ & p) - zGround_ + z0_)/max(z0_, scalar(0.001)))
+ *log(((zDir_ & p) - zGround_ + z0_)/z0_)
);
return flowDir_*Un;
--
GitLab