Aims. In interstellar clouds, molecular hydrogens are formed from atomic hydrogen on grain surfaces. An atomic hydrogen hops around until it finds another one with which it combines. This necessarily implies that the average recombination time, or equivalently, the effective grain surface area depends on the relative numbers of atomic hydrogen influx rate and the number of sites on the grain. Our aim is to discover this dependency. Methods. We perform a numerical simulation to study the recombination of hydrogen on grain surfaces in a variety of cloud conditions. We use a square lattice (with a periodic boundary condition) of various sizes on two types of grains, namely, amorphous carbon and olivine. Results. We find that the steady state results of our simulation match very well with those obtained from a simpler analytical considerations provided the "effective" grain surface area is written as ∼ S α , where S is the actual physical grain area and α is a function of the flux of atomic hydrogen, which is determined from our simulation. We carry out the simulation for various astrophysically relevant accretion rates. For high accretion rates, small grains tend to become partly saturated with H and H 2 , and the subsequent accretion will be partly inhibited. For very low accretion rates, the number of sites to be swept before a molecular hydrogen can form is too large compared to the actual number of sites on the grain, implying that α is greater than unity.