The exceptional performance of the silicon sensors in the radiation environment has led to their extensive application in high energy physics. Even so, the future experiments foresee these sensors to be exposed to higher radiation levels. Radiation induces a change in the macroscopic properties of the sensor, thus, severely affecting the sensor performance and ultimately becoming the limiting factor for its operation. With an aim to extend the radiation hardness capabilities of the silicon sensors for the future experiments there has been a growing interest in sensors with novel designs and unique characteristic of intrinsic charge multiplication. \newline
However, it is important to understand the effect of radiation damage on these sensors, before employing them in the main detector system. The RD50 collaboration extensively employs TCAD simulation tools for an in-depth understanding and structural optimization of the newly proposed sensor technologies, complementing the measurement results. The simulation package, by the finite element method, solves the Poisson equation coupled with the current continuity equations in the main device, typically with the drift-diffusion current model. There is also the possibility of co-simulating an electronic circuit connected to the main device using the incorporated SPICE package, for example to predict single event effects or read-out related waveforms. The simulation tools also provide an insight into the sensor operation both in the non-irradiated and the irradiated scenario in order to predict the voltage dependence of the leakage current, charge collection and electrical field distribution as function of irradiation. This has required the development of a radiation damage model within the simulation tools such that the measurements are well reproduced. The details of the radiation damage modelling using two commercial TCAD tools $–$ Silvaco and Synopsys, are discussed in this work.