A pseudo-transient finite element method for the resolution of viscoelastic fluid flow problems by the method of characteristics

doi:10.1016/0377-0257(94)80074-X

Journal of Non-Newtonian Fluid Mechanics

Volume 55, Issue 3, December 1994, Pages 283-305

https://doi.org/10.1016/0377-0257(94)80074-X Get rights and content

Abstract

A pseudo-transient finite element method for the computation of two-dimensional viscoelastic flows using the method of characteristics is presented. The momentum and continuity equations are first solved as a saddle-point problem using the Uzawa algorithm with the viscoelastic stress treated as a fixed body force. The constitutive equation is then integrated with known kinematics using the method of characteristics. Numerical results involving the four-to-one abrupt circular contraction and the flow around a sphere in a tube are presented to assess the validity and the robustness of the proposed approach.

References (27)

F.G. Basombrio
J. Non-Newtonian Fluid Mech.
(1991)
J.M. Marchal et al.
J. Non-Newtonian Fluid Mech.
(1987)
D. Rajagopalan et al.
J. Non-Newtonian Fluid Mech.
(1990)
H. Giesekus
J. Non-Newtonian Fluid Mech.
(1982)
M.J. Crochet et al.
J. Non-Newtonian Fluid Mech.
(1982)
J.M. Grygiel et al.
Comp. Meth. Appl. Mech. Eng.
(1991)
P. Hansbo
Comp. Meth. Appl. Mech. Eng.
(1992)
A. Fortin et al.
J. Non-Newtonian Fluid Mech.
(1992)
M.J. Crochet et al.
J. Non-Newtonian Fluid Mech.
(1992)
T.J.R. Hughes et al.
Comp. Meth. Appl. Mech. Eng.
(1987)

R.B. King et al.

Comp. Meth. Appl. Mech. Eng.

(1987)

J.M. Winget et al.

Comp. Meth. Appl. Mech. Eng.

(1985)

J.P. Benqué et al.

Int. J. Numer. Meth. Fluids

(1988)

Cited by (9)

Stable Finite Element Discretizations for Viscoelastic Flow Models
2011, Handbook of Numerical Analysis
This chapter describes the way a general macroscopic viscoelastic fluid model can be reformulated, in terms of the conformation tensor, as a Riccati differential equation. It also discusses the way this reformulation can be used to establish the positive definiteness of the conformation tensor and key numerical methods based on the Eulerian–Lagrangian method, which discretizes the momentum equation and constitutive equations by solving the nonlinear ordinary differential equations that define the characteristics related to the transport part of the equation. The way the resulting discrete system can be effectively solved iteratively by combining multigrid and parallel computing techniques is discussed in the chapter. The chapter reviews the basic properties of the flow maps, the generalized Lie derivatives, and the algebraic Riccati differential equations. It introduces the connection between the algebraic Riccati differential equations and the macroscopic constitutive relations for viscoelastic fluids. It discusses the properties of various macroscopic models for viscoelastic fluids. The existing numerical schemes designed for simulating viscoelastic fluids at high Weissenberg number regimes is reviewed in the chapter.
A new finite element scheme using the Lagrangian framework for simulation of viscoelastic fluid flows
2002, Journal of Non-Newtonian Fluid Mechanics
A new numerical scheme for simulation of viscoelastic fluid flows was designed, making use of finite element algorithms generally regarded as advantageous for tackling the problem. This includes the Lagrangian approach for the solution of viscoelastic constitutive equation using the co-deformational frame of reference with a possibility of analytically solving the equation along the particles trajectories, which in turn allowed eluding the solution of any system of linear equations for the stress. Then, the full ellipticity of the momentum conservation equation was utilised thanks to a possibility of accurate determination of the stress tensor independently of the velocity field at the current stage of computation. The needed independent stress was calculated at each time step on the basis of the past deformation history, which in turn was determined on the basis of the past velocity fields, all incorporated into a modified Euler time stepping algorithm. Owing to explicit inclusion of the full viscous term from the viscoelastic model into the momentum conservation equation, no stress splitting was necessary. The trajectory feet tracking was done accurately using a semi-analytic solution of the displacement gradient evolution equation and a weak formulation of the kinematics equation, the latter at the expense of solving an extra symmetric system of linear equations.
The error expressed in the form of the Sobolev norms was determined using a comparison with available analytical solution for UCM fluid in the transient regime or numerically obtained steady-state stress values for the PTT fluid in Couette flow. The implementation of the PTT fluid model was done by modifying the relative displacement gradient tensor so that a new convective frame was defined.
The stability of the algorithm was assessed using the well-known benchmark problem of a sphere sedimenting in a tube with viscoelastic fluid. The stable numerical results were obtained at high Weissenberg numbers, with the limit of convergence Wi=6.6, exceeding any previously reported values. The robustness of the code was proven by simulation of the Weissenberg effect (the rod-climbing phenomenon) with the use of PTT fluid.
Segregated formulations in computational aspects of complex viscoelastic flows
1999, Rheology Series
This chapter examines the segregated formulations adopted in numerical simulations of viscoelastic flows and discusses the possibilities of emerging algorithms. It presents governing equations for viscoelastic flows in conjunction with basic decoupled techniques for the extra-stress. Several strategies handling the hyperbolic nature of the constitutive equations are discussed. The segregated formulations of the Navier-Stokes equation are treated and the possibilities of emerging segregated algorithms are considered. For an incompressible viscoelastic flow, the governing equations are the continuity, momentum and constitutive equations. The streamline upwind (SU) scheme for the control volume finite element method (SUCV) is proposed. When the FEM is considered in conjunction with the SU technique, it gives satisfactory results in analyzing viscoelastic flows, adopting SUCV into the existing algorithm based on the FVM seems to be shows potential.
Mixed finite element methods for viscoelastic flow analysis: A review
1998, Journal of Non-Newtonian Fluid Mechanics
The progress made during the past decade in the application of mixed finite element methods to solve viscoelastic flow problems using differential constitutive equations is reviewed. The algorithmic developments are discussed in detail. Starting with the classical mixed formulation, the elastic viscous stress splitting (EVSS) method as well as the related discrete EVSS and the so-called EVSS-G method are discussed among others. Furthermore, stabilization techniques such as the streamline upwind Petrov–Galerkin (SUPG) and the discontinuous Galerkin (DG) are reviewed. The performance of the numerical schemes for both smooth and non-smooth benchmark problems is discussed. Finally, the capabilities of viscoelastic flow solvers to predict experimental observations are reviewed.
A split Lagrangian-Eulerian method for simulating transient viscoelastic flows
1995, Journal of Non-Newtonian Fluid Mechanics
A novel numerical method for simulating time-dependent flow of viscoelastic fluids derived from dumbbell models is described. The constitutive equation is solved in a co-deforming frame, where the natural time-derivative is the upper-convected derivative. Mesh reconnection is achieved using a variant of Delaunay triangulation. The velocity and pressure are found via a finite element solution of the momentum equations. The method is tested by applying it to the benchmark problem of a sphere falling along the axis of a cylindrical tube.
Numerical simulation of the flow and fiber orientation in reinforced thermoplastic injection molded products
1997, International Polymer Processing

View all citing articles on Scopus

¹: Corresponding author.

²: CERSIM, Département de génie chimique. Université Lavai, Que., G1K 7P4, Canada.

View full text

A pseudo-transient finite element method for the resolution of viscoelastic fluid flow problems by the method of characteristics

Abstract

J. Non-Newtonian Fluid Mech.

J. Non-Newtonian Fluid Mech.

J. Non-Newtonian Fluid Mech.

J. Non-Newtonian Fluid Mech.

J. Non-Newtonian Fluid Mech.

Comp. Meth. Appl. Mech. Eng.

Comp. Meth. Appl. Mech. Eng.

J. Non-Newtonian Fluid Mech.

J. Non-Newtonian Fluid Mech.

Comp. Meth. Appl. Mech. Eng.

Comp. Meth. Appl. Mech. Eng.

Comp. Meth. Appl. Mech. Eng.

Int. J. Numer. Meth. Fluids