This work reports the numerical investigation of the secondary breakup of non-Newtonian droplets at different Weber W e g and Ohnesorge O h l 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 W e g and O h l numbers on secondary breakup for Newtonian, shear-thinning, and shear-thickening fluids. O h l varies from 10 - 3 to 10 1, and, correspondingly, W e g varies from 10 1 to 10 3. 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 y ˙, and coefficient of pressure C p are studied and compared with the internal flow theory. A generalized relation for critical Weber number W e g, c r 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. © 2022 Author(s).