MonodisperseSmoothInclinedChutePatternInflow.cpp File Reference

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

Classes
class	SmoothChute

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

Function Documentation

main()

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

{
    SmoothChute problem;
    problem.setName("FlowOverSmoothInclinedChannelPatternInflow");
    problem.setTimeMax(0.001);
    problem.setChuteAngleAndMagnitudeOfGravity(30.0,9.81);
    
    problem.insertTime = 0.01;
    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;
 
    particleSpecies->setDensity(2550.0); // sets the species type-0 density
    particleSpecies->setStiffness(1000.0);// sets the spring stiffness
    particleSpecies->setSlidingStiffness(2.0/7.0*particleSpecies->getStiffness());
    particleSpecies->setDissipation(0.0);// sets the dissipation
    particleSpecies->setSlidingFrictionCoefficient(0.12);
 
    wallParticleSpecies->setStiffness(979.38);
    wallParticleSpecies->setDissipation(particleSpecies->getDissipation());
    wallParticleSpecies->setSlidingStiffness(2.0/7.0*wallParticleSpecies->getStiffness());
    wallParticleSpecies->setSlidingFrictionCoefficient(0.22);
 
    //For this problem, we want a datafile and a restartfile as output, and
    // they should be printed every 100 time steps.
    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");
 
    problem.setTimeStep(1e-4);
    problem.solve(argc, argv);
    
}

References e(), NO_FILE, ONE_FILE, and problem.