Rheological properties of suspensions of the green microalga Chlorella vulgaris at various volume fractions

A systematic study of the rheological properties of solutions of non-motile microalgae (Chlorella vulgaris CCAP 211-19) in a wide range of volume fractions is presented. As the volume fraction is gradually increased, several rheological regimes are observed. At low volume fractions (but yet beyond the Einstein diluted limit), the suspensions display a Newtonian rheological behaviour and the volume fraction dependence of the viscosity can be well described by the Quemada model (Quemada, Eur Phys J Appl Phys 1:119–127, 1997). For intermediate values of the volume fraction, a shear thinning behaviour is observed and the volume fraction dependence of the viscosity can be described by the Simha model (Simha, J Appl Phys 23:1020–1024, 1952). For the largest values of the volume fraction investigated, an apparent yield stress behaviour is observed. Increasing and decreasing stress ramps within this range of volume fractions indicate a thixotropic behaviour as well. The rheological behaviour observed within the high concentration regime bears similarities with the measurements performed by Heymann and Aksel (Phys Rev E 75:021505, 2007) on polymethyl methacrylate suspensions: irreversible flow behaviour (upon increasing/decreasing stresses) and dependence of the flow curve on the characteristic time of forcing (the averaging time per stress values). All these findings indicate a behaviour of the microalgae suspensions similar to that of suspensions of rigid particles. A deeper insight into the physical mechanisms underlying the shear thinning and the apparent yield stress regime is obtained by an in situ analysis of the microscopic flow of the suspension under shear. The shear thinning regime is associated to the formation of cell aggregates (flocs). Based on the Voronoi analysis of the correlation between the cell distribution and cell sizes, we suggest that the repulsive electrostatic interactions are responsible for this microscale organisation. The apparent yield stress regime originates in the formation of large-scale cell aggregates which behave as rigid plugs leading to a maximally random jammed state.     

Dynamic simulation of non-spherical particulate suspensions

Particle-level simulation has been employed to investigate rheology and microstructure of non-spherical particulate suspensions in a simple shear flow. Non-spherical particles in Newtonian fluids are modeled as three-dimensional clusters of neutrally buoyant, non-Brownian spheres linked together by Hookean-type constraint force. Rotne–Prager correction to velocity disturbance has been employed to account for far-field hydrodynamic interactions. An isolated rod-like particle in simple shear flow exhibits a periodic orientation distribution, commonly referred to as Jeffery orbit. Lubrication-like repulsive potential between clusters have been included in simulation of rod-like suspensions at various aspect ratios over dilute to semi-dilute volume fractions. Shear viscosity evaluated by orientation distribution qualitatively agrees with one obtained by direct computation of shear stress.     

A model for the shear viscosity of non-colloidal suspensions with Newtonian matrix fluids

We present a model for the shear viscosity of non-colloidal suspensions with Newtonian matrix fluids. The model is based on the original idea first presented by Brinkman (Applied Sci Research A1:27-34. 1947) for the viscous force exerted by a flowing fluid on a dense swarm of spherical particles. In particular, we consider an inertialess suspension in which the mean flow is driven by a pressure difference, and simultaneously, the suspension is subject to simple shear. Assuming steady state, incompressibility and taking into account a resistance force which is generated due to the presence of the particles in the flow, the three-dimensional governing equations for the mean flow around a single spherical particle are solved analytically. Self-consistency of the model provides a relationship between the resistance parameter and the volume fraction of the solid phase. A volume, or an ensemble, averaging of the total stress gives the bulk properties and an expression for the relative (bulk) viscosity of the suspension. The viscosity expression reduces to the Einstein limit for dilute suspensions and agrees well with empirical formulas from the literature in the semi-dilute and concentrated regimes. Since the model is based on a single particle and its average interaction with the other particles is isotropic, no normal stress differences can be predicted. A possible method of addressing this problem is provided in the paper.     

Shear rheometry and visualization of glass fiber suspensions

The nonlinear rheological behavior of short glass fiber suspensions has been investigated in this work by rotational rheometry and flow visualization. A Newtonian and a Boger fluid (BF) were used as suspending media. The suspensions exhibited shear thinning in the semidilute regime and weaker shear thinning in the transition to the concentrated one. Normal stresses and relative viscosity were higher for the BF suspensions than for the Newtonian ones presumably due to enhanced hydrodynamic interactions resulting from BF elasticity. In addition, relative viscosity of the suspensions increased rapidly with fiber content, suggesting that the rheological behavior in the concentrated regime is dominated by mechanical contacts between fibers. Visualization of individual fibers and their interactions under flow allowed the detection of aggregates, which arise from adhesive contacts. The orientation states of the fibers were quantified by a second order tensor and fast Fourier transforms of the flow field images. Fully oriented states occurred for shear rates around 20 s^− 1. Finally, the energy required to orient the fibers was higher in step forward than in reversal flow experiments due to a change in the spatial distribution of fibers, from isotropic to planar oriented, during the forward experiments.     

On the behaviour of gelled fibre suspensions in steady shear

The shear rheological properties of suspensions of gelled agar fibres in a low viscosity Newtonian matrix fluid were investigated. Two classes of fibres, low aspect ratio fibres and high aspect ratio fibres with an aspect ratio of the order of 10 and 100 respectively were included in the investigations. For all fibre phase volumes investigated, from as low as 0.01 upwards, the flow curves are characterised by an apparent yield stress followed by shear-thinning which was independent of the fibre aspect ratio. Based on our analysis of the flow curves, we conclude that the high aspect ratio fibres behave like flexible threads in contrast to the low aspect ratio fibres whose high shear relative viscosity is successfully described by a relation for long rigid rods. These findings are supported by flow visualisation using an optical shearing stage coupled to a light microscope.     

Der Einfluß der elastischen Eigenschaften des Trägermediums auf das Fließverhalten von Kugel- und Fasersuspensionen

The influence of the elastic properties of the suspending medium on the flow of viscoelastic glass bead and glass fibre suspensions through flat orifices was investigated. The results are discussed by contrasting the flow behaviour of the viscoelastic suspensions with that of corresponding suspensions in a Newtonian suspending medium. For the suspensions in a Newtonian oil linear relationships were always found between the pressure loss and the effective velocity gradient in orifice flow. Thus it can be concluded that in this case the influence of the filler on the flow behaviour is independent of the imposed strain. Increasing the filler content or using more anisotropic particles led to higher viscosities and thus to larger pressure losses. It is well known that viscoelastic polyisobutene solutions show strain rate dependent flow behaviour. Due to the increasing influence of elasticity with increasing strain rate, the apparent flow curves can be divided into characteristic regions with different slopes. The addition of filler to such solutions altered the shape of the flow curves and it was found that the onset of “flow hardening” occurred at lower imposed strains. In addition, characteristic changes in the hardening behaviour and flow stability were observed; these were most pronounced for the fibre suspensions, even at low concentrations. For the fibre suspensions, these phenomena could be related to the influence of the fibres on the undisturbed flow field near the orifice, leading, in general, to higher strain velocities between the fibres. On the other hand, enhanced extensional strains are induced at the ends of the fibres and, at the same time, shear flow occurs along the fibres. These two effects counteract each other with respect to the hardening behaviour of the polymer. Similar effects may also occur in the glass bead suspensions, although they would certainly be less pronounced.     

A parametric study of droplet deformation through a microfluidic contraction: Shear thinning liquids

Numerical simulations of a droplet passing through an axisymmetric microfluidic contraction are presented, focusing on systems where one of the two liquids present is shear thinning. The simulations are performed using a transient Volume of Fluid (VOF) algorithm. When the droplet is shear thinning and the surrounding phase Newtonian, droplets deform in a similar way to Newtonian droplets that have a viscosity equal to the average viscosity of the shear thinning fluid while it is within the contraction. When the surrounding phase is shear thinning and the droplet Newtonian, droplets deform in a similar way to droplets contained within a Newtonian liquid that has a viscosity that is lower than that of the droplet. In both cases the behaviour of the shear thinning fluid can be broadly described in terms of a ‘characteristic’ Newtonian viscosity: However, determining the exact value of this viscosity without performing a full shear thinning simulation is not possible.     

A self-similar behavior for the relative viscosity of concentrated suspensions of rigid spheroids

A viscosity model for suspensions of rigid particles with predictive capability over a wide range of particle volume fraction and shear conditions is of interest to quantify the transport of suspensions in fluid flow models. We study the shear viscosity of suspensions and focus on the effect of particle aspect ratio and shear conditions on the rheological behavior of suspensions of rigid bi-axially symmetric ellipsoids (spheroids). We propose a framework that forms the basis to microscopically parameterize the evolution of the suspension microstructures and its effect on the shear viscosity of suspensions. We find that two state variables, the intrinsic viscosity in concentrated limit and the self-crowding factor, control the state of dispersion of the suspension. A combination of these two variables is shown to be invariant with the imposed shear stress (or shear rate) and depends only on the particle aspect ratio. This self-similar behavior, tested against available experimental and numerical data, allows us to derive a predictive model for the relative viscosity of concentrated suspensions of spheroids subjected to low (near zero) strain rates. At higher imposed strain rates, one needs to constrain one of the state variables independently to constrain the state of dispersion of the suspension and its shear dynamic viscosity. Alternatively, the obtained self-similar behavior provides the means to estimate the state variables from the viscosity measurements made in the laboratory, and to relate them to microstructure rearrangements and evolution occurring during deformation.     

Extensional rheology of a shear-thickening cornstarch and water suspension

A filament-stretching rheometer is used to measure the extensional viscosity of a shear-thickening suspension of cornstarch in water. The experiments are performed at a concentration of 55 wt.%. The shear rheology of these suspensions demonstrates a strong shear-thickening behavior. The extensional rheology of the suspensions demonstrates a Newtonian response at low extension rates. At moderate strain rates, the fluid strain hardens. The speed of the strain hardening and the extensional viscosity achieved increase quickly with increasing extension rate. Above a critical extension rate, the extensional viscosity goes through a maximum and the fluid filaments fail through a brittle fracture at a constant tensile stress. The glassy response of the suspension is likely the result of jamming of particles or clusters of particles at these high extension rates. This same mechanism is responsible for the shear thickening of these suspensions. In capillary breakup extensional rheometry, measurement of these suspensions demonstrates a divergence in the extensional viscosity as the fluid stops draining after a modest strain is accumulated.     

Shear behaviour of biopolymer suspensions with spheroidal and cylindrical particles

The rheological behaviour of suspensions is influenced by many parameters, one of which is the particle shape. For rigid particle suspensions a number of studies demonstrate the effects of the particle aspect ratio. Indeed, fibres are widely used as rheology modifiers in different materials such as synthetic polymers. This work is concerned with testing the hypothesis that regularly shaped particles with aspect ratios larger than one that are made of gelled biopolymers could be used as rheology modifiers for biopolymer solutions. Biopolymers, and mixtures thereof are a widely used ingredient in foods and other products with structure functionality. Tailoring rheology modifiers by morphology offers an alternative to using different biopolymers. It is demonstrated how biopolymer suspensions with regular spheroidal, or cylindrical particle shapes can be produced by gelling the droplet phase of a liquid two phase biopolymer mixture in a shear field. Biopolymers were chosen such that gelation is initiated by cooling. Shear-cooling at constant stresses leads to the formation of ellipsoidal particles. Cylindrical particles can be generated by stepping up the shear stress prior to gelation, i.e., stretching the droplet phase into fibrils, and trapping the shape prior to break-up through gelation. Morphologies and steady shear rheological data for suspensions of the two biopolymers gellan and κ-carrageenan with an internal phase volume of 0.2 are reported. The influence of particle shape on relative viscosity is pronounced. At high shear stresses particle orientation leads to decreased viscosity with increasing particle aspect ratio. In the low shear region, higher aspect ratio suspensions show higher viscosities. Additionally, the material properties, including the interfacial tension, which influence the suspension morphology are reported. Received: 3 March 2000 Accepted: 22 August 2000     

Reological behaviour of concentrated monodisperse suspensions as a function of preshear conditions and temperature: an experimental study

 The influence of preshearing on the rheological behaviour of model suspensions was investigated with a stress-controlled cone-and-plate rheometer. The used matrix fluids showed Newtonian behaviour over the whole range of applied shear stresses. Highly monodisperse spherical glass spheres with various particle diameters were used as fillers. By applying steady preshearing at a low preshear stress, where a diffusion of particles can be expected, it was found for all model suspensions investigated at volume fractions ranging from 0.20 to 0.35 that the time-temperature superposition in the steady shear and in the dynamic mode holds within the chosen temperature range. Furthermore, all presheared model suspensions displayed a high and a low frequency range which are either separated by a shoulder or by a plateau value of G′ at intermediate frequencies. It could clearly be demonstrated that the low frequency range strongly depends on the preshear conditions. Hence, the features observed in the low frequency range can be attributed to a structure formation of a particulate network. In the high frequency range a frequency-dependent behaviour was observed which obeys the classical behaviour of Newtonian fluids (G′∝ω², G′′∝ω). The resulting temperature shift factors from the dynamic and the steady shear mode are identical and independent of the volume fraction and the particle size of the filler. Received: 29 November 2000 Accepted: 12 July 2001     

Rheology of concentrated suspensions containing mixtures of spheres and fibres

Optimising flow properties of concentrated suspensions is an important issue common for many industries. The rheology of concentrated suspensions has therefore been studied intensively both experimentally and theoretically. Most studies have focused on monodisperse and polydisperse suspensions of either spheres or fibres. In practice, most suspensions contain particles that are polydisperse both in size and shape. A mixing rule for such systems is expected to be a powerful tool for engineers and product designers. Therefore in this work, suspensions of spheres, fibres and mixtures thereof were characterised using rotational shear rheometry and in-line image analyses. Thereby, total solids volume concentration and fibre fraction was varied. Results from transient and steady-state shear rheometry are discussed with respect to concentration, fibre fraction, and shear induced microstructure. Experimentally obtained viscosity data were accurately fitted using the model proposed by Farris (T Soc Rheol 12:281, 1968) for mixtures of monodisperse non-interacting spheres of different sizes.Originally presented at the Annual European Rheology Conference 2003, AERC 2003     

On the elastic properties of model suspensions as investigated by creep recovery measurements in shear

 The elastic properties of model suspensions with spherical monodisperse hydrophilic glass spheres that were dispersed in a Newtonian liquid were determined in creep and creep recovery measurements in shear with a magnetic bearing torsional creep rheometer. The creep and creep recovery measurements were performed depending on the applied level of shear stresses ranging from 0.19 Pa to 200 Pa. Since the recoverable creep compliances of the chosen suspending medium (i.e. a low molecular weight polyisobutylene) were far below the lower limit of the resolution of the creep rheometer it can be considered to behave as purely viscous. By applying a large shear stress in the creep tests the investigated suspensions with a volume fraction of Φ _t =0.35 behave as Newtonian liquids, too. For these suspensions no significant recoverable creep compliances could be detected, as well. In contrast to the Newtonian state of suspensions at high shear stresses, where a shear induced ordering of the particles can be expected, a non-Newtonian behaviour arises by applying a very low shear stress in the creep test. In this state large recoverable creep compliances were detected for the suspensions. The magnitude of the recoverable creep compliances of the suspensions exceeded the largest creep compliances of polymer melts that are reported in the literature by more than two decades. From the results obtained by creep recovery measurements with a magnetic bearing torsional creep rheometer it can clearly be concluded that the particle structure present in the chosen model suspension gives rise to a pronounced elasticity. Received: 21 November 2000 Accepted: 12 July 2001     

Experiments on the effective viscosity of concentrated suspensions of solid spheres

The effective viscosity of concentrated suspensions of solid spherical particles was determined experimentally, under zero shear conditions, by means of a convection experiment. For solid concentrations in the range 30–50 per cent the results verify extrapolations to low shear rates of previously reported conventional viscosity measurements, and suggest that, in this range of concentrations and at low shear, a suspension of neutrally bouyant spheres behaves effectively as an equivalent Newtonian fluid.     

A numerical simulation of suspension flow using a constitutive model based on anisotropic interparticle interactions

We report a Brownian configuration field implementation of a recent constitutive equation for suspensions, reported by Phan-Thien et al. 1999. The numerical method is a hybrid technique, which combines a modification of the Brownian configuration field method described by Hulsen et al. 1997 and the adaptive viscosity split stress formulation proposed by Sun et al. 1996. The implementation is used to examine the flow past a sphere in a tube. The relative viscosity derived from the drag force/sedimentation velocity agrees well with a well-known empiricism. In addition, the ratio of the pressure force to the drag on the sphere seems to be weakly dependent on the volume fraction, and is somewhat higher than Brenner's results of 1962, which were derived for Newtonian fluids. Received: 5 April 1999/Accepted: 27 September 1999     

Überlegungen zur Scherviskosität von Suspensionen aufgrund einer Dimensionsanalyse

Dimensional analysis of the motion of solid particles suspended in a fluid phase shows that the macroscopic relative shear viscosity of suspensions generally depends not only on the volume concentration and particle shape but also on two Reynolds numbers and a dimensionless sedimentation number. These dimensionless numbers are formed using parameters characterizing the structure and motion of the suspension at the microscopic level. The analysis was based on the assumptions that the dispersed particles are rigid and sufficiently large that Brownian motion may be neglected, that the continuous fluid phase is Newtonian and that the interactions between particles and between particles and fluid phase are only hydrodynamic. The Reynolds numbers describe the influence of the inertial forces at the microscopic level, and the sedimentation number the influence of gravity. The dimensionless numbers can be neglected if their values are much smaller than one. For each of the dimensionless numbers both the shear rate and the particle size influence the shear viscosity. Thus sedimentation number is large for low shear rates, whereas the Reynolds numbers are large for high shear rates. The viscosity function for one suspension can be transformed into the viscosity function for another suspension with geometrically similar particles but of a different size. The scale-up rules are derived from the requirement that the relevant dimensionless numbers must be constant. The influence of non-hydrodynamic effects at the microscopic level on the shear viscosity can be detected by deviations from the derived scale-up rules.     

Droplet deformation of a strongly shear thinning dense suspension of polymeric micelles

We investigated the deformation of a strong shear thinning droplet undergoing simple shear flow in a Newtonian liquid. The droplet was an aqueous solution of poly(ethylene oxide) end capped with an alkyl group that forms spherical micelles in aqueous solution. At high concentrations and below a critical temperature, the jammed micelles showed strong shear thinning behaviour, and neither a yield stress nor a Newtonian viscosity was observed. At small shear rates, the droplet rotated and aligned in the flow, but did not deform or only very weakly. At high shear rates, the droplet deformation increased with increasing shear rate. The deformed droplet did not relax after the shear was stopped except for a modest rounding of the edges. For each shear rate, an apparent viscosity, η _ad, of the equivalent Newtonian droplet was calculated assuming affine deformation. η _ad showed a power law dependence on the capillary number Ca with an exponent of − 1.8 and was larger than the shear viscosity of the micelle suspension at the same shear rates. The results were explained by the existence of a strong gradient of the viscosity inside the droplet leading to a very low viscosity fluid layer near the droplet/matrix interface.     

A generalized equation for rheology of emulsions and suspensions of deformable particles subjected to simple shear at low Reynolds number

We present analyses to provide a generalized rheological equation for suspensions and emulsions of non-Brownian particles. These multiparticle systems are subjected to a steady straining flow at low Reynolds number. We first consider the effect of a single deformable fluid particle on the ambient velocity and stress fields to constrain the rheological behavior of dilute mixtures. In the homogenization process, we introduce a first volume correction by considering a finite domain for the incompressible matrix. We then extend the solution for the rheology of concentrated system using an incremental differential method operating in a fixed and finite volume, where we account for the effective volume of particles through a crowding factor. This approach provides a self-consistent method to approximate hydrodynamic interactions between bubbles, droplets, or solid particles in concentrated systems. The resultant non-linear model predicts the relative viscosity over particle volume fractions ranging from dilute to the the random close packing in the limit of small deformation (capillary or Weissenberg numbers) for any viscosity ratio between the dispersed and continuous phases. The predictions from our model are tested against published datasets and other constitutive equations over different ranges of viscosity ratio, volume fraction, and shear rate. These comparisons show that our model, is in excellent agreement with published datasets. Moreover, comparisons with experimental data show that the model performs very well when extrapolated to high capillary numbers (C a?1). We also predict the existence of two dimensionless numbers; a critical viscosity ratio and critical capillary numbers that characterize transitions in the macroscopic rheological behavior of emulsions. Finally, we present a regime diagram in terms of the viscosity ratio and capillary number that constrains conditions where emulsions behave like Newtonian or Non-Newtonian fluids.     

Rheology of oil-in-water emulsions

The effect of interfacial tension on the steady-flow and dynamic viscoelastic behavior of emulsions are studied experimentally. At very low inter-facial tensions and low volume fractions, the viscosity decreases with increasing shear rate and becomes constant at high shear rates. The high-shear-rate Newtonian viscosity is not affected by interfacial tension, but the transition from pseudoplastic to Newtonian flow shifts to lower shear rates as the interfacial tension decreases. At an interfacial tension of 5 × 10^–3 Nm^–1, the viscosity decreases, passes through a minimum, and then increases as the shear rate is increased. The dilatant behavior may be attributed to elastic responses of interfaces during collision of drops. At high volume fractions, the emulsions show remarkable elasticity resulting from the interfacial energy associated with deformation of liquid films. The modulus and viscosity are proportional to interfacial tension and inversely proportional to drop size.     

Resistance to uniaxial extensional flow of fibre suspensions

The tensile stress due to resistance to uniaxial extensional flow of fibre suspensions in Newtonian and non-Newtonian fluids has been measured using the filament stretching technique. It has been found that addition of fibres to a Newtonian fluid increases the extensional viscosity. The steady state results agree with Bachelors theory and the stress growth behaviour is qualitatively predicted by the theory of Dinh and Armstrong. Experimental results from this work have also shown that the behaviour of a fibre suspension in viscoelastic fluid is qualitatively described by Fans equation. The added fibres increase the extensional stress growth coefficient of the viscoelastic fluid at low strain but have marginal effect on the fluid after the onset of strain-hardening.