  We calculate the scattering cross section between two $0^{++}$ glueballs in
$SU(2)$ Yang-Mills theory on lattice at $\beta = 2.1, 2.2, 2.3, 2.4$, and 2.5
using the indirect (HAL QCD) method. We employ the cluster-decomposition error
reduction technique and use all space-time symmetries to improve the signal. In
the use of the HAL QCD method, the centrifugal force was subtracted to remove
the systematic effect due to nonzero angular momenta of lattice discretization.
From the extracted interglueball potential we determine the low energy glueball
effective theory by matching with the one-glueball exchange process. We then
calculate the scattering phase shift, and derive the relation between the
interglueball cross section and the scale parameter $\Lambda$ as $\sigma_{\phi
\phi} = (2 - 51) \Lambda^{-2}$ (stat.+sys.). From the observational constraints
of galactic collisions, we obtain the lower bound of the scale parameter, as
$\Lambda > 60$ MeV. We also discuss the naturalness of the Yang-Mills theory as
the theory explaining dark matter.
