From 6f4908ed000d223b797825b308a5a7ae6bbc8a63 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Mon, 14 Nov 2016 15:03:08 +0000 Subject: [PATCH] TurbulenceModels::ReynoldsStress: remove the spherical part of the wall-function Reynolds stress 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 --- .../ReynoldsStress/ReynoldsStress.C | 13 +++++-------- 1 file changed, 5 insertions(+), 8 deletions(-) diff --git a/src/TurbulenceModels/turbulenceModels/ReynoldsStress/ReynoldsStress.C b/src/TurbulenceModels/turbulenceModels/ReynoldsStress/ReynoldsStress.C index ab0830a20..016407476 100644 --- a/src/TurbulenceModels/turbulenceModels/ReynoldsStress/ReynoldsStress.C +++ b/src/TurbulenceModels/turbulenceModels/ReynoldsStress/ReynoldsStress.C @@ -89,16 +89,13 @@ void Foam::ReynoldsStress::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)); } } }