We calculate the scattering cross section between two 0 þþ glueballs in SUð2Þ Yang-Mills theory on a lattice at β ¼ 2.1, 2.2, 2.3, 2.4, and 2.5 using the indirect (HAL QCD) method. We employ the clusterdecomposition 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 the 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 Λ as σ ϕϕ ¼ ð2-51ÞΛ −2 (stat þ sys). From the observational constraints of galactic collisions, we obtain the lower bound of the scale parameter as Λ > 60 MeV. We also discuss the naturalness of the Yang-Mills theory as the theory explaining dark matter.