Numerical modeling of a Global Navigation Satellite System in a general relativistic framework