We have extended Monte Carlo simulations of hopping transport in completely disordered 2D conductors to the process of external charge relaxation. In this situation, a conductor of area L × W shunts an external capacitor C with initial charge Q i . At low temperatures, the charge relaxation process stops at some "residual" charge value corresponding to the effective threshold of the Coulomb blockade of hopping. We have calculated the r.m.s. value Q R of the residual charge for a statistical ensemble of capacitor-shunting conductors with random distribution of localized sites in space and energy and random Q i , as a function of macroscopic parameters of the system. Rather unexpectedly, Q R has turned out to depend only on some parameter combination: X 0 ≡ LW ν 0 e 2 /C for negligible Coulomb interaction and X χ ≡ LW κ 2 /C 2 for substantial interaction. (Here ν 0 is the seed density of localized states, while κ is the dielectric constant.) For sufficiently large conductors, both functions Q R /e = F (X) follow the power law F (X) = DX −β , but with different exponents: β = 0.41 ± 0.01 for negligible and β = 0.28 ± 0.01 for significant Coulomb interaction. We have been able to derive this law analytically for the former (most practical) case, and also explain the scaling (but not the exact value of the exponent) for the latter case. In conclusion, we discuss possible applications of the sub-electron charge transfer for "grounding" random background charge in single-electron devices.