FreeCooling3DDemo.cpp File Reference

#include <iostream>
#include <Species/LinearViscoelasticSpecies.h>
#include "Mercury3D.h"
#include "Particles/BaseParticle.h"
#include "Walls/InfiniteWall.h"
#include "MercuryTime.h"

Go to the source code of this file.

Classes
class	FreeCooling3DDemoProblem
	todo{This code is not working as is wanted} More...

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

Function Documentation

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

Start off my solving the default problem

Definition at line 86 of file FreeCooling3DDemo.cpp.

References SpeciesHandler::addObject(), FreeCooling3DDemoProblem::N, DPMBase::setGravity(), MercuryBase::setHGridCellOverSizeRatio(), MercuryBase::setHGridMaxLevels(), MercuryBase::setHGridUpdateEachTimeStep(), Files::setName(), Files::setSaveCount(), DPMBase::setSystemDimensions(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::setXMax(), DPMBase::setYMax(), DPMBase::setZMax(), DPMBase::solve(), DPMBase::speciesHandler, Time::tic(), and Time::toc().

{
    std::cout << "In this file 10^3 particles with the same velocity are placed "
    "in a tri-axial box. This makes them collide with the walls and eachother. "
    "Afterwards the same run is performed with hgrid on. It tests the working "
    "(and speedup) of the hgrid." << std::endl;
    
    Time time;
    time.tic();
    FreeCooling3DDemoProblem problem;
    auto species = new LinearViscoelasticSpecies;
    problem.speciesHandler.addObject(species);
    species->setDensity(2000);
    species->setDissipation(0.005);
    species->setStiffness(1e3);

    problem.setXMax(0.1);
    problem.setYMax(0.1);
    problem.setZMax(0.1);
    problem.N=1000;
    problem.setName("FreeCooling3DDemo");
    problem.setGravity(Vec3D(0.0,0.0,0.0));
    problem.setTimeStep(5e-5);
    problem.setSaveCount(4);
    problem.setTimeMax(0.08);
    problem.setSystemDimensions(3);
    
    problem.setHGridMaxLevels(1);
    problem.setHGridCellOverSizeRatio(1.2);
    problem.setHGridUpdateEachTimeStep(false);    
    problem.solve();
    
    std::cout << "Total time to run this simulation: " << time.toc() << "s" << std::endl;
}