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

updated V1.0 RotaryAirLock

This commit is contained in:
ramin1728 2024-12-28 13:09:06 +03:30
parent 6179f6ea3b
commit 278afc8ab6
8 changed files with 216 additions and 160 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
tutorials/sphereGranFlow/RotaryAirLockValve/ReadMe.md
View File

@ -3,6 +3,7 @@ The problem is to simulate a Rotary Air-Lock Valve. The external diameter of rot
* **28000** particles with **5 mm** diameter are inserted into the valve with the rate of **4000 particles/s**. * **28000** particles with **5 mm** diameter are inserted into the valve with the rate of **4000 particles/s**.
* The rotor starts its ortation at t = 1.25 s at the rate of 2.1 rad/s. * The rotor starts its ortation at t = 1.25 s at the rate of 2.1 rad/s.
<html> <html>
<body> <body>
<div align="center"><b> <div align="center"><b>

84
tutorials/sphereGranFlow/RotaryAirLockValve/caseSetup/interaction
View File

@ -6,78 +6,56 @@ objectName interaction;
objectType dicrionary; objectType dicrionary;
fileFormat ASCII; fileFormat ASCII;
/*---------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------*/
materials (sphereMat wallMat); // a list of materials names
// a list of materials names densities (1000 2500); // density of materials [kg/m3]
materials (sphereMat wallMat);
// density of materials [kg/m3]
densities (1000 2500);
contactListType sortedContactList; contactListType sortedContactList;
contactSearch
{
method NBS; // method for broad search
updateInterval 10;
sizeRatio 1.1;
cellExtent 0.55;
adjustableBox Yes;
}
model model
{ {
contactForceModel nonLinearNonLimited; contactForceModel nonLinearNonLimited;
rollingFrictionModel normal; rollingFrictionModel normal;
/* /*
Property (sphereMat-sphereMat sphereMat-wallMat Property (sphereMat-sphereMat sphereMat-wallMat
wallMat-wallMat);
wallMat-wallMat);
*/ */
// Young modulus [Pa] Yeff (1.0e6 1.0e6 // Young modulus [Pa]
Yeff (1.0e6 1.0e6
1.0e6); 1.0e6);
// Shear modulus [Pa] Geff (0.8e6 0.8e6 // Shear modulus [Pa]
Geff (0.8e6 0.8e6
0.8e6); 0.8e6);
// Poisson's ratio [-] nu (0.25 0.25 // Poisson's ratio [-]
nu (0.25 0.25 0.25);
0.25);
// coefficient of normal restitution en (0.7 0.8 // coefficient of normal restitution
en (0.7 0.8 1.0);
1.0);
// coefficient of tangential restitution et (1.0 1.0 // coefficient of tangential restitution
et (1.0 1.0 1.0);
1.0);
// dynamic friction mu (0.3 0.35 // dynamic friction
mu (0.3 0.35 0.35);
0.35);
// rolling friction mur (0.1 0.1 // rolling friction
mur (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 10 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;
}
}

37
tutorials/sphereGranFlow/RotaryAirLockValve/caseSetup/particleInsertion
View File

@ -6,45 +6,38 @@ objectName particleInsertion;
objectType dicrionary; objectType dicrionary;
fileFormat ASCII; fileFormat ASCII;
/*---------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------*/
active yes; // is insertion active -> yes or no
// is insertion active? checkForCollision No; // is checked -> yes or no
active yes;
// not implemented for yes
collisionCheck No;
/* /*
one region is considered for inserting particles. one region is considered for inserting particles.
*/ */
topRegion topRegion
{ {
regionType box; // type of insertion region
// type of insertion region timeControl simulationTime;
type boxRegion;
// insertion rate (particles/s) rate 4000; // insertion rate (particles/s)
rate 4000;
// Start time of Particles insertion (s) startTime 0; // Start time of Particles insertion (s)
startTime 0;
// End time of Particles insertion (s) endTime 7; // End time of Particles insertion (s)
endTime 7;
// Time Interval between each insertion (s) insertionInterval 0.025; // Time Interval between each insertion (s)
interval 0.025;
// Coordinates of BoxRegion (m,m,m) boxInfo // Coordinates of BoxRegion (m,m,m)
boxRegionInfo
{ {
min ( 0.48 0.58 0.01 ); // (m,m,m) min ( 0.48 0.58 0.01 ); // (m,m,m)
max ( 0.64 0.59 0.05 ); // (m,m,m)
max ( 0.64 0.59 0.05 ); // (m,m,m)
} }
setFields setFields
{ {
// initial velocity of inserted particles velocity realx3 (0.0 -0.6 0.0); // initial velocity of inserted particles
velocity realx3 (0.0 -0.6 0.0);
} }
mixture mixture

13
tutorials/sphereGranFlow/RotaryAirLockValve/caseSetup/shapes Normal file
View File

@ -0,0 +1,13 @@
/* -------------------------------*- C++ -*--------------------------------- *\
| phasicFlow File |
| copyright: www.cemf.ir |
\* ------------------------------------------------------------------------- */
objectName particleInsertion;
objectType dicrionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
names (sphere); // names of shapes
diameters (0.005); // diameter of shapes
materials (sphereMat); // material names for shapes

64
tutorials/sphereGranFlow/RotaryAirLockValve/settings/domainDict Executable file
View File

@ -0,0 +1,64 @@
/* -------------------------------*- C++ -*--------------------------------- *\
| phasicFlow File |
| copyright: www.cemf.ir |
\* ------------------------------------------------------------------------- */
objectName domainDict;
objectType dictionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
globalBox // Simulation domain: every particles that goes outside this domain will be deleted
{
min (0.397538 0.178212 0.00);
max (0.725537 0.600214 0.06);
}
decomposition
{
direction z;
}
boundaries
{
neighborListUpdateInterval 50; /* Determines how often (how many iterations) do you want to
rebuild the list of particles in the neighbor list
of all boundaries in the simulation domain */
updateInterval 10; // Determines how often do you want to update the new changes in the boundary
neighborLength 0.004; // The distance from the boundary plane within which particles are marked to be in the boundary list
left
{
type exit; // other options: periodict, reflective
}
right
{
type exit; // other options: periodict, reflective
}
bottom
{
type exit; // other options: periodict, reflective
}
top
{
type exit; // other options: periodict, reflective
}
rear
{
type exit; // other options: periodict, reflective
}
front
{
type exit; // other options: periodict, reflective
}
}

51
tutorials/sphereGranFlow/RotaryAirLockValve/settings/geometryDict
View File

@ -6,30 +6,21 @@ objectName geometryDict;
objectType dictionary; objectType dictionary;
fileFormat ASCII; fileFormat ASCII;
/*---------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------*/
motionModel rotatingAxis; // motion model: rotating object around an axis
// motion model: rotating object around an axis rotatingAxisInfo // information for rotatingAxisMotion motion model
motionModel rotatingAxisMotion;
// information for rotatingAxisMotion motion model
rotatingAxisMotionInfo
{ {
rotAxis rotAxis
{ {
p1 (0.561547 0.372714 0.000); // first point for the axis of rotation
// first point for the axis of rotation p2 (0.561547 0.372714 0.010); // second point for the axis of rotation
p1 (0.561547 0.372714 0.000);
// second point for the axis of rotation omega 2.1; // rotation speed (rad/s)
p2 (0.561547 0.372714 0.010);
// rotation speed (rad/s) startTime 1.25; // Start time of Geometry Rotating (s)
omega 2.1;
// Start time of Geometry Rotating (s) endTime 7; // End time of Geometry Rotating (s)
startTime 1.25;
// End time of Geometry Rotating (s)
endTime 7;
} }
} }
@ -37,30 +28,24 @@ surfaces
{ {
gear gear
{ {
// type of the wall type stlWall; // type of the wall
type stlWall;
// file name in stl folder file gear.stl; // file name in stl folder
file gear.stl;
// material name of this wall material wallMat; // material name of this wall
material wallMat;
// motion component name motion rotAxis; // motion component name
motion rotAxis;
} }
surfaces
surfaces
{ {
// type of the wall type stlWall; // type of the wall
type stlWall;
// file name in stl folder file surfaces.stl; // file name in stl folder
file surfaces.stl;
// material name of this wall material wallMat; // material name of this wall
material wallMat;
// motion component name motion none; // motion component name
motion none; }
} }

69
tutorials/sphereGranFlow/RotaryAirLockValve/settings/particlesDict
View File

@ -2,42 +2,73 @@
| phasicFlow File | | phasicFlow File |
| copyright: www.cemf.ir | | copyright: www.cemf.ir |
\* ------------------------------------------------------------------------- */ \* ------------------------------------------------------------------------- */
objectName geometryDict; objectName particlesDict;
objectType dictionary; objectType dictionary;
fileFormat ASCII; fileFormat ASCII;
/*---------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------*/
setFields setFields
{ {
/*
Default value for fields defined for particles:
These fields should always be defined for simulations with spherical particles
*/
defaultValue defaultValue
{ {
// linear velocity (m/s) velocity realx3 (0 0 0); // linear velocity (m/s)
velocity realx3 (0 0 0);
// linear acceleration (m/s2) acceleration realx3 (0 0 0); // linear acceleration (m/s2)
acceleration realx3 (0 0 0);
// rotational velocity (rad/s) rVelocity realx3 (0 0 0); // rotational velocity (rad/s)
rVelocity realx3 (0 0 0);
// name of the particle shape shapeName word sphere; // name of the particle shape
shapeName word sphere;
} }
selectors selectors
{} {
shapeAssigne
{
selector stridedRange; // other options: box, cylinder, sphere, randomPoints
stridedRangeInfo
{
begin 0; // begin index of points
end 50000; // end index of points
stride 3; // stride for selector
}
fieldValue // fields that the selector is applied to
{
shapeName word sphere; // sets shapeName of the selected points to largeSphere
}
}
}
} }
// positions particles positionParticles // positions particles
positionParticles
{ {
method ordered; // other options: random and empty
// creates the required fields with zero particles (empty). orderedInfo
method empty; {
diameter 0.005; // diameter of particles
// maximum number of particles in the simulation numPoints 50000; // number of particles in the simulation
maxNumberOfParticles 50000;
// perform initial sorting based on morton code? axisOrder (z x y); // axis order for filling the space with particles
mortonSorting Yes; }
regionType cylinder; // other options: box and sphere
cylinderInfo // cylinder information for positioning particles
{
p1 (0.0 0.0 0.003); // begin point of cylinder axis
p2 (0.0 0.0 0.22); // end point of cylinder axis
radius 0.117; // radius of cylinder
}
} }

41
tutorials/sphereGranFlow/RotaryAirLockValve/settings/settingsDict
View File

@ -6,43 +6,34 @@ objectName geometryDict;
objectType dictionary; objectType dictionary;
fileFormat ASCII; fileFormat ASCII;
/*---------------------------------------------------------------------------*/ /*---------------------------------------------------------------------------*/
run rotatingValve; run rotatingValve;
// time step for integration (s) dt 0.00001; // time step for integration (s)
dt 0.00001;
// start time for simulation startTime 0; // start time for simulation
startTime 0;
// end time for simulation endTime 7; // end time for simulation
endTime 7;
// time interval for saving the simulation saveInterval 0.05; // time interval for saving the simulation
saveInterval 0.05;
// maximum number of digits for time folder timePrecision 6; // maximum number of digits for time folder
timePrecision 6;
// gravity vector (m/s2) g (0 -9.8 0); // gravity vector (m/s2)
g (0 -9.8 0);
/* /*
Simulation domain every particles that goes outside this domain is deleted. Simulation domain every particles that goes outside this domain is deleted
*/ */
domain includeObjects (diameter); // save necessary (i.e., required) data on disk
{
min (0.397538 0.178212 0.00);
max (0.725537 0.600214 0.06);
}
// integration method // exclude unnecessary data from saving on disk
integrationMethod AdamsBashforth3; excludeObjects (rVelocity.dy1 pStructPosition.dy1 pStructVelocity.dy1);
// report timers? integrationMethod AdamsBashforth2; // integration method
timersReport Yes;
// time interval for reporting timers writeFormat ascii; // data writting format (ascii or binary)
timersReportInterval 0.01;
timersReport Yes; // report timers: Yes or No
timersReportInterval 0.01; // time interval for reporting timers