Electric effects on the rheology of insulating oils in electrodes with flocked fabric

The rheological behavior of insulating oils is studied in nonuniform electric fields which are generated by an electrode covered with flocked fabric. Although the oils show no electrorheological effects in uniform fields between metal electrodes with smooth surfaces, the flocked fabric leads to a striking increase of viscosity in steady shear. The viscosity enhancement increases with decreasing zero-field viscosity and decreasing conductivity of oils. In the limit of zero shear rate, the oils with low conductivity behave as solids with yield stress. When a very small quantity of fine particles is introduced into electrified oils without shear, a rapid and large-scale motion of particles is observed between the tips of fibers and the plate electrode. The local motion of fluids in high electric fields is referred to as electrohydrodynamic (EHD) convection. Periodic patterns of circulation flow are formed in static oils. The electric energy which is dissipated during the circulation motion contributes to holding the periodic flow in static oils. When the stress is very low, the periodic patterns are not broken down. The yield stress corresponds to the force required to rupture the domain structures of EHD convection. In shear fields, the additional energy may be required to change the periodic patterns of EHD convection. The striking increase of viscosity in steady shear can be attributed to the interactions between EHD convection and external shear. Received: 31 August 1998 Accepted: 17 February 1999     

Electrorheology of dielectric liquids

Dielectric liquids that show striking electrorheological (ER) effects are formulated by controlling the conductivity. Although the viscosity is increased on the application of a d.c. field, the flow of electrified fluids is Newtonian in the plain electrodes with smooth surfaces. When the liquids are sandwiched between the electrodes with flocked fabrics, the viscosity behavior is converted from Newtonian to shear-thinning flow. In electric fields, the convective flow is induced over the system due to the electrohydrodynamic(EHD) effect. The interactions between EHD convection and external shear give rise to the additional energy dissipation and in turn the increase in viscosity. The ER effects of simple liquids are very attractive in application to new fluid devices.     

A “zero-parameter” constitutive relation for simple shear viscoelasticity

Based on the Cox-Merz rule and Eyring's expression for the nonlinear shear viscosity, a Wagner-type constitutive relation with no nontrivial adjustable parameters is proposed for simple shear viscoelasticity. The predictions for a number of non-steady shear flows are worked out analytically. It is shown that most features of shear viscoelasticity are reproduced by the model.     

Viscoelasticity of a dielectric fluid in nonuniform electric fields generated by electrodes with flocked fabrics

The rheological behavior of a dielectric fluid is studied in nonuniform electric fields which are generated by an electrode covered with flocked fabrics. Although no electrorheological (ER) effects are observed in uniform fields between metal electrodes with smooth surfaces, striking increases in viscosity and elastic response are induced by the electrode with flocked fabrics. The presence of flocked fabrics does not have a significant effect on the fluid rheology without electric fields. The ER behavior and current density are influenced by the fiber length even at a constant field strength. When a very small amount of fine particles is introduced in the electrified fluid without shear, we can see the rapid and large-scale motion of particles between the tips of fibers and plate electrode. In high DC fields, the Coulomb force acting on a free charge often gives rise to the secondary motion of fluid. The local motion of fluid in high electric fields is refereed to as electrohydrodynamic (EHD) convection. The additional energy may be required to change the periodic patterns of EHD convection by forced shear. Therefore, the ER effect demonstrated by the modification of electrode with flocked fabrics can be attributed to a combined effect of EHD convection and external shear. Received: 10 March 1998 Accepted: 1 June 1998     

Shear-thickening flow of suspensions of carbon nanofibers in aqueous PVA solutions

The suspensions of carbon nanofibers in aqueous poly(vinyl alcohol) solutions were prepared in the presence of spherical carbon black particles, and the steady-shear viscosity and dynamic viscoelasticity were measured for complex suspensions. Although the single suspensions of carbon black are highly stable, the flocculation of carbon nanofibers is promoted by the addition of carbon black particles. The complex suspensions show remarkable shear thickening in the steady-flow and strain hardening in oscillatory shear with large amplitude. The nonlinear responses strongly depend on the carbon black concentration, whereas the dynamic viscoelasticity at low strains in the linear ranges is not significantly influenced. As the highly elastic effects arise from the long-range motion of particles, the possible mechanism may be the orientation of nanofibers in strong shear fields. The suspensions show the time-dependent behavior of viscosity when the time-scale of measurements is shorter than that of structural recovery to the isotropic states.     

Rheology of an emulsion of viscoelastic drops in steady shear

Steady shear rheology of a dilute emulsion with viscoelastic inclusions is numerically investigated using direct numerical simulations. Batchelor's formulation for rheology of a viscous emulsion is extended for a viscoelastic system. Viscoelasticity is modeled using the Oldroyd-B constitutive equation. A front-tracking finite difference code is used to numerically determine the drop shape, and solve for the velocity and stress fields. The effective stress of the viscoelastic emulsion has three different components due to interfacial tension, viscosity difference (not considered here) and the drop phase viscoelasticity. The interfacial contributions – first and second normal stress differences and shear stresses – vary with Capillary number in a manner similar to those of a Newtonian system. However the shear viscosity decreases with viscoelasticity at low Capillary numbers, and increases at high Capillary numbers. The first normal stress difference due to interfacial contribution decreases with increasing drop phase viscoelasticity. The first normal stress difference due to the drop phase viscoelasticity is found to have a complex dependence on Capillary and Deborah numbers, in contrast with the linear mixing rule. Drop phase viscoelasticity does not contribute significantly to effective shear viscosity of the emulsion. The total first normal stress difference shows an increase with drop phase viscoelasticity at high Capillary numbers. However at low Capillary numbers, a non-monotonic behavior is observed. The results are explained by examining the stress field and the drop shape.     

Surfactant spreading on a thin weakly viscoelastic film

A mathematical model is presented for surfactant-driven thin weakly viscoelastic film flows on a flat, impermeable plane. The Oldroyd-B constitutive relation is used to model the viscoelastic fluid. Lubrication theory and a perturbation expansion in powers of the Weissenberg number (We) are employed, which give rise to non-linear coupled evolution equations governing the transport of insoluble surfactant and thin liquid film thickness. Spreading on a Newtonian film is recovered to leading order and corrections to viscoelasticity are obtained at order We. These equations are solved numerically over a wide range of viscosity ratio (ratio of solvent viscosity to the sum of solvent and polymeric viscosities), pre-existing surfactant level and Peclet number (Pe). The effect of viscoelasticity on surfactant transport and fluid flow is investigated and the mechanisms underlying this effect are explored. Shear stress, streamwise normal stress and the temporal rate of change of extra shear stress generated from gradients in surfactant concentration dominate thin viscoelastic film flows whereas only shear stresses play a role in Newtonian thin film flows. Our results also reveal that, for weak viscoelasticity, the influence of viscosity ratio on the evolution of surfactant concentration and film thickness can be significant and varies considerably, depending on the concentration of pre-existing surfactant and surfactant surface diffusivity.     

Nonlinear dependence of viscosity in modeling the rate-dependent response of natural and high damping rubbers in compression and shear: Experimental identification and numerical verification

The rate-dependent behavior of filled natural rubber (NR) and high damping rubber (HDR) is investigated in compression and shear regimes. In order to describe the viscosity-induced rate-dependent effects, a constitutive model of finite strain viscoelasticity founded on the basis of the multiplicative decomposition of the deformation gradient tensor into elastic and inelastic parts is proposed. The total stress is decomposed into an equilibrium stress and a viscosity-induced overstress by following the concept of the Zener model. To identify the constitutive equation for the viscosity from direct experimental observations, an analytical scheme that ascertains the fundamental relation between the inelastic strain rate and the overstress tensor of the Mandel type by evaluating simple relaxation test results is proposed. Evaluation of the experimental results using the proposed analytical scheme confirms the necessity of considering both the current overstress and the current deformation as variables to describe the evolution of the rate-dependent phenomena. Based on this experimentally based motivation, an evolution equation using power laws is proposed to represent the effects of internal variables on viscosity phenomena. The proposed evolution equation has been incorporated in the finite strain viscoelasticity model in a thermodynamically consistent way. Simulation results for simple relaxation tests, multi-step relaxation tests and monotonic tests at different strain rates using the developed model show an encouraging correlation with the experiments conducted on HDR and NR in both compression and shear regimes. Finally, an approach to extend the proposed evolution equation for rate-dependent cyclic processes is proposed. The simulation results are critically compared with the experiments.     

Interaction between Taylor bubbles rising in stagnant non-Newtonian fluids

The interaction between Taylor bubbles rising in stagnant non-Newtonian solutions was studied. Aqueous solutions of carboxymethylcellulose (CMC) and polyacrylamide (PAA) polymers were used to study the effect of different rheological properties: shear viscosity and viscoelasticity. The solutions studied covered a range of Reynolds numbers between 10 and 714, and Deborah numbers up to 14. The study was performed with pairs of Taylor bubbles rising in a vertical column (0.032 m internal diameter) filled with stagnant liquid. The velocities of the leading and trailing bubbles were measured by sets of laser diodes/photocells placed along the column. The velocity of the trailing bubble was analysed together with the liquid velocity profile in the wake of a single rising bubble (Particle Image Velocimetry data obtained from the literature). For the less concentrated CMC solutions, with moderate shear viscosity and low viscoelasticity, the interaction between Taylor bubbles was similar to that found in Newtonian fluids. For the most concentrated CMC solution, which has high shear viscosity and moderate viscoelasticity, a negative wake forms behind the Taylor bubbles, inhibiting coalescence since the bubbles maintain a minimum distance of about 1D between them. For the PAA solutions, with moderate shear viscosity but higher viscoelasticity than the CMC solutions, longer wake lengths are seen, which are responsible for trailing bubble acceleration at greater distances from the leading bubble. Also in the PAA solutions, the long time needed for the fluid to recover its initial shear viscosity after the passage of the first bubble makes the fluid less resistant to the trailing bubble flow. Hence, the trailing bubble can travel at a higher velocity than the leading bubble, even at distances above 90D.     

Bimodal model of slurry viscosity with application to coal-slurries. Part 1. Theory and experiment

The rheological behavior of stable slurries is shown to be characterized by a bimodal model that represents a slurry as made up of a coarse fraction and a colloidal size fine fraction. According to the model, the two fractions behave independently of each other, and the non-Newtonian behavior of the viscosity is solely caused by the colloidal fraction, while the coarse fraction increases the viscosity level through hydrodynamic interactions. Data from experiments run with colloidal coal particles of about 2–3 µm average size dispersed in water show the viscosity of these colloidal suspensions to exhibit a highly shearrate-dependent behavior and, in the high shear limit, to match very closely the viscosity of suspensions of uniform size rigid spheres although the coal volume fraction must be determined semi-empirically. Different amounts of coarse coal particles are added to the colloidal suspension and the viscosity of the truly bimodal slurries measured as a function of shear rate. In agreement with the bimodal model, the measured shear viscosities show the coarse fraction to behave independently of the colloidal fraction and its contribution to the viscosity rise to be independent of the shear rate. It is shown that the shear rate exerted on the colloidal fraction is higher than that applied by the viscometer as a result of hydrodynamic interactions between the coarse particles, and that it is this effective higher shear rate which is necessary to apply in the correlations. For determining the coal volume fraction a relatively simple and quite accurate measurement technique is developed for determining the density and void fraction of coarse porous particles; the technique directly relates volume fraction to mass fraction.     

Drag reduction in the turbulent pipe flow of polymers

The paper concerns an experimental study of the fully developed turbulent pipe flow of several different aqueous polymer solutions: 0.25%, 0.3% and 0.4% carboxymethylcellulose (CMC), 0.2% xanthan gum (XG), a 0.09%/0.09% CMC/XG blend, 0.125% and 0.2% polyacrylamide (PAA). The flow data include friction factor vs. Reynolds number, mean velocity and near-wall shear rate distributions, and axial velocity fluctuation intensity u′ at a fixed radial location as a laminar/turbulent transition indicator. For each fluid we also include measurements of shear viscosity, first normal-stress difference and extensional viscosity. At high shear rates we find that the degree of viscoelasticity increases with concentration (0.3% CMC is an exception) for a given polymer, and in the sequence XG, CMC/XG, CMC, PAA, whilst at low shear rates the ranking changes to CMC, CMC/XG, XG, PAA. The extensional viscosity ranking is XG/CMC, XG, CMC, PAA at high strain rates and the same as that for the viscoelasticity at low shear rates. We find that the observed drag-reduction behaviour is consistent for most part with the viscoelastic and extensional-viscosity behaviour at the low shear and strain rates typical of those occurring in the outer zone of the buffer region.Although laminar/turbulent transition is practically indiscernible from the friction factor vs. Reynolds number plots, particularly for PAA and XG, the u′ level provides a very clear indicator and it is found that the transition delay follows much the same trend with elasticity/extensional viscosity as the drag reduction.     

Electrorheology of dilute suspensions induced by hydrodynamic instability

The viscosity behavior in electric fields was measured for dilute suspensions of p-[perfluoro(2-isopropyl-1,3-dimethyl-1-butenyl)oxy]benzoic acid particles (PFNA) in silicone oils. The application of electric fields causes a viscosity increase in a wide range of shear rates. Since the electrorheological (ER) effect is much stronger at low shear rates, the flow becomes shear-thinning. However, contrary to conventional ER suspensions which are reversibly converted between Newtonian fluids and Bingham solids, the PFNA suspensions are fluids even in electric fields. When the particle concentration is increased to 5 wt.%, the ER effect reaches saturation. Further increase does not contribute to additional viscosity enhancement. These results cannot be explained through the chain formation mechanism established for conventional systems. After the ER experiments, the bob surface of the rheometer is covered with several stripes of aggregated particles. Although the strength of electric and shear fields is constant in the rheometer, the periodic structure may be formed in the flow of electrified suspensions. When a dielectric liquid is subjected to high electric fields, the secondary motion of liquid can be induced by the Coulomb force acting on free charge. The electrohydrodynamic (EHD) convection is responsible for the periodic distribution of particles concentration. The ER effect of PFNA suspensions may be generated by a combined effect of EHD convection and external shear.     

Viscosity of suspensions in polymeric solutions

Viscosities of suspended particles in polymeric solutions depend upon dissolved polymer concentration, volume fraction of particles and shear rate. In this analysis of viscosity data, relative viscosity was defined as the ratio of suspension viscosity to solution viscosity at the same shear stress rather than shear rate. These relative viscosities reached asymptotic values at high shear stress for all concentrations of dissolved polymer and for all particle loadings. At a given particle loading, the asymptotic values of relative viscosity were nearly independent of the concentration of dissolved polymer. Realtive viscosities were correlated with volume fraction by the one-constant equation of Maron and Pierce.     

Effect of variations in counterion to surfactant ratio on rheology and microstructures of drag reducing cationic surfactant systems

Rheology, drag reduction and cryo-TEM experiments were performed on Arquad 16–50/NaSal and Ethoquad O/12/NaSal surfactant systems at different counterion-to-surfactant ratios and at constant low surfactant concentrations, 5 mM, appropriate for drag reduction. The molar ratio of counterion-to-surface was varied from 0.6 to 2.5. All the surfactant systems described here are viscoelastic and drag reducing. The viscoelasticity and drag reducing effectiveness increase with increase in counterion/surfactant ratio. Network are present in the solutions with high ratio, and they are viscoelastic. However, shear is needed to induce network formation for solutions at low ratio. Cryo-TEM images confirm the existence of thread-like micelles which form entanglement networks, and show that the micellar network becomes denser with increasing counterion/surfactant ratio in one surfactant series. Both increase in the counterion/surfactant ratio and increase in the shear rate result in shorter relaxation times. For some of these systems, abrupt increase in viscosity is observed at certain shear rates which are time effects affecting microstructure rearrangements rather than formation of shear induced structures.     

Influence of viscoelasticity on drop deformation and orientation in shear flow: Part 1. Stationary states

The influence of matrix and droplet viscoelasticity on the steady deformation and orientation of a single droplet subjected to simple shear is investigated microscopically. Experimental data are obtained in the velocity–vorticity and velocity–velocity gradient plane. A constant viscosity Boger fluid is used, as well as a shear-thinning viscoelastic fluid. These materials are described by means of an Oldroyd-B, Giesekus, Ellis, or multi-mode Giesekus constitutive equation. The drop-to-matrix viscosity ratio is 1.5. The numerical simulations in 3D are performed with a volume-of-fluid algorithm and focus on capillary numbers 0.15 and 0.35. In the case of a viscoelastic matrix, viscoelastic stress fields, computed at varying Deborah numbers, show maxima slightly above the drop tip at the back and below the tip at the front. At both capillary numbers, the simulations with the Oldroyd-B constitutive equation predict the experimentally observed phenomena that matrix viscoelasticity significantly suppresses droplet deformation and promotes droplet orientation. These two effects saturate experimentally at high Deborah numbers. Experimentally, the high Deborah numbers are achieved by decreasing the droplet radius with other parameters unchanged. At the higher capillary and Deborah numbers, the use of the Giesekus model with a small amount of shear-thinning dampens the stationary state deformation slightly and increases the angle of orientation. Droplet viscoelasticity on the other hand hardly affects the steady droplet deformation and orientation, both experimentally and numerically, even at moderate to high capillary and Deborah numbers.     

The influence of fluid properties on electrohydrodynamic heat transfer enhancement in liquids under viscous and electrically dominated flow conditions

Fluid property effects on electrohydrodynamic (EHD) heat transfer enhancement were investigated. Heat transfer, pressure drop, electrical power requirements, and the transition between the viscous dominated and electrically dominated flow regimes as a function of fluid properties were examined using three cooling oils having widely varying physical properties. Low viscosity and low electrical conductivity gave the greatest heat transfer enhancement for a given electrical power input. The required electrical power to achieve a specified heat transfer enhancement was greater for working fluids that had a small charge relaxation time, defined as the ratio of the electrical permittivity to the electrical conductivity. These results correlate well with available experimental and analytical data. A theoretical prediction of the effect of fluid properties and forced flow rate on the onset of EHD enhancement was experimentally verified. The onset of significant EHD heat transfer enhancement occurs most readily in low viscosity liquids at low Reynolds number flows for a given electrical power input.     

A numerical investigation of electrohydrodynamic (EHD) effects on bubble deformation under pseudo-nucleate boiling conditions

In this article, the electrohydrodynamic (EHD) effects on nucleate boiling are studied by developing a numerical modelling of EHD effect on bubble deformation in pseudo-nucleate boiling conditions. The volume of fluid (VOF) method is employed to track the interface between the gas–liquid two phases; the user-defined code is written and added to the commercial software FLUENT to solve the electric field and the corresponding electric body force. On this basis, the model is applied to study the EHD effects on heat transfer and fluid flows. An initial air bubble surrounded by liquid CCl₄ and attached to a horizontal superheated wall under the action of electric field is studied. The results of the EHD effect on bubble shape evolution are compared with those of available experiments showing good agreement. The mechanism of EHD enhancement of heat transfer and the EHD induced phenomena including bubble elongation and detachment are analyzed in detail.     

Discontinuous shear thickening in the presence of polymers adsorbed on the surface of calcium carbonate particles

In the presence of dispersant molecules currently used in cement industry and based on polyethylene oxide (PEO), we found a strong discontinuous shear thickening (DST) at high volume fraction in suspensions of calcium carbonate particles. The transition was reversible and the critical shear rate and shear stress for which this instability appears are reported versus the volume fraction of particles. A model of repulsive forces between polymers, taking into account the thickness of the polymer layer and the density of adsorption on the surface of the particles, can explain the differences of critical stresses observed between these three dispersant molecules. In particular, it explains why a small polymer densely adsorbed can be more efficient to repel the transition at higher stress than a larger molecule less densely adsorbed. Above the transition, we find that the suspension presents a special kind of stick-slip instability with even the presence of a negative shear rate under constant applied stress. A model is proposed which well predicts this regime by taking into account both the inertia of the apparatus and the viscoelasticity of the suspension.     

On uniaxial extension of an elasto-viscous cylinder

Summary This paper deals with slow processes of extension of an elasto-visoous cylinder under conditions where the theory of linear viscoelasticity is applicable. Kinematic dependences are given and four nonstationary problems are solved concerning the extension of a cylinder when one of the following parameters is constant: strain rate, extension rate, stress, and tensile force.Experiments were performed on a low-molecular polyisobutylene at 25° and constant extension rates. The changes in total deformation and stresses, and after unloading, in highly elastic and irreversible deformation, were determined. The effect of the forces of surface tension on the elastic recovery of the samples was taken into account. The formulas of the linear theory of viscoelasticity are shown to apply within the range of extension rates studied. The viscosity and rubbery elasticity modulus were determined and found to agree well with the results obtained during steady shear flows.     

Flow of granular materials modeled as a non-linear fluid

Recently, Massoudi (2011a) derived a generalized form of a constitutive relation related to Reiner's fluid model for wet sand, where not only the effects of volume fraction are incorporated in the rheological properties of the fluid, but also the shear viscosity depends on the shear rate. In this paper, we use this model to study the fully developed flow of granular-like materials down an inclined plane. The governing equations are made dimensionless and numerical solutions are presented for the various dimensionless parameters.