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

V-blender updated.

This commit is contained in:
Ramin Khodabandeh 2024-05-28 18:03:19 +04:30
parent 8490668d51
commit 8d14850a1e
7 changed files with 202 additions and 198 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

73
tutorials/sphereGranFlow/V-blender/caseSetup/interaction Normal file → Executable file
View File

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

110
tutorials/sphereGranFlow/V-blender/caseSetup/particleInsertion Normal file → Executable file
View File

@ -6,12 +6,9 @@ objectName particleInsertion;
objectType dicrionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
active Yes; // is insertion active -> yes or no
// is insertion active?
active yes;
// not implemented for yes
collisionCheck No;
checkForCollision No; // is checked -> yes or no
/*
two layers of particles are packed one-by-one using tho insertion steps.
@ -19,84 +16,79 @@ 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;
regionType cylinder; // type of insertion region
cylinderRegionInfo
cylinderInfo
{
// 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);
// 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); // (m,m,m)
p2 (0.1150615 0.12 0.4811585); // (m,m,m)
radius 0.1; // radius of cylinder (m)
}
timeControl simulationTime;
insertionInterval 0.025; // seconds
rate 10000; // insertion rate (particles/s)
startTime 0.0; // start time of small particles insertion (s)
endTime 1; // end time of small particles insertion (s)
interval 0.025; // time interval of small particles insertion (s)
setFields
{
// initial velocity of inserted particles
velocity realx3 (1.2 0.0 -1.2);
velocity realx3 (1.2 0.0 -1.2); // initial velocity of inserted particles
}
mixture
{
// first layer of inserted particles
smallSphere 1;
smallSphere 1; // first layer of inserted particles
}
}
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;
regionType cylinder; // type of insertion region
cylinderRegionInfo
cylinderInfo
{
// 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);
// coordinates of center of both ends of the insertion cylinder on
// the right side of the V-blender (m,m,m)
p1 ( 0.7562545 0.12 0.50079); // (m,m,m)
p2 ( 0.7362545 0.12 0.48079); // (m,m,m)
radius 0.1; // radius of cylinder (m)
}
timeControl simulationTime;
insertionInterval 0.025; // seconds
rate 10000; // insertion rate (particles/s)
startTime 1.5; // start time of large particles insertion (s)
endTime 2.5; // end time of large particles insertion (s)
interval 0.025; // time interval of large particles insertion (s)
setFields
{
// initial velocity of inserted particles
velocity realx3 (-1.2 0.0 -1.2);
velocity realx3 (-1.2 0.0 -1.2); // initial velocity of inserted particles
}
mixture
{
// second layer of inserted particles
largeSphere 1;
largeSphere 1; // second layer of inserted particles
}
}

11
tutorials/sphereGranFlow/V-blender/caseSetup/sphereShape → tutorials/sphereGranFlow/V-blender/caseSetup/shapes Normal file → Executable file
View File

@ -6,12 +6,11 @@ objectName sphereDict;
objectType sphereShape;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
names (smallSphere largeSphere); // names of particles
diameters (0.01 0.0101); // diameter of particles
materials (lightMat lightMat); // material of the particles
// names of particles
names (smallSphere largeSphere);
// diameter of particles
diameters (0.01 0.0101);
// material names for particles
materials (lightMat lightMat);

59
tutorials/sphereGranFlow/V-blender/settings/domainDict Executable file
View File

@ -0,0 +1,59 @@
/* -------------------------------*- 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.1 -0.4 0);
max (0.86 0.6 0.6);
}
boundaries
{
// 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
neighborListUpdateInterval 50;
// Determines how often do you want to update the new changes in the boundary
updateInterval 0.0001;
// The distance from the boundary plane within which particles are marked to be in the boundary list
neighborLength 0.004;
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
}
}

63
tutorials/sphereGranFlow/V-blender/settings/geometryDict Normal file → Executable file
View File

@ -6,49 +6,38 @@ objectName geometryDict;
objectType dictionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
motionModel rotatingAxis; // motion model can be rotatingAxis or stationary or vibrating
// 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
rotatingAxisInfo // information for rotatingAxisMotion motion model
{
rotAxis
{
// first point for the axis of rotation
p1 (0.128228 0.116446 0.297901);
p1 (0.128228 0.116446 0.297901); // first point for the axis of rotation
p2 (0.722596 0.116459 0.297901); // second point for the axis of rotation
// second point for the axis of rotation
p2 (0.722596 0.116459 0.297901);
omega 3.14; // rotation speed (rad/s)
// rotation speed (rad/s)
omega 3.14;
startTime 3; // start time of rotation
// start time of rotation
startTime 3;
// end time of rotation
endTime 10;
endTime 10; // end time of rotation
}
}
surfaces
{
body
{
type stlWall; // type of the wall
file body.stl; // file name in stl folder
material wallMat; // material name of this wall
motion rotAxis; // motion component name
}
}

39
tutorials/sphereGranFlow/V-blender/settings/particlesDict Normal file → Executable file
View File

@ -6,7 +6,6 @@ objectName particlesDict;
objectType dictionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
setFields
{
/*
@ -17,32 +16,36 @@ setFields
defaultValue
{
// linear velocity (m/s)
velocity realx3 (0 0 0);
// linear acceleration (m/s2)
acceleration realx3 (0 0 0);
velocity realx3 (0 0 0); // linear velocity (m/s)
acceleration realx3 (0 0 0); // linear acceleration (m/s2)
rVelocity realx3 (0 0 0); // rotational velocity (rad/s)
// rotational velocity (rad/s)
rVelocity realx3 (0 0 0);
// name of the particle shape
shapeName word smallSphere;
shapeName word smallSphere; // name of the particle shape
}
selectors
{}
}
// positions particles
positionParticles
positionParticles // positions particles
{
// no particle is initially in the simulation
method empty;
method empty; // other options: ordered and random
// maximum number of particles in the simulation
maxNumberOfParticles 25000;
maxNumberOfParticles 25000; // maximum number of particles in the simulation
// perform initial sorting based on morton code?
mortonSorting Yes;
regionType box; // other options: cylinder and sphere
boxInfo // box for positioning particles
{
min (-0.08 -0.08 0.015); // lower corner point of the box
max ( 0.08 0.08 0.098); // upper corner point of the box
}
}

45
tutorials/sphereGranFlow/V-blender/settings/settingsDict Normal file → Executable file
View File

@ -6,45 +6,32 @@ objectName settingsDict;
objectType dictionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
run rotatingVblender;
run rotatingVblender;
dt 0.00001; // time step for integration
// time step for integration (s)
dt 0.00001;
startTime 0.0; // start time for simulation
// start time for simulation
startTime 0;
endTime 10.0; // end time for simulation
// end time for simulation
endTime 10;
saveInterval 0.05; // time interval for saving the simulation
// time interval for saving the simulation
saveInterval 0.05;
timePrecision 6; // maximum number of digits for time folder
// maximum number of digits for time folder
timePrecision 6;
g (0 0 -9.8); // gravity vector (m/s2)
// gravity vector (m/s2)
g (0 0 -9.8);
// save necessary (i.e., required) data on disk
includeObjects (diameter);
// exclude unnecessary data from saving on disk
excludeObjects (rVelocity.dy1 pStructPosition.dy1 pStructVelocity.dy1);
/*
Simulation domain
every particles that goes outside this domain is deleted.
*/
integrationMethod AdamsBashforth2; // integration method
domain
{
min (-0.1 -0.4 0);
max (0.86 0.6 0.6);
}
writeFormat ascii; // data writting format (ascii or binary)
// integration method
integrationMethod AdamsBashforth2;
timersReport Yes; // report timers
timersReportInterval 0.01; // time interval for reporting timers
writeFormat ascii;
// report timers?
timersReport Yes;
// time interval for reporting timers
timersReportInterval 0.01;