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

viewFactorsGen: Handle baffles

Browse Source 
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1680
This commit is contained in:
Henry
2015-04-30 22:25:40 +01:00
parent df3e05beb9
commit 9b64e50ca2
1 changed files with 20 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

31
applications/utilities/preProcessing/viewFactorsGen/viewFactorsGen.C
View File

@ -69,6 +69,7 @@ Description
using namespace Foam; using namespace Foam;
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
triSurface triangulate triSurface triangulate
( (
@ -208,19 +209,27 @@ scalar calculateViewFactorFij
{ {
vector r = i - j; vector r = i - j;
scalar rMag = mag(r); scalar rMag = mag(r);
scalar dAiMag = mag(dAi);
scalar dAjMag = mag(dAj);
vector ni = dAi/dAiMag; if (rMag > SMALL)
vector nj = dAj/dAjMag; {
scalar cosThetaJ = mag(nj & r)/rMag; scalar dAiMag = mag(dAi);
scalar cosThetaI = mag(ni & r)/rMag; scalar dAjMag = mag(dAj);
return vector ni = dAi/dAiMag;
( vector nj = dAj/dAjMag;
(cosThetaI*cosThetaJ*dAjMag*dAiMag) scalar cosThetaJ = mag(nj & r)/rMag;
/(sqr(rMag)*constant::mathematical::pi) scalar cosThetaI = mag(ni & r)/rMag;
);
return
(
(cosThetaI*cosThetaJ*dAjMag*dAiMag)
/(sqr(rMag)*constant::mathematical::pi)
);
}
else
{
return 0;
}
} }