Added Rotating Drum with Baffles Tutorial

tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.txt

@@ -0,0 +1,10 @@
+To run this case please follow the below Steps:
+1- Open a terminal in the RotatingDrumWithBaffles
+2- Enter the code ./runThisCase to run the Case
+3- After finishing the run time please type pFlowtoVTK to convert data to VTK format
+4- After Post-Producting and viewing the case please enter ./cleanThisCase to delete data
+
+
+For more questions please contact:
+omidava.khosravi@gmail.com
+<3

tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/interaction

@@ -0,0 +1,72 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+
+objectName interaction;
+objectType dicrionary;
+
+materials   (lightMat heavyMat wallMat); // a list of materials names
+densities   (1000 1500.0 2500);      // density of materials [kg/m3]
+
+contactListType   sortedContactList;
+
+model
+{
+   contactForceModel nonLinearLimited;
+   rollingFrictionModel normal;
+
+   /*
+   Property (lightMat-lightMat   lightMat-heavyMat    lightMat-wallMat
+                                 heavyMat-heavyMat    heavyMat-wallMat
+                                                      wallMat-wallMat );
+   */
+
+   Yeff (1.0e6 1.0e6 1.0e6  // Young modulus [Pa]
+               1.0e6 1.0e6
+                     1.0e6);
+   
+   Geff (0.8e6 0.8e6 0.8e6  // Shear modulus [Pa]
+               0.8e6 0.8e6
+                     0.8e6);
+
+   nu    (0.25 0.25  0.25  // Poisson's ratio [-]
+               0.25  0.25
+                     0.25);
+
+   en (0.97  0.97    0.85   // coefficient of normal restitution
+             0.97    0.85
+                     1.00);
+
+   et (1.0   1.0  1.0   // coefficient of tangential restitution 
+             1.0  1.0
+                  1.0);
+
+   mu (0.65   0.65 0.35   // dynamic friction 
+              0.65 0.35
+                   0.35);
+
+   mur (0.1  0.1 0.1    // rolling friction 
+             0.1 0.1
+                 0.1);
+        
+}
+
+contactSearch
+{
+   method         NBS;           // method for broad search particle-particle
+   wallMapping    cellMapping;   // method for broad search particle-wall 
+
+   NBSInfo
+   {
+      updateFrequency 10;        // each 20 timesteps, update neighbor list 
+      sizeRatio      1.1;        // bounding box size to particle diameter (max)
+   }
+
+   cellMappingInfo
+   {
+      updateFrequency 10;        // each 20 timesteps, update neighbor list  
+      cellExtent     0.6;        // bounding box for particle-wall search (> 0.5)
+   }
+
+}

tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion

@@ -0,0 +1,71 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+
+objectName particleInsertion;
+objectType dicrionary;
+
+
+active  yes;          // is insertion active?
+
+collisionCheck No;   // not implemented for yes
+
+/*
+five layers of particles are packed one-by-one using 5 insertion steps. 
+*/
+
+layer0
+{
+   type	          cylinderRegion; // type of insertion region
+   rate 	  5000;  // insertion rate (particles/s)
+   startTime 	  0;      // (s)
+   endTime   	  1;    // (s)
+   interval       0.025;  //s  
+
+   cylinderRegionInfo 
+   {
+   	p2 (-0.15 0.25 0.05); // (m,m,m)
+   	p1 (-0.15 0.24	0.05); // (m,m,m)
+   	radius 0.035;  // radius of cylinder (m)
+   }
+   
+   setFields
+   {
+     velocity    realx3 (0.0 -0.6 0.0); // initial velocity of inserted particles 
+   }
+   
+   mixture
+   {
+      lightSphere 1; // mixture composition of inserted particles 
+      
+   }
+}
+
+layer1
+{
+   type         cylinderRegion;
+   rate         7500; // (particles/s)
+   startTime    0; // (s)
+   endTime      1; // (s)
+   interval     0.025; //s  
+
+   cylinderRegionInfo 
+   {
+      p2 (-0.23	 0.25	0.05); // (m,m,m)
+      p1 (-0.23	 0.24	0.05); // (m,m,m)
+      radius 0.035;  // radius of cylinder (m)
+   }
+   
+   setFields
+   {
+     velocity    realx3 (0.0 -0.6 0.0);
+   }
+   
+   mixture
+   {
+      heavySphere 1; // only heavySphere
+      
+   }
+}
+

tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/sphereShape

@@ -0,0 +1,11 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+
+objectName 	sphereDict;
+objectType 	sphereShape;
+
+names 		(lightSphere heavySphere); 	// names of shapes 
+diameters 	(0.003 0.005);				// diameter of shapes 
+materials	(lightMat heavyMat);		// material names for shapes 

tutorials/sphereGranFlow/RotatingDrumWithBaffles/cleanThisCase

@@ -0,0 +1,7 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+
+ls | grep -P "^(([0-9]+\.?[0-9]*)|(\.[0-9]+))$" | xargs -d"\n" rm -rf
+rm -rf VTK 
+
+#------------------------------------------------------------------------------

tutorials/sphereGranFlow/RotatingDrumWithBaffles/runThisCase

@@ -0,0 +1,21 @@
+#!/bin/sh
+cd ${0%/*} || exit 1    # Run from this directory
+echo "\n<--------------------------------------------------------------------->"
+echo "1) Creating particles"
+echo "<--------------------------------------------------------------------->\n"
+particlesPhasicFlow
+
+echo "\n<--------------------------------------------------------------------->"
+echo "2) Creating geometry"
+echo "<--------------------------------------------------------------------->\n"
+geometryPhasicFlow
+
+echo "\n<--------------------------------------------------------------------->"
+echo "3) Running the case"
+echo "<--------------------------------------------------------------------->\n"
+sphereGranFlow
+
+
+
+
+#------------------------------------------------------------------------------

tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict

@@ -0,0 +1,63 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+
+objectName 	geometryDict;
+objectType 	dictionary;
+
+// motion model: all surfaces are fixed
+motionModel rotatingAxisMotion; 
+
+surfaces
+{
+	Shell
+	{
+		type 	 stlWall;  	// type of the wall
+		file 	 shell.stl;	// file name in stl folder		
+		material wallMat;   // 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.4 0.10 0.0);	// first point of the wall
+		p2 ( 0.0 0.10 0.0);   // second point
+		p3 ( 0.0 0.35 0.0);   // third point
+		p4 (-0.4 0.35 0.0);   // fourth point 
+		material wallMat;         // material name of the wall  
+		motion rotAxis;			// motion component name 
+	}
+
+	/*
+	This is a plane wall at the front end of cylinder
+	*/
+	wall2
+	{
+		type planeWall;			// type of the wall
+		p1 (-0.4 0.10 0.1);	// first point of the wall
+		p2 ( 0.0 0.10 0.1);   // second point
+		p3 ( 0.0 0.35 0.1);   // third point
+		p4 (-0.4 0.35 0.1);   // fourth point 
+		material wallMat;         // material name of the wall  
+		motion rotAxis;			// motion component name 
+	}
+}
+	// information for rotatingAxisMotion motion model 
+rotatingAxisMotionInfo
+{
+	rotAxis 
+	{
+		p1 (-0.1974  0.2269  0);	// first point for the axis of rotation 
+		p2 (-0.1974  0.2269  0.1);	// second point for the axis of rotation
+		omega 2.38733; 		// rotation speed (rad/s) => 15 rpm
+		startTime 2;
+		endTime 9.5;
+	}
+}
+
+

tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict

@@ -0,0 +1,31 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+
+objectName 	particlesDict;
+objectType 	dictionary;
+
+setFields
+{
+	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	lightSphere; 	// name of the particle shape 
+	}
+
+	selectors
+	{}
+}
+
+// positions particles 
+positionParticles
+{
+	method empty;     		// creates the required fields with zero particles (empty). 
+
+	maxNumberOfParticles 	50000; // maximum number of particles in the simulation
+	mortonSorting 			Yes;    // perform initial sorting based on morton code?   
+
+}

tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/settingsDict

@@ -0,0 +1,36 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	settingsDict;
+objectType 	dictionary;;
+
+run 			layerdSiloFilling;
+
+dt 				0.00001; 	// time step for integration (s)
+
+startTime 		0; 			// start time for simulation 
+
+endTime 		10;	 		// end time for simulation 
+
+saveInterval 	0.05; 		// 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.328 0.095 -0.025);
+	max (-0.068  0.355 0.125);
+}
+
+integrationMethod 		AdamsBashforth3; 	// integration method 
+
+timersReport 			Yes;  	// report timers?
+
+timersReportInterval   	0.01;	// time interval for reporting timers

tutorials/sphereGranFlow/RotatingDrumWithBaffles/stl/note

@@ -0,0 +1 @@
+/* The STL files should be Placed in this Folder *\