Numerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries

  • Published : 2007.12.31

Abstract

Here we demonstrate complex transient behavior of viscoelastic liquid described numerically with the Leonov model in straight and contraction channel flow domains. Finite element and implicit Euler time integration methods are employed for spatial discretization and time marching. In order to stabilize the computational procedure, the tensor-logarithmic formulation of the constitutive equation with SUPG and DEVSS algorithms is implemented. For completeness of numerical formulation, the so called traction boundaries are assigned for flow inlet and outlet boundaries. At the inlet, finite traction force in the flow direction with stress free condition is allocated whereas the traction free boundary is assigned at the outlet. The numerical result has illustrated severe forward-backward fluctuations of overall flow rate in inertial straight channel flow ultimately followed by steady state of forward flow. When the flow reversal occurs, the flow patterns exhibit quite complicated time variation of streamlines. In the inertialess flow, it takes much more time to reach the steady state in the contraction flow than in the straight pipe flow. Even in the inertialess case during startup contraction flow, quite distinctly altering flow patterns with the lapse of time have been observed such as appearing and vanishing of lip vortices, coexistence of multiple vortices at the contraction comer and their merging into one.

Keywords

References

  1. Fattal, R. and R. Kupferman, 2004, Constitutive laws of the matrix-logarithm of the conformation tensor, J. Non-Newtonian Fluid Mech. 123, 281-285 https://doi.org/10.1016/j.jnnfm.2004.08.008
  2. Fattal, R. and R. Kupferman, 2005, Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation, J. Non-Newtonian Fluid Mech. 126, 23-37 https://doi.org/10.1016/j.jnnfm.2004.12.003
  3. Gupta, M., 1997, Viscoelastic modeling of entrance flow using multimode Leonov model, Int. J. Numer. Meth. Fluids 24, 493-517 https://doi.org/10.1002/(SICI)1097-0363(19970315)24:5<493::AID-FLD502>3.0.CO;2-W
  4. Hulsen, M. A., 2004, Keynote presentation in Internatioan Congress on Rheology 2004, Seoul, Korea
  5. Hulsen, M. A., R. Fattal and R. Kupferman, 2005, Flow of viscoelastic fluids past a cylinder at high Weissenberg number: Stabilized simulations using matrix logarithms, J. Non-Newtonian Fluid Mech. 127, 27-39 https://doi.org/10.1016/j.jnnfm.2005.01.002
  6. Kwon, Y., 2004, Finite element analysis of planar 4:1 contraction flow with the tensor-logarithmic formulation of differential constitutive equations, Korea-Australia Rheology J. 16, 183-191
  7. Leonov, A. I., 1976, Nonequilibrium thermodynamics and rheology of viscoelastic polymer media, Rheol. Acta 15, 85-98 https://doi.org/10.1007/BF01517499
  8. Sato, T. and M. Richardson, 1994, Explicit numerical simulation of time-dependent viscoelastic flow problems by a finite element/finite volume method, J. Non-Newtonian Fluid Mech. 51, 249-275 https://doi.org/10.1016/0377-0257(94)85019-4
  9. Webster, M. F., H. R. Tamaddon-Jahromi and M. Aboubacar, 2004, Transient viscoelastic flows in planar contractions, J. Non-Newtonian Fluid Mech. 118, 83-101 https://doi.org/10.1016/j.jnnfm.2004.03.001
  10. Yoon, S. and Y. Kwon, 2005, Finite element analysis of viscoelastic flows in a domain with geometric singularities, Korea-Australia Rheology J. 17, 99-110