tutorial: layeredSiloFilling

tutorials/sphereGranFlow/layeredSiloFilling/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    cellsSimple;   // 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)
+   }
+
+   cellsSimpleInfo
+   {
+      updateFrequency 10;        // each 20 timesteps, update neighbor list  
+      cellExtent     0.6;        // bounding box for particle-wall search (> 0.5)
+   }
+
+}

tutorials/sphereGranFlow/layeredSiloFilling/caseSetup/particleInsertion

@@ -0,0 +1,152 @@
+/* -------------------------------*- 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 	      15000;  // insertion rate (particles/s)
+   startTime 	  0;      // (s)
+   endTime   	  0.5;    // (s)
+   interval       0.025;  //s  
+
+   cylinderRegionInfo 
+   {
+   	radius 	0.09;       // radius of cylinder (m)
+   	p1 ( 0.0  0.0 0.1 ); // (m,m,m)
+   	p2 ( 0.0  0.0 0.11); // (m,m,m)
+   }
+   
+   setFields
+   {
+     velocity    realx3 (0.0 0.0 -0.6); // initial velocity of inserted particles 
+   }
+   
+   mixture
+   {
+      lightSphere 1; // mixture composition of inserted particles 
+      
+   }
+}
+
+layer1
+{
+   type         cylinderRegion;
+   rate         15000; // (particles/s)
+   startTime    0.7; // (s)
+   endTime      1.2; // (s)
+   interval     0.025; //s  
+
+   cylinderRegionInfo 
+   {
+      radius   0.09;
+      p1 ( 0.0  0.0 0.16 ); // (m,m,m)
+      p2 ( 0.0  0.0 0.17); // (m,m,m)
+   }
+   
+   setFields
+   {
+     velocity    realx3 (0.0 0.0 -0.6);
+   }
+   
+   mixture
+   {
+      heavySphere 1; // only heavySphere
+      
+   }
+}
+
+layer2
+{
+   type         cylinderRegion;
+   rate         15000; // (particles/s)
+   startTime    1.4; // (s)
+   endTime      1.9; // (s)
+   interval     0.025; //s  
+
+   cylinderRegionInfo 
+   {
+      radius   0.09;
+      p1 ( 0.0  0.0 0.2 ); // (m,m,m)
+      p2 ( 0.0  0.0 0.21); // (m,m,m)
+   }
+   
+   setFields
+   {
+     velocity    realx3 (0.0 0.0 -0.6);
+   }
+   
+   mixture
+   {
+      lightSphere 1; // only lightSphere
+      
+   }
+}
+
+layer3
+{
+   type         cylinderRegion;
+   rate         15000; // (particles/s)
+   startTime    2.1; // (s)
+   endTime      2.6; // (s)
+   interval     0.025; //s  
+
+   cylinderRegionInfo 
+   {
+      radius   0.09;
+      p1 ( 0.0  0.0 0.28 ); // (m,m,m)
+      p2 ( 0.0  0.0 0.29); // (m,m,m)
+   }
+   
+   setFields
+   {
+     velocity   realx3 (0.0 0.0 -0.6);
+   }
+   
+   mixture
+   {
+      heavySphere 1;
+      
+   }
+}
+
+layer4
+{
+   type         cylinderRegion;
+   rate         15000; // (particles/s)
+   startTime    2.8; // (s)
+   endTime      3.3; // (s)
+   interval     0.025; //s  
+
+   cylinderRegionInfo 
+   {
+      radius   0.09;
+      p1 ( 0.0  0.0 0.37 ); // (m,m,m)
+      p2 ( 0.0  0.0 0.38); // (m,m,m)
+   }
+   
+   setFields
+   {
+     velocity    realx3 (0.0 0.0 -0.6);
+   }
+   
+   mixture
+   {
+      lightSphere 1;
+      
+   }
+}
+

tutorials/sphereGranFlow/layeredSiloFilling/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.007 0.007);				// diameter of shapes 
+materials	(lightMat heavyMat);		// material names for shapes