CGDisplacementFieldUnitTest.cpp File Reference

#include "Mercury3D.h"
#include "Species/LinearViscoelasticSpecies.h"
#include "CG/TimeAveragedCG.h"
#include "CG/Fields/DisplacementField.h"

Classes
class	CGDisplacementFieldUnitTest
	Tests if the displacement momentum density is correctly calculated. More...

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

Function Documentation

main()

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

{
    //declare the DPM problem and set the name
    CGDisplacementFieldUnitTest problem;
    problem.setTimeMax(1);
    problem.setTimeStep(0.1);
    problem.setFileType(FileType::NO_FILE);
 
    //define different coarse-graining objects (which should all result in the same mean values)
    auto cg0 = problem.cgHandler.copyAndAddObject(TimeAveragedCG<CGCoordinates::O,CGFunctions::Lucy,CGFields::DisplacementField>());
 
    // Initially, the previous position is not set, so the displacement would not make sense at time step 0
    cg0->setTimeMin(0.2);
    cg0->selectSpecies(0);
    
    //run the simulation
    problem.solve(argc, argv);
 
    Vec3D displacementMomentum = cg0->evaluateAverage().getDisplacementMomentum();
 
    // Calculate the expected displacement momentum
    BaseParticle* p0 = problem.particleHandler.getObject(0);
    Vec3D diff = problem.getMax() - problem.getMin();
    Mdouble simulationVolume = diff.X * diff.Y * diff.Z;
    Vec3D expectedDisplacementMomentum = p0->getMass() * p0->getVelocity() / simulationVolume;
 
    logger.assert_always((displacementMomentum-expectedDisplacementMomentum).getLength() < 1e-10, "Displacement momentum is not correct");
}

References e(), logger, NO_FILE, p0, problem, Vec3D::X, Vec3D::Y, and Vec3D::Z.