Skip to main content
SpringerLink
Log in

On the numerical stability of mixed finite-element methods for viscoelastic flows governed by differential constitutive equations

  • Published:
Theoretical and Computational Fluid Dynamics Aims and scope Submit manuscript

Abstract

Mixed finite-element methods for computation of viscoelastic flows governed by differential constitutive equations vary by the polynomial approximations used for the velocity, pressure, and stress fields, and by the weighted residual methods used to discretize the momentum, continuity, and constitutive equations. This paper focuses on computation of the linear stability of the planar Couette flow as a test of the numerical stability for solution of the upper-convected Maxwell model. Previous theoretical results prove this inertialess flow to be always stable, but that accurate calculation is difficult at high De because eigenvalues with fine spatial structure and high temporal frequency approach neutral stability. Computations with the much used biquadratic finite-element approximations for velocity and deviatoric stress and bilinear interpolation for pressure demonstrate numerical instability beyond a critical value of De for either the explicitly elliptic momentum equation (EEME) or elastic-viscous split-stress (EVSS) formulations, applying Galerkin's method for solution of the momentum and continuity equations, and using streamline upwind Petrov-Galerkin (SUPG) method for solution of the hyperbolic constitutive equation. The disturbance that causes the instability is concentrated near the stationary streamline of the base flow. The removal of this instability in a slightly modified form of the EEME formulation suggests that the instability results from coupling the approximations to the variables. A new mixed finite-element method, EVSS-G, is presented that includes smooth interpolation of the velocity gradients in the constitutive equation that is compatible with bilinear interpolation of the stress field. This formulation is tested with SUPG, streamline upwinding (SU), and Galerkin least squares (GLS) discretization of the constitutive equation. The EVSS-G/SUPG and EVSS-G/SU do not have the numerical instability described above; linear stability calculations for planar Couette flow are stable to values of De in excess of 50 and converge with mesh and time step. Calculations for the steady-state flow and its linear stability for a sphere falling in a tube demonstrate the appearance of linear instability to a time-periodic instability simultaneously with the apparent loss of existence of the steady-state solution. The instability appears as finely structured secondary cells that move from the front to the back of the sphere.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Similar content being viewed by others

References

  • Apelian, M.R., Armstrong, R.C., and Brown, R.A. (1988) Impact of constitutive equation and singularity on the calculation of stick-slip flow: the modified upper-convected Maxwell fluid. J. Non-Newtonian Fluid Mech., 27, 299–321.

    Google Scholar 

  • Baijens, F.T.P. (1992) Application of Galerkin-least-squares related methods to the mixed formulations of steady viscoelastic flow. J. Non-Newtonian Fluid Mech., submitted.

  • Baranger, J., and Sandri, D. (1991a) Approximation par finis d'ecoulements de fluides viscoelastiques: Existence de solutions approchees er majoration d'erreur. I. Constraintes discontinues, C. R. Acad. Sci. Paris Sér. I, 312, 541–544.

    Google Scholar 

  • Baranger, J., and Sandri, D. (1991b) Approximation par finis d'ecoulements de fluides viscoelastiques: Existence de solutions approchees er majoration d'erreur. II. Constraintes continues. C.R. Acad. Sci. Paris Sér. I, 313, 111–114.

    Google Scholar 

  • Beris, A.N., Armstrong, R.C., and Brown, R.A. (1987) Finite-element/spectral calculations of the flow of a Maxwell fluid between eccentric rotating cylinders. J. Non-Newtonian Fluid Mech., 22, 129–167.

    Google Scholar 

  • Bisgaard, C. (1983) Velocity fields around spheres and bubbles investigated by laser-Doppler anemometry, J Non-Newtonian Fluid Mech., 12, 282–302.

    Google Scholar 

  • Brezzi, F., and Fortin, M. (1991) Mixed and Hybrid Finite Element Methods. Series in Computational Mechanics, Springer-Verlag, New York.

    Google Scholar 

  • Brooks, A.N., and Hughes, T.J.R. (1982) Streamline upwind/Petrov Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier—Stokes equations. Comput. Methods Appl. Mech. Engrg., 32, 199–259.

    Google Scholar 

  • Brown, R.A., Armstrong, R.C. Beris, A.N., and Yeh, P.-W. (1986) Finite element analysis of complex viscoelastic flows. Comput. Methods Appl. Mech. Engrg., 58, 201–226.

    Google Scholar 

  • Burdette, S., Coates, P.J., Armstrong, R.C., and Brown, R.A. (1989) Finite element calculation of viscoelastic flow through an axisymmetric corrugated tube using the explicitly elliptic momentum equation formulation (EEME). J. Non-Newtonian Fluid Mech., 33, 1–23.

    Google Scholar 

  • Carey, G.F., and Oden, J.T. (1986) Finite Elements, Vols. 4 and 6. Prentice-Hall, Englewood-Cliffs, NJ.

    Google Scholar 

  • Chilcott, M.D., and Rallison, J.M. (1988) Creeping flow of dilute polymer solutions past cylinders and spheres. J. Non-Newtonian Fluid Mech., 29, 381–432.

    Google Scholar 

  • Coates, P.J., Armstrong, R.C., and Brown, R.A. (1992) Calculation of steady-state viscoleastic flow through axisymmetric contractions using the EEME formulation. J. Non-Newtonian Fluid Mech., 42, 141–188.

    Google Scholar 

  • Crochet, M.J., Davies, A.R., and Walters, K. (1984) Numerical Simulation of Non-Newtonian Flow. Elsevier, Amsterdam.

    Google Scholar 

  • Dupret, F., Marchal, J.M., and Crochet, M.J. (1985) On the consequence of discretization errors in the numerical calculation of viscoelastic flow. J. Non-Newtonian Fluid Mech., 18, 173–186.

    Google Scholar 

  • Franca, L.P., Frey, S.L., and Hughes, T.J.R. (1992) Stabilized finite element methods: I. Application to the advective—diffusive model. Comput. Methods Appl. Mech. Engrg., 95, 253–276.

    Google Scholar 

  • Gorodstov, V.A., and Leonov, A.I. (1967) On a linear instability of plane parallel Couette flow of viscoelastic fluid. J. Appl. Math. Mech., 31, 289–299.

    Google Scholar 

  • Hassager, O., and Bisgaard, C. (1983) A Lagrangian finite element method for the simulation of flow of non-Newtonian liquids. J. Non-Newtonian Fluid Mech., 12, 153–164.

    Google Scholar 

  • Hughes, T.J.R., Franca, L.P., and Hulbert, G.M. (1989) A new finite element formulation for computational fluid dynamics. VIII. The Galerkin least-squares method for advective—diffusive equations. Comput. Methods Appl. Mech. Engrg., 73, 173–189.

    Google Scholar 

  • Johnson, C., Nävert, U., and Pitkaranta, J. (1984) Finite element methods for linear hyperbolic problems. Comput. Methods Appl. Mech. Engrg., 45, 285–312.

    Google Scholar 

  • Joseph, D.D., Renardy, M., and Saut, J.C. (1985) Hyperbolicity and change of type in the flow of viscoelastic fluids. Arch. Rational Mech. Anal., 87, 213–251.

    Google Scholar 

  • Keiller, R.A. (1992) Numerical instability of time dependent flows. J. Non-Newtonian Fluid Mech., 43, 229–246.

    Google Scholar 

  • King, R.K., Apelian, M.R., Armstrong, R.C., and Brown, R.A. (1988) Numerically stable finite-element methods for solution of steady viscoelastic flows. J. Non-Newtonian Fluid Mech., 29, 147–216

    Google Scholar 

  • Leborgne, G. (1992) Simulation numerique d'ecoulements de fluides viscoelasticques de type Oldroyd: utilisation d'element finis de degre 1. Doctoral thesis, I'Ecole Poltechnique, Paris.

    Google Scholar 

  • Lee, K., and Finlayson, B.A. (1986) Stability of plane Poiseuille flow and Couette flow of a Maxwell fluid. J. Non-Newtonian Fluid Mech., 21, 65–78.

    Google Scholar 

  • Lunsmann, W.J., Genieser, L., Armstrong, R.C., and Brown, R.A. (1993) Finite element analysis of steady viscoelastic flow around a sphere in a tube: calculations with constant viscosity models. J. Non-Newtonian Fluid Mech., in press.

  • Marchal, J.M., and Crochet, M.J. (1987) A new mixed finite element method for calculating viscoelastic flow. J. Non-Newtonian Fluid Mech., 26, 77–114.

    Google Scholar 

  • Mendelson, M.A., Yeh, P.-W., Brown, R.A., and Armstrong, R.C. (1982) Approximation error in finite element calculation of viscoelastic fluid flows. J. Non-Newtonian Fluid Mech., 10, 31–54.

    Google Scholar 

  • Northey, P., Armstrong, R.C., and Brown, R.A. (1990) Finite element calculations of time-dependent two-dimensional flows of an upper-convected Maxwell fluid using the EEME formulation. J. Non-Newtonian Fluid Mech., 36, 109–134.

    Google Scholar 

  • Northey, P., Armstrong, R.C., and Brown R.A. (1992) Finite amplitude time-periodic states in viscoelastic Taylor—Couette flow of a UCM fluid. J. Non-Newtonian Fluid Mech., 42, 117–140.

    Google Scholar 

  • Pilitsis, S., and Beris, A.N. (1989) Calculations of steady-state viscoelastic flow in an undulating tube. J. Non-Newtonian Fluid Mech., 31, 231–287.

    Google Scholar 

  • Rajagopalan, D., Armstrong, R.C., and Brown, R.A. (1990a) Finite element analysis of viscoelastic flow of a multi-mode Maxwell fluid: application of EEME formulation. J. Non-Newtonian Fluid Mech., 36, 135–158.

    Google Scholar 

  • Rajagopalan, D., Armstrong, R.C., and Brown, R.A. (1990b) Finite element methods for calculation of viscoelastic flow for constitutive equations with a Newtonian viscosity. J. Non-Newtonian Fluid Mech., 36, 159–192.

    Google Scholar 

  • Renardy, M. (1985) Existence of slow steady flows of viscoelastic fluids with differential constitutive equations. Z. Angew. Math. Mech., 65, 449–451.

    Google Scholar 

  • Renardy, M., and Renardy, Y. (1986) Linear stability of plane Couette flow of an upper convected Maxwell fluid. J. Non-Newtonian Fluid Mech., 22, 23–33.

    Google Scholar 

  • Renardy, M., Hrusa, W.J., and Nohel, J.A. (1987) Mathematical Problems in Viscoelasticity. Longmans, New York.

    Google Scholar 

  • Rutkevich, I.M. (1970) The propogation of small perturbations in a viscoelastic fluid. J. Appl. Math. Mech., 34, 35–50.

    Google Scholar 

  • Rutkevich, I.M. (1972) On the thermodynamic interpretation of the evolutionary conditions of the equations of the mechanics of finitely deformable viscoelastic media of Maxwell type. J. Appl. Math. Mech., 36, 283–295.

    Google Scholar 

  • Zheng, R., Phan-Thien, N., and Tanner, R.I. (1990) On the flow past a sphere in a cylinderical tube: limiting Weissenberg number. J. Non-Newtonian Fluid Mech., 36, 27–49.

    Google Scholar 

Download references

Author information

Authors and Affiliations

  1. Department of Chemical Engineering, Massachusetts Institute of Technology, 02139, Cambridge, MA, U.S.A.

    R. A. Brown, M. J. Szady, P. J. Northey & R. C. Armstrong

Authors
  1. R. A. Brown
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. M. J. Szady
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. P. J. Northey
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. R. C. Armstrong
    View author publications

    You can also search for this author in PubMed Google Scholar

Additional information

Communicated by D.D. Joseph and M. Renardy

Financial support for this research was given by the National Science Foundation, the Office of Naval Research, and the Defense Research Projects Agency. Computational resources were supplied by a grant from the Pittsburgh National Supercomputer Center and by the MIT Supercomputer Facility.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Brown, R.A., Szady, M.J., Northey, P.J. et al. On the numerical stability of mixed finite-element methods for viscoelastic flows governed by differential constitutive equations. Theoret. Comput. Fluid Dynamics 5, 77–106 (1993). https://doi.org/10.1007/BF00311812

Download citation

  • Issue Date:

  • DOI: https://doi.org/10.1007/BF00311812

Keywords

Search

Navigation