This work reports the numerical investigation of the secondary breakup of non-Newtonian droplets at different Weber Weg and Ohnesorge Ohl numbers. As part of this work, an in-house coupled level set volume of fluid solver is developed based on OpenFOAM libraries. It uses improved curvature calculation techniques like smoothening and the closest point search method. Flow is assumed to be axisymmetric. Approximately 95 different cases were simulated to investigate the effect of Weg and Ohl numbers on secondary breakup for Newtonian, shear-thinning, and shear-thickening fluids. Ohl varies from 10−3 to 101, and, correspondingly, Weg varies from 101 to 103. The non-Newtonian rheology is modeled as a power-law fluid, and the power-law index n ranged from 0.2–1.8. The present work describes the flow field near the droplet and the effects of non-Newtonian parameters and viscosity on the flow field. The various aspects of droplet dynamics like droplet deformation ratio y, deformation rate ẏ, and coefficient of pressure Cp are studied and compared with the internal flow theory. A generalized relation for critical Weber number Weg,cr is proposed for both Newtonian and non-Newtonian fluids and is shown in a phase diagram plot to map the different regimes of secondary droplet breakup.