I want to add the Beetstra drag force for a polydisperse system. Implementing the drag force is easy, but the problem is with the "Mean Sauter Diameter" that needs to be in the Reynolds calculation. For the Mean Sauter Diameter, we need to access the solid volume fraction of each particle and the total of them, and the diameter of each particle.
Please give some guidance in this regard