Classes | Functions
CGLiquidsUnitTest.cpp File Reference
#include <Walls/InfiniteWall.h>
#include <Species/LinearViscoelasticFrictionLiquidMigrationWilletSpecies.h>
#include <Particles/LiquidFilmParticle.h>
#include <CG/Fields/LiquidMigrationFields.h>
#include "Mercury3D.h"
#include "CG/CG.h"

Classes

class  TwoParticles
 

Functions

int main ()
 

Function Documentation

◆ main()

int main ( )
43 {
44  logger(INFO," Simulates a particle-particle collision.\n"
45  " Checks cg \n"
46  " - at one instant in time\n"
47  " - standard fields (density and stress)\n"
48  " - for O, Z, and XYZ coordinates\n"
49  " - for both Gauss and Lucy kernel functions\n");
50 
51  TwoParticles dpm;
52  dpm.setTimeMax(0);
53  auto o = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::O,CGFunctions::Gauss,CGFields::LiquidMigrationFields>());
54  auto z = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::Z,CGFunctions::Gauss,CGFields::LiquidMigrationFields>());
55  z->setWidth(0.16);
56  z->setN(100);
57  auto xz = dpm.cgHandler.copyAndAddObject(CG<CGCoordinates::XZ,CGFunctions::Gauss,CGFields::LiquidMigrationFields>());
58  xz->setWidth(0.16);
59  xz->setN(100);
60  dpm.solve();
61 
62  //A few checks
63  logger(INFO,"Checking a few cg parameters");
64 
65  //Checks density for unresolved cg, rho = M/V = 2/16 = 0.125
66  helpers::check(o->getPoint(0).getLiquidFilmVolume(),1,1e-15, "Average liquidFilmVolume");
67  helpers::check(o->getPoint(0).getLiquidBridgeVolume(),1,1e-15, "Average liquidBridgeVolume");
68  helpers::check(z->evaluateAverage().getLiquidFilmVolume(),1,3e-2, "Average liquidFilmVolume for z-resolved stats");
69  helpers::check(z->evaluateAverage().getLiquidBridgeVolume(),1,1e-2, "Average liquidBridgeVolume for z-resolved stats");
70  helpers::check(xz->evaluateAverage().getLiquidFilmVolume(),1,3e-2, "Average liquidFilmVolume for xz-resolved stats");
71  helpers::check(xz->evaluateAverage().getLiquidBridgeVolume(),1,1e-2, "Average liquidBridgeVolume for xz-resolved stats");
72 
73  return 0;
74 }
e
Array< double, 1, 3 > e(1./3., 0.5, 2.)
logger
Logger< MERCURYDPM_LOGLEVEL > logger("MercuryKernel")
Definition of different loggers with certain modules. A user can define its own custom logger here.
BaseHandler::copyAndAddObject
std::enable_if<!std::is_pointer< U >::value, U * >::type copyAndAddObject(const U &object)
Creates a copy of a Object and adds it to the BaseHandler.
Definition: BaseHandler.h:360
CG
Evaluates time-resolved continuum fields and writes the data into a stat file.
Definition: CG.h:55
DPMBase::cgHandler
CGHandler cgHandler
Object of the class cgHandler.
Definition: DPMBase.h:1479
DPMBase::setTimeMax
void setTimeMax(Mdouble newTMax)
Sets a new value for the maximum simulation duration.
Definition: DPMBase.cc:864
DPMBase::solve
void solve()
The work horse of the code.
Definition: DPMBase.cc:4334
TwoParticles
Definition: CGLiquidsUnitTest.cpp:14
INFO
#define INFO(i)
Definition: mumps_solver.h:54
helpers::check
void check(double real, double ideal, double error, std::string errorMessage)
Definition: TestHelpers.cc:16

References DPMBase::cgHandler, helpers::check(), BaseHandler< T >::copyAndAddObject(), e(), INFO, logger, DPMBase::setTimeMax(), and DPMBase::solve().