Giesekus Constitutive Model for Thermoviscoelastic Fluids based on Ordered Rate Constitutive Theories

K. S. Surana; D. Nunez; J. N. Reddy

doi:10.6000/1929-5995.2013.02.04.5

Authors

K. S. Surana The University of Kansas
D. Nunez The University of Kansas
J. N. Reddy Texas A&M University

DOI:

https://doi.org/10.6000/1929-5995.2013.02.04.5

Keywords:

Contravariant, covariant, Jaumann, upper convective, lower convective, least squares

Abstract

This paper presents derivation of Giesekus constitutive model in Eulerian description based on ordered rate constitutive theories for thermoviscoelastic fluids for compressible and incompressible cases in contra-, co-variant and Jaumann bases. The ordered rate constitutive theories for thermoviscoelastic fluids of orders (m, n) consider convected time derivative of order m of the deviatoric Cauchy stress tensor in a chosen basis (i.e. co-, contra-variant or Jaumann) as dependent variable in the development of constitutive theories for the stress tensor. Its argument tensors consist of density, temperature, convected time derivatives of the deviatoric Cauchy stress tensor of up to order m-1 and convected time derivative of up to order n of the conjugate strain tensor. In addition, constitutive theory for the heat vector compatible with the constitutive theory for the deviatoric stress tensor is also presented in co-, contra-variant and Jaumann bases. It is shown that the Giesekus constitutive model is a subset of the rate constitutive theory of orders m = n = 1. It is also shown that the deviatoric Cauchy stress tensor (contra-, co-variant or Jaumann basis) naturally results as dependent variable in the constitutive theory, and that currently used Giesekus constitutive model in deviatoric polymer Cauchy stress tensor is not derivable based on axioms and principles of the constitutive theory in continuum mechanics. Numerical studies are presented for fully developed flow between parallel plates for a dense polymeric liquid using the Giesekus constitutive model derived in this paper as well as currently used model.

Author Biographies

K. S. Surana, The University of Kansas

Mechanical Engineering Department

Deane E. Ackers Distinguished Professor of Mechanical Engineering

D. Nunez, The University of Kansas

Department of Mechanical Engineering

Post-Doctoral Research Scientist

J. N. Reddy, Texas A&M University

Mechanical Engineering Department

Distinguished Professor and Wyatt Chair in Mechanical Engineering

References

Maxwell JC. On the dynamical theory of gases. Philos Trans R Soc Lond A 1867; A157: p. 49-88. DOI: https://doi.org/10.1098/rstl.1867.0004

Bird RB, Armstrong RC, Hassager O. Dynamics of polymeric liquids, Volume 1: fluid mechanics. 2nd ed. New York: Wiley 1987.

Bird RB, Armstrong RC, Hassager O. Dynamics of polymeric liquids, Volume 2: kinetic theory. 2nd ed. New York: Wiley 1987.

Oldroyd JG. On the formulation of rheological equations of state. Proc R Soc Lond A Math Phys Sci 1950; A200: p. 523-541. DOI: https://doi.org/10.1098/rspa.1950.0035

Giesekus H. A simple constitutive equation for polymer fluids based on the concept of deformationdependent tensorial mobility. J Nonnewton Fluid Mech 1982; 11: p. 69-110. DOI: https://doi.org/10.1016/0377-0257(82)85016-7

Phan-Thien N, Tanner RI. A nonlinear network viscoelastic model. J Rheol 1978; 22: p. 259-283. DOI: https://doi.org/10.1122/1.549481

Eringen AC. Mechanics of continua. New York: Wiley 1967.

Surana KS, Nunez D, Reddy JN, Romkes A. Rate constitutive theory for ordered thermofluids. Continuum Mech Thermodyn 2013; 25: p. 625-662. DOI: https://doi.org/10.1007/s00161-012-0257-6

Surana KS, Nunez D, Reddy JN, Romkes A. Rate constitutive theory for ordered thermoelastic solids. Ann solid struct mech 2012; 3: p. 27-54. DOI: https://doi.org/10.1007/s12356-012-0027-y

Brown R. A brief account of microscopical observations. Edinburgh J Sci 1828; 5: p. 358-371.

Brown R. Additional remarks on active molecules. Edinburgh J Sci 1829; 1: p. 314-319.

Udaykumar HS, Tran L, Belk DM, Vanden KJ. An Eulerian method for computation of multimaterial impact with ENO shock-capturing and sharp interfaces. J Comput Phys 2003; 186: p. 136-177. DOI: https://doi.org/10.1016/S0021-9991(03)00027-5

Tran L, Udaykumar HS. A particle-level set-based sharp interface cartesian grid method for impact, penetration, and void collapse. J Comput Phys 2004; 193: p. 469-510. DOI: https://doi.org/10.1016/j.jcp.2003.07.023

Rivlin RS, Ericksen JL. Stress-deformation relations for isotropic materials. Arch Ration Mech Anal 1955; 4: p. 323-425. DOI: https://doi.org/10.1512/iumj.1955.4.54011

Rivlin RS. Further remarks on the stressdeformation relations for isotropic materials. Arch Ration Mech Anal 1955; 4: p. 681-702. DOI: https://doi.org/10.1512/iumj.1955.4.54025

White FM. Fluid mechanics. 7th ed. New York: McGraw-Hill 2010.

Surana KS, MA Y, Reddy JN, Romkes A. The rate constitutive equations and their validity for progressively increasing deformation. Mech Adv Mater Struc 2010; 17: p. 509-533. DOI: https://doi.org/10.1080/15376494.2010.509195

Giesekus H. Die Rheologische Zustandsgleichung Elasto-Viskoser Fl ¨ussigkeiten - Insbesondere Von Weissenberg-Fl¨ussigkeiten - f ¨ ur Allgemeine und Station ¨are Fließvorg¨ange. J Appl Math Mech 1962; 42: p. 32-61. DOI: https://doi.org/10.1002/zamm.19620420105

Quinzani LM, Armstrong RC, Brown RA. Use of coupled birefringence and LDV studies of flow through a planar contraction to test constitutive equations for concentrated polymer solutions. J Rheol 1995; 39: p. 1201-1227. DOI: https://doi.org/10.1122/1.550725

Bell B, Surana KS. The p-version least squaresfinite element formulation for two-dimensional, incompressible, non-Newtonian isothermal and nonisothermal fluid flow. Int J Numer Methods Fluids 1994; 18: p. 127-162. DOI: https://doi.org/10.1002/fld.1650180202

Bell B, Surana KS. The p-version space-time coupled least squares finite element formulation for two-dimensional unsteady incompressible, Newtonian fluid flow. Adv Finite Elem Anal Fluid Dyn 1993; 171: p. 75-87.

Winterscheidt D, Surana KS. p-Version least squares finite element formulation for twodimensional, incompressible fluid flow. Int J Numer Methods Fluids 1994; 18: p. 43-69. DOI: https://doi.org/10.1002/fld.1650180104

Surana KS, Allu S, Tenpas PW, Reddy JN. The k- version of finite element method in gas dynamics: higher order global differentiability numerical solutions. Int J Numer Methods Eng 2006; 69: p. 1109- 1157. DOI: https://doi.org/10.1002/nme.1801

Surana KS, Reddy JN, Allu S. The k-version of finite element method for IVPs: mathematical and computational framework. Int J Comput Methods Eng Sci Mech 2007; 8: p. 123-136. DOI: https://doi.org/10.1080/15502280701252321

Surana KS, Allu S, Reddy JN, Tenpas PW. Least squares finite element processes in h; p; k mathematical and computational framework for a nonlinear conservation law. Int J Numer Methods Fluids 2008; 57: p. 1545-1568. DOI: https://doi.org/10.1002/fld.1695

Rajagopalan D, Armstrong RC, Brown RA. Finite element methods for calculation of steady viscoelastic flow using constitutive equations with a Newtonian viscosity. J Nonnewton Fluid Mech 1990; 36: p. 159-192. DOI: https://doi.org/10.1016/0377-0257(90)85008-M

Rajagopalan D, Phillips RJ, Armstrong RC, Brown RA, Bose A. The influence of viscoelasticity on the existence of steady solutions in two-dimensional rimming flow. J Fluid Mech 1992; 235: p. 611-642. DOI: https://doi.org/10.1017/S0022112092001241