exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Gauss grading bug fix, Zeljko Tukovic

Browse source
This commit is contained in:
Hrvoje Jasak 2012-04-23 15:44:42 +01:00
parent 190c4f8eae
commit 523f2eb35a
1 changed files with 2 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/finiteArea/finiteArea/gradSchemes/gaussFaGrad/gaussFaGrad.C
View file

@ -22,8 +22,6 @@ License
along with OpenFOAM; if not, write to the Free Software Foundation, along with OpenFOAM; if not, write to the Free Software Foundation,
Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
Description
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "gaussFaGrad.H" #include "gaussFaGrad.H"
@ -69,12 +67,12 @@ gaussGrad<Type>::grad
GeometricField<GradType, faPatchField, areaMesh>& gGrad = tgGrad(); GeometricField<GradType, faPatchField, areaMesh>& gGrad = tgGrad();
gGrad -= vsf*fac::edgeIntegrate(vsf.mesh().Le());
// Remove component of gradient normal to surface (area) // Remove component of gradient normal to surface (area)
const areaVectorField& n = vsf.mesh().faceAreaNormals(); const areaVectorField& n = vsf.mesh().faceAreaNormals();
// const areaScalarField& K = vsf.mesh().faceCurvatures();
gGrad -= n*(n & gGrad); gGrad -= n*(n & gGrad);
// gGrad -= K*(n*vsf);
gGrad.correctBoundaryConditions(); gGrad.correctBoundaryConditions();
gGrad.rename("grad(" + vsf.name() + ')'); gGrad.rename("grad(" + vsf.name() + ')');