diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/RASProperties b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/RASProperties
deleted file mode 100644
index 379c5ea5..00000000
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/RASProperties
+++ /dev/null
@@ -1,24 +0,0 @@
-/*--------------------------------*- 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    "constant";
-    object      RASProperties;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-RASModel        laminar;    //kEpsilon;
-
-turbulence      off;        //on;
-
-printCoeffs     no;         //yes;
-
-// ************************************************************************* //
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties
index 0801f008..61a80cfa 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties
@@ -87,7 +87,7 @@ momCoupleModels
     implicitCouple
 );
 
-turbulenceModelType "RASProperties";//"LESProperties";//
+turbulenceModelType "turbulenceProperties";//"LESProperties";//
 
 //===========================================================================//
 // sub-model properties
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/polyMesh/boundary b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/polyMesh/boundary
deleted file mode 100644
index 3d85114f..00000000
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/polyMesh/boundary
+++ /dev/null
@@ -1,55 +0,0 @@
-/*--------------------------------*- C++ -*----------------------------------*\
-| =========                 |                                                 |
-| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
-|  \\    /   O peration     | Version:  2.3.x                                 |
-|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
-|    \\/     M anipulation  |                                                 |
-\*---------------------------------------------------------------------------*/
-FoamFile
-{
-    version     2.0;
-    format      ascii;
-    class       polyBoundaryMesh;
-    location    "constant/polyMesh";
-    object      boundary;
-}
-// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-
-5
-(
-    top
-    {
-        type            wall;
-        inGroups        1(wall);
-        nFaces          52;
-        startFace       504;
-    }
-    bottom
-    {
-        type            wall;
-        inGroups        1(wall);
-        nFaces          52;
-        startFace       556;
-    }
-    side-walls
-    {
-        type            wall;
-        inGroups        1(wall);
-        nFaces          104;
-        startFace       608;
-    }
-    inlet
-    {
-        type            patch;
-        nFaces          16;
-        startFace       712;
-    }
-    outlet
-    {
-        type            patch;
-        nFaces          16;
-        startFace       728;
-    }
-)
-
-// ************************************************************************* //
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/turbulenceProperties b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/turbulenceProperties
index 17a36ec3..e8fa087c 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/turbulenceProperties
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/turbulenceProperties
@@ -15,7 +15,28 @@ FoamFile
 }
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
-simulationType  RASModel;
+simulationType  laminar;
+
+ /*laminar
+{
+    RASModel laminar;
+    
+    turbulence off;
+    
+    printCoeffs off;
+} */
+
+/*
+ * simulationType  RAS;
+
+RAS
+{
+    RASModel        kOmegaSST;
+
+    turbulence      on;
+
+    printCoeffs     on;
+} */
 
 
 // ************************************************************************* //
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict
index f9ad5f5a..6e380e13 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict
@@ -53,7 +53,7 @@ maxDeltaT       0.1;
 
 
 // ************************************************************************* //
-libs (  
+/* libs (  
         "libsimpleFunctionObjects.so"
         "libfieldFunctionObjects.so"
      );
@@ -142,7 +142,7 @@ functions
     } 
               
 
-} 
+} */
 
 // ************************************************************************* //
  
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes
index 5a354281..4404d9c3 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes
+++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/fvSchemes
@@ -37,9 +37,10 @@ divSchemes
     div(phi,epsilon)   Gauss upwind;
     div(U)             Gauss linear;
     div(phi,Yi_h)      Gauss upwind;
-    div((muEff*dev2(T(grad(U)))))           Gauss linear;
+    // div((muEff*dev2(T(grad(U)))))           Gauss linear;
     div((viscousTerm*dev2(grad(U).T())))    Gauss linear;
     div((thermo:mu*dev(grad(U).T())))       Gauss linear;
+    div((nuEff*dev2(T(grad(U)))))           Gauss linear;
 }
 
 laplacianSchemes