From 99da1d5a209b25264762a34f1985f7dae9174f6f Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:47:19 +0330
Subject: [PATCH 01/10] Create toteblender folder
---
tutorials/sphereGranFlow/toteblender/toteblender | 1 +
1 file changed, 1 insertion(+)
create mode 100644 tutorials/sphereGranFlow/toteblender/toteblender
diff --git a/tutorials/sphereGranFlow/toteblender/toteblender b/tutorials/sphereGranFlow/toteblender/toteblender
new file mode 100644
index 00000000..e8689b24
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/toteblender
@@ -0,0 +1 @@
+hh
From 46385b695315db53133ff2ec81f88bf17931ac1c Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:53:46 +0330
Subject: [PATCH 02/10] Add tote blender files via Upload
---
.../sphereGranFlow/toteblender/cleanThisCase | 7 +++++++
.../sphereGranFlow/toteblender/runThisCase | 21 +++++++++++++++++++
2 files changed, 28 insertions(+)
create mode 100644 tutorials/sphereGranFlow/toteblender/cleanThisCase
create mode 100644 tutorials/sphereGranFlow/toteblender/runThisCase
diff --git a/tutorials/sphereGranFlow/toteblender/cleanThisCase b/tutorials/sphereGranFlow/toteblender/cleanThisCase
new file mode 100644
index 00000000..8a0ab919
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/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
+
+#------------------------------------------------------------------------------
diff --git a/tutorials/sphereGranFlow/toteblender/runThisCase b/tutorials/sphereGranFlow/toteblender/runThisCase
new file mode 100644
index 00000000..c48d71fe
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/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
+
+
+
+
+#------------------------------------------------------------------------------
From 8666e87f83ef00fb86acbca990dc05a528b78bb9 Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:54:33 +0330
Subject: [PATCH 03/10] Delete test file
---
tutorials/sphereGranFlow/toteblender/toteblender | 1 -
1 file changed, 1 deletion(-)
delete mode 100644 tutorials/sphereGranFlow/toteblender/toteblender
diff --git a/tutorials/sphereGranFlow/toteblender/toteblender b/tutorials/sphereGranFlow/toteblender/toteblender
deleted file mode 100644
index e8689b24..00000000
--- a/tutorials/sphereGranFlow/toteblender/toteblender
+++ /dev/null
@@ -1 +0,0 @@
-hh
From bc3409d1540d75386833a31880d189ad14520052 Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:55:33 +0330
Subject: [PATCH 04/10] create test
---
tutorials/sphereGranFlow/toteblender/caseSetup/test | 1 +
1 file changed, 1 insertion(+)
create mode 100644 tutorials/sphereGranFlow/toteblender/caseSetup/test
diff --git a/tutorials/sphereGranFlow/toteblender/caseSetup/test b/tutorials/sphereGranFlow/toteblender/caseSetup/test
new file mode 100644
index 00000000..8b137891
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/caseSetup/test
@@ -0,0 +1 @@
+
From 1b0667576bde9bb0bc2e5042c0f997bb4e011a26 Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:56:03 +0330
Subject: [PATCH 05/10] Add caseSetup files via Upload
---
.../toteblender/caseSetup/interaction | 61 +++++++++++++++++++
.../toteblender/caseSetup/particleInsertion | 13 ++++
.../toteblender/caseSetup/sphereShape | 13 ++++
3 files changed, 87 insertions(+)
create mode 100644 tutorials/sphereGranFlow/toteblender/caseSetup/interaction
create mode 100644 tutorials/sphereGranFlow/toteblender/caseSetup/particleInsertion
create mode 100644 tutorials/sphereGranFlow/toteblender/caseSetup/sphereShape
diff --git a/tutorials/sphereGranFlow/toteblender/caseSetup/interaction b/tutorials/sphereGranFlow/toteblender/caseSetup/interaction
new file mode 100644
index 00000000..f97d2f4f
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/caseSetup/interaction
@@ -0,0 +1,61 @@
+/* -------------------------------*- C++ -*--------------------------------- *\
+| phasicFlow File |
+| copyright: www.cemf.ir |
+\* ------------------------------------------------------------------------- */
+objectName interaction;
+objectType dicrionary;
+/* ------------------------------------------------------------------------- */
+// a list of materials names
+materials (prop1);
+// density of materials [kg/m3]
+densities (1000.0);
+
+contactListType sortedContactList;
+
+model
+{
+ contactForceModel nonLinearNonLimited;
+ rollingFrictionModel normal;
+ /*
+ Property (prop1-prop1);
+ */
+// Young modulus [Pa]
+ Yeff (1.0e6);
+// Shear modulus [Pa]
+ Geff (0.8e6);
+// Poisson's ratio [-]
+ nu (0.25);
+// coefficient of normal restitution
+ en (0.7);
+// coefficient of tangential restitution
+ et (1.0);
+// dynamic friction
+ mu (0.3);
+// rolling friction
+ mur (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;
+ }
+
+}
\ No newline at end of file
diff --git a/tutorials/sphereGranFlow/toteblender/caseSetup/particleInsertion b/tutorials/sphereGranFlow/toteblender/caseSetup/particleInsertion
new file mode 100644
index 00000000..d8f022c6
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/caseSetup/particleInsertion
@@ -0,0 +1,13 @@
+/* -------------------------------*- C++ -*--------------------------------- *\
+| phasicFlow File |
+| copyright: www.cemf.ir |
+\* ------------------------------------------------------------------------- */
+objectName particleInsertion;
+objectType dicrionary;
+/* ------------------------------------------------------------------------- */
+// is insertion active?
+active no;
+// not implemented for yes
+collisionCheck No;
+
+
diff --git a/tutorials/sphereGranFlow/toteblender/caseSetup/sphereShape b/tutorials/sphereGranFlow/toteblender/caseSetup/sphereShape
new file mode 100644
index 00000000..ff2df8d1
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/caseSetup/sphereShape
@@ -0,0 +1,13 @@
+/* -------------------------------*- C++ -*--------------------------------- *\
+| phasicFlow File |
+| copyright: www.cemf.ir |
+\* ------------------------------------------------------------------------- */
+objectName sphereDict;
+objectType sphereShape;
+/* ------------------------------------------------------------------------- */
+// names of shapes
+names (sphere1);
+// diameter of shapes (m)
+diameters (0.004);
+// material names for shapes
+materials (prop1);
\ No newline at end of file
From 4ecc667a10ad36fd72dd44ad15cd493c54bb53af Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:56:41 +0330
Subject: [PATCH 06/10] Delete test
---
tutorials/sphereGranFlow/toteblender/caseSetup/test | 1 -
1 file changed, 1 deletion(-)
delete mode 100644 tutorials/sphereGranFlow/toteblender/caseSetup/test
diff --git a/tutorials/sphereGranFlow/toteblender/caseSetup/test b/tutorials/sphereGranFlow/toteblender/caseSetup/test
deleted file mode 100644
index 8b137891..00000000
--- a/tutorials/sphereGranFlow/toteblender/caseSetup/test
+++ /dev/null
@@ -1 +0,0 @@
-
From 05328a395171d805e57c20194d8caf2e1cc5e406 Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:57:13 +0330
Subject: [PATCH 07/10] Create test
---
tutorials/sphereGranFlow/toteblender/settings/test | 1 +
1 file changed, 1 insertion(+)
create mode 100644 tutorials/sphereGranFlow/toteblender/settings/test
diff --git a/tutorials/sphereGranFlow/toteblender/settings/test b/tutorials/sphereGranFlow/toteblender/settings/test
new file mode 100644
index 00000000..6a69f920
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/settings/test
@@ -0,0 +1 @@
+f
From 0dcff307a6181f3b64f0f6834a2c9aca4443c95a Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:57:44 +0330
Subject: [PATCH 08/10] Add settings files via Upload
---
.../toteblender/settings/geometryDict | 151 ++++++++++++++++++
.../toteblender/settings/particlesDict | 59 +++++++
.../toteblender/settings/settingsDict | 34 ++++
3 files changed, 244 insertions(+)
create mode 100644 tutorials/sphereGranFlow/toteblender/settings/geometryDict
create mode 100644 tutorials/sphereGranFlow/toteblender/settings/particlesDict
create mode 100644 tutorials/sphereGranFlow/toteblender/settings/settingsDict
diff --git a/tutorials/sphereGranFlow/toteblender/settings/geometryDict b/tutorials/sphereGranFlow/toteblender/settings/geometryDict
new file mode 100644
index 00000000..49d486bd
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/settings/geometryDict
@@ -0,0 +1,151 @@
+/* -------------------------------*- C++ -*--------------------------------- *\
+| phasicFlow File |
+| copyright: www.cemf.ir |
+\* ------------------------------------------------------------------------- */
+objectName geometryDict;
+objectType dictionary;
+/* ------------------------------------------------------------------------- */
+// motion model: rotating object around an axis
+motionModel rotatingAxisMotion;
+
+surfaces
+{
+ enterGate
+ {
+ // type of wall
+ type planeWall;
+ // coords of wall
+ p1 (-0.05 -0.05 0.3);
+ p2 (-0.05 0.05 0.3);
+ p3 ( 0.05 0.05 0.3);
+ p4 (0.05 -0.05 0.3);
+ // material of wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ cylinderinlet
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.28);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.3);
+ // radius at p1
+ radius1 0.03;
+ // radius at p2
+ radius2 0.03;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ coneShelltop
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.2);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.28);
+ // radius at p1
+ radius1 0.1;
+ // radius at p2
+ radius2 0.03;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ cylinderShell
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.1);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.2);
+ // radius at p1
+ radius1 0.1;
+ // radius at p2
+ radius2 0.1;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ coneShelldown
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.02);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.1);
+ // radius at p1
+ radius1 0.03;
+ // radius at p2
+ radius2 0.1;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+ /*
+ This is a plane wall at the exit of silo
+ */
+
+ cylinderoutlet
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.0);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.02);
+ // radius at p1
+ radius1 0.03;
+ // radius at p2
+ radius2 0.03;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+ exitGate
+ {
+ type planeWall;
+ p1 (-0.05 -0.05 0);
+ p2 (-0.05 0.05 0);
+ p3 ( 0.05 0.05 0);
+ p4 (0.05 -0.05 0);
+ material prop1;
+ motion rotAxis;
+ }
+
+}
+// information for rotatingAxisMotion motion model
+rotatingAxisMotionInfo
+{
+ rotAxis
+ {
+ p1 (-0.1 0.0 0.15); // first point for the axis of rotation
+ p2 (0.1 0.0 0.15); // second point for the axis of rotation
+ omega 3; // rotation speed (rad/s)
+ }
+}
diff --git a/tutorials/sphereGranFlow/toteblender/settings/particlesDict b/tutorials/sphereGranFlow/toteblender/settings/particlesDict
new file mode 100644
index 00000000..49a26e76
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/settings/particlesDict
@@ -0,0 +1,59 @@
+/* -------------------------------*- 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 sphere1;
+ }
+
+ selectors
+ {}
+}
+
+// positions particles
+positionParticles
+{
+// ordered positioning
+ method positionOrdered;
+// maximum number of particles in the simulation
+ maxNumberOfParticles 40000;
+// perform initial sorting based on morton code?
+ mortonSorting Yes;
+// box for positioning particles
+ box
+ {
+// lower corner point of the box
+ min (-0.06 -0.06 0.08);
+// upper corner point of the box
+ max (0.06 0.06 0.18);
+ }
+
+ positionOrderedInfo
+ {
+// minimum space between centers of particles
+ diameter 0.004;
+// number of particles in the simulation
+ numPoints 20000;
+// axis order for filling the space with particles
+ axisOrder (z y x);
+ }
+}
diff --git a/tutorials/sphereGranFlow/toteblender/settings/settingsDict b/tutorials/sphereGranFlow/toteblender/settings/settingsDict
new file mode 100644
index 00000000..28b9a552
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/settings/settingsDict
@@ -0,0 +1,34 @@
+/* -------------------------------*- C++ -*--------------------------------- *\
+| phasicFlow File |
+| copyright: www.cemf.ir |
+\* ------------------------------------------------------------------------- */
+objectName settingsDict;
+objectType dictionary;;
+/*---------------------------------------------------------------------------*/
+run toteBlender;
+// 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.1;
+// 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.5 -0.5 -0.5);
+ max (0.5 0.5 0.5);
+}
+// integration method
+integrationMethod AdamsBashforth2;
+// report timers?
+timersReport Yes;
+// time interval for reporting timers
+timersReportInterval 0.01;
From d4dc0d7e80663a3a3d7de210ce800ffdc1da1710 Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:58:15 +0330
Subject: [PATCH 09/10] Delete test
---
tutorials/sphereGranFlow/toteblender/settings/test | 1 -
1 file changed, 1 deletion(-)
delete mode 100644 tutorials/sphereGranFlow/toteblender/settings/test
diff --git a/tutorials/sphereGranFlow/toteblender/settings/test b/tutorials/sphereGranFlow/toteblender/settings/test
deleted file mode 100644
index 6a69f920..00000000
--- a/tutorials/sphereGranFlow/toteblender/settings/test
+++ /dev/null
@@ -1 +0,0 @@
-f
From 04d7ecfe9e91155945c541ec28ef6fae0bb746f9 Mon Sep 17 00:00:00 2001
From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com>
Date: Sun, 26 Mar 2023 18:58:42 +0330
Subject: [PATCH 10/10] Add ReadMe file via Upload
---
.../sphereGranFlow/toteblender/ReadMe.md | 244 ++++++++++++++++++
1 file changed, 244 insertions(+)
create mode 100644 tutorials/sphereGranFlow/toteblender/ReadMe.md
diff --git a/tutorials/sphereGranFlow/toteblender/ReadMe.md b/tutorials/sphereGranFlow/toteblender/ReadMe.md
new file mode 100644
index 00000000..4a9ff9bf
--- /dev/null
+++ b/tutorials/sphereGranFlow/toteblender/ReadMe.md
@@ -0,0 +1,244 @@
+# Problem Definition
+The problem is to simulate a double pedestal tote blender with the diameter **0.03 m** and **0.1 m** respectively, the length **0.3 m**, rotating at **28 rpm**. This blender is filled with **20000** Particles. The timestep for integration is **0.00001 s**. There is one type of Particle in this blender that are being inserted during simulation to fill the drum.
+* **20000** particles with **4 mm** diameter, at the rate of 20000 particles/s for 1 sec.
+
+
+
+
+ a view of the tote-blender while rotating
+
+
+
+
+in caseSetup/particleInsertion file
+
+
+```C++
+// positions particles
+positionParticles
+{
+// ordered positioning
+ method positionOrdered;
+// maximum number of particles in the simulation
+ maxNumberOfParticles 40000;
+// perform initial sorting based on morton code?
+ mortonSorting Yes;
+// box for positioning particles
+ box
+ {
+// lower corner point of the box
+ min (-0.06 -0.06 0.08);
+// upper corner point of the box
+ max (0.06 0.06 0.18);
+ }
+```
+ ## Interaction between particles
+ In `caseSetup/interaction` file, material names and properties and interaction parameters are defined: interaction between the particles of rotating drum. Since we are defining 1 material for simulation, the interaction matrix is 1x1 (interactions are symetric).
+```C++
+ // a list of materials names
+materials (prop1);
+// density of materials [kg/m3]
+densities (1000.0);
+
+contactListType sortedContactList;
+
+model
+{
+ contactForceModel nonLinearNonLimited;
+ rollingFrictionModel normal;
+ /*
+ Property (prop1-prop1);
+ */
+// Young modulus [Pa]
+ Yeff (1.0e6);
+// Shear modulus [Pa]
+ Geff (0.8e6);
+// Poisson's ratio [-]
+ nu (0.25);
+// coefficient of normal restitution
+ en (0.7);
+// coefficient of tangential restitution
+ et (1.0);
+// dynamic friction
+ mu (0.3);
+// rolling friction
+ mur (0.1);
+
+}
+```
+## Settings
+### Geometry
+In the `settings/geometryDict` file, the geometry and axis of rotation is defined for the drum. The geometry is composed of a cylinder inlet and outlet, cone shell top and down, a cylinder shell and enter and exit Gate.
+```C++
+surfaces
+{
+ enterGate
+ {
+ // type of wall
+ type planeWall;
+ // coords of wall
+ p1 (-0.05 -0.05 0.3);
+ p2 (-0.05 0.05 0.3);
+ p3 ( 0.05 0.05 0.3);
+ p4 (0.05 -0.05 0.3);
+ // material of wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ cylinderinlet
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.28);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.3);
+ // radius at p1
+ radius1 0.03;
+ // radius at p2
+ radius2 0.03;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ coneShelltop
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.2);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.28);
+ // radius at p1
+ radius1 0.1;
+ // radius at p2
+ radius2 0.03;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ cylinderShell
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.1);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.2);
+ // radius at p1
+ radius1 0.1;
+ // radius at p2
+ radius2 0.1;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+
+ coneShelldown
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.02);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.1);
+ // radius at p1
+ radius1 0.03;
+ // radius at p2
+ radius2 0.1;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+ /*
+ This is a plane wall at the exit of silo
+ */
+
+ cylinderoutlet
+ {
+ // type of the wall
+ type cylinderWall;
+ // begin point of cylinder axis
+ p1 (0.0 0.0 0.0);
+ // end point of cylinder axis
+ p2 (0.0 0.0 0.02);
+ // radius at p1
+ radius1 0.03;
+ // radius at p2
+ radius2 0.03;
+ // number of divisions
+ resolution 36;
+ // material name of this wall
+ material prop1;
+ // motion component name
+ motion rotAxis;
+ }
+ exitGate
+ {
+ type planeWall;
+ p1 (-0.05 -0.05 0);
+ p2 (-0.05 0.05 0);
+ p3 ( 0.05 0.05 0);
+ p4 (0.05 -0.05 0);
+ material prop1;
+ motion rotAxis;
+ }
+
+}
+```
+### Rotating Axis Info
+In this part of `geometryDict` the information of rotating axis and speed of rotation are defined. Unlike the previous cases, the rotation of this blender starts at time=**0 s**.
+```C++
+rotatingAxisMotionInfo
+{
+ rotAxis
+ {
+ p1 (-0.1 0.0 0.15); // first point for the axis of rotation
+ p2 (0.1 0.0 0.15); // second point for the axis of rotation
+ omega 3; // rotation speed (rad/s)
+ }
+}
+```
+## Performing Simulation
+To perform simulations, enter the following commands one after another in the terminal.
+
+Enter `$ particlesPhasicFlow` command to create the initial fields for particles.
+Enter `$ geometryPhasicFlow` command to create the Geometry.
+At last, enter `$ sphereGranFlow` command to start the simulation.
+After finishing the simulation, you can use `$ pFlowtoVTK` to convert the results into vtk format storred in ./VTK folder.
\ No newline at end of file