zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

TurbulenceModels::ReynoldsStress: remove the spherical part of the wall-function Reynolds stress

Browse Source 
The spherical part of the Reynolds stress is included in the pressure so
that the wall boundary condition for the pressure is zeroGradient.

Resolves bug-report http://bugs.openfoam.org/view.php?id=2325
This commit is contained in:
Henry Weller
2016-11-14 15:03:08 +00:00
parent e540e333f0
commit 6f4908ed00
1 changed files with 5 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/TurbulenceModels/turbulenceModels/ReynoldsStress/ReynoldsStress.C
View File

@ -89,16 +89,13 @@ void Foam::ReynoldsStress<BasicTurbulenceModel>::correctWallShearStress
forAll(curPatch, facei)
{
// Calculate near-wall velocity gradient
tensor gradUw
const tensor gradUw
= (faceAreas[facei]/magFaceAreas[facei])*snGradU[facei];
// Calculate near-wall shear-stress tensor
tensor tauw = -nutw[facei]*2*dev(symm(gradUw));
// Reset the shear components of the stress tensor
Rw[facei].xy() = tauw.xy();
Rw[facei].xz() = tauw.xz();
Rw[facei].yz() = tauw.yz();
// Set the wall Reynolds-stress to the near-wall shear-stress
// Note: the spherical part of the normal stress is included in
// the pressure
Rw[facei] = -nutw[facei]*2*dev(symm(gradUw));
}
}
}