TriangleWallCubeSelfTest.cpp File Reference

#include "Mercury3D.h"
#include "Walls/TriangulatedWall.h"
#include "Walls/InfiniteWall.h"
#include "Species/LinearViscoelasticSpecies.h"
#include "Logger.h"

Functions
int	main ()
	Tests the implementation of TriangulatedWall. More...

Function Documentation

main()

int main

(

)

Tests the implementation of TriangulatedWall.

{
    //write input files
    helpers::writeToFile("TriangulatedWallSimple.vtk","# vtk DataFile Version 2.0\n"
     "Defines a domain outside of a tetrahedron\n"
     "ASCII\n"
     "DATASET UNSTRUCTURED_GRID\n"
     "POINTS 4 double\n"
     "0 0 0\n"
     "1 0 0\n"
     "0 1 0\n"
     "0 0 1\n"
     "\n"
     "CELLS 4 16\n"
     "3 0 1 2\n"
     "3 0 2 3\n"
     "3 0 3 1\n"
     "3 1 3 2\n"
     "\n"
     "CELL_TYPES 4\n"
     "5\n"
     "5\n"
     "5\n"
     "5");
 
    helpers::writeToFile("TriangulatedWallInverted.vtk","# vtk DataFile Version 2.0\n"
     "Defines a domain inside of a tetrahedron\n"
     "ASCII\n"
     "DATASET UNSTRUCTURED_GRID\n"
     "POINTS 4 double\n"
     "0 0 0\n"
     "1 0 0\n"
     "0 1 0\n"
     "0 0 1\n"
     "\n"
     "CELLS 4 16\n"
     "3 0 2 1\n"
     "3 0 3 2\n"
     "3 0 1 3\n"
     "3 1 2 3\n"
     "\n"
     "CELL_TYPES 4\n"
     "5\n"
     "5\n"
     "5\n"
     "5");
 
    helpers::writeToFile("TriangulatedWall.vtk","# vtk DataFile Version 2.0\n"
     "Cube Normals (cross of first two vec's) into the wall\n"
     "ASCII\n"
     "DATASET UNSTRUCTURED_GRID\n"
     "POINTS 8 double\n"
     "0 0 0\n"
     "1 0 0\n"
     "1 1 0\n"
     "0 1 0\n"
     "0 0 1\n"
     "1 0 1\n"
     "1 1 1\n"
     "0 1 1\n"
     "\n"
     "CELLS 8 32\n"
     "3 0 1 2\n"
     "3 0 7 4\n"
     "3 2 6 7\n"
     "3 0 2 7\n"
     "3 1 0 4\n"
     "3 1 4 7\n"
     "3 6 2 1\n"
     "3 1 7 6\n"
     "\n"
     "CELL_TYPES 8\n"
     "5\n"
     "5\n"
     "5\n"
     "5\n"
     "5\n"
     "5\n"
     "5\n"
     "5");
 
    logger(INFO,"To get VTK output, uncomment the relevant lines in the main function.");
 
 
    Mercury3D dpm;
    dpm.setName("TriangleWallCubeSelfTest");
    dpm.setMin({-.2,-.2,-.2});
    dpm.setMax({1.2,1.2,1.2});
 
    dpm.setGravity(10*Vec3D(0,0,-1));
    dpm.setTimeStep(1e-4);
    dpm.setTimeMax(1e4*dpm.getTimeStep());
    //setTimeMax(getTimeStep());
    dpm.setSaveCount(100);
    dpm.setXBallsAdditionalArguments("-v0 -solidf");
    dpm.fStatFile.setFileType(FileType::NO_FILE);
    dpm.restartFile.setSaveCount(1e5);
 
    Mdouble radius = 0.08;
    auto s = dpm.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    s->setDensity(1);
    s->setCollisionTimeAndRestitutionCoefficient(20.0*dpm.getTimeStep(),0.1,s->getMassFromRadius(radius));
 
    //create domain walls
    InfiniteWall i;
    i.setSpecies(dpm.speciesHandler.getLastObject());
    i.set({0,0,-1},dpm.getMin());
    dpm.wallHandler.copyAndAddObject(i);
    i.set({0,-1,0},dpm.getMin());
    dpm.wallHandler.copyAndAddObject(i);
    i.set({-1,0,0},dpm.getMin());
    dpm.wallHandler.copyAndAddObject(i);
    i.set({0,1,0},dpm.getMax());
    dpm.wallHandler.copyAndAddObject(i);
    i.set({1,0,0},dpm.getMax());
    dpm.wallHandler.copyAndAddObject(i);
 
    //read in walls
    dpm.wallHandler.readTriangleWall("TriangulatedWall.vtk",s);
 
    //introduce particles
    SphericalParticle p(s);
 
    //random insertion
    for (unsigned i=0; i<300; i++)
    {
        Vec3D r;
        r.X = dpm.random.getRandomNumber(dpm.getXMin(),dpm.getXMax());
        r.Y = dpm.random.getRandomNumber(dpm.getYMin(),dpm.getYMax());
        r.Z = dpm.random.getRandomNumber(dpm.getZMax(),dpm.getZMax()+1.0);
        p.setPosition(r);
        p.setRadius(dpm.random.getRandomNumber(radius,1.1*radius));
        if (dpm.checkParticleForInteraction(p))
            dpm.particleHandler.copyAndAddObject(p);
    }
 
    logger(INFO,"Inserted % particles",dpm.particleHandler.getNumberOfObjects());
    //uncomment to get VTK output
    dpm.setParticlesWriteVTK(true);
    dpm.wallHandler.setWriteVTK(FileType::ONE_FILE);
    dpm.solve();
    return 0;
}

References MercuryBase::checkParticleForInteraction(), BaseHandler< T >::copyAndAddObject(), DPMBase::fStatFile, BaseHandler< T >::getLastObject(), DPMBase::getMax(), DPMBase::getMin(), ParticleHandler::getNumberOfObjects(), RNG::getRandomNumber(), DPMBase::getTimeStep(), DPMBase::getXMax(), DPMBase::getXMin(), DPMBase::getYMax(), DPMBase::getYMin(), DPMBase::getZMax(), constants::i, INFO, logger, NO_FILE, ONE_FILE, DPMBase::particleHandler, DPMBase::random, WallHandler::readTriangleWall(), DPMBase::restartFile, File::setFileType(), DPMBase::setGravity(), DPMBase::setMax(), DPMBase::setMin(), DPMBase::setName(), DPMBase::setParticlesWriteVTK(), BaseInteractable::setPosition(), BaseParticle::setRadius(), DPMBase::setSaveCount(), File::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), WallHandler::setWriteVTK(), DPMBase::setXBallsAdditionalArguments(), DPMBase::solve(), DPMBase::speciesHandler, DPMBase::wallHandler, helpers::writeToFile(), Vec3D::X, Vec3D::Y, and Vec3D::Z.