Merge pull request #83 from omid-khosravi/main

Rotary Air-Lock Valve

tutorials/sphereGranFlow/RotaryAirLockValve/caseSetup/interaction

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

tutorials/sphereGranFlow/RotaryAirLockValve/caseSetup/particleInsertion

@@ -0,0 +1,55 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */
+objectName particleInsertion;
+objectType dicrionary;
+fileFormat  ASCII;
+/*---------------------------------------------------------------------------*/
+
+// is insertion active?
+active  yes;
+
+// not implemented for yes
+collisionCheck No;
+
+/*
+one layer of particles is added using an insertion step. 
+*/
+//Insertion Layer
+layer0
+{
+
+	// type of insertion region
+	type  boxRegion;
+
+	// insertion rate (particles/s)
+	rate 	  4000;
+
+	// Start time of Particles insertion (s)
+	startTime 	  0;
+
+	// End time of Particles insertion (s)	
+	endTime   	  7;
+
+	// Time Interval of Particles insertion (s)
+	interval       0.025;
+
+	// Coordinates of BoxRegion (m,m,m)
+	boxRegionInfo 
+	{
+	min ( 0.48 0.58 0.01 ); // (m,m,m)
+	max ( 0.64 0.59 0.05 ); // (m,m,m)
+	}
+
+	setFields
+	{
+	// initial velocity of inserted particles
+	velocity    realx3 (0.0 -0.6 0.0);  
+	}
+
+	mixture
+	{
+	sphere 1;  
+	}
+}

tutorials/sphereGranFlow/RotaryAirLockValve/caseSetup/sphereShape

@@ -0,0 +1,19 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+objectName particleInsertion;
+objectType dicrionary;
+fileFormat  ASCII;
+/*---------------------------------------------------------------------------*/
+objectName 	sphereDict;
+objectType 	sphereShape;
+
+// names of shapes
+names 		(sphere);
+
+// diameter of shapes 
+diameters 	(0.005);
+
+// material names for shapes 
+materials	(sphereMat);

tutorials/sphereGranFlow/RotaryAirLockValve/settings/geometryDict

@@ -0,0 +1,66 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	geometryDict;
+objectType 	dictionary;
+fileFormat  ASCII;
+/*---------------------------------------------------------------------------*/
+
+// motion model: rotating object around an axis
+motionModel rotatingAxisMotion; 
+
+// information for rotatingAxisMotion motion model 
+rotatingAxisMotionInfo
+{
+	rotAxis 
+	{
+
+		// first point for the axis of rotation
+		p1 (0.561547 0.372714 0.000);
+
+		// second point for the axis of rotation
+		p2 (0.561547 0.372714 0.010);
+
+		// rotation speed (rad/s)
+		omega 2.1;
+
+		// Start time of Geometry Rotating (s)
+		startTime 1.25;
+		
+		// End time of Geometry Rotating (s)
+		endTime 7;
+	}
+}
+
+surfaces
+{
+	gear
+	{
+		// type of the wall
+		type 	 stlWall;
+
+		// file name in stl folder
+		file 	 gear.stl;
+
+		// material name of this wall
+		material wallMat;
+
+		// motion component name 
+		motion 	 rotAxis;
+	}
+	surfaces
+	{
+		// type of the wall
+		type 	 stlWall;
+
+		// file name in stl folder
+		file 	 surfaces.stl;
+
+		// material name of this wall
+		material wallMat;
+
+		// motion component name 
+		motion 	 none;
+}
+

tutorials/sphereGranFlow/RotaryAirLockValve/settings/particlesDict

@@ -0,0 +1,43 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	geometryDict;
+objectType 	dictionary;
+fileFormat  ASCII;
+/*---------------------------------------------------------------------------*/
+
+setFields
+{
+	defaultValue 
+	{
+		// linear velocity (m/s)
+		velocity 	realx3 	(0 0 0);
+
+		// 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	sphere;
+	}
+
+	selectors
+	{}
+}
+
+// positions particles 
+positionParticles
+{
+
+	// creates the required fields with zero particles (empty).
+	method empty;
+
+	// maximum number of particles in the simulation
+	maxNumberOfParticles 	50000;
+
+	// perform initial sorting based on morton code?
+	mortonSorting 		Yes;
+}

tutorials/sphereGranFlow/RotaryAirLockValve/settings/settingsDict

@@ -0,0 +1,48 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	geometryDict;
+objectType 	dictionary;
+fileFormat  ASCII;
+/*---------------------------------------------------------------------------*/
+
+run		rotatingValve;
+
+// time step for integration (s)
+dt			0.00001;
+
+// start time for simulation
+startTime		0;
+
+// end time for simulation
+endTime			7;
+
+// time interval for saving the simulation
+saveInterval	0.05;
+
+// maximum number of digits for time folder
+timePrecision	6;
+
+// gravity vector (m/s2) 
+g		(0 -9.8 0);
+
+/*
+	Simulation domain every particles that goes outside this domain is deleted.
+*/
+
+domain 
+{
+	min (0.397538  0.178212 0.00);
+	max (0.725537  0.600214 0.06);
+}
+
+// integration method 
+integrationMethod	AdamsBashforth3;
+
+// report timers?
+timersReport		Yes;
+
+// time interval for reporting timers
+timersReportInterval		  0.01;	                
+