zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: Fix the calculation on max delta normal to the wall on each cell.It

Browse Source 
loops over all faces
This commit is contained in:
sergio
2012-03-22 10:56:01 +00:00
parent d2e63f4293
commit 06bbf06c5b
1 changed files with 9 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
src/turbulenceModels/incompressible/LES/SpalartAllmarasIDDES/IDDESDelta/IDDESDelta.C
View File

@ -74,16 +74,20 @@ void Foam::IDDESDelta::calcDelta()
{ {
scalar deltaMaxTmp = 0.0; scalar deltaMaxTmp = 0.0;
const labelList& cFaces = cells[cellI]; const labelList& cFaces = cells[cellI];
const point& faceCentre = faceCentres[cFaces[0]];
const vector nCell = n[cellI]; const vector nCell = n[cellI];
forAll(cFaces, cFaceI) forAll(cFaces, cFaceI)
{ {
label faceI = cFaces[cFaceI]; label faceI = cFaces[cFaceI];
const point& faceCentreTwo = faceCentres[faceI]; const point& faceCentreI = faceCentres[faceI];
scalar tmp = (faceCentre - faceCentreTwo) & nCell; forAll(cFaces, cFaceJ)
if (tmp > deltaMaxTmp)
{ {
deltaMaxTmp = tmp; label faceJ = cFaces[cFaceJ];
const point& faceCentreJ = faceCentres[faceJ];
scalar tmp = (faceCentreJ - faceCentreI) & nCell;
if (tmp > deltaMaxTmp)
{
deltaMaxTmp = tmp;
}
} }
} }
faceToFacenMax[cellI] = deltaMaxTmp; faceToFacenMax[cellI] = deltaMaxTmp;