MonodisperseSmoothInclinedChute.cpp File Reference

#include "Chute.h"
#include <Species/LinearViscoelasticFrictionSpecies.h>
#include <Walls/InfiniteWall.h>

Classes
class	SmoothChute

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

Function Documentation

main()

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

{
    logger(INFO, "Chute flow over a smooth bottom");
 
    //set the problem parameters: name, simulation time and time step
    SmoothChute problem;
    problem.setName("MonodisperseSmoothInclinedChute");
    problem.setTimeMax(100);
    problem.setTimeStep(1e-4);
 
    //set the direction of the gravity as 30 degrees, and magnitude 9.81 m/s^2
    problem.setChuteAngleAndMagnitudeOfGravity(30.0, 9.81);
    
    problem.numberOfParticlesToBeInserted = 286;
    problem.insertTime = 0.006;
    problem.insertTimeInterval = problem.insertTime;
    //
    auto particleSpecies = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
    auto wallSpecies = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
    auto wallParticleSpecies = problem.speciesHandler.getMixedObject(wallSpecies, particleSpecies);
    
    // Mdouble effectiveMass = 0.5*particleSpecies->getMassFromRadius(0.5*0.003);
    // std::cout << "Mass" << 2.0*effectiveMass << std::endl;
    // helpers::KAndDisp kAndDispPP = helpers::computeKAndDispFromCollisionTimeAndRestitutionCoefficientAndEffectiveMass(0.005,0.97,effectiveMass);
    //helpers::KAndDisp kAndDispPW = helpers::computeKAndDispFromCollisionTimeAndRestitutionCoefficientAndEffectiveMass(0.005,0.95,effectiveMass);
    
    //std::cout << "Stifness" << kAndDispPP.k << std::endl;
 
    //set parameters for the particle species: density, stiffness, dissipation and sliding friction coefficient
    particleSpecies->setDensity(2550.0);
    particleSpecies->setStiffness(1000.0);
    particleSpecies->setSlidingStiffness(2.0 / 7.0 * particleSpecies->getStiffness());
    particleSpecies->setDissipation(0.0);
    particleSpecies->setSlidingFrictionCoefficient(0.12);
 
    //set the parameters for the wall species
    wallParticleSpecies->setStiffness(979.38);
    wallParticleSpecies->setDissipation(particleSpecies->getDissipation());
    wallParticleSpecies->setSlidingStiffness(2. / 7.0 * wallParticleSpecies->getStiffness());
    wallParticleSpecies->setSlidingFrictionCoefficient(0.22);
 
    //set the output parameters: how many time steps should be skipped before a next time step is written to the files,
    //and which files should be written.
    problem.setSaveCount(100);
    problem.dataFile.setFileType(FileType::ONE_FILE);
    problem.restartFile.setFileType(FileType::ONE_FILE);
    problem.fStatFile.setFileType(FileType::NO_FILE);
    problem.eneFile.setFileType(FileType::NO_FILE);
    //
    //problem.autoNumber();
    problem.setXBallsAdditionalArguments("-solidf -v0");
 
    //run the simulation with all parameters given before
    problem.solve(argc, argv);
    
    return 0;
}

References BaseHandler< T >::copyAndAddObject(), DPMBase::dataFile, DPMBase::eneFile, DPMBase::fStatFile, SpeciesHandler::getMixedObject(), INFO, SmoothChute::insertTime, SmoothChute::insertTimeInterval, logger, NO_FILE, SmoothChute::numberOfParticlesToBeInserted, ONE_FILE, DPMBase::restartFile, Chute::setChuteAngleAndMagnitudeOfGravity(), File::setFileType(), DPMBase::setName(), DPMBase::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::setXBallsAdditionalArguments(), DPMBase::solve(), and DPMBase::speciesHandler.