Yes / A recent development of the Molecular Stress Function constitutive model, the Hierarchical Multi-Mode Molecular Stress Function (HMMSF) model has been shown to fit a large range of rheometrical data with accuracy, for a large range of polymer melts. We develop a 3D simulation of the HMMSF model and compare it to experimental data for the flow of Lupolen 1840H LDPE through an abrupt 3D contraction flow. We believe this to be the first finite element implementation of the HMMSF model. It is shown that the model gives a striking agreement with experimental vortex opening angles, with very good agreement to full-field birefringence measurements, over a wide range of flow rates. A method to give fully-developed inlet boundary conditions is implemented (in place of using parabolic inlet boundary conditions), which gives a significantly improved match to birefringence measurements in the inlet area, and in low stress areas downstream from the inlet. Alternative constitutive model parameters are assessed following the principle that extensional rheometer data actually provides a ‘lower bound’ for peak extensional viscosity. It is shown that the model robustly maintains an accurate fit to vortex opening angle and full-field birefringence data, provided that both adjustable parameters are kept such that both shear and extensional data are well fitted.

Yes / The Molecular Stress Function (MSF) constitutive model with convective constraint release mechanism has been shown to accurately fit a large range of viscometric data, and also shown to give strong vortex growth in flows of LDPE through planar and axisymmetric contractions. This work compares simulation and experimental results for 3-D flows of Lupolen 1840H LDPE through a contraction slit; 3-D effects are introduced by using a slit with a low upstream aspect ratio of 5:3. Comparisons are made with vortex opening angles obtained from streak photography, and also with stress birefringence measurements. The comparisons are made with two versions of the convective constraint release (CCR) mechanism. The simulated vortex angles for one version of the CCR mechanism are found to approach what is seen experimentally. The best-fit value for the stress optical coefficient was found to vary between CCRs and to decrease with flow rate. This is partially explained by different centreline elongational rates with the two CCRs, which in turn is related to different opening angles. A 3-D simulation is compared to the corresponding 2-D simulation. It is shown that both velocity vectors and birefringence show only small changes to around 60% of the distance to the side wall.