rotatingDrumSmall updated.

tutorials/sphereGranFlow/rotatingDrumSmall/caseSetup/interaction

@@ -6,17 +6,33 @@ objectName interaction;
 objectType dicrionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-
 materials                        (prop1);          // a list of materials names
+
 densities                       (1000.0);          // density of materials [kg/m3]
  
 contactListType        sortedContactList;
 
+contactSearch
+{
+   method                            NBS;          // method for broad search particle-particle
+   
+   updateInterval                     10;
+
+   sizeRatio                         1.1;
+
+   cellExtent                       0.55;
+
+   adjustableBox                     Yes;
+} 
+
 model
 {
    contactForceModel nonLinearNonLimited;
+
    rollingFrictionModel           normal;
 
+   // Property (solid-solid Properties)
+
    Yeff                          (1.0e6);          // Young modulus [Pa]
 
    Geff                          (0.8e6);          // Shear modulus [Pa]
@@ -30,24 +46,6 @@ model
    mu                              (0.3);          // dynamic friction 
    
    mur                             (0.1);          // rolling friction     
-        
 }
 
-contactSearch
-{
-   method         NBS;           // method for broad search particle-particle
-   wallMapping    cellMapping;   // method for broad search particle-wall 
 
-   NBSInfo
-   {
-      updateFrequency 20;        // each 20 timesteps, update neighbor list 
-      sizeRatio      1.1;        // bounding box size to particle diameter (max)
-   }
-
-   cellMappingInfo
-   {
-      updateFrequency 20;        // each 20 timesteps, update neighbor list  
-      cellExtent     0.7;        // bounding box for particle-wall search (> 0.5)
-   }
-
-}

tutorials/sphereGranFlow/rotatingDrumSmall/caseSetup/particleInsertion

@@ -1,14 +0,0 @@
-/* -------------------------------*- C++ -*--------------------------------- *\ 
-|  phasicFlow File                                                            | 
-|  copyright: www.cemf.ir                                                     | 
-\* ------------------------------------------------------------------------- */ 
-objectName particleInsertion;
-objectType dicrionary;
-fileFormat  ASCII;
-/*---------------------------------------------------------------------------*/
-
-active  no;          // is insertion active?
-
-collisionCheck No;   // not implemented for yes
-
-

tutorials/sphereGranFlow/rotatingDrumSmall/caseSetup/shapes

@@ -6,7 +6,10 @@ objectName 	sphereDict;
 objectType 	sphereShape;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-
 names 		  (sphere1); 	      // names of shapes 
+
 diameters 	    (0.004);	      // diameter of shapes 
+
 materials	    (prop1);	      // material names for shapes 
+
+

tutorials/sphereGranFlow/rotatingDrumSmall/settings/domainDict

@@ -0,0 +1,58 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	domainDict;
+objectType 	dictionary;
+fileFormat  ASCII;
+/*---------------------------------------------------------------------------*/
+globalBox          // Simulation domain: every particles that goes outside this domain will be deleted
+{
+	min    (-0.12 -0.12 0);
+	
+	max (0.12   0.12 0.11);
+}
+
+boundaries
+{   
+    // Determines how often (how many iterations) do you want to 
+	// rebuild the list of particles in the neighbor list 
+	// of all boundaries in the simulation domain 
+	neighborListUpdateInterval 50;
+	
+	// Determines how often do you want to update the new changes in the boundary
+	updateInterval         0.0001;
+
+	// The distance from the boundary plane within which particles are marked to be in the boundary list 
+	neighborLength          0.004;
+
+	left
+	{
+		type  exit;          // other options: periodict, reflective 	
+	}
+
+	right
+	{
+		type  exit;          // other options: periodict, reflective 
+	}
+
+	bottom
+	{
+		type  exit;          // other options: periodict, reflective 
+	}
+
+	top
+	{
+		type  exit;          // other options: periodict, reflective 
+	}
+
+	rear
+	{
+		type  exit;          // other options: periodict, reflective 
+	}
+
+	front
+	{
+		type  exit;          // other options: periodict, reflective 
+	}
+}

tutorials/sphereGranFlow/rotatingDrumSmall/settings/geometryDict

@@ -6,66 +6,86 @@ objectName 	geometryDict;
 objectType 	dictionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
+motionModel     rotatingAxis; 		   // motion model can be rotatingAxis or stationary or vibrating
 
-// motion model: rotating object around an axis
+rotatingAxisInfo					   // information for rotatingAxisMotion motion model
-motionModel rotatingAxisMotion; 
+{
+	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)
+	}
+}
 
 surfaces
 {
 	/*
-	    A cylinder with begin and end radii 0.12 m and axis points at (0 0 0) 
+	    A cylinder with begin and end radii 0.12 m and axis points at (0 0 0) and (0 0 0.1)
-	and (0 0 0.1)
 	*/
+
 	cylinder
 	{
-		type cylinderWall;  // type of the wall
+		type    cylinderWall;  		   // other options: cuboidWall and planeWall
+
 		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
+		type       planeWall;		   // other options: cuboidWall and cylinderWall 
+
 		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;
+		type       planeWall;		   // other options: cuboidWall and cylinderWall 
-		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;
-	}
 
+		p1 (-0.12 -0.12 0.1);          // first point of the wall
 
-}
+		p2 ( 0.12 -0.12 0.1);          // second point
 
-// information for rotatingAxisMotion motion model 
+		p3 ( 0.12  0.12 0.1);          // third point
-rotatingAxisMotionInfo
+
-{
+		p4 (-0.12  0.12 0.1);          // fourth point
-	rotAxis 
+
-	{
+		material       prop1;          // material name of the wall
-		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
+		motion       rotAxis;          // motion component name 
-		omega 1.214; 		// rotation speed (rad/s)
 	}
 }
+
+

tutorials/sphereGranFlow/rotatingDrumSmall/settings/particlesDict

@@ -6,7 +6,6 @@ objectName 	particlesDict;
 objectType 	dictionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-
 setFields
 {
 	/*
@@ -18,8 +17,11 @@ setFields
 	defaultValue 
 	{
 		velocity 	   realx3  (0 0 0); 		 // linear velocity (m/s)
+
 		acceleration   realx3  (0 0 0); 		 // linear acceleration (m/s2)
+
 		rVelocity 	   realx3  (0 0 0); 		 // rotational velocity (rad/s)
+
 		shapeName      word	   sphere1; 		 // name of the particle shape 
 	}
 
@@ -27,24 +29,26 @@ setFields
 	{}
 }
 
-// positions particles 
+positionParticles					             // positions particles
-positionParticles
 {
-	method positionOrdered;     // ordered positioning
+	method                      ordered;         // can be ordered or random or empty
 
-	maxNumberOfParticles 40000; // maximum number of particles in the simulation
+	orderedInfo
-	mortonSorting Yes;          // perform initial sorting based on morton code?   
-
-	box  // box for positioning particles 
 	{
-		min (-0.08 -0.08 0.015);    // lower corner point of the box 
+		diameter                 0.004; 	     // diameter of particles
-		max ( 0.08  0.08 0.098);	// upper corner point of the box 
-	}
 
-	positionOrderedInfo
-	{
-		diameter 0.004; 	// minimum space between centers of particles
 		numPoints                30000; 		 // number of particles in the simulation 
+
 		axisOrder              (z y x);  		 // axis order for filling the space with particles
 	}
+
+	regionType                     box;          // can be box or cylinder or sphere   
+
+	boxInfo                                      // box for positioning particles 
+	{
+		min        (-0.08 -0.08 0.015);          // lower corner point of the box 
+		
+		max        ( 0.08  0.08 0.098);	         // upper corner point of the box 
+	}	
 }
+

tutorials/sphereGranFlow/rotatingDrumSmall/settings/settingsDict

@@ -6,35 +6,30 @@ objectName 	settingsDict;
 objectType 	dictionary;
 fileFormat  ASCII;
 /*---------------------------------------------------------------------------*/
-
 run 			rotatingDrumSmall;
 
 dt 				          0.00001; 	        // time step for integration (s)
 
-startTime 		0; 			// start time for simulation 
+startTime 		              0.0; 		    // start time for simulation 
 
-endTime 		10; 		// end time for simulation 
+endTime 		             10.0; 		    // end time for simulation 
 
-saveInterval 	0.1; 		// time interval for saving the simulation
+saveInterval 	             0.01; 		    // time interval for saving the simulation
 
 timePrecision                   6;			// maximum number of digits for time folder 
 
 g 				       (0 -9.8 0);          // gravity vector (m/s2) 
 
-/*
+// save necessary (i.e., required) data on disk
-	Simulation domain 
+includeObjects (diameter);                                                 
-	every particles that goes outside this domain is deleted.
+// exclude unnecessary data from saving on disk
-*/
+excludeObjects (rVelocity.dy1 pStructPosition.dy1 pStructVelocity.dy1);    
-domain 
-{
-	min (-0.12 -0.12 0);
-	max (0.12   0.12 0.11);
-}
 
 integrationMethod       AdamsBashforth2; 	      // integration method 
 
-writeFormat				ascii;
+writeFormat				          ascii;          // data writting format (ascii or binary)
 
-timersReport 			Yes;  	// report timers?
+timersReport 			            Yes;  	      // report timers
 
 timersReportInterval   	           0.01;	      // time interval for reporting timers
+