Instationary Eulerian viscoelastic flow simulations using time separable Rivlin–Sawyers constitutive equations

• Published in: Journal of non-Newtonian fluid mechanics Vol. 89; no. 1; pp. 209 - 228
• Main Authors: Peters, E.A.J.F., Hulsen, M.A., van den Brule, B.H.A.A.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.02.2000; Elsevier
• Subjects: Computational methods in fluid dynamics; Cross-disciplinary physics: materials science; rheology; DEVSS/DG; Engineering Sciences; Eulerian description; Exact sciences and technology; Falling sphere in tube; Fluid dynamics; Fundamental areas of phenomenology (including applications); Integral models; Mechanics; Non-newtonian fluid flows; Physics; Rheology; Structural mechanics; Techniques and apparatus; Transient viscoelastic flow; Viscosity measurements; Eulerian description; Falling sphere in tube; DEVSS/DG; Integral models; Transient viscoelastic flow; Constitutive equation; Viscoelastic fluid; Falling sphere viscosimeter; Refinement method; Digital simulation; Maxwell fluid; Unsteady flow; Euler scheme; Drag coefficient; Mesh generation; Discretization method
• ISSN: 0377-0257; 1873-2631
• DOI: 10.1016/S0377-0257(99)00026-9

A time dependent method for solving integral constitutive equations of the Rivlin–Sawyers type is introduced. The deformation history is represented by a finite number of deformation fields. Using these fields the stress integral is approximated as a finite sum. When the flow evolves the deformation fields are convected and deformed. The approach presented in this paper is the first Eulerian method that can handle integral equations in a time dependent way. The method is validated by using the upper-convected Maxwell (UCM) benchmark of a sphere moving in a tube. We show that the method converges with mesh and time step refinement and that the results are accurate, comparable to the results obtained with the differential equivalent of the UCM model. To demonstrate that complicated linear spectra are easily incorporated, results of a Rouse model simulation of 100 modes are presented. We also compare results on a falling sphere problem to the results obtained by a Lagrangian method as reported by Rasmussen and Hassager [H.K. Rasmussen, O. Hassager, On the sedimentation velocity of spheres in a polymeric liquid, Chem. Eng. Sci. 51 (1996) 1431–1440]. The model being employed is the PSM model, for which no differential equivalent exists.