The end correction for power-law fluids in the capillary rheometer

The finite element scheme developed by Nickell, Tanner and Caswell is used to compute the entry and exit losses for creeping flow of power-law fluids in a capillary rheometer. The predicted entry losses for a Newtonian fluid agree well with available experimental and theoretical results. The entry losses for inelastic power-law fluids increased with decreasing flow behaviour index and show an increasing deviation from available upper bound results as the flow behaviour index in the power-law decreases.The exit losses are found to be finite for inelastic power-law fluids and increase as the flow behaviour index decreases. The predicted die swell for Newtonian fluids agrees well with the available experimental data while the influence of shear thinning is to reduce the die swell.The end correction which is the sum of the entry and exit losses relative to twice the viscometric wall shear stress varies from 0.834 for n = 1 to 2.917 for n = 1/6. This figure reaches a very high value as n tends to zero. The experimental variation in the Couette correction factor in capillary rheometry is explained in terms of the shear thinning characteristics of the fluid. It is concluded that the exit flow is not viscometric, contrary to a common assumption.     

Film thickness in blade coating of viscous and viscoelastic liquids

Coating of viscous and viscoelastic liquids is examined both theoretically and experimentally. A rigid blade, accurately positioned over a rotating roll, provides an experimental system in which coating thickness is measured as a function of geometric parameters. A perturbation solution to the Navier—Stokes equations yields a lubrication theory which shows agreement with the data to an extent depending on the specific geometry.The effect of a non-Newtonian viscosity is explored by adopting a purely viscous power-law model. The lubrication equations are solved by the method of Horowitz and Steidler [1], and predict an increase in coating thickness relative to the Newtonian case. Data for viscoelastic fluids show both an increase and a decrease in coating thickness compared with Newtonian liquids depending on the relative magnitude of shear thinning and elastic effects.     

Extrudate swell in some dilute elastic solutions

An experimental programme has been carried out to obtain data for the isothermal swelling of dilute viscoelastic fluids based on the polyisobutylene-polybutene system. The fluid was thus a constant viscosity one and so the observed swell was predominantly due to elastic effects alone. To establish the experimental rig and procedure, preliminary results for the creeping flow of a Newtonian fluid, based on a polybutene, were acquired. The Newtonian results, showing a creeping flow swell of 13.8%, demonstrated that the rig and procedure were satisfactory for the formation of free, isothermal jets. The viscoelastic results were quite unexpected, showing that dilute elastic solutions can exhibit significant swelling despite the fact that the recoverable shear is small (S _R < 0.1). These results qualitatively confirm a recent numerical study [1].     

The translation of a Newtonian droplet in a 4-constant Oldroyd fluid

The viscoelastic correction to the Hadamard—Rybczynski terminal velocity was calculated for the retarded motion translation of a finite viscosity Newtonian droplets in a 4-constant Oldroyd fluid. The droplet translational velocity can be enhanced or hindered relative to the Hadamard—Rybczinski value depending upon the degree of shear thinning, and elongational and memory effects in the viscoelastic fluid. Specifically, it is found that for a large droplet viscosity, for which surface mobility is negligible, shear thinning and fluid memory enhance translational motion. At the other extreme, when the droplet surface is very mobile (i.e. for gas bubbles) the translational velocity can be enhanced or hindered relative to the Hadamard—Rybczynski value depending on the relative influence of the memory and the elongational properties of the viscoelastic phase.     

Non-Newtonian flow in microporous structures under the electroviscous effect

Single phase non-Newtonian microporous flow combined with the electroviscous effect is investigated in the pore-scale under conditions of various rheological properties and electrokinetic parameters. The lattice Boltzmann method is employed to solve both the electric potential field and flow velocity field. The simulation of commonly used power-law non-Newtonian flow shows that the electroviscous effect on the flow depends on both the fluid rheological behavior and pore surface area ratio significantly. For the shear thinning fluid with power-law exponent n < 1, the fluid viscosity near the wall is smaller and the electrovicous effect plays a more important role compared to the Newtonian fluid and shear thickening fluid. The high pore surface area ratio in the porous structure increases the electroviscous force and the induced flow resistance becomes important even to the flow of Newtonian and shear thickening fluids.     

A note on the translation of a thin rod inside a cylinder

The problem of a thin rod moving longitudinally along the axis of symmetry of a cylindrical vessel is examined for Newtonian and non-Newtonian liquids. For non-Newtonian fluids, the inelastic power-law type solution predicts the experimental results particularly well. On account of wall effects, the induced pressure gradients are much greater for a Newtonian fluid than for a viscoelastic fluid. In fact, in the latter case, they may be considered negligible when the radius of the inner cylinder is small compared to the one of the outer cylinder.     

Highly elastic solutions for Oldroyd-B and Phan-Thien/Tanner fluids with a finite volume/element method: planar contraction flows

A hybrid finite volume/element method is analysed through the computation of creeping flows of viscoelastic fluids in plane 4:1 sharp and rounded-corner contraction geometries. Simulations are presented for three models: a constant viscosity Oldroyd-B fluid, and Phan-Thien/Tanner (PTT) shear thinning fluids of exponential and linear approximation form. A Taylor–Galerkin/pressure-correction scheme is implemented as the base time-stepping framework. The momentum equations are solved by a finite element method, whilst the constitutive equations are solved by a finite volume approach. Mesh convergence is analysed via refinement around the contraction to capture boundary layers and flow structure. Pressure drop is shown to increase with flow rate for a fixed fluid. For the Oldroyd-B model, singular behaviour is reported in the main stress component as one approaches the corner in the rounded, as with the sharp geometry. Velocity components display an asymptotic trend with a positive slope. Higher values of Weissenberg numbers (We) are reached with these finite volume schemes compared to their finite element counterparts, attributing this to superior accuracy properties.     

Viscous dissipation with fluid inertia in oscillatory shear flow

For liquids with high viscosity and low thermal conductivity, viscous dissipation can cause appreciable errors in rheological property measurements. Here, the influences of both viscous dissipation and fluid inertia on the property measurements in oscillatory sliding plate rheometry are investigated. For Newtonian fluids, Bird (1965) solved the combined problem analytically, but only for high frequencies. Here his solution is extended to any frequencies. Also, the equations of motion and energy are solved for linear viscoelastic fluids, and new analytical solutions for the velocity and temperature profiles are given. In both Newtonian and linear viscoelastic fluids, the temperature rise in the gap increases with frequency. The location of the maximum temperature shifts from the mid-plane at low frequency towards the moving wall at high frequency. The fluid inertia increases the viscous dissipation in both fluids. By solving the combined problem, this paper simplifies rheometer design by providing one unified criterion for avoiding measurement errors. Operating limits are presented graphically for minimizing the effects of both fluid inertia and viscous dissipation in oscillatory sliding plate rheometry.     

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.     

Dispersion of linear gravity waves on a viscoelastic fluid in an horizontal canal

A characteristic equation is derived that describes the spatial decay of linear surface gravity waves on Maxwell fluids. Except at small frequencies, the derived dispersion relation is different from the temporal decay dispersion relation which is normally studied within fluid mechanics. The implications for waves on viscous Newtonian fluids using the two different dispersion relations is briefly discussed. The wave number is measured experimentally as function of the frequency in a horizontal canal. Seven Newtonian fluids and four viscoelastic liquids with constant viscosity have been used in the experiments. The spatial decay theory for Newtonian fluids fits well to the experimental data. The model and experiments are used to determine limits for the Maxwell fluid time numbers for the four viscoelastic liquids. As a result of low viscosity it was not possible within this study to obtain these time numbers from oscillatory experiments. Therefore, a comparison of surface gravity wave experiments with theory is applicable as a method to evaluate memory times of low viscosity viscoelastic fluids.     

A parametric study of droplet deformation through a microfluidic contraction: Shear thinning liquids

Numerical simulations of a droplet passing through an axisymmetric microfluidic contraction are presented, focusing on systems where one of the two liquids present is shear thinning. The simulations are performed using a transient Volume of Fluid (VOF) algorithm. When the droplet is shear thinning and the surrounding phase Newtonian, droplets deform in a similar way to Newtonian droplets that have a viscosity equal to the average viscosity of the shear thinning fluid while it is within the contraction. When the surrounding phase is shear thinning and the droplet Newtonian, droplets deform in a similar way to droplets contained within a Newtonian liquid that has a viscosity that is lower than that of the droplet. In both cases the behaviour of the shear thinning fluid can be broadly described in terms of a ‘characteristic’ Newtonian viscosity: However, determining the exact value of this viscosity without performing a full shear thinning simulation is not possible.     

Turbulent channel flow of dilute polymeric solutions: Drag reduction scaling and an eddy viscosity model

Direct numerical simulation of viscoelastic turbulent channel flows up to the maximum drag reduction (MDR) limit has been performed. The simulation results in turn have been used to develop relationships between the flow and fluid rheological parameters, i.e. maximum chain extensibility, Reynolds number, Re_τ, and Weissenberg number, We_τ and percent drag reduction (%DR) as well as the slope increment of the mean velocity profile. Moreover, based on the trends observed in the mean velocity profile and the overall momentum balance three different regimes of drag reduction (DR), namely, low drag reduction (LDR; 0 ≤ %DR ≤ 20), high drag reduction (HDR; 20 ≤ %DR ≤ 52) and MDR (52 ≤ %DR ≤ 74) have been identified and mathematical expressions for the eddy viscosity in these regimes are presented. It is found that both in LDR and HDR regimes the eddy viscosity varies with the distance from the channel wall. However, in the MDR regime the ratio of the eddy viscosity to the Newtonian one tends to a very small value around 0.1 within the channel. Based on these expressions a procedure that relies on the DNS predictions of the budgets of momentum and viscoelastic shear stress is developed for evaluating the mean velocity profile.     

The roles of inertia and shear-thinning in flow of an inelastic liquid through an axisymmetric sudden contraction

The roles of luid inertia and shear-rate dependent viscosity in determining the velocity field in an axisymmetric sudden contraction are assessed by finite-element analysis for a generalized Newtonian fluid with viscosity function given by a Carreau equation. Acting alone, either increasing shear-thinning of the viscosity or increasing fluid inertia suppresses the upstream vortex that surrounds the opening to the small tube. For creeping flows, shear thinning does not produce concavities and off-centered maxima in the axial velocity profile just inside the small tube, even at high Carreau numbers where the velocity field approaches the limiting form for a power-law fluid. Peaks in the axial velocity away from the center of the tube were found only for moderate and high Reynolds numbers and were enhanced by shear thinning, which decreased the viscosity and consequently increased the “local” Reynolds number near the wall of the small tube. The effect of steep velocity gradients near this surface on the accuracy of the finite-element approximations is discussed.     

Newtonian,power law,and infinite shear flow characteristics of concentrated slurries using percolation theory concepts

There is a strong interest today in concentrated particulate-filled dispersion and slurries in both polymeric and Newtonian fluids. This paper reviews and extends theoretical approaches using percolation theory concepts to characterize the rheological behavior of fluids filled with particulate solids. First, a previously proposed limiting, zero shear viscosity model based on percolation theory concepts is reviewed. This model has been primarily tested with rigid fillers in a Newtonian carrier and polymeric fluids. Second, all Newtonian fluid-based slurries that have a high concentration of filler become pseudoplastic, shear-thinning slurries at some threshold shear rate. A new theory is reviewed and new data are evaluated that correlate the power law constant, n, to cluster formation of the fillers suspended in the fluids in shear flow. Slurry systems reported here cover a size range from 58 nm to 200 μm. Third, this cluster percolation-based rheological analysis is then extended to a newly proposed model for the calculation of the ratio of infinite shear, η_∞, to the zero shear viscosity, η₀. Using literature data, it is demonstrated that measurements of the viscosity ratio, η_∞/η₀, correlate with the power law through the use of an energy dissipation-based model for Bingham rheological fluids.     

Performance modeling and analysis of blood flow in elastic arteries

Nomenclature　　τ　　wallshearstressγshearrateτy yieldstressηc Cassonviscosityktheconsistencyindexnnon_Newtonianindexτp shearstressofthepthelementωangularvelocityRvessel’sradiusCwavespeedM　　magneticparameter (Hartmannnumber)u,w　velocitycomponentinther_andz_directions,respectivelyP　　pressureα　　unsteadinessparameter k , R meanparametersTp relaxationtimeofthepthelementρ densityIntroductionTheimportancetoatherogenesisofarterialflowphenomenasuchasflowseparation ,recirculationands…     

Development of a Low-Reynolds-number k-ω Model for FENE-P Fluids

A low-Reynolds-number k-ω model for Newtonian fluids has been developed to predict drag reduction of viscoelastic fluids described by the FENE-P model. The model is an extension to viscoelastic fluids of the model for Newtonian fluids developed by Bredberg et al. (Int J Heat Fluid Flow 23:731–743, 2002). The performance of the model was assessed using results from direct numerical simulations for fully developed turbulent channel flow of FENE-P fluids. It should only be used for drag reductions of up to 50 % (low and intermediate drag reductions), because of the limiting assumption of turbulence isotropy leading to an under-prediction of k, but compares favourably with results from k-ε models in the literature based on turbulence isotropy.     

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.     

GPU accelerated simulations of three-dimensional flow of power-law fluids in a driven cube

Newtonian fluid flow in two- and three-dimensional cavities with a moving wall has been studied extensively in a number of previous works. However, relatively a fewer number of studies have considered the motion of non-Newtonian fluids such as shear thinning and shear thickening power law fluids. In this paper, we have simulated the three-dimensional, non-Newtonian flow of a power law fluid in a cubic cavity driven by shear from the top wall. We have used an in-house developed fractional step code, implemented on a Graphics Processor Unit. Three Reynolds numbers have been studied with power law index set to 0.5, 1.0 and 1.5. The flow patterns, viscosity distributions and velocity profiles are presented for Reynolds numbers of 100, 400 and 1000. All three Reynolds numbers are found to yield steady state flows. Tabulated values of velocity are given for the nine cases studied, including the Newtonian cases.     

A FENE-P k–ε turbulence model for low and intermediate regimes of polymer-induced drag reduction

A new low-Reynolds-number k–ε turbulence model is developed for flows of viscoelastic fluids described by the finitely extensible nonlinear elastic rheological constitutive equation with Peterlin approximation (FENE-P model). The model is validated against direct numerical simulations in the low and intermediate drag reduction (DR) regimes (DR up to 50%). The results obtained represent an improvement over the low DR model of Pinho et al. (2008) [A low Reynolds number k–ε turbulence model for FENE-P viscoelastic fluids, Journal of Non-Newtonian Fluid Mechanics, 154, 89–108]. In extending the range of application to higher values of drag reduction, three main improvements were incorporated: a modified eddy viscosity closure, the inclusion of direct viscoelastic contributions into the transport equations for turbulent kinetic energy (k) and its dissipation rate, and a new closure for the cross-correlations between the fluctuating components of the polymer conformation and rate of strain tensors (NLT_ij). The NLT_ij appears in the Reynolds-averaged evolution equation for the conformation tensor (RACE), which is required to calculate the average polymer stress, and in the viscoelastic stress work in the transport equation of k. It is shown that the predictions of mean velocity, turbulent kinetic energy, its rate of dissipation by the Newtonian solvent, conformation tensor and polymer and Reynolds shear stresses are improved compared to those obtained from the earlier model.