GranularJet.cpp File Reference

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <Species/LinearViscoelasticSlidingFrictionSpecies.h>
#include "Funnel.h"

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

Function Documentation

main()

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

{
    Funnel problem;
    
    // Particle properties
    problem.setFixedParticleRadius(.3e-3);
    problem.setInflowParticleRadius(.3e-3);
    
    //Contact properties
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSlidingFrictionSpecies());
    species->setSlidingFrictionCoefficient(0.5);
    species->setDensity(1442.0); //sand
    species->setCollisionTimeAndRestitutionCoefficient(4e-4, 1.0 - (0.12 * (50.0 / 15.0)), species->getMassFromRadius(
            problem.getInflowParticleRadius()));//eps=0.6
    species->setSlidingStiffness(species->getStiffness() * 2.0 / 7.0);
    species->setSlidingDissipation(species->getDissipation() * 2.0 / 7.0);
    
    logger(INFO, "Setting k to % and disp to % with radius: %",
           species->getStiffness(), species->getDissipation(), problem.getInflowParticleRadius());
    
    // Default problem parameters
    problem.setChuteLength(0.25);
    problem.setChuteWidth(0.15);
    problem.setChuteAngle(26.7);
    problem.setRoughBottomType(MONOLAYER_DISORDERED);
    
    //funr; // Funnel radius.
    //funO[2]; // (x,y) location of the center of the top of the funnel.
    //funa; // Angle of the funnel
    //funH; // Heigth of the funnel
    //funHf; // Falling heigth
    //funD; // Funnel diameter at the downside of the funnel.
    //funnz; //Number of particles along the heigth of the funnel
    //funfr; // Filling ratio of the funnel
    //fundiag; //The diagonal of the filling region
    //funrmax; //The maximum range for r
    problem.set_funa(60.0); // Angle of the funnel
    problem.set_funD(0.015); // Funnel diameter at the downside of the funnel.
    problem.set_funHf(0.05); // Falling heigth
    problem.set_funnz(25.0); //Number of particles along the heigth of the funnel
    //problem.set_funO(0, 0.5*(problem.getYMax()+problem.getYMin())); // (x,y) location of the center of the top of the funnel.
    problem.set_funfr(0.3); // Filling ratio of the funnel
    
    problem.setInflowVelocity(0);
    problem.setInflowVelocityVariance(0.01);
    problem.setMaxFailed(1);
    
    logger(INFO, "Chute inflow height: % Chute inflow velocity: % Chute inflow velocity variance: %",
           problem.getInflowHeight(), problem.getInflowVelocity(), problem.getInflowVelocityVariance());
    
    //Discretization parameters
    //problem.setHGridMaxLevels(1);
    //problem.setHGridNumberOfBucketsToPower(1e6); //automated
    double mass = species->getMassFromRadius(
            0.5 * (problem.getMinInflowParticleRadius() + problem.getMaxInflowParticleRadius()));
    problem.setTimeStep(0.02 * species->getCollisionTime(mass));
    problem.setTimeMax(.5e-2);
    problem.setSaveCount(10);
    
    logger(INFO, " dt = % tmax = %", problem.getTimeStep(), problem.getTimeMax());
    //  logger(INFO, "Maximum allowed speed of particles: %", problem.getMaximumVelocity()); // speed
    //  allowed before particles move through each other!
    
    
    problem.readArguments(argc, argv);
    problem.set_funO(-0.9 * problem.get_funHf() * sin(problem.getChuteAngle()), 0.5 * (problem.getYMax() +
                                                                                       problem.getYMin())); // (x,y) location of the center of the top of the funnel.
    problem.setName_();
    problem.solve();
}

References e(), INFO, logger, MONOLAYER_DISORDERED, problem, and sin().