FreeCoolingAdhesiveDemo.cpp File Reference

#include <iostream>
#include "Mercury2D.h"
#include "Particles/BaseParticle.h"
#include "Species/LinearViscoelasticSlidingFrictionReversibleAdhesiveSpecies.h"
#include "Walls/InfiniteWall.h"

Go to the source code of this file.

Classes
class	my_problem_HGRID
	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 80 of file FreeCoolingAdhesiveDemo.cpp.

References BaseHandler< T >::copyAndAddObject(), Files::dataFile, Files::fStatFile, BaseParticle::getInteractionRadius(), BaseHandler< T >::getObject(), BaseParticle::getRadius(), my_problem_HGRID::N, NO_FILE, ONE_FILE, DPMBase::particleHandler, File::setFileType(), DPMBase::setGravity(), MercuryBase::setHGridMaxLevels(), MercuryBase::setHGridUpdateEachTimeStep(), Files::setName(), Files::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::setXMax(), DPMBase::setYMax(), DPMBase::solve(), and DPMBase::speciesHandler.

{
    std::cout<<"In this file 32^2 particles with the same velocity are placed "
    "in a bi-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;
    
    my_problem_HGRID problem;
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticSlidingFrictionReversibleAdhesiveSpecies());
    species->setDensity(2000);
    species->setDissipation(0.01);
    species->setStiffness(1e4);
    species->setSlidingFrictionCoefficient(0.1);
    species->setDensity(2000);
    species->setAdhesionStiffness(0.05* species->getStiffness());
    species->setAdhesionForceMax(1.1e-4*0.05* species->getStiffness());

    problem.N=25;
    problem.setName("free_cooling_adhesive");
    problem.setXMax(0.0018);
    problem.setYMax(0.0018);
    problem.setGravity(Vec3D(0.0,0.0,0.0));
    problem.setTimeStep(3e-6);
    problem.setSaveCount(50);
    problem.setTimeMax(0.05);
    problem.setHGridMaxLevels(1);
    //problem.setHGridCellOverSizeRatio(1.2);
    problem.setHGridUpdateEachTimeStep(true);

        
    problem.dataFile.setFileType(FileType::ONE_FILE);
    problem.fStatFile.setFileType(FileType::NO_FILE);
    problem.solve();
    std::cout << problem.particleHandler.getObject(0)->getRadius() << std::endl;
    std::cout << problem.particleHandler.getObject(0)->getInteractionRadius() << std::endl;
    return 0;
}