diff --git a/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C b/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
index a732ec4be3..8450bd912b 100644
--- a/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
+++ b/applications/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
@@ -1442,8 +1442,8 @@ void Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::partialElimination
         {
             for (label j = i + 1; j < phases.size(); j++)
             {
-                KdByAs[i][j] /= KdByAs[i][i];
-                KdByAfs[i][j] /= KdByAfs[i][i];
+                KdByAs[j][i] /= KdByAs[i][i];
+                KdByAfs[j][i] /= KdByAfs[i][i];
                 for (label k = i + 1; k < phases.size(); ++ 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++)
             {
-                phiKdfs[i][j] /= phiKdfs[i][i];
+                phiKdfs[j][i] /= phiKdfs[i][i];
                 for (label k = i + 1; k < phases.size(); ++ k)
                 {
                     phiKdfs[j][k] -= phiKdfs[j][i]*phiKdfs[i][k];