This paper presents a novel Voronoi-TCAD co-simulation framework to investigate the impact of extrinsic and intrinsic variation sources on the performance and variability of 2D-layered TFTs. By focusing on grain boundaries, trap density, grain size, and layer-to-layer interactions, we provide insights of degradation mechanisms, and suggest process optimizations for enhanced device performance and lower variability of 2DTFTs.