impact.cpp File Reference

#include <sstream>
#include <iostream>
#include <iomanip>
#include <cmath>
#include "scr/ChuteWithHopper.h"
#include <sys/types.h>
#include <sys/stat.h>

Classes
class	GranularJet

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

Function Documentation

main()

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

{
    GranularJet problem;
    
    // Problem parameters
    problem.setName("impact");
    //Should be 10 for full length problem, but here I keep it low for a test case
    problem.setTimeMax(25);
    
    
    // Particle properties
    problem.setDensity(2400.0);
    
    problem.setInflowParticleRadius(0.5,0.5);
    problem.speciesHandler.getObject(0)->setCollisionTimeAndRestitutionCoefficient(4e-3, 0.6);
    problem.speciesHandler.getObject(0)->setSlidingDissipation(problem.get_dissipation());
    problem.speciesHandler.getObject(0)->setSlidingFrictionCoefficient(0.8);
    problem.setFixedParticleRadius(0.0);
    problem.setRoughBottomType(MONOLAYER_ORDERED);
    problem.setAlignBase(false);
    
    // Chute properties
    problem.setChuteAngle(0.0);
    problem.setChuteLength(100.0);
    problem.setChuteWidth(100.0);
    problem.setMaxFailed(6);
    //problem.makeChutePeriodic();
    problem.setHopperDimension(2);
    problem.setHopperLift(20);
    double ExitHeight = 8.0, ExitLength = 1.0 * ExitHeight, hopperAngle_ = 30.0, hopperLength_ = 6.0 * ExitLength;
    problem.setIsHopperCentred(true);
    problem.setHopperLowerFillingHeight(.75);
    problem.setAlignBase(false);
    problem.set_Hopper(ExitLength,ExitHeight,hopperAngle_,hopperLength_);
    problem.MassFlowFactor=20;
    //solve
    cout << "Maximum allowed speed of particles: " << problem.getMaximumVelocity() << endl; // speed allowed before particles move through each other!
    problem.setTimeStep();
    //problem.setSaveCount(helpers::getSaveCountFromNumberOfSavesAndTimeMaxAndTimeStep(100,problem.getTimeMax(),problem.getTimeStep()));
    problem.setSaveCount(500);
    cout << "dt=" << problem.getTimeStep() << endl;
    
    problem.setXBallsAdditionalArguments("-sort -v0 -solidf -v0 -oh -200 -p 20");
    
    problem.auto_number();
    problem.readArguments(argc, argv);
    problem.solve();
    problem.write(cout);
    //problem.HGRID_base::write(cout);
    
    
    //cout << problem << endl;
    problem.writeRestartFile();
}

References e(), MONOLAYER_ORDERED, and problem.