From 06bbf06c5b489f90ec3e66213acb3e406c09b300 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Thu, 22 Mar 2012 10:56:01 +0000
Subject: [PATCH] BUG: Fix the calculation on max delta normal to the wall on
each cell.It loops over all faces
---
.../SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C | 14 +++++++++-----
1 file changed, 9 insertions(+), 5 deletions(-)
diff --git a/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C b/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C
index 4cd85a7cf60..5a680bc56b1 100644
--- a/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C
+++ b/src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C
@@ -74,16 +74,20 @@ void Foam::IDDESDelta::calcDelta()
{
scalar deltaMaxTmp = 0.0;
const labelList& cFaces = cells[cellI];
- const point& faceCentre = faceCentres[cFaces[0]];
const vector nCell = n[cellI];
forAll(cFaces, cFaceI)
{
label faceI = cFaces[cFaceI];
- const point& faceCentreTwo = faceCentres[faceI];
- scalar tmp = (faceCentre - faceCentreTwo) & nCell;
- if (tmp > deltaMaxTmp)
+ const point& faceCentreI = faceCentres[faceI];
+ forAll(cFaces, cFaceJ)
{
- deltaMaxTmp = tmp;
+ label faceJ = cFaces[cFaceJ];
+ const point& faceCentreJ = faceCentres[faceJ];
+ scalar tmp = (faceCentreJ - faceCentreI) & nCell;
+ if (tmp > deltaMaxTmp)
+ {
+ deltaMaxTmp = tmp;
+ }
}
}
faceToFacenMax[cellI] = deltaMaxTmp;
--
GitLab