Merge pull request #72 from omid-khosravi/main

Rotating Drum with Baffles Tutorial

tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.txt

@@ -0,0 +1,10 @@
+To run this case please follow the below Steps:
+1- Open a terminal in the RotatingDrumWithBaffles
+2- Enter the code ./runThisCase to run the Case
+3- After finishing the run time please type pFlowtoVTK to convert data to VTK format
+4- After Post-Producting and viewing the case please enter ./cleanThisCase to delete data
+
+
+For more questions please contact:
+omidava.khosravi@gmail.com
+<3

tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/interaction

@@ -0,0 +1,76 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+objectName interaction;
+objectType dicrionary;
+/*----------------------------------------------------------------------------*/
+// a list of materials names
+materials   (lightMat heavyMat wallMat); 
+// density of materials [kg/m3]
+densities   (1000 1500 2500);     
+
+contactListType   sortedContactList;
+
+model
+{
+   contactForceModel nonLinearLimited;
+   rollingFrictionModel normal;
+
+   /*
+   Property (lightMat-lightMat   lightMat-heavyMat    lightMat-wallMat
+                                 heavyMat-heavyMat    heavyMat-wallMat
+                                                      wallMat-wallMat );
+   */
+// Young modulus [Pa]
+   Yeff (1.0e6 1.0e6 1.0e6  
+               1.0e6 1.0e6
+                     1.0e6);
+// Shear modulus [Pa]   
+   Geff (0.8e6 0.8e6 0.8e6  
+               0.8e6 0.8e6
+                     0.8e6);
+// Poisson's ratio [-]
+   nu    (0.25 0.25  0.25  
+               0.25  0.25
+                     0.25);
+// coefficient of normal restitution
+   en (0.97  0.97    0.85   
+             0.97    0.85
+                     1.00);
+// coefficient of tangential restitution
+   et (1.0   1.0  1.0    
+             1.0  1.0
+                  1.0);
+// dynamic friction 
+   mu (0.65   0.65 0.35   
+              0.65 0.35
+                   0.35);
+// rolling friction 
+   mur (0.1  0.1 0.1    
+             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 20 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/RotatingDrumWithBaffles/caseSetup/particleInsertion

@@ -0,0 +1,86 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+objectName particleInsertion;
+objectType dicrionary;
+/*---------------------------------------------------------------------------*/
+
+// is insertion active?
+active  yes;          
+// not implemented for yes
+collisionCheck No;
+
+/*
+Two layers of particles are packed one-by-one using 1 insertion steps. 
+*/
+
+// Right Layer Region
+   layerrightregion 
+{
+// type of insertion region
+   type	          cylinderRegion;
+// insertion rate (particles/s)    
+   rate 	  12500;
+// Start time of LightParticles insertion (s)
+   startTime 	  0; 
+// End time of LightParticles insertion (s)      
+   endTime   	  1;
+// Time Interval of LightParticles insertion (s)
+   interval       0.025; 
+
+   cylinderRegionInfo 
+   {
+// Coordinates of cylinderRegion (m,m,m)
+   	p2 (-0.15 0.25 0.05);
+   	p1 (-0.15 0.24	0.05);
+// radius of cylinder (m)
+   	radius 0.035;
+   }
+   
+   setFields
+   {
+// initial velocity of inserted particles
+     velocity    realx3 (0.0 -0.6 0.0);  
+   }
+   
+   mixture
+   {
+// mixture composition of inserted particles
+      lightSphere 1;    
+   }
+}
+// Left Layer Region
+   layerleftregion
+{
+   type         cylinderRegion;
+//Left Region Particles Insertion Rate (particles/s)
+   rate         7500;
+// Start time of LightParticles insertion (s)
+   startTime    0;
+// End time of LightParticles insertion (s)  
+   endTime      1;
+// Time Interval of LightParticles insertion (s)
+   interval     0.025;
+
+   cylinderRegionInfo 
+   {
+// Coordinates of cylinderRegion (m,m,m)
+      p2 (-0.23	 0.25	0.05); 
+      p1 (-0.23	 0.24	0.05);
+// radius of cylinder (m)      
+      radius 0.035;
+   }
+   
+   setFields
+   {
+     velocity    realx3 (0.0 -0.6 0.0);
+   }
+   
+   mixture
+   {
+      heavySphere 1; // only heavySphere
+      
+   }
+}
+

tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/sphereShape

@@ -0,0 +1,13 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+objectName 	sphereDict;
+objectType 	sphereShape;
+/*---------------------------------------------------------------------------*/
+// names of shapes 
+names 		(lightSphere heavySphere);
+// diameter of shapes (m) 	
+diameters 	(0.004 0.005);
+// material names for shapes 			
+materials	(lightMat heavyMat);		

tutorials/sphereGranFlow/RotatingDrumWithBaffles/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 
+
+#------------------------------------------------------------------------------

tutorials/sphereGranFlow/RotatingDrumWithBaffles/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
+
+
+
+
+#------------------------------------------------------------------------------

tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict

@@ -0,0 +1,78 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	geometryDict;
+objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
+
+// motion model: all surfaces are fixed
+motionModel rotatingAxisMotion; 
+
+surfaces
+{
+	body
+	{
+		// type of the wall
+		type 	 stlWall;
+		// file name in stl folder  	
+		file 	 Body.stl;
+		// material name of this wall		
+		material wallMat;
+		// motion component name   
+		motion 	 rotAxis;		 
+	}
+	/*	This is a Cylinder Wall at the rear of cylinder	*/
+	CylinderRear1
+	{
+		// type of the wall
+		type cylinderWall;
+		// first point for the axis of rotation			
+		p1 (-0.1974  0.2269  -0.001);
+		// second point for the axis of rotation	 
+		p2 (-0.1974  0.2269   0.0);
+		// Radius of p1	
+		radius1 0.0001;
+		// Radius of p2
+		radius2 0.12;
+		// material name of the wall
+		material wallMat;
+		// motion component name          
+		motion rotAxis;			 
+	}
+	/*	This a cylinder Wall at the front of Cylinder */
+	CylinderFront2
+	{
+		// type of the wall
+		type cylinderWall;
+		// first point for the axis of rotation			
+		p1 (-0.1974  0.2269   0.0989);
+		// second point for the axis of rotation	 
+		p2 (-0.1974  0.2269   0.0990);	
+		// Radius of p1
+		radius1 0.0001;
+		// Radius of p2
+		radius2 0.12;
+		// material name of the wall
+		material wallMat; 
+		// motion component name          
+		motion rotAxis;			 
+	}
+}
+	// information for rotatingAxisMotion motion model 
+rotatingAxisMotionInfo
+{
+	rotAxis 
+	{
+		// first point for the axis of rotation
+		p1 (-0.1974  0.2269  0);
+		// second point for the axis of rotation	 
+		p2 (-0.1974  0.2269  0.1);
+		// rotation speed (rad/s) => 15 rpm	
+		omega 2.38733;
+		// Start time of Geometry Rotating 		
+		startTime 2;
+		// End time of Geometry Rotating
+		endTime 9.5;
+	}
+}

tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict

@@ -0,0 +1,36 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	particlesDict;
+objectType 	dictionary;
+/*---------------------------------------------------------------------------*/
+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	lightSphere; 	
+	}
+
+	selectors
+	{}
+}
+
+// positions particles 
+positionParticles
+{
+	// creates the required spaces with zero particles (empty).
+	method empty;     		 
+	// maximum number of particles in the simulation
+	maxNumberOfParticles 	20050; 
+	// perform initial sorting based on morton code?
+	mortonSorting 	       	  Yes;       
+
+}

tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/settingsDict

@@ -0,0 +1,36 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+
+objectName 	settingsDict;
+objectType 	dictionary;;
+/*---------------------------------------------------------------------------*/
+run 			RotatingDrumwithBaffles;
+
+// time step for integration (s)
+	dt 			 0.00001;
+// start time for simulation 	
+	startTime 	   	   0;
+// end time for simulation 
+	endTime 		  10;	 		
+// 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.328 0.095 -0.025);
+	max (-0.068  0.355 0.125);
+}
+// integration method 
+	integrationMethod 		AdamsBashforth3; 	
+// report timers?
+	timersReport 			Yes;  	
+// time interval for reporting timers
+	timersReportInterval   	0.01;	

tutorials/sphereGranFlow/RotatingDrumWithBaffles/stl/note

@@ -0,0 +1 @@
+/* The STL files should be Placed in this Folder *\