We present a calculation methodology of classical Hall mobility based on the non-equilibrium Green's function method, and calculate the Hall mobility in two-dimensional systems without and with potential disorders. By using the Büttiker probe method, inelastic scattering is taken into account. The Hall mobility in systems without disorders agrees with the expectation by the scattering rate, validating our methodology. Then, the Hall mobility in systems with disorders is simulated, and the behavior is discussed by comparing with the results by Fermi's golden rule, a scattering theory, and an analysis based on an average effective mass calculated considering the disorders. Reasonable agreement between the Hall mobility by the non-equilibrium Green's function method and the estimation by the average effective mass is obtained, suggesting that the average effective mass can be used to estimate the Hall mobility in systems with disorders.