diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/p b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/p
index 634fa811..4e5c08c5 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/p
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/p
@@ -40,7 +40,8 @@ boundaryField
 
     outlet
     {
-        type            zeroGradient;
+        type            fixedValue;
+        value           $internalField;
     }
 }
 
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init
index 98f6916a..f8f0b06c 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init
@@ -28,7 +28,7 @@ pair_style  gran model hertz tangential history # Hertzian without cohesion
 pair_coeff  * *
 
 # timestep, gravity
-timestep    0.000001
+timestep    1e-6
 fix         gravi all gravity 0.0 vector 0.0 -1.0 0.0
 
 # walls
@@ -43,16 +43,15 @@ fix     zwalls2 all wall/gran model hertz tangential history primitive type 1 zp
 #set                 atom * diameter 0.07 density 2267 vx 0 vy 0 vz 0
 
 # create single particle in a specific spot
-#create_atoms 1 single 1.25 .40 0.31  units box
-#set atom 2 diameter 0.07 density 2267 vx 0 vy 0 vz 0
-#group part2 id 2  
+#create_atoms 1 single 0.25 .25 0.25  units box
+#set atom 1 diameter 0.07 density 2267 vx 0 vy 0 vz 0
 
 # particle distributions and insertion
-region  bc block 0. 1.5 0. 0.5 0. 0.5 units box
+region  bc block 0.25 1.25 0. 0.5 0. 0.5 units box
 fix     pts1 all particletemplate/sphere 1 atom_type 1 density constant 2267 radius constant 0.035 
 fix     pdd1 all particledistribution/discrete 1 1 pts1 1.0
 
-fix     ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. 0. insert_every once overlapcheck yes all_in yes particles_in_region 100 region bc
+fix     ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. 0. insert_every once overlapcheck yes all_in yes particles_in_region 10 region bc
 
 # apply nve integration to all particles that are inserted as single particles
 #fix     integr all nve/sphere
diff --git a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
index f623aa40..37efd16b 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
+++ b/tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
@@ -31,7 +31,7 @@ pair_style  gran model hertz tangential history # Hertzian without cohesion
 pair_coeff  * *
 
 # timestep, gravity
-timestep    0.000001
+timestep    1e-6
 fix         gravi all gravity 0.0 vector 0.0 -1.0 0.0
 
 # walls
@@ -51,7 +51,7 @@ fix     cfd3 all couple/cfd/chemistry n_species 2 species_names O2 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.5e3 rdef hertzpct 0.2
-fix     cfd4 all chem/shrink speciesA O2 molMassA 31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 25 rmin 0.005
+fix     cfd4 all chem/shrink speciesA O2 molMassA 31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 2.5e3 rmin 0.01
 
 # apply nve integration to all particles that are inserted as single particles
 fix     integr all nve/sphere