Numerical simulation of a viscoelastic fluid with anisotropic heat conduction

For the nonisothermal flow of a viscoelastic fluid we have taken into account temperature dependency of the relaxation times and the viscosities in the constitutive equation for the stress. In the energy equation the heat flux is specified by Fourier's law, where anisotropic heat conduction has been taken into account. Furthermore one has to specify which part of the stress work is dissipated and which part is stored as elastic energy. The equations are solved with a finite element method for the balance equations and a streamline integration method for the constitutive equation. The influence of the Deborah number, the Péclet number and the cooling temperature are examined in a flow through a 4 to 1 contraction

Response to "Comment on: `Thermodynamics of viscoelastic fluids: the temperature equation'"

No abstrac

A modification of the convective constraint release mechanism in the molecular stress function model giving enhanced vortex growth

The molecular stress function model with convective constraint release (MSF with CCR) constitutive model [J. Rheol. 45 (2001), 1387] is capable of fitting all viscometric data for IUPAC LDPE, with only two adjustable parameters (with difference found only on reported ¿steady-state¿ elongational viscosities). The full MSF with CCR model is implemented in a backwards particle-tracking implementation, using an adaptive method for the computation of relative stretch that reduces simulation time many-fold, with insignificant loss of accuracy. The model is shown to give improved results over earlier versions of the MSF (without CCR) when compared to well-known experimental data from White and Kondo [J. non-Newt. Fluid Mech., 3 (1977), 41]; but still to under-predict contraction flow opening angles. The discrepancy is traced to the interaction between the rotational dissipative function and the large stretch levels caused by the contraction flow. A modified combination of dissipative functions in the constraint release mechanism is proposed, which aims to reduce this interaction to allow greater strain hardening in a mixed flow. The modified constraint release mechanism is shown to fit viscometric rheological data equally well, but to give opening angles in the complex contraction flow that are much closer to the experimental data from White and Kondo. It is shown (we believe for the first time) that a constitutive model demonstrates an accurate fit to all planar elongational, uniaxial elongational and shear viscometric data, with a simultaneous agreement with this well-known experimental opening angle data. The sensitivity of results to inaccuracies caused by representing the components of the deformation gradient tensor to finite precision is examined; results are found to be insensitive to even large reductions in the precision used for the representation of components. It is shown that two models that give identical response in elongational flow, and a very similar fit to available shear data, give significantly different results in flows containing a mix of deformation modes. The implication for constitutive models is that evaluation against mixed deformation mode flow data is desirable in addition to evaluation against viscometric measurements

Stress singularities and the formation of birefringent strands in stagnation flows of dilute polymer solutions

We consider stagnation point flow away from a wall for creeping flow of dilute polymer solutions. For a simplified flow geometry, we explicitly show that a narrow region of strong polymer extension (a birefringent strand) forms downstream of the stagnation point in the UCM model and extensions, like the FENE-P model. These strands are associated with the existence of an essential singularity in the stresses, which is induced by the fact that the stagnation point makes the convective term in the constitutive equation into a singular point. We argue that the mechanism is quite general, so that all flows that have a separatrix going away from the stagnation point exhibit some singular behaviour. These findings are the counterpart for wall stagnation points of the recently discovered singular behaviour in purely elongational flows: the underlying mechanism is the same while the different nature of the singular stress behaviour reflects the different form of the velocity expansion close to the stagnation point.Comment: 15 pages, 6 figure