A three-dimensional mixing length model is proposed for modeling local Reynolds stress behavior in rectangular ducts of arbitrary aspect ratio. The model is applicable to both developing and fully-developed flows, and can be applied to other 90-degree corner flows with mild streamwise pressure gradients. Comparisons between theory and experiment show that all components of the Reynolds stress tensor are modeled reasonably well, both in the vicinity of a corner and in two-dimensional regions away from the corner.