Merge pull request #73 from saharsiahtiri/main

Updated and added V-blender tutorial

tutorials/sphereGranFlow/V-blender/caseSetup/interaction

@@ -0,0 +1,74 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+objectName interaction;
+objectType dicrionary;
+/* ------------------------------------------------------------------------- */
+
+// a list of materials names
+materials      (wallMat lightMat);
+
+// density of materials [kg/m3]
+densities      (1000 1000);
+
+contactListType   sortedContactList; 
+
+model
+{
+   contactForceModel nonLinearNonLimited;
+   rollingFrictionModel normal;
+
+   // 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.97 0.85
+            0.97);
+
+   // dynamic friction
+   mu (0.65 0.35 
+            0.65);
+
+   // 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 20 timesteps, update neighbor list
+      updateFrequency 20;
+      
+      // bounding box size to particle diameter (max) 
+      sizeRatio      1.1;
+   }
+
+   cellMappingInfo
+   {
+      // each 20 timesteps, update neighbor list
+      updateFrequency 20;
+      
+      // bounding box for particle-wall search (> 0.5)
+      cellExtent     0.7;
+   }
+
+}

tutorials/sphereGranFlow/V-blender/caseSetup/particleInsertion

@@ -0,0 +1,101 @@
+/* -------------------------------*- 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 tho insertion steps. 
+*/
+
+rightregion
+{
+   // type of insertion region
+   type		  cylinderRegion;
+   
+   // insertion rate (particles/s)
+   rate 	  10000;
+   
+   // start time of small particles insertion (s)
+   startTime 	  0.0;
+   
+   // end time of small particles insertion (s)
+   endTime   	  1;
+   
+   // time interval of small particles insertion (s)
+   interval       0.025;  
+
+   cylinderRegionInfo 
+   {
+        // radius of cylinder (m)
+   	radius 	0.1;
+   	
+   	// coordinates of center of both ends of the insertion cylinder on 
+   	// the right side of the V-blender (m,m,m)
+   	p1 (0.0950615  0.12 0.5011585);
+   	p2 (0.1150615  0.12 0.4811585);
+   }
+   
+   setFields
+   {
+     // initial velocity of inserted particles
+     velocity    realx3 (1.2 0.0 -1.2); 
+   }
+   
+   mixture
+   {
+      // first layer of inserted particles
+      smallSphere 1;
+      
+   }
+}
+
+leftregion
+{
+   // type of insertion region
+   type         cylinderRegion;
+   
+   // insertion rate (particles/s)
+   rate         10000;
+   
+   // start time of large particles insertion (s)
+   startTime    1.5;
+   
+   // end time of large particles insertion (s)
+   endTime      2.5;
+   
+   // time interval of large particles insertion (s)
+   interval     0.025;
+
+   cylinderRegionInfo 
+   {
+      // radius of cylinder (m)
+      radius   0.1;
+      
+      // coordinates of center of both ends of the insertion cylinder on
+      // the left side of the V-blender (m,m,m)      
+      p1 ( 0.7562545  0.12 0.50079);
+      p2 ( 0.7362545  0.12 0.48079);
+   }
+   
+   setFields
+   {
+     // initial velocity of inserted particles
+     velocity    realx3 (-1.2 0.0 -1.2);
+   }
+   
+   mixture
+   {
+      // second layer of inserted particles
+      largeSphere 1;
+      
+   }
+}

tutorials/sphereGranFlow/V-blender/caseSetup/sphereShape

@@ -0,0 +1,16 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */ 
+objectName 	sphereDict;
+objectType 	sphereShape;
+/* ------------------------------------------------------------------------- */
+
+// names of particles
+names 		(smallSphere largeSphere);
+
+// diameter of particles
+diameters 	(0.01 0.0101);
+
+// material names for particles
+materials	(lightMat lightMat); 

tutorials/sphereGranFlow/V-blender/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/V-blender/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/V-blender/settings/geometryDict

@@ -0,0 +1,53 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	geometryDict;
+objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
+
+// motion model: rotating object around an axis
+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;
+	}
+
+	
+}
+
+// information for rotatingAxisMotion motion model 
+rotatingAxisMotionInfo
+{
+	rotAxis 
+	{
+	        // first point for the axis of rotation
+		p1 (0.128228 0.116446 0.297901);
+		
+		// second point for the axis of rotation
+		p2 (0.722596 0.116459 0.297901);
+		
+		// rotation speed (rad/s)
+		omega 3.14;
+		
+		// start time of rotation
+		startTime 3;
+		
+		// end time of rotation
+		endTime 10;
+	}
+}

tutorials/sphereGranFlow/V-blender/settings/particlesDict

@@ -0,0 +1,47 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	particlesDict;
+objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
+
+setFields
+{
+	/*
+	Default value for fields defined for particles
+	These fields should always be defined for simulations with 
+	spherical particles.
+	*/
+
+	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	smallSphere;
+	}
+
+	selectors
+	{}
+}
+
+// positions particles 
+positionParticles
+{
+        // no particle is initially in the simulation
+	method empty;
+
+        // maximum number of particles in the simulation
+	maxNumberOfParticles 25000;
+	
+	// perform initial sorting based on morton code?
+	mortonSorting Yes;
+}

tutorials/sphereGranFlow/V-blender/settings/settingsDict

@@ -0,0 +1,47 @@
+/* -------------------------------*- C++ -*--------------------------------- *\ 
+|  phasicFlow File                                                            | 
+|  copyright: www.cemf.ir                                                     | 
+\* ------------------------------------------------------------------------- */  
+objectName 	settingsDict;
+objectType 	dictionary;
+/* ------------------------------------------------------------------------- */
+
+run 		rotatingVblender;
+
+// 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 0 -9.8);
+
+/*
+	Simulation domain 
+	every particles that goes outside this domain is deleted.
+*/
+
+domain 
+{
+	min (-0.1 -0.4 0);
+	max (0.86 0.6 0.6);
+}
+
+// integration method
+integrationMethod 	AdamsBashforth2;
+
+// report timers?
+timersReport 		Yes;
+
+// time interval for reporting timers
+timersReportInterval   	0.01;