updated V-blender tutorial

tutorials/sphereGranFlow/V-blender/caseSetup/interaction

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

tutorials/sphereGranFlow/V-blender/caseSetup/particleInsertion

@@ -2,71 +2,100 @@
 |  phasicFlow File                                                            | 
 |  copyright: www.cemf.ir                                                     | 
 \* ------------------------------------------------------------------------- */ 
-
 objectName particleInsertion;
 objectType dicrionary;
+/* ------------------------------------------------------------------------- */
 
+// is insertion active?
+active  yes;
 
-active  yes;          // is insertion active?
+// not implemented for yes
-
+collisionCheck No;
-collisionCheck No;   // not implemented for yes
 
 /*
-one layer of particles are packed one-by-one using tho insertion steps. 
+two layers of particles are packed one-by-one using tho insertion steps. 
 */
 
-layer0
+rightregion
 {
-   type		  cylinderRegion; // type of insertion region
+   // type of insertion region
-   rate 	  5000;  // insertion rate (particles/s)
+   type		  cylinderRegion;
-   startTime 	  0.0;   // (s)
+   
-   endTime   	  1;     // (s)
+   // insertion rate (particles/s)
-   interval       0.025; //s  
+   rate 	  10000;
+   
+   // start time of small particles insertion (s)
+   startTime 	  0.0;
+   
+   // end time of small particles insertion (s)
+   endTime   	  1;
+   
+   // time interval of small particles insertion (s)
+   interval       0.025;  
 
    cylinderRegionInfo 
    {
-   	radius 	0.1;       // radius of cylinder (m)
+        // radius of cylinder (m)
-   	p1 ( 0.0950615  0.12 0.5011585 ); // (m,m,m)
+   	radius 	0.1;
-   	p2 ( 0.1150615  0.12 0.4811585); // (m,m,m)
+   	
+   	// coordinates of center of both ends of the insertion cylinder on 
+   	// the right side of the V-blender (m,m,m)
+   	p1 (0.0950615  0.12 0.5011585);
+   	p2 (0.1150615  0.12 0.4811585);
    }
    
    setFields
    {
-     velocity    realx3 (1.2 0.0 -1.2); // initial velocity of inserted particles 
+     // initial velocity of inserted particles
+     velocity    realx3 (1.2 0.0 -1.2); 
    }
    
    mixture
    {
-      lightSphere 1; // mixture composition of inserted particles 
+      // first layer of inserted particles
+      smallSphere 1;
       
    }
 }
 
-layer1
+leftregion
 {
+   // type of insertion region
    type         cylinderRegion;
-   rate         5000; // (particles/s)
+   
-   startTime    0.0; // (s)
+   // insertion rate (particles/s)
-   endTime      1; // (s)
+   rate         10000;
-   interval     0.025; //s  
+   
+   // start time of large particles insertion (s)
+   startTime    1.5;
+   
+   // end time of large particles insertion (s)
+   endTime      2.5;
+   
+   // time interval of large particles insertion (s)
+   interval     0.025;
 
    cylinderRegionInfo 
    {
+      // radius of cylinder (m)
       radius   0.1;
-      p1 ( 0.7562545  0.12 0.50079); // (m,m,m)
+      
-      p2 ( 0.7362545  0.12 0.48079); // (m,m,m)
+      // coordinates of center of both ends of the insertion cylinder on
+      // the left side of the V-blender (m,m,m)      
+      p1 ( 0.7562545  0.12 0.50079);
+      p2 ( 0.7362545  0.12 0.48079);
    }
    
    setFields
    {
+     // initial velocity of inserted particles
      velocity    realx3 (-1.2 0.0 -1.2);
    }
    
    mixture
    {
-      heavySphere 1; // only heavySphere
+      // second layer of inserted particles
+      largeSphere 1;
       
    }
 }
-
-

tutorials/sphereGranFlow/V-blender/caseSetup/sphereShape

@@ -2,10 +2,15 @@
 |  phasicFlow File                                                            | 
 |  copyright: www.cemf.ir                                                     | 
 \* ------------------------------------------------------------------------- */ 
-
 objectName 	sphereDict;
 objectType 	sphereShape;
+/* ------------------------------------------------------------------------- */
 
-names 		(lightSphere heavySphere); 	// names of shapes 
+// names of particles
-diameters 	(0.01 0.0101);				// diameter of shapes 
+names 		(smallSphere largeSphere);
-materials	(lightMat heavyMat);		// material names for shapes 
+
+// diameter of particles
+diameters 	(0.01 0.0101);
+
+// material names for particles
+materials	(lightMat lightMat); 

tutorials/sphereGranFlow/V-blender/settings/geometryDict

@@ -2,9 +2,9 @@
 |  phasicFlow File                                                            | 
 |  copyright: www.cemf.ir                                                     | 
 \* ------------------------------------------------------------------------- */  
-
 objectName 	geometryDict;
 objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
 
 // motion model: rotating object around an axis
 motionModel rotatingAxisMotion; 
@@ -12,12 +12,19 @@ motionModel rotatingAxisMotion;
 surfaces
 {
 	
-	rod
+	body
 	{
-		type 	 stlWall;   // type of the wall
+	        // type of the wall
-		file 	 cone.stl;	// file name in stl folder
+		type 	 stlWall;
-		material prop1;   // material name of this wall
+		
-		motion 	 rotAxis;	// motion component name 
+		// file name in stl folder
+		file 	 body.stl;
+		
+		// material name of this wall
+		material wallMat;
+		
+		// motion component name
+		motion 	 rotAxis;
 	}
 
 	
@@ -28,10 +35,19 @@ rotatingAxisMotionInfo
 {
 	rotAxis 
 	{
-		p1 (0.128228 0.116446 0.297901);	// first point for the axis of rotation 
+	        // first point for the axis of rotation
-		p2 (0.722596 0.116459 0.297901);	// second point for the axis of rotation
+		p1 (0.128228 0.116446 0.297901);
-		omega 3.14; 		// rotation speed (rad/s).
+		
-		startTime 1.5;
+		// second point for the axis of rotation
+		p2 (0.722596 0.116459 0.297901);
+		
+		// rotation speed (rad/s)
+		omega 3.14;
+		
+		// start time of rotation
+		startTime 3;
+		
+		// end time of rotation
 		endTime 10;
 	}
 }

tutorials/sphereGranFlow/V-blender/settings/particlesDict

@@ -2,9 +2,9 @@
 |  phasicFlow File                                                            | 
 |  copyright: www.cemf.ir                                                     | 
 \* ------------------------------------------------------------------------- */  
-
 objectName 	particlesDict;
 objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
 
 setFields
 {
@@ -16,10 +16,17 @@ setFields
 
 	defaultValue 
 	{
-		velocity 	realx3 	(0 0 0); // linear velocity (m/s)
+		// linear velocity (m/s)
-		acceleration 	realx3 	(0 0 0); // linear acceleration (m/s2)
+		velocity 	realx3 	(0 0 0);
-		rotVelocity 	realx3 	(0 0 0); // rotational velocity (rad/s)
+		
-		shapeName 	word	lightSphere; // name of the particle shape 
+		// linear acceleration (m/s2)
+		acceleration 	realx3 	(0 0 0);
+		
+		// rotational velocity (rad/s)
+		rotVelocity 	realx3 	(0 0 0);
+		
+		// name of the particle shape
+		shapeName 	word	smallSphere;
 	}
 
 	selectors
@@ -29,10 +36,12 @@ setFields
 // positions particles 
 positionParticles
 {
-	method empty;     // ordered positioning
+        // no particle is initially in the simulation
-
+	method empty;
-	maxNumberOfParticles 80000; // maximum number of particles in the simulation
-	mortonSorting Yes;          // perform initial sorting based on morton code?   
 
+        // maximum number of particles in the simulation
+	maxNumberOfParticles 25000;
 	
+	// perform initial sorting based on morton code?
+	mortonSorting Yes;
 }

tutorials/sphereGranFlow/V-blender/settings/settingsDict

@@ -3,34 +3,45 @@
 |  copyright: www.cemf.ir                                                     | 
 \* ------------------------------------------------------------------------- */  
 objectName 	settingsDict;
-objectType 	dictionary;;
+objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
 
-run 			rotatingDrumSmall;
+run 		rotatingVblender;
 
-dt 			0.00001; 	// time step for integration (s)
+// time step for integration (s)
+dt 			0.00001;
 
-startTime 		0; 		// start time for simulation 
+// start time for simulation
+startTime 		0;
 
-endTime 		10; 		// end time for simulation 
+// end time for simulation
+endTime 		10;
 
-saveInterval 	0.05; 		        // time interval for saving the simulation
+// time interval for saving the simulation
+saveInterval    	0.05;
 
-timePrecision   6;			// maximum number of digits for time folder 
+// maximum number of digits for time folder
+timePrecision    	6; 
 
-g 				(0 0 -9.8); // gravity vector (m/s2) 
+// gravity vector (m/s2)
+g 			(0 0 -9.8);
 
 /*
 	Simulation domain 
 	every particles that goes outside this domain is deleted.
 */
+
 domain 
 {
 	min (-0.1 -0.4 0);
 	max (0.86 0.6 0.6);
 }
 
-integrationMethod 		AdamsBashforth2; 	// integration method 
+// integration method
+integrationMethod 	AdamsBashforth2;
 
-timersReport 			Yes;  	// report timers?
+// report timers?
+timersReport 		Yes;
 
-timersReportInterval   	0.01;	// time interval for reporting timers
+// time interval for reporting timers
+timersReportInterval   	0.01;