Transition from vortex to wall driven turbulence production in the Taylor–Couette system with a rotating inner cylinder

Axisymmetrically stable turbulent Taylor vortices between two concentric cylinders are studied with respect to the transition from vortex to wall driven turbulent production. The outer cylinder is stationary and the inner cylinder rotates. A low Reynolds number turbulence model using the k‐ω formulation, facilitates an analysis of the velocity gradients in the Taylor–Couette flow. For a fixed inner radius, three radius ratios 0.734, 0.941 and 0.985 are employed to identify the Reynolds number range at which this transition occurs. At relatively low Reynolds numbers, turbulent production is shown to be dominated by the outflowing boundary of the Taylor vortex. As the Reynolds number increases, shear driven turbulence (due to the rotating cylinder) becomes the dominating factor. For relatively small gaps turbulent flow is shown to occur at Taylor numbers lower than previously reported. Copyright © 2002 John Wiley & Sons, Ltd.     

Discrete simulation and micromechanical analysis of two-dimensional saturated granular media

In this study, a novel approach to incorporate the pore water pressure in the discrete element method (DEM) to comprehensively model saturated granular media was developed. A numerical model was constructed based on the DEM by implanting additional routines in the basic DEM code; pore water pressure calculations were used with a two-dimensional (2D) model to simulate the undrained behavior of saturated granular media. This model coupled the interaction of solid particles and the pore fluid in saturated granular media. Finally, several 2D undrained shear tests were simulated. The test results showed that the model could predict the response of the saturated granular soil to shear loading. The effect of initial compaction was investigated. Biaxial tests on dense and loose specimens were conducted, and the effect of the initial density on the change in shear strength and the volume change of the system was investigated. The overall behavior of loose and dense specimens was phenomenologically similar to the real granular material. Constant volume tests were simulated, and the results were compared to those from the coupled model. Induced anisotropy was micromechanically investigated by studying the contact force orientation. The change in anisotropy depended on the modeling scheme. However, the overall responses of the media obtained using the coupled and constant volume methods were similar.     

Discrete simulation and micromechanical analysis of two-dimensional saturated granular media

In this study, a novel approach to incorporate the pore water pressure in the discrete element method （DEM） to comprehensively model saturated granular media was developed. A numerical model was constructed based on the DEM by implanting additional routines in the basic DEM code; pore water pressure calculations were used with a two-dimensional （2D） model to simulate the undrained behavior of satu- rated granular media. This model coupled the interaction of solid particles and the pore fluid in saturated granular media. Finally, several 2D undrained shear tests were simulated. The test results showed that the model could predict the response of the saturated granular soil to shear loading. The effect of initial compaction was investigated. Biaxial tests on dense and loose specimens were conducted, and the effect of the initial density on the change in shear strength and the volume change of the system was inves- tigated. The overall behavior of loose and dense specimens was phenomenologically similar to the real granular material. Constant volume tests were simulated, and the results were compared to those from the coupled model. Induced anisotropy was micromechanically investigated by studying the contact force orientation. The change in anisotropy depended on the modeling scheme. However, the overall responses of the media obtained usinz the couoled and constant volume methods were similar.     

Drag reduction using surfactants in a rotating cylinder geometry

Turbulent Couette flow between two circular cylinders has been used for drag reduction experiments using surfactants. In the experiments presented here, only the outer cylinder rotates, the inner cylinder remains at rest and accurate measurements of the torque at the inner cylinder are measured. Water is used as a reference fluid. A drag reducing surfactant called Arquad S-50 (Akzo Nobel Surface Chemistry LLC, Chicago, Ill., USA) (5 mM)+NaSal (12.5 mM) was used as the drag reduction agent. This surfactant can reduce the drag up to 70% (a Reynolds number of about 70,000–150,000) as measured by pressure drop in a pipe flow. Experiments in Couette flow also show drag reduction in the turbulent range. Two arrangements were used, (1) one small trip-wire on the inner cylinder, and (2) four larger trip-wires on the outer cylinder. These trips reduce the critical Reynolds number for transition from laminar to turbulent flow. In case (1), we obtained 18% drag reduction at 5,000<Re<15,000 and in case (2), we obtained an average reduction of about 20% at 2,000<Re<10,000, increasing up to 30% at Re=15,000. The paper also discusses two important problems. First, the shear rate is not constant in the radial gap in circular Couette flow. For non-Newtonian fluids, where the molecular viscosity is a function of the shear rate, this effect must be considered. Second, which viscosity should be used in the Reynolds number? For pipe flow measurements, most authors use the viscosity of the solvent (generally water and Newtonian). For measurements in the Couette flow, we use a different approach, which is described in this paper. We conclude that Couette flow is a useful method for drag reduction investigations. Its advantage is the much smaller geometry in comparison to those of conventional test facilities such as wind tunnels, water, or oil channels or in tubes.     

Description of local dilatancy and local rotation of granular assemblies by microstretch modeling

This study investigates the microstretch continuum modeling of granular assemblies while accounting for both the dilatant and rotational degrees of freedom of a macroelement. By introducing the solid volume fraction and the gyration radius of a granular system, the balance equations of the microstretch continuum are transformed into a new formulation of evolution equations comprising six variables: the solid volume fraction, the gyration radius, the velocity field, the averaged angular velocity, the rate of gyration radius, and the internal energy. The bulk microinertia density, the averaged angular velocity, and the microgyration tensor at a macroscopic point are obtained in terms of discrete physical quantities. The bulk part and the rotational part of the microgyration tensor are proposed as the two indices to measure the local dilatancy and local rotation of granular assemblies. It is demonstrated in the numerical simulation that the two indices can be used to identify the shear band evolution in a granular system under a biaxial compression.     

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.     

Integration approach of the Couette inverse problem of powder type self-compacting concrete in a wide-gap concentric cylinder rheometer

For powder type self-compacting concrete (SCC) mixes, commonly used in Belgium, a shear thickening (Herschel–Bulkley) flow behaviour of the fresh mixes is quite often observed.A longstanding problem in rheometry is the so-called “Couette inverse problem”, where one tries to derive the flow curve from the torque measurements T(N) in a (wide-gap) concentric cylinder (Couette) rheometer, with T the torque registered at the inner, stationary cylinder and N the rotational velocity of the outer, rotating, cylinder.In this paper, the Couette inverse problem is approached by means of the integration method in order to convert T(N) into for a wide-gap (R_o/R_i = 1.45) concentric cylinder rheometer. The approach consists in the decoupling of the flow resistance and the power-law flow behaviour after exceeding the flow resistance. The integration approach is validated by experimental verification with different powder type SCC mixtures. By means of illustration, the results of one limestone powder type SCC mixture with different superplasticizer contents are shown in this paper.     

散体介质SHPB被动围压试验体应力计算的理论修正方法

SHPB被动围压试验为探究散体介质在爆炸和冲击荷载作用下的力学行为提供了一个行之有效的方法。针对相关试验设计和计算中存在的弊端和不足,借助经典板壳理论将SHPB被动围压试验中用于约束散体介质的刚性套筒简化为受均匀带状内压作用的圆柱形壳体。理论计算了套筒径向位移、环向应变与均匀带状内压及套筒几何、力学参数的关系,得到了套筒径向位移、环向应变沿其轴向的分布规律;分析了套筒长度、厚度、内外径以及均匀带状内压宽度之间等无量纲几何参数对计算结果的影响;将理论计算结果与试验和数值模拟结果进行对比,验证了理论计算结果的正确性。本文中提出的理论修正方法可为指导散体介质SHPB被动围压试验提供参考。     

Stability of flow between two rotating cylinders in the presence of a constant heat flux at the outer cylinder and radial temperature gradient – wide gap problem

The stability of Couette flow of a viscous incompressible fluid between two concentric rotating cylinders in the presence of a radial temperature gradient due to a constant heat flux at the outer cylinder is studied. The critical values of `a' (the wave number) and Ta (the Taylor number) are listed in a table and some critical Taylor numbers are shown graphically. It is shown that as the heat flux is increased the flow becomes more unstable for all values of μ calculated, where μ is the ratio of the angular velocity of the outer cylinder to that of the inner cylinder. Received on 04 March 1997     

Normal stress effects in the gravity driven flow of granular materials

In this paper, we study the fully developed gravity-driven flow of granular materials between two inclined plates. We assume that the granular materials can be represented by a modified form of the second grade fluid where the viscosity depends on the shear rate and volume fraction and the normal stress coefficients depend on the volume fraction. We also propose a new isotropic (spherical) part of the stress tensor which can be related to the compactness of the (rigid) particles. This new term ensures that the rigid solid particles cannot be compacted beyond a point, namely when the volume fraction has reached the critical/maximum packing value. The numerical results indicate that the newly proposed stress tensor has obvious and physically meaningful effects on both the velocity and the volume fraction fields.     

Processing the Couette viscometry data using a Bingham approximation in shear rate calculation

This paper presents an approach to computing the shear flow curve from torque–rotational velocity data in a Couette rheometer. The approximation techniques in shear rate calculation are generally dictated by the radius ratio between coaxial cylinders and the rheological behaviour of fluid tested. Here, the approach consists in analysing the sheared material as a Bingham fluid and computing an average shear rate when the fluid in the cylindrical gap is partially and fully sheared. We focus in particular on the applicability of the Bingham approximation in shear rate calculation. First, the approach is assessed by examining synthetic data generated with Newtonian, non-Newtonian and yield stress materials with known properties, varying the gap radius ratio. The results, which are compared with commonly used techniques in shear rate calculation, prove the relevance of the proposed approach. Finally, its efficiency is examined by applying it to process Couette data of yield stress fluids taken from published works.     

Shear localisation in stirred yoghurt

The shear flow behaviour of stirred yoghurt in the cone-and-plate and cylindrical Couette geometries was studied using diffusing wave spectroscopy (DWS) and nuclear magnetic resonance (NMR) velocimetry. Differences between the transmission and backscattering DWS correlations suggest the formation of a high shear rate band near the surface of a moving cone of a cone-and-plate geometry at low shear rates. At higher shear rates, homogeneous shear flow is indicated. NMR velocimetry unambiguously demonstrated that a high shear rate band forms at the moving inner wall of a cylindrical Couette geometry at low shear rates. At intermediate shear rate, a high shear rate band is formed at the stationary outer wall and plug-like flow is observed mid-gap. At higher shear rates, homogeneous shear flow is observed. Slip is seen at both walls. The three flow regimes appear to correlate loosely with transitions in the pseudo-steady-state flow curve and may reflect a break-up of the protein aggregates observed with confocal microscopy.     

Migration of a sphere suspended in viscoelastic liquids in Couette flow: experiments and simulations

The intrinsically coupled effects of the curvature of the flow-field and of the viscoelastic nature of suspending medium on the cross-stream lateral migration of a single non-Brownian sphere in wide-gap Couette flow are studied. Quantitative videomicroscopy experiments using a counterrotating device are compared to the results of 3D finite element simulations. To evaluate the effects of differences in rheological properties of the suspending media, fluids have been selected which highlight specific constitutive features, including a reference Newtonian fluid, a single relaxation time wormlike micellar surfactant solution, a broad spectrum shear-thinning elastic polymer solution and a constant viscosity, highly elastic Boger fluid. As expected for conditions corresponding to Stokes flow, migration is absent in the Newtonian fluid. In the wormlike micellar solution and the shear-thinning polymer solution, spheres are observed to migrate in the direction of decreasing shear rate gradient, i.e. the outer cylinder, except when the sphere is initially released close to the inner cylinder, in which case the migration is towards it. The migration is enhanced by faster relative angular velocities of the cylinders. Shear-thinning reduces the migration velocity, showing an opposite behavior as compared to previous results in planar shear flow. In the Boger fluid, within experimental error no migration could be observed, likely due to the large solvent contribution to the overall viscosity. For small Deborah numbers the migration results are well described by an heuristic argument based on a local stress balance.     

Cessation of Newtonian circular and plane Couette flows with wall slip and non-zero slip yield stress

We solve analytically the cessation flows of a Newtonian fluid in circular and plane Couette geometries assuming that wall slip occurs provided that the wall shear stress exceeds a critical threshold, the slip yield stress. In steady-state, slip occurs only beyond a critical value of the angular velocity of the rotating inner cylinder in circular Couette flow or of the speed of the moving upper plate in plane Couette flow. Hence, in cessation, the classical no-slip solution holds if the corresponding wall speed is below the critical value. Otherwise, slip occurs only initially along both walls. Beyond a first critical time, slip along the fixed wall ceases, and beyond a second critical time slip ceases also along the initially moving wall. Beyond this second critical time no slip is observed and the decay of the velocity is faster. The velocity decays exponentially in all regimes and the decay is reduced with slip. The effects of slip and the slip yield stress are discussed.     

Numerical study of the generalized cylindrical Couette flow of rarefied gas

We study the cylindrical Couette flow of a rarefied gas between two cylinders in the generalized setup in which the inner of which not only rotates but also slides along its axis. The analysis is based on the numerical solution of the S-model kinetic equation. The influence of ratio of cylinder radiuses, velocities of the inner cylinder and Knudsen number on shear stresses, mass-flow rates as well as macroscopic parameters is investigated in the broad range of Knudsen numbers.     

Resolution sensitivity of momentum‐exchange and immersed boundary methods for solid–fluid interaction in the lattice Boltzmann method

In the lattice Boltzmann method (LBM), the mechanism of fluid–solid interaction can be effectively captured by appropriately enforcing the no‐slip conditions in shear direction, and bounce‐back of the non‐equilibrium distribution portion in the normal direction at fluid–solid interfaces. Among various solid–fluid interaction schemes being proposed for LBM in recent decades, two simple fluid–solid interaction methods—the momentum exchange algorithm (MEA) and the immersed boundary scheme (IBS)—were developed based on the above concept. In this paper, MEA and IBS are implemented in a D2Q9 LBGK system and applied to measure the wall correction factors of drag force upon a stationary circular particle midway in the Poiseuille channel flow at very low Reynolds number and drag coefficients at low to moderate Reynolds numbers. MEA and IBS are also employed to compare the fluid‐induced torque over the cylinder in the Taylor–Couette flow, and the steady velocity of a particle settling under the influence of gravity inside a tube. The above experiments show that IBS seems to be more accurate and less demanding on lattice resolution. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulation of Taylor Couette flow of Bingham fluids

The Bingham fluid flow between two concentric cylinders is studied using numerical simulation. The cylinders are assumed to rotate independently, and with an imposed axial sliding. The flow field is decomposed with linearity arguments of the base circular Couette shear flow and corresponding deviation field. The numerical methods are based on the expression of the deviation field in terms of complete sets of orthogonal functions and Chebyshev series. The Galerkin projection method is used with the pressure term being eliminated. The Adams Bashforth scheme is adopted for time marching. The results show that the vortices are squeezed toward the inner cylinder due to the effect of yield stress. When the outer cylinder is held stationary, the yield stress plays a role in weakening the vortex flow. However, for the co-rotation situation, the vortex flow is initially strengthened with an increase of yield stress, and then weakened as the yield stress is raised large enough. The annular unyielded regions emerge and stick to the outer cylinder. In case of Taylor Couette flow with an imposed axial sliding, a spiral vortex flow is visible with spiral unyielded region being obtained.     

A method for correction of the wall-slip effect in a Couette rheometer

A simple method for correction of the wall-slip effect in a Couette rheometer was derived. The method requires only two series of measurements (two flow curves) performed in two measuring sets of any dimensions, and therefore it may be applied for the results obtained in each rheometer with a standard cup and bob set. The method was checked for experimental data and also verified theoretically for a hypothetical liquid. H height of cylinder - M torque - r distance from axis - R _i ,R ₀ radius of inner and outer cylinder - R _m average radius defined by Eq. (7) - u slip velocity - shear rate - shear rate for no-slip conditions - Newtonian viscosity - angular speed - angular speed of the rotating cylinder