Abstract
We present a numerical scheme for a previously unexploited formulation of the equations for unsteady viscoelastic flow. The formulation aligns the polymer stress along particle paths/streamlines, utilising the characteristic curves associated with the hyperbolic part of the constitutive equations. We illustrate the approach for the Oldroyd-B model in the benchmark 4:1 contraction for moderate elasticity numbers. We show that the scheme is able to accurately capture the re-entrant corner singularity for the polymer stresses and the pressure, the latter variable being inaccurately determined by schemes using the traditional formulation in terms of Cartesian polymer stresses. A space-step restriction for stability is derived, which can be numerically limiting in certain recirculation regions. This contrasts with the equivalent space-step restriction for the formulation in Cartesian stresses, which is limiting in flow regions of high velocity gradients, for example, at sharp corners in contraction flows.
Original language | English |
---|---|
Pages (from-to) | 462-489 |
Number of pages | 28 |
Journal | Journal of Computational Physics |
Volume | 388 |
Early online date | 23 Mar 2019 |
DOIs | |
Publication status | Published - 1 Jul 2019 |
Funding
The authors would like to thank the financial support given by FAPESP grants no. 2013/07375-0 , 2015/50094-7 , 2016/00456-2 and 2017/04471-9 , and The Royal Society Newton International Exchanges grant no. 2015/NI150225 . C.M. Oishi would like to thank the support provided by CNPq grant 307459/2016-0 . Appendix A
Keywords
- Natural Stress Formulation
- Numerical simulation
- Sharp corner flows
- Unsteady viscoelastic flows
ASJC Scopus subject areas
- Physics and Astronomy (miscellaneous)
- Computer Science Applications