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

homogenizationSilo-PeriodicBoundary Updated.

This commit is contained in:
ramin1728 2025-02-20 20:05:05 +03:30
parent b8bfd04dfa
commit 78b32b5027
8 changed files with 151 additions and 142 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
tutorials/sphereGranFlow/conveyorBelt/settings/particlesDict
View File

@ -18,13 +18,13 @@ setFields
defaultValue
{
velocity realx3 (0 0 0); // linear velocity (m/s)
velocity realx3 (0 0 0); // linear velocity (m/s)
acceleration realx3 (0 0 0); // linear acceleration (m/s2)
acceleration realx3 (0 0 0); // linear acceleration (m/s2)
rVelocity realx3 (0 0 0); // rotational velocity (rad/s)
rVelocity realx3 (0 0 0); // rotational velocity (rad/s)
shapeName word lightSphere; // name of the particle shape
shapeName word lightSphere; // name of the particle shape
}
selectors
@ -33,9 +33,9 @@ setFields
positionParticles // positions particles
{
method empty; // other options: ordered and random
method empty; // other options: ordered and random
regionType box; // other options: cylinder and sphere
regionType box; // other options: cylinder and sphere
boxInfo // box region for positioning particles
{

20
tutorials/sphereGranFlow/homogenizationSilo-PeriodicBoundary/caseSetup/interaction
View File

@ -6,30 +6,30 @@ objectName interaction;
objectType dicrionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
materials (lightMat heavyMat wallMat); // a list of materials names
materials (lightMat heavyMat wallMat); // a list of materials names
densities (1000 1500.0 2500); // density of materials [kg/m3]
densities (1000 1500.0 2500); // density of materials [kg/m3]
contactListType sortedContactList;
contactListType sortedContactList;
contactSearch
{
method NBS; // method for broad search particle-particle
method NBS; // method for broad search particle-particle
updateInterval 10;
updateInterval 10;
sizeRatio 1.1;
sizeRatio 1.1;
cellExtent 0.55;
cellExtent 0.55;
adjustableBox No;
adjustableBox No;
}
model
{
contactForceModel nonLinearLimited;
contactForceModel nonLinearLimited;
rollingFrictionModel normal;
rollingFrictionModel normal;
/*

149
tutorials/sphereGranFlow/homogenizationSilo-PeriodicBoundary/caseSetup/particleInsertion
View File

@ -6,33 +6,33 @@ objectName particleInsertion;
objectType dicrionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
active Yes; // is insertion active -> yes or no
active Yes; // is insertion active -> yes or no
/*
one layers of particles are packed
*/
layer0
{
timeControl simulationTime;
regionType cylinder; // type of insertion region
regionType cylinder; // type of insertion region
rate 5100; // insertion rate (particles/s)
rate 5100; // insertion rate (particles/s)
startTime 0; // (s)
startTime 0; // (s)
endTime 0.5; // (s)
endTime 0.5; // (s)
insertionInterval 0.025; // s
insertionInterval 0.025; // (s)
cylinderInfo
{
radius 0.09; // radius of cylinder (m)
radius 0.09; // radius of cylinder (m)
p1 (0.0 0.0 0.1); // (m,m,m)
p2 (0.0 0.0 0.11); // (m,m,m)
p1 (0.0 0.0 0.1); // (m,m,m)
p2 (0.0 0.0 0.11); // (m,m,m)
}
setFields
@ -42,173 +42,180 @@ layer0
mixture
{
parType1 1; // mixture composition of inserted particles
parType1 1; // mixture composition of inserted particles
}
}
layer1
{
timeControl simulationTime;
timeControl simulationTime;
regionType cylinder; // type of insertion region
regionType cylinder; // type of insertion region
rate 5100; // insertion rate (particles/s)
rate 5100; // insertion rate (particles/s)
startTime 0.7; // (s)
startTime 0.7; // (s)
endTime 1.2; // (s)
endTime 1.2; // (s)
insertionInterval 0.025; // s
insertionInterval 0.025; // s
cylinderInfo
{
radius 0.09;
p1 (0.0 0.0 0.16 ); // (m,m,m)
p2 (0.0 0.0 0.17); // (m,m,m)
radius 0.09;
p1 (0.0 0.0 0.16); // (m,m,m)
p2 (0.0 0.0 0.17); // (m,m,m)
}
setFields
{
velocity realx3 (0.0 0.0 -0.6);
velocity realx3 (0.0 0.0 -0.6);
}
mixture
{
parType2 1; // only parType2
parType2 1; // only parType2
}
}
layer2
{
timeControl simulationTime;
timeControl simulationTime;
regionType cylinder; // type of insertion region
regionType cylinder; // type of insertion region
rate 5100; // insertion rate (particles/s)
rate 5100; // insertion rate (particles/s)
startTime 1.4; // (s)
startTime 1.4; // (s)
endTime 1.9; // (s)
endTime 1.9; // (s)
insertionInterval 0.025; // s
insertionInterval 0.025; // (s)
cylinderInfo
{
radius 0.09;
p1 ( 0.0 0.0 0.2 ); // (m,m,m)
p2 ( 0.0 0.0 0.21); // (m,m,m)
radius 0.09;
p1 (0.0 0.0 0.2); // (m,m,m)
p2 (0.0 0.0 0.21); // (m,m,m)
}
setFields
{
velocity realx3 (0.0 0.0 -0.6);
velocity realx3 (0.0 0.0 -0.6);
}
mixture
{
lightSphere 1; // only lightSphere
lightSphere 1; // only lightSphere
}
}
layer3
{
timeControl simulationTime;
timeControl simulationTime;
regionType cylinder; // type of insertion region
regionType cylinder; // type of insertion region
rate 5100; // insertion rate (particles/s)
rate 5100; // insertion rate (particles/s)
startTime 2.1; // (s)
startTime 2.1; // (s)
endTime 2.6; // (s)
endTime 2.6; // (s)
insertionInterval 0.025; // s
insertionInterval 0.025; // (s)
cylinderInfo
{
radius 0.09;
p1 ( 0.0 0.0 0.28 ); // (m,m,m)
p2 ( 0.0 0.0 0.29); // (m,m,m)
radius 0.09;
p1 (0.0 0.0 0.28); // (m,m,m)
p2 (0.0 0.0 0.29); // (m,m,m)
}
setFields
{
velocity realx3 (0.0 0.0 -0.6);
velocity realx3 (0.0 0.0 -0.6);
}
mixture
{
heavySphere 1;
heavySphere 1;
}
}
layer4
{
timeControl simulationTime;
timeControl simulationTime;
regionType cylinder; // type of insertion region
regionType cylinder; // type of insertion region
rate 5100; // insertion rate (particles/s)
rate 5100; // insertion rate (particles/s)
startTime 2.8; // (s)
startTime 2.8; // (s)
endTime 3.3; // (s)
endTime 3.3; // (s)
insertionInterval 0.025; // s
insertionInterval 0.025; // (s)
cylinderInfo
{
radius 0.09;
p1 ( 0.0 0.0 0.37 ); // (m,m,m)
p2 ( 0.0 0.0 0.38); // (m,m,m)
radius 0.09;
p1 (0.0 0.0 0.37); // (m,m,m)
p2 (0.0 0.0 0.38); // (m,m,m)
}
setFields
{
velocity realx3 (0.0 0.0 -0.6);
velocity realx3 (0.0 0.0 -0.6);
}
mixture
{
lightSphere 1;
lightSphere 1;
}
}
layer5
{
timeControl simulationTime;
timeControl simulationTime;
regionType cylinder; // type of insertion region
regionType cylinder; // type of insertion region
rate 5100; // insertion rate (particles/s)
rate 5100; // insertion rate (particles/s)
startTime 3.4; // (s)
startTime 3.4; // (s)
endTime 3.9; // (s)
endTime 3.9; // (s)
insertionInterval 0.025; // s
insertionInterval 0.025; // (s)
cylinderInfo
{
radius 0.09;
p1 ( 0.0 0.0 0.38 ); // (m,m,m)
p2 ( 0.0 0.0 0.39); // (m,m,m)
radius 0.09;
p1 (0.0 0.0 0.38); // (m,m,m)
p2 (0.0 0.0 0.39); // (m,m,m)
}
setFields
{
velocity realx3 (0.0 0.0 -0.6);
velocity realx3 (0.0 0.0 -0.6);
}
mixture
{
heavySphere 1;
heavySphere 1;
}
}

6
tutorials/sphereGranFlow/homogenizationSilo-PeriodicBoundary/caseSetup/shapes
View File

@ -6,10 +6,10 @@ objectName sphereDict;
objectType sphereShape;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
names (parType1 parType2); // names of shapes
names (parType1 parType2); // names of shapes
diameters (0.00885 0.0089); // diameter of shapes
diameters (0.00885 0.0089); // diameter of shapes
materials (lightMat heavyMat); // material names for shapes
materials (lightMat heavyMat); // material names for shapes

4
tutorials/sphereGranFlow/homogenizationSilo-PeriodicBoundary/settings/domainDict
View File

@ -6,9 +6,7 @@ objectName domainDict;
objectType dictionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
// Simulation domain: every particles that goes outside this domain will be deleted
globalBox
globalBox // Simulation domain: every particles that goes outside this domain will be deleted
{
min (-0.11 -0.11 -0.15);

62
tutorials/sphereGranFlow/homogenizationSilo-PeriodicBoundary/settings/geometryDict
View File

@ -6,19 +6,21 @@ objectName geometryDict;
objectType dictionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
motionModel rotatingAxis;
motionModel rotatingAxis; // motion model can be rotatingAxis or stationary or vibrating
rotatingAxisInfo
{
// for opening the gate of silo between time 4.1 and 5.1 s
gateMotion
gateMotion // for opening the gate of silo between time 4.1 and 5.1 s
{
p1 (-0.04 -0.04 -0.1);
p2 (-0.04 -0.04 0.0);
omega 3.14;
startTime 4.1;
endTime 5.1;
p1 (-0.04 -0.04 -0.1);
p2 (-0.04 -0.04 0.0);
omega 3.14;
startTime 4.1;
endTime 5.1;
}
}
@ -26,36 +28,36 @@ surfaces
{
cylinderShell
{
type cylinderWall; // other options: cuboidWall and planeWall
type cylinderWall; // other options: cuboidWall and planeWall
p1 (0.0 0.0 0.0); // begin point of cylinder axis
p1 (0.0 0.0 0.0); // begin point of cylinder axis
p2 (0.0 0.0 0.4); // end point of cylinder axis
p2 (0.0 0.0 0.4); // end point of cylinder axis
radius1 0.1; // radius at p1
radius1 0.1; // radius at p1
radius2 0.1; // radius at p2
radius2 0.1; // radius at p2
resolution 36; // number of divisions
resolution 36; // number of divisions
material wallMat; // material name of this wall
material wallMat; // material name of this wall
}
coneShell
{
type cylinderWall; // other options: cuboidWall and planeWall
type cylinderWall; // other options: cuboidWall and planeWall
p1 (0.0 0.0 -0.1); // begin point of cylinder axis
p1 (0.0 0.0 -0.1); // begin point of cylinder axis
p2 (0.0 0.0 0.0); // end point of cylinder axis
p2 (0.0 0.0 0.0); // end point of cylinder axis
radius1 0.04; // radius at p1
radius1 0.04; // radius at p1
radius2 0.1; // radius at p2
radius2 0.1; // radius at p2
resolution 36; // number of divisions
resolution 36; // number of divisions
material wallMat; // material name of this wall
material wallMat; // material name of this wall
}
/*
@ -64,19 +66,19 @@ surfaces
exitGate
{
type planeWall; // other options: cuboidWall and cylinderWall
type planeWall; // other options: cuboidWall and cylinderWall
p1 (-0.04 -0.04 -0.1); // first point of the wall
p1 (-0.04 -0.04 -0.1); // first point of the wall
p2 ( 0.04 -0.04 -0.1); // second point of the wall
p2 (0.04 -0.04 -0.1); // second point of the wall
p3 ( 0.04 0.04 -0.1); // third point of the wall
p3 (0.04 0.04 -0.1); // third point of the wall
p4 (-0.04 0.04 -0.1); // fourth point of the wall
p4 (-0.04 0.04 -0.1); // fourth point of the wall
material wallMat; // material name of the wall
material wallMat; // material name of the wall
motion gateMotion;
motion gateMotion;
}
}

8
tutorials/sphereGranFlow/homogenizationSilo-PeriodicBoundary/settings/particlesDict
View File

@ -16,13 +16,13 @@ setFields
defaultValue
{
velocity realx3 (0 0 0); // linear velocity (m/s)
velocity realx3 (0 0 0); // linear velocity (m/s)
acceleration realx3 (0 0 0); // linear acceleration (m/s2)
acceleration realx3 (0 0 0); // linear acceleration (m/s2)
rVelocity realx3 (0 0 0); // rotational velocity (rad/s)
rVelocity realx3 (0 0 0); // rotational velocity (rad/s)
shapeName word parType1; // name of the particle shape
shapeName word parType1; // name of the particle shape
}
selectors

32
tutorials/sphereGranFlow/homogenizationSilo-PeriodicBoundary/settings/settingsDict
View File

@ -6,36 +6,38 @@ objectName settingsDict;
objectType dictionary;
fileFormat ASCII;
/*---------------------------------------------------------------------------*/
run homogenizationSilo;
run homogenizationSilo;
dt 0.00001; // time step for integration (s)
dt 0.00001; // time step for integration (s)
startTime 0.0; // start time for simulation
startTime 0.0; // start time for simulation
endTime 20; // end time for simulation
endTime 20; // end time for simulation
saveInterval 0.05; // time interval for saving the simulation
saveInterval 0.05; // time interval for saving the simulation
timePrecision 4; // maximum number of digits for time folder
timePrecision 4; // maximum number of digits for time folder
g (0 0 -9.8); // gravity vector (m/s2)
g (0 0 -9.8); // gravity vector (m/s2)
// overrides the default behavior
includeObjects (diameter);
includeObjects (diameter);
// exclude unnecessary data from saving on disk
excludeObjects (rVelocity.dy1 rVelocity.dy2 rVelocity.dy3
pStructPosition.dy1 pStructPosition.dy2 pStructPosition.dy3
pStructVelocity.dy1 pStructVelocity.dy2 pStructVelocity.dy3);
integrationMethod AdamsBashforth4; // integration method
excludeObjects (rVelocity.dy1 rVelocity.dy2 rVelocity.dy3
pStructPosition.dy1 pStructPosition.dy2 pStructPosition.dy3
pStructVelocity.dy1 pStructVelocity.dy2 pStructVelocity.dy3);
writeFormat binary; // data writting format (ascii or binary)
integrationMethod AdamsBashforth4; // integration method
timersReport Yes; // report timers
writeFormat binary; // data writting format (ascii or binary)
timersReportInterval 0.1; // time interval for reporting timers
timersReport Yes; // report timers
timersReportInterval 0.1; // time interval for reporting timers