A numerical modeling recipe has been established to perform a parametric simulation on the transport properties of charged species in dye sensitized solar cells (DSSCs). The main outputs of the simulation are the charged species densities and current voltage curve which determine the performance parameters of DSSC. But a realistic estimation of the performance parameters strongly depends on the number of physical parameters involved. To do so, we have attempted to engage the most recent theories in the present work. The results of the simulation allow the optimization of DSSC performance as a function of photoanode thickness. The highest power conversion efficiency belongs to the device with photoanode thickness of 20 mu m. The model also highlights the importance of substrate/electrolyte interface in recombination of carriers. This model can be considered to be a tool to give a comprehensive view of DSSC operation in order to reduce the experimental workload aimed at optimizing DSSC performance.