Maurice Couette,one of the founders of rheology

This article presents a biography of Maurice Couette, whose name is associated with a type of flow, of viscometer, and with a correction method for end effects in capillary flows. His life and work are described, with special mention being made of the cylinder apparatus that he designed. The relevance of his work to present day rheology is stressed.     

An analytical solution for the velocity and shear rate distribution of non-ideal Bingham fluids in concentric cylinder viscometers

An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/_y - c = 2µ - d = ₀ ²–2c_y - f() = (–₀)²+2c(–_y) - r Radius - r _b Bob radius - r _c Cup radius - r _y Yield radius - r ₀ Stationary surface radius - r Rotating Stationary radius - Y ₀ Shear rate parameter = /µ Greek letters Shear rate - = (r _y /r _b )²– 1 - µ Bingham viscosity - µ₀ Initial differential viscosity - µ µ₀-µ - Rotation speed - Angular velocity - Shear stress - _b Bob shear stress - _B Bingham stress - _y (True) yield stress - ₀ Stress parameter = _B +µ Y ₀ - _B - _y      

Obtaining the shear stress shear rate relationship and yield stress of liquid foods from Couette viscometry data

A method based on Tikhonov regularisation is used to convert the Couette viscometry data of a number of liquid foods into shear stress vs shear rate relationships. For liquid foods that have a yield stress and if the viscometry data cover the appropriate range of shear stress, Tikhonov regularisation is used to compute simultaneously the yield stress. The versatility of Tikhonov regularisation is demonstrated by applying it to process data from wide as well as narrow gap Couette viscometers. The results generated by Tikhonov regularisation are compared against that obtained by methods currently employed by rheologists.     

Couette flow of a Bingham plastic in a channel with equally porous parallel walls

In the present paper the flow of a Bingham fluid between two parallel porous walls is studied. One of the walls moves with constant velocity parallel to the other, which is fixed, while a longitudinal pressure gradient exists, as well as a transverse flow field due the porosity of the walls. An exact analytical solution is given for the u-velocity field, which has four different forms depending on the values of the three dimensionless parameters, which are the Bingham, Couette and Reynolds numbers.     

Start-up Flow of a Bingham Fluid in a Pipe

We present an analytical solution of axisymmetric motion for a Bingham fluid initially at rest subjected to a constant pressure gradient applied suddenly. Using the Laplace transform, we obtain expressions which allow the calculation of the instantaneous velocity, plug radius and rate of flow as a function of time. We also give a relation for the shear stress in the plug and in the region where the behaviour of the fluid is Newtonian.     

Rheometric and rheo-optical investigation on the effect of the aliphatic chain length of the surfactant on the shear thickening of dilute worm-like micellar solutions

In this study, a new way of understanding the shear-thickening phenomenon in self-assembled solutions is introduced. The near- and out-of-equilibrium behavior is investigated in four aqueous micellar solutions containing surfactants of the same family: the alkyltrimethylammonium bromide mixed with an ionic salt at equimolar concentration. Thus, the four molecules of surfactants have the same polar head but different aliphatic chain length containing 12, 14, 16, and 18 carbon atoms. In aqueous solutions, the attractive forces between the surfactant molecules depend on the length of the aliphatic chain, thus, varying this parameter will have a definite influence on the aggregation number, the shape and dimension of the micelles. According to our results, the evolution of the low shear viscosity is affected by the chain length. The rheometric measurements performed in steady and time-dependent flows show that the emergence, the range, and the amplitude of the shear thickening also depend on this parameter. However, in the domain of high shear rates, after reaching the maximum of viscosity, all flow curves superimpose irrespective of the chain length. The rheo-optic measurements confirm the apparition of the shear-induced structure (SIS); this new phase is locally oriented and the chain length affects strongly the micelles orientation and the birefringence intensity. These results undoubtedly demonstrate that the chain length plays an important role in the behavior near equilibrium and under shear flow of the micellar systems. Paper presented at the third Annual European Rheology Conference (AERC 2006) on April 27–29, 2006, Crete, Greece.     

An implicit solution of bi-penalty approximation with orthogonality projection for the numerical simulation of Bingham fluid flow

1.IntroductionItisamajordiffct.encefi-omtheNewtonnuidflowthattheBinghammodelofNonNewtonfluidflowischaracterizedbytwoparameters:ayieldstressandaviscosity.WhenthestressoftheBinghalnmaterialbelowtheyieldstress,materialisrigidotherwisethequasiNewtolliannowresultstll:'71.Hence,therearesomeofthefloating"rigidcores"involvedintheBinghamfluidfloworsomeofthe'rigidcores"attachedtotheboundaries,inwhichthelocationsalldshapesofthese"rigidcores"maychangeforthetransientBinghamfluid,flow.ThisBingllammodelh…     

Peristaltic flow of a Bingham fluid in a channel

We study the peristaltic transport of a Bingham fluid in a channel with small aspect ratio whose walls behave as a periodic traveling wave. The governing equations in the unyielded phase are obtained writing the integral formulation for the momentum balance. As shown in Fusi et al. (2015), this approach allows to overcome the so-called “lubrication paradox” which may arise in the thin film approximation. We consider the case in which the inlet flux is prescribed and the one in which the flow is driven by a given pressure drop. In both cases the solution of the problem is determined solving a nonlinear integral equation for the yield surface. We perform some numerical simulations to illustrate the behavior of the yield surface, assuming that the traveling wave describing the peristaltic motion has a sinusoidal shape.     

Spots and turbulent domains in a model of transitional plane Couette flow

A review of the globally subcritical transition to turbulence in shear flows is presented, with an emphasis on the cases of plane and circular Couette flows (pCf and cCf, respectively). A Swift–Hohenberg-like model is next proposed to interpret the behavior of plane Couette flow in the vicinity of its global stability threshold. We present results of numerical simulations supporting this proposal and helping us to raise good questions about the growth and decay of intermittent turbulent domains in this precise context, and more generally about the coexistence of laminar flow and turbulence in other spatio-temporally intermittent flows. PACS 47.27.-i, 47.54.-r, 05.45.-a     

Flow development in compression of a finite amount of a Bingham plastic

The flow and shape evolution during the compression of a finite amount of a Bingham plastic is investigated by means of numerical simulations. The problem relates to the popular compression test used for the rheological characterization of non-Newtonian fluids. The flow is modelled in Lagrangian coordinates using the Papanastasiou regularization for the Bingham plastic and a mixed-Galerkin finite element method. Simulations have been performed for compression under both constant load and constant velocity. Results for various Reynolds and Bingham numbers are presented and discussed.     

Dispersion in a Bingham plastic fluid: Effects of peripheral layer

The dispersion of a soluble matter in a plastic fluid flowing through a tube and a channel has been analysed by taking into account the variations of viscosity, diffusivity and yield stress. It has been shown that in the special case of a Bingham fluid, surrounded by a peripheral layer of a Newtonian fluid, the effective dispersion coefficient with which the solute disperses across a plane moving with the mean speed of the flow decreases with the viscosity of the peripheral layer fluid but increases as the molecular diffusion coefficient of this layer decreases. Further, the effective dispersion coefficient also decreases as the yield stress of the Bingham fluid increases.     

Plane Couette flow of binary gaseous mixture in the whole range of the Knudsen number

The Couette flow of binary gaseous mixtures is studied on the basis of the McCormack model of the Boltzmann equation, which was solved numerically by the discrete velocity method. The calculations were carried out for three mixtures of noble gases: neon–argon, helium–argon, and helium–xenon. The stress tensor and bulk velocity of both species were calculated for several values of the gas rarefaction in the range from 0.01 to 40 for three values of the molar concentrations: 0.1,0.5 and 0.9. The numerical solution together with an analytical solution based on the slip boundary condition cover the whole range of the gas rarefaction. It was showed that the Couette flow is weakly affected by the intermolecular interaction law.     

On the modeling of fluidity loss phenomena in Couette and Poiseuille flows of elastic liquids

Some effects of the possible relaxation transition from viscoelastic liquid state to highly elastic solid state were theoretically and numerically investigated in the shear situations, within the approach proposed in papers [1, 2, 5, 16]. It was found that for a single Maxwellian model the constitutive equations developed in [1, 2, 5] are not valid at elevated shear stresses. Some new aspects of the possible rheological behavior of elastic liquids in subcritical (before transition) and supercritical (after transition) regimes were demonstrated. The mechanism of fluidity loss studied in this paper could serve as a possible trigger mechanism for the melt flow instabilities.     

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.     

From global to local bifurcations in a forced Taylor–Couette flow

The unfolding due to imperfections of a gluing bifurcation occurring in a periodically forced Taylor–Couette system is analyzed numerically. In the absence of imperfections, a temporal glide-reflection Z₂ symmetry exists, and two global bifurcations occur within a small region of parameter space: a heteroclinic bifurcation between two saddle two-tori and a gluing bifurcation of three-tori. As the imperfection parameter increase, these two global bifurcations collide, and all the global bifurcations become local (fold and Hopf bifurcations). This severely restricts the range of validity of the theoretical picture in the neighborhood of the gluing bifurcation considered, and has significant implications for the interpretation of experimental results. PACS 47.20.Ky, 47.20.Lz, 47.20.Ft     

Sensitivity of shear rate in artificial grafts using automatic differentiation

An accurate numerical simulation of blood requires the solution of incompressible Navier–Stokes equations coupled with specific constitutive models. We consider a generalized Newtonian fluid model in which viscosity depends on shear rate, accounting for the shear‐thinning behavior of blood. Previous work on the design of an artificial graft indicated that there is an influence of the fluid model on the solution of the partial differential equation‐constrained shape optimization problem. Therefore, we carry out a sensitivity analysis of the actual implementation of the flow solver using automatic differentiation (AD). We compare the sensitivities of shear rate with respect to viscosity for different configurations and validate the truncation‐error‐free sensitivities obtained from AD with those based on divided differencing and, if available, with analytic derivatives. Copyright © 2009 John Wiley & Sons, Ltd.     

On the stability of the shear flow of a viscoelastic fluid with slip along the fixed wall

In this paper we solve the time-dependent shear flow of an Oldroyd-B fluid with slip along the fixed wall. We use a non-linear slip model relating the shear stress to the velocity at the wall and exhibiting a maximum and a minimum. We assume that the material parameters in the slip equation are such that multiple steady-state solutions do not exist. The stability of the steady-state solutions is investigated by means of a one-dimensional linear stability analysis and by numerical calculations. The instability regimes are always within or coincide with the negative-slope regime of the slip equation. As expected, the numerical results show that the instability regimes are much broader than those predicted by the linear stability analysis. Under our assumptions for the slip equation, the Newtonian solutions are stable everywhere. The interval of instability grows as one moves from the Newtonian to the upper-convected Maxwell model. Perturbing an unstable steady-state solution leads to periodic solutions. The amplitude and the period of the oscillations increase with elasticity.     

缓凝砂浆在单螺杆挤出机中的流动分析

在预应力钢筋上连续包裹一层缓凝砂浆是新型缓粘结预应力混凝土体系得以工程应用的关键。为了弄清包覆过程,了解物料性能、加工条件和螺杆参数等输入变量对包覆输出结果的影响,基于缓凝预应力筋包覆试验,运用Polyflow软件,对宾汉姆体的缓凝砂浆在螺旋槽内流动和变形情况进行CFD(computational fluid dynamics)模拟。得到了可视化的速度场和剪切速率场,模拟了分析螺杆转速、螺杆几何参数对物料流变特性的影响。模拟结果可视,且直观地验证了包覆试验现象,揭示了包覆机理,为包覆装置可靠作业和结构优化提供了理论依据。