HourGlass3DDemo.cpp File Reference

#include "Mercury3D.h"
#include "Walls/IntersectionOfWalls.h"
#include "Walls/AxisymmetricIntersectionOfWalls.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"

Classes
class	HourGlass

Functions
int	main (int argc, char *argv[])

Function Documentation

main()

int main	(	int argc	,
		char *	argv[]
	)

{
    logger(INFO, "Hourglass Simulation");
    // note: the scaling of the hourglass is based on Stefan Luding's demo code MDCLR/DEMO/W7
 
    //all parameters should be set in the main file
    //here, we use SI units (generally, other, scaled units are possible)
 
    //create an instance of the class and name it
    HourGlass HG;
    HG.setName("HourGlass3DDemo");
 
 
    //these parameters are needed in setupInitialConditions()
 
    //specify geometry
    //specify dimensions of the hourglass
    const Mdouble width = 10e-2; // 10cm
    const Mdouble height = 60e-2; // 60cm
    //set domain accordingly (domain boundaries are not walls!)
    HG.setXMin(0.0);
    HG.setXMax(width);
    HG.setYMin(0.0);
    HG.setYMax(width);
    HG.setZMin(0.0);
    HG.setZMax(height);
 
    //specify how big the wedge of the contraction should be
    const Mdouble contractionWidth = 2.5e-2; //2.5cm
    const Mdouble contractionHeight = 5e-2; //5cm
    HG.contractionWidth = contractionWidth;
    HG.contractionHeight = contractionHeight;
 
    //specify particle properties
 
    //these parameters are needed in setupInitialConditions()
    HG.minParticleRadius = 6e-3; // 6mm
    HG.maxParticleRadius = 10e-3; //10mm
 
    //specify body forces
    HG.setGravity(Vec3D(0.0, 0.0, -9.8));
 
    //make the species of the particle and wall
    LinearViscoelasticFrictionSpecies species;
    species.setDensity(2000);
    //specify contact properties
    //normal forces
    species.setStiffness(1e5);
    species.setDissipation(0.63);
    //tangential (sliding) forces
    species.setSlidingFrictionCoefficient(0.5);
    species.setSlidingStiffness(1.2e4);
    species.setSlidingDissipation(0.16);
    //tangential (rolling) torques
    species.setRollingFrictionCoefficient(0.2);
    species.setRollingStiffness(1.2e4);
    species.setRollingDissipation(6.3e-2);
    //normal (torsion/spin) torques
    species.setTorsionFrictionCoefficient(0.1);
    species.setTorsionStiffness(1.2e4);
    species.setSlidingDissipation(6.3e-2);
    HG.speciesHandler.copyAndAddObject(species);
 
    //test normal forces
    const Mdouble minParticleMass = species.getDensity() * 4.0 / 3.0 * constants::pi * mathsFunc::cubic(HG.minParticleRadius);
    //Calculates collision time for two copies of a particle of given dissipation_, k, effective mass
    logger(INFO, "minParticleMass = %", minParticleMass);
    //Calculates collision time for two copies of a particle of given dissipation_, k, effective mass
    const Mdouble tc = species.getCollisionTime(minParticleMass);
    logger(INFO, "tc  = %", tc);
    //Calculates restitution coefficient for two copies of given dissipation_, k, effective mass
    const Mdouble r = species.getRestitutionCoefficient(minParticleMass);
    logger(INFO, "restitution coefficient = %", r);
 
    //set other simulation parameters
    HG.setTimeStep(tc / 50.0);
    HG.setTimeMax(10.0);
    HG.setSaveCount(500);
    HG.setXBallsAdditionalArguments("-v0 -solidf");
    HG.N = 600; //number of particles
    //uncomment next two line 2 to create paraview files
    HG.setParticlesWriteVTK(true);
    HG.wallHandler.setWriteVTK(FileType::ONE_FILE);
 
    HG.solve(argc, argv);
    return 0;
}

References HourGlass::contractionHeight, HourGlass::contractionWidth, BaseHandler< T >::copyAndAddObject(), mathsFunc::cubic(), ParticleSpecies::getDensity(), INFO, logger, HourGlass::maxParticleRadius, HourGlass::minParticleRadius, HourGlass::N, ONE_FILE, constants::pi, ParticleSpecies::setDensity(), DPMBase::setGravity(), DPMBase::setName(), DPMBase::setParticlesWriteVTK(), DPMBase::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), WallHandler::setWriteVTK(), DPMBase::setXBallsAdditionalArguments(), DPMBase::setXMax(), DPMBase::setXMin(), DPMBase::setYMax(), DPMBase::setYMin(), DPMBase::setZMax(), DPMBase::setZMin(), DPMBase::solve(), DPMBase::speciesHandler, and DPMBase::wallHandler.