adding forces only to DEM integration (not to momentum exchange terms)

Hi!

Words of warning:
Please note, that not all force models will have this functionality!
More important, current solvers do not support explicit source terms (feel free to add it to the solver - everything necessary is allready there).

please see e.g. Archimedes.C

// add force only to DEM calc, not to momentum exchange terms:
114 if(treatDEM_) for(int j=0;j<3;j++) DEMForces[index][j] += force[j];

// or add force to an explicit source term "f"
115 else if(treatExplicit_) for(int j=0;j<3;j++) expForces[index][j] += force[j];

// or add force to an implicit source term "Ksl"
116 else for(int j=0;j<3;j++) impForces[index][j] += force[j];

so you can choose only one of the options! If not defined it will add it to an implicit source term.

you can activate it via:
94 ArchimedesProps
95 {
96 densityFieldName "rho";
97 gravityFieldName "g";
98 treatDEM;
99 }

Cheers,
Chris