Asymmetry parameters of the phase function for isolated and densely packed spherical particles with multiple internal inclusions in the geometric optics limit