xfygogo
/
phasicFlow
mirror of https://github.com/PhasicFlow/phasicFlow.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Update Tote Blender ReadMe 

# Tote Blender ReadMe
* ASCII file format explanation.
* Code blocks were fixed.
* Particles number were changed from 20000 to 24000.
This commit is contained in:
Omid Khosravi 2023-04-20 11:35:38 +04:30
parent 61f5288e18
commit 1771704270
1 changed files with 188 additions and 110 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

202
tutorials/sphereGranFlow/toteblender/ReadMe.md
View File

@ -1,6 +1,6 @@
# Problem Definition # 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 blender. 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 **24000** 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 blender.
* **20000** particles with **4 mm** diameter, at the rate of 20000 particles/s for 1 sec. َAfter settling particles, this blender starts to rotate at t=**1s**. * **24000** particles with **5 mm** diameter, at the rate of 24000 particles/s for 1 sec. َAfter settling particles, this blender starts to rotate at t=**1s**. For better and faster performace in simulations where the number of particles is very large, the format of the files is saved as **ASCII**.
<html> <html>
<body> <body>
@ -19,12 +19,14 @@ As it has been explained in the previous cases, the simulation case setup is bas
## Defining particles ## Defining particles
Then in the `caseSetup/sphereShape` the diameter and the material name of the particles are defined. Then in the `caseSetup/sphereShape` the diameter and the material name of the particles are defined.
```C++ ```C++
// names of shapes // name of shapes
names (sphere1); names (sphere1);
// diameter of shapes (m) // diameter of shapes (m)
diameters (0.004); diameters (0.005);
// material names for shapes
materials (prop1); // material name for shapes
materials (solidProperty);
``` ```
## Particle Insertion ## Particle Insertion
In this case we have a region for ordering particles. These particles are placed in this blender. For example the script for the inserted particles is shown below. In this case we have a region for ordering particles. These particles are placed in this blender. For example the script for the inserted particles is shown below.
@ -37,38 +39,46 @@ in <b>caseSetup/particleInsertion</b> file
// positions particles // positions particles
positionParticles positionParticles
{ {
// ordered positioning // ordered positioning
method positionOrdered; method positionOrdered;
// maximum number of particles in the simulation
maxNumberOfParticles 40000; // maximum number of particles in the simulation
// perform initial sorting based on morton code? maxNumberOfParticles 25001;
// perform initial sorting based on morton code?
mortonSorting Yes; mortonSorting Yes;
// cylinder for positioning particles
// cylinder for positioning particles
cylinder cylinder
{ {
// Coordinates of top cylinderRegion (m,m,m) // Coordinates of top cylinderRegion (m,m,m)
p1 (0.05 0.0 0.12); p1 (0.0 0.0 0.09);
p2 (0.05 0.0 0.22);
// radius of cylinder p2 (0.0 0.0 0.21);
radius 0.066;
// radius of cylinder
radius 0.09;
} }
positionOrderedInfo positionOrderedInfo
{ {
// minimum space between centers of particles // minimum space between centers of particles
diameter 0.003; diameter 0.005;
// number of particles in the simulation
numPoints 20000; // number of particles in the simulation
// axis order for filling the space with particles numPoints 24000;
axisOrder (z y x);
// axis order for filling the space with particles
axisOrder (x y z);
} }
} }
``` ```
## Interaction between particles ## 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). In `caseSetup/interaction` file, material names and properties and interaction parameters are defined: interaction between the particles of Tote Blender. Since we are defining 1 material for simulation, the interaction matrix is 1x1 (interactions are symetric).
```C++ ```C++
// a list of materials names // a list of materials names
materials (prop1); materials (solidProperty);
// density of materials [kg/m3] // density of materials [kg/m3]
densities (1000.0); densities (1000.0);
@ -77,48 +87,62 @@ contactListType sortedContactList;
model model
{ {
contactForceModel nonLinearNonLimited; 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);
rollingFrictionModel normal;
/*
Property (solidProperty-solidProperty);
*/
// 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 ## Settings
### Geometry ### 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. In the `settings/geometryDict` file, the geometry and axis of rotation is defined for the blender. The geometry is composed of a cylinder inlet and outlet, cone shell top and down, a cylinder shell and enter and exit Gate.
```C++ ```C++
surfaces surfaces
{ {
topGate
topGate topGate
{ {
// type of wall // type of wall
type cylinderWall; type cylinderWall;
// begin point of cylinder axis // begin point of cylinder axis
p1 (0.0 0.0 0.299); p1 (0.0 0.0 0.3);
// end point of cylinder axis // end point of cylinder axis
p2 (0.0 0.0 0.3); p2 (0.0 0.0 0.301);
// radius at p1 // radius at p1
radius1 0.03; radius1 0.03;
// radius at p2 // radius at p2
radius2 0.0001; radius2 0.0001;
// material of wall // material of wall
material solidProperty; material solidProperty;
// motion component name // motion component name
motion axisOfRotation; motion axisOfRotation;
} }
@ -127,18 +151,25 @@ surfaces
{ {
// type of the wall // type of the wall
type cylinderWall; type cylinderWall;
// begin point of cylinder axis // begin point of cylinder axis
p1 (0.0 0.0 0.28); p1 (0.0 0.0 0.28);
// end point of cylinder axis // end point of cylinder axis
p2 (0.0 0.0 0.3); p2 (0.0 0.0 0.3);
// radius at p1 // radius at p1
radius1 0.03; radius1 0.03;
// radius at p2 // radius at p2
radius2 0.03; radius2 0.03;
// number of divisions // number of divisions
resolution 36; resolution 36;
// material name of this wall // material name of this wall
material prop1; material solidProperty;
// motion component name // motion component name
motion axisOfRotation; motion axisOfRotation;
} }
@ -147,18 +178,25 @@ surfaces
{ {
// type of the wall // type of the wall
type cylinderWall; type cylinderWall;
// begin point of cylinder axis // begin point of cylinder axis
p1 (0.0 0.0 0.2); p1 (0.0 0.0 0.2);
// end point of cylinder axis // end point of cylinder axis
p2 (0.0 0.0 0.28); p2 (0.0 0.0 0.28);
// radius at p1 // radius at p1
radius1 0.1; radius1 0.1;
// radius at p2 // radius at p2
radius2 0.03; radius2 0.03;
// number of divisions // number of divisions
resolution 36; resolution 36;
// material name of this wall // material name of this wall
material prop1; material solidProperty;
// motion component name // motion component name
motion axisOfRotation; motion axisOfRotation;
} }
@ -167,72 +205,109 @@ surfaces
{ {
// type of the wall // type of the wall
type cylinderWall; type cylinderWall;
// begin point of cylinder axis // begin point of cylinder axis
p1 (0.0 0.0 0.1); p1 (0.0 0.0 0.1);
// end point of cylinder axis // end point of cylinder axis
p2 (0.0 0.0 0.2); p2 (0.0 0.0 0.2);
// radius at p1 // radius at p1
radius1 0.1; radius1 0.1;
// radius at p2 // radius at p2
radius2 0.1; radius2 0.1;
// number of divisions // number of divisions
resolution 36; resolution 36;
// material name of this wall // material name of this wall
material prop1; material solidProperty;
// motion component name // motion component name
motion axisOfRotation; motion axisOfRotation;
} }
coneShelldown coneShelldown
{ {
// type of the wall // type of the wall
type cylinderWall; type cylinderWall;
// begin point of cylinder axis // begin point of cylinder axis
p1 (0.0 0.0 0.02); p1 (0.0 0.0 0.02);
// end point of cylinder axis // end point of cylinder axis
p2 (0.0 0.0 0.1); p2 (0.0 0.0 0.1);
// radius at p1 // radius at p1
radius1 0.03; radius1 0.03;
// radius at p2 // radius at p2
radius2 0.1; radius2 0.1;
// number of divisions // number of divisions
resolution 36; resolution 36;
// material name of this wall // material name of this wall
material prop1; material solidProperty;
// motion component name // motion component name
motion axisOfRotation; motion axisOfRotation;
} }
/*
This is a plane wall at the exit of silo
*/
bottomCylinder bottomCylinder
{ {
// type of the wall // type of the wall
type cylinderWall; type cylinderWall;
// begin point of cylinder axis // begin point of cylinder axis
p1 (0.0 0.0 0.0); p1 (0.0 0.0 0.0);
// end point of cylinder axis // end point of cylinder axis
p2 (0.0 0.0 0.02); p2 (0.0 0.0 0.02);
// radius at p1 // radius at p1
radius1 0.03; radius1 0.03;
// radius at p2 // radius at p2
radius2 0.03; radius2 0.03;
// number of divisions // number of divisions
resolution 36; resolution 36;
// material name of this wall // material name of this wall
material prop1; material solidProperty;
// motion component name // motion component name
motion axisOfRotation; motion axisOfRotation;
} }
exitGate exitGate
{ {
type planeWall;
p1 (-0.05 -0.05 0); // type of the wall
p2 (-0.05 0.05 0); type cylinderWall;
p3 ( 0.05 0.05 0);
p4 (0.05 -0.05 0); // begin point of cylinder axis
material prop1; p1 (0.0 0.0 -0.001);
// end point of cylinder axis
p2 (0.0 0.0 0.0);
// radius at p1
radius1 0.03;
// radius at p2
radius2 0.0001;
// number of divisions
resolution 36;
// material name of this wall
material solidProperty;
// motion component name
motion axisOfRotation; motion axisOfRotation;
} }
@ -247,10 +322,13 @@ rotatingAxisMotionInfo
axisOfRotation axisOfRotation
{ {
p1 (-0.1 0.0 0.15); // first point for the axis of rotation 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 p2 ( 0.1 0.0 0.15); // second point for the axis of rotation
omega 1.5708; // rotation speed ==> 15 rad/s omega 1.5708; // rotation speed ==> 15 rad/s
// Start time of Geometry Rotating (s) // Start time of Geometry Rotating (s)
startTime 1; startTime 0.5;
// End time of Geometry Rotating (s) // End time of Geometry Rotating (s)
endTime 9.5; endTime 9.5;
} }