FluxBoundaryUnitTest.cpp File Reference

#include "Mercury3D.h"
#include "Species/LinearViscoelasticSpecies.h"
#include "Boundaries/FluxBoundary.h"

Classes
class	FluxBoundaryUnitTest
	Tests if FluxBoundary correctly measures the flux of particles crossing the boundary. 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
    FluxBoundaryUnitTest problem;
    problem.setTimeMax(1);
    problem.setTimeStep(0.1);
    problem.setFileType(FileType::NO_FILE);
    
    //run the simulation
    problem.solve(argc, argv);
 
    // Check the value of the fluxBoundary
    FluxBoundary* fluxBoundary = static_cast<FluxBoundary*>(problem.boundaryHandler.getObject(0));
 
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedNet() == 1, "Wrong number of particles crossed the flux boundary");
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedForw() == 2, "Wrong number of particles crossed the flux boundary (Forward)");
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedForw() == 2, "Wrong number of particles crossed the flux boundary (Backward)");
 
    // Get and ckeck volume of particles
    Mdouble particleVolume = problem.particleHandler.getObject(0)->getVolume();
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedNet()*particleVolume == fluxBoundary->getVolumeOfParticlesCrossedNet(), "Wrong volume of particles crossed the flux boundary");
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedForw()*particleVolume == fluxBoundary->getVolumeOfParticlesCrossedForw(), "Wrong volume of particles crossed the flux boundary (Forward)");
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedBack()*particleVolume == fluxBoundary->getVolumeOfParticlesCrossedBack(), "Wrong volume of particles crossed the flux boundary (Backward)");
 
    // Get and check mass of particles
    Mdouble particleMass = problem.particleHandler.getObject(0)->getMass();
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedNet()*particleMass == fluxBoundary->getMassOfParticlesCrossedNet(), "Wrong mass of particles crossed the flux boundary");
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedForw()*particleMass == fluxBoundary->getMassOfParticlesCrossedForw(), "Wrong mass of particles crossed the flux boundary (Forward)");
    logger.assert_always(fluxBoundary->getNumberOfParticlesCrossedBack()*particleMass == fluxBoundary->getMassOfParticlesCrossedBack(), "Wrong mass of particles crossed the flux boundary (Backward)");
}

References FluxBoundary::getMassOfParticlesCrossedBack(), FluxBoundary::getMassOfParticlesCrossedForw(), FluxBoundary::getMassOfParticlesCrossedNet(), FluxBoundary::getNumberOfParticlesCrossedBack(), FluxBoundary::getNumberOfParticlesCrossedForw(), FluxBoundary::getNumberOfParticlesCrossedNet(), FluxBoundary::getVolumeOfParticlesCrossedBack(), FluxBoundary::getVolumeOfParticlesCrossedForw(), FluxBoundary::getVolumeOfParticlesCrossedNet(), logger, NO_FILE, and problem.