Flow of non-Newtonian fluids through porous media at high Reynolds numbers is often encountered in chemical, pharmaceutical and food, as well as petroleum and groundwater engineering, and in many other industrial applications. Under the majority of operating conditions typically explored, the dependence of pressure drops on flow rate is non-linear and the development of models capable of describing accurately this dependence, in conjunction with non-trivial rheological behaviors, is of paramount importance. In this work, pore-scale single-phase flow simulations conducted on synthetic two-dimensional porous media are performed via computational fluid dynamics for both Newtonian and non-Newtonian fluids and the results are used for the extension and validation of the Darcy-Forchheimer law, herein proposed for shear-thinning fluid models of Cross, Ellis and Carreau. The inertial parameter beta is demonstrated to be independent of the viscous properties of the fluids. The results of flow simulations show the superposition of two contributions to pressure drops: one, strictly related to the non-Newtonian properties of the fluid, dominates at low Reynolds numbers, while a quadratic one, arising at higher Reynolds numbers, is dependent on the porous medium properties. The use of pore-scale flow simulations on limited portions of the porous medium is here proposed for the determination of the macroscale-averaged parameters (permeability K, inertial coefficient beta and shift factor alpha), which are required for the estimation of pressure drops via the extended Darcy-Forchheimer law. The method can be applied for those fluids which would lead to critical conditions (high pressures for low permeability media and/or high flow rates) in laboratory tests.