In the framework of the CERN-RD50 Collaboration, the adoption of p-type substrates has been proposed as a suitable mean to improve the radiation hardness of silicon detectors up to fluencies of 1times10 16 n/cm 2. In this work two numerical simulation models will be presented for p-type and n-type silicon detectors, respectively. A comprehensive analysis of the variation of the effective doping concentration (N eff), the leakage current density and the charge collection efficiency as a function of the fluence has been performed using the Synopsys T-CAD device simulator. The simulated electrical characteristics of irradiated detectors have been compared with experimental measurements extracted from the literature, showing a very good agreement. The predicted behaviour of p-type silicon detectors after irradiation up to 10 16 n/cm 2 shows better results in terms of charge …