TY - JOUR

T1 - Effect of artificial stress diffusivity on the stability of numerical calculations and the flow dynamics of time-dependent viscoelastic flows

AU - Sureshkumar, R.

AU - Beris, Antony N.

N1 - Funding Information:
We would like to acknowledge the financial support by NSF division of Fluid Mechanics and Particulate Flows under grant No. CTS-9114508. We are also grateful to the Pittsburgh Supercomputing Center for providing the computational resources used in this research.

PY - 1995/10

Y1 - 1995/10

N2 - In this work, we investigate the effect of the introduction of a stress diffusive term into the classical Oldroyd-B constitutive equation on the numerical stability of time-dependent viscoelastic flow calculations. The channel Poiseuille flow at Re ≫ 1 and O(1) We is chosen as a test problem. Through a linear stability analysis, we demonstrate that the introduction of a small amount of (dimensionless) diffusivity, typically of the order 10-3, does not affect the critical eigenmodes of the viscoelastic Orr-Sommerfeld problem appreciably. However, a diffusive term of that magnitude is shown to have a significant influence on the singular eigenmodes of the classical Oldroyd-B model, associated with the continuum spectra. A finite amplitude perturbation is constructed as a linear superposition of the eigenvectors corresponding to the most unstable eigenvalues of the problem. This is superimposed on the steady Poiseuille flow solution to provide the initial conditions for time-dependent simulations. The numerical algorithm involves a fully spectral spatial discretization and a semi-implicit second order integration in time. For the Oldroyd-B fluid, depending on the magnitude of the initial perturbation, numerical instabilities set in at relatively short times while the components of the conformation tensor increase monotonically in magnitude with time. Introduction of a diffusive term into this model is shown to stabilize the calculations remarkably, and for a three-dimensional simulation with Re = 5000 and We = 1, no instabilities were observed even at very large times. The effect of the magnitude of the diffusivity on the stability and the flow dynamics is addressed through a direct comparison of the results with those obtained for the Oldroyd-B model.

AB - In this work, we investigate the effect of the introduction of a stress diffusive term into the classical Oldroyd-B constitutive equation on the numerical stability of time-dependent viscoelastic flow calculations. The channel Poiseuille flow at Re ≫ 1 and O(1) We is chosen as a test problem. Through a linear stability analysis, we demonstrate that the introduction of a small amount of (dimensionless) diffusivity, typically of the order 10-3, does not affect the critical eigenmodes of the viscoelastic Orr-Sommerfeld problem appreciably. However, a diffusive term of that magnitude is shown to have a significant influence on the singular eigenmodes of the classical Oldroyd-B model, associated with the continuum spectra. A finite amplitude perturbation is constructed as a linear superposition of the eigenvectors corresponding to the most unstable eigenvalues of the problem. This is superimposed on the steady Poiseuille flow solution to provide the initial conditions for time-dependent simulations. The numerical algorithm involves a fully spectral spatial discretization and a semi-implicit second order integration in time. For the Oldroyd-B fluid, depending on the magnitude of the initial perturbation, numerical instabilities set in at relatively short times while the components of the conformation tensor increase monotonically in magnitude with time. Introduction of a diffusive term into this model is shown to stabilize the calculations remarkably, and for a three-dimensional simulation with Re = 5000 and We = 1, no instabilities were observed even at very large times. The effect of the magnitude of the diffusivity on the stability and the flow dynamics is addressed through a direct comparison of the results with those obtained for the Oldroyd-B model.

KW - Artificial stress diffusivity

KW - Flow dynamics

KW - Numerical calculations

KW - Time-dependent viscoelastic flow

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

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

U2 - 10.1016/0377-0257(95)01377-8

DO - 10.1016/0377-0257(95)01377-8

M3 - Article

AN - SCOPUS:0029379964

SN - 0377-0257

VL - 60

SP - 53

EP - 80

JO - Journal of Non-Newtonian Fluid Mechanics

JF - Journal of Non-Newtonian Fluid Mechanics

IS - 1

ER -