Rheo-PIV of a yield-stress fluid in a capillary with slip at the wall

An analysis of the yielding and flow behavior of a model yield-stress fluid, 0.2 wt% Carbopol gel, in a capillary with slip at the wall has been carried out in the present work. For this, a study of the flow kinematics in a capillary rheometer was performed with a two-dimensional particle image velocimetry (PIV) system. Besides, a stress-controlled rotational rheometer with a vane rotor was used as an independent way to measure the yield stress. The results in this work show that in the limit of resolution of the PIV technique, the flow behavior agrees with the existence of a yield stress, but there is a smooth solid?Cliquid transition in the capillary flow curve, which complicates the determination of the yield stress from rheometrical data. This complication, however, is overcome by using the solely velocity profiles and the measured wall shear stresses, from which the yield-stress value is reliably determined. The main details of the kinematics in the presence of slip were all captured during the experiments, namely, a purely plug flow before yielding, the solid?Cliquid transition, as well as the behavior under flow, respectively. Finally, it was found that the slip velocity increases in a power-law way with the shear stress.     

Stick-slip flow of high density polyethylene in a transparent slit die investigated by laser Doppler velocimetry

The oscillating flow instability of a molten linear high-density polyethylene is carefully studied using a single screw extruder equipped with a transparent slit die. Experiments are performed using laser Doppler velocimetry in order to obtain the local velocities field across the entire die width. At low flow rate, the extrusion is stable and steady state velocity profiles are obtained. During the instability, the velocity oscillates between two steady state limits, suggesting a periodic stick-slip transition mechanism. At high flow rate, the flow is mainly characterized by a pronounced wall slip. We show that wall slip occurs all along the die land. An investigation of the slip flow conditions shows that wall slip is not homogeneous in a cross section of the slit die, and that pure plug flow occurs only for very high flow rates. A numerical computation of the profile assuming wall slip boundary conditions is done to obtain the true local wall slip velocity. It confirms that slip velocities are of the same order of magnitude as those measured with a capillary rheometer.     

Enhanced melt strength and stretching of linear low-density polyethylene extruded under strong slip conditions

The influence of extrusion under strong slip conditions on the extensional properties of linear low-density polyethylene was studied in this work. The material was extruded at two different temperatures under strong slip and no slip conditions, and was subsequently subjected to uniaxial elongational flow by means of a Rheotens device. Strong slip was evident through the elimination of sharkskin distortions and the stick-slip instability, as well as by the electrification of the extrudates. The extrudate swell was smaller in the presence of slip when comparing with no slip conditions at constant apparent shear rate, but it was found to be a unique function of the shear stress if comparison was performed at constant stress. The draw ratio and melt strength of the filaments obtained under slip conditions were larger compared to those without slip. In addition, draw resonance was postponed to higher draw ratios during the extrusion with strong slip at constant apparent shear rate. It is suggested that slip of the polymer at the die wall decreases the shear stress in the bulk, and therefore, restricts the disentanglement and orientation of macromolecules during flow, which subsequently produces the increase in draw ratio and melt strength during stretching.     

Particle imaging velocimetry of a wet aqueous foam with an underlying liquid film

 We investigate the utility of particle imaging velocimetry (PIV) for performing kinematic measurements in wet aqueous foam with a liquid film beneath it. The flow velocities are measured near the walls of a square cross-section horizontal duct. The flow velocities are useful for validating the rheological models. We show that there is a discrepancy between the velocity profiles in the wet foam and the Bingham plastic model of flow. The velocity measurements reveal a more complex flow pattern, which may be analysed following three different regimes: a plug flow, a shear flow in a vertical plane and a three-dimensional shear flow. The transition between the plug flow and the shear flows may be explained by a shear-induced migration of bubbles. Received: 25 April 2000 / Accepted: 26 February 2001     

Slip at the interface between immiscible polymer melts I: method to measure slip

Slip at the interface between immiscible polymer melts remains poorly understood. A method that relies solely on rheological measurements to obtain the interfacial slip velocity uses the slip-induced deviation in the flow variables. To use the method, accurate estimates of the flow variables under the assumption of no-slip are necessary. Although such estimates can be easily derived under some cases, in general, this is not straightforward. Therefore, methods to determine the interfacial slip velocity without using estimates for the flow variables under no-slip conditions are desirable. In this work, we focus on investigations of slip at the interface between two immiscible polymer melts undergoing two-phase coaxial flow. To enable such investigations, we have adapted the Mooney method, usually used to investigate wall slip, to investigate polymer/polymer interfacial slip. Using this method, we have measured the slip velocity at the interface between polypropylene and polystyrene as a function of the interfacial stress. To determine the validity of the modified Mooney method, we also determine the slip velocity using the slip-induced deviation in the flow variables. To enable this determination, we use polypropylene and polystyrene with almost identical shear rate-dependent viscosities over a range of shear rates. The slip velocity obtained from the modified Mooney method displayed excellent agreement with that determined using the deviation from no-slip. In agreement with prior work, the dependence of the slip velocity on the interfacial stress is a power-law. Our investigation spans a sufficiently wide range of interfacial stress to enable the direct observation of two power-law regimes and also the transition between the two regimes. We also find that the power-law exponent of approximately 3 at low stresses decreases to approximately 2 at high stresses.     

Extrusion of a compressible Newtonian fluid with periodic inflow and slip at the wall

We explore a mechanism of extrusion instability, based on the combination of nonlinear slip and compressibility. We consider the time-dependent compressible Newtonian extrudate swell problem with slip at the wall. Steady-state solutions are unstable in regimes where the shear stress is a decreasing function of the velocity at the wall. Compressibility provides the means for the alternate storage and release of elastic energy, and, consequently, gives rise to periodic solutions. The added novelty in the present work is the assumption of periodic volumetric flow rate at the inlet of the die. This leads to more involved periodic responses and to free surface oscillations similar to those observed experimentally with the stick-slip instability. To numerically simulate the flow, we use finite elements in space and a fully-implicit scheme in time.Dedicated to the memory of Prof. Tasos Papanastasiou     

A linear model of the stick-slip phenomena in polymer flow in rheometers

A linear approach was employed for the qualitative theoretical study of the stick-slip phenomena in polymer flows in rheometers. For this purpose, the familiar three-constant rheological equation and a linear wall friction law (i.e. wall stress — slip velocity dependence) were used with some additional hypotheses about the onset of the stick-slip behaviour. The friction law used was derived from a crude molecular approach. On the basis of these equations the inertialess stickslip behaviour of a viscoelastic liquid flowing through a capillary at a constant flow rate was considered. To be able to describe some transient phenomena in this problem, inertial effects (as an example) were taken into account. Furthermore, the distortions on the extrudate surface due to the slip phenomena inside the capillary were described theoretically within the framework of a linear approximation. In the final part of the study the possibility of rapid stochastization was discussed for rotational stick-slip flow of polymers in a cone-plate rheometer.     

Rheology of shear thickening suspensions and the effects of wall slip in torsional flow

The rheological characterisation of concentrated shear thickening materials suspensions is challenging, as complicated and occasionally discontinuous rheograms are produced. Wall slip is often apparent and when combined with a shear thickening fluid the usual means of calculating rim shear stress in torsional flow is inaccurate due to a more complex flow field. As the flow is no longer “controlled”, a rheological model must be assumed and the wall boundary conditions are redefined to allow for slip. A technique is described where, by examining the angular velocity response in very low torque experiments, it is possible to indirectly measure the wall slip velocity. The suspension is then tested at higher applied torques and different rheometer gaps. The results are integrated numerically to produce shear stress and shear rate values. This enables the measurement of true suspension bulk flow properties and wall slip velocity, with simple rheological models describing the observed complex rheograms.     

In vitro post-stenotic flow quantification and validation using echo particle image velocimetry (Echo PIV)

Echo particle image velocimetry (Echo PIV) presents itself as an attractive in vivo flow quantification technique to traditional approaches. Promising results have been acquired; however, limited quantification and validation is available for post-stenotic flows. We focus here on the comprehensive evaluation of in vitro downstream stenotic flow quantified by Echo PIV and validated in relation to digital particle image velocimetry (DPIV). A Newtonian blood analog was circulated through a closed flow loop and quantified immediately downstream of a 50 % axisymmetric blockage at two Reynolds numbers (Re) using time-averaged Echo PIV and DPIV. Centerline velocities were in good agreement at all Re; however, Echo PIV measurements presented with elevated standard deviation (SD) at all measurements points. SD was improved using increased line density (LD); however, frame rate or field of view (FOV) is compromised. Radial velocity profiles showed close agreement with DPIV with the largest disparity in the shear layer and near-wall recirculation. Downstream recirculation zones were resolved by Echo PIV at both Re; however, magnitude and spatial coverage was reduced compared to DPIV that coincided with reduced contrast agent penetration beyond the shear layer. Our findings support the use of increased LD at a cost to FOV and highlight reduced microbubble penetration beyond the shear layer. High local SD at near-wall measurements suggests that further refinement is required before proceeding to in vivo quantification studies of wall shear stress in complex flow environments.     

Effect of pressure-dependent slip on flow curve multiplicity

Various microstructural pictures for slip at polymer/solid interfaces lead to relations which have a region where multiple values of slip velocity are predicted for the same shear stress. This leads to the expectation of multivalued flow curves, which has been verified in specific cases by numerous researchers. We study the effect of pressure dependence on flow curve multiplicity using a simple multivalued slip relation to model the phenomena of hysteresis and spurt flow in polymer extrusion. A continuation technique is used to trace out the boundaries of the region of flow curve multiplicity as pressure drop and die length to diameter (L/D) ratio are changed. Results for Newtonian, shear thinning and viscoelastic constitutive equations show that, despite the multivalued nature of the slip model, multiplicity (and thus hysteresis) is absent at high L/D.  For the sake of completeness, we also carry out time-dependent simulations at constant piston speed taking fluid compressibility into account. These simulations show that oscillations in the pressure drop and exit volumetric flow rate result only if the system is operated in the multiplicity region of the steady state flow curve, in agreement with the results of similar simulations by researchers using various multivalued slip models without pressure dependence. The results demonstrate that a multivalued slip model does not guarantee multiplicity in the flow curve for the constant pressure drop operation, nor oscillations for constant piston speed operation. Received: 18 August 1997 Accepted: 30 March 1998     

PIV measurements of flow in drying polymer solutions during solvent casting

An experimental method based on confocal microscopy and particle image velocimetry (PIV) is used to characterize the flow in a polymer solution during solvent casting. The flow inside a 200-μm-thick film of a poly(vinyl alcohol) (PVA) solution is visualized near a vertical wall of a mold using confocal microscopy of seed particles during solvent evaporation at 25, 35, and 45°C, and the corresponding velocity vector fields are determined from projections of the confocal images. Flow toward the vertical wall is observed inside the film as well as a slower Marangoni-type counter flow at the film surface during the initial phase of solvent evaporation, resulting from a polymer concentration gradient along the film due to a local variation in evaporation rate. Total volume of the polymer solution in the observation volume as well as solvent evaporation rate are determined as a function of time, both revealing close correlation to average horizontal velocity data from PIV. The PIV measurements show significant differences in the flow velocity fields at different temperatures. The PIV measurements correlate with the solvent evaporation rates as well as the final polymer thicknesses on the vertical wall of the mold. Surface tension and viscosity measurements are taken for different concentrations of PVA solution.     

Determination of the wall slip velocity in the flow of a SBR compound

Rubber compounds are known to exhibit slip at the wall in particular flow conditions. The slip velocity is usually determined by using the classical Mooney method. The rheological behavior of a styrene butadiene rubber (SBR) compound was studied with three different rheometers. Biconical rotational, capillary and slit die rheometers were used to define the true viscous behavior of the compound and the slip velocity. It was shown that it was impossible to apply the Mooney method to our experimental data. New characterizations were thus developed for both capillary and slit die experiments. They were based on the dependency of the slip velocity on the local flow gap. Contrarily to the Mooney method, they provided physically acceptable results and led to a power-law relationship between wall slip, wall shear stress and local geometry of the flow.     

Rheology of pulp suspensions using ultrasonic Doppler velocimetry

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

Stick-slip transition in capillary flow of linear polyethylene: 3. Surface conditions

Effects of surface topology and energy on the stick-slip transition were studied in capillary flow of highly entangled polyethylene (PE) melts. Surface roughness was shown to increase the critical stress of the stick-slip transition because of the increased resistance to interfacial disentanglement. Lowering the surface energy of a smooth die wall by treatment with a fluorocarbon completely eliminates the stick-slip transition and produces massive interfacial slip at PE/wall boundary down to a stress level of 0.05 MPa. On the other hand, considerable roughness on the same low energy surface can produce a stick hydrodynamic boundary condition and restore the stick-slip transition despite the weak PE/wall interfacial interactions. Additionally, a slip-slip transition was found in the die with a nearly non-adsorbing wall that appears to involve a secondary chain-debonding process.     

Specific rheology – morphology relationships for some blends containing LCPs

For the blend melts of isotropic polysulfone (PSF) and LC polyester (PES), differing in viscosity more than 10 times, the flow curves with maxima were observed in cone and plate geometry. The low shear rate branch is located near the PSF flow curve, and the high shear rate branch is close to the PES flow curve. At high strains, the formation of the ring-like morphology of the blend sample, accompanied by appearance of maximum on flow curve, was registered. The scaling analysis of the reasons for the ring morphology formation was based on stretching of the large, low-viscous LC droplet, embedded to the high-viscous polymer matrix, in a homogeneous shear field. It was shown that, if the critical Taylor radius is not exceeded, the droplet may form the closed torus. Under strong flows, the PSF melt manifests the “spurt effect”, consisting of a drastic increase of the shear rate when the critical value of the shear stress is reached. The pattern of the blend flow curves with maxima may be explained by a vanishing PSF input to the total shear stress, inherent for blends, while the PES melt continues to be in a liquid state and, consequently, is responsible for the blend viscosity at the high shear rates. The presence of regular heterogeneities in the blend in the form of LC rings may initiate the rupture of the entanglements network of the matrix PSF (close to LC rings) under strong shear flows. The appearance of the low-viscous “cracks” at the critical shear stress will diminish the contribution of the PSF to the blends rheological response. Received: 20 April 1999 Accepted: 28 January 2000     

Turbulent free convection in a vertical converging channel

This paper presents the results of experimental and numerical investigations of the problem of turbulent natural convection in a converging-plate vertical channel. The channel has two isothermally heated inclined walls and two adiabatic vertical side walls. The parameters involved in this study are the channel geometry represented by the channel width at exit, the inclination of the heated walls and the temperature difference between the heated walls and the ambient. The investigation covered modified Rayleigh numbers up to 10⁸ in the computational study and up to 9.3 × 10⁶ in the experimental work. The experimental measurements focused on the velocity field and were carried out using a PIV system and included measurements of the mean velocity profiles as well as the root-mean-square velocity and shear stress profiles. The experiments were conducted for an inclination angle of 30°, a gap width of 10 mm and two temperature differences (∆T=25.4°C and 49.8°C). The velocity profiles in the lower part of the channel indicated the presence of two distinct layers. The first layer is adjacent to the heated plate and driven by buoyancy forces while the second layer extends from the point of maximum velocity to the channel center plane and driven mainly by shear forces. The velocity profile at the upper portion of the channel has shown the merging of the two boundary layers growing over the two heated walls. The measured values of the Reynolds shear stress and root mean square of the horizontal and vertical velocity fluctuation components have reached their maximum near the wall while having smaller values in the core region. The computational results have shown that the average Nusselt number increases approximately linearly with the increase of the modified Rayleigh number when plotted on log–log scale. The variation of the local Nusselt number indicated infinite values at the channel inlet (leading edge effect) and high values at the channel exit (trailing edge effect). For a fixed value of the top channel opening, the increase of the inclination angle tended to reduce flow velocity at the inlet section while changing the flow structure near the heated plates in such a way to create boundary-layer type flow. The maximum value of the average Nusselt number occurs when θ = 0 and decreases with the increase of the inclination angle. On the other hand, the increase of the channel width at exit for the same inclination angle caused a monotonic increase in the flow velocity at the channel inlet.     

Optimization of multiplane μPIV for wall shear stress and wall topography characterization

Multiplane μPIV can be utilized to determine the wall shear stress and wall topology from the measured flow over a structured surface. A theoretical model was developed to predict the measurement error for the surface topography and shear stress, based on a theoretical analysis of the precision in PIV measurements. The main parameters that affect the accuracy of the measurement are identified. The effect of different parameter settings is studied by means of Monte Carlo simulations, and the results are compared with an experimental test case. The results are used to determine the recommended parameter settings for this measurement approach.     

Simultaneous determination of shear and extensional viscosities and wall slip by on-line rheometry with parameter fitting: an application to EPDM rubber

Shear and extensional viscosities and wall slip are determined simultaneously under extrusion processing conditions using an on-line rheometer. Because it is not possible to independently control flow rate and temperature, classical methods for interpretation of capillary data cannot be used with on-line rheometry. This limitation is overcome using computational optimization to fit parameters in a flow model. This consists of three parts, representing shear viscosity, extensional viscosity, and wall slip. Three-parameter, power law forms, based on local instantaneous deformation rates and including temperature dependence, are used for each, and analytic solutions applied for entry flow and flow in the capillary. For entry flow, the Cogswell–Binding approach is used, and for developed flow in the capillary a solution incorporating wall slip is derived. The rheometer, with interchangeable capillaries, is mounted in place of the die on a rubber profile extrusion line. Pressure drops and temperatures for extrusion of an EPDM rubber through 2 mm diameter capillaries of length 0, 2, 3, 4, and 5 mm are logged and flow rates determined for a range of extruder speeds (5 to 20 rpm). Pressures ranged from 60 to 75 bar and temperatures from 86 to 116 °C. Mean flow velocity in the capillaries was between 5 × 10⁻³ and 5 × 10⁻¹ m s⁻¹. The nine material parameters are optimized for best fit of the analytic pressure drops to experimental data, using about 100 data points, with the Levenberg–Marquardt method. It is concluded that flow is dominated by extension and wall slip. Shear flow appears to play little part. The slip model indicates that slip velocity increases much more rapidly than the wall shear stress (in the range 0.5–1 MPa) and decreases with temperature for a given stress level. Results for the (uniaxial) extensional viscosity represent an engineering approximation to this complex phenomenon at the high strains (approximately 200) and high extension rates (up to 800 s⁻¹) applying in the extrusion. Results indicate a slight extension hardening and a decrease with temperature. Results are put into the context of the available studies in the literature, which, particularly with regard to wall-slip and extensional flow, consider conditions far removed from those applying in industrial extrusion. The present methods provide a powerful means for flow characterization under processing conditions, providing data suitable for use in computer simulations of extrusion and optimization of die design.     

Apparent wall slip velocity measurements in free surface flow of concentrated suspensions

In this work we have experimentally measured the apparent wall slip velocity in open channel flow of neutrally buoyant suspension of non-colloidal particles. The free surface velocity profile was measured using the tool of particle imaging velocimetry (PIV) for two different channels made of plane and rough walls. The rough walled channel prevents wall slip, whereas the plane wall showed significant wall slip due to formation of slip layer. By comparing the velocity profiles from these two cases we were able to determine the apparent wall slip velocity. This method allows characterization of wall slip in suspension of large sized particles which cannot be performed in conventional rheometers. Experiments were carried out for concentrated suspensions of various particle volume concentrations and for two different sizes of particles. It was observed that wall slip velocity increases with particle size and concentration but decreases with increase in the viscosity of suspending fluid. The apparent wall slip velocity coefficients are in qualitative agreement with the earlier measurements. The effect of wall slip on free surface corrugation was also studied by analyzing the power spectral density (PSD) of the refracted light from the free surface. Our results indicate that free surface corrugation is a bulk flow response and it does not arise from boundary problem such as development of slip layer.     

Non-oscillating/Oscillating Shear Layer over a Large Deep Cavity at Low-Subsonic Speeds

The flow over a deep cavity at low subsonic velocity is considered in the present paper. The cavity length-to-depth aspect ratio is L/H = 0.2. Single hot-wire measurements characterized the incident turbulent boundary layer and show the influence of the cavity on the streamwise statistic components just downstream from the cavity. The streamwise mean and fluctuating velocity profiles are affected by the cavity. PIV measurements reveal the presence for ejection-like events responsible of local perturbations of the skewness and the flatness coefficients. Time-resolved PIV technic is also used to characterize phase properties of shear layer oscillating cycle. It is shown that for deep cavity with first Rossiter mode, only one vortical structure is formed at the cavity leading edge. Then, it grows while convecting downstream along the shear layer. A well-defined ejection process begins after the vortex impact near the cavity downstream corner. A cylinder device placed spanwisely near the cavity leading edge eliminates the resonance and highly modifies the behavior of the shear layer flow. In fact, the shear layer could be divided into upper and lower parts with different structure aspects.