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.     

Diffusion effects on the interfacial tension of immiscible polymer blends

Most methods of measuring the interfacial tension between two immiscible polymers are based on the analysis of the shape that a drop of one polymer immersed in the other one exhibits under the action of flow or gravity. In such a situation, the small, yet nonzero mutual solubility between the two polymers acts toward mass transfer between the drop and the surrounding fluid. In this work, diffusion effects on the interfacial tension of the pair polyisobutylene/polydimethylsiloxane have been investigated by drop deformation under shear flow. When the drop was made of polyisobutylene, drop size decreased with time due to diffusion. Drop shrinkage was associated with a significant increase in interfacial tension, until an apparent plateau value was reached. The effect was attributed to a selective migration of molecular weights, which would act to enrich the drop with higher molar mass material. To support such an interpretation, drop viscosity was evaluated by drop shape analysis and it was actually found to increase with time. In some cases, the ratio between drop and continuous phase viscosity became higher than the critical value for drop breakup in shear flow. Upon inverting the phases (i.e., when the drop was made of polydimethylsiloxane), no significant transient effects were observed. In the light of these results, the problem of what are the correct values of interfacial tension and viscosity ratio for a polymer blend of a certain composition will also be discussed. Received: 25 January 1999 Accepted: 24 May 1999     

Determination of interfacial tension by the retraction method of highly deformed drop

The traditional retraction of the deformed drop method (DDRM) to determine the interfacial tension is reformulated to relax the limit the small deformation assumption. The kernel of the new formalism is the calculation of the velocity gradient on the vertex of the ellipsoidal drop. Two models were used for such calculations: the Jackson and Tucker model [J Rheol 47:659–682] and the Yu and Bousmina model [J Rheol 47:1011–1039]. The method can be used either in the retraction of shear deformed drop, or in the retraction of elongated drops produced by the breakup of a long thread. Comparison with experimental results of the literature showed that conversely to the classical DDRM, good accuracy is obtained when the new modeling for the determination of interfacial tension is used both under small and large deformations.     

Shear and squeeze rheometry of suspensions of magnetic polymerized chains

We present the first experimental results on the magnetorheology of suspensions of non-Brownian magnetic ellipsoidal particles. These particles are made of spherical iron particles linked by polymers and are called polymerized chains. Steady shear, oscillatory shear, and oscillatory squeeze rheological tests have been performed. The rheological properties of the suspension of polymerized chains have been compared with those of the suspension of spherical iron particles. In shear flow, both suspensions develop nearly the same yield stress, while in squeeze flow, the yield stress is several times higher for the suspension of polymerized chains. We show that the squeezing force of a suspension of spherical particles is an increasing function of the magnetic field intensity at low magnetic fields but decreases dramatically at higher fields. Surprisingly, this phenomenon, attributed to cavitation or air entrainment, does not occur in the suspension of polymerized chains.     

Viscosity of a Newtonian fluid calculated from the deformation of droplets covered with a surfactant under a linear shear flow

The viscosity of small fluid droplets covered with a surfactant is determined using drop deformation techniques. This method, proposed by Hu and Lips, is here extended to the case of the presence of a surface-active adsorpted at the liquid–liquid interface, to consider more general scenarios. In these experiments, a droplet is sheared by another immiscible fluid of known viscosity, both Newtonian liquids. From the steady-state deformation and retraction mechanisms, the droplet viscosity is calculated using an equation derived from the theories of Taylor and Rallison. Although these theories were expressed for surfactant-free interfaces, they can be applied when a surfactant is present in the system if the sheared droplet reaches reliable steady-state deformations and the surfactant attains its equilibrium adsorption concentration. These determinations are compared to bulk viscosities measured in a rheometer for systems with different viscosity ratios and surfactant concentrations. Very good agreement between both determinations is found for drops more viscous than the continuous phase.     

Strain recovery of model immiscible blends: effects of added compatibilizer

The effect of added compatibilizer on the strain recovery of model immiscible blends after cessation of shear was studied. Blends were composed of polyisobutylene drops (up to 30% by weight) in a polydimethylsiloxane matrix, with viscosity ratio (viscosity of the drops relative to the matrix viscosity) ranging from 0.3 to 1.7. Up to 1% by weight of a PIB-PDMS diblock copolymer was added as compatibilizer. The ultimate recovery recorded after reaching steady-shear conditions increased significantly due to added compatibilizer. Furthermore, the compatibilizer also slowed down the kinetics of the recovery; however, unlike uncompatibilized blends, the recovery could no longer be captured by a single retardation time. The largest increase in ultimate recovery due to compatibilizer occurred at the lowest viscosity ratio. In contrast, the greatest slowing down of the recovery due to compatibilizer occurred at the highest viscosity ratio. The rheological data by themselves are insufficient to reach a definitive conclusion about the mechanism of compatibilizer action. The results are consistent with the effects of flow-induced gradients in compatibilizer concentration. An alternative constitutive modeling approach that captures compatibilizer effects in terms of an interfacial dilational elasticity can reproduce the recovery curves qualitatively, but some predictions of the model contradict experimental results.     

Strain recovery of model immiscible blends without compatibilizer

Strain recovery after the cessation of shear was studied in model immiscible blends composed of polyisobutylene drops (10–30% by weight) in a polydimethylsiloxane matrix. Blends of viscosity ratio (viscosity of the drops relative to the matrix viscosity) ranging from 0.3 to 1.7 were studied. Most of the strain recovery was attributable to interfacial tension, and could be well-described by just two parameters: the ultimate recovery and a single retardation time. Both these parameters were found to increase with the capillary number of the drops prior to cessation of shear. For blends that had reached steady shear conditions, the ultimate recovery decreased with increasing viscosity ratio, whereas the retardation time increased with increasing viscosity ratio. The retardation time was well-predicted, but the ultimate recovery was over-predicted by a linear viscoelastic model developed previously by Vinckier et al. (Rheol Acta 38:65–72, 1999).     

Effect of glycerol addition on phase inversion in horizontal dispersed oil-water pipe flows

The effect of interfacial tension on the phase inversion process during horizontal pipe flow of an oil-aqueous solution was investigated. Interfacial tension was varied by adding small amounts of glycerol in the water phase. At these glycerol concentrations the density and viscosity of the aqueous phase changed by 1% or less. Exxsol™ D140 (5.5 mPa s, 828 kg m⁻³) was used as the oil phase. The experiments were carried out in a 38 mm ID acrylic test pipe. The phase continuity and appearance of phase inversion were investigated using conductivity (wire and ring) probes and an Electrical Resistance Tomographic (ERT) system. The ERT also provided diagrams of the phase distribution in a pipe cross section. Drop size distribution was monitored using a dual impedance probe. It was found that starting from a water continuous flow with increasing oil fraction at constant mixture velocity the mixture inverted initially in the middle of the pipe (measured at 19 mm from the top pipe wall) while a higher oil fraction was required for inversion at the top (measured at 4 mm from the top pipe wall) and finally the rest of the pipe. The addition of glycerol did not affect the phase fraction where the initial inversion occurred but caused an increase in the oil fraction needed to complete the inversion. The drop size measurements were used to explain this behaviour. Pressure drop was found to decrease with increasing oil fraction but this trend reversed when inversion spread to the pipe wall and the oil continuous phase came in contact with it.     

An extended unresolved CFD-DEM coupling method for simulation of fluid and non-spherical particles

This study develops an extended unresolved CFD-DEM coupling method for simulation of the fluid–solid flow with non-spherical particles. The limitation of fluid grid size is discussed, by simulating the settling of a cylinder in a Newtonian fluid based on the resolved and unresolved CFD-DEM coupling method. Then, the calculation of porosity and the fluid–particle relative velocity based on the particle shape enlargement method for simulation of non-spherical particles is proposed. The availability of the particle shape enlargement method for the simulation of non-spherical particles with different sphericity is discussed in this work, by comparing it with the results from the equivalent diameter enlargement method. The limitation of the equivalent diameter enlargement method for non-spherical particles is revealed from the simulation results. Several typical cases are employed to elaborate and verify the extended unresolved CFD-DEM method based on particle shape enlargement method, by presenting a good consistency with the experimental results. It proves that the extended unresolved CFD-DEM method is suitable for different CFD grid size ratios, and consolidates that it is a universal calculation method for CFD-DEM coupling simulation.     

Multi-level DEM study on silo discharge behaviors of non-spherical particles

The silo discharge of non-spherical particles has been widely practiced in engineering processes, yet the understanding of multi-level mechanisms during solid transportation is still lacking. In this study, a high-fidelity super-ellipsoid Discrete Element Method (DEM) model is established to investigate the discharge behaviors of non-spherical particles with different size distributions. After the comprehensive model validations, we investigated the effects of particle shape (aspect ratio and particle sharpness) on the particle level discharge behaviors. The discharge rates of the ellipsoid particles used in the current work are larger than the spherical particles due to the larger solid fraction. The discharge rates of the cuboid-like particles are determined by the combined effect of the solid fraction and the contact force. Parcel level data show that the translational movements of the ellipsoid particles are more ordered, which is supported by the global level data. Strong correlations exist between the particle level and parcel level data, especially the ellipsoid particles and the large particles in the polydispersed cases.     

GPU-enhanced DEM analysis of flow behaviour of irregularly shaped particles in a full-scale twin screw granulator

During twin screw granulation (TSG), small particles, which generally have irregular shapes, agglomerate together to form larger granules with improved properties. However, how particle shape impacts the conveying characteristics during TSG is not explored nor well understood. In this study, a graphic processor units (GPUs) enhanced discrete element method (DEM) is adopted to examine the effect of particle shape on the conveying characteristics in a full scale twin screw granulator for the first time. It is found that TSG with spherical particles has the smallest particle retention number, mean residence time, and power consumption; while for TSG with hexagonal prism (Hexp) shaped particles the largest particle retention number is obtained, and TSG with cubic particles requires the highest power consumption. Furthermore, spherical particles exhibit a flow pattern closer to an ideal plug flow, while cubic particles present a flow pattern approaching a perfect mixing. It is demonstrated that the GPU-enhanced DEM is capable of simulating the complex TSG process in a full-scale twin screw granulator with non-spherical particles.     

Rheological properties of suspensions of TiO2 particles in polymer solutions. 1. Shear viscosity

This paper presents results on the rheological behaviour of suspensions of two kinds of TiO₂ particles in two different polymer solutions. The particles differ in their hydrophilic or hydrophobic properties. The dispersing media are a solution of high molecular weight polyisobutylene in decalin and a solution of a low molecular weight polybutene in decalin. The concentrations of polymer are adjusted in order to get the same zero shear viscosity. The shear viscosity measurements display an apparent yield stress in some cases. The existence and the values of the yield stress depend on the volume fraction of solid particles and on the type of particles. The evolution of the intrinsic viscosity and of the maximum packing fraction vs the shear rate is interpreted in terms of evolution of the size and of the shape of aggregates of particles under shear. The effect of temperature on the development of the yield stress is also discussed. The results are completed by microscopic observations.     

NUMERICAL SIMULATION OF ORIENTATION DISTRIBUTION FUNCTION OF CYLINDRICAL PARTICLE SUSPENSIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｙｌｉｎｄｒｉｃａｌｐａｒｔｉｃｌｅｓｕｓｐｅｎｓｉｏｎｓｉｎｖｏｌｖｅｓｏｍｅｄｉｆｆｉｃｕｌｔｐｒｏｂｌｅｍｓｉｎｔｈｅｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈｏｆｍｕｌｔｉｐｈａｓｅｆｌｏｗａｎｄｎｏｎ＿Ｎｅｗｔｏｎｉａｎｆｌｏｗ ,ａｎｄｈａｖｅｗｉｄｅａｐｐｌｉｃａｔｉｏｎｉｎｔｈｅｉｎｄｕｓｔｒｉａｌｅｎｇｉｎｅｅｒｉｎｇ .Ｆｏｒｅｘａｍｐｌｅ,ｄｕｒｉｎｇｔｈｅｐｒｏｃｅｓｓｉｎｇｏｆｓｏｍｅｍａｔｅｒｉａｌｓ,ｐｒｉｍｉｔｉｖｅｍａｔｅｒｉａｌｐｒｏｐｅ…     

Rheological behaviour of two-phase flow of solid particles in a gas

The present investigation was concerned with the rheological behaviour of dilute suspensions of solid particles in a gas in a vertical cocurrent flow moving upwards. Starting from the experimentally determined dependence of the pressure drop on the concentration of solid particles and the Reynolds number of the carrier medium in the steady flow region, the rheological parameters were estimated using pseudo-shear diagrams. Air was the carrier medium and the dispersed phase was one of six fractions of polypropylene powder and five fractions of glass ballotini. The results show that the investigated two-phase systems have pseudoplastic character which becomes more pronounced with increases in concentration, equivalent diameter and density of solid particles in the flowing suspension. C _d coefficient of particle resistance - d _e equivalent diameter of particles - D column diameter - Fr Froude number - g gravitational acceleration - K rheological parameter - L length - n rheological parameter - p _t pressure drop due to friction - p _m total pressure drop - p _ag pressure drop due to acceleration of the gas phase - p _as pressure drop due to acceleration of the solid phase - p _g hydrostatic pressure of the gas phase - p _s specific effective weight of the dispersed phase - r radius - Re Reynolds number - Re _p Reynolds number of a particle - Re _G generalized Reynolds number - Re _G1 generalized Reynolds number relating to the end of the laminar flow region - Re _G2 generalized Reynolds number relating to the beginning of the turbulent flow region - w _z axial component of velocity - u _t steady free-fall velocity of a single particle - w average velocity - w _g average velocity of the gas phase - w _s average velocity of the dispersed phase of solid particles - relative mass fraction of solid particles - x _s volume fraction of solid particles - _g coefficient of pressure drop due to friction - µ dynamic viscosity - _g density of the gas phase - _m density of the suspension - _s density of solid particles - _ds density of the dispersed phase - _w shear stress at the wall     

Modulation of homogeneous turbulence seeded with finite size bubbles or particles

The dynamics of homogeneous, isotropic turbulence seeded with finite sized particles or bubbles is investigated in a series of numerical simulations, using the force-coupling method for the particle phase and low wavenumber forcing of the flow to sustain the turbulence. Results are given on the modulation of the turbulence due to massless bubbles, neutrally buoyant particles and inertial particles of specific density 1.4 at volumetric concentrations of 6%. Buoyancy forces due to gravity are excluded to emphasize finite size and inertial effects for the bubbles or particles and their interactions with the turbulence. Besides observing the classical entrapment of bubbles and the expulsion of inertial particles by vortex structures, we analyze the Lagrangian statistics for the velocity and acceleration of the dispersed phase. The turbulent fluctuations are damped at mid-range wavenumbers by the bubbles or particles while the small-scale kinetic energy is significantly enhanced. Unexpectedly, the modulation of turbulence depends only slightly on the dispersion characteristics (bubble entrapment in vortices or inertial sweeping of the solid particles) but is closely related to the stresslet component (finite size effect) of the flow disturbances. The pivoting wavenumber characterizing the transition from damped to enhanced energy content is shown to vary with the size of the bubbles or particles. The spectrum for the energy transfer by the particle phase is examined and the possibility of representing this, at large scales, through an additional effective viscosity is discussed.     

Shear viscosity behavior of highly concentrated suspensions at low and high shear-rates

A method is proposed for determining the shear viscosity behavior of highly concentrated suspensions at low and high shear-rates through the use of a formulation that is a function of three parameters signifying the effects of particle size distribution. These parameters are the intrinsic viscosity [], a parametern that reflects the level of particle association at the initiation of motion and the maximum packing concentration _m. The formulation reduces to the modified Eilers equation withn = 2 for high shear rates. An analytical method was used for the calculation of maximum packing concentration which was subsequently correlated with the experimental values to account for the surface induced interaction of particles with the fluid. The calculated values of viscosities at low and high shear-rates were found to be in good agreement with various experimental data reported in literature. A brief discussion is also offered on the reliability of the methods of measuring the maximum packing concentration. _r = /₀ relative viscosity of the suspension - volumetric concentration of solids - k _n coefficient which characterizes a specific effect of particle interactions - _m maximum packing concentration - _r,0 relative viscosity at low shear-rates - [] intrinsic viscosity - n, n parameter that reflects the level of particle interactions at low and high shear-rates, respectively - _r, relative viscosity at high shear-rates - (_m)_s, (_m)_i, (_m)_l packing factors for small, intermediate and large diameter classes - v _s, v_i, v_l volume fractions of small, intermediate and large diameter classes, respectively - _si, _sl coefficient to be used in relating a smaller to an intermediate and larger particle group, respectively - _is, _il coefficient to be used in relating an intermediate to a smaller and larger particle group, respectively - _ls, _li coefficient to be used in relating a larger to a smaller and intermediate particle group, respectively - _m0 maximum packing concentration for binary mixtures - _m,e measured maximum packing concentration - _m,c calculated maximum packing concentration