Cello Polydisperse drag

[Nimajhi]

The implementation of the Sauter diameter is based on a correction that the cello drag force has.
Obviously there are bottlenecks in this code to increase speed compared to monodisperse mode.

But for other polydisperse drag forces there are different coefficients as you can see in the picture. So what is your suggestion for adding other models in the same format or different?

[PhasicFlow]

First tell me what y and x are standing for in all the equations?
y: ratio of particle diameter to average particle diameter in the cell?
x: fraction of particle of type i in the cell?

[Nimajhi]

yeah, exactly

[PhasicFlow]

I think you first need to evaluate the r.h.s. of the equations for each class/size of particle in each cell and then apply it to all particles.

For doing this, you need to have particle classification field. When in the base class, the average particle diameter is being calculated, this field should also be updated. Then having this field, you know which size class the particle belongs to and then the volume fraction of the particle in that cell is used in the calculation of the r.h.s of equation.

[PhasicFlow]

I added a new method to couplingSystem, shapeDiameters(), which returns a vector of shape diameters. the size is equal to number of particle diameters in the system and it contains the diameter values of all the shapes in the system. Using it, you can classify them based on their size as you are calculating average diameter.