Unsteady-state development of plane couette flow for viscoelastic fluids

Unsteady-state development of plane Couette flow for viscoelastic fluids is analyzed using a constitutive equation that can be obtained from molecular theory, in which the molecules are regarded as finitely extensible dumbbells. Typical features of the flow situation are as follows: (i) For a fluid with moderate elasticity, not only stress overshoot but also velocity overshoot are predicted. (ii) For suitable combinations of elasticity and gap width, and for some time intervals stress propagation and reflection phenomena are predicted. (iii) After a sufficient time has elapsed, the stress state behaves similarly to that corresponding to the start-up of a steady simple shear flow.     

Squeezing a viscoelastic fluid from a cone

The paper reports an exact kinematics for the squeezing flow from a cone of a general viscoelastic fluid. To obtain numerical values for the stresses, a network model that allows stress overshoot and shear-thinning in the start-up of a shear flow is adopted. Both these features are important in this flow. For the special case of an Oldroyd-B fluid it is shown that there is a limiting Weissenberg number above which at least one component of the stresses increases unboundedly with time.     

Squeezing a viscoelastic liquid from a wedge: an exact solution

The paper reports an exact solution for the squeezing flow from a wedge of a general viscoelastic liquid. To obtain numerical values for the field variables, a network model that allows stress overshoot and shear-thinning in the start-up of a shear flow is adopted. It is found that both these features are important in this transient flow; stress overshoot is responsible for a stiffer response of the fluid (compared to the inelastic case) at moderate time —at large time, shear-thinning dominates and the fluid behaves like an inelastic fluid. On the other hand, the Oldroyd-B fluid always predicts a softer response than the Newtonian one. Furthermore, there is a limiting Weissenberg number above which one component of the stresses of the Oldroyd-B fluid increases unboundedly with time. This limiting Weissenberg number is approximately $\text{sol2}\text{3}$.     

Squeeze-film flow of a viscoelastic fluid a lubrication model

In this paper the constant-speed squeezing flow of a model viscoelastic fluid which allows for a Carreau-type viscosity and a stress overshoot in a start up of a shear flow is considered. Based on the knowledge of the full numerical solution by a Boundary Element method a set of lubrication equations is constructed which we believe to retain all the essential physics. The lubrication equations consist of four partial differential equations in time and one spatial coordinate. These are more tractable than the full set of governing equations (partial differential equations in time and two spatial coordinates). The lubrication equations are solved by a Finite Difference scheme, and the results are compared favourably with the full numerical solutions. The mechanism for load enhancement is also discussed.     

Entropy effects in viscoelastic fluid theory

A derivation of viscoelastic fluid theory is given which includes the effect of the entropy dependence upon deformation history. The theory is then specialized to model second-order effects. The resulting theory predicts the observed effects of shear thinning and stress overshoot. Comparison is made with experimental data.     

Unsteady compressible boundary layer flow over a circular cone near a plane of symmetry

An analysis has been performed to study the unsteady laminar compressible boundary layer governing the hypersonic flow over a circular cone at an angle of attack near a plane of symmetry with either inflow or outflow in the presence of suction. The flow is assumed to be steady at time t=0 and at t>0 it becomes unsteady due to the time-dependent free stream velocity which varies arbitrarily with time. The nonlinear coupled parabolic partial differential equations under boundary layer approximations have been solved by using an implicit finite-difference method. It is found that suction plays an important role in stabilising the fluid motion and in obtaining unique solution of the problem. The effect of the cross flow parameter is found to be more pronounced on the cross flow surface shear stress than on the streamwise surface shear stress and surface heat transfer. Beyond a certain value of the cross flow parameter overshoot in the cross flow velocity occurs and the magnitude of this overshoot increases with the cross flow parameter. The time variation of the streamwise surface shear stress is more significant than that of the cross flow surface shear stress and surface heat transfer. The suction and the total enthalpy at the wall exert strong influence on the streamwise and cross flow surface shear stresses and the surface heat transfer except that the effect of suction on the cross flow surface shear stress is small.     

Kinematics and normal stress differences of a viscoelastic fluid undergoing wiggle flow

An experimental program was carried out to determine the laminar regime kinematics and normal stress differences of a viscoelastic fluid in wiggle flow employing non-contact measurement techniques. The viscoelastic fluid was a 5% by weight solution of polyisobutylene dissolved in Primol 355, a high purity mineral oil.The kinematics were determined by Laser-Doppler Anemometry and compared with the data obtained for a Newtonian fluid, Primol 355, under identical flow conditions. It was found that the normalized axial velocity versus axial position curves along the centerline for both fluids superimposed at very low flow rates, an experimental verification that a viscoelastic fluid behaves like a Newtonian fluid under very low shear rates. However, at higher flow rates the behaviour of the viscoelastic fluid curves changed appreciably whereas the Newtonian fluid curves did not change at all. Thus, the effect of flow rate on viscoelastic fluid behaviour was also experimentally established.The normal stress differences were determined using a stress-birefringence apparatus. Data obtained along the centerline clearly exhibited a delayed growth of stress which should be attributed to the expected memory effects in viscoelastic fluid flow.     

Computation of rheological properties of suspensions of rigid rods: stress growth after inception of steady shear flow

The orientation distribution and stress growth for a suspension of rigid rods (or dumbbells) in a Newtonian solvent are calculated for inception of steady shear flow. Galerkin's method, with spherical harmonics as trial functions, is used in the spatial coordinates to obtain a system of ordinary differential equations in time which is solved by the spectral method. The method is applicable over a wide range of dimensionless shear rates (Peclet numbers) and has been coded with standard system-solvent and eigensystem packages. For sufficiently large Peclet numbers, the results give the well known rigid-dumbbell prediction of an overshoot in the shear viscosity and normal stress differences. This overshoot is then followed by an undershoot. An explicit analytical approximation for the fluid stresses is presented which is reasonably accurate for Peclet numbers less than unity.     

Plane surface suddenly set in motion in a viscoelastic fluid with fractional Maxwell model

The fractional calculus approach in the constitutive relationship model of viscoelastic fluid is introduced. The flow near a wall suddenly set in motion is studied for a non-Newtonian viscoelastic fluid with the fractional Maxwell model. Exact solutions of velocity and stress are obtained by using the discrete inverse Laplace transform of the sequential fractional derivatives. It is found that the effect of the fractional orders in the constitutive relationship on the flow field is significant. The results show that for small times there are appreciable viscoelastic effects on the shear stress at the plate, for large times the viscoelastic effects become weak. The project supported by the National Natural Science Foundation of China (10002003), Foundation for University Key Teacher by the Ministry of Education, Research Fund for the Doctoral Program of Higher Education     

Compliance effects on the torsional flow of a viscoelastic fluid

The effects of transducer compliance on transient stress measurements in torsional flows of a viscoelastic fluid are investigated theoretically. The analysis is based on the torsional flow of an upper-convected Maxwell fluid between a rotating and ‘stationary’ disk, which is allowed to twist and displace axially as a result of the stresses exerted on the disk by the fluid. An approximate analytical solution to the governing equations is obtained using a standard perturbation method. Results of the analysis are used to examine how the fluid velocity is altered by the motion of the stationary disk and to gain insight on how transient stress measurements are affected by transducer compliance. The analysis shows that compliance effects increase with applied shear rate and that the effects of torsional and axial compliance are coupled in measurements of the shear stress and first normal stress difference.     

Rheological behaviour of spray-dried egg yolk/xanthan gum aqueous dispersions

This paper provides information on the microstructure of, and reports particle size distributions and rheological results for, aqueous dispersions of spray-dried egg yolk and xanthan gum prepared on a laboratory scale using two types of homogenisers. Laser light scattering results demonstrated that higher energy input during homogenisation yielded a dispersion with a lower average particle size and a wider polydispersity, slightly influencing the linear dynamic viscoelastic functions due to the low concentration of egg yolk particles. These dispersions exhibited weak gel properties at the composition studied. The mechanical spectrum and the corresponding relaxation spectrum were dominated by the xanthan/gum-water matrix which controls the structure of the continuous phase. This fact explained the lack of any wall depletion effects. Several controlled-rate and controlled-stress rotational rheometers and a capillary rheometer were used to obtain information on flow properties. The shear rate dependence of steady state viscosity was determined through twelve decades, and was fitted using the Carreau equation. The kinetics of structural recovery after steady-state shear was studied by start-up at the inception of shear and flow interrupted experiments under controlled shear history. The results were analysed in terms of the ratio of a time-dependent amount of overshoot to the amount of overshoot of the original sample, using the addition of two first order equations. Additionally, combined steady state flow properties at fixed shear stress/oscillatory shear experiments were also used. The increase of the storage modulus with time, checking a linear viscoelastic response, tracked the structural recovery after steady shear. Laser light scattering of sheared samples helped gain a better understanding of the role of egg yolk particles on the rheology of these dispersions. Received: 6 February 2000 Accepted: 5 September 2000     

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.     

Observability of viscoelastic fluids

We apply the observability rank condition to study the observability of various viscoelastic fluids under imposed shear or extensional flows. In this paper the observability means the ability of determining the viscoelastic stress from the time history of the observations of the first normal stress difference. We consider four viscoelastic models: the upper convected Maxwell (UCM) model, the Phan–Thien–Tanner (PTT) model, the Johnson–Segalman (JS) model and the Giesekus model. Our study reveals that all of the four models have observability for all stress components almost everywhere under shear flow whereas under extensional flow most of the models have no observability for the shear stress component. More specifically, for UCM and JS models under imposed shear flow, the observations of the first normal stress difference allow the reconstruction of all components of viscoelastic stress. For UCM and JS models under extensional flow, the two normal stress components can be determined from the measurements of the first normal stress difference; the shear stress component does not affect the evolution of the normal stress components and consequently it cannot be extracted from the observations. Under shear flow, the PTT and Giesekus models have observability almost everywhere. That is, all components of the viscoelastic stress can be determined from the observations when the vector formed by the components of viscoelastic stress does not lie on a certain surface. Under extensional flow, the PTT model has observability almost everywhere for normal stress components whereas the Giesekus model has observability almost everywhere for all stress components. We also run simulations using the unscented Kalman filter (UKF) to reconstruct the viscoelastic stress from observations without and with noises. The UKF yields accurate and robust estimates for the viscoelastic stress both in the absence and in the presence of observation noises.     

Characteristic turbulent structure of a modified drag-reduced surfactant solution flow via dosing water from channel wall

The modification of the near-wall structure is very important for the control of wall turbulence. To ascertain the effect of near-wall modulation on the viscoelastic drag-reduced flow, the modified characteristics of a surfactant solution channel flow were investigated experimentally. The modulation was conducted on the boundary of the channel flow by injecting water from the whole surface of one side of the channel wall. The diffusion process of the injected water was observed by using the planar laser-induced fluorescence technique. The velocity statistics and characteristic structure including the spatial distributions of instantaneous streamwise velocity, swirling strength, and Reynolds shear stress were analyzed based on the velocity vectors acquired in the streamwise wall-normal plane by using the particle imaging velocimetry technique. The results indicated that the disturbance of the injected water was constricted within a finite range very near the dosing wall, and the Reynolds shear stress was increased in this region. However, the eventual drag reduction rate was found to be increased due to a relatively large decrement of viscoelastic shear stress in this near-wall region. Moreover, the flow structure under this modulation presented obvious regional characteristics. In the unstable disturbed region, the mixing of high-speed and low-speed fluids and the motions of ejection and sweep occurred actively. Many clockwise vortex cores were also found to be generated. This characteristic structure was similar to that in the ordinary turbulence of Newtonian fluid. Nevertheless, outside this disturbed region, the structure still maintained the characteristics of the drag-reduced flow with non-Newtonian viscoelastic additives. These results proved that the injected Newtonian fluid associated with the modified stress distribution creates a diverse characteristic structure and subsequent enhanced drag reduction. This investigation can provide the experimental basis for further study of turbulence control.     

Theoretical and experimental investigation of the flow of a viscoelastic fluid in the gap between two rotating disks

A theory of the nonlinear viscoelastic behavior of polymer fluids has been constructed in [1]. The theory was used in [2] to investigate the motion of a nonlinear viscoelastic medium under steady and unsteady deformation rates in simple shear flow, and a comparison was made with experiment. The experiments in [2], which were performed on a cone-plate Weissenberg rheogoniometer, indicate that this arrangement is unsuitable for measurements of normal stresses under unsteady conditions in fluids with a fairly high viscosity. Below, we will show the suitability of using a disk-disk Weissenberg rheogoniometer to measure normal stresses in this case for unsteady conditions (transition to steady flow and stress relaxation). In this regard, a theoretical study of the flow of a viscoelastic fluid in the gap between rotating disks is needed. Note that in this case new information will be obtained from a comparison with simple uniform shear flow, since in the flow of a polymer between two disks all three normal stress components contribute to the axial force, while in the gap between a cone and a plate only the first normal stress difference contributes to the normal force.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 25–30, March–April, 1976.     

Nonlinear viscoelasticity of PP/PS/SEBS blends

The nonlinear viscoelastic behavior of polypropylene/polystyrene (PP/PS) blends compatibilized or not with the linear triblock copolymer (styrene-ethylene-/butylene-styrene, SEBS) was investigated. Start-up of steady-shear at rates from 0.1 to 10 s^–1 was carried out using a controlled strain rotational rheometer and a sliding plate rheometer for strain histories involving one or several shear rates. The shear stress and first normal shear stress difference were measured as functions of time, and the morphologies of the samples before and after shearing were determined. For each strain history except that involving a single shear rate of 0.1 s^–1 the blends showed typical non-linear viscoelastic behavior: a shear stress overshoot/undershoot, depending on the history, followed by a steady state for each step. The first normal stress difference increased monotonically to a steady-state value. The values of the stresses increased with the addition of SEBS. The shear stress overshoot and undershoot and the times at which they occurred depended strongly on the strain history, decreasing for a subsequent shear rate step performed in the same direction as the former, and the time at which stress undershoot occurred increased for a subsequent shear rate step performed in the opposite direction, irrespective of the magnitude of the shear rate. This behavior was observed for all the blends studied. The time of overshoot in a single-step shear rate experiment is inversely proportional to the shear rate, and the steady-state value of N₁ scaled linearly with shear rate, whereas the steady-state shear stress did not. The average diameter of the dispersed phase decreased for all strain histories when the blend was not compatibilized. When the blend was compatibilized, the average diameter of the dispersed phase changed only during the stronger flows. Experimental data were compared with the predictions of a model formulated using ideas of Doi and Ohta (1991), Lacroix et al. (1998) and Bousmina et al. (2001). The model correctly predicted the behavior of the uncompatibilized blends for single-step shear rates but not that of the compatibilized blends, nor did it predict morphologies after shearing.     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method

An efficient immersed boundary-lattice Boltzmann method (IB-LBM) is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows. Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor. A boundary thickening-based direct forcing IB method is adopted to solve the particle–fluid interactions with high accuracy for non-slip boundary conditions. A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computational parameters. An asynchronous calculation strategy is adopted to further improve the computing efficiency. The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods. The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed. Subsequently, 504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh–Taylor instabilities was observed. At last, simulations of a dense flow involving 11001 particles, the largest number of particles to date, were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid. The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time. These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.     

Negative wake in the uniform flow past a cylinder

The upstream/downstream streamline shift and the associated negative wake generation (streamwise velocity overshoot in the wake) in a viscoelastic flow past a cylinder are studied in this paper, for the Oldroyd-B, UCM, PTT, and FENE-CR fluids, using the Discrete Elastic Viscous Split Stress Vorticity (DEVSS-ω) scheme (Dou HS, Phan-Thien N (1999). The flow of an Oldroyd-B fluid past a cylinder in a channel: adaptive viscosity vorticity (DAVSS-ω) formulation. J Non-Newtonian Fluid Mech 87:47–73). The numerical algorithm is a parallelized unstructured Finite Volume Method (FVM), running under a distributed computing environment through the Parallel Virtual Machine (PVM) library. It is demonstrated that both the normal stress and its gradient are responsible for the negative wake generation and streamline shifting. Fluid extensional rheology plays an important role in the generation of the negative wake. The negative wake can occur in flows where the fluid extensional viscosity does not increase rapidly with strain rate. The formation of the negative wake does not depend on whether the streamlines undergo an upstream or a downstream shift. Shear-thinning viscosity weakens the velocity overshoot and while shear-thinning first normal stress coefficient enhances the velocity overshoot. Wall proximity is not necessary for the velocity overshoot; however, it enhances the strength of the negative wake. For the Oldroyd-B fluid, the ratio of the solvent viscosity to the zero-shear viscosity plays an important role in the streamline shift. In addition, mesh dependent behaviour of normal stresses along the centreline at high De in most cylinder/sphere simulations is due to the convection of normal stress from the cylinder to the wake, which results in the maximum of the normal stress being located off the centreline by a short distance at high De.