FreeCoolingDemo.cpp File Reference

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

Go to the source code of this file.

Classes
class	FreeCoolingDemoProblem
	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
	)

Now do the same problem with the HGrid turned on, simple done by create a new problem from the old problem.

Definition at line 79 of file FreeCoolingDemo.cpp.

References FreeCoolingDemoProblem::N, DPMBase::setGravity(), MercuryBase::setHGridCellOverSizeRatio(), MercuryBase::setHGridMaxLevels(), MercuryBase::setHGridUpdateEachTimeStep(), Files::setName(), Files::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), and DPMBase::solve().

{
    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;   

    FreeCoolingDemoProblem problem;
    problem.setName("FreeCoolingDemo");

    problem.N=1024;
    problem.setGravity(Vec3D(0.0,0.0,0.0));
    problem.setTimeStep(5e-5);
    problem.setSaveCount(10);
    problem.setTimeMax(0.05);

    problem.setHGridMaxLevels(1);
    problem.setHGridCellOverSizeRatio(1.2);
    problem.setHGridUpdateEachTimeStep(false);
    problem.solve();
}