research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

BUG: Adding Qext to the correct position in the matrix

Browse source 
http://www.openfoam.org/mantisbt/view.php?id=1309
This commit is contained in:
sergio 2014-06-10 10:14:53 +01:00
parent 4d7c1cc2c0
commit cdd3e495bb
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/thermophysicalModels/radiationModels/radiationModel/viewFactor/viewFactor.C
View file

@ -531,12 +531,12 @@ void Foam::radiation::viewFactor::calculate()
if (i==j) if (i==j)
{ {
C[i][j] = invEj - (invEj - 1.0)*Fmatrix_()[i][j]; C[i][j] = invEj - (invEj - 1.0)*Fmatrix_()[i][j];
q[i] += (Fmatrix_()[i][j] - 1.0)*sigmaT4; q[i] += (Fmatrix_()[i][j] - 1.0)*sigmaT4 - QrExt[j];
} }
else else
{ {
C[i][j] = (1.0 - invEj)*Fmatrix_()[i][j]; C[i][j] = (1.0 - invEj)*Fmatrix_()[i][j];
q[i] += Fmatrix_()[i][j]*sigmaT4 - QrExt[j]; q[i] += Fmatrix_()[i][j]*sigmaT4;
} }
} }
@ -580,11 +580,11 @@ void Foam::radiation::viewFactor::calculate()
if (i==j) if (i==j)
{ {
q[i] += (Fmatrix_()[i][j] - 1.0)*sigmaT4; q[i] += (Fmatrix_()[i][j] - 1.0)*sigmaT4 - QrExt[j];
} }
else else
{ {
q[i] += Fmatrix_()[i][j]*sigmaT4 - QrExt[j]; q[i] += Fmatrix_()[i][j]*sigmaT4;
} }
} }
} }