From 8b90aa522000cf57db5e1d21c94a67ca9dd93862 Mon Sep 17 00:00:00 2001 From: Sergio Ferraris Date: Wed, 23 Oct 2013 10:43:41 +0100 Subject: [PATCH] ENH: Adding FPE stabilization to Edash for nutkAtmRoughWall --- .../nutkAtmRoughWallFunctionFvPatchScalarField.C | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkAtmRoughWallFunction/nutkAtmRoughWallFunctionFvPatchScalarField.C b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkAtmRoughWallFunction/nutkAtmRoughWallFunctionFvPatchScalarField.C index 2406f30d0a..9e724596bf 100644 --- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkAtmRoughWallFunction/nutkAtmRoughWallFunctionFvPatchScalarField.C +++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkAtmRoughWallFunction/nutkAtmRoughWallFunctionFvPatchScalarField.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -63,10 +63,10 @@ tmp nutkAtmRoughWallFunctionFvPatchScalarField::calcNut() const scalar uStar = Cmu25*sqrt(k[faceCellI]); scalar yPlus = uStar*y[faceI]/nuw[faceI]; - scalar Edash = (y[faceI] + z0_[faceI])/z0_[faceI]; + scalar Edash = (y[faceI] + z0_[faceI])/(z0_[faceI] + 1e-4); nutw[faceI] = - nuw[faceI]*(yPlus*kappa_/log(max(Edash, 1+1e-4)) - 1); + nuw[faceI]*(yPlus*kappa_/log(max(Edash, 1 + 1e-4)) - 1); if (debug) {