benchmarks update

benchmarks/rotatingDrum_4MParticles/caseSetup/interaction

@@ -0,0 +1,59 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+
+objectName interaction;
+objectType dicrionary;
+
+materials   (glassMat wallMat); // a list of materials names
+densities   (2500.0 2500);      // density of materials [kg/m3]
+
+contactListType   sortedContactList;
+
+model
+{
+   contactForceModel nonLinearLimited;
+   rollingFrictionModel normal;
+
+   Yeff (1.0e6 1.0e6  // Young modulus [Pa]
+               1.0e6);
+   
+   Geff (0.8e6 0.8e6  // Shear modulus [Pa]
+               0.8e6);
+
+   nu    (0.25 0.25  // Poisson's ratio [-]
+               0.25);
+
+   en (0.97  0.85   // coefficient of normal restitution
+             1.00);
+
+   et (1.0   1.0    // coefficient of tangential restitution 
+             1.0);
+
+   mu (0.65  0.65   // dynamic friction 
+             0.65);
+
+   mur (0.1  0.1    // rolling friction 
+             0.1);
+        
+}
+
+contactSearch
+{
+   method         NBS;           
+   wallMapping    cellsSimple;   
+
+   NBSInfo
+   {
+      updateFrequency 10;        // each 20 timesteps, update neighbor list 
+      sizeRatio      1.05;       // 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)
+   }
+
+}

benchmarks/rotatingDrum_4MParticles/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
+
+

benchmarks/rotatingDrum_4MParticles/caseSetup/sphereShape

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

benchmarks/rotatingDrum_4MParticles/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 
+
+#------------------------------------------------------------------------------

benchmarks/rotatingDrum_4MParticles/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
+
+
+
+
+#------------------------------------------------------------------------------

benchmarks/rotatingDrum_4MParticles/settings/geometryDict

@@ -0,0 +1,63 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+
+objectName 	geometryDict;
+objectType 	dictionary;
+
+motionModel rotatingAxisMotion; 
+
+surfaces
+{
+	
+	cylinder
+	{
+		type cylinderWall;  
+		p1 (0.0 0.0 0.0);   
+		p2 (0.0 0.0 1.6);  
+		radius1 0.2;		
+		radius2 0.2;		
+		resolution 24;      
+		material wallMat;   
+		motion rotAxis;
+	}
+
+	
+	wall1
+	{
+		type planeWall;			
+		p1 (-0.2 -0.2 0.0);	
+		p2 ( 0.2 -0.2 0.0);   
+		p3 ( 0.2  0.2 0.0);   
+		p4 (-0.2  0.2 0.0);   
+		material wallMat;     
+		motion rotAxis;		
+	}
+
+	/*
+	This is a plane wall at the front end of cylinder
+	*/
+	wall2
+	{
+		type planeWall;
+		p1 (-0.2 -0.2 1.6);
+		p2 ( 0.2 -0.2 1.6);
+		p3 ( 0.2  0.2 1.6);
+		p4 (-0.2  0.2 1.6);
+		material wallMat;
+		motion rotAxis;
+	}
+		
+}
+
+// information for rotatingAxisMotion motion model 
+rotatingAxisMotionInfo
+{
+	rotAxis 
+	{
+		p1 (0.0 0.0 0.0);	
+		p2 (0.0 0.0 1.0);	
+		omega 1.256; 		// rotation speed (rad/s) => 12 rpm
+	}
+}

benchmarks/rotatingDrum_4MParticles/settings/particlesDict

@@ -0,0 +1,44 @@
+/* -------------------------------*- 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	glassBead; // name of the particle shape 
+	}
+
+	selectors
+	{}
+}
+
+positionParticles
+{
+	method positionOrdered;     		
+
+	maxNumberOfParticles 	4000001; 
+	mortonSorting 			Yes;    
+
+	cylinder  // box for positioning particles 
+	{
+		p1 ( 0.0  0.0 0.01);    // lower corner point of the box 
+		p2 ( 0.0  0.0 1.59);	// upper corner point of the box 
+		radius 0.195;
+	}
+	
+	positionOrderedInfo
+	{
+		diameter 0.003; 	// minimum space between centers of particles
+		numPoints 4000000; 	// number of particles in the simulation 
+		axisOrder (z x y);  // axis order for filling the space with particles
+	}
+}

benchmarks/rotatingDrum_4MParticles/settings/settingsDict

@@ -0,0 +1,32 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	settingsDict;
+objectType 	dictionary;;
+
+run 			rotatingDrum_1;
+
+dt 				0.00001; 	// time step for integration (s)
+
+startTime 		0; 			// start time for simulation 
+
+endTime 		10; 		// end time for simulation 
+
+saveInterval 	0.2; 		// time interval for saving the simulation
+
+timePrecision   5;			// maximum number of digits for time folder 
+
+g 				(0 -9.8 0); // gravity vector (m/s2) 
+
+domain 
+{
+	min (-0.2 -0.2 -0.0);
+	max ( 0.2  0.2  1.6);
+}
+
+integrationMethod 		AdamsBashforth3; 	// integration method 
+
+timersReport 			Yes;  	
+
+timersReportInterval   	0.01;