The dataset include mean flow data of square duct flow up to friction Reynolds number 1000 obtained from direct numerical simulation. The dataset also include coefficients obtained from Pope's generalized eddy viscosity framework, for different truncation of the polynomial expansion.