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.     

Application of Rheological Models in Prediction of Turbulent Slurry Flow

The paper deals with fully developed steady turbulent flow of slurry in a circular straight and smooth pipe. The Kaolin slurry consists of very fine solid particles, so the solid particles concentration, and density, and viscosity are assumed to be constant across the pipe. The mathematical model is based on the time averaged momentum equation. The problem of closure was solved by the Launder and Sharma k-ε turbulence model (Launder and Sharma, Lett Heat Mass Transf 1:131–138, 1974) but with a different turbulence damping function. The turbulence damping function, used in the mathematical model in the present paper, is that proposed by Bartosik (1997). The mathematical model uses the apparent viscosity concept and the apparent viscosity was calculated using two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley. The main aim of the paper is to compare measurements and predictions of the frictional head loss and velocity distribution, taking into account two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley, if the Kaolin slurry possesses low, moderate, and high yield stress. Predictions compared with measurements show an observable advantage of the Herschel–Bulkley rheological model over the Bingham model particularly if the bulk velocity decreases.     

The influence of flow confinement on the rheological properties of complex fluids

The rheological properties of plastic and viscoplastic complex fluids have been investigated using constant velocity squeeze flow rheometry at small gaps in order to examine the influence of the resulting flow confinement. The fluids investigated were aqueous carbopol suspensions and three commercial products (Tesco English mustard, Tesco value lemon curd and Heinz tomato ketchup (HTK)). The bulk rheological properties were measured using parallel plate rheometry. Rough plates were used to eliminate wall slip in both rheometric configurations. The commercial products are Herschel–Bulkley fluids and, for gaps less than a critical value of ~200 μm, the yield stresses tended to large values. There was a corresponding trend to small values for the flow consistencies except for HTK, which has a negligibly small bulk value. A possible explanation is that the micro-structural correlation lengths of these fluids are of the order of that of the critical separation. The yield stresses of the carbopol suspensions also increased sharply at the critical gap but then decreased gradually with decreasing gap separation. The transition may correspond to the formation of a jammed state with the subsequent plastic flow at smaller gap separations arising from localised shearing.     

Analysis on creeping channel flows of compressible fluids subject to wall slip

Creeping channel flows of compressible fluids subject to wall slip are widely encountered in industries. This paper analyzes such flows driven by pressure in planar as well as circular channels. The analysis elucidates unsteady flows of Newtonian fluids subject to the Navier slip condition, followed by steady flows of viscoplastic fluids, in particular, Herschel–Bulkley fluids and their simplifications including power law and Newtonian fluids, that slip at wall with a constant coefficient or a coefficient inversely proportional to pressure. Under the lubrication assumption, analytical solutions are derived, validated, and discussed over a wide range of parameters. Analysis based on the derived solutions indicates that unsteadiness alters cross-section velocity profiles. It is demonstrated that compressibility of the fluids gives rise to a concave pressure distribution in the longitudinal direction, whereas wall slip with a slip coefficient that is inversely proportional to pressure leads to a convex pressure distribution. Energy dissipation resulting from slippage can be a significant portion in the total dissipation of such a flow. A distinctive feature of the flow is that, in case of the pressure-dependent slip coefficient, the slip velocity increases rapidly in the flow direction and the flow can evolve into a pure plug flow at the exit.     

The influence of temperature on rheological properties of blood mixtures with different volume expanders—implications in numerical arterial hemodynamics simulations

During the complicated cardiac surgery on a non-beating heart with cardiopulmonary bypass, protection of the heart is accomplished by injecting cold cardioplegic solutions. In most forms of circulatory shock, it is necessary to immediately restore the circulating volume. Intravenous solutions of volume expanders, such as hydroxyethyl starch and dextrans, are used to increase the volume of fluid in the circulating blood. In this work, blood samples of six donors were obtained and used to prepare mixtures with different volume expanders in concentrations ranging from 10 to 50 vol./vol.%. The flow curves of all mixtures in the temperature range from 4°C to 37°C were constructed and fitted to the Herschel–Bulkley model, in order to extract the shear thinning and yield stress parameters. To assess the influence of the observed changes in the rheological properties of blood on the hemodynamics in arterial vasculature, a realistic three-dimensional rigid-wall computational model was constructed from MRI images of the right carotid bifurcation obtained in vivo from a healthy male volunteer. The time-varying flow field was numerically computed using the Newtonian model as well as the Herschel–Bulkley model with the Papanastasiou regularization. The numerical simulations indicate only moderate changes in the time-averaged flow field that become accentuated when the instantaneous flow field is considered. We also found that although the influence of temperature, hematocrit, and volume expanders on hemodynamics is significant, this can primarily be attributed to the changes in the nominal viscosity of the flow medium.     

Thixotropy, yielding and ultrasonic Doppler velocimetry in pulp fibre suspensions

This paper reports thixotropy in concentrated pulp fibre suspensions and studies their transient flow behaviour using conventional rheometry coupled with a velocimetry technique. Specifically, an ultrasonic Doppler velocimeter is used in conjunction with a rate-controlled rheometer to deduce the local velocity profiles of pulp fibre suspensions. Pulp suspensions are found to exhibit a plateau in their flow curves where a slight increase in the shear stress generates a jump in the corresponding shear rate, implying the occurrence of shear banding. The velocity profiles were found to be discontinuous in the vicinity of the yielding radius where the Herschel–Bulkley model failed to predict the flow. Shear history and the time of rest prior to the measurement were found to play a significant role on the rheology and the local velocity profiles of pulp suspensions.     

Permanent Waves in Slow Free-Surface Flow of a Herschel–Bulkley fluid

Unsteady flow of a viscoplastic fluid on an inclined plane is examined. The fluid is described by the three-parameter Herschel–Bulkley constitutive equation. The set of equations governing the flow is presented, recovering earlier results for a Bingham fluid and steady uniform motion. A permanent wave solution is then derived, and the relation between wave speed and flow depth is discussed. It is shown that more types of gravity currents are possible than in a Newtonian fluid; these include some cases of flows propagating up a slope. The speed of permanent waves is derived and the possible surface profiles are illustrated as functions of the flow behavior index.     

Squeeze-flow and vane rheometry of a gas–liquid foam

The rheology and slip of a dry shaving foam are investigated using squeeze-flow and rotating-vane methods. Constant-force squeeze flow between planar surfaces is used to study the effect of surface roughness on slip and to obtain the yield stress. Non-slip vane measurements are used to obtain the linear shear viscosity and elasticity at small strains, and the yield stress and strain at large strains. Data are compared with the small-strain Maxwell and Kelvin–Voigt linear-viscoelastic models. An apparent dependence of the yield stress and elasticity on the rotational speed of the vane is shown to result from time-dependent rheological parameters as the foam ages. The effect of viscosity in the pre-yield region may give an erroneous identification of yield.     

The rheological characterisation of bread dough using capillary rheometry

An Australian hard wheat flour–water dough has been characterised using parallel plate and capillary rheometers over an extensive range of apparent shear rates (10^− 3–10³ s^− 1) relevant to process conditions. Torsional measurements showed that the shear viscosity of the dough increased with strain to a maximum value and then decreased, suggesting a breakdown of the dough structure. Both torsional and capillary experiments revealed the shear-thinning behaviour of the dough. The wall slip phenomenon in capillary rheometry was investigated and found to be diameter dependent and occurred at a critical shear stress of approximately 5–10 kPa. A two-regime power law behaviour was observed, with the power law index approximately 0.3 in the low shear rate range increasing to 0.67 in the high shear rate range. Pressure fluctuation was observed in the capillary data and increased with shear rate, in particular, at shear rates approaching 10⁴ s^− 1. The results demonstrate that capillary rheometry is a viable means of rheologically testing dough at high shear rates provided pressure fluctuation is carefully monitored and capillary rheometry corrections, including wall slip, are accounted for.     

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     

Axial laminar flow of viscoplastic fluids in a concentric annulus subject to wall slip

The flow of non-Newtonian fluids in annular geometries is an important problem, especially for the extrusion of polymeric melts and suspensions and for oil and gas exploration. Here, an analytical solution of the equation of motion for the axial flow of an incompressible viscoplastic fluid (represented by the Hershel–Bulkley equation) in a long concentric annulus under isothermal, fully developed, and creeping conditions and subject to true or apparent wall slip is provided. The simplifications of the analytical model for Hershel–Bulkley fluid subject to wall slip also provide the analytical solutions for the axial annular flows of Bingham plastic, power-law, and Newtonian fluids with and without wall slip at one or both surfaces of the annulus.     

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.     

Droplet dynamics in sub-critical complex flows

A newly designed eccentric cylinder device has been used to study the deformation and orientation of single Newtonian droplets immersed in an immiscible Newtonian liquid in a controlled complex flow field. Optical microscopy coupled with image acquisition analysis allows monitoring the dynamics of droplets flowing in the gap between the eccentric cylinders. Throughout the experiments, the flow intensity was kept below the critical conditions for droplet break-up. The experimental results are compared with predictions which are obtained using the transient form of the phenomenological model of Maffettone and Minale (J Non-Newtonian Fluid Mech 78:227–241, 1998; J Non-Newtonian Fluid Mech 84:105–106, 1999), incorporating a flow type parameter that accounts for the relative amount of elongational effects in the flow field and adapting the capillary number to mixed flows. For all the sub-critical flows studied here, good agreement was found between model predictions and experimental data, providing, for the first time, a quantitative assessment of drop shape predictions in complex flows.     

Hybrid LES–RANS technique based on a one-equation near-wall model

In order to reduce the high computational effort of wall-resolved large-eddy simulations (LES), the present paper suggests a hybrid LES–RANS approach which splits up the simulation into a near-wall RANS part and an outer LES part. Generally, RANS is adequate for attached boundary layers requiring reasonable CPU-time and memory, where LES can also be applied but demands extremely large resources. Contrarily, RANS often fails in flows with massive separation or large-scale vortical structures. Here, LES is without a doubt the best choice. The basic concept of hybrid methods is to combine the advantages of both approaches yielding a prediction method, which, on the one hand, assures reliable results for complex turbulent flows, including large-scale flow phenomena and massive separation, but, on the other hand, consumes much fewer resources than LES, especially for high Reynolds number flows encountered in technical applications. In the present study, a non-zonal hybrid technique is considered (according to the signification retained by the authors concerning the terms zonal and non-zonal), which leads to an approach where the suitable simulation technique is chosen more or less automatically. For this purpose the hybrid approach proposed relies on a unique modeling concept. In the LES mode a subgrid-scale model based on a one-equation model for the subgrid-scale turbulent kinetic energy is applied, where the length scale is defined by the filter width. For the viscosity-affected near-wall RANS mode the one-equation model proposed by Rodi et al. (J Fluids Eng 115:196–205, 1993) is used, which is based on the wall-normal velocity fluctuations as the velocity scale and algebraic relations for the length scales. Although the idea of combined LES–RANS methods is not new, a variety of open questions still has to be answered. This includes, in particular, the demand for appropriate coupling techniques between LES and RANS, adaptive control mechanisms, and proper subgrid-scale and RANS models. Here, in addition to the study on the behavior of the suggested hybrid LES–RANS approach, special emphasis is put on the investigation of suitable interface criteria and the adjustment of the RANS model. To investigate these issues, two different test cases are considered. Besides the standard plane channel flow test case, the flow over a periodic arrangement of hills is studied in detail. This test case includes a pressure-induced flow separation and subsequent reattachment. In comparison with a wall-resolved LES prediction encouraging results are achieved.      

Estimation of Chaotic and Regular (Stick–Slip and Slip–Slip) Oscillations Exhibited by Coupled Oscillators with Dry Friction

In this paper, we present a novel approach to quantify regular or chaotic dynamics of either smooth or non-smooth dynamical systems. The introduced method is applied to trace regular and chaotic stick–slip and slip–slip dynamics. Stick–slip and slip–slip periodic and chaotic trajectories are analyzed (for the investigated parameters, a stick–slip dynamics dominates). Advantages of the proposed numerical technique are given.     

Experimental and numerical investigation of a viscoplastic Carbopol gel injected into a prototype 3D mold cavity

A numerical model for free-surface flow of a viscoplastic liquid into a cavity is presented. This flow is regarded as a basic model of injection molding, which is a widely used processing technology. Model experiments of the injection process are performed with a water-based gel with shear-thinning behavior. The filling process is visualized by tracing the free surface of the gel within the cavity. Filling times of the cavity are deduced from the experimental observations. The filling process is also analyzed by means of numerical simulation.The flow equations are integrated according to the finite-volume method. The volume-of-fluid method is employed in order to describe the flow of two incompressible, immiscible phases, the phase interface is resolved by the method of geometric reconstruction or alternatively by the method of surface compression. The Herschel–Bulkley model is used in order to describe the shear-thinning behavior of the gel and the effects of a yielding point. The governing equations of the flow are solved by means of the commercial code Fluent as well as the Open Source code OpenFOAM.The results of the numerical simulations are analyzed in detail. They are compared with the experimental findings. Cavity filling times in the experiments and the simulations are in good agreement. Different patterns of the filling flow depending on the injection parameters are evident in the experiments and the simulations. They are characterized and arranged with respect to the similarity parameters of the flow configuration. Again, the results of the simulation are found to agree well with the experimental observations.     

An outline of the non-linear viscoelastic behaviour of wheat flour dough in shear

Creep and creep recovery, stress relaxation and small- and large-amplitude oscillatory shear experiments have been used to study the steady-state flow behaviour and the transient viscoelastic response of wheat flour dough in shear over large ranges of time, stress and strain. The results are discussed with reference to the limited body of reliable literature data. Dough does display a linear viscoelastic domain. The complex character of its non-linear viscoelastic properties is essentially due to the extremely low shear rate limit of the initial Newtonian plateau and to the onset of time-dependent flow behaviour above a certain strain threshold, which explain qualitatively the discrepancies observed in certain cases on a part of the range of the rheological variables explored, despite global self-consistency of the results. Comparison of gluten and dough linear viscoelastic properties shows that dough cannot be viewed simply as a concentrated suspension of starch granules in the hydrated viscoelastic gluten matrix.Paper presented at the second Annual European Rheology Conference (AERC 2005) held in Grenoble, France on April 21–23, 2005.     

On the pressure and stress singularities induced by steady flows of incompressible viscous fluids

Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, to date there has been relatively little explicit identification of stress singularities caused by fluid flows. In this study, stress and pressure singularities induced by steady flows of viscous incompressible fluids are asymptotically identified. This is done by taking advantage of an earlier result that the Navier-Stokes equations are locally governed by Stokes flow in angular corners. Findings for power singularities are confirmed by developing and using an analogy with solid mechanics. This analogy also facilitates the identification of flow-induced log singularities. Both types of singularity are further confirmed for two global configurations by applying convergence-divergence checks to numerical results. Even though these flow-induced stress singularities are analogous to singularities in solid mechanics, they nonetheless render a number of structural configurations singular that were not previously appreciated as such from identifications within solid mechanics alone.     

Numerical study of the combined effects of inertia,slip, and compressibility in extrusion of yield stress fluids

The axisymmetric extrudate swell flow of a compressible Herschel–Bulkley fluid with wall slip is solved numerically. The Papanastasiou-regularized version of the constitutive equation is employed, together with a linear equation of state relating the density of the fluid to the pressure. Wall slip is assumed to obey Navier’s slip law. The combined effects of yield stress, inertia, slip, and compressibility on the extrudate shape and the extrudate swell ratio are analyzed for representative values of the power-law exponent. When the Reynolds number is zero or low, swelling is reduced with the yield stress and eventually the extrudate contracts so that the extrudate swell ratio reaches a minimum beyond which it starts increasing asymptotically to unity. Slip suppresses both swelling and contraction in this regime. For moderate Reynolds numbers, the extrudate may exhibit necking and the extrudate swell ratio initially increases with yield stress reaching a maximum; then, it decreases till a minimum corresponding to contraction, and finally, it converges asymptotically to unity. In this regime, slip tends to eliminate necking and may initially cause further swelling of the extrudate, which is suppressed if slip becomes stronger. Compressibility was found to slightly increase swelling, this effect being more pronounced for moderate yield stress values and wall slip.     

Simultaneous determination of shear and extensional viscosities and wall slip by on-line rheometry with parameter fitting: an application to EPDM rubber

Shear and extensional viscosities and wall slip are determined simultaneously under extrusion processing conditions using an on-line rheometer. Because it is not possible to independently control flow rate and temperature, classical methods for interpretation of capillary data cannot be used with on-line rheometry. This limitation is overcome using computational optimization to fit parameters in a flow model. This consists of three parts, representing shear viscosity, extensional viscosity, and wall slip. Three-parameter, power law forms, based on local instantaneous deformation rates and including temperature dependence, are used for each, and analytic solutions applied for entry flow and flow in the capillary. For entry flow, the Cogswell–Binding approach is used, and for developed flow in the capillary a solution incorporating wall slip is derived. The rheometer, with interchangeable capillaries, is mounted in place of the die on a rubber profile extrusion line. Pressure drops and temperatures for extrusion of an EPDM rubber through 2 mm diameter capillaries of length 0, 2, 3, 4, and 5 mm are logged and flow rates determined for a range of extruder speeds (5 to 20 rpm). Pressures ranged from 60 to 75 bar and temperatures from 86 to 116 °C. Mean flow velocity in the capillaries was between 5 × 10⁻³ and 5 × 10⁻¹ m s⁻¹. The nine material parameters are optimized for best fit of the analytic pressure drops to experimental data, using about 100 data points, with the Levenberg–Marquardt method. It is concluded that flow is dominated by extension and wall slip. Shear flow appears to play little part. The slip model indicates that slip velocity increases much more rapidly than the wall shear stress (in the range 0.5–1 MPa) and decreases with temperature for a given stress level. Results for the (uniaxial) extensional viscosity represent an engineering approximation to this complex phenomenon at the high strains (approximately 200) and high extension rates (up to 800 s⁻¹) applying in the extrusion. Results indicate a slight extension hardening and a decrease with temperature. Results are put into the context of the available studies in the literature, which, particularly with regard to wall-slip and extensional flow, consider conditions far removed from those applying in industrial extrusion. The present methods provide a powerful means for flow characterization under processing conditions, providing data suitable for use in computer simulations of extrusion and optimization of die design.