code: def time_col_turbulent(EnergyDis, ConcAl, ConcClay, coag, material,
                         DiamTarget, DIM_FRACTAL):
      
      return((1/6) * (6/np.pi)**(1/9) * EnergyDis**(-1/3) * DiamTarget**(2/3)
             * frac_vol_floc_initial(ConcAl, ConcClay, coag, material)**(-8/9)
             * (DiamTarget / material.Diameter)**((8*(DIM_FRACTAL-3)) / 9)
             )
docs: Calculate single collision time for turbulent flow mediated collisions.
  
      Calculated as a function of floc size.