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.     

A procedure to determine the viscosity function from experimental data of capillary flow

The accurate calculation of the viscosity η as function of the shear rate &γdot; from capillary viscometry is still a matter of debate in the literature. In fact, this problem involves the inversion of an integral equation, which leads to multiple solutions due to the unavoidable noise present in the experimental data. The purpose of this work is to develop an efficient procedure to determine the viscosity function from experimental data of capillary flow without presenting the difficulties inherent in other methods discussed previously in the literature. The system identification procedure is used here to estimate the parameters of a viscosity model, which is appropriately selected for the fluid under study through preliminary calculations involving the apparent shear rate – shear stress data. Once the model is chosen by satisfying criteria for the fit goodness and its parameters are evaluated, a smooth and continuous function η(γdot;) is obtained in the range of experimental shear rates. The procedure proposed is also applicable to fluids in shear flow that present two Newtonian plateaus, as it is typically found in macromolecular dilute solutions. The mean value theorem of continuous functions is used to reduce significantly the computational time. Received: 15 November 1999 Accepted: 7 November 2000     

Obtaining the steady shear rheological properties and apparent wall slip velocity data of a water-in-oil emulsion from gap-dependent parallel plate viscometry data

A two-stage Tikhonov regularisation procedure has been used to obtain rheological properties for a high internal phase emulsion from gap-dependent steady-state parallel plate shear data. This method is beneficial in that it can convert the steady shear data into rheological property functions. The built-in regularisation parameters of the method are able to keep noise amplification under control. The two-stage method is able to obtain not only the shear stress–shear rate function but also the apparent slip velocity as a function of wall shear stress. The method is such that it obtains the rheological functions over the maximum range of shear rate covered by the data. The results obtained using the new method are compared to those obtained using the vane geometry with good agreement being observed.     

A general method for obtaining shear stress and normal stress functions from parallel disk rheometry data

The problems of converting the torque and normal force versus rim shear rate data generated by parallel disk rheometers into shear stress and normal stress difference as functions of shear rate are formulated as two independent integral equations of the first kind. Tikhonov regularization is used to obtain approximate solutions of these equations. This way of handling parallel disk rheometer data has the advantage that it is independent of the rheological constitutive equation and noise amplification is kept under control by the user-specified parameter in Tikhonov regularization. If the fluid under test exhibits a yield stress, Tikhonov regularization computation will simultaneously give an estimate of the yield stress. The performance of this method is demonstrated by applying it to a number of data sets taken from the published literature and to laboratory measurements conducted specifically for this investigation.     

Anomalous predictions of pressure drop and heat transfer in ducts of arbitrary cross-section with modified power law fluids

 The apparent viscosities of purely viscous non-Newtonian fluids are shear rate dependent. At low shear rates, many of such fluids exhibit Newtonian behaviour while at higher shear rates non-Newtonian, power law characteristics exist. Between these two ranges, the fluid's viscous properties are neither Newtonian or power law. Utilizing an apparent viscosity constitutive equation called the “Modified Power Law” which accounts for the above behavior, solutions have been obtained for forced convection flows. A shear rate similarity parameter is identified which specifies both the shear rate range for a given fluid and set of operating conditions and the appropriate solution for that range. The results of numerical solutions for the friction factor–Reynolds number product and for the Nusselt number as a function of a dimensionless shear rate parameter have been presented for forced fully developed laminer duct flows of different cross-sections with modified power law fluids. Experimental data is also presented showing the suitability of the “Modified Power Law” constitutive equation to represent the apparent viscosity of various polymer solutions. Received on 21 August 2000     

Elongational and shear rheology of carbon nanotube suspensions

Rheological behavior of concentrated suspensions of chemical vapor deposition carbon nanotubes in uniaxial elongation and simple shear is studied experimentally and theoretically. Nanotubes are suspended in viscous host liquids—castor oil or its blends with n-decane. The elongational measurements are performed by analyzing self-thinning (due to surface tension effect) liquid threads of nanotube suspensions. A quasi-one-dimensional model is used to describe the self-thinning process, whereas corrections accounting for thread nonuniformity and necking are introduced a posteriori. The effects of nanotube concentration and aspect ratio, viscosity of the suspending liquid, and initial diameter of the self-thinning thread in uniaxial elongation are elucidated. The results for uniaxial elongation are compared with those for simple shear. The correspondence in the results of the shear and elongational measurements is addressed and interpreted. The results conform to the Herschel–Bulkley rheological constitutive equation (i.e., power law fluids with yield stress). However, the yield stress in elongation is about 40% higher than in simple shear flow, which suggests that the original Herschel–Bulkley model need modification with the yield stress being a function of the second invariant of the deviatoric stress tensor. The present effort is the first to study capillary self-thinning of Herschel–Bulkley liquids, which are exemplified here by suspensions of carbon nanotubes.     

Numerical breakdown at high Weissenberg number in non-Newtonian contraction flows

 Planar contraction flows of non-Newtonian fluids with integral constitutive models are studied to investigate the problem of numerical breakdown at high Weissenberg or Debrorah numbers. Spurious shear stress extrema are found on the wall downstream of the re-entrant corner for both sharp and rounded corners. Moreover, a non-monotonic relation between shear stress and strain rate is found when the Deborah number limit is approached, which correlates with these shear extrema. This strongly suggests that non-monotonicity between shear stress and strain rate may be responsible for the Deborah number limit problem in contraction flow simulations. This non-monotonicity is caused by the inaccuracy of the quadrature, using constitutive equations that do not have shear stress maxima when exactly evaluated. This conclusion agrees with recent analytical findings by others that inaccuracy of the integration along the streamlines – either by numerical integration or asymptotic approximation – makes the problem ill-conditioned, with spurious growth occurring on the wall downstream of the re-entrant corner. Received: 5 March 1999/Accepted: 1 September 1999     

Nonlinear oscillating-cup viscometry

The general algorithm of calculations for oscillating-cup viscometry of fluids with viscous, elastic and plastic components is developed for damped and forced modes. The working equations efficient for practical applications are presented and related problems of data processing are analyzed. The case of forced oscillations of viscometer filled with rheostable fluids is treated in detail.     

Stability of a liquid jet into incompressible gases and liquids: Part 2. Effects of the irrotational viscous pressure

In this paper we investigate the effects of an irrotational, viscous pressure on the stability of a liquid jet into gases and liquids. The analysis extends our earlier work (part 1) in which the stability of the viscous jet was studied assuming that the motion and pressure are irrotational and the viscosity enters through the jump in the viscous normal stress in the normal stress balance at the interface. The liquid jet is always unstable; at high Weber numbers the instability is dominated by capillary instability; at low W the instability is dominated by Kelvin–Helmholtz (KH) waves generated by pressures driven by the discontinuous velocity. In the irrotational analysis the viscosity is important but the effects of shear are neglected. In fact a discontinuous velocity is not compatible with the continuity of the tangential components of velocity and shear stress so that KH instability is not properly posed for exact study using the no-slip condition but some of the effects of viscosity can be ascertained using viscous potential flow. The theory is called viscous potential flow (VPF). Here we develop another irrotational theory in which the discontinuities in the irrotational tangential velocity and shear stress are eliminated in the global energy balance by selecting viscous contributions to the irrotational pressure. These pressures generate a hierarchy of potential flows in powers of the viscosity, but only the first one, linear in viscosity, in the irrotational viscous stress, is thought to have physical significance. The tangential velocity and shear stress in an irrotational study cannot be made continuous, but the effects of the discontinuous velocity and stress in the mechanical energy balance can be removed “in the mean.” This theory with the additional viscous pressure is called VCVPF, viscous correction of VPF. VCVPF is VPF with the additional pressures. The theory here cannot be compared with an exact solution, which would not allow the discontinuous velocity and stress. In other problems, like capillary instability, in which VCVPF can be compared with an exact solution, the agreements are uniformly excellent in the wave number when one of the fluids is gas and in good but not uniform, agreement when both fluids are liquids.     

High resolution COD image analysis for health monitoring of reinforced concrete structures through inverse analysis

Direct and inverse problems of a fracture mechanics based RC beam model are solved. Solution of the direct problem that maps crack bridging stresses into crack opening displacements (COD) is straightforward, but the inverse problem is ill-posed, and better solved by the theory of inverse problems. This paper exploits the Tikhonov regularization method to solve the inverse problem, and estimates the force and location of rebar in buried concrete from CODs. Bending tests are carried out on model RC beams in the laboratory to demonstrate the applicability of the method. During the tests, a microscopic camera snaps high resolution digital pictures of cracked concrete surface. The images are analyzed by a software to measure surface CODs that are input into the inverse problem. The practical CODs inevitably include noise due to experimental error, which makes the inverse problem ill-posed, and necessitates regularization. In the current inverse analysis by the Tikhonov regularization method, bridging stress profiles, i.e. variation of the crack bridging stress along the crack length, has been figured out. Results are compared with those from other theoretical methods of analysis as well as with the readings from strain gauges. The method is a suitable non-destructive means for existing structures in cases where the section information is inadequate, or damages/repairs have altered the designed cross-section.     

Measuring static yield stress of electrorheological fluids using the slotted plate device

An electrorheological (ER) response is defined as the dramatic change in rheological properties of a suspension of small particles due to the application of a large electric field transverse to the direction of flow. ER fluids are typically composed of nonconducting or semiconducting particles dispersed in a nonconducting continuous phase. A sufficiently large electric field will cause ER fluids to solidify, giving rising to a yield stress. Many applications in torque and stress transfer devices were proposed employing the reversible yielding behavior of ER fluids. Successful applications depend on a large yield stress of ER fluids and therefore accurate measurements of the yield stress of ER fluids are required. Reported experimental yield stresses of ER fluids have been dynamic yield stresses obtained by extrapolating the shear stress–shear rate data to zero-shear rate. It would be very helpful to the understanding of ER behaviors and the applications of ER fluids to be able to measure the static yield stress of ER fluids accurately. The slotted plate technique has been shown to be a successful method to determine the static yield stress of suspensions. The values obtained via the slotted plate method are static yield stress as the platform is designed for extremely low-speed motion. In this study, we modified the slotted plate device for the application of large electric fields and measured the static yield stress of TiO₂ ER fluids under various electric fields. The measured static yield stress values are also compared with the static yield stress values from a commercial rheometer.     

Characterization of yield stress fluids with concentric cylinder viscometers

The methods normally employed for shear rate calculations from concentric cylinder viscometer data generally are not applicable for fluids with a yield stress. In cylindrical systems with large radius ratios, as usually is the case with suspensions, the yield stress induces two possible flow regimes in the annulus. Unless the yield value is exceeded everywhere in the gap only part of the fluid can be sheared while the remaining region behaves like a solid plug. A correct calculation of the shear rate must take into account the presence of a variable effective gap width determined by the extent of the sheared layer. For time-independent yield stress fluids, a two-step procedure, which does not require any specific flow model, is proposed for analysing the experimental torque-speed data. Under the partially sheared condition, the shear rate can be computed exactly, whereas for the fully sheared flow the Krieger and Elrod approximation is satisfactory. The method is assessed by examining both semi-ideal data generated with a Casson fluid with known properties, and experimental data with an industrial suspension. A more complicated problem associated with characterization of time-dependent yield stress fluids is also identified and discussed. An approximate procedure is used to illustrate the dependence of the shear rate on time of shear in constant-speed experiments.     

Wall slip of mayonnaises in viscometers

The shear flow of mayonnaise is generally characterized by an apparent yield stress, shear thinning in steady flow, stress overshoots upon inception of flow and other time-dependent effects. These observations are usually understood to be the result of structural rearrangement within the material. Additionally and separately, the possibility that emulsions may exhibit apparent wall slip on a microscopic scale at a solid-liquid boundary has been reported by some researchers. Thus, observed rheological behavior is likely to be the result of the interplay between these two phenomena. In the present work, it is demonstrated that when measurements are sought to be made on mayonnaise using rotational viscometers visible wall slip occurs, rendering such instruments ineffective for the purpose of making viscosity measurements even at shear rates as low as 10^–3s^–1. The factors that influence the onset and extent of slip are investigated with the help of parallel plate viscometers, and it is concluded that the observed “yielding” of mayonnaise is actually an artifact of the onset of macroscopic slip. Slip effects are also found in capillary flow but are ameliorated with increasing shear rate. To circumvent these problems, it is proposed that extensional viscometry be employed for determining the flow behavior of mayonnaises. Received: 18 August 1997 Accepted: 1 April 1998     

The motion of long bubbles through viscoelastic fluids in capillary tubes

The penetration of long gas bubble through a viscoelastic fluid in a capillary tube has been studied in order to investigate the influence of viscoelastic material properties on the hydrodynamic coating thickness and local flow kinematics. Experiments are conducted for three tailored ideal elastic (Boger) fluids, designed to exhibit similar steady shear properties but substantially different elastic material functions. This allows for the isolation of elastic and extensional material effects on the bubble penetration process. The shear and extensional rheology of the fluid is characterized using rotational and filament stretching rheometers (FSR). The fluids are designed such that the steady-state extensional viscosity measured by the FSR at a Deborah number (De) greater than 1 differs over three orders of magnitude (Trouton ratio = 10³–10⁶). The experiment set up to measure the hydrodynamic coating thickness is designed to provide accurate data over a wide range of capillary numbers (0.01 < Ca < 100). The results indicate that the coating thickness in this process increases with an increase in the extensionally thickening nature of the fluid. Experiments are also conducted using several different capillary tube diameters (0.1 < D < 1 cm), in order to compare responses at similar Ca but different flow De. Suitable scaling methods and nonlinear viscoelastic constitutive equations are explored to characterize the displacement process for polymeric fluids. Bubble tip shapes at different De are recorded using a CCD camera, and measured using an edge detection algorithm. The influence of the mixed flow field on the bubble tip shape is examined. Particle tracking velocimetry experiments are conducted to compare the influence of viscoelastic properties on the velocity field in the vicinity of the bubble tip. Local shear and extension rates are calculated in the vicinity of the bubble tip from the velocity data. The results provide quantitative information on the influence of elastic and extensional properties on the bubble penetration process in gas-assisted injection molding. The bubble shape and velocity field information provides a basis for evaluating the performance of constitutive equations in mixed flow. Received: 19 January 1999 Accepted: 30 June 1999     

Tikhonovs regularization method for ill-posed problems

Frequently the determination of material characteristic functions, such as the molecular mass distribution of a polymeric sample or the relaxation spectrum of a viscoelastic fluid, leads to an ill-posed problem. When Tikhonov regularization is applied to such a problem the problem of an appropriate choice of the regularization parameter arises. Well-known methods to determine this parameter, such as the discrepancy principle, and a method based on the minimization of the predictive mean-square signal error are compared with a self-consistence method. Monte Carlo simulations have been carried out for the determination of the relaxation spectrum from small amplitude oscillatory shear flow data. The self-consistence method has proven to be much more robust and reliable.     

Super-stable parallel flows of multiple visco-plastic fluids

We consider the stability of a multi-layer plane Poiseuille flow of two Bingham fluids. It is shown that this two-fluid flow is frequently more stable than the equivalent flow of either fluid alone. This phenomenon of super-stability results only when the yield stress of the fluid next to the channel wall is larger than that of the fluid in the centre of the channel, which need not have a yield stress. Our result is in direct contrast to the stability of analogous flows of purely viscous generalised Newtonian fluids, for which short wavelength interfacial instabilities can be found at relatively low Reynolds numbers. The results imply the existence of parameter regimes where visco-plastic lubrication is possible, permitting transport of an inelastic generalised Newtonian fluid in the centre of a channel, lubricated at the walls by a visco-plastic fluid, travelling in a stable laminar flow at higher flow rates than would be possible for the single fluid alone.     

Yield stress of structured fluids measured by squeeze flow

Various structured fluids were placed between the parallel circular plates of a squeeze-flow rheometer and squeezed by a force F until the fluid thickness h was stationary. Fluid thickness down to a few microns could be measured. Most fluids showed two kinds of dependence of f on h according to an experimentally-determined thickness h ^*. If h > h ^* then F varied in proportion to h ⁻¹ as predicted by Scott (1931) for a fluid with a shear yield stress τ₀. The magnitude of τ₀ from squeeze-flow data in this region was compared with the yield stress measured by the vane method. For some fluids τ₀ measured by squeeze flow was less than the vane yield stress, suggesting that the yield stress of fluid in contact with the plates was less than the bulk yield stress. If h < h ^* then F varied approximately as h ^−5/2 and the squeeze-flow data in this region analysed with Scott's relationship gave a yield stress which increased as the fluid thickness decreased. This previously unreported effect may result from unconnected regions of large yield stress in the fluid of size similar to h ^* which are not sensed by the vane and which become effective in squeeze flow only when h < h ^*. Received: 13 December 1999/Accepted: 4 January 2000     

Large amplitude oscillatory shear of pseudoplastic and elastoviscoplastic materials

We explore the utility of strain-controlled large amplitude oscillatory shear (LAOS) deformation for identifying and characterizing apparent yield stress responses in elastoviscoplastic materials. Our approach emphasizes the visual representation of the LAOS stress response within the framework of Lissajous curves with strain, strain rate, and stress as the coordinate axes, in conjunction with quantitative analysis of the corresponding limit cycle behavior. This approach enables us to explore how the material properties characterizing the yielding response depend on both strain amplitude and frequency of deformation. Canonical constitutive models (including the purely viscous Carreau model and the elastic Bingham model) are used to illustrate the characteristic features of pseudoplastic and elastoplastic material responses under large amplitude oscillatory shear. A new parameter, the perfect plastic dissipation ratio, is introduced for uniquely identifying plastic behavior. Experimental results are presented for two complex fluids, a pseudoplastic shear-thinning xanthan gum solution and an elastoviscoplastic invert-emulsion drilling fluid. The LAOS test protocols and the associated material measures provide a rheological fingerprint of the yielding behavior of a complex fluid that can be compactly represented within the domain of a Pipkin diagram defined by the amplitude and timescale of deformation.     

Rheo-PIV of a yield-stress fluid in a capillary with slip at the wall

An analysis of the yielding and flow behavior of a model yield-stress fluid, 0.2 wt% Carbopol gel, in a capillary with slip at the wall has been carried out in the present work. For this, a study of the flow kinematics in a capillary rheometer was performed with a two-dimensional particle image velocimetry (PIV) system. Besides, a stress-controlled rotational rheometer with a vane rotor was used as an independent way to measure the yield stress. The results in this work show that in the limit of resolution of the PIV technique, the flow behavior agrees with the existence of a yield stress, but there is a smooth solid?Cliquid transition in the capillary flow curve, which complicates the determination of the yield stress from rheometrical data. This complication, however, is overcome by using the solely velocity profiles and the measured wall shear stresses, from which the yield-stress value is reliably determined. The main details of the kinematics in the presence of slip were all captured during the experiments, namely, a purely plug flow before yielding, the solid?Cliquid transition, as well as the behavior under flow, respectively. Finally, it was found that the slip velocity increases in a power-law way with the shear stress.