dec13_A254_Hi0075_RC06_MU05.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;
    
    // Problem parameters
    problem.setName("dec13_A254_Hi0075_RC06_MU05");
 
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSlidingFrictionSpecies());
    species->setSlidingFrictionCoefficient(0.5);
    species->setDensity(1442.0);
 
    //problem.set_HGRID_max_levels(2);
    //problem.set_HGRID_num_buckets(1e6);
    // Particle properties
 
    problem.setFixedParticleRadius(300e-6);
    problem.setInflowParticleRadius(300e-6);
 
    species->setCollisionTimeAndRestitutionCoefficient(4e-4, 0.6, species->getMassFromRadius(problem.getFixedParticleRadius()));
 
 
    //problem.setStiffness(1e5*4/3*pi*problem.getInflowParticleRadius()*9.81*problem.getDensity());
    //problem.set_disp(50*sqrt(9.81/(2*problem.getInflowParticleRadius());
    logger(INFO, "Setting k to % and disp to %", species->getStiffness(), species->getDissipation());
 
    species->setSlidingStiffness(species->getStiffness()*2.0/7.0);
    species->setSlidingDissipation(species->getDissipation()*2.0/7.0);
 
    double mass = species->getMassFromRadius(0.5 * (problem.getMinInflowParticleRadius() + problem.getMaxInflowParticleRadius()));
    problem.setTimeStep(0.02 * species->getCollisionTime(mass));
    problem.setTimeMax(3);
    problem.setSaveCount(100);
    
    problem.setChuteLength(0.6);
    problem.setChuteWidth(0.25);
    problem.setChuteAngle(25.4);
    problem.setRoughBottomType(MONOLAYER_DISORDERED);
    
    
    double funx = problem.getXMin() + 0.5 * (problem.getXMax() - problem.getXMin());
    
    
    problem.set_funa(60.0);
    problem.set_funD(0.015);
    problem.set_funHf(0.075 + (problem.getXMax() - funx) * sin(problem.getChuteAngle()));
    problem.set_funnz(50.0);
    problem.set_funO(-funx, 0.5 * (problem.getYMax() - problem.getYMin()));
    problem.set_funfr(0.3);
    
    problem.setInflowVelocity(0);
    problem.setInflowVelocityVariance(0.01);
    problem.setMaxFailed(1);
    
    
    
    //solve
    //cout << "Maximum allowed speed of particles: " << problem.getMaximumVelocity() << endl; // speed allowed before particles move through each other!
    
    
    logger(INFO, "dt=%", problem.getTimeStep());
    problem.solve();
    problem.writeRestartFile();
}

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