Comparison of a quasi-Newtonian fluid with a viscoelastic fluid in planar contraction flow

The flow of a 5.0 wt.% solution of polyisobutylene in tetradecane through a planar 4 : 1 contraction exhibiting a shear thinning viscosity is simulated using the flow-type sensitive quasi-Newtonian fluid model. The shear viscosity is fitted by the Giesekus model, which, with the chosen parameters, leads to an extension thickening elongational viscosity. The stress and velocity fields of the numerical simulations are compared with the experimental results of Quinzani et al. [J. Non-Newtonian Fluid Mech. 52 (1994) 1–36] and the numerical results of the viscoelastic simulation using the Giesekus model of Azaiez et al. [J. Non-Newtonian Fluid Mech. 62 (1996) 253–277]. It can be shown that the quasi-Newtonian fluid qualitatively predicts the essential features of the flow in the vicinity of the contraction.     

A method for calculating rheological and morphological properties of constant-volume polymer blend models in inhomogeneous shear fields

We show how to formulate two-point boundary-value problems in order to compute fully-developed laminar channel and tube flow profiles for viscoelastic fluid models. The formulation is applied to Couette and pressure-driven flows separately, or a combination of both. The application of this methodology is illustrated analytically for the Upper-Convected Maxwell Model, and it is applied computationally for the Phan-Thien/Tanner and Giesekus Models. Numerical solutions exist for the last two models [J.Y. Yoo, H.C. Choi, On the steady simple shear flows of the one-mode Giesekus fluid, Rheol. Acta 28 (1989) 13–24; P.J. Oliveira, F.T. Pinho, Analytical solution for fully developed channel and pipe flow of Phan-Thien–Tanner fluids, J. Fluid Mech. 387 (1999) 271–280; M.A. Alves, F.T. Pinho, P.J. Oliveira, Study of steady pipe and channel flows of a single-mode Phan-Thien–Tanner fluid, J. Non-Newtonian Fluid Mech. 101 (2001) 55–76], allowing verification of the computational technique. Subsequently, the computational algorithm is applied to the constant-volume polymer blend models of Maffettone and Minale [P.L. Maffettone, M. Minale, Equation of change for ellipsoidal drops in viscous flow, J. Non-Newtonian Fluid Mech. 84 (1999) 105–106 (Erratum), J. Non-Newtonian Fluid Mech. 78 (1998) 227–241] and Dressler and Edwards [M. Dressler, B.J. Edwards, The influence of matrix viscoelasticity on the rheology of polymer blends, Rheol. Acta 43 (2004) 257–282; M. Dressler, B.J. Edwards, Rheology of polymer blends with matrix-phase viscoelasticity and a narrow droplet size distribution, J. Non-Newtonian Fluid Mech. 120 (2004) 189–205]. Rheological and morphological properties of the model blends are thus obtained as functions of the spatial position within the channel, applied pressure drop, and shear rate at the wall.     

Dynamics of viscoelastic liquid filaments: Low capillary number flows

Many applications of viscoelastic free surface flows requiring formation of drops from small nozzles, e.g., ink-jet printing, micro-arraying, and atomization, involve predominantly extensional deformations of liquid filaments. The capillary number, which represents the ratio of viscous to surface tension forces, is small in such processes when drops of water-like liquids are formed. The dynamics of extensional deformations of viscoelastic liquids that are weakly strain hardening, i.e., liquids for which the growth in the extensional viscosity is small and bounded, are here modeled by the Giesekus, FENE-P, and FENE-CR constitutive relations and studied at low capillary numbers using full 2D numerical computations. A new computational algorithm using the general conformation tensor based constitutive equation [M. Pasquali, L.E. Scriven, Theoretical modeling of microstructured liquids: a simple thermodynamic approach, J. Non-Newtonian Fluid Mech. 120 (2004) 101–135] to compute the time dependent viscoelastic free surface flows is presented. DEVSS-TG/SUPG mixed finite element method [M. Pasquali, L.E. Scriven, Free surface flows of polymer solutions with models based on conformation tensor, J. Non-Newtonian Fluid Mech. 108 (2002) 363–409] is used for the spatial discretization and a fully implicit second-order predictor–corrector scheme is used for the time integration. Inertia, capillarity, and viscoelasticity are incorporated in the computations and the free surface shapes are computed along with all the other field variables in a fully coupled way. Among the three models, Giesekus filaments show the most drastic thinning in the low capillary number regime. The dependence of the transient Trouton ratio on the capillary number in the Giesekus model is demonstrated. The elastic unloading near the end plates is investigated using both kinematic [M. Yao, G.H. McKinley, B. Debbaut, Extensional deformation, stress relaxation and necking failure of viscoelastic filaments, J. Non-Newtonian Fluid Mech. 79 (1998) 469–501] and energy analyses. The magnitude of elastic unloading, which increases with growing elasticity, is shown to be the largest for Giesekus filaments, thereby suggesting that necking and elastic unloading are related.     

Transient 3D flow of polymer solutions: A Lagrangian computational method

We describe a computational method for the numerical simulation of three-dimensional transient flows of polymer solutions that extends the work of Harlen et al. [O.G. Harlen, J.M. Rallison, P. Szabó, A split Lagrangian–Eulerian method for simulating transient viscoelastic flows, J. Non-Newtonian Fluid Mech. 60 (1995) 81–104]. The method uses a Lagrangian computation of the stress together with an Eulerian computation of the velocity field. Adaptive mesh reconnection based on Delaunay tetrahedra is used to ensure well-shaped elements. Additional shape-quality improvement procedures are developed to improve the algorithm. We validate the method for the benchmark problem of a rigid sphere falling in a cylindrical pipe. Inertia is neglected. We compare results for the axisymmetric case with previous work (using a FENE model), and then consider the off-axis non-axisymmetric case. In the latter case, we find that as the sphere falls, it drifts across the pipe, a phenomenon previously observed in experiments but not fully explained. The physical mechanisms that cause the time-dependent drift are identified, and a simple model based on the normal stresses in the fluid is shown to predict the magnitude of the drift velocity.We also consider a second benchmark problem involving a constriction in an axisymmetric pipe. Numerical difficulties associated with ill-shaped elements near the concave boundary arise for higher Weissenberg numbers. The merits and drawbacks of the new numerical method, and its applicability to various flow geometries are discussed.     

An adaptive remeshing strategy for viscoelastic fluid flow simulations

In the last few years, we have developed a fairly general adaptive finite element solution procedure which can be applied to a large variety of problems. In this paper, this strategy is briefly recalled and applied to the solution of two-dimensional viscoelastic fluid flow problems. A log-conformation formulation recently introduced by Fattal and Kupferman [R. Fattal, R. Kupferman, Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation, J. Non-Newtonian Fluid Mech. 126 (2005) 23-37] was implemented in order to improve the convergence properties of the numerical scheme. We confirm some results obtained in Hulsen, Fattal and Kupferman [M. Hulsen, R. Fattal, R. Kupferman, Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithm, J. Non-Newtonian Fluid Mech. 127 (2005) 27-39] and in some instances, we show that mesh adaptation allows to almost automatically reproduce accurate results obtained on very fine structured meshes.     

Multiscale simulation of viscoelastic free surface flows

A micro–macro approach based on combining the Brownian configuration fields (BCF) method [M.A. Hulsen, A.P.G. van Heel, B.H.A.A. van den Brule, Simulation of viscoelastic flow using Brownian configuration fields, J. Non-Newtonian Fluid Mech. 70 (1997) 79–101] with an Arbitrary Lagrangian–Eulerian (ALE) Galerkin finite element method, using elliptic mesh generation equations coupled with time-dependent conservation equations, is applied to study slot coating flows of polymer solutions. The polymer molecules are represented by dumbbells with both linear and non-linear springs; hydrodynamic interactions between beads are incorporated. Calculations with infinitely extensible (Hookean) and pre-averaged finitely extensible (FENE-P) dumbbell models are performed and compared with equivalent closed-form macroscopic models in a conformation tensor based formulation [M. Pasquali, L.E. Scriven, Free surface flows of polymer solutions with models based on the conformation tensor, J. Non-Newtonian Fluid Mech. 108 (2002) 363–409]. The BCF equation for linear dumbbell models is solved using a fully implicit time integration scheme which is found to be more stable than the explicit Euler scheme used previously to compute complex flows. We find excellent agreement between the results of the BCF based formulation and the macroscopic conformation tensor based formulation. The computations using the BCF approach are stable at much higher Weissenberg numbers, (where λ is the characteristic relaxation time of polymer, and is the characteristic rate of strain) compared to the purely macroscopic conformation tensor based approach, which fail beyond a maximum Wi. A novel computational algorithm is introduced to compute complex flows with non-linear microscopic constitutive models (i.e. non-linear FENE dumbbells and dumbbells with hydrodynamic interactions) for which no closed-form constitutive equations exist. This algorithm is fast and computationally efficient when compared to both an explicit scheme and a fully implicit scheme involving the solution of the non-linear equations with Newton’s method for each configuration field.     

Direct simulation of flexible fibers

A numerical simulation of multiple flexible fibers in suspension in Newtonian simple shear flow is presented. The method used is similar to those of previous recent simulation works by Fan et al. [J. Non-Newtonian Fluid Mech. 74 (1998) 113] and Yamane et al. [J. Non-Newtonian Fluid Mech. 54 (1994) 405], however, the method has been modified to allow a small amount of bending and torsion in the fibers. A restoring moment acts to straighten the fibers as they interact in the flow.It is demonstrated that this simulation can be used to extract basic rheological information about the suspension including fiber orientations and suspension viscosity. The viscosity of semi-concentrated to concentrated flexible fiber suspensions are shown to increase by a magnitude of the order 7–10% greater than the equivalent rigid fiber suspension tested. This is in qualitative agreement with previous experimental work by Goto et al. [Rheologica Acta 25 (1986) 119] and Blakeney [J. Colloid Interface Sci. 22 (1966) 324]. The implication is that any constitutive relation involving particulate suspensions described by orientation vectors may quantitatively underestimate suspension viscosity, particularly for fibers of large aspect ratio, or low Young’s modulus, whereby the tendency to flex is greater [Rheologica Acta 25 (1986) 119]. If particulate deformation were accounted for (by whatever means) in the existing constitutive relationship, predictions of bulk suspension parameters such as viscosity should be noticeably improved. A method is developed to modify an existing rigid-fiber viscosity to an equivalent flexible fiber viscosity, hence improving viscosity prediction ability.     

A rising bubble in a polymer solution

We propose a boundary integral method to study the shape of a bubble rising under gravity in a dilute polymer solution. Constitutive properties are modelled using a FENE model [M.D. Chilcott, J.M. Rallison, J. Non-Newtonian Fluid Mech. 29 (1988) 381] with a pure surface tension interface. We employ a birefringent strand representation [O.G. Harlen, J.M. Rallison, M.D. Chilcott, High-Deborah-number flows of dilute polymer, J. Non-Newtonian Fluid Mech.34 (1990) 319–349] of the wake to simulate the shape and the time-dependent motion of the bubble. Steady and non-steady solutions reproduce qualitatively the bubble deformation seen in experiment with a small region of very high curvature near the rear stagnation point of the bubble. We find a limit point for steady axisymmetric solutions if the polymer concentration is increased or the surface tension is decreased. Rise speed jump discontinuities were not found.     

Nonequilibrium stretching dynamics of dilute and entangled linear polymers in extensional flow

We propose an extension of the FENE-CR model for dilute polymer solutions [M.D. Chilcott, J.M. Rallison, Creeping flow of dilute polymer solutions past cylinders and spheres, J. Non-Newtonian Fluid Mech. 29 (1988) 382–432] and the Rouse-CCR tube model for linear entangled polymers [A.E. Likhtman, R.S. Graham, Simple constitutive equation for linear polymer melts derived from molecular theory: Rolie–Poly equation, J. Non-Newtonian Fluid Mech. 114 (2003) 1–12], to describe the nonequilibrium stretching dynamics of polymer chains in strong extensional flows. The resulting models, designed to capture the progressive changes in the average internal structure (kinked state) of the polymer chain, include an ‘effective’ maximum contour length that depends on local flow dynamics. The rheological behavior of the modified models is compared with various results already published in the literature for entangled polystyrene solutions, and for the Kramers chain model (dilute polymer solutions). It is shown that the FENE-CR model with an ‘effective’ maximum contour length is able to describe correctly the hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation also observed and computed by Doyle et al. [P.S. Doyle, E.S.G. Shaqfeh, G.H. McKinley, S.H. Spiegelberg, Relaxation of dilute polymer solutions following extensional flow, J. Non-Newtonian Fluid Mech. 76 (1998) 79–110] and Li and Larson [L. Li, R.G. Larson, Excluded volume effects on the birefringence and stress of dilute polymer solutions in extensional flow, Rheol. Acta 39 (2000) 419–427] using Brownian dynamics simulations of bead–spring model. The Rolie–Poly model with an ‘effective’ maximum contour length exhibits a less pronounced hysteretic behavior in stress versus birefringence in start-up of uniaxial extensional flow and subsequent relaxation.     

A fictitious domain/distributed Lagrange multiplier method for the particulate flow of Oldroyd-B fluids: A positive definiteness preserving approach

In this article we present a numerical method for simulating the sedimentation of circular particles in a two-dimensional channel filled with an Oldroyd-B fluid. We have combined a fictitious domain/distributed Lagrange multiplier method with a factorization approach from Lozinski and Owens [J. Non-Newtonian Fluid Mech. 112 (2003) 161] via an operator splitting technique. The new scheme preserves the positive definiteness of the conformation tensor at the discrete level. The method is validated by performing a convergence study which shows that the results are independent of the mesh and time step sizes. Our results show that when the elasticity number (E) is less than a critical value (which depends upon the blockage ratio), two particles will sediment in the channel-like particles in Newtonian fluids; when the elasticity number is greater than the critical value, chains are formed for the case of two particles sedimenting in an Oldroyd-B fluid and the center line is aligned with the falling direction. These results agree with those presented in [P.Y. Huang, H.H. Hu, and D.D. Joseph, J. Fluid Mech. 362 (1998) 297]. For the cases of three and six particles, when the elasticity number is greater than a critical value and the viscoelastic Mach number is less than one, chains are also formed and move to the center of the channel.     

Effect of purely elastic flow instability on molecular conformation and drag

Linear stability analysis has shown that viscoelastic creeping flow of an Oldroyd-B liquid through a sinusoidal channel is unstable to stationary, wall-localized and short wavelength perturbations [B. Sadanandan, R. Sureshkumar, Global linear stability analysis of non-separated viscoelastic flow through a periodically constricted channel, J. Non-Newtonian Fluid Mech. 122 (2004) 55]. In this work, time-dependent simulations are performed to determine the nonlinear evolution of finite amplitude disturbances in the post-critical flow regime. It is shown that a nonlinear transition, which is facilitated by a supercritical pitchfork bifurcation, establishes a finite amplitude state (FAS) in which the average polymer stretch is highly modulated. The maximum normal stress, observed at the channel nip, can increase by up to approximately 100% when the Weissenberg number, defined as the ratio of the fluid relaxation time to an inverse characteristic shear rate, is increased by only 10% beyond its critical value. This is attributed to the amplification of configurational perturbations by the base flow shear rate, which attains its maximum at the channel nip. The effect of finite chain extensibility on the critical condition and nonlinear instability is investigated using the FENE-CR model. The stabilizing effect of finite extensibility can be expressed through a renormalization of the Weissenberg number by accounting for the screening effect of the nonlinear force law on the transmission of configurational perturbations to polymeric stress. The principal features of the FAS are qualitatively model-independent. The FAS exhibits a small, but numerically perceptible increase in the friction factor as compared to the base flow. The implication of the findings on the experimentally observed flow resistance enhancement phenomenon in viscoelastic creeping flows through converging/diverging geometries is discussed.     

Viscoelastic stresses in the stagnation flow of a dilute polymer solution

In this paper, we consider viscoelastic stresses T₁₁, T₁₂ and T₂₂ arising in the stagnation flow of a dilute polymer solution; in particular, we consider an upper convected Maxwell (UCM) fluid. We present exact solutions to the coupled partial differential equations describing the viscoelastic stresses and deduce the results for the stress T₂₂ of Becherer et al. [P. Becherer, A.N. Morozov, W. van Saarloos, Scaling of singular structures in extensional flow of dilute polymer solutions, J. Non-Newtonian Fluid Mech. 153 (2008) 183–190]. As we considered the viscoelastic stresses over two spatial variables, we are able to study the effect of variable boundary data at the inflow. As such, our results are applicable to a wider range of fluid flow problems.     

An extended FENE dumbbell model theory for concentration dependent shear-induced anisotropy in dilute polymer solutions: addenda

Schneggenburger et al. [C. Schneggenburger, M. Kröger, S. Hess, An extended FENE dumbbell theory for concentration dependent shear-induced anisotropy in dilute polymer solutions, J. Non-Newtonian Fluid Mech. 62 (1996) 235] extended the original FENE dumbbell kinetic theory to describe concentration dependent shear-induced anisotropy in dilute polymer solutions by a mean-field approach. Besides providing an erratum to the above-mentioned paper and two revised figures we present related analytic results for steady shear and uniaxial elongational flow. Within the same framework we further consider a modified FENE potential and briefly discuss its implications.     

A stable unstructured finite volume method for parallel large-scale viscoelastic fluid flow calculations

A new stable unstructured finite volume method is presented for parallel large-scale simulation of viscoelastic fluid flows. The numerical method is based on the side-centered finite volume method where the velocity vector components are defined at the mid-point of each cell face, while the pressure term and the extra stress tensor are defined at element centroids. The present arrangement of the primitive variables leads to a stable numerical scheme and it does not require any ad-hoc modifications in order to enhance the pressure–velocity–stress coupling. The log-conformation representation proposed in [R. Fattal, R. Kupferman, Constitutive laws for the matrix–logarithm of the conformation tensor, J. Non-Newtonian Fluid Mech. 123 (2004) 281–285] has been implemented in order improve the limiting Weissenberg numbers in the proposed finite volume method. The time stepping algorithm used decouples the calculation of the polymeric stress by solution of a hyperbolic constitutive equation from the evolution of the velocity and pressure fields by solution of a generalized Stokes problem. The resulting algebraic linear systems are solved using the FGMRES(m) Krylov iterative method with the restricted additive Schwarz preconditioner for the extra stress tensor and the geometric non-nested multilevel preconditioner for the Stokes system. The implementation of the preconditioned iterative solvers is based on the PETSc library for improving the efficiency of the parallel code. The present numerical algorithm is validated for the Kovasznay flow, the flow of an Oldroyd-B fluid past a confined circular cylinder in a channel and the three-dimensional flow of an Oldroyd-B fluid around a rigid sphere falling in a cylindrical tube. Parallel large-scale calculations are presented up to 523,094 quadrilateral elements in two-dimension and 1,190,376 hexahedral elements in three-dimension.     

Stress communication and filtering of viscoelastic layers in oscillatory shear

We revisit a classical topic: response functions of viscoelastic layers in large amplitude oscillatory shear. Motivated by questions concerning protective biological layers, we focus on boundary stresses in a parallel plate geometry with imposed oscillatory strain or stress. These features are gleaned from resolution and analysis of coupled standing waves of deformation and stress. We identify a robust non-monotone variation in boundary stress signals with respect to all experimental controls: viscoelastic moduli of the layer, layer thickness, and driving frequency. This structure of peaks and valleys in boundary values of shear and normal stress indicates redundant mechanisms for stress communication (by tuning to the peaks) and stress filtering (by tuning to the valleys). In this paper, we first restrict to a single-mode non-linear Maxwell model for the viscoelastic layer, where analysis renders a transparent explanation of the phenomena. We then consider a Giesekus constitutive model of the layer, where analysis is supplanted by numerical simulations of coupled non-linear partial differential equations. Parametric studies of wall stress values from standing waves confirm persistence of the Maxwell model phenomena. The analysis and simulations rely on and extend our recent studies of shear waves in a micro parallel plate rheometer [S.M. Mitran, M.G. Forest, L. Yao, B. Lindley, D. Hill, Extenstions of the Ferry shear wave model for active linear and nonlinear microrheology, J. Non-Newtonian Fluid Mech. 154 (2008) 120–135; D.B. Hill, B. Lindley, M.G. Forest, S.M. Mitran, R. Superfine, Experimental and modeling protocols from a micro-parallel plate rheometer, UNC Preprint, 2008].     

Numerical stability of the method of Brownian configuration fields

We investigate numerical aspects of the Brownian configuration fields method, and in particular its numerical stability as the Weissenberg number increases. Our results show the method to be immune to the type of instability leading to numerical blowup in the simulation of macroscopic models. We discuss this finding in the light of the stability criterion proposed in Fattal et al. [R. Fattal, R. Kupferman, Time-dependent simulation of viscoelastic flows at high Weissenberg using the log-conformation representation, J. Non Newtonian Fluid Mech. 126 (2005) 23–37].     

Droplet dynamics in sub-critical complex flows

A newly designed eccentric cylinder device has been used to study the deformation and orientation of single Newtonian droplets immersed in an immiscible Newtonian liquid in a controlled complex flow field. Optical microscopy coupled with image acquisition analysis allows monitoring the dynamics of droplets flowing in the gap between the eccentric cylinders. Throughout the experiments, the flow intensity was kept below the critical conditions for droplet break-up. The experimental results are compared with predictions which are obtained using the transient form of the phenomenological model of Maffettone and Minale (J Non-Newtonian Fluid Mech 78:227–241, 1998; J Non-Newtonian Fluid Mech 84:105–106, 1999), incorporating a flow type parameter that accounts for the relative amount of elongational effects in the flow field and adapting the capillary number to mixed flows. For all the sub-critical flows studied here, good agreement was found between model predictions and experimental data, providing, for the first time, a quantitative assessment of drop shape predictions in complex flows.     

Thermodynamic model formulation for viscoplastic solids as general equations for non-equilibrium reversible–irreversible coupling

Thermodynamic models for viscoplastic solids are often formulated in the context of continuum thermodynamics and the dissipation principle. The purpose of the current work is to show that models for such material behavior can also be formulated in the form of a General Equation for Non-Equilibrium Reversible–Irreversible Coupling (GENERIC), see, e.g., Grmela and Öttinger (Phys Rev E, 56:6620–6632, 1997), Öttinger and Grmela (Phys Rev E, 56:6633–6655, 1997), Grmela (J Non-Newtonian Fluid Mech, 165:980–986, 2010). A GENERIC combines Hamiltonian-dynamics-based modeling of time-reversible processes with Onsager–Casimir-based modeling of time-irreversible processes. The result is a model for the approach of non-equilibrium systems to thermodynamic equilibrium. Originally developed to model complex fluids, it has recently been applied to anisotropic inelastic solids in Eulerian (Hütter and Tervoort, in J Non-Newtonian Fluid Mech, 152:45–52, 2008; Hütter and Tervoort, in J Non-Newtonian Fluid Mech, 152:53–65, 2008; Hütter and Tervoort, in Adv Appl Mech, 42:254–317, 2008) and Lagrangian (Hütter and Svendsen, in J Elast 104:357–368, 2011) settings, as well as to damage mechanics. For simplicity, attention is focused in the current work on the case of thermoelastic viscoplasticity. Central to this formulation is a GENERIC-based form for the viscoplastic flow rule. A detailed comparison with the formulation based on continuum thermodynamics and the dissipation principle is given.     

Existence and Stability of Solutions for Steady Flows of Fibre Suspension Flows

We establish existence, uniqueness, convergence and stability of solutions to the equations of steady flows of fibre suspension flows. The existence of a unique steady solution is proven by using an iterative scheme. One of the restrictions imposed on the data confirms a well known fact proven in Galdi and Reddy (J Non-Newtonian Fluid Mech 83:205–230, 1999), Munganga and Reddy (Math Models Methods Appl Sci 12:1177–1203, 2002) and Munganga et al. (J Non-Newtonian fluid Mech 92:135–150, 2000) that the particle number N _p must be less than 35/2. Exact solutions are calculated for Couette and Poiseuille flows. Solutions of Poiseuille flows are shown to be more accurate than those of Couette flow when a time perturbation is considered.     

Corner flow of a suspension of rigid rods

Stress singularities in the neighbourhood of sharp corners can be a source of severe problems in the numerical simulation of non-Newtonian flows leading to loss of convergence with grid refinement (G.G. Lipscombe, R. Kennings and M.M. Denn, J. Non-Newtonian Fluid Mech., 24 (1987) 85 [1]). For Newtonian flows the nature of this singularity is given by the analysis of Dean and Montagnon (W.R. Dean and P.E. Montagnon, Phil. Trans. R. Soc. London, Ser. A., 308 (1949) 199 [2]) in terms of similarity solutions. In this paper we extend this similarity analysis to a suspension of rigid rods. In the limit of nearly full extension the FENE constitutive model has the same behaviour as such a suspension. Our analysis predicts the possibility of lip vortices but their behaviour is somewhat inconsistent with those observed experimentally.