Se plantea este estudio para determinar un modelo matemático que describa la cinética de degradación del paracetamol para su aplicación en reactores fotocatalíticos a gran escala, teniendo en cuenta, aspectos de diseño, estudio de emisión, y de cinética de reacción. A partir de un análisis de literatura de los modelos Generalizado, Zalazar, Ballari y Langmuir-Hinshelwood modificado, para elegir los más adecuados en función de su desempeño con otras sustancias, su error calculado, y parámetros ajustables; se estableció que los modelos adecuados para simular fueron L-H modificado y Generalizado. El valor de la tasa volumétrica global de absorción de fotones (OVRPA) para dar inicio a la simulación fue estimada con gráficos de la tasa volumétrica de absorción de fotones (VRPA) [3]. Se usó el módulo de simulación y evaluación de reactores fotocatalíticos computacional Photoreac V1.0 [2], con base a datos experimentales [1] para analizar la viabilidad de los modelos usando diferentes concentraciones (41,64; 87,6; 149,8 ppm). Se emplearon dos criterios de selección: análisis de dispersión de parámetros, y de residuales en una formula ponderada a: 70% para dispersión de parámetros y un 30% para análisis de residuales; resultando el modelo L-H modificado como el mejor.