From 8ac2a1c7ee674ce23fd888843e7adbd90d7894ac Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Tue, 10 Jun 2014 10:14:53 +0100
Subject: [PATCH] BUG: Adding Qext to the correct position in the matrix
http://www.openfoam.org/mantisbt/view.php?id=1309
---
.../radiationModel/viewFactor/viewFactor.C | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/thermophysicalModels/radiationModels/radiationModel/viewFactor/viewFactor.C b/src/thermophysicalModels/radiationModels/radiationModel/viewFactor/viewFactor.C
index 5cfb6b11c10..c382f3ed788 100644
--- a/src/thermophysicalModels/radiationModels/radiationModel/viewFactor/viewFactor.C
+++ b/src/thermophysicalModels/radiationModels/radiationModel/viewFactor/viewFactor.C
@@ -531,12 +531,12 @@ void Foam::radiation::viewFactor::calculate()
if (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
{
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)
{
- q[i] += (Fmatrix_()[i][j] - 1.0)*sigmaT4;
+ q[i] += (Fmatrix_()[i][j] - 1.0)*sigmaT4 - QrExt[j];
}
else
{
- q[i] += Fmatrix_()[i][j]*sigmaT4 - QrExt[j];
+ q[i] += Fmatrix_()[i][j]*sigmaT4;
}
}
}
--
GitLab