Radiation hardness is a critical design concern for present and future silicon detectors in high energy physics. Tracking systems at the CERN Large Hadron Collider (LHC) are expected to operate for ten years and to receive fast hadron fluences equivalent to 1015cm-2 1-MeV neutrons. Recently, low temperature operating conditions have been suggested as a means of suppressing the negative effects of radiation damage on detector charge collection properties. To investigate this effect, simulations have been carried out using the ISE-TCAD DESSIS device simulator. The so-called "three-level model" has been used. A comprehensive analysis of the influence of the V2, CiOi and V2O capture cross sections on the effective doping concentration (Neff) as a function of temperature and fluence has been carried out. The capture cross sections have been varied in the range 10-18 - 10-12 cm2. The simulated results are compared with charge collection spectra obtained with 1064-nm laser pulses on devices irradiated with 23-GeV protons as a function of detector bias voltage. To validate the model, a wide range of temperature and fluence has been studied using a one-dimensional (1-D) simplified structure. Thousands of simulation results have been cross checked with the experimental data. The data between 190 K (the lower limit for simulations due to computational difficulties) and 290 K are well reproduced for all of the fluences considered. We conclude that the three-level model can be successfully used to predict irradiated detector behavior down to a temperature of at least 190 K.