diff --git a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
index eaa64aebf5..55478e59aa 100644
--- a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
+++ b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
@@ -670,17 +670,20 @@ void Foam::MULES::limitSum(SurfaceScalarFieldList& phiPsiCorrs)
 
     forAll(phiPsiCorrs[0].boundaryField(), patchi)
     {
-        UPtrList<scalarField> phiPsiCorrsPatch(phiPsiCorrs.size());
-        forAll(phiPsiCorrs, phasei)
+        if (phiPsiCorrs[0].boundaryField()[patchi].coupled())
         {
-            phiPsiCorrsPatch.set
-            (
-                phasei,
-                &phiPsiCorrs[phasei].boundaryField()[patchi]
-            );
-        }
+            UPtrList<scalarField> phiPsiCorrsPatch(phiPsiCorrs.size());
+            forAll(phiPsiCorrs, phasei)
+            {
+                phiPsiCorrsPatch.set
+                (
+                    phasei,
+                    &phiPsiCorrs[phasei].boundaryField()[patchi]
+                );
+            }
 
-        limitSum(phiPsiCorrsPatch);
+            limitSum(phiPsiCorrsPatch);
+        }
     }
 }