Nonlinear rheology of concentrated spherical silica suspensions: 3. Concentration dependence

Nonlinear rheology was examined for concentrated suspensions of spherical silica particles (with radius of 40 nm) in viscous media, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture and pure ethylene glycol. The particles were randomly and isotropically dispersed in the media in the quiescent state, and their effective volume fraction φ_eff ranged from 0.36 to 0.59. For small strains, the particles exhibited linear relaxation of the Brownian stress σ_B due to their diffusion. For large step strains γ, the nonlinear relaxation modulus G(t,γ) exhibited strong damping and obeyed the time-strain separability. This damping was related to γ-insensitivity of strain-induced anisotropy in the particle distribution that resulted in decreases of σ_B/γ. The damping became stronger for larger φ_eff. This φ_eff dependence was related to a hard-core volume effect, i.e., strain-induced collision of the particles that is enhanced for larger φ_eff. Under steady/transient shear flow, the particles exhibited thinning and thickening at low and high γ˙, respectively. The thinning behavior was well described by a BKZ constitutive equation using the G(t,γ) data and attributable to decreases of a Brownian contribution, σ_B/γ˙. The thickening behavior, not described by this equation, was related to dynamic clustering of the particles and corresponding enhancement of the hydrodynamic stress at high γ˙. In this thickening regime, the viscosity growth η₊ after start-up of flow was scaled with a strain γ˙t. Specifically, critical strains γ_d and γ_s for the onset of thickening and achievement of the steadily thickened state were independent of γ˙ but decreased with increasing φ_eff. This φ_eff dependence was again related to the hard-core volume effect, flow-induced collision of the particles enhanced for larger φ_eff. Received: 26 June 1998 Accepted: 9 December 1998     

Viscosity of concentrated noncolloidal bidisperse suspensions

At the same solid volume fraction (Φ) the relative viscosity (η _r ) of a concentrated noncolloidal bidisperse suspension of hard spherical particles is lower than that of a monodisperse suspension. In this paper a semi-analytical viscosity model of noncolloidal bidisperse suspensions is derived using an integration method. In this model the random loose packing density obtained by computer simulation is taken as the limit of solid volume fraction Φ _m which depends upon both the diameter ratio (λ) of large to small particles and the volume fraction of large particles (ξ=Φ _l /Φ). This model shows that at high solid volume fraction, Φ > 0.40, both λ and ξ significantly influence η _r . For example, at Φ=0.5, it predicts that for monodisperse suspensions η _r =70, while for bidisperse suspensions (λ=2 and ξ=0.7) η _r =40. Comparison shows that, at high solid volume fraction (0.4–0.5), the relative viscosity predicted by this model is in good agreement with that predicted by the work of Shapiro and Probstein (1992) and of Patlazhan (1993), but is higher than that predicted by the work of others. Received: 27 February 2001 Accepted: 25 April 2001     

An empirical model predicting the viscosity of highly concentrated, bimodal dispersions with colloidal interactions

The relationship between particle size distribution and viscosity of concentrated dispersions is of great industrial importance, since it is the key to get high solids dispersions or suspensions. The problem is treated here experimentally as well as theoretically for the special case of strongly interacting colloidal particles. An empirical model based on a generalized Quemada equation is used to describe η as a function of volume fraction for mono- as well as multimodal dispersions. The pre-factor η˜ accounts for the shear rate dependence of η and does not affect the shape of the η vs φ curves. It is shown here for the first time that colloidal interactions do not show up in the maximum packing parameter and φ_max can be calculated from the particle size distribution without further knowledge of the interactions among the suspended particles. On the other hand, the exponent ɛ is controlled by the interactions among the particles. Starting from a limiting value of 2 for non-interacting either colloidal or non-colloidal particles, ɛ generally increases strongly with decreasing particle size. For a given particle system it then can be expressed as a function of the number average particle diameter. As a consequence, the viscosity of bimodal dispersions varies not only with the size ratio of large to small particles, but also depends on the absolute particle size going through a minimum as the size ratio increases. Furthermore, the well-known viscosity minimum for bimodal dispersions with volumetric mixing ratios of around 30/70 of small to large particles is shown to vanish if colloidal interactions contribute significantly. Received: 7 June 2000/Accepted: 12 February 2001     

Rheology of carbon black suspensions. II. Well dispersed system

Linear viscoelastic properties were investigated for the suspensions of carbon black (CB) particles having covalently-fused aggregate structures of an average diameter a=120 nm. The suspending medium, an alkyd resin (AR), had a high affinity toward the CB particles, and the aggregates of these particles were well dispersed to form no higher-order agglomerates. Consequently, the suspensions obeyed the time-temperature superposition and their Arrhenius-type activation energy was identical to that of the medium. From comparison of the zero shear viscosity η₀ for the CB suspensions and hard-sphere silica suspensions, an effective volume fraction φ_eff of the CB particles was found to be 2.7 times larger than the bare volume fraction of the particles. The CB particles exhibited a slow relaxation process, and the terminal relaxation time of this process was close to the Peclet time (Brownian diffusion time) evaluated from the aggregate size a and a high frequency viscosity. Furthermore, the terminal relaxation mode distribution of the CB suspensions was well scaled with an intensity factor H_t that was evaluated from the φ_eff in a way utilized for the hard-sphere silica suspensions. These results demonstrated that the slow relaxation in the CB suspensions was dominated by the Brownian diffusion of the CB aggregates, as similar to the situation for the silica suspensions.     

Non-isobaric Marangoni boundary layer flow for Cu,Al<Subscript>2</Subscript>O<Subscript>3</Subscript> and TiO<Subscript>2</Subscript> nanoparticles in a water based fluid

In this paper, a non-isobaric Marangoni boundary layer flow that can be formed along the interface of immiscible nanofluids in surface driven flows due to an imposed temperature gradient, is considered. The solution is determined using a similarity solution for both the momentum and energy equations and assuming developing boundary layer flow along the interface of the immiscible nanofluids. The resulting system of nonlinear ordinary differential equations is solved numerically using the shooting method along with the Runge-Kutta-Fehlberg method. Numerical results are obtained for the interface velocity, the surface temperature gradient as well as the velocity and temperature profiles for some values of the governing parameters, namely the nanoparticle volume fraction φ (0≤φ≤0.2) and the constant exponent β. Three different types of nanoparticles, namely Cu, Al₂O₃ and TiO₂ are considered by using water-based fluid with Prandtl number Pr =6.2. It was found that nanoparticles with low thermal conductivity, TiO₂, have better enhancement on heat transfer compared to Al₂O₃ and Cu. The results also indicate that dual solutions exist when β<0.5. The paper complements also the work by Golia and Viviani (Meccanica 21:200–204, 1986) concerning the dual solutions in the case of adverse pressure gradient.     

Linear viscoelastic behavior of perfluorooctyl sulfonate micelles: effects of counter-ions

Linear viscoelastic behavior was investigated for aqueous solutions of perfluorooctyl sulfonate (C₈F₁₇SO⁻ ₃; abbreviated as FOS) micelles having a mixture of tetraethylammonium (N⁺(C₂H₅)₄; TEA) and lithium (Li⁺) ions as the counter-ions. The solutions had the same FOS concentration (0.1 mol l⁻¹) and various Li⁺ fractions in the counter-ions, φ_Li = 0−0.6, and the FOS micelles in these solutions formed threads which further organized into dendritic networks. At T ≤ 15 °C, the terminal relaxation time τ and the viscosity η, governed by thermal scission of the networks, increased with increasing φ_Li up to 0.55. A further increase of φ_Li resulted in decreases of τ and η and in broadening of the relaxation mode distribution. These rheological changes are discussed in relation to the role of TEA ions in thermal scission: Previous NMR studies revealed that only a fraction of TEA ions were tightly bound to the FOS micellar surfaces and these bound ions stabilized the thread/network structures. The concentration of non-bound TEA ions, C_TEA ^*, decreased and finally vanished on increasing φ_Li up to φ_Li ^* ≅ 0.6, and the concentration of the bound TEA ions significantly decreased on a further increase of φ_Li. The non-bound TEA ions appeared to catalyze the thermal scission of the FOS threads, and the observed increases of τ and η for φ_Li < 0.55 were attributed to the decrease of C_TEA ^*. On the other hand, the decreases of τ and η as well as the broadening of the mode distribution, found for φ_Li > 0.55 (where C_TEA ^* ≅ 0), were related to destabilization of the FOS threads/networks due to a shortage of the bound TEA ions and to the existence of concentrated Li⁺ ions. Viscoelastic data of pure FOSTEA and FOSTEA/FOSLi/TEACl solutions lent support to these arguments for the role of TEA ions in the relaxation of FOSTEA/FOSLi solutions. Received: 12 October 1999/Accepted: 1 November 1999     

Recoverable strains and retardation times of monodisperse suspensions of silicon dioxide spheres in poly(dimethylsiloxane)

Steady-state viscosities η, steady-state recoverable strains γ _rs and characteristic retardation time τ _1/2 were measured for suspensions of monodisperse silicon dioxide (SiO₂) spheres in poly(dimethylsiloxane) (PDMS) with various volume fractions Φ of the suspended spheres at various creep stresses σ ₀. Two different regions are found in plots of η/η _m vs γ _rs, where η/η _m denotes the relative viscosity of the suspensions. In one region, η/η _m is proportional to γ _rs, while γ _rs is independent of η/η _m in the other region. In both regions, τ _1/2 is the functions of the shear strain rate in the steady-state of creep test independently of Φ. The origin of the elasticity is related to the ‘maximally distorted’ cages recovered owing to the repulsive interaction between the SiO₂ spheres and recovery of the cages in the shear-induced clusters of the suspended spheres.     

The role of the mathematical mean value theorem (MVT) in rheometry: an easy way to convert apparent flow curves into correct ones

The mean value theorem of integral calculus guarantees that the apparent viscosity η _a can easily be converted into the correct viscosity η. For ordinary liquids there is a direct identity between η _a and η but the apparent shear rate (or apparent shear stress) has to be shifted to the representative shear rate γ˙^ (or representative shear stress τ^). A model free approximation scheme is introduced which implies a constant shift factor. The corresponding approximation for η is acceptable for liquids most commonly encountered. For plastic fluids the relation between η and η _a is more complex since it involves a function depending upon α; the yield stress relative to the maximum stress within the viscometer. Using the same approximation scheme as before the shift factor will involve α as well. The corresponding approximation of η is shown to be acceptable for the whole range of α. Received: 7 February 2000/Accepted: 15 February 2000     

Edge fracture in cone-plate and parallel plate flows

Edge fracture is an instability of cone-plate and parallel plate flows of viscoelastic liquids and suspensions, characterised by the formation of a `crack' or indentation at a critical shear rate on the free surface of the liquid. A study is undertaken of the theoretical, experimental and computational aspects of edge fracture. The Tanner-Keentok theory of edge fracture in second-order liquids is re-examined and is approximately extended to cover the Criminale-Ericksen-Filbey (CEF) model. The second-order theory shows that the stress distribution on the semi-circular crack is not constant, requiring an average to be taken of the stress; this affects the proportionality constant, K in the edge fracture equation −N _2c = KΓ/a, where N _2c is the critical second normal stress difference, Γ is the surface tension coefficient and a is the fracture diameter. When the minimum stress is used, K = 2/3 as found by Tanner and Keentok (1983). Consideration is given to the sources of experimental error, including secondary flow and slip (wall effect). The effect of inertia on edge fracture is derived. A video camera was used to record the inception and development of edge fracture in four viscoelastic liquids and two suspensions. The recorded image was then measured to obtain the fracture diameter. The edge fracture phenomenon was examined to find its dependence on the physical dimensions of the flow (i.e. parallel plate gap or cone angle), on the surface tension coefficient, on the critical shear rate and on the critical second normal stress difference. The critical second normal stress difference was found to depend on the surface tension coefficient and the fracture diameter, as shown by the theory of Tanner and Keentok (1983); however, the experimental data were best fitted by the equation −N _2c = 1.095Γ/a. It was found that edge fracture in viscoelastic liquids depends on the Reynolds number, which is in good agreement with the inertial theory of edge fracture. Edge fracture in lubricating grease and toothpaste is broadly consistent with the CEF model of edge fracture. A finite volume method program was used to simulate the flow of a viscoelastic liquid, obeying the modified Phan-Thien-Tanner model, to obtain the velocity and stress distribution in parallel plate flow in three dimensions. Stress concentrations of the second normal stress difference (N ₂) were found in the plane of the crack; the velocity distribution shows a secondary flow tending to aid crack formation if N ₂ is negative, and a secondary flow tending to suppress crack formation if N ₂ is positive. Received: 4 January 1999 Accepted: 19 May 1999     

Influence of molecular parameters on the stress dependence of viscous and elastic properties of polypropylene melts in shear

The stress dependencies of the steady-state viscosity η and, particularly, that of the steady-state elastic compliance J _e of various linear isotactic polypropylenes (PP) and one long-chain branched PP are investigated using creep-recovery tests. The creep stresses applied range from 2 to 10,000 Pa. In order to discuss the stress-dependent viscosity η and elastic compliance J _e with respect to the influence of the weight average molar mass M _w and the polydispersity factor M _w/M _n the PP are characterized by SEC–MALLS. For the linear PP, linear steady-state elastic compliances J_e⁰J_{\rm e}^0 in the range of 10^− 5–10^− 3 Pa^− 1 are obtained depending on the molar mass distribution. J_e⁰J_{\rm e}^0 of the LCB-PP is distinctly higher and comes to lie at around 10^− 2 Pa^− 1. J_e⁰J_{\rm e}^0 is found to be independent of M _w but strongly dependent on polydispersity. η and J _e decrease with increasing stress. For the linear PP, J _e as a function of the stress τ is temperature independent. The higher M _w/M _n the stronger is the shear thinning of η and the more pronounced is the stress dependence of J _e. For the LCB-PP, the strongest stress dependence of η and J _e is observed. Furthermore, for all PP J _e reacts more sensitively to an increasing stress than η. A qualitative explanation for the stronger stress dependence of J _e compared to η is given by analyzing the contribution of long relaxation times to the viscosity and elasticity.     

Simulation of the rheological properties of suspensions of oblate spheroidal particles in a Newtonian fluid

A simulation algorithm was developed to predict the rheological properties of oblate spheroidal suspensions. The motion of each particle is described by Jeffery’s solution, which is then modified by the interactions between the particles. The interactions are considered to be short range and are described by results from lubrication theory and by approximating locally the spheroid surface by an equivalent spherical surface. The simulation is first tested on a sphere suspension, results are compared with known experimental and numerical data, and good agreement is found. Results are then presented for suspensions of oblate spheroids of two mean aspect ratios of 0.3 and 0.2. Results for the relative viscosity η _r, normal stress differences N ₁ and N ₂ are reported and compared with the few available results on oblate particle suspensions in a hydrodynamic regime. Evolution of the orientation of the particles is also observed, and a clear alignment with the flow is found to occur after a transient period. A change of sign of N ₁ from negative to positive as the particle concentration is increased is observed. This phenomenon is more significant as the particle aspect ratio increases. It is believed to arise from a change in the suspension microstructure as the particle alignment increases.     

Thermal convection in a porous layer: Effects of anisotropy and surface boundary conditions

The theory describing the onset of convection in a homogeneous porous layer bounded above and below by isothermal surfaces is extended to consider an upper boundary which is partly permeable. The general boundary condition p + λ ∂p/∂n = constant is applied at the top surface and the flow is investigated for various λ in the range 0 ⩽ λ < ∞. Estimates of the magnitude and horizontal distribution of the vertical mass and heat fluxes at the surface, the horizontally-averaged heat flux (Nusselt number) and the fraction of the fluid which recirculates within the layer are found for slightly supercritical conditions. Comparisons are made with the two limiting cases λ → ∞, where the surface is completely impermeable, and λ = 0, where the surface is at constant pressure. Also studied are the effects of anisotropy in permeability, ξ = K _H /K _V , and anisotropy is thermal conductivity, η = k _H /k _V , both parameters being ratios of horizontal to vertical quantities. Quantitative results are given for a wide variety of the parameters λ, ξ and η. In the limit ξ/η → 0 there is no recirculation, all fluid being converted out of the top surface, while in the limit ξ/η → ∞ there is full recirculation.     

Rheological behavior of cellulose/monohydrate of n-methylmorpholine n-oxide solutions Part 1: Liquid state

Steady-state and dynamic experiments have been performed on solutions containing cellulose dissolved in monohydrate of N-methylmorpholine N-oxide (NMMO). The dependence of the zero-shear viscosity η₀, and of the terminal relaxation time τ _c , on concentration, average degree of polymerization (DP) and temperature are discussed. The behavior of this semi-rigid, polymolecular polymer in solution differs from that of flexible monodisperse ones. The slope of the plot of log(η₀) versus, on the one hand, log(c) at fixed molecular weight (DP)=600, and, on the other hand, log(DP) at fixed concentration (c=5%w/w) are equal to 4.6 and 5 respectively, instead of 3.4 in the concentrated region. Experimental data for the shear modulus were fitted using the classical Doi-Edwards equation with a log normal distribution of relaxation time. This distribution is compared to the distribution of DP. Received: 25 February1997 Accepted: 30 December 1997     

Time Periodic Solutions to the One-Dimensional Nonlinear Wave Equation

This paper is concerned with the time periodic solutions to the one-dimensional nonlinear wave equation with either variable or constant coefficients. By adjusting the basis of L ² function space, we can circumvent the difficulties caused by η _u  = 0 and obtain the existence of a weak periodic solution, which was posed as an open problem by Baubu and Pavel in (Trans Am Math Soc 349:2035–2048, 1997). Finally, an application to the forced Sine-Gordon equation is presented to illustrate the utility of this technique.     

Rheological properties and dispersion stability of magnetorheological (MR) suspensions

In the present article, the rheological responses and dispersion stability of magnetorheological (MR) fluids were investigated experimentally. Suspensions of magnetite and carbonyl iron particles were prepared as model MR fluids. Under an external magnetic field (H ₀) and a steady shear flow, the yield stress depends upon H ₀ ^3/2. The Yield stress depended on the volume fraction of the particle (φ) linearly only at low concentration and increased faster at high fraction. Rheological behavior of MR fluids subjected to a small-strain oscillatory shear flow was investigated as a function of the strain amplitude, frequency, and the external magnetic field. In order to improve the stability of MR fluid, ferromagnetic Co-γ-Fe₂O₃ and CrO₂ particles were added as the stabilizing and thickening agent in the carbonyl iron suspension. Such needle-like particles seem to play a role in the steric repulsion between the relatively large carbonyl iron particles, resulting in improved stability against rapid sedimentation of dense iron particles. Furthermore, the additive-containing MR suspensions exhibited larger yield stress, especially at higher magnetic field strength. Received: 4 April 2000 Accepted: 6 November 2000     

Rheology of carbon black suspensions. III. Sol-gel transition system

Linear viscoelastic properties of carbon black (CB) suspensions with various CB volume fractions () in a rosin-modified phenol resin type varnish (Varnish-1) were investigated at various temperatures (T). The CB/Varnish-1 suspensions exhibited a sol-gel transition on an increase in , and the _gel value at the gelation point decreased with increasing T. This T dependence of _gel, being opposite to the dependence seen for usual gelling systems, can be related to a phenol resin type polymeric component included in the Varnish-1. At low T, this polymeric component appeared to be rather well solvated in the Varnish-1 thereby allowing the gelation due to bare attraction between the CB particles at large . In contrast, at high T, the polymeric component appeared to have been less solvated, as evidence from a moderate failure of the time-temperature superposition of pure Varnish-1 and a decrease of its elasticity (in a shifted frequency scale) with increasing T. This less solvated polymeric component would have been adsorbed on the CB particles, thereby allowing the agglomeration of the particles at small _gel at high T.     

Viscous properties of calcium carbonate filled polymer melts

Summary The viscous properties of calcium carbonate filled polyethylene and polystyrene melts were examined. The relative vircosity _r defined in the previous paper gave an asymtptotic value( _r)_l in the range of the shear stress below 10⁵ dyne/cm².( _r)_l of the calcium carbonate filled system was higher than that of the glass beads or glass balloons filled system at the same volume fraction of the filler. Maron-Pierce equation with ₀ = 0.44 was able to approximate the( _r)_l — relationship. However, it was deduced here that the high value of( _r)_l of calcium carbonyl filled system was due to the apparent increase of and this increase was attributed to the fixed polymer layer formed on the powder particle. By assuming the particle as a sphere with a diameter of 2 µm, the thickness of the fixed polymer layer was estimated as about 0.17 µm. The yield stress estimated from the Casson's plots increased exponentially with.

---

Zusammenfassung Es wurden die viskosen Eigenschaften von Polyäthylen-und Polystyrol-Schmelzen untersucht, die mit Kalziumkarbonat-Teilchen gefüllt waren. Für die relative Viskosität _r, wie sie in einer vorangegangenen Veröffentlichung definiert worden war, ergab sich bei Schubspannungen unterhalb 10⁵ dyn/cm² ein asymptotischer Wert( _r)_l. Dieser war bei den mit Kalziumkarbonat gefüllten Schmelzen höher als bei Schmelzen, die bis zur gleichen Volumenkonzentration mit Glaskugeln oder Glasballons gefüllt waren. Die ( _r) _l —-Abhängigkeit ließ sich durch eine Gleichung nachMaron und Pierce mit ₀ = 0,44 beschreiben. Es wurde jedoch geschlossen, daß der hohe( _r)_l-Wert der mit Kalziumkarbonat gefüllten Schmelzen auf eine scheinbare Zunahme von zurückzuführen ist, verursacht durch eine feste Polymerschicht auf der Teilchenoberfläche. Unter Annahme kugelförmiger Teilchen mit einem Durchmesser von 2 µm ließ sich die zugeordnete Schichtdicke zu 0,17 µm abschätzen. Die mittels der Casson-Beziehung geschätzte Fließspannung ergab eine exponentielle-Abhängigkeit.

---

---

With 7 figures and 1 table     

On a Crystalline Variational Problem, Part II:¶BV Regularity and Structure of Minimizers on Facets

For a nonsmooth positively one-homogeneous convex function φ:ℝ ⁿ → [0,+∞[, it is possible to introduce the class ?_φ (ℝ ⁿ ) of smooth boundaries with respect to φ, to define their φ-mean curvature κ_φ, and to prove that, for E∈?_φ (ℝ ⁿ ), κ_φ∈L ^∞(δE) [9]. Based on these results, we continue the analysis on the structure of δE and on the regularity properties of κ_φ. We prove that a facet F of δE is Lipschitz (up to negligible sets) and that κ_φ has bounded variation on F. Further properties of the jump set of κ_φ are inspected: in particular, in three space dimensions, we relate the sublevel sets of κ_φ on F to the geometry of the Wulff shape ?_φ≔{φ≤ 1 }. Accepted October 11, 2000?Published online 14 February, 2001     

Rheological analysis of highly concentrated w/o emulsions

A series of highly concentrated lipophilic cosmetic emulsions were analysed, in order to determine their rheological and textural properties, as a function of their microstructure. The originality of this study lies in the methodology used, especially the shear-stress scanning analysis. The results of a very powerful and comprehensive dynamic rheological analysis suggest the existence of two critical volume fraction values: besides the “close-packed” value φ _c , a “slack-packed” value φ₀, close to 0.60 could be demonstrated. It has been shown that the close-packed structure is stable under shear; in constrast, the slack-packed configuration, defined as φ₀<φ<φ _c is unstable under shear. A comparison with theoretical models, especially that of Princen, showed good agreement and allowed the close-packed value φ _c to be defined more precisely as 0.67. The gap between 0.67 and 0.74 is probably indicative of a highly polydisperse distribution, as confirmed by microscopic analysis. Flow experiments confirmed the validity of Princen‘s model. Received: 20 February 1997 Accepted: 20 January 1998     

Flocculation and rheology of kaolinite/quartz suspensions

Rheological properties of suspensions of Na-kaolinite and colloidal quartz (Min-U-Sil) at constant overall volume concentration of 2% are determined with a Weissenberg Rheogoniometer using a combined Couette and cone-and-plate geometry. The results are interpreted in terms of the flocculation behaviour of the constituent particles in the presence of high salt concentrations (0.1–0.75m NaCl) at pH 6, 7 and 8. In these chemical environments these suspensions are pseudoplastic for much of the range of mixture compositions becoming Newtonian for suspensions containing only quartz. These properties reflect the dominant influence of interactions between kaolinite particles on the flocculation behaviour of the mixture.Nomenclature a ₁, a₂ radii of spheres - A Hamaker constant - b radius of cylinder - C volumetric solids concentration - e ₀ electronic charge - H ₀ shortest distance between surfaces - I ionic strength - J collision frequency per unit volume - k Boltzmann constant - l length of cylinder - m number in eq. (2) - N particle number concentration - S (H ₀ + b)/b in eq. (4) - T absolute temperature - U electrophoretic mobility - V _A van der Waals attractive energy - V _R coulombic energy - V _T total energy of interaction - X H ₀/2a₁ in eq. (5) - Y a ₂/a₁ in eq. (5) - thickness of plate - shear rate - permittivity - zeta potential - k Debye-Hückel parameter - µ dynamic viscosity - µ _pl plastic viscosity - v valency of counter ion - shear stress - _B Bingham stress - ₁, ₂ dimensionless potentials - ₁, ₂ surface potentials