diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/couplingProperties b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/couplingProperties
index 80e777cf..ebfa060a 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/couplingProperties
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/couplingProperties
@@ -31,7 +31,7 @@ FoamFile
 
 modelType "A"; // A or B
 
-couplingInterval 10;//1000;
+couplingInterval 100;//1000;
 
 voidFractionModel divided;//centre;//
 
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict
index d519b0ca..eee46edd 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/system/controlDict
@@ -25,11 +25,11 @@ stopAt          endTime;
 
 endTime         0.1;
 
-deltaT          0.001;
+deltaT          0.0005;
 
 writeControl    adjustableRunTime;
 
-writeInterval   0.0025;
+writeInterval   0.001;
 
 purgeWrite      0;
 
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init
index 3d63a5c9..dd86f1fd 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init
@@ -14,7 +14,7 @@ processors      2 2 1
 region          reg block 0. 1.5 0. 0.5 0. 0.5 units box
 create_box      1 reg
 
-neighbor        0.001 bin
+neighbor        0.0005 bin
 neigh_modify    delay 0 
 
 # Material properties required for granular pair styles
@@ -28,7 +28,7 @@ pair_style  gran model hertz tangential history # Hertzian without cohesion
 pair_coeff  * *
 
 # timestep, gravity
-timestep    0.0001 #0.00001
+timestep    0.00001
 fix         gravi all gravity 0.0 vector 0.0 0.0 -1.0
 
 # walls
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
index f9570ce0..58d050e4 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
@@ -15,7 +15,7 @@ processors      2 2 1
 # read the restart file
 read_restart    ../DEM/post/restart/liggghts.restart
 
-neighbor        0.001 bin
+neighbor        0.0005 bin
 neigh_modify    delay 0 binsize 0.01 ##0.01
 
 
@@ -31,7 +31,7 @@ pair_style  gran model hertz tangential history # Hertzian without cohesion
 pair_coeff  * *
 
 # timestep, gravity
-timestep    0.0001 #0.00001
+timestep    0.00001
 fix         gravi all gravity 0.0 vector 0.0 -1.0 0.0
 
 # walls
@@ -47,7 +47,7 @@ fix     cfd all couple/cfd couple_every 1 mpi
 fix     cfd2 all couple/cfd/force
 
 # this should invoke chemistry
-fix     cfd3 all couple/cfd/chemistry n_species 2 species_names O2 CO2
+fix     cfd3 all couple/cfd/chemistry n_species 3 species_names O2 CO CO2
 
 # this should shrink the particle
 fix     cfd4 all chem/shrink speciesA O2 molMassA +31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 2.5e-7