Computational fluid dynamics (CFD) software was used to construct two-dimensional axisymmetric simulations of the turbulent flow inside 1.5, 5, 9.15 and 10 m tanks containing rotating water. The water rotation was induced mathematically by either placing a rotating cover on the surface of the water or locating a tangential water inlet near the top of the side wall. Central and side drains evacuated the water from the tank when using an inlet. The rotating cover produced a forced vortex in the tanks similar to the one observed in actual rearing tanks when the majority of the flow is leaving via the side drain. The predicted flow structure in tanks with an inlet can be divided into three regions: a narrow region along the side wall where the rotating water is moving down toward the floor, a thin boundary layer along the floor where the rotating water is moving radially inward, and a large region in the bulk of the tank where radial and axial velocities are small and the tangential velocity is independent of elevation but varies with radial position. Calculations were performed with the k-ε (RNG) and Reynolds stress turbulence models. Converged solutions were easier to obtain with the k-ε (RNG) model but only the Reynolds stress model could predict the strong vortex which is experimentally observed to form in the central portion of the tank as the flow leaving via the center drains is increased. The simulations predict that the thickness of the floor boundary layer is proportional to the tank diameter and that the radial velocity within the floor boundary layer is maximum at an elevation equal to about 10% of the boundary layer thickness. At any given radial position, the maximum value of the radial velocity next to the floor is between 15 and 45% of the tangential velocity in the bulk of the tank. Predicted shear stresses along bounding surfaces were used to obtain correlations for the side-wall and floor friction coefficients in terms of the Reynolds number. These correlations were in turn substituted in an overall moment balance to obtain an analytical model for predicting the maximum tangential velocity Vθw near the side wall of multi-drain rearing tanks. The predictions of the analytical model are in excellent agreement with Vθw values reported in the literature for tanks with diameters ranging from 1.5–15 m and indicate that the drag from inlet structures is non-negligible. The results demonstrate that the axisymmetric CFD model can simulate the main features of the rotating flow inside circular tanks and provide valuable boundary-layer information that is difficult to obtain experimentally.