DEM Forces

Submitted by dusweave on Fri, 10/15/2021 - 20:48

Hi,

I'm having somewhat of an issue understanding the forces acting on particles from the fluid in full CFDEM simulation.
When look at the codes for drag forces the particle volume is calculated by using
ds = 2*particleCloud.radius(index);
Vs = ds*ds*ds*M_PI/6;

Where other forces such as pressure gradient, viscosity, and Archimedes all use
particleVolume().

Tracing particleVolume() function -> forceModel.H ->cfdemCloud.h -> cfdemCloudI.h. This is where I find that the code uses the void fraction model to calculate the total amount of solids inside of each cell. Therefore in my mind we are calculating total amount of solids inside of a cell using particleVolume() instead of merely single particle volume. Am I correct that these other forces average out all the forces on particles inside of a single cell and only drag has individual particle forces?
If this is the case then where is the code that divides the force by the number of pieces to get the force for each particle?

Any help with this will be greatly appreciated.