Very recently, a hepatitis B epidemic model with saturated incidence rate has been proposed and analyzed in Khan et al. (Chaos Solitons Fractals 124:1-9, 2019). The local asymptotic stability of the disease endemic equilibrium (DEE) point of model has been established theoretically but its global asymptotic stability has not been studied. In this paper, we present a mathematically rigorous analysis for the global asymptotic stability of the DEE point of the model. More precisely, we prove that if the DEE point exists, then it is not only locally asymptotically stable but also globally asymptotically stable. Furthermore, we present an alternative proof for the global stability of the disease-free equilibrium point. The main result is that we obtain the complete global stability of the hepatitis B virus model. Besides, we construct nonstandard finite difference (NSFD) schemes preserving the essential qualitative properties of the continuous model. These properties include the positivity of solutions and the stability of the model. Dynamical properties of the proposed schemes are rigorously investigated by mathematical analyses and numerical simulations. Finally, numerical simulations are performed to confirm the validity of the theoretical results as well as the advantages of the NSFD schemes. Numerical simulations also indicate that the NSFD schemes are appropriate and effective to solve the continuous model. Meanwhile, the standard finite difference schemes such as the Euler scheme, the classical fourth-order Runge-Kutta scheme cannot preserve the essential properties of the continuous model. Consequently, they can generate numerical approximations which are completely different from the solutions of the model.