Computational modelling of turbulent Taylor–Couette flow for bearing chamber applications: A comparison of unsteady Reynolds-averaged Navier–Stokes models

Abstract

The capability to accurately model fluid flow within rotating Taylor–Couette systems has a primary role in informing computational investigations of rotating machinery. There is considerable uncertainty regarding selection of modelling approach, including a suitable turbulence model, that can accurately resolve turbulence within such complex flows while remaining computationally feasible for industrially relevant applications. This paper presents a numerical comparison of axisymmetric and three-dimensional unsteady Reynolds-averaged Navier–Stokes (URANS) turbulence models within ANSYS Fluent. The CFD geometries are representative of ones for which there are published experimental measurements. For the Taylor–Couette study, investigation into inner cylinder start-up procedure, based on previous published findings, confirmed that the final state of the flow is highly dependent on the initial conditions and acceleration rate. Once Taylor vortices form and stabilise, they are not disrupted by small steps in inner cylinder speed, allowing computationally efficient accelerations. Investigations into applying rotational periodicity were unsuccessful, resulting in a significantly reduced core velocity. Axisymmetric predictions provided reasonable agreement with experimental data only at low rotation rates. A good prediction of the velocity flow field was obtained for three-dimensional simulations of the full 360° domain with differences of less than 5% for radial velocities. Among the URANS models, the standard k-ω model and baseline Reynolds stress model (BSL-RSM) provided the closest agreement to published experimental data. In the paper, the developed Taylor–Couette turbulence modelling methodology is extended to a bearing chamber geometry. Analysis of the secondary vortex flow field is compared both qualitatively and quantitatively to published bearing chamber experimental measurements. Overall, whilst a good agreement is still found using the standard k-ω turbulence model, discrepancies arise with the BSL-RSM. However, for this more complex bearing chamber environment compared to a Taylor–Couette flow, the shear stress transport k-ω turbulence model provided the closest agreement and is recommended for future bearing chamber modelling.