The turbulent flow in an enclosed corotating disk pair: axisymmetric numerical simulation and Reynolds stress modelling