High Weissenberg number simulation of die swell for differential type viscoelastic model by the streamline-upwind finite element method

Yoshifumi Kuwano, Toshihisa Kajiwara, Kazumori Funatsu

Research output: Contribution to journalArticle

Abstract

In order to predict the die swell seen in the actual polymer processing, the planar, the capillary and the annular die swell simulations have been performed by the streamline-upwind finite element method with the subelements for stress components, which was shown effective to calculate up to high Weissenberg number (We) for the entry flow simulation in an earlier study. The calculation using the Giesekus model, which is the differential type viscoelastic model, was found feasible over hundreds of We in the planar and the capillary die swell simulations as long as the primary normal stress difference was not so large. The shape of free surface at high We under the condition of no gravitation once showed the maximum swell and became an equilibrium one after shrinking back a little. This tendency became remarkable for the model with larger We and larger primary normal stress difference. Through the examination of the velocity profile, it was found that the velocity near the free surface was accelerated during the swell after extrusion and was larger than the inside velocity in the neighborhood of the position which showed the maximum swell. Since the accelerated outer fluid dragged the inside fluid, the swelling ratio was supposed to take an equilibrium value after shrinking back a little. Also, as the primary normal stress difference became larger, the axial position of the maximum swell approached the die. This may be due to that the model with the large primary normal stress difference predicts faster swell, because the elastic recovery force after extrusion is large. The axial position of the maximum swell became distant from the die with increasing We. We interpreted that the axial position of the maximum swell shifted downstream as the representative relaxation time was longer, or the velocity became larger. On the other hand, the calculation became impossible for We > 3 in the annular die swell simulation. In order to examine this reason, we performed the calculation for the planar die swell in two different analysis regions, i.e., one being the whole flow region with two singular points and the other, the half flow region with one singular point in consideration of flow symmetry. The calculation in the whole flow region with two singular points was unsuccessful at high We. It seems that the presence of two singular points in the analysis region made it impossible to perform the annular die swell up to high We simulation.

Original languageEnglish
Pages (from-to)755-756
Number of pages2
JournalJournal of Rheology
Volume38
Issue number3
DOIs
Publication statusPublished - Jan 1 1994

Fingerprint

finite element method
Finite element method
simulation
Extrusion
Fluids
fluids
Flow simulation
entry
swelling
Relaxation time
Swelling
Gravitation
Polymers
tendencies
velocity distribution
examination
relaxation time
recovery
gravitation
Recovery

All Science Journal Classification (ASJC) codes

  • Materials Science(all)
  • Condensed Matter Physics
  • Mechanics of Materials
  • Mechanical Engineering

Cite this

High Weissenberg number simulation of die swell for differential type viscoelastic model by the streamline-upwind finite element method. / Kuwano, Yoshifumi; Kajiwara, Toshihisa; Funatsu, Kazumori.

In: Journal of Rheology, Vol. 38, No. 3, 01.01.1994, p. 755-756.

Research output: Contribution to journalArticle

@article{47756042bf2a45889ca22cd48df7016e,
title = "High Weissenberg number simulation of die swell for differential type viscoelastic model by the streamline-upwind finite element method",
abstract = "In order to predict the die swell seen in the actual polymer processing, the planar, the capillary and the annular die swell simulations have been performed by the streamline-upwind finite element method with the subelements for stress components, which was shown effective to calculate up to high Weissenberg number (We) for the entry flow simulation in an earlier study. The calculation using the Giesekus model, which is the differential type viscoelastic model, was found feasible over hundreds of We in the planar and the capillary die swell simulations as long as the primary normal stress difference was not so large. The shape of free surface at high We under the condition of no gravitation once showed the maximum swell and became an equilibrium one after shrinking back a little. This tendency became remarkable for the model with larger We and larger primary normal stress difference. Through the examination of the velocity profile, it was found that the velocity near the free surface was accelerated during the swell after extrusion and was larger than the inside velocity in the neighborhood of the position which showed the maximum swell. Since the accelerated outer fluid dragged the inside fluid, the swelling ratio was supposed to take an equilibrium value after shrinking back a little. Also, as the primary normal stress difference became larger, the axial position of the maximum swell approached the die. This may be due to that the model with the large primary normal stress difference predicts faster swell, because the elastic recovery force after extrusion is large. The axial position of the maximum swell became distant from the die with increasing We. We interpreted that the axial position of the maximum swell shifted downstream as the representative relaxation time was longer, or the velocity became larger. On the other hand, the calculation became impossible for We > 3 in the annular die swell simulation. In order to examine this reason, we performed the calculation for the planar die swell in two different analysis regions, i.e., one being the whole flow region with two singular points and the other, the half flow region with one singular point in consideration of flow symmetry. The calculation in the whole flow region with two singular points was unsuccessful at high We. It seems that the presence of two singular points in the analysis region made it impossible to perform the annular die swell up to high We simulation.",
author = "Yoshifumi Kuwano and Toshihisa Kajiwara and Kazumori Funatsu",
year = "1994",
month = "1",
day = "1",
doi = "10.1122/1.550571",
language = "English",
volume = "38",
pages = "755--756",
journal = "Journal of Rheology",
issn = "0148-6055",
publisher = "Society of Rheology",
number = "3",

}

TY - JOUR

T1 - High Weissenberg number simulation of die swell for differential type viscoelastic model by the streamline-upwind finite element method

AU - Kuwano, Yoshifumi

AU - Kajiwara, Toshihisa

AU - Funatsu, Kazumori

PY - 1994/1/1

Y1 - 1994/1/1

N2 - In order to predict the die swell seen in the actual polymer processing, the planar, the capillary and the annular die swell simulations have been performed by the streamline-upwind finite element method with the subelements for stress components, which was shown effective to calculate up to high Weissenberg number (We) for the entry flow simulation in an earlier study. The calculation using the Giesekus model, which is the differential type viscoelastic model, was found feasible over hundreds of We in the planar and the capillary die swell simulations as long as the primary normal stress difference was not so large. The shape of free surface at high We under the condition of no gravitation once showed the maximum swell and became an equilibrium one after shrinking back a little. This tendency became remarkable for the model with larger We and larger primary normal stress difference. Through the examination of the velocity profile, it was found that the velocity near the free surface was accelerated during the swell after extrusion and was larger than the inside velocity in the neighborhood of the position which showed the maximum swell. Since the accelerated outer fluid dragged the inside fluid, the swelling ratio was supposed to take an equilibrium value after shrinking back a little. Also, as the primary normal stress difference became larger, the axial position of the maximum swell approached the die. This may be due to that the model with the large primary normal stress difference predicts faster swell, because the elastic recovery force after extrusion is large. The axial position of the maximum swell became distant from the die with increasing We. We interpreted that the axial position of the maximum swell shifted downstream as the representative relaxation time was longer, or the velocity became larger. On the other hand, the calculation became impossible for We > 3 in the annular die swell simulation. In order to examine this reason, we performed the calculation for the planar die swell in two different analysis regions, i.e., one being the whole flow region with two singular points and the other, the half flow region with one singular point in consideration of flow symmetry. The calculation in the whole flow region with two singular points was unsuccessful at high We. It seems that the presence of two singular points in the analysis region made it impossible to perform the annular die swell up to high We simulation.

AB - In order to predict the die swell seen in the actual polymer processing, the planar, the capillary and the annular die swell simulations have been performed by the streamline-upwind finite element method with the subelements for stress components, which was shown effective to calculate up to high Weissenberg number (We) for the entry flow simulation in an earlier study. The calculation using the Giesekus model, which is the differential type viscoelastic model, was found feasible over hundreds of We in the planar and the capillary die swell simulations as long as the primary normal stress difference was not so large. The shape of free surface at high We under the condition of no gravitation once showed the maximum swell and became an equilibrium one after shrinking back a little. This tendency became remarkable for the model with larger We and larger primary normal stress difference. Through the examination of the velocity profile, it was found that the velocity near the free surface was accelerated during the swell after extrusion and was larger than the inside velocity in the neighborhood of the position which showed the maximum swell. Since the accelerated outer fluid dragged the inside fluid, the swelling ratio was supposed to take an equilibrium value after shrinking back a little. Also, as the primary normal stress difference became larger, the axial position of the maximum swell approached the die. This may be due to that the model with the large primary normal stress difference predicts faster swell, because the elastic recovery force after extrusion is large. The axial position of the maximum swell became distant from the die with increasing We. We interpreted that the axial position of the maximum swell shifted downstream as the representative relaxation time was longer, or the velocity became larger. On the other hand, the calculation became impossible for We > 3 in the annular die swell simulation. In order to examine this reason, we performed the calculation for the planar die swell in two different analysis regions, i.e., one being the whole flow region with two singular points and the other, the half flow region with one singular point in consideration of flow symmetry. The calculation in the whole flow region with two singular points was unsuccessful at high We. It seems that the presence of two singular points in the analysis region made it impossible to perform the annular die swell up to high We simulation.

UR - http://www.scopus.com/inward/record.url?scp=84957276272&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84957276272&partnerID=8YFLogxK

U2 - 10.1122/1.550571

DO - 10.1122/1.550571

M3 - Article

VL - 38

SP - 755

EP - 756

JO - Journal of Rheology

JF - Journal of Rheology

SN - 0148-6055

IS - 3

ER -