From aadacc5812c7c39be87d98246ac546bb414ba848 Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Wed, 6 Oct 2010 17:41:40 +0100
Subject: [PATCH] BUG: corrected laminar clause in kappt wall function
---
...ppatJayatillekeWallFunctionFvPatchScalarField.C | 14 ++++++--------
1 file changed, 6 insertions(+), 8 deletions(-)
diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/kappatWallFunctions/kappatJayatillekeWallFunction/kappatJayatillekeWallFunctionFvPatchScalarField.C b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/kappatWallFunctions/kappatJayatillekeWallFunction/kappatJayatillekeWallFunctionFvPatchScalarField.C
index 7aa9fa2fe94..6efef7364b0 100644
--- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/kappatWallFunctions/kappatJayatillekeWallFunction/kappatJayatillekeWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/kappatWallFunctions/kappatJayatillekeWallFunction/kappatJayatillekeWallFunctionFvPatchScalarField.C
@@ -231,19 +231,17 @@ void kappatJayatillekeWallFunctionFvPatchScalarField::updateCoeffs()
scalar P = Psmooth(Prat);
scalar yPlusTherm = this->yPlusTherm(P, Prat);
- // Evaluate new effective thermal diffusivity
- scalar kappaEff = 0.0;
- if (yPlus < yPlusTherm)
+ // Update turbulent thermal conductivity
+ if (yPlus > yPlusTherm)
{
- kappaEff = Pr*yPlus;
+ scalar nu = nuw[faceI];
+ scalar kt = nu*(yPlus/(Prt_/kappa_*log(E_*yPlusTherm) + P) - 1/Pr);
+ kappatw[faceI] = max(0.0, kt);
}
else
{
- kappaEff = nuw[faceI]*yPlus/(Prt_/kappa_*log(E_*yPlusTherm) + P);
+ kappatw[faceI] = 0.0;
}
-
- // Update turbulent thermal diffusivity
- kappatw[faceI] = max(0.0, kappaEff - nuw[faceI]/Pr);
}
fixedValueFvPatchField<scalar>::updateCoeffs();
--
GitLab