Skip to content
Snippets Groups Projects
Commit 8ac2a1c7 authored by sergio's avatar sergio Committed by Andrew Heather
Browse files

BUG: Adding Qext to the correct position in the matrix

parent 01972cf5
No related branches found
No related tags found
No related merge requests found
...@@ -531,12 +531,12 @@ void Foam::radiation::viewFactor::calculate() ...@@ -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() ...@@ -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;
} }
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment