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 -