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     

Coriolis effects on the stability of pulsed flows in a Taylor–Couette geometry

Results steming from the linear stability of time-periodic flows in a Taylor–Couette geometry with cylinders oscillating in phase or out-of-phase are presented. Our analysis takes into account the gap size effects and investigates the influence of a superimposed mean angular rotation of the whole system.In case of no mean rotation, the finite gap geometry is found to affect the shape of the stability diagrams (critical Taylor number versus the frequency parameter) which consist of two distinct branches as opposed to being continuous in the narrow gap approximation. In particular, in the out-of-phase configuration a new branch for low frequencies was found, thus enabling better agreement with available experimental results.When cylinders are co-rotating and subject to rotation effects, our calculations provide the evolution of the critical Taylor number versus the rotation number for two values of the frequency. The stability curves are found to be in qualitative agreement with available experimental data revealing a maximum of instability for a rotation number of about 0.3.In the high rotation regime, enhancement of the critical Taylor number is investigated through an asymptotic analysis and the value of the rotation number at which restabilization occurs is found to depend on the frequency parameter.A restabilization of the flow also occurs when the rotation number and the gap size are of the same order, a phenomenon already pointed out in the case of steady flows and attributed to the near cancellation of Coriolis and centrifugal effects. Our investigation proves that the same mechanism still holds for time-periodic flows.     

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     

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.     

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.     

The effects of chain conformation in the microfluidic entry flow of polymer–surfactant systems

An associative polymer–surfactant system has been used to observe the effects of chain conformation in the entry flow through a microfabricated planar 16:1:16 contraction–expansion geometry. The well-studied system of the flexible polymer poly(ethylene oxide) (PEO) and anionic surfactant sodium dodecyl sulfate (SDS) was used. Dilute polymer solutions with increasing SDS concentration were characterized in steady and dynamic shear, as well as capillary breakup extensional rheology. Based on this characterization, the primary quantitative difference is an increase in zero-shear viscosity as a result of the PEO chain expansion brought on by association of SDS surfactant micelles. However, these quantitatively similar solutions were observed to exhibit much more qualitatively different flow patterns via fluorescent streak imaging in the entry flow. In contrast to previous work on PEO solutions, the PEO–SDS systems were observed to transition to a steady viscoelastic flow regime characterized by stable lip vortices at much lower elasticity and Weissenberg numbers. The resulting insight gained regarding the utility of microfluidic flows in elucidating effects of subtle conformational changes further illustrates the potential for using microfabricated devices as rheometric tools for measuring the properties of dilute and weakly viscoelastic fluids.     

Taylor–Couette flow with independently rotating end plates

Results are presented from a combined numerical and experimental study of steady bifurcation phenomena in a modified Taylor–Couette geometry where the end plates of the flow domain are allowed to rotate independently of the inner cylinder. The ends rotate synchronously and the ratio between the rate of rotation of the ends _e and the inner cylinder _i defines a control parameter :=_e/_i. Stationary ends favour inward motion along the end walls whereas rotating walls promote outward flow. We study the exchange between such states and focus on two-cell flows, which are found in the parameter range between =0 and =1 for =2. Hence is used as an unfolding parameter. A cusp bifurcation is uncovered as the organizing centre for the stability exchange between the two states. Symmetry breaking bifurcations, which lead to flows that break the mid-plane symmetry are also revealed. Overall, excellent agreement is found between numerical and experimental results. PACS 47.20, 47.11, 47.54     

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.     

Modelling the three-dimensional flow of a semi-dilute polymer solution in microfluidics—on the effect of aspect ratio

The flow of polymer solutions in microfluidic devices is inherently three-dimensional, especially in the non-linear flow regime, and often results in flow phenomena that might not even be encountered in macro-devices. Using a multi-mode Phan-Thien–Tanner model, three-dimensional (3-D) simulations of a semi-dilute polyethylene oxide (PEO) solution through 8:1 planar contraction micro-channels with various depths have been carried out to systematically study the effect of the aspect ratio on the flow fields. Vortex dynamics in the upstream flow section and excess pressure drop are quantified in detail. A transition from a salient-corner vortex mechanism to a lip vortex mechanism is observed as the aspect ratio is varied from 1 to 1/4, which corresponds to the elasticity numbers El = 36.3 to 48.3. The numerical results show that varying the aspect ratio of microfluidic channels has similar effects to varying other parameters, such as fluid properties, which influence the elasticity number. Thus, our results support the view that vortex growth mechanism is determined by the elasticity number, which is fixed for a given fluid and geometry. The principle is of significance to the design of new microfluidic chips for a wide range of applications.     

An approximate solution for the Couette–Poiseuille flow of the Giesekus model between parallel plates

An approximate analytical solution is derived for the Couette–Poiseuille flow of a nonlinear viscoelastic fluid obeying the Giesekus constitutive equation between parallel plates for the case where the upper plate moves at constant velocity, and the lower one is at rest. Validity of this approximation is examined by comparison to the exact solution during a parametric study. The influence of Deborah number (De) and Giesekus model parameter (α) on the velocity profile, normal stress, and friction factor are investigated. Results show strong effects of viscoelastic parameters on velocity profile and normal stress. In addition, five velocity profile types were obtained for different values of α, De, and the dimensionless pressure gradient (G).     

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.     

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.     

Channel, tube, and Taylor–Couette flow of complex viscoelastic fluid models

We show how to formulate two-point boundary value problems to compute laminar channel, tube, and Taylor–Couette flow profiles for some complex viscoelastic fluid models of differential type. The models examined herein are the Pom-Pom Model [McLeish and Larson 42:81–110, (1998)] the Pompon Model [Öttinger 40:317–321, (2001)] and the Two Coupled Maxwell Modes Model (Beris and Edwards 1994). For the two-mode Upper-Convected Maxwell Model, we calculate analytical solutions for the three flow geometries and use the solutions to validate the numerical methodology. We illustrate how to calculate the velocity, pressure, conformation tensor, backbone orientation tensor, backbone stretch, and extra stress profiles for various models. For the Pom-Pom Model, we find that the two-point boundary value problem is numerically unstable, which is due to the aphysical non-monotonic shear stress vs shear rate prediction of the model. For the other two models, we compute laminar flow profiles over a wide range of pressure drops and inner cylinder velocities. The volumetric flow rate and the nonlinear viscoelastic material properties on the boundaries of the flow geometries are determined as functions of the applied pressure drop, allowing easy analysis of experimentally measurable quantities.     

Direct numerical simulation of Taylor–Couette turbulent flow controlled by a traveling wave-like blowing and suction

In wall turbulence, a traveling wave-like control is known to decrease the skin-friction drag and induce the relaminarization phenomenon. Because it is noteworthy to investigate the control effect in other canonical flows, direct numerical simulations of fully developed turbulent Taylor–Couette flows are performed. The Reynolds number, based on the wall velocity of a rotating inner cylinder and the radius of a centerline between cylinders, is set to 84,000. The traveling wave-like blowing and suction is imposed on the inner or outer cylinder wall, and the control effect is parametrically investigated. In the inner cylinder control, the torque reduction is obtained when the wave travels in the co-rotating direction with the inner cylinder, and its wavespeed is faster than the rotation. In the outer cylinder control, in contrast, the torque reduction is obtained when the wave propagates in the opposite direction. While the control is imposed on one side wall (i.e., inner or outer cylinder), the control affects the entire flow region. The Taylor vortex remains, while the traveling wave affects its strength. The three-component decomposition analysis shows that the traveling wave creates the coherent contribution on the torque, while the random contribution on it is reduced. Accordingly, a major factor of the torque reduction in the Taylor–Couette flow is the reduction of the random contribution. In addition, for the faster wavespeed cases with the small wavenumber (i.e., the long wavelength), the drag reduction larger than 60% is obtained and the relaminarization occurs in these cases.     

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.     

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.