diff --git a/tutorials/postprocessPhasicFlow/segregation/settings/settingsDict b/tutorials/postprocessPhasicFlow/segregation/settings/settingsDict
index e6e2d9b9..7ee75f0d 100644
--- a/tutorials/postprocessPhasicFlow/segregation/settings/settingsDict
+++ b/tutorials/postprocessPhasicFlow/segregation/settings/settingsDict
@@ -5,7 +5,7 @@
 objectName 	settingsDict;
 objectType 	dictionary;;
 
-run 			rotatingDrumSmall;
+run 			postprocessSegregation;
 
 dt 				0.00001; 	// time step for integration (s)
 
@@ -15,7 +15,7 @@ endTime 		10; 		// end time for simulation
 
 saveInterval 	0.1; 		// time interval for saving the simulation
 
-timePrecision   6;			// maximum number of digits for time folder 
+timePrecision   4;			// maximum number of digits for time folder 
 
 g 				(0 -9.8 0); // gravity vector (m/s2) 
 
diff --git a/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/interaction b/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/interaction
new file mode 100755
index 00000000..50ab5159
--- /dev/null
+++ b/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/interaction
@@ -0,0 +1,52 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+
+objectName interaction;
+objectType dicrionary;
+
+materials      (prop1);    // a list of materials names
+densities      (1000.0);   // density of materials [kg/m3]
+
+contactListType   sortedContactList; 
+
+model
+{
+   contactForceModel nonLinearNonLimited;
+   rollingFrictionModel normal;
+
+   Yeff  (1.0e6);       // Young modulus [Pa]
+
+   Geff  (0.8e6);       // Shear modulus [Pa]
+
+   nu    (0.25);        // Poisson's ratio [-]
+
+   en    (0.7);         // coefficient of normal restitution
+
+   et    (1.0);         // coefficient of tangential restitution 
+
+   mu    (0.3);         // dynamic friction 
+
+   mur   (0.1);         // rolling friction 
+        
+}
+
+contactSearch
+{
+   method         NBS;           // method for broad search particle-particle
+   wallMapping    cellsSimple;   // method for broad search particle-wall 
+
+   NBSInfo
+   {
+      updateFrequency 10;        // each 10 timesteps, update neighbor list 
+      sizeRatio      1.1;        // bounding box size to particle diameter (max)
+   }
+
+   cellsSimpleInfo
+   {
+      updateFrequency 10;        // each 10 timesteps, update neighbor list  
+      cellExtent     0.6;        // bounding box for particle-wall search (> 0.5)
+   }
+
+}
\ No newline at end of file
diff --git a/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/particleInsertion b/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/particleInsertion
new file mode 100755
index 00000000..eec7b7f9
--- /dev/null
+++ b/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/particleInsertion
@@ -0,0 +1,14 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+
+objectName particleInsertion;
+objectType dicrionary;
+
+
+active  no;          // is insertion active?
+
+collisionCheck No;   // not implemented for yes
+
+
diff --git a/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/sphereShape b/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/sphereShape
new file mode 100755
index 00000000..15060c08
--- /dev/null
+++ b/tutorials/sphereGranFlow/binarySystemOfParticles/caseSetup/sphereShape
@@ -0,0 +1,11 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+
+objectName 	sphereDict;
+objectType 	sphereShape;
+
+names 		(smallSphere largeSphere); 	// names of shapes 
+diameters 	(0.003 	0.005);	// diameter of shapes 
+materials	(prop1  prop1);	// material names for shapes 
diff --git a/tutorials/sphereGranFlow/binarySystemOfParticles/settings/geometryDict b/tutorials/sphereGranFlow/binarySystemOfParticles/settings/geometryDict
new file mode 100644
index 00000000..7a2667f7
--- /dev/null
+++ b/tutorials/sphereGranFlow/binarySystemOfParticles/settings/geometryDict
@@ -0,0 +1,70 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+
+objectName 	geometryDict;
+objectType 	dictionary;
+
+// motion model: rotating object around an axis
+motionModel rotatingAxisMotion; 
+
+surfaces
+{
+	/*
+	    A cylinder with begin and end radii 0.12 m and axis points at (0 0 0) 
+	and (0 0 0.1)
+	*/
+	cylinder
+	{
+		type cylinderWall;  // type of the wall
+		p1 (0.0 0.0 0.0);   // begin point of cylinder axis
+		p2 (0.0 0.0 0.1);   // end point of cylinder axis
+		radius1 0.12;		// radius at p1
+		radius2 0.12;		// radius at p2
+		resolution 24;      // number of divisions
+		material prop1;     // material name of this wall
+		motion rotAxis;		// motion component name 
+	}
+
+	/*
+	This is a plane wall at the rear end of cylinder
+	*/
+	wall1
+	{
+		type planeWall;			// type of the wall
+		p1 (-0.12 -0.12 0.0);	// first point of the wall
+		p2 ( 0.12 -0.12 0.0);   // second point
+		p3 ( 0.12  0.12 0.0);   // third point
+		p4 (-0.12  0.12 0.0);   // fourth point 
+		material prop1;         // material name of the wall  
+		motion rotAxis;			// motion component name 
+	}
+
+	/*
+	This is a plane wall at the front end of cylinder
+	*/
+	wall2
+	{
+		type planeWall;
+		p1 (-0.12 -0.12 0.1);
+		p2 ( 0.12 -0.12 0.1);
+		p3 ( 0.12  0.12 0.1);
+		p4 (-0.12  0.12 0.1);
+		material prop1;
+		motion rotAxis;
+	}
+		
+	
+}
+
+// information for rotatingAxisMotion motion model 
+rotatingAxisMotionInfo
+{
+	rotAxis 
+	{
+		p1 (0.0 0.0 0.0);	// first point for the axis of rotation 
+		p2 (0.0 0.0 1.0);	// second point for the axis of rotation
+		omega 1.214; 		// rotation speed (rad/s)
+	}
+}
\ No newline at end of file
diff --git a/tutorials/sphereGranFlow/binarySystemOfParticles/settings/particlesDict b/tutorials/sphereGranFlow/binarySystemOfParticles/settings/particlesDict
new file mode 100644
index 00000000..c3e05448
--- /dev/null
+++ b/tutorials/sphereGranFlow/binarySystemOfParticles/settings/particlesDict
@@ -0,0 +1,68 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+
+objectName 	particlesDict;
+objectType 	dictionary;
+
+// positions particles 
+positionParticles
+{
+	method positionOrdered;     // ordered positioning
+
+	maxNumberOfParticles 30001; // maximum number of particles in the simulation
+	mortonSorting Yes;          // perform initial sorting based on morton code?   
+
+	cylinder  // cylinder region for positioning particles 
+	{
+		p1 (0.0 0.0 0.003); 	// begin point of cylinder axis
+		p2 (0.0 0.0 0.097);		// end point of cylinder axis
+		radius 		0.117;		// radius of cylinder 
+	}
+		
+	positionOrderedInfo
+	{		
+		diameter 0.005; 	// minimum space between centers of particles
+		numPoints 30000; 	// number of particles in the simulation 
+		axisOrder (z x y);  // axis order for filling the space with particles
+	}
+}
+
+setFields
+{
+	/*
+	    Default value for fields defined for particles
+	These fields should always be defined for simulations with 
+	spherical particles.*/
+
+	defaultValue 
+	{
+		velocity 		realx3 	(0 0 0); // linear velocity (m/s)
+		acceleration 	realx3 	(0 0 0); // linear acceleration (m/s2)
+		rotVelocity 	realx3 	(0 0 0); // rotational velocity (rad/s)
+		shapeName 		word	smallSphere; // name of the particle shape 
+	}
+
+	selectors
+	{
+		shapeAssigne
+		{
+			selector 	rangeAll; // type of point selector
+			rangeAllInfo
+			{
+				begin 	0;			// begin index of points
+				end 	30000;		// end index of points 
+				stride  3;			// stride for selector 
+			}
+			fieldValue  // fields that the selector is applied to 
+			{
+				/*
+				sets shapeName of the selected points to largeSphere*/
+				shapeName 	word 	largeSphere;
+			}
+		}
+	}
+}
+
+
diff --git a/tutorials/sphereGranFlow/binarySystemOfParticles/settings/settingsDict b/tutorials/sphereGranFlow/binarySystemOfParticles/settings/settingsDict
new file mode 100644
index 00000000..5e633b90
--- /dev/null
+++ b/tutorials/sphereGranFlow/binarySystemOfParticles/settings/settingsDict
@@ -0,0 +1,36 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	settingsDict;
+objectType 	dictionary;;
+
+run 			binarySystemofParticles;
+
+dt 				0.00001; 	// time step for integration (s)
+
+startTime 		0; 			// start time for simulation 
+
+endTime 		10; 		// end time for simulation 
+
+saveInterval 	0.1; 		// time interval for saving the simulation
+
+timePrecision   6;			// maximum number of digits for time folder 
+
+g 				(0 -9.8 0); // gravity vector (m/s2) 
+
+/*
+	Simulation domain 
+	every particles that goes outside this domain is deleted.
+*/
+domain 
+{
+	min (-0.12 -0.12 0);
+	max (0.12   0.12 0.1);
+}
+
+integrationMethod 		AdamsBashforth2; 	// integration method 
+
+timersReport 			Yes;  	// report timers?
+
+timersReportInterval   	0.01;	// time interval for reporting timers