diff --git a/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C b/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C
index 1ce0a83871..54df8c5bad 100644
--- a/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C
+++ b/src/finiteVolume/cfdTools/general/porosityModel/DarcyForchheimer/DarcyForchheimer.C
@@ -203,9 +203,21 @@ void Foam::porosityModels::DarcyForchheimer::correct
     if (UEqn.dimensions() == dimForce)
     {
         const volScalarField& rho = mesh_.lookupObject<volScalarField>(rhoName);
-        const volScalarField& mu = mesh_.lookupObject<volScalarField>(muName);
 
-        apply(Udiag, Usource, V, rho, mu, U);
+        if (mesh_.foundObject<volScalarField>(muName))
+        {
+            const volScalarField& mu =
+                mesh_.lookupObject<volScalarField>(muName);
+
+            apply(Udiag, Usource, V, rho, mu, U);
+        }
+        else
+        {
+            const volScalarField& nu =
+                mesh_.lookupObject<volScalarField>(nuName);
+
+            apply(Udiag, Usource, V, rho, rho*nu, U);
+        }
     }
     else
     {