From 5e9f6e150c6d711d03c69eaa28aae51d8b44a493 Mon Sep 17 00:00:00 2001 From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com> Date: Fri, 17 Feb 2023 20:20:39 +0330 Subject: [PATCH 1/8] Added Picture to Tutorial --- tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md index 7e0d3258..d100fcfc 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md @@ -3,6 +3,13 @@ The problem is to simulate a Rotating Drum with **6** Baffles with the diameter * **12500** Particles with **4 mm** diameter * **7500** Particles with **5mm** diameter + + + + + + + ## Setting up the Case As it has been explained in the previous Cases, these Tutorials are based on text-based scripts. There are three parts in this case to study `caseSetup`, `setting` and `stl`. ## Particle Insertion From 4a591edb876ca9c548b42e66cb3d70602ae2b64c Mon Sep 17 00:00:00 2001 From: Omid Khosravi <123903048+omid-khosravi@users.noreply.github.com> Date: Fri, 17 Feb 2023 20:25:07 +0330 Subject: [PATCH 2/8] Added the Interaction Code block to Readme --- .../RotatingDrumWithBaffles/ReadMe.md | 37 ++++++++++++++++++- 1 file changed, 36 insertions(+), 1 deletion(-) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md index d100fcfc..c258eae9 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md @@ -50,7 +50,42 @@ diameters (0.004 0.005); materials (lightMat heavyMat); ``` In this Case we have two types of Particle with 4mm and 5mm diameters. -At the end of `caseSetup`, the interaction between the particles and the Shell of Rotating Drum is defined. You can see the Coefficients of the Interactions between the particles and shell of Rotating Drum in `interaction`. +At the end of `caseSetup`, the interaction between the particles and the Shell of Rotating Drum is defined. You can see the Coefficients of the Interactions between the particles and shell of Rotating Drum in `interaction`. Because we have 3 kind of interactions between these Particles and the Drum, we need to define a 3*3 Matrix. +```C++ + /* + 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); +``` ## Settings ### Geometry In the Settings folder the Specifications of our Rotating Drum and the information of rotating axis are brought. In this case we use two solid cylinders to keep our rotating drum isolated. This is to prevent particles, from being thrown out. From 3bcfc31fd23a197b71377b16e95f41c1c0da2593 Mon Sep 17 00:00:00 2001 From: Hamidreza Norouzi <72734524+hamidrezanorouzi@users.noreply.github.com> Date: Sun, 19 Feb 2023 23:55:44 +0330 Subject: [PATCH 3/8] image-source replaced --- tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md index c258eae9..df4fc2a3 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md @@ -5,7 +5,7 @@ The problem is to simulate a Rotating Drum with **6** Baffles with the diameter - + From aa149bcef79ba52856e11e71f15072217dd373d0 Mon Sep 17 00:00:00 2001 From: Hamidreza Norouzi <72734524+hamidrezanorouzi@users.noreply.github.com> Date: Mon, 20 Feb 2023 12:29:22 +0330 Subject: [PATCH 4/8] Update ReadMe.md --- .../RotatingDrumWithBaffles/ReadMe.md | 144 ++++++++++++------ 1 file changed, 99 insertions(+), 45 deletions(-) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md index df4fc2a3..54d15d41 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md @@ -1,20 +1,43 @@ # Problem Definition -The problem is to simulate a Rotating Drum with **6** Baffles with the diameter **0.24m** and the length **0.1m** rotating at **15 rad/s**. This Rotating Drum is filled with **20000** Particles.The timestep for integration is **0.00001 s**. There are 2 types of Particles in this Rotating Drum: -* **12500** Particles with **4 mm** diameter -* **7500** Particles with **5mm** diameter +The problem is to simulate a rotating drum with the diameter **0.24 m**, the length **0.1 m** and **6** Baffles, rotating at **15 rpm**. This drum is filled with **20000** Particles.The timestep for integration is **0.00001 s**. There are 2 types of Particles in this drum each are beining inserted during simulation to fill the drum. +* **12500** Particles with **4 mm** diameter, at the rate of 12500 particles/s for 1 sec. +* **7500** Particles with **5mm** diameter, at the rate of 7500 particles/s for 1 sec. - - +
+ a view of the drum while rotating +
+
+ +
-## Setting up the Case -As it has been explained in the previous Cases, these Tutorials are based on text-based scripts. There are three parts in this case to study `caseSetup`, `setting` and `stl`. +# Setting up the Case +As it has been explained in the previous cases, the simulation case setup is based on text-based scripts. Here, the simulation case setup are sotred in three folders: `caseSetup`, `setting` and `stl` (see the above folders). + +## Defining small and large particles +Then in the `caseSetup/sphereShape` the diameter and the material name of the particles are defined. Two sizes are defined: 4 and 5 mm. +```C++ +// names of shapes +names (smallSphere largeSphere); +// diameter of shapes (m) +diameters (0.004 0.005); +// material names for shapes +materials (lightMat heavyMat); +``` + + ## Particle Insertion -In this case we have two region for inserting our particles. In the both region we define rate of Insertion, start and end time of Insertion, coordinates of Insertion and radius of Insertion. -An example for the Right Layer Region of insertion of Particles is shown below. +In this case we have two region for inserting our particles. In the both region we define rate of insertion, start and end time of insertion, information for the volume of the space throught which particles are being inserted. The insertion phase in the simulation is performed between times 0 and 1 seconds. +For example, for the insertion region for inserting light particles is shown below. + +
+in caseSetup/particleInsertion file +
+ + ```C++ // Right Layer Region layerrightregion @@ -40,18 +63,15 @@ layerrightregion } } ``` -Then in the `sphereShape` the diameter and the material of our Particles are defined. -```C++ -// names of shapes -names (lightSphere heavySphere); -// diameter of shapes (m) -diameters (0.004 0.005); -// material names for shapes -materials (lightMat heavyMat); -``` -In this Case we have two types of Particle with 4mm and 5mm diameters. -At the end of `caseSetup`, the interaction between the particles and the Shell of Rotating Drum is defined. You can see the Coefficients of the Interactions between the particles and shell of Rotating Drum in `interaction`. Because we have 3 kind of interactions between these Particles and the Drum, we need to define a 3*3 Matrix. +## Interaction between particles and walls +In `caseSetup/interaction` file, material names and properties and interaction parameters are defined: interaction between the particles and the shell of rotating drum. Since we are defining 3 materials for simulation, the interaction matrix is 3x3, while we are only required to enter upper-triangle elements (interactions are symetric). + ```C++ +// a list of materials names +materials (lightMat heavyMat wallMat); +// density of materials [kg/m3] +densities (1000 1500 2500); + /* Property (lightMat-lightMat lightMat-heavyMat lightMat-wallMat heavyMat-heavyMat heavyMat-wallMat @@ -88,30 +108,63 @@ At the end of `caseSetup`, the interaction between the particles and the Shell o ``` ## Settings ### Geometry -In the Settings folder the Specifications of our Rotating Drum and the information of rotating axis are brought. In this case we use two solid cylinders to keep our rotating drum isolated. This is to prevent particles, from being thrown out. -For example the codes for the rear cylinder is brought below. +In the `settings/geometryDict` file, the geometry and axis of rotation is defined for the drum. The geometry is composed of a body, front and rear ends. + ```C++ -/*This is a Cylinder Wall at the rear of cylinder */ -CylinderRear1 +surfaces { - // 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; + 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 */ + rearEnd + { + // 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 */ + frontEnd + { + // 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; + } } ``` ### Rotating Axis Info -In this part of `geometryDict` the information of `rotating axis` and `velocity` of this Rotating Drum is defined. Also in purpose to settle down Particles after they were inserted we use a `startTime` and `endTime` function. This shows the start time of rotation. +In this part of `geometryDict` the information of rotating axis and speed of rotation are defined. The start of rotation is at 2 s. The first 2 seconds of simulation is for allowing particles to settle donw in the drum. + ```C++ rotatingAxisMotionInfo { @@ -130,9 +183,10 @@ rotatingAxisMotionInfo } } ``` -## Starting Simulation -To start Simulation we have to create our Particles at first. -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. +## 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. From 6fdf0d43f6c75b04fb6e477cb8da122f79590b22 Mon Sep 17 00:00:00 2001 From: Hamidreza Norouzi <72734524+hamidrezanorouzi@users.noreply.github.com> Date: Mon, 20 Feb 2023 12:30:34 +0330 Subject: [PATCH 5/8] Update geometryDict --- .../RotatingDrumWithBaffles/settings/geometryDict | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict index 5da00ec3..bb372bce 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict @@ -23,7 +23,7 @@ surfaces motion rotAxis; } /* This is a Cylinder Wall at the rear of cylinder */ - CylinderRear1 + rearEnd { // type of the wall type cylinderWall; @@ -41,7 +41,7 @@ surfaces motion rotAxis; } /* This a cylinder Wall at the front of Cylinder */ - CylinderFront2 + frontEnd { // type of the wall type cylinderWall; From d61c46f8794f9cc1bb69fb1fc62fe1d34b7d61ea Mon Sep 17 00:00:00 2001 From: Hamidreza Norouzi <72734524+hamidrezanorouzi@users.noreply.github.com> Date: Mon, 20 Feb 2023 12:31:12 +0330 Subject: [PATCH 6/8] Update particlesDict --- .../RotatingDrumWithBaffles/settings/particlesDict | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict index a7f58c61..4be2a5c3 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict @@ -16,7 +16,7 @@ setFields // rotational velocity (rad/s) rotVelocity realx3 (0 0 0); // name of the particle shape - shapeName word lightSphere; + shapeName word smallSphere; } selectors From aece11657d75d24b49fc4d126f2b9ded84f1cb0f Mon Sep 17 00:00:00 2001 From: Hamidreza Norouzi <72734524+hamidrezanorouzi@users.noreply.github.com> Date: Mon, 20 Feb 2023 12:32:10 +0330 Subject: [PATCH 7/8] Update sphereShape --- .../RotatingDrumWithBaffles/caseSetup/sphereShape | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/sphereShape b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/sphereShape index 1b94f4f0..714e1b9a 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/sphereShape +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/sphereShape @@ -6,7 +6,7 @@ objectName sphereDict; objectType sphereShape; /*---------------------------------------------------------------------------*/ // names of shapes -names (lightSphere heavySphere); +names (smallSphere largeSphere); // diameter of shapes (m) diameters (0.004 0.005); // material names for shapes From f8995445d1cdff6334c541c7e0addba3ab507695 Mon Sep 17 00:00:00 2001 From: Hamidreza Norouzi <72734524+hamidrezanorouzi@users.noreply.github.com> Date: Mon, 20 Feb 2023 12:32:47 +0330 Subject: [PATCH 8/8] Update particleInsertion --- .../RotatingDrumWithBaffles/caseSetup/particleInsertion | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion index 8a16d732..68712a02 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion @@ -47,7 +47,7 @@ Two layers of particles are packed one-by-one using 1 insertion steps. mixture { // mixture composition of inserted particles - lightSphere 1; + smallSphere 1; } } // Left Layer Region @@ -79,7 +79,7 @@ Two layers of particles are packed one-by-one using 1 insertion steps. mixture { - heavySphere 1; // only heavySphere + largeSphere 1; // only heavySphere } }