Abstract Numerical calculations of the flow of viscoelastic fluids through an abrupt contraction are prone to important errors in the evaluation of the stress and velocity fields near the re-entrant corner. We introduce a new type of rectangular finite element for calculating viscoelastic flows. Hermitian shape functions are used for interpolating the velocity field, while the stress shape functions are able to represent every derivative of the velocity field. The domain of convergence in Deborah number for the flow of an Oldroyd-B fluid is extended by a factor of four by comparison with earlier calculations. The numerical results show important differences between plane and axisymmetric flows in a four-to-one abrupt contraction.
[1]
K. Walters,et al.
Flow characteristics associated with abrupt changes in geometry in the case of highly elastic liquids
,
1986
.
[2]
D. V. Boger,et al.
Further observations of elastic effects in tubular entry flows
,
1986
.
[3]
A. Beris,et al.
Finite element calculation of viscoelastic flow in a journal bearing: II. Moderate eccentricity
,
1986
.
[4]
B. M. Fraeijs de Veubeke,et al.
A conforming finite element for plate bending
,
1968
.
[5]
M. Crochet,et al.
On the consequence of discretization errors in the numerical calculation of viscoelastic flow
,
1985
.