no_seg.cpp File Reference

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <Species/LinearViscoelasticFrictionSpecies.h>
#include "Chute.h"

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

Function Documentation

main()

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

A quasi-2D inclined plane with inflow conditions on the left boundary, and deletion of particles when they exit the domain.

{
    //Print description
    logger(INFO, "\nDescription: A quasi-2D inclined plane with inflow conditions on the left boundary, and deletion of"
                 " particles when they exit the domain.");
    
    // Problem parameters
    Chute problem;
    problem.setName("chute_demo");
    problem.setTimeMax(0.5);
    
    // Particle properties
    problem.setFixedParticleRadius(0.001);
    problem.setInflowParticleRadius(0.001);
    
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
    double mass = 0.5*species->getMassFromRadius(0.5*(problem.getMinInflowParticleRadius() + problem.getMaxInflowParticleRadius()));
    species->setCollisionTimeAndRestitutionCoefficient(2.5e-3,0.8,mass);
    species->setSlidingDissipation(0.0);
    problem.setTimeStep(0.02 * 2.5e-3);
 
    
    // Chute properties
    problem.setChuteAngle(30.0);
    problem.setXMax(0.1);
    problem.setYMax(2.0*problem.getInflowParticleRadius());
    
    // Inflow properties
    problem.setInflowHeight(0.1);
    problem.setInflowVelocity(0.1);
    problem.setInflowVelocityVariance(0.02);
    
    // Output properties
    problem.setSaveCount(helpers::getSaveCountFromNumberOfSavesAndTimeMaxAndTimeStep(1, problem.getTimeMax(), problem.getTimeStep()));
    problem.dataFile.setSaveCount(100*problem.restartFile.getSaveCount());
    problem.eneFile.setSaveCount(100*problem.restartFile.getSaveCount());
//    problem.setSaveCount(helpers::getSaveCountFromNumberOfSavesAndTimeMaxAndTimeStep(1,problem.getTimeMax(),problem.getTimeStep())); //minimize output to the last time step
//    problem.set_number_of_saves_data(100); //allow enough data output so the evolution can be viewed in xballs
//    problem.set_number_of_saves_ene(100);
 
    //solve
    problem.solve(argc, argv);
}

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