A finite element method for liquid with memory

A. C. Papanastasiou, L. E. Scriven, Chris Macosko

Research output: Contribution to journalArticlepeer-review

36 Scopus citations


Steady, two-dimensional flows of viscoelastic liquid with memory are analyzed by means of a nonlinear integral constitutive equation, novel streamlined finite elements, Galerkin's method of weighted residuals, and Newton iteration. By their relation to streamlines the new elements allow, for the first time, full Newton iteration of the algebraic equation set to which the governing integrodifferential system reduces. Streamlines are computed simultaneously with velocity components and pressure, the primary unknowns; to track the history of liquid particles the deformation equations are solved analytically in a Protean system of coordinates that conforms to the streamlines. In the cases studied the Newton iteration converges quadratically to a creeping flow state of Weissenberg number (product of upstream wall shear rate and relaxation time of the liquid) up to about 20 in channel flow, 10 in film flow, and 2 in die-swell flow. Shear-thinning shrinks the domain of convergence; it also reduces die-swell. A slip boundary condition in the vicinity of the contact line extends the domain of convergence.

Original languageEnglish (US)
Pages (from-to)271-288
Number of pages18
JournalJournal of Non-Newtonian Fluid Mechanics
Issue number3
StatePublished - 1987

Fingerprint Dive into the research topics of 'A finite element method for liquid with memory'. Together they form a unique fingerprint.

Cite this