Merge branch 'main' of github.com:PhasicFlow/phasicFlow into main

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 *\

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;