This paper presents a groundbreaking numerical approach, the fractional differential quadrature method (FDQM), to simulate the complex dynamics of organic polymer solar cells. The method, which leverages polynomial-based differential quadrature and Cardinal sine functions coupled with the Caputo-type fractional derivative, offers a significant improvement in accuracy and efficiency over traditional methods. By employing a block-marching technique, we effectively address the time-dependent nature of the governing equations. The efficacy of the proposed method is validated through rigorous numerical simulations and comparisons with existing analytical and numerical solutions. Each scheme’s computational characteristics are tailored to achieve high accuracy, ensuring an error margin on the order of 10−8 or less. Additionally, a comprehensive parametric study is conducted to investigate the impact of key parameters on device performance. These parameters include supporting conditions, time evolution, carrier mobilities, charge carrier densities, geminate pair distances, recombination rate constants, and generation efficiency. The findings of this research offer valuable insights for optimizing and enhancing the performance of organic polymer solar cell devices.