6
I have a simple program like this:
length = 10;
o = {0, 0, 0};
dp = 0.25 (*m*);
rho = 2345 (*kg/m^3*);
rho = dp^3*rho;
region = Cuboid[o, o + length];
cavity = Table[Ball[RandomPoint[region], RandomReal[{0.6 dp, 1.5 dp}]], 4];
region = Fold[RegionDifference, region, cavity];
RegionPlot3D[region,...
