diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md index d01d31e2..8d5c799c 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/ReadMe.md @@ -15,7 +15,7 @@ The problem is to simulate a rotating drum with a diameter of 0.24 m, a length o # 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 sorted in three folders: `caseSetup`, `setting` and `stl` (see the above folders). +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 sorted in three folders: `caseSetup`, `setting` and `stl`. ## Defining small and large particles Then in the `caseSetup/shapes` the diameter and the material name of the particles are defined. Two sizes are defined: 4 and 5 mm. @@ -43,7 +43,8 @@ in caseSetup/particleInsertion file layerrightregion { // type of insertion region - type cylinderRegion; + timeControl simulationTime; + regionType cylinder; // insertion rate (particles/s) rate 12500; // Start time of LightParticles insertion (s) @@ -53,7 +54,7 @@ layerrightregion // Time Interval of LightParticles insertion (s) interval 0.025; - cylinderRegionInfo + cylinderInfo { // Coordinates of cylinderRegion (m,m,m) p2 (-0.15 0.25 0.05); @@ -93,10 +94,6 @@ densities (1000 1500 2500); 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 @@ -166,6 +163,7 @@ surfaces 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++ +motionModel rotatingAxis; rotatingAxisInfo { rotAxis @@ -186,7 +184,7 @@ rotatingAxisInfo ## Performing Simulation To run simulations, type the following commands in the terminal one at a time. -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 stored in ./VTK folder. \ No newline at end of file +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 stored in ./VTK folder. \ No newline at end of file diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/interaction b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/interaction index 27d71289..7572adb9 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/interaction +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/interaction @@ -14,54 +14,50 @@ contactListType sortedContactList; contactSearch { - method NBS; // method for broad search + method NBS; // method for broad search - updateInterval 10; + updateInterval 10; - sizeRatio 1.1; + sizeRatio 1.1; - cellExtent 0.55; + cellExtent 0.55; - adjustableBox Yes; + adjustableBox Yes; } model { - contactForceModel nonLinearLimited; + contactForceModel nonLinearLimited; - rollingFrictionModel normal; + rollingFrictionModel normal; - /* - Property (lightMat-lightMat lightMat-heavyMat lightMat-wallMat - heavyMat-heavyMat heavyMat-wallMat - wallMat-wallMat); - */ + /* + Property (lightMat-lightMat lightMat-heavyMat lightMat-wallMat + heavyMat-heavyMat heavyMat-wallMat + wallMat-wallMat); + */ - Yeff (1.0e6 1.0e6 1.0e6 - 1.0e6 1.0e6 - 1.0e6); // Young modulus [Pa] + Yeff (1.0e6 1.0e6 1.0e6 + 1.0e6 1.0e6 + 1.0e6); // Young modulus [Pa] - Geff (0.8e6 0.8e6 0.8e6 - 0.8e6 0.8e6 - 0.8e6); // Shear modulus [Pa] + Geff (0.8e6 0.8e6 0.8e6 + 0.8e6 0.8e6 + 0.8e6); // Shear modulus [Pa] - nu (0.25 0.25 0.25 - 0.25 0.25 - 0.25); // Poisson's ratio [-] + nu (0.25 0.25 0.25 + 0.25 0.25 + 0.25); // Poisson's ratio [-] - en (0.97 0.97 0.85 - 0.97 0.85 - 1.00); // coefficient of normal restitution - - et (1.0 1.0 1.0 - 1.0 1.0 - 1.0); // coefficient of tangential restitution + en (0.97 0.97 0.85 + 0.97 0.85 + 1.00); // coefficient of normal restitution - mu (0.65 0.65 0.35 - 0.65 0.35 - 0.35); // dynamic friction + mu (0.65 0.65 0.35 + 0.65 0.35 + 0.35); // dynamic friction - mur (0.1 0.1 0.1 - 0.1 0.1 - 0.1); // rolling friction + mur (0.1 0.1 0.1 + 0.1 0.1 + 0.1); // rolling friction } diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion index 73ebc90d..6ba98b98 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/caseSetup/particleInsertion @@ -15,70 +15,69 @@ active Yes; // is insertion active -> Yes o layerrightregion // Right Layer Region { - timeControl simulationTime; + timeControl simulationTime; - regionType cylinder; // type of insertion region + regionType cylinder; // type of insertion region - rate 12500; // Right Region Particles Insertion Rate (particles/s) + rate 12500; // Right Region Particles Insertion Rate (particles/s) - startTime 0; // Start time of LightParticles insertion (s) + startTime 0; // Start time of LightParticles insertion (s) - endTime 1; // End time of LightParticles insertion (s) + endTime 1; // End time of LightParticles insertion (s) - insertionInterval 0.025; // Time Interval of LightParticles insertion (s) + insertionInterval 0.025; // Time Interval of LightParticles insertion (s) - cylinderInfo - { + cylinderInfo + { + p2 (-0.15 0.25 0.05); // Top of cylinderRegion (m,m,m) - p2 (-0.15 0.25 0.05); // Top of cylinderRegion (m,m,m) + p1 (-0.15 0.24 0.05); // Bottom of cylinderRegion (m,m,m) - p1 (-0.15 0.24 0.05); // Bottom of cylinderRegion (m,m,m) - - radius 0.035; // radius of cylinder (m) - } + radius 0.035; // radius of cylinder (m) + } - setFields - { - velocity realx3 (0.0 -0.6 0.0); // initial velocity of inserted particles - } + setFields + { + velocity realx3 (0.0 -0.6 0.0); // initial velocity of inserted particles + } - mixture - { - smallSphere 1; // mixture composition of inserted particles - } + mixture + { + smallSphere 1; // mixture composition of inserted particles + } } layerleftregion // Left Layer Region { - timeControl simulationTime; + timeControl simulationTime; - regionType cylinder; // type of insertion region + regionType cylinder; // type of insertion region - rate 7500; // Left Region Particles Insertion Rate (particles/s) + rate 7500; // Left Region Particles Insertion Rate (particles/s) - startTime 0; // Start time of LightParticles insertion (s) + startTime 0; // Start time of LightParticles insertion (s) - endTime 1; // End time of LightParticles insertion (s) + endTime 1; // End time of LightParticles insertion (s) - insertionInterval 0.025; // Time Interval of LightParticles insertion (s) + insertionInterval 0.025; // Time Interval of LightParticles insertion (s) - cylinderInfo - { - p2 (-0.23 0.25 0.05); // Top of cylinderRegion (m,m,m) + cylinderInfo + { + p2 (-0.23 0.25 0.05); // Top of cylinderRegion (m,m,m) - p1 (-0.23 0.24 0.05); // Bottom of cylinderRegion (m,m,m) + p1 (-0.23 0.24 0.05); // Bottom of cylinderRegion (m,m,m) - radius 0.035; // radius of cylinder (m) - } + radius 0.035; // radius of cylinder (m) + } - setFields - { - velocity realx3 (0.0 -0.6 0.0); // initial velocity of inserted particles - } + setFields + { + velocity realx3 (0.0 -0.6 0.0); // initial velocity of inserted particles + } - mixture - { - largeSphere 1; // only heavySphere - } + mixture + { + largeSphere 1; // only heavySphere + } } diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/domainDict b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/domainDict index cd36ab02..74c86ad3 100755 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/domainDict +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/domainDict @@ -8,16 +8,16 @@ fileFormat ASCII; /*---------------------------------------------------------------------------*/ globalBox // Simulation domain: every particles that goes outside this domain will be deleted { - min (-0.328 0.095 -0.025); // lower corner point of the box + min (-0.328 0.095 -0.025); // lower corner point of the box - max (-0.068 0.355 0.125); // upper corner point of the box + max (-0.068 0.355 0.125); // upper corner point of the box } boundaries { - left + left { - type exit; // other options: periodic, reflective + type exit; // other options: periodic, reflective } right diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict index 99132c9e..a7502d1c 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/geometryDict @@ -10,7 +10,7 @@ motionModel rotatingAxis; rotatingAxisInfo { - rotAxis + rotAxis { p1 (-0.1974 0.2269 0); // first point for the axis of rotation diff --git a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict index 6a182c2b..b4d5bf80 100644 --- a/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict +++ b/tutorials/sphereGranFlow/RotatingDrumWithBaffles/settings/particlesDict @@ -16,7 +16,7 @@ setFields spherical particles. */ - defaultValue + defaultValue { velocity realx3 (0 0 0); // linear velocity (m/s)