From 923f4023bbe4326450cb8ec6cbe74ed53b549660 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 4 Feb 2015 22:11:15 +0000
Subject: [PATCH] kLowReWallFunction: Limit k to avoid failure of the
turbulence model due to division by k
---
.../kLowReWallFunction/kLowReWallFunctionFvPatchScalarField.C | 3 +++
1 file changed, 3 insertions(+)
diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/kqRWallFunctions/kLowReWallFunction/kLowReWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/kqRWallFunctions/kLowReWallFunction/kLowReWallFunctionFvPatchScalarField.C
index 8749bc02a41..29a6f0e2b41 100644
--- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/kqRWallFunctions/kLowReWallFunction/kLowReWallFunctionFvPatchScalarField.C
+++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/kqRWallFunctions/kLowReWallFunction/kLowReWallFunctionFvPatchScalarField.C
@@ -215,6 +215,9 @@ void kLowReWallFunctionFvPatchScalarField::updateCoeffs()
kw[faceI] *= sqr(uTau);
}
+ // Limit kw to avoid failure of the turbulence model due to division by kw
+ kw = max(kw, SMALL);
+
fixedValueFvPatchField<scalar>::updateCoeffs();
// TODO: perform averaging for cells sharing more than one boundary face
--
GitLab