diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.C b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.C
index b71bb4f8d1c5134a8b08b2579c68d0d42308583c..c51a7d7e8cd84f921e194bc9ad081f68702a4b9b 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -61,7 +61,8 @@ JohnsonJackson
     Fr_("Fr", dimensionSet(1, -1, -2, 0, 0), coeffDict_.lookup("Fr")),
     eta_("eta", dimless, coeffDict_.lookup("eta")),
     p_("p", dimless, coeffDict_.lookup("p")),
-    phi_("phi", dimless, coeffDict_.lookup("phi"))
+    phi_("phi", dimless, coeffDict_.lookup("phi")),
+    alphaDeltaMax_("alphaDeltaMax", dimless, coeffDict_.lookup("alphaDeltaMax"))
 {
     phi_ *= constant::mathematical::pi/180.0;
 }
@@ -88,7 +89,7 @@ frictionalPressure
 
     return
         Fr_*pow(max(alpha1 - alphaMinFriction, scalar(0)), eta_)
-       /pow(max(alphaMax - alpha1, scalar(5.0e-2)), p_);
+       /pow(max(alphaMax - alpha1, alphaDeltaMax_), p_);
 }
 
 
@@ -106,7 +107,7 @@ frictionalPressurePrime
         eta_*pow(max(alpha1 - alphaMinFriction, scalar(0)), eta_ - 1.0)
        *(alphaMax-alpha1)
       + p_*pow(max(alpha1 - alphaMinFriction, scalar(0)), eta_)
-    )/pow(max(alphaMax - alpha1, scalar(5.0e-2)), p_ + 1.0);
+    )/pow(max(alphaMax - alpha1, alphaDeltaMax_), p_ + 1.0);
 }
 
 
@@ -134,6 +135,8 @@ bool Foam::kineticTheoryModels::frictionalStressModels::JohnsonJackson::read()
     phi_.read(coeffDict_);
     phi_ *= constant::mathematical::pi/180.0;
 
+    alphaDeltaMax_.read(coeffDict_);
+
     return true;
 }
 
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.H b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.H
index 0e5f7c9a84d4923456d031d3fbe2a0fb5acceb59..1d42177bf0532559ccb23098d9121fadb622d476 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/phaseCompressibleTurbulenceModels/kineticTheoryModels/frictionalStressModel/JohnsonJackson/JohnsonJacksonFrictionalStress.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2014 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -69,6 +69,9 @@ class JohnsonJackson
         //- Angle of internal friction
         dimensionedScalar phi_;
 
+        //- Limit of alphaMax - alpha1 (~5e-2 or higher)
+        dimensionedScalar alphaDeltaMax_;
+
 
 public:
 
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/turbulenceProperties.particles b/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/turbulenceProperties.particles
index 5ade08c60002f89cbfbe632cb0469ae739076dff..7d74e9d4baed2f833e266f569db49d7aa353eadb 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/turbulenceProperties.particles
+++ b/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/turbulenceProperties.particles
@@ -45,6 +45,7 @@ RAS
             eta                     2;
             p                       5;
             phi                     28.5;
+            alphaDeltaMax           0.05;
         }
     }
 
diff --git a/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/turbulenceProperties.particles b/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/turbulenceProperties.particles
index 39b5549d3234fd6884b88f5a04e256edf3df950d..d36ee69392bc156830cf934fd1f0d9490dd968fd 100644
--- a/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/turbulenceProperties.particles
+++ b/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/turbulenceProperties.particles
@@ -44,6 +44,7 @@ RAS
             eta                     2;
             p                       5;
             phi                     28.5;
+            alphaDeltaMax           0.05;
         }
     }