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

multiphaseEuler: Fix indexing bug in partial elimination

Browse Source
This commit is contained in:
Will Bainbridge
2023-08-16 16:03:32 +01:00
parent 7ab2321a07
commit e3b53f2727
1 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
View File

@ -1442,8 +1442,8 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::partialElimination
{ {
for (label j = i + 1; j < phases.size(); j++) for (label j = i + 1; j < phases.size(); j++)
{ {
KdByAs[i][j] /= KdByAs[i][i]; KdByAs[j][i] /= KdByAs[i][i];
KdByAfs[i][j] /= KdByAfs[i][i]; KdByAfs[j][i] /= KdByAfs[i][i];
for (label k = i + 1; k < phases.size(); ++ k) for (label k = i + 1; k < phases.size(); ++ k)
{ {
KdByAs[j][k] -= KdByAs[j][i]*KdByAs[i][k]; KdByAs[j][k] -= KdByAs[j][i]*KdByAs[i][k];
@ -1584,7 +1584,7 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::partialEliminationf
{ {
for (label j = i + 1; j < phases.size(); j++) for (label j = i + 1; j < phases.size(); j++)
{ {
phiKdfs[i][j] /= phiKdfs[i][i]; phiKdfs[j][i] /= phiKdfs[i][i];
for (label k = i + 1; k < phases.size(); ++ k) for (label k = i + 1; k < phases.size(); ++ k)
{ {
phiKdfs[j][k] -= phiKdfs[j][i]*phiKdfs[i][k]; phiKdfs[j][k] -= phiKdfs[j][i]*phiKdfs[i][k];