research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

atmBoundaryLayer: Removed the hard-coded limit on z0

Browse source 
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1638
This commit is contained in:
Henry 2015-03-27 15:20:18 +00:00
parent c08a8c97a0
commit 1449a67364
3 changed files with 5 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayerInletEpsilon/atmBoundaryLayerInletEpsilonFvPatchScalarField.C
View file

@ -110,7 +110,7 @@ atmBoundaryLayerInletEpsilonFvPatchScalarField
{ {
Ustar_[i] = Ustar_[i] =
kappa_*Uref_ kappa_*Uref_
/(log((Href_ + z0_[i])/max(z0_[i] , scalar(0.001)))); /(log((Href_ + z0_[i])/z0_[i]));
} }
z_ /= mag(z_); z_ /= mag(z_);

8
src/TurbulenceModels/turbulenceModels/RAS/derivedFvPatchFields/atmBoundaryLayerInletVelocity/atmBoundaryLayerInletVelocityFvPatchVectorField.C
View file

@ -116,7 +116,7 @@ atmBoundaryLayerInletVelocityFvPatchVectorField
{ {
Ustar_[i] = Ustar_[i] =
kappa_*Uref_ kappa_*Uref_
/(log((Href_ + z0_[i])/max(z0_[i] , scalar(0.001)))); /(log((Href_ + z0_[i])/z0_[i]));
} }
const vectorField& c = patch().Cf(); const vectorField& c = patch().Cf();
@ -129,11 +129,7 @@ atmBoundaryLayerInletVelocityFvPatchVectorField
{ {
Un[i] = Un[i] =
(Ustar_[i]/kappa_) (Ustar_[i]/kappa_)
*log *log((coord[i] - zGround_[i] + z0_[i])/z0_[i]);
(
(coord[i] - zGround_[i] + z0_[i])
/max(z0_[i], scalar(0.001))
);
} }
else else
{ {

4
src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/atmBoundaryLayer/atmBoundaryLayer.C
View file

@ -73,7 +73,7 @@ atmBoundaryLayer::atmBoundaryLayer(const vectorField& p, const dictionary& dict)
flowDir_ /= mag(flowDir_); flowDir_ /= mag(flowDir_);
zDir_ /= mag(zDir_); zDir_ /= mag(zDir_);
Ustar_ = kappa_*Uref_/(log((Zref_ + z0_)/max(z0_, scalar(0.001)))); Ustar_ = kappa_*Uref_/(log((Zref_ + z0_)/z0_));
} }
@ -136,7 +136,7 @@ tmp<vectorField> atmBoundaryLayer::U(const vectorField& p) const
scalarField Un scalarField Un
( (
(Ustar_/kappa_) (Ustar_/kappa_)
*log(((zDir_ & p) - zGround_ + z0_)/max(z0_, scalar(0.001))) *log(((zDir_ & p) - zGround_ + z0_)/z0_)
); );
return flowDir_*Un; return flowDir_*Un;