diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/U b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/U
index 962f2631..9ac61d7e 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/U
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/U
@@ -17,7 +17,7 @@ FoamFile
 
 dimensions      [0 1 -1 0 0 0 0];
 
-internalField   uniform (0.1525 0 0);
+internalField   uniform (0.0218 0 0);
 
 boundaryField
 {
@@ -40,7 +40,7 @@ boundaryField
     inlet
     {
         type            fixedValue;
-        value           uniform (0.1525 0 0);
+        value           uniform (0.0218 0 0);
     }
 
     outlet
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/transportProperties b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/transportProperties
index 8e0fd139..125570c0 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/transportProperties
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/transportProperties
@@ -17,7 +17,7 @@ FoamFile
 
 transportModel  Newtonian;
 
-nu              nu [ 0 2 -1 0 0 0 0 ] 1.525e-05;    // ḱinematic viscosity of O2 = 1.525e-05; // kinematic viscosity of air (o2+n2) = 1.48e-05;
+nu              nu [ 0 2 -1 0 0 0 0 ] 1.5327e-05;    // ḱinematic viscosity of O2 = 1.525e-05; // kinematic viscosity of air (o2+n2) = 1.48e-05;
 
 kf		kf [ 1 1 -3 -1 0 0 0 ] 0;           // comes from energy model
 
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict
index e2692b8f..02e1df80 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict
@@ -29,7 +29,7 @@ deltaT          0.0005;
 
 writeControl    adjustableRunTime;
 
-writeInterval   0.025;
+writeInterval   0.25;
 
 purgeWrite      0;
 
@@ -47,28 +47,28 @@ runTimeModifiable true;
 
 adjustTimeStep  yes;
 
-maxCo           1.0;
+maxCo           0.5;
 
-maxDeltaT       1;
+maxDeltaT       0.1;
 
 
 // ************************************************************************* //
-libs (
+/*libs (
     "libsampling.so"
     "libutilityFunctionObjects.so"
     "libOpenFOAM.so"
 //"libsimpleSwakFunctionObjects.so"
 //"libswakFunctionObjects.so"
 //"libgroovyBC.so"
-);
-
+);*/
+/*
 functions
 {  
   testing
    {
       type                 writeRegisteredObject;
       functionObjectLibs   ("libIOFunctionObjects.so");
-/*      outputControl        outputTime;*/
+      outputControl        outputTime;
       outputControl        timeStep;
       outputInterval       1;
       objectNames
@@ -113,5 +113,5 @@ functions
         ); 
     } 
 }
-
+*/
 // ************************************************************************* //
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSchemes b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSchemes
index f51c504c..5a354281 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSchemes
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSchemes
@@ -29,8 +29,8 @@ divSchemes
 {
     default         none;
 
-    div(phi,U)         Gauss upwind;
-    div(phid,p)        Gauss upwind;
+    div(phi,U)         Gauss linear; //Gauss upwind;
+    div(phid,p)        Gauss linear; //Gauss upwind;
     div(phi,K)         Gauss linear;
     div(phi,h)         Gauss upwind;
     div(phi,k)         Gauss upwind;
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSolution b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSolution
index 03dd528a..076b7976 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSolution
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/fvSolution
@@ -69,12 +69,12 @@ solvers
         relTol          0;
     }
 
-    "(Yi|CO|CO2|O2)"
+    "(Yi|CO2|O2)"
     {
         $h;
     }
 
-    "(Yi|CO|CO2|O2)Final"
+    "(Yi|CO2|O2)Final"
     {
         $Yi;
         tolerance       1e-06;
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
index 19aa7891..1733e389 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
@@ -43,7 +43,7 @@ fix     zwalls1 all wall/gran model hertz tangential history primitive type 1 zp
 fix     zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.5
 
 # cfd coupling
-fix     cfd all couple/cfd couple_every 1 mpi
+fix     cfd all couple/cfd couple_every 100 mpi
 fix     cfd2 all couple/cfd/force
 
 # this should invoke chemistry
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/shrink.dat b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/shrink.dat
new file mode 100644
index 00000000..a4ade2e2
--- /dev/null
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/shrink.dat
@@ -0,0 +1,203 @@
+# Time-averaged data for fix shrink
+# TimeStep c_radreduce
+0 0.035
+50 0.035
+100 0.035
+150 0.035
+200 0.035
+250 0.035
+300 0.035
+350 0.035
+400 0.035
+450 0.035
+500 0.035
+550 0.035
+600 0.035
+650 0.035
+700 0.035
+750 0.035
+800 0.035
+850 0.035
+900 0.035
+950 0.035
+1000 0.035
+1050 0.035
+1100 0.035
+1150 0.035
+1200 0.035
+1250 0.035
+1300 0.035
+1350 0.035
+1400 0.035
+1450 0.035
+1500 0.035
+1550 0.035
+1600 0.035
+1650 0.035
+1700 0.035
+1750 0.035
+1800 0.035
+1850 0.035
+1900 0.035
+1950 0.035
+2000 0.035
+2050 0.035
+2100 0.035
+2150 0.035
+2200 0.035
+2250 0.035
+2300 0.035
+2350 0.035
+2400 0.035
+2450 0.035
+2500 0.035
+2550 0.035
+2600 0.035
+2650 0.035
+2700 0.035
+2750 0.035
+2800 0.035
+2850 0.035
+2900 0.035
+2950 0.035
+3000 0.035
+3050 0.035
+3100 0.035
+3150 0.035
+3200 0.035
+3250 0.035
+3300 0.035
+3350 0.035
+3400 0.035
+3450 0.035
+3500 0.035
+3550 0.035
+3600 0.035
+3650 0.035
+3700 0.035
+3750 0.035
+3800 0.035
+3850 0.035
+3900 0.035
+3950 0.035
+4000 0.035
+4050 0.035
+4100 0.035
+4150 0.035
+4200 0.035
+4250 0.035
+4300 0.035
+4350 0.035
+4400 0.035
+4450 0.035
+4500 0.035
+4550 0.035
+4600 0.035
+4650 0.035
+4700 0.035
+4750 0.035
+4800 0.035
+4850 0.035
+4900 0.035
+4950 0.035
+5000 0.035
+5050 0.035
+5100 0.035
+5150 0.035
+5200 0.035
+5250 0.035
+5300 0.035
+5350 0.035
+5400 0.035
+5450 0.035
+5500 0.035
+5550 0.035
+5600 0.035
+5650 0.035
+5700 0.035
+5750 0.035
+5800 0.035
+5850 0.035
+5900 0.035
+5950 0.035
+6000 0.035
+6050 0.035
+6100 0.035
+6150 0.035
+6200 0.035
+6250 0.035
+6300 0.035
+6350 0.035
+6400 0.035
+6450 0.035
+6500 0.035
+6550 0.035
+6600 0.035
+6650 0.035
+6700 0.035
+6750 0.035
+6800 0.035
+6850 0.035
+6900 0.035
+6950 0.035
+7000 0.035
+7050 0.035
+7100 0.035
+7150 0.035
+7200 0.035
+7250 0.035
+7300 0.035
+7350 0.035
+7400 0.035
+7450 0.035
+7500 0.035
+7550 0.035
+7600 0.035
+7650 0.035
+7700 0.035
+7750 0.035
+7800 0.035
+7850 0.035
+7900 0.035
+7950 0.035
+8000 0.035
+8050 0.035
+8100 0.035
+8150 0.035
+8200 0.035
+8250 0.035
+8300 0.035
+8350 0.035
+8400 0.035
+8450 0.035
+8500 0.035
+8550 0.035
+8600 0.035
+8650 0.035
+8700 0.035
+8750 0.035
+8800 0.035
+8850 0.035
+8900 0.035
+8950 0.035
+9000 0.035
+9050 0.035
+9100 0.035
+9150 0.035
+9200 0.035
+9250 0.035
+9300 0.035
+9350 0.035
+9400 0.035
+9450 0.035
+9500 0.035
+9550 0.035
+9600 0.035
+9650 0.035
+9700 0.035
+9750 0.035
+9800 0.035
+9850 0.035
+9900 0.035
+9950 0.035
+10000 0.035