A non-linear k-epsilon model for turbulent shear flows