Droplet formation in a flow-focusing device was modeled using the open source CFD package, OpenFOAM, with the VOF model for two-phase flow. Predictions using theinterFoamsolver and a power-law non-Newtonian model were first validated against experimental data in the literature. Then, the formation of Newtonian fluid droplets in power-law fluids was modeled during tubing, squeezing, dripping and jetting. The effects of the continuous phase rheological parameters on the droplet formation were investigated by changing the power law index (n) and the consistency coefficient (K). The results show that the droplet length and the spacing between two droplets decrease asnorKincrease. However, the formation frequency and droplet velocity in the main channel increase asnorKincrease. The results also show thatnhas a greater effect thanKon the droplet formation. A method was developed to calculate the capillary number of the power-law continuous phase in the squeezing and dripping regimes including the influences ofnandK. For a given dispersed phase flow rate, the formation frequency is inversely proportional to the droplet volume. A scaling law was also developed to predict the formation frequency since the droplet volume is found to vary linearly with the non-dimensional droplet length. The present work is useful for controlling droplet formation and designing microfluidic devices in areas where non-Newtonian fluids are used as the continuous phase.