Spark plasma sintering (SPS) is one of the most widely used sintering techniques that utilizes pulsed direct current together with uniaxial pressure to consolidate a wide variety of materials. The unique mechanisms of SPS enable it to sinter powder compacts at a lower temperature and in a shorter time than the conventional hot pressing, hot isostatic pressing and vacuum sintering process. One of the limitations of SPS is the presence of temperature gradients inside the sample, which could result in non-uniform physical and microstructural properties. Detailed study of the temperature and current distributions inside the sintered sample is necessary to minimize the temperature gradients and achieve desired properties. In the present study, a coupled thermal-electric model was developed using finite element codes in ABAQUS software to investigate the temperature and current distributions inside the conductive and non-conductive samples. An integrated experimental-numerical methodology was implemented to determine the system contact resistances accurately. The developed sintering model was validated by a series of experiments, which showed good agreements with simulation results. The temperature distribution inside the sample depends on some process parameters such as sample and tool geometry, punch and die position, applied current and thermal insulation around the die. The role of these parameters on sample temperature distribution was systematically analyzed. The findings of this research could prove very useful for the reliable production of large size sintered samples with controlled and tailored properties.