TY - JOUR
T1 - Finite element analysis of stability of two-dimensional viscoelastic flows to three-dimensional perturbations
AU - Smith, M. D.
AU - Armstrong, R. C.
AU - Brown, R. A.
AU - Sureshkumar, R.
N1 - Funding Information:
We would like to thank Howard E. Covert for many helpful conversations regarding the Arnoldi method and Usamah Ahmad Al-Mubaiyedh for eigenspectra for the Couette problem. This work was supported primarily by the ERC Program of the National Science Foundation under Award Number EEC-9731680.
PY - 2000/10
Y1 - 2000/10
N2 - We present numerical methods for analysis of the stability of two-dimensional steady viscoelastic flows to small-amplitude, two-dimensional and three-dimensional disturbances based on finite element calculations of the steady base flow and the perturbation. Direct time integration of the linearized equations of motion and iterative calculation of the most dangerous components of the eigenspectrum are tested. Finite element discretizations based on the DEVSS-G finite element discretization with Newton's method used to compute steady-state solutions. Two different time integration schemes are tested for computing the time evolution of general, random disturbances: a θ-method operator-splitting scheme and a fourth-order Runge-Kutta method. For both time integrators, time stepping is decoupled into a solution of a modified Stokes problem and an evaluation of the time-dependent constitutive equation. The overall efficiency of both methods is extremely high, as is the potential for implementation on parallel computers. An algorithm also is presented for calculating eigenvalues with the largest real part that combines time integration of the linearized equations with a Krylov subspace method to accelerate the calculation of the eigenvalues. Although this method does not dramatically reduce the computational cost over use of time integration alone, it does provide a more complete analysis of the eigenspectrum. For both direct time integration and the hybrid time integration/Krylov calculation, the stability results for cylindrical Couette flow show quantitative agreement with the eigenvalues calculated by using other methods of analysis [M. Avgousti, A.N. Beris, J. Non-Newtonian Fluid Mech. 50 (1993) 225-251]. Contrary to the results in our previous paper [R. Sureshkumar, M.D. Smith, R.C. Armstrong, R.A. Brown, J. Non-Newtonian Fluid Mech. 82 (1999) 57-104], we find that the flow of an Oldroyd-B fluid through a closely-spaced cylinder array is stable to two-dimensional perturbations. However, allowing the perturbations to be three-dimensional and considering an isolated cylinder does not alter the conclusions of our earlier study [R. Sureshkumar, M.D. Smith, R.C. Armstrong, R.A. Brown, J. Non-Newtonian Fluid Mech. 82 (1999) 57-104] of the two-dimensional stability of widely-spaced arrays of cylinders; the flow around an isolated cylinder is computed to be stable for all values of the Weissenberg number obtainable with these calculations, We≤0.75. (C) 2000 Elsevier Science B.V. All rights reserved.
AB - We present numerical methods for analysis of the stability of two-dimensional steady viscoelastic flows to small-amplitude, two-dimensional and three-dimensional disturbances based on finite element calculations of the steady base flow and the perturbation. Direct time integration of the linearized equations of motion and iterative calculation of the most dangerous components of the eigenspectrum are tested. Finite element discretizations based on the DEVSS-G finite element discretization with Newton's method used to compute steady-state solutions. Two different time integration schemes are tested for computing the time evolution of general, random disturbances: a θ-method operator-splitting scheme and a fourth-order Runge-Kutta method. For both time integrators, time stepping is decoupled into a solution of a modified Stokes problem and an evaluation of the time-dependent constitutive equation. The overall efficiency of both methods is extremely high, as is the potential for implementation on parallel computers. An algorithm also is presented for calculating eigenvalues with the largest real part that combines time integration of the linearized equations with a Krylov subspace method to accelerate the calculation of the eigenvalues. Although this method does not dramatically reduce the computational cost over use of time integration alone, it does provide a more complete analysis of the eigenspectrum. For both direct time integration and the hybrid time integration/Krylov calculation, the stability results for cylindrical Couette flow show quantitative agreement with the eigenvalues calculated by using other methods of analysis [M. Avgousti, A.N. Beris, J. Non-Newtonian Fluid Mech. 50 (1993) 225-251]. Contrary to the results in our previous paper [R. Sureshkumar, M.D. Smith, R.C. Armstrong, R.A. Brown, J. Non-Newtonian Fluid Mech. 82 (1999) 57-104], we find that the flow of an Oldroyd-B fluid through a closely-spaced cylinder array is stable to two-dimensional perturbations. However, allowing the perturbations to be three-dimensional and considering an isolated cylinder does not alter the conclusions of our earlier study [R. Sureshkumar, M.D. Smith, R.C. Armstrong, R.A. Brown, J. Non-Newtonian Fluid Mech. 82 (1999) 57-104] of the two-dimensional stability of widely-spaced arrays of cylinders; the flow around an isolated cylinder is computed to be stable for all values of the Weissenberg number obtainable with these calculations, We≤0.75. (C) 2000 Elsevier Science B.V. All rights reserved.
KW - Oldroyd-B fluid
KW - Runge-Kutta method
KW - Viscoelastic flows
UR - http://www.scopus.com/inward/record.url?scp=0034307617&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=0034307617&partnerID=8YFLogxK
U2 - 10.1016/S0377-0257(00)00124-5
DO - 10.1016/S0377-0257(00)00124-5
M3 - Article
AN - SCOPUS:0034307617
SN - 0377-0257
VL - 93
SP - 203
EP - 244
JO - Journal of Non-Newtonian Fluid Mechanics
JF - Journal of Non-Newtonian Fluid Mechanics
IS - 2-3
ER -