4 word modelType = particleCloud.modelType();
7 if (modelType==
"Bfull"){
8 Info <<
"solving volume averaged Navier Stokes equations of type B\n"<< endl;
12 forAll(particleCloud.forceModels(),i)
14 if(particleCloud.forceModels()[i]==
"Archimedes")
18 FatalError <<
"Archimedes model not found!\n" << abort(FatalError);
22 forAll(particleCloud.forceModels(),i)
24 if(particleCloud.forceModels()[i]==
"gradPForce")
28 FatalError <<
"gradPForce model not found!\n" << abort(FatalError);
32 forAll(particleCloud.forceModels(),i)
34 if(particleCloud.forceModels()[i]==
"viscForce")
38 FatalError <<
"viscForce model not found!\n" << abort(FatalError);
40 }
else if(modelType==
"B"){
41 Info <<
"solving volume averaged Navier Stokes equations of type B\n"<< endl;
45 forAll(particleCloud.forceModels(),i)
47 if(particleCloud.forceModels()[i]==
"Archimedes")
51 FatalError <<
"Archimedes model not found!\n" << abort(FatalError);
57 forAll(particleCloud.forceModels(),i)
59 if(particleCloud.forceModels()[i]==
"gradPForce" || particleCloud.forceModels()[i]==
"viscForce")
63 FatalError <<
"do not use gradPForce and viscForce with model type B!\n" << abort(FatalError);
65 }
else if (modelType==
"A"){
66 Info <<
"solving volume averaged Navier Stokes equations of type A\n"
67 <<
"check if Archimedes lift is necessary (depends on pressure formulation)\n"<< endl;
71 forAll(particleCloud.forceModels(),i)
73 if(particleCloud.forceModels()[i]==
"Archimedes")
77 Info <<
"WARNING!!! Are you sure to use Archimedes lift, not necessary if gradPForce includes this effect." << endl;
81 forAll(particleCloud.forceModels(),i)
83 if(particleCloud.forceModels()[i]==
"gradPForce")
87 FatalError <<
"gradPForce model not found!\n" << abort(FatalError);
91 forAll(particleCloud.forceModels(),i)
93 if(particleCloud.forceModels()[i]==
"viscForce")
97 FatalError <<
"viscForce model not found!\n" << abort(FatalError);
98 }
else if (modelType==
"none")
99 Warning <<
"You chose model type -none- you might get erroneous results!" << endl;
101 FatalError <<
"no suitable model type specified:" << modelType <<
"\n" << abort(FatalError);