diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phasePair/phasePair/phasePair.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phasePair/phasePair/phasePair.C
index 77d113817..701e1eee8 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phasePair/phasePair/phasePair.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/phasePair/phasePair/phasePair.C
@@ -25,6 +25,7 @@ License
 
 #include "phasePair.H"
 #include "phaseSystem.H"
+#include "surfaceTensionModel.H"
 
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
@@ -165,7 +166,11 @@ Foam::tmp<Foam::volScalarField> Foam::phasePair::EoH2() const
 
 Foam::tmp<Foam::volScalarField> Foam::phasePair::sigma() const
 {
-    return phase1().fluid().sigma(phasePair(phase1(), phase2()));
+    return
+        phase1().fluid().lookupSubModel<surfaceTensionModel>
+        (
+            phasePair(phase1(), phase2())
+        ).sigma();
 }
 
 
diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C
index 6d269cd03..8cbe0d25a 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C
@@ -30,6 +30,7 @@ License
 #include "driftModel.H"
 #include "nucleationModel.H"
 #include "phaseSystem.H"
+#include "surfaceTensionModel.H"
 #include "fvmDdt.H"
 #include "fvcDdt.H"
 #include "fvmSup.H"
@@ -1152,13 +1153,11 @@ Foam::diameterModels::populationBalanceModel::sigmaWithContinuousPhase
     const phaseModel& dispersedPhase
 ) const
 {
-    const phasePairKey key
-    (
-        dispersedPhase.name(),
-        continuousPhase_.name()
-    );
-
-    return fluid_.sigma(key);
+    return
+        fluid_.lookupSubModel<surfaceTensionModel>
+        (
+            phasePair(dispersedPhase, continuousPhase_)
+        ).sigma();
 }