Chute/segregation.cpp File Reference

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cmath>
#include "scr/Chute.h"
#include <sys/types.h>
#include <sys/stat.h>

Classes
class	SegregationPeriodic
	This class does segregation problems in a periodic chute. More...

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

Function Documentation

main()

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

{
    SegregationPeriodic problem;
    
    // Problem parameters
    problem.setName("segregation");
    problem.setTimeMax(100);
    
    
    // Particle properties
    problem.setDensity(2400.0);
    
    problem.setInflowParticleRadius(0.3e-3,0.60e-3);
    problem.speciesHandler.getObject(0)->setCollisionTimeAndRestitutionCoefficient(4e-4, 0.6);
    problem.speciesHandler.getObject(0)->setSlidingDissipation(problem.get_dissipation());
    problem.speciesHandler.getObject(0)->setSlidingFrictionCoefficient(0.8);
    problem.setFixedParticleRadius(0.3e-3);
    problem.setRoughBottomType(MONOLAYER_DISORDERED);
    
    
    // Chute properties
    problem.setChuteAngle(25.0);
    problem.setChuteLength(50.0e-3);
    problem.setChuteWidth(3e-3);
    problem.setMaxFailed(6);
    problem.makeChutePeriodic();
    
    //solve
    cout << "Maximum allowed speed of particles: " << problem.getMaximumVelocity() << endl; // speed allowed before particles move through each other!
    problem.setTimeStep(); 
    //This is based on the fact in general you get too much data, so prob at worst you want to turn it into a 20 at 60fps (which is its self overkill)
    problem.setSaveCount(helpers::getSaveCountFromNumberOfSavesAndTimeMaxAndTimeStep(20*60,getTimeMax(),getTimeStep()));
    //problem.setSaveCount(1);
    cout << "dt=" << problem.getTimeStep() << endl;
    
    problem.auto_number();
 
    
    //This set to colouring based of size and small vectors
    problem.setXBallsColourMode(7);
    problem.setXBallsVectorScale(1);
    problem.setXBallsAdditionalArguments("-sort -v0 -solidf");
    
    //problem.setSaveCount(1000);
    
    problem.solve();
    //problem.HGRID_base::write(cout);
    
    
    //cout << problem << endl;
    problem.writeRestartFile();
}

References e(), helpers::getSaveCountFromNumberOfSavesAndTimeMaxAndTimeStep(), MONOLAYER_DISORDERED, and problem.