Computation of the Turbulent Flow in a Square Duct Using a Cubic Low-Reynolds Stress Model