From df43c8cb1de8047ce5a3bb265dc5f2f7622a7c52 Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Thu, 28 Jun 2018 18:02:38 +0200
Subject: [PATCH] [kepsilon] Do not limit uStarNominal
---
dumux/freeflow/rans/twoeq/kepsilon/problem.hh | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/dumux/freeflow/rans/twoeq/kepsilon/problem.hh b/dumux/freeflow/rans/twoeq/kepsilon/problem.hh
index 604730549d..1c236bd050 100644
--- a/dumux/freeflow/rans/twoeq/kepsilon/problem.hh
+++ b/dumux/freeflow/rans/twoeq/kepsilon/problem.hh
@@ -254,12 +254,11 @@ public:
//! \brief Returns the nominal wall shear stress velocity (accounts for poor approximation of viscous sublayer)
const Scalar uStarNominal(unsigned int elementID) const
{
- using std::max;
using std::pow;
using std::sqrt;
unsigned int matchingPointID = matchingPointID_[asImp_().wallElementID_[elementID]];
return pow(cMu(), 0.25)
- * sqrt(max(storedTurbulentKineticEnergy_[matchingPointID],1e-8));
+ * sqrt(storedTurbulentKineticEnergy_[matchingPointID]);
}
/*!
--
GitLab