diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2 b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2
index ff10ed20..14c11b22 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2
@@ -38,7 +38,8 @@ boundaryField
 
     inlet
     {
-        type            zeroGradient;
+        type            fixedValue;
+        value           uniform 0.0;
     }
 
     outlet
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2 b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2
index 4f1a34f6..d4c9aabd 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2
@@ -38,7 +38,8 @@ boundaryField
 
     inlet
     {
-	type		zeroGradient;
+        type            fixedValue;
+        value           uniform 1.0;
     }
 
     outlet
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/T b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/T
index dc46d191..18f141a1 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/T
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/T
@@ -38,7 +38,8 @@ boundaryField
 
     inlet
     {
-        type            zeroGradient;
+        type            fixedValue;
+        value           uniform 1293.15;
     }
 
     outlet
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/rho b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/rho
index c2e2d9f1..14ee0ebb 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/rho
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/rho
@@ -37,7 +37,8 @@ boundaryField
 
     inlet
     {
-        type            zeroGradient;
+        type            fixedValue;
+        value           uniform 0.38;
     }
 
     outlet
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvOptions b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvOptions
new file mode 100644
index 00000000..1b6c818a
--- /dev/null
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvOptions
@@ -0,0 +1,28 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  2.3.0                                 |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "system";
+    object      fvOptions;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+source1
+{
+    type        limitTemperature;
+    active          yes;
+    limitTemperatureCoeffs
+    {
+        active          yes;
+        selectionMode   all;
+        Tmin     500;
+        Tmax     1500;
+    }
+}
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes
index 4404d9c3..204ce281 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes
@@ -29,7 +29,7 @@ divSchemes
 {
     default         none;
 
-    div(phi,U)         Gauss linear; //Gauss upwind;
+    div(phi,U)         Gauss upwind;
     div(phid,p)        Gauss linear; //Gauss upwind;
     div(phi,K)         Gauss linear;
     div(phi,h)         Gauss upwind;
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSolution
index 9832d2cc..a21941dd 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSolution
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSolution
@@ -50,7 +50,7 @@ solvers
     pFinal
     {
         $p;
-        tolerance       1e-06;
+        tolerance       1e-07;
         relTol          0;
     }
 
@@ -65,7 +65,7 @@ solvers
     "(U|h|R|k|epsilon)Final"
     {
         $U;
-        tolerance       1e-05;
+        tolerance       1e-07;
         relTol          0;
     }
 
@@ -77,7 +77,7 @@ solvers
     "(Yi|CO2|O2)Final"
     {
         $Yi;
-        tolerance       1e-06;
+        tolerance       1e-08;
         relTol          0;
     }
 }