Slow viscoelastic radial flow between parallel disks

A perturbation analysis is presented for the steady-state radial flow of a third-order fluid between two parallel disks. The results include previous perturbation analyses in which various other rheological models were used. The pressure drop needed to maintain the radial flow is less than that for the Newtonian creeping-flow solution because of fluid inertia and shear-thinning viscosity, whereas the normal stresses have the opposite effect. Possible use of the “radial flow viscometer” for experimental evaluation of second-order constants is also discussed. Finally, molecular stretching in the flow system is examined using the elastic dumbbell model for a polymer molecule.     

Numerical simulation of blade coating with short dwell and roll application coaters

Paper and board are often coated at high speeds with a mineral-based aqueous suspension in order to improve their printing properties. This suspension is usually called coating colour. The flow behaviour of the coating colour in the cavity of the short dwell coater (SDC) and in the vicinity of the blade tip when paper is coated with a stiff blade has been analysed using the finite element method. The models used to simulate the flow incorporated free surfaces and shear-thinning colours. The Newtonian case was in some cases also included in the modelling. The viscosity level and the shear-thinning character of the coating colour had a significant influence on the flow in the SDC cavity, although the overall behaviour was to a large extent governed by the speed of the coater. The pressure distribution along the paper surface in the SDC cavity was also analysed.In agreement with earlier reported results, increasing the machine speed raised the pressure level in the colour close to the blade tip. The rheological properties of the coating colour also affected the flow field and the pressure distribution in this region; e.g. it was found here that a colour with a high viscosity level at low shear rates developed a high pressure level close to the blade tip. The most interesting result revealed by the analysis was that changes in the configuration close to the blade tip (converging flow between the blade tip and the paper and compression of the substrate under the blade tip) and boundary conditions at the blade tip surface had a very significant effect on the pressure distribution. The predictions of the numerical simulation were to some extent compared with experience from practical coating trials.     

Migration of a sphere suspended in viscoelastic liquids in Couette flow: experiments and simulations

The intrinsically coupled effects of the curvature of the flow-field and of the viscoelastic nature of suspending medium on the cross-stream lateral migration of a single non-Brownian sphere in wide-gap Couette flow are studied. Quantitative videomicroscopy experiments using a counterrotating device are compared to the results of 3D finite element simulations. To evaluate the effects of differences in rheological properties of the suspending media, fluids have been selected which highlight specific constitutive features, including a reference Newtonian fluid, a single relaxation time wormlike micellar surfactant solution, a broad spectrum shear-thinning elastic polymer solution and a constant viscosity, highly elastic Boger fluid. As expected for conditions corresponding to Stokes flow, migration is absent in the Newtonian fluid. In the wormlike micellar solution and the shear-thinning polymer solution, spheres are observed to migrate in the direction of decreasing shear rate gradient, i.e. the outer cylinder, except when the sphere is initially released close to the inner cylinder, in which case the migration is towards it. The migration is enhanced by faster relative angular velocities of the cylinders. Shear-thinning reduces the migration velocity, showing an opposite behavior as compared to previous results in planar shear flow. In the Boger fluid, within experimental error no migration could be observed, likely due to the large solvent contribution to the overall viscosity. For small Deborah numbers the migration results are well described by an heuristic argument based on a local stress balance.     

Viscoelastic properties of magnetorheological fluids

We have studied the rheological properties of some magnetorheological fluids (MRF). MRF are known to exhibit original rheological properties when an external magnetic field is applied, useful in many applications such as clutches, damping devices, pumps, antiseismic protections, etc. While exploiting parameters such as magnetic field intensity, particle concentration and the viscosity of the suspending fluid, we highlighted the importance of each one of these parameters on rheology in the presence of a magnetic field. We made this study by conducting rheological experiments in dynamic mode at very low strain which facilitates the comprehension of the influence of the structure on MRF rheology. Our results confirmed the link between the magnetic forces which ensure the cohesion of the particles in aggregates, and the elastic modulus. Moreover, we found that the loss modulus varies with the frequency in a similar manner than the elastic modulus. The system, even with the smallest deformations, was thus not purely elastic but dissipates also much energy. Moreover, we demonstrated that this dissipation of energy was not due to the matrix viscosity. Actually, we attributed viscous losses to particle movements within aggregates.     

Steady shear viscosity of short fibre suspensions in thermoplastics

The rheological behaviour of a polyethylene, two polyamides and a silicone oil filled with different fibre contents are studied in capillary rheometry. The viscosity increase induced by the fibres is important for the silicone oil, and negligible for the polyethylene. The polyamide is intermediate. The same classification stands for the pressure loss in the convergent channel upstream from the capillary. A constitutive equation based on a cell model which takes into account the shear-thinning behaviour of the matrix is built. The predictions of the model are in correct agreement with the measurements. Received: 22 September 1997 Accepted: 5 March 1998     

Numerical simulation of shear-thinning flow problems in mixing vessels

In many applications of chemical and process engineering numerous important flow situations appear, in which the shear-thinning properties of the fluid dominate the normal stress effects and time-dependent elastic effects. A mixed finite element method for steady incompressible flow is presented taking as a basis the model of a generalized Newtonian fluid. An augmented Lagrangian functional is constructed corresponding to the equations of motion and the continuity constraint together with appropriate boundary conditions. The treatment of the resulting nonlinear system of equations by the Newton-Raphson scheme is made for general models of the viscosity function. The linear subproblems are solved by Uzawa's algorithm. The method is applied to the numerical simulation of various mixing problems in cylindrical unbaffled vessels. The computations were performed for a real polymer liquid (CMC in water), where the viscosity data were found experimentally and were fitted to a suitable mathematical model. Several numerical results are discussed and as far as possible compared with experimental data.This work was supported by Deutsche Forschungsgemeinschaft.     

3D particle image velocimetry of the flow field around a sphere sedimenting near a wall: Part 1. Effects of Weissenberg number

The flow fields surrounding a sphere sedimenting through a liquid near a vertical wall are characterized using 3D stereoscopic particle-image velocimetry (PIV) experiments. Three different fluids, a Newtonian reference fluid, a constant (shear) viscosity Boger fluid, and a shear-thinning elastic fluid, are used to determine the effects of both elasticity and shear-thinning on the flow field. All three fluids have similar zero shear viscosities. The Weissenberg number is manipulated by varying the diameter and the composition of the ball. Significant differences are found for the different types of fluid, demonstrating both the influence of elasticity and shear-thinning on the velocity fields. In addition, the impact of the wall on the flow field is qualitatively different for each fluid. We find that the flow behind the sphere is strongly dependent on the fluid properties as well as the elasticity. Also, the presence of a negative wake is found for the shear-thinning fluid at high Weissenberg number (Wi > 1).     

Effect of poly(ethylene oxide) on the rheological behavior of silica suspensions

The steady-shear viscosity, dynamic viscoelasticity, and sedimentation behavior were measured for silica suspensions dispersed in aqueous solutions of poly(ethylene oxide) (PEO). For suspensions prepared with polymer solutions in which the transient network is developed by entanglements, the viscosity at a given shear rate decreases, shows a minimum, and then increases with increasing particle concentration. Because the suspensions are sterically stabilized under the conditions where the particle surfaces are fully covered with by a thick layer of adsorbed polymer, the viscosity decrease can be attributed to the reduction of network density in solution. But under the low coverage conditions, the particles are flocculated by bridging and this leads to a viscosity increase with shear-thinning profiles. The polymer chains with high molecular weights form flexible bridges between particles. The stress-dependent curve of storage modulus measured by a stress amplitude sweep shows an increase prior to a drastic drop due to structural breakdown. The increase in elastic responses may arise from the restoring forces of extended bridges with high deformability. The effect of PEO on the rheological behavior of silica suspensions can be explained by a combination of concentration reduction of polymer in solution and flocculation by bridging.     

Delamination mechanism maps for a strong elastic coating on an elastic–plastic substrate subjected to contact loading

Hard wear resistant coatings that are subjected to contact loading sometimes fail because the coating delaminates from the substrate. In this report, systematic finite element computations are used to model coating delamination under contact loading. The coating and substrate are idealized as elastic and elastic–plastic solids, respectively. The interface between coating and substrate is represented using a cohesive zone law, which can be characterized by its strength and fracture toughness. The system is loaded by an axisymmetric, frictionless spherical indenter. We observe two failure modes: shear cracks may nucleate just outside the contact area if the indentation depth or load exceeds a critical value; in addition, tensile cracks may nucleate at the center of the contact when the indenter is subsequently removed from the surface. Delamination mechanism maps are constructed which show the critical indentation depth and force required to initiate both shear and tensile cracks, as functions of relevant material properties. The fictitious viscosity technique for avoiding convergence problems in finite element simulations of crack nucleation and growth on cohesive interfaces allows us to explore a wider parametric space that a conventional cohesive model cannot handle. Numerical results have also been compared to analytical analyses of asymptotic limits using plate bending and membrane stretching theories, thus providing guidelines for interpreting the simulation results.     

Viscoelastic flow of polymer solutions around a periodic,linear array of cylinders: comparisons of predictions for microstructure and flow fields

Numerical simulation is used to investigate the flow of polymer solutions around a periodic, linear array of cylinders by using three constitutive equations derived from kinetic theory of dilute polymer solutions: the Giesekus model; the finitely extensible, nonlinear elastic dumbbell model with Peterlin's approximation (FENE-P); and the FENE dumbbell model of Chilcott–Rallison (CR). In the Giesekus model, intramolecular forces are described by a Hookean spring, whereas a finitely extensible spring whose modulus is given by the Warner approximation is used in both the FENE-P and CR models. Hydro dynamic drag on the beads is taken to be anisotropic for the Giesekus model and isotropic for the other two models. The CR and FENE-P models differ subtly in their approximate treatment of the nonlinear force law. The three models exhibit very similar rheological behavior in viscometric flow and steady elongational flow, with the notable exception that the viscosity for the CR model is shear-rate independent. Finite element simulations are performed by using two different formulations: the elastic-viscous split-stress gradient (EVSS-G) method and a new variant of this formulation, the discrete EVSS-G (DEVSS-G) formulation, in which the elliptic stabilization term is added only to the discrete version of the momentum equation, and the constitutive equation is solved directly in terms of the polymer contribution to the stress tensor. Calculations are performed for all models up to a Weissenberg number We, where the configuration tensor 〈QQ〉 loses positive definiteness. However, by locally refining the mesh in the gap region, the positive definiteness of 〈QQ〉 is recovered. The flow and stress fields predicted by the three constitutive equations are qualitatively similar. A `birefringent strand' of highly stretched polymer molecules, which appears to emanate from the rear stagnation point in the cylinder, strengthens as We is increased. Not surprisingly, the molecular extension computed for the Giesekus model is considerably larger than that of the two FENE spring models. The drag force on the cylinders differs for the FENE-P and CR models, because of the difference in the shear-thinning viscosity resulting from the different approximations used in these models.     

Extracting extensional properties through excess pressure drop estimation in axisymmetric contraction and expansion flows for constant shear viscosity,extension strain-hardening fluids

In this study, hyperbolic contraction–expansion flow (HCF) devices have been investigated with the specific aim of devising new experimental measuring systems for extensional rheological properties. To this end, a hyperbolic contraction–expansion configuration has been designed to minimize the influence of shear in the flow. Experiments have been conducted using well-characterized model fluids, alongside simulations using a viscoelastic White–Metzner/FENE-CR model and finite element/finite volume analysis. Here, the application of appropriate rheological models to reproduce quantitative pressure drop predictions for constant shear viscosity fluids has been investigated, in order to extract the relevant extensional properties for the various test fluids in question. Accordingly, experimental evaluation of the hyperbolic contraction–expansion configuration has shown rising corrected pressure drops with increasing elastic behaviour (De=0～16), evidence which has been corroborated through numerical prediction. Moreover, theoretical to predicted solution correspondence has been established between extensional viscosity and first normal stress difference. This leads to a practical means to measure extensional viscosity for elastic fluids, obtained through the derived pressure drop data in these HCF devices.     

Rheological properties of oil paints and their flow instabilities in blade coating

Painters express their ideas and emotions by making use of oil paints, tools, and painting techniques. Their artworks appear on the surface, which depend on colors, tools, and the sensuous touch of the painters. From the scientific viewpoint, it is determined by the rheological properties of oil paints and their flow behavior while the external force is applied to oil paints. It means that the artists know through learning or by instinct that the oil paints have different rheological properties according to the color, and that various surfaces can be generated by changing the operating conditions. In this study, we measured the rheological properties of oil paints and investigated their flow behavior when the flow was applied. Flow instability was visualized, and the coating window was analyzed, which was supposed to mimic the painting process of the artists. We could understand that oil paints are similar in ingredients, rheological properties, and flow characteristics to industrial coating process. One important thing to note is that the painting is a process to pursue heterogeneity that is generated by the touches of the artists, while the coating process in industry is a process to pursue homogeneity so as to get defect-free surface with high productivity.     

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.     

Rheological properties of disperse systems at low shear stresses

Three-dimensional network structures can be built up in disperse systems due to long-range colloidal interactions between the dispersed particles. The rheological behaviour of such coagulation structures has been studied by means of creep and recovery experiments at low shear stresses, i.e. by measuring the shear strain as a function of time under constant stress and after removal of stress. Measurements of this type give insight into the elastic and viscous deformations and the retardation times necessary to reach equilibrium or steady-state conditions.Results obtained with dispersions of pigments in polymer solutions and with monodisperse polymer latexes indicate the existence of an equilibrium state at low shear stresses with a predominant elastic deformation and a high viscosity suggesting that the disperse systems investigated do not behave exactly as rigid gels but apparently exhibit a dynamic equilibrium of structural break-down and formation under applied stress. This behaviour is approximately described by a 4-parameter-model with an instantaneous and a steady-state compliance, one retardation time, and a viscosity.At higher shear stresses thixotropic structural break-down occurs resulting in a transition from the rheological behaviour described here to a liquid-like state with a comparatively low viscosity. In this stress range the viscoelastic properties become strongly time-dependent.These measurements give evidence of the presence of two types of deformation: an instantaneous, purely elastic deformation attributable to the unperturbed coagulation structure and the creep-recovery behaviour of an elastic liquid apparently related to the breaking and re-forming of bonds.     

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau rheological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (ρ₁/ρ_g up to 10³) and high viscosity (η₁/η_g up to 10⁴). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids.     

Rheology and birefringence of Fomblin YR at very high shear rates

The rheological and structural properties of perfluoropolyether (PFPE) lubricant films including viscosity, shear stress, and birefringence were measured at relatively low to extremely high shear rates using a rotational optical rheometer. The viscosity of various films with different thicknesses exhibit Newtonian behavior up to a shear rate 1 × 10⁴ s⁻¹, with a transition to shear-thinning behavior obvious at higher shear rates. Birefringence of these films was also measured for the first time, and these results indicate chain alignment with shear in the shear-thinning regime. The shear rate at which alignment occurs is similar to that of the onset of shear thinning. This correlation between chain alignment and shear thinning provides direct evidence that the ability of PFPEs to lubricate hard drives at high shear rates is a direct consequence of the ability of the applied shear field to align the molecules on a molecular level.     

Leveling of thixotropic liquids

In this paper the effect of thixotropy in the hydrodynamic behavior of thin films is studied. The simple problem of leveling on a horizontal substrate is considered. The rheological properties of the material are assumed to evolve over time due purely to changes in its internal structure. These changes are modeled in terms of a single structural variable. Neither elastic nor yielding effects are taken into account. More specifically, two distinct rheological models are considered: the simple model proposed by Moore and the more complex model proposed by Baravanian et al. These models exhibit a large range of variation for the liquid viscosity across the film thickness. After deriving the hydrodynamic equations governing leveling flows with the standard assumptions required by the lubrication approximation and running time-dependent numerical simulations, the nonlinear leveling history of the liquid can be predicted as a function of the initial microstructural state, rheological parameters, and initial disturbance of the liquid free surface. The main effort of this work is devoted to devising approximation schemes which lead to significant simplifications of the governing equations and their numerical computations. By approximating the inverse of viscosity as a monotonic function between its substrate and free-surface values, excellent agreement is found for the film amplitude, irrespective of the values of the rheological parameters of both models. Finally, a linear analysis yields a generalization of the Orchard’s law of leveling for Newtonian liquids to take into account the effect of thixotropy.     

Operability windows in viscoelastic slot coating flows using a simplified viscoelastic-capillary model

The viscoelastic-capillary model to predict approximately coating windows for the stable operations of viscoelastic coating liquids is derived using a lubrication approximation in slot coating processes. Pressure distributions and velocity profiles for viscoelastic liquids based on the Oldroyd-B and Phan-Thien and Tanner (PTT) models are solved in the coating bead region considering the Couette-Poiseuille flow feature and the pressure jumps at upstream and downstream menisci. Practical operating limits for the uniform coating of rheologically different liquids that are free from leaking and bead break-up defects are constructed under various conditions, incorporating the position of the upstream meniscus as an important indicator while determining limits. The shift of the uniform operating range shows different patterns for the Oldroyd-B liquid with a constant shear viscosity and the PTT liquid with a shear-thinning nature in comparison with the Newtonian case. The windows predicted by the simplified model are corroborated with experimental observations for one Newtonian and two viscoelastic liquids.     

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.