Non-Newtonian fluid flows experience turbulent regime in some industrial applications. Several approaches have been proposed for numerical simulation of turbulent flows that each one has specific features. RANS turbulence models have reasonable computational costs, while include several sources of uncertainties affecting simulation results. In addition, developed RANS models for non-Newtonian fluids are modified versions of available models for Newtonian fluids, therefore, they cannot provide reliable estimation for viscoplastic stress term. On the contrary, DNS delivers accurate results but with high computational costs. Consequently, use of DNS data for estimation of uncertainty in RANS models can provide better decision making for engineers based on RANS results. In the present study, a turbulence model based on for power-law non-Newtonian fluid is developed and employed for simulation of flow in a pipe. Then, an efficient method is proposed for quantification of available model-form uncertainty. Moreover, it is assumed that uncertainties originating from various sources are combined together in calculation of Reynolds stress as well as viscoplastic stress. Deviation of the stresses, computed using RANS turbulence model, from DNS data are modeled through Gaussian Random Field. Thereafter, Karhunen-Loeve expansion is employed for uncertainty propagation in simulation process. Finally, the effects of these uncertainties on RANS results are shown in velocity field demonstrating the fact that the presented approach is accurate enough for statistical modeling of model-form uncertainty in RANS turbulence models.