Numerical study of the steady-state uniform flow past a rotating cylinder

Padrino, J.C. and Joseph, D.D. (2006) Numerical study of the steady-state uniform flow past a rotating cylinder. Journal of Fluid Mechanics, 557, 191-223. (doi:10.1017/S0022112006009682).

Abstract

Results from the numerical simulation of the two-dimensional incompressible unsteady Navier-Stokes equations for streaming flow past a rotating circular cylinder are presented in this study. The numerical solution of the equations of motion is conducted with a commercial computational fluid dynamics package which discretizes the equations applying the control volume method. The numerical set-up is validated by comparing results for a Reynolds number based on the free stream of Re = 200 and dimensionless peripheral speed of q̃ = 3, 4 and 5 with results from the literature. After the validation stage, various pairs of Re q̃ are specified in order to carry out the numerical experiments. These values are Re = 200 with q̃ = 4 and 5; Re = 400 with q̃ = 4, 5 and 6, and Re = 1000 with q̃ = 3. In all these cases, gentle convergence to fully developed steady state is reached. From the numerical vorticity distribution, the position of the outer edge of the vortical region is determined as a function of the angular coordinate. This position is found by means of a reasonable criterion set to define the outmost curve around the cylinder where the vorticity magnitude reaches a certain cut-off value. By considering the average value of this profile, a uniform vortical region thickness is specified for every pair of Re and q̃. Next, the theoretical approach of Wang & Joseph (2006a; see the companion paper) and the numerical results are used to determine two different values of the effective vortical region thickness for every pair of Re and q̃. One effective thickness δ_D/a is obtained from the match between the additional drag on the outer edge of the vortical region according to the viscous correction of viscous potential flow (VCVPF) and the corresponding numerical profile while the other thickness δ_L/a is determined from the match between the pressure lift on the cylinder obtained from Wang & Joseph (2006a)'s simple modification of the boundary-layer analysis due to Glauert (Proc. R. Soc. Lond. A, vol. 242, 1957, p. 108) and the numerical value of the pressure lift coefficient. The values of δ_D/a and δ_L/a are used in the computation of various parameters associated with the flow, namely, the torque on the rotating cylinder, the circulatory velocity at the edge of the vortical region, which links the cylinder's angular velocity with the circulation of the irrotational flow of the viscous fluid outside this region, and the viscous dissipation. Predictions from the approaches of Glauert (1957) and Wang & Joseph (2006a) are also included for comparison. The values of both effective thicknesses, δ_D/a and δ_L/a, are found to have the same order of magnitude. Then, we show that choosing δ_D/ a as a unique effective thickness, the modification of Glauert's boundary-layer analysis and the VCVPF approach as proposed by Wang & Joseph (2006a) produce results which are in better general agreement with the values from numerical simulation than those from Glauert's solution.

This record has no associated files available for download.

Contributors

Author: J.C. Padrino

Download statistics