Identification of complex flows in Taylor–Couette counter-rotating cavities

The transition in confined rotating flows is a topical problem with many industrial and fundamental applications. The purpose of this study is to investigate the Taylor–Couette flow in a finite-length cavity with counter-rotating walls, for two aspect ratios L=5 or L=6. Two complex regimes of wavy vortex and spirals are emphasized for the first time via direct numerical simulation, by using a three-dimensional spectral method. The spatio-temporal behavior of the solutions is analyzed and compared to the few data actually available.     

Stress gradient-induced migration effects in the Taylor–Couette flow of a dilute polymer solution

We present an investigation of the phenomenon of stress-induced polymer migration for dilute polymer solutions in the Taylor–Couette device, consisting of two infinitely long, concentric cylinders rotating at constant angular velocities. The underlying physical model is represented by the dilute limit of a two-fluid Hamiltonian system involving two components: one (the polymer) is viscoelastic and obeys the Oldroyd-B constitutive equation, and the other (the solvent) is viscous Newtonian. The two components are considered to be in thermal, but not mechanical equilibrium, interacting with each other through an isotropic drag coefficient tensor. This allows for stress-induced diffusion of polymer chains. The governing equations consist of the continuity and the momentum equations for the bulk velocity, the constitutive model for the polymer chain conformation tensor and the diffusion equation for the polymer concentration. The diffusion equation contains an extra source term, which is proportional to gradients in the polymer stress, so that polymer concentration gradients can develop even in the absence of externally imposed fluxes in the presence of stress inhomogeneities. The solution to the steady-state purely azimuthal flow is obtained first using a spectral collocation method and an adaptive mesh formulation to track the steep changes of the concentration in the flow domain. The calculations show the development of strong polymer migration towards the inner cylinder with increasing Deborah number (De) in agreement with experimental observations. The migration is enhanced for increasing values of the gap thickness resulting in concentration changes by several orders of magnitude in the area between the inner and outer cylinder walls. The extent of the migration also depends strongly on the ratio of the solvent to the polymer viscosity. In addition to a strongly inhomogeneous polymer concentration, significant deviations from the homogenous flow are also observed in the velocity profile. Next, results are reported from a linear stability analysis around the steady-state solution against axisymmetric disturbances corresponding to various wavenumbers in the axial direction. The calculations show that the steady-state solution remains stable up to moderate values of the Deborah number, explaining why some of the predicted stress-induced migration effects should be experimentally observable. The role of the Peclet number (Pe) on the stability of the system is elucidated.     

Turbulent Couette–Taylor flows with endwall effects: A numerical benchmark

The accurate prediction of fluid flow within rotating systems has a primary role for the reliability and performance of rotating machineries. The selection of a suitable model to account for the effects of turbulence on such complex flows remains an open issue in the literature. This paper reports a numerical benchmark of different approaches available within commercial CFD solvers together with results obtained by means of in-house developed or open-source available research codes exploiting a suitable Reynolds Stress Model (RSM) closure, Large Eddy Simulation (LES) and a direct numerical simulation (DNS). The predictions are compared to the experimental data of Burin et al. (2010) in an original enclosed Couette–Taylor apparatus with endcap rings. The results are discussed in details for both the mean and turbulent fields. A particular attention has been turned to the scaling of the turbulent angular momentum G with the Reynolds number Re. By DNS, G is found to be proportional to Re^α, the exponent α = 1.9 being constant in our case for the whole range of Reynolds numbers. Most of the approaches predict quite well the good trends apart from the k–ω SST model, which provides relatively poor agreement with the experiments even for the mean tangential velocity profile. Among the RANS models, even though no approach appears to be fully satisfactory, the RSM closure offers the best overall agreement.     

Numerical study of eccentric Couette–Taylor flows and effect of eccentricity on flow patterns

In this study, the differential quadrature (DQ) method was used to simulate the eccentric Couette–Taylor vortex flow in an annulus between two eccentric cylinders with rotating inner cylinder and stationary outer cylinder. An approach combining the SIMPLE (semi-implicit method for pressure-linked equations) and DQ discretization on a non-staggered mesh was proposed to solve the time-dependent, three-dimensional incompressible Navier–Stokes equations in the primitive variable form. The eccentric steady Couette–Taylor flow patterns were obtained from the solution of three-dimensional Navier–Stokes equations. The reported numerical results for steady Couette flow were compared with those from Chou [1], and San and Szeri [2]. Very good agreement was achieved. For steady eccentric Taylor vortex flow, detailed flow patterns were obtained and analyzed. The effect of eccentricity on the eccentric Taylor vortex flow pattern was also studied.     

Nonnormal Effects in the Taylor–Couette Problem

This work is devoted to the study of transient growth of perturbations in the Taylor–Couette problem due to linear nonnormal mechanisms. The study is carried out for a particular small gap case and is mostly focused on the linearly stable regime of counter-rotation. The exploration covers a wide range of inner and outer angular speeds as well as axial and azimuthal modes. Significant transient growth is found in the regime of stable counter-rotation. The numerical results are in agreement with former analyses based on energy methods and other independent numerical studies. The optimal energy transient growth factor appears to be consistent with experimental observations. This study might shed some light on the subcritical transition to turbulence which is found experimentally in Taylor–Couette flow when the cylinders rotate in opposite directions. Received 13 February 2001 and accepted 29 March 2002 Published online: 2 October 2002 RID="*" ID="*" This work was supported by the UK EPSRC under Grant GR/M30890. The author thanks Nick Trefethen for fruitful discussions. RID="*" ID="*" Present address: Departament de Fisica Aplicada, Univ. Politecnica de Catalunya, 08034 Barcelona, Spain (alvar@fa.upc.es) Communicated by H.J.S. Fernando     

Flow pattern and stability of pseudoplastic axial Taylor–Couette flow

The effect of an axial flow on the stability of the Taylor–Couette flow is explored for pseudoplastic fluids. The fluid is assumed to follow the Carreau–Bird model and mixed boundary conditions are imposed while the axial flow can be independent of rotational flow. The four-dimensional low-order dynamical system, resulted from Galerkin projection of the conservation of mass and momentum equations, includes additional non-linear terms in the velocity components originated from the shear-dependent viscosity. In absence of axial flow the base flow loses its radial flow stability to the vortex structure at a lower critical Taylor number, as the pseudoplasticity effects increases. The emergence of the vortices corresponds to the onset of a supercritical bifurcation which is also seen in the flow of a linear fluid. However, unlike the Newtonian case, pseudoplastic Taylor vortices lose their stability as the Taylor number reaches a second critical number corresponding to the onset of a Hopf bifurcation. Existence of an axial flow, induced by a pressure gradient appears to further advance each critical point on the bifurcation diagram. Complete flow field together with viscosity maps are given for stability regions in the bifurcation diagram.     

Resonances in the Codimension-2 Bifurcations in the Couette–Taylor Problem

The investigation of codimension-2 bifurcations, in particular in systems with cylindric symmetry, enables us to deduce new types of secondary regimes branching-off from the symmetric regimes. This investigation also allows us the unique possibility of a rigorous treatment of chaotic solutions to Navier–Stokes and other nonlinear PDE’s. The central manifold approach combined with the reduction to the normal form lead to the so-called amplitude systems. These ODE systems describe the nonlinear interaction between the neutral modes, and always include several nonlinear terms due to so-called intrinsic resonances. However, sometimes additional resonances appear. In this paper we present the complete list of all possible resonances in dynamic systems with cylindric symmetry and the corresponding forms of the amplitude equations. Further, we present the results of extensive numerical investigation of the resonant codimension-2 bifurcations in the Couette–Taylor problem, thus creating an intriguing subject for further investigation.     

Stability of plane Poiseuille–Couette flows of a piezo-viscous fluid

We examine stability of fully developed isothermal unidirectional plane Poiseuille–Couette flows of an incompressible fluid whose viscosity depends linearly on the pressure as previously considered in Hron et al. [J. Hron, J. Málek, K.R. Rajagopal, Simple flows of fluids with pressure-dependent viscosities, Proc. R. Soc. Lond. A 457 (2001) 1603–1622] and Suslov and Tran [S.A. Suslov, T.D. Tran, Revisiting plane Couette–Poiseuille flows of a piezo-viscous fluid, J. Non-Newtonian Fluid Mech. 154 (2008) 170–178]. Stability results for a piezo-viscous fluid are compared with those for a Newtonian fluid with constant viscosity. We show that piezo-viscous effects generally lead to stabilisation of a primary flow when the applied pressure gradient is increased. We also show that the flow becomes less stable as the pressure and therefore the fluid viscosity decrease downstream. These features drastically distinguish flows of a piezo-viscous fluid from those of its constant-viscosity counterpart. At the same time the increase in the boundary velocity results in a flow stabilisation which is similar to that observed in Newtonian fluids with constant viscosity.     

Simulations of radiative effects on the Rayleigh–Taylor instability using the CRASH code

Future experiments at the National Ignition Facility will be able to generate diagnosable Rayleigh–Taylor instability growth in the presence of locally generated, high radiation fluxes. This interplay of radiative energy transfer and hydrodynamic instability is relevant to many astrophysical systems, such as core-collapse red supergiant supernovae. Previous simulations of high-energy-density Rayleigh–Taylor instabilities in the presence of a hot environment near a radiative shock demonstrate behavior that differs from that found in non-radiative cases. However, these simulations considered only 1D or single wavelength cases. Here we report simulations of an entire experimental system using the CRASH code. These simulations lead to modified predictions, attributed to the effects of radial energy losses.     

Bubble capture and migration in Couette–Taylor flow

 Bubble capture and migration under the effect of organized structures in weak turbulent Couette–Taylor flow between two concentric cylinders, the inner one rotating, has been investigated. Bubbles generated at the free surface for large enough angular velocities are sucked into the flow by the upper organized structures. Then they migrate progressively from top to bottom by jumping from cell to cell. With an upper solid stationary wall instead of the free surface, injected bubbles are trapped by the coherent vortices beyond a critical Taylor number. However, in this situation there is no migration mechanism carrying the bubbles from top to bottom. This particular migration and capture process, able to act against the forces of buoyancy, has been investigated by perturbing the flow by adding a vertical plate protruding from the inner surface of the solid stationary wall. The perturbation so introduced causes the deformation of the upper coherent structures and reinstalls the migration of the bubbles. Received: 27 October 1997/Accepted: 21 May 1998     

Viscous effects on the non-classical Rayleigh–Taylor instability of spherical material interfaces

The viscous and conductivity effects on the instability of a rapidly expanding material interface produced by a spherical shock tube are investigated through the employment of a high-order WENO scheme. The instability is influenced by various mechanisms, which include (a) classical Rayleigh–Taylor (RT) effects, (b) Bell–Plesset or geometry/curvature effects, (c) the effects of impulsively accelerating the interface, (d) compressibility effects, (e) finite thickness effects, and (f) viscous effects. Henceforth, the present instability studied is more appropriately referred to as non-classical RT instability to distinguish it from classical RT instability. The linear regime is examined and the development of the viscous three-dimensional perturbations is obtained by solving a one-dimensional system of partial differential equations. Numerical simulations are performed to illustrate the viscous effects on the growth of the disturbances for various conditions. The inviscid analysis does not show the existence of a maximum amplification rate. The present viscous analysis, however, shows that the growth rate increases with increasing the wave number, but there exists a peak wavenumber beyond which the growth rate decreases with increasing the wave number due to viscous effects.     

On the formation of Taylor–G?rtler vortices in RMF-driven spin-up flows

This paper is concerned with a liquid metal flow driven by a rotating magnetic field inside a stationary cylinder. We consider especially the secondary meridional flow during the time when the fluid spins up from rest. The developing flow is investigated experimentally and by direct numerical simulations. The vertical profiles of the axial velocity are measured by means of the ultrasound Doppler velocimetry. Evolving instabilities in the form of Taylor–G?rtler vortices have been observed just above the instability threshold (Ta ≥ 1.5· Ta ^cr). The rotational symmetry may survive over a distinct time even if a first Taylor–G?rtler vortex pair has been formed as closed rings along the cylinder perimeter. The transition to a three-dimensional flow in the side layers results from the advection or a precession and splitting of the Taylor–G?rtler vortex rings. The predictable behaviour of the Taylor–G?rtler vortices disappears with increasing magnetic field strength. The numerical simulations agree very well with the flow measurements.     

Influence of a homogeneous axial magnetic field on Taylor–Couette flow of ferrofluids with low particle–particle interaction

Experimental results concerning the stability of Couette flow of ferrofluids under magnetic field influence are presented. The fluid cell of the Taylor–Couette system is subject to a homogeneous axial magnetic field and the axial flow profiles are measured by ultrasound Doppler velocimetry. It has been found that an axial magnetic field stabilizes the Couette flow. This effect decreases with a rotating outer cylinder. Moreover, it could be observed that lower axial wave numbers are more stable at a higher axial magnetic field strength. Since the used ferrofluid shows a negligible particle–particle interaction, the observed effects are considered to be solely based on the hindrance of free particle rotation.     

Large-eddy simulation of counter-rotating Taylor–Couette flow: The effects of angular velocity and eccentricity

In the present work, turbulent flow in the annulus of a counter-rotating Taylor-Couette (CRTC) system is studied using large-eddy simulation. The numerical methodology employed is validated, for both the mean and second-order statistics, with the direct numerical simulation (DNS) data available in the literature, for a range of Reynolds numbers from 500 to 4000. Thereafter, turbulent flow occurring in this system at Reynolds numbers of 8000 and 16000 are studied, and the results obtained are analyzed using mean and second-order statistics, vortical structures, velocity vector plots and power energy spectra. Further, the spatio-temporal variation of azimuthal velocity, extracted near the inner cylinder, shows the existence of herringbone like patterns similar to that observed in the previous studies. The effect of eccentricity of the inner cylinder with respect to the outer cylinder is studied, on the turbulent flow in the CRTC system, for two different eccentricity ratios of 0.2 and 0.5 and for two different Reynolds numbers of 1500 and 4000. The results of the eccentric CRTC are analyzed using contours of pressure, mean and second-order statistics, velocity vectors, vortical structures, and turbulence anisotropy maps. It is observed from the eccentric CRTC simulations that the smaller-gap region seems to contain higher amplitude fluctuations and more vortical structures when compared with the larger-gap region. The mean turbulent kinetic energy contours do not change qualitatively with the Reynolds number, however, quantitatively a higher turbulent kinetic energy is observed in the higher Reynolds number case of 4000.     

Numerical study of Saffman–Taylor instability in immiscible nonlinear viscoelastic flows

In this paper, a numerical solution for Saffman–Taylor instability of immiscible nonlinear viscoelastic-Newtonian displacement in a Hele–Shaw cell is presented. Here, a nonlinear viscoelastic fluid pushes a Newtonian fluid and the volume of fluid method is applied to predict the formation of two phases. The Giesekus model is considered as the constitutive equation to describe the nonlinear viscoelastic behavior. The simulation is performed by a parallelized finite volume method (FVM) using second order in both the spatial and the temporal discretization. The effect of rheological properties and surface tension on the immiscible Saffman–Taylor instability are studied in detail. The destabilizing effect of shear-thinning behavior of nonlinear viscoelastic fluid on the instability is studied by changing the mobility factor of Giesekus model. Results indicate that the fluid elasticity and capillary number decrease the intensity of Saffman–Taylor instability.     

Effects of longitudinal grooves on the Couette–Poiseuille flow

The effects of longitudinal grooves on the flow resistance in a channel where the flow is driven by movement of one of the walls and modified by a streamwise pressure gradient have been studied. The reducedorder geometrymodel has been used to extract geometric features that are hydraulically relevant. Three distinct zones leading to the reduced resistance have been identified, depending on the flow pressure gradient and the groove wave number. Two of these zones correspond to grooves with long wavelengths and one to grooves with short wavelengths. Optimization has been used to determine shapes that provide the largest flow rate. In the case of the long-wavelength grooves, the optimal shapes depend on the constraints. These shapes are well approximated by a certain universal trapezoid for grooves that have the same height and depth. There exists an optimum depth which, combined with the corresponding shape, defines the optimal geometry in the case of the unequal-depth grooves; this shape is well approximated by a Gaussian function. No optimal shape exists for the short-wavelength grooves if the groove amplitude is sufficiently small; the shortest admissible wavelength dominates system performance under such conditions. The most effective groove wave number does exist for higher grooves, but the optimal shape cannot be determined due to numerical limitations.     

A note on the far-asymptotics of Helmholtz–Kirchhoff flows

In this sequel to a rather recent paper on the classical problem of Helmholtz–Kirchhoff flows by Vic. V. Sychev (TsAGI Sci J 41(5):531–533, 2010), the representation of the flow far from the body and its specific implications discussed in that study are revisited. Here the concise derivation of these findings resorts to well-known Levi-Cività’s method and, alternatively, only fundamental properties of analytic functions and thin-airfoil theory. As particularly of interest when the well-known Kirchhoff parabola degenerates to an infinitely long cusp, integration constants debated controversially so far and important for the understanding and computation of those flows are specified by the integral conservation of momentum. Also, the parametric modification towards flows encompassing stagnant-fluid regions of finite extent and the previously unnoticed impact of higher-order terms on the associated high-Reynolds-number flows are addressed.