diff --git a/dumux/freeflow/rans/twoeq/kepsilon/problem.hh b/dumux/freeflow/rans/twoeq/kepsilon/problem.hh
index 604730549df2f7128b129796b3d263a738a8f234..1c236bd050be44c2633016939869c7316b2f2ea2 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]);
     }
 
     /*!