SegregationNew.cpp File Reference

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cmath>
#include "Species/LinearViscoelasticFrictionSpecies.h"
#include "Chute.h"
#include "Walls/InfiniteWall.h"
#include "Boundaries/PeriodicBoundary.h"
#include "Boundaries/CubeInsertionBoundary.h"
#include <CG/CG.h>

Classes
class	Chutebelt
	if you restart this code the third argument will be used as the number of large particles to add and the forth the number of small. More...

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

Function Documentation

main()

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

Todo:

change this to the mean particle size (not small particle size)

{
    //Print description
    logger(INFO, "\nDescription: A quasi-2D moving-bed channel with walls on the left and right boundary.");
    
    // Problem parameters
    Chutebelt problem;
 
    problem.autoNumber();
    problem.setName("Segregation");
    problem.setTimeMax(2000.0);
    problem.setTimeStep(1. / (200.0 * 50.0));
    
    problem.set_radiusLarge(2.0);
    problem.set_particle_number_volRatio(1.0); //volume ratio of large to small
    problem.set_particle_numbers(2500, 130); //Should be 5000 and 130
 
    problem.setChuteAngleAndMagnitudeOfGravity(25.0, 1.0);
    problem.set_beltSpeed(0.0);
    
    // Chute properties : Simply remove the first line to add side walls.
    problem.makeChutePeriodic();
    problem.setXMax(150.0);
    problem.setYMax(5.0);
    
    problem.setZMin(0.0);
    problem.setZMax(250.0);
 
    problem.setFixedParticleRadius(0.75);
    
    //Swap the next two lines to swap between the different type of rought bottoms.
    problem.setRoughBottomType(MULTILAYER);
    //problem.setRoughBottomType(MONOLAYER_DISORDERED);
 
    problem.setSaveCount(helpers::getSaveCountFromNumberOfSavesAndTimeMaxAndTimeStep(2000, problem.getTimeMax(), problem.getTimeStep()));
    
    problem.setXBallsAdditionalArguments("-cmode 7");
    
    problem.readArguments(argc, argv);
    
    if (argc > 4)
    {
        problem.num_restart_large=atoi(argv[3]);
        problem.num_restart_small=atoi(argv[4]);
    }
    else
    {
        problem.num_restart_large=0;
        problem.num_restart_small=0;
    }
    
    logger(INFO, "% %", problem.num_restart_small, problem.num_restart_large);
 
 
    problem.solve();
 
    
}

References helpers::getSaveCountFromNumberOfSavesAndTimeMaxAndTimeStep(), INFO, logger, MULTILAYER, and problem.