CGParticleWallUnitTest.cpp File Reference

#include <Walls/InfiniteWall.h>
#include "Mercury3D.h"
#include "CG/CG.h"
#include "Species/LinearViscoelasticSpecies.h"

Classes
class	ParticleWall

Functions
int	main ()

Function Documentation

main()

int main

(

)

Todo:

TW should there be a discrepancy between distance and distance2?

{
    logger(INFO," Simulates a particle-wall collision.\n"
            " Checks cg \n"
            "   - at one instant in time\n"
            "   - standard fields (density and stress)\n"
            "   - for O, Z, and XYZ coordinates\n"
            "   - for both Gauss and Lucy kernel functions\n");
 
    ParticleWall dpm;
    dpm.setName("CGParticleWallUnitTest");
    dpm.setTimeStep(1e-4);
    dpm.setTimeMax(0);
    dpm.setFileType(FileType::NO_FILE);
 
    //Define coarse graining objects
    Mdouble width = 0.075;
 
    //Unresolved
    auto O = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::O>());
 
    //Resolved in z, at particle center
    auto Z = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::Z>(width,1));
    Z->setMin({-1,-1,0.95});
    Z->setMax({1,1,0.95});
 
    //Resolved in z with Gauss cg function
    auto G = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::Z,CGFunctions::Gauss>(width,1));
    G->setMin({-1,-1,0.95});
    G->setMax({1,1,0.95});
 
    //Resolved in z, on contact line
    auto C = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::Z>(width,1));
    C->setMin({-1,-1,0.25});
    C->setMax({1,1,0.25});
 
    //Resolved in xyz, at particle center
    auto XYZ = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::XYZ>(width,1));
    XYZ->setMin({-0.95,-1,0.95});
    XYZ->setMax({1.05,1,0.95});
 
    dpm.solve();
 
    //A few checks
    logger(INFO,"Checking a few cg parameters");
 
    //Checks if the unresolved cg returns the right value for stress, sigma_zz = f_12*r_12/V
    Vec3D contactPoint = dpm.interactionHandler.getObject(0)->getContactPoint();
    Mdouble overlap = dpm.interactionHandler.getObject(0)->getOverlap();
    Mdouble distance = dpm.interactionHandler.getObject(0)->getDistance();
    Mdouble distance2 = Vec3D::getDistance(contactPoint,dpm.particleHandler.getObject(0)->getPosition());
    Mdouble distance3 = distance+0.5*overlap;
    logger(INFO,"Distance % % % Overlap % Contact % Distance %",distance,distance2,distance3,overlap,contactPoint);
    helpers::check(O->getPoint(0).getContactStress().trace(),2e5*overlap*distance2/16,1e-11, "Average stress");
 
    //Checks if z-resolved cg returns the right value for stress at the center of the contact,
    //sigma_zz = f_12*r_12/r_12z/A, assuming the cutoff is smaller than the branch vector
    Mdouble distanceZ = dpm.particleHandler.getObject(0)->getPosition().Z;
    helpers::check(C->getPoint(0).getContactStress().trace(),2e5*overlap*distance/distanceZ/4,1e-11, "Z-resolved stress at contact point");
 
 
    return 0;
}

References DPMBase::cgHandler, helpers::check(), BaseHandler< T >::copyAndAddObject(), e(), G, BaseInteraction::getContactPoint(), BaseInteraction::getDistance(), Vec3D::getDistance(), BaseHandler< T >::getObject(), BaseInteraction::getOverlap(), BaseInteractable::getPosition(), INFO, DPMBase::interactionHandler, logger, NO_FILE, DPMBase::particleHandler, DPMBase::setFileType(), DPMBase::setName(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::solve(), XYZ, Z, and Vec3D::Z.