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 BaseHandler< T >::copyAndAddObject(), Chute::getChuteAngle(), Chute::getFixedParticleRadius(), Chute::getMaxInflowParticleRadius(), Chute::getMinInflowParticleRadius(), DPMBase::getTimeStep(), DPMBase::getXMax(), DPMBase::getXMin(), DPMBase::getYMax(), DPMBase::getYMin(), INFO, logger, MONOLAYER_DISORDERED, Funnel::set_funa(), Funnel::set_funD(), Funnel::set_funfr(), Funnel::set_funHf(), Funnel::set_funnz(), Funnel::set_funO(), Chute::setChuteAngle(), Chute::setChuteLength(), Chute::setChuteWidth(), Chute::setFixedParticleRadius(), Chute::setInflowParticleRadius(), Chute::setInflowVelocity(), Chute::setInflowVelocityVariance(), Chute::setMaxFailed(), DPMBase::setName(), Chute::setRoughBottomType(), DPMBase::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), mathsFunc::sin(), DPMBase::solve(), DPMBase::speciesHandler, and DPMBase::writeRestartFile().