STLRotationSelfTest.cpp File Reference

#include <Mercury3D.h>
#include <Species/LinearViscoelasticSpecies.h>
#include <Walls/TriangleWall.h>

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

Function Documentation

main()

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

{
    // Problem setup
    Mercury3D problem;
    // Set name of output files
    problem.setName("STLRotationSelfTest");
    // Set domain size
    problem.setMax({450,50,50});
    problem.setMin({-450,-50,-50});
    // Set time step, final time and how often to output
    problem.setTimeStep(0.01);
    problem.setTimeMax(0.1);
    problem.setSaveCount(1);
    // Turn off most output files except data and wall-vtu files
    problem.fStatFile.setFileType(FileType::NO_FILE);
    problem.eneFile.setFileType(FileType::NO_FILE);
    problem.restartFile.setFileType(FileType::NO_FILE);
    problem.wallHandler.setWriteVTK(true);
    // Introduce a material (no properties set, as no collisions happen here)
    problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    // introduce an outer wall
    problem.wallHandler.readTriangleWall("Casing.stl",problem.speciesHandler.getObject(0));
    // introduce an inner, rotating wall
    Mdouble scaleFactor = 1;
    Vec3D shift = {0,0,0};
    Vec3D velocity = {0,0,0};
    Vec3D angularVelocity = {2.0*constants::pi,0,0};
    problem.wallHandler.readTriangleWall("Screw.stl",problem.speciesHandler.getObject(0),
            scaleFactor,shift,velocity,angularVelocity);
    // start solving in time
    problem.solve();
    logger.assert_always(problem.wallHandler.getSize()==4316,"Didn't read the right number of walls");
    logger(INFO,"Load %Wall_*.vtu in paraview to see the wall geometry",problem.getName());
    return 0;
}

References INFO, logger, NO_FILE, constants::pi, problem, and Jeffery_Solution::velocity().