A numerical study of the flow of Bingham-like fluids in two-dimensional vane and cylinder rheometers using a smoothed particle hydrodynamics (SPH) based method

In this paper, a Lagrangian formulation of the Navier–Stokes equations, based on the smoothed particle hydrodynamics (SPH) approach, was applied to determine how well rheological parameters such as plastic viscosity can be determined from vane rheometer measurements. First, to validate this approach, a Bingham/Papanastasiou constitutive model was implemented into the SPH model and tests comparing simulation results to well established theoretical predictions were conducted. Numerical simulations for the flow of fluids in vane and coaxial cylinder rheometers were then performed. A comparison to experimental data was also made to verify the application of the SPH method in realistic flow geometries. Finally, results are presented from a parametric study of the flow of Bingham fluids with different yield stresses under various applied angular velocities of the outer cylindrical wall in the vane and coaxial cylinder rheometers. The stress, strain rate and velocity profiles, especially in the vicinity of the vane blades, were computed. By comparing the calculated stress and flow fields between the two rheometers, the validity of the assumption that the vane could be approximated as a cylinder for measuring the rheological properties of Bingham fluids at different shear rates was tested.     

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

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

The yield surface of viscoelastic and plastic fluids in a vane viscometer

Vane viscometers are often used to investigate the low shear rate properties of plastic fluids. The shear stress is determined by assuming that the material is held in the space between the vane blades so that it behaves like a rigid cylinder. Experimental evidence supports this assumption and the aim of the present study is to model numerically the yield process in a vane rheometer using viscoelastic and plastic fluids. The finite element method has been used to model the behavior of Herschel-Bulkley (Bingham), Casson and viscoelastic (Maxwell type) fluids. The penalty function approach for the pressure approximation and a rotating reference frame are used together with fine meshes containing more than 1300 elements. The results show that for Herschel-Bulkley (Bingham), and Casson fluids a rotating rigid cylinder of fluid is trapped inside the periphery of the vane, the shear stress is uniformly distributed over the surface of the cylinder. Finally a modified second order fluid is used to simulate the viscoelastic behaviour, anticipated to be an intermediate between the elastic deformation and the plastic flow, to provide a more realistic simulation of the yield process about a vane. In this case, contrast with the concentration of the elastic strain rate at the blade tips, a nearly uniform distribution of the plastic shear rate is still found. This implies that the plastic shear always distributes uniformly during the entire yielding process. Evidently the assumption of uniform shear on a rotating cylinder of material occluded in the blades of a vane is a valid and useful model for many types of fluid possessing a yield stress.     

Rheometry on magnetorheological (MR) fluids

Test fixtures of a commercial concentric cylinder rheometer (Physica Rheolab MC 20) were modified to enable measurements under magnetic inductions up to 0.5 Tesla in a shear rate range of 0.1 up to 1000 s ^–1 and temperatures 0° to 150°C. In the 2 x90°-cups only two 90° sectors of the stationary part of the double concentric cylinder arrangement are submitted to the magnetic field which is created outside the test tools by an electromagnet. A prototype of a 360°-cup contains the electromagnet within the cup and avoids the correction necessary for the sector geometry. Measurements are shown for a carbonyl iron MR fluid and two nano MR fluids. An encouraging comparison of the viscosity function and MR effect (shear stress changes due to the field) measured by using the various cups is presented. The detailed investigation of the magnetic field distribution in the tools yields a distinct radial field gradient and also stray fields that make the quantification of the effective field in the gap difficult. The change of the field when the gap is filled with MR fluid is addressed. MR effects up to 13 000 Pa have been found, the limited torque range of the rheometer making it necessary to use relatively small gap dimensions which introduce errors due to edge effects. Shear rates up to 40000 s^–1 as typical for the application in dampers were investigated by a piston-driven capillary rheometer making use of a thermostated rectangular slit with superimposed magnetic field. A satisfactory agreement of the magnetorheological data with the concentric cylinder results is found in the overlapping shear rate range.     

Optimization of the vane geometry

The use of nonstandard geometries like the vane is essential to measure the rheological characteristics of complex fluids such as non-Newtonian fluids or particle dispersions. For this geometry which is of Couette type, there is no analytical simple model defining the relation between the shear stress and the torque or relating the angular velocity to the shear rate. This study consists on calibrating a nonstandard vane geometry using a finite volume method with the Ansys Fluent software. The influence of geometrical parameters and rheological characteristics of the complex fluids are considered. First, the Newtonian fluid flow in a rotative vane geometry was simulated and a parametric model is derived therefrom. The results show an excellent agreement between the calculated torque and the measured one. They provide the possibility to define equivalent dimensions by reference to a standard geometry with concentric cylinders where the relationships between shear stress (resp. shear rate) and the torque (resp. the angular rotation) are classical. Non-Newtonian fluid flows obeying a power law rheology with different indices were then simulated. The results of these numerical simulations are in very good agreement with the preceding Newtonian-based model in some ranges of indices. The absolute difference still under 5 % provided the index is below 0.45. Finally, this study provides a calibration protocol in order to use nonstandard vane geometries with various heights, gaps, and distance to the cup bottom for measuring the rheology of complex fluids like shear thinning fluids and concentrated suspensions.     

Development of a new mixing rheometer for studying rheological behaviour of concentrated energetic suspensions

The overall objective is to present a procedure based on a Couette analogy to quantitatively analyse torque/rotor speed data and extract viscosity/shear-rate curves using a non-conventional geometry. Diphasic flows of energetic concentrated suspensions of melt-cast insensitive explosives exhibit particular rheological properties. The characterization of these complex fluids may be a challenging task when conventional rheometers are used. Placing these dense suspensions in a classic cylindrical geometry may lead to a partial destruction of the internal fluid structure. To prevent that, a “RheoXF” a mixer-type rheometer has been developed: it consists of a mixing device with quite a complex geometry rotating in a cylindrical tank. To evaluate the rheological constants (virtual radius, virtual shear rate and stress constants) of thus mixing rheometer, we used five Newtonian fluids. After this calibration, the rheological characterizations were carried out on five formulations. The unique parameter which changes in these formulations is the batch's origin of a secondary explosive: the 3-nitro-1,2,4-triazole-5-one. These energetic particles differ by their morphology, maximum packing density and may be by their process synthesis. After having determined pseudoplastic parameters, a correlation has been made with the evolution of maximum packing density values calculated with De Larrard model.     

A double wall-ring geometry for interfacial shear rheometry

The rheological properties of complex fluid interfaces are of prime importance in a number of technological and biological applications. Whereas several methods have been proposed to measure the surface rheological properties, it remains an intrinsically challenging problem due to the small forces and torques involved and due to the intricate coupling between interfacial and bulk flows. In the present work, a double wall-ring geometry to measure the viscoelastic properties of interfaces in shear flows is presented. The geometry can be used in combination with a modern rotational rheometer. A numerical analysis of the flow field as a function of the surface viscoelastic properties is presented to evaluate the non-linearities in the surface velocity profile at a low Boussinesq number. The sensitivity of the geometry, as well as its applicability, are demonstrated using some reference Newtonian and viscoelastic fluids. Oscillatory and steady shear measurements on these reference complex fluid interfaces demonstrate the intrinsic sensitivity, the accuracy, and the dynamic range of the geometry when used in combination with a sensitive rheometer.     

Numerical investigation of plate edge and slot size effects in low yield stress measurements with a slotted plate device

There is a need for accurate yield stress measurements, especially in the case of low yield stress complex materials such as biological samples. This task cannot be accomplished with conventional rotational rheometers due to significant wall slip effects and the necessity to operate the device at very low shear rates, often beyond the limit that such rheometers can achieve. In this paper, we focused on the slotted plate method proposed recently for low yield stress measurements. Using computational fluid dynamics, we studied the effects of plate geometry on the measurement accuracy of the slotted plate method. Results of this study indicate that both wall slip effects and pressure drag force can be substantially reduced by adopting a thin plate with sharp front and rear edges, high slot area ratio, and large number of slots. If the plate has 30° triangular edges, a slot area ratio of 80%, and 12 slots, the slotted plate method overpredicts the yield stress of a 0.09 wt.% Carbopol dispersion (yield stress of 9.17 Pa) by only 8.4% under no-slip conditions and underpredicts the yield stress by 12.3% under free-slip conditions. Similar results were obtained for human saliva characterized by a very low yield stress (0.073 Pa).     

Rheometry for large-particulated fluids: analysis of the ball measuring system and comparison to debris flow rheometry

For large-particulated fluids encountered in natural debris flow, building materials, and sewage treatment, only a few rheometers exist that allow the determination of yield stress and viscosity. In the present investigation, we focus on the rheometrical analysis of the ball measuring system as a suitable tool to measure the rheology of particulated fluids up to grain sizes of 10 mm. The ball measuring system consists of a sphere that is dragged through a sample volume of approximately 0.5 l. Implemented in a standard rheometer, torques exerted on the sphere and the corresponding rotational speeds are recorded within a wide measuring range. In the second part of this investigation, six rheometric devices to determine flow curve and yield stress of fluids containing large particles with maximum grain sizes of 1 to 25 mm are compared, considering both rheological data and application in practical use. The large-scale rheometer of Coussot and Piau, the building material learning viscometer of Wallevik and Gjorv, and the ball measuring system were used for the flow curve determination and a capillary rheometer, the inclined plane test, and the slump test were used for the yield stress determination. For different coarse and concentrated sediment–water mixtures, the flow curves and the yield stresses agree well, except for the capillary rheometer, which exhibits much larger yield stress values. Differences are also noted in the measuring range of the different devices, as well as for the required sample volume that is crucial for application.     

A modified rheometer for in-situ radial and tangential SAXS studies on shear-induced alignment

A tooth rheometer designed to investigate the influence of large amplitude oscillatory shear on the macroscopic orientation of complex fluids is described. The rheometer allows in-situ two-directional small-angle scattering measurements, making it possible to observe in particular the orientation process of hierarchically structured block copolymers in dynamics. Experiments performed with a linear poly(vinyl methyl ether) show that, despite a considerable gap/width ratio, the rheometer can be used for adequate rheological measurements. In addition we found an unexpected dependence of the dynamic parameters on the rheometer gap. The small amount of sample required favors the use of this rheometer for measuring the rheological characteristics of new compounds.     

Rheology of pulp suspensions using ultrasonic Doppler velocimetry

Conventional rheometry coupled with local velocity measurements (ultrasonic Doppler velocimetry) are used to study the flow behaviour of various commercial pulp fibre suspensions at fibre mass concentrations ranging from 1 to 5 wt.%. Experimental data obtained using a stress-controlled rheometer by implementing a vane in large cup geometry exhibits apparent yield stress values which are lower than those predicted before mainly due to existence of apparent slip. Pulp suspensions exhibit shear-thinning behaviour up to a high shear rate value after which Newtonian behaviour prevails. Local velocity measurements prove the existence of significant wall slippage at the vane surface. The velocimetry technique is also used to study the influence of pH and lignin content on the flow behaviour of pulp suspensions. The Herschel–Bulkley constitutive equation is used to fit the local steady-state velocity profiles and to predict the steady-state flow curves obtained by conventional rheometry. Consistency between the various sets of data is found for all suspensions studied, including apparent yield stress, apparent wall slip and complete flow curves.     

Flow instability and shear localization in a drilling mud

To have a better knowledge of problems occurring with drilling fluids in complex wells, we carried out a detailed rheological analysis of a typical drilling mud at low shear rates using both conventional rheometry and MRI velocimetry. We show the existence of a viscosity bifurcation effect: Below a critical stress value, the mud tends to completely stop flowing, whereas beyond this critical stress, it reaches an apparent shear rate larger than a finite (critical) value, and no stable flows can be obtained between this critical shear rate value and zero. These results are confirmed by MRI velocity profiles, which exhibit a slope break at the interface between the solid and the liquid phases inside the Couette geometry. Moreover, this viscosity bifurcation is a transient phenomenon, the progressive development of which can be observed by MRI. A further examination of MRI data shows that, in the transient regime, the shear rate does not vary monotonously in the rheometer gap and is particularly large along the outer (rough) cylinder, which might be at the origin of the development of a region of constant shear rate in the apparent flow curve.     

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.     

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.     

Rotational rheometry of liquid metal systems: Measurement geometry selection and flow curve analysis

In the present study, a rotational measurement technique was used to evaluate viscosities of liquid metals and metallic alloys. Three types of measurement geometries in a high temperature rotational rheometer were evaluated: cone and plate, DIN coaxial, and double concentric cylinder (DCC). The DCC geometry proved to be the most effective. An analytical solution has been presented to evaluate the viscosity as a function of shear rate for DCC geometry. The flow curves and shear viscosities of pure Al, pure Zn and Sn_95.8Ag_3.28Cu_0.92 solder alloy have been evaluated as a function of shear rate and melt superheat temperature. It is proposed that liquid metal systems are non-Newtonian and strongly shear thinning in flow behavior.     

Yield stress of structured fluids measured by squeeze flow

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

Electrorheological response of SnO<Subscript>2</Subscript> and Y<Subscript>2</Subscript>O<Subscript>3</Subscript> nanoparticles in silicon oil

The electrorheological properties (ER) of some fluids containing particles change extensively under the external electrical field. This phenomenon is applicable in many industries and equipments, such as clutches and motor driven rotor, which would transfer the spin to a drive shaft through a thin layer of electrorheological fluid. In this investigation, the effects of external electrical field on ER properties of non-Newtonian fluids (silicon oil) with the addition of SnO₂ and Y₂O₃ nanoparticles were studied. The ER properties were measured for a wide range of SnO₂ and Y₂O₃ nanoparticle concentrations and DC electrical voltages using concentric cylinder rotary rheometer. Based on the results, ER properties of nanofluids, e.g., apparent viscosity, shear stress, and yield stress, were enhanced by applying electrical field and increasing SnO₂ and Y₂O₃ concentrations.

An investigation of squeeze flow as a viable technique for determining the yield stress

A systematic study of squeeze flow (SF) was performed on different concentrations of Carbopol with varying yield stresses. A sample of constant volume was placed between two parallel plates and a series of constant force steps applied, following the plate separation as a function of time. Precise rheological measurements of the model yield stress fluids were performed in addition to the well-controlled SF tests. These rheological measurements were used in conjunction with the SF equations to determine the time-dependent plate separation, allowing a direct comparison of theory and experiment throughout the entire test. The limiting height achieved during constant force SF reveals information about the yield stress of the fluid as predicted by the theory. It appears that by carefully controlling the experimental conditions of the squeeze test one can obtain yield stress values that agree with the rheological measurements within 10%. Additionally, the validity of the lubricational theory was tested; not only for the determination of the yield stress but throughout the flow as well.     

Fluid inertia effects on a Controlled Stress Rheometer in its oscillatory mode

This paper is concerned with the effect of fluid inertia on experimental oscillatory stress data taken from a Controlled Stress Rheometer. A linear viscoelastic theory is developed which includes the effect of fluid inertia for cone and plate, parallel plate and concentric cylinder geometries. This theory is used to interpret dynamic data for both a slightly elastic and a highly viscoelastic fluid. It is shown that intertial effects are very small for both a cone and plate and parallel plate geometries. Inertial effects, however, can be important in the concentric cylinder geometry.     

Study of magnetorheological fluids at high shear rates

The tunable rheological properties of magnetorheological (MR) materials at high shear rates are studied using a piston-driven flow-mode-type rheometer. The proposed method provides measurement of the apparent viscosity and yield stress of MR fluids for a shear rate range of 50 to 40,000 s⁻¹. The rheological properties of a commercial MR fluid, as well as a newly developed MR polymeric gel, and a ferrofluid-based MR fluid are investigated. The results for apparent viscosity and dynamic and static shear stresses under different applied magnetic fields are reported.