From a166237d087d9ed36106fa9f0614a0dd27a772e5 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Fri, 17 Mar 2017 11:26:24 +0000
Subject: [PATCH] reactingMultiphaseEulerFoam: Ensure that an unordered phase
 pair gets generated whenever a BlendedInterfacialModel is created.

Resolves bug-report https://bugs.openfoam.org/view.php?id=2472
---
 .../phaseSystem/phaseSystemTemplates.C           | 16 ++++++++++++++++
 1 file changed, 16 insertions(+)

diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseSystem/phaseSystemTemplates.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseSystem/phaseSystemTemplates.C
index 01f36c881..ffa5af34c 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseSystem/phaseSystemTemplates.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phaseSystem/phaseSystemTemplates.C
@@ -131,6 +131,22 @@ void Foam::phaseSystem::generatePairsAndSubModels
                 )
             )
         );
+
+        if (!phasePairs_.found(key))
+        {
+            phasePairs_.insert
+            (
+                key,
+                autoPtr<phasePair>
+                (
+                    new phasePair
+                    (
+                        phaseModels_[key.first()],
+                        phaseModels_[key.second()]
+                    )
+                )
+            );
+        }
     }
 }