From 84569f0f12c48a259b205afb5372e9f59d78b9fe Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Tue, 11 Jun 2013 10:43:45 +0100
Subject: [PATCH] BUG: Updated accumulations in wall functions
---
.../epsilonWallFunctionFvPatchScalarField.C | 4 ++--
.../omegaWallFunction/omegaWallFunctionFvPatchScalarField.C | 4 ++--
.../omegaWallFunction/omegaWallFunctionFvPatchScalarField.C | 4 ++--
3 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
index 5e78e48c973..1b3e4ce3d14 100644
--- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
@@ -237,9 +237,9 @@ void epsilonWallFunctionFvPatchScalarField::calculate
scalar w = cornerWeights[faceI];
- epsilon[cellI] = w*Cmu75*pow(k[cellI], 1.5)/(kappa_*y[faceI]);
+ epsilon[cellI] += w*Cmu75*pow(k[cellI], 1.5)/(kappa_*y[faceI]);
- G[cellI] =
+ G[cellI] +=
w
*(mutw[faceI] + muw[faceI])
*magGradUw[faceI]
diff --git a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
index fff165e15e3..929ca245982 100644
--- a/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
@@ -244,9 +244,9 @@ void omegaWallFunctionFvPatchScalarField::calculate
scalar omegaLog = sqrt(k[cellI])/(Cmu25*kappa_*y[faceI]);
- omega[cellI] = w*sqrt(sqr(omegaVis) + sqr(omegaLog));
+ omega[cellI] += w*sqrt(sqr(omegaVis) + sqr(omegaLog));
- G[cellI] =
+ G[cellI] +=
w
*(mutw[faceI] + muw[faceI])
*magGradUw[faceI]
diff --git a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
index 554c18542d9..6b8504a4e9b 100644
--- a/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
+++ b/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
@@ -242,9 +242,9 @@ void omegaWallFunctionFvPatchScalarField::calculate
scalar omegaLog = sqrt(k[cellI])/(Cmu25*kappa_*y[faceI]);
- omega[cellI] = w*sqrt(sqr(omegaVis) + sqr(omegaLog));
+ omega[cellI] += w*sqrt(sqr(omegaVis) + sqr(omegaLog));
- G[cellI] =
+ G[cellI] +=
w
*(nutw[faceI] + nuw[faceI])
*magGradUw[faceI]
--
GitLab