We present a simulation framework based on the combination of Monte Carlo transport and Landauer-Buttiker formalism to model and optimize ballistic rectifiers exploiting the long mean-free-path of graphene. The influence on the responsivity (i.e. the induced voltage normalized by the impinging power) of the main geometrical parameters, of bias and of different scattering mechanisms is critically analyzed.