Simulation of nonisothermal and time-dependent viscoelastic flows
Name
47351879-MIT.pdf
Description
Full printable version
Size
19.02 MB
Format
Adobe PDF
Checksum (MD5)
e86b9b7dc6373ba7d9b85a6bf2da8e80
Author(s)
Smith, Mark Davis, 1969-
Advisor(s)
Robert C. Armstrong and Robert A. Brown.
Date Issued
2000
Publisher
Massachusetts Institute of Technology
Abstract
The need to predict the interactions between the microstructure of polymeric fluids and the macroscopic flow field in polymer processing operations has lead to the development of many numerical methods for the simulation of viscoelastic flows. The two main obstacles to the use of these tools in the quantitative analysis of polymer processing operations are the need for accurate constitutive equations and the computational efficiency of the numerical methods. The goal of the first part of this thesis is to develop constitutive equations appropriate for modeling nonisothermal viscoelastic flows and to examine the influence of temperature variations on the stress in the fluid. The pseudotime method developed by Crochet and Naghdi (1978) and by Sarti (1977) is used to obtain thermodynamically consistent constitutive equations for the internal energy, the stress, and the heat flux vector. Fourier"s Jaw is assumed for the thermal conductivity, and the internal energy is chosen to depend only on the temperature of the fluid. The non isothermal analog to the Giesekus constitutive equation for the polymer contribution to the stress is then de1ived with the pseudotime approach. The constitutive equations are combined with conservation equations for mass. momentum, and energy to model the nonisothermal viscoelastic extrudate swell problem. For contraction-shaped dies, the final thickness of the extrudate is shown to depend on both the thermal and the strain histories of the fluid emerging, at the die exit. In this thesis second part of this thesis, efficient time integration algorithms are developed within the framework of examining the linear stability of isothermal viscoelastic two dimensional steady flows in complex geometries to three-dimensional perturbations. In this approach, a steady base state is found by discretizing in two dimensions with the discrete elastic-viscous split stress gradient (DEVSS-G) finite element formulation of the momentum-continuity pair combined with either the streamwise upwind Petrov-Galerkin (SUPG) or discontinuous Galerkin (DG) finite element discretization of the hyperbolic constitutive equation for the polymer contribution to the stress tensor. The resulting set of nonlinear algebraic equations is solved by using the Newton-Raphson method. The stability of these states is determinined by time integration of the evolution equations for infinitessimal perturbations governed by the equations of motion and the constitutive equation for the stress linearized about the base state. Two time integration methods are investigated: the {}-method operator splitting scheme and a fourth-order Runge-Kutta method. Both schemes reduce to a solution to 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 parts 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 time integration alone. it does provide a more complete analysis of the eigenspectrum. The linear stability of the flow of an Oldroyd-B model is analyzed with both the time integration methods and the hybrid time integration/Krylov algorithm for flow through several model geometries. First, the circular Couette problem is examined to benchmark the methods, and the stability results presented here are in good agreement with those obtained by other methods of analysis. The second model geometry is flow around an isolated cylinder confined in a channel. This investigation is motivated by the observation that the flow of a Boger fluid in this geometry undergoes a transition from a two-dimensional flow to a three-dimensional flow at a critical value of the Weissenberg number (Mc Kinley, 1993 ). Contrary to experimental observations, the flow with the Oldroyd-B model is predicted to be stable for all values of the Weissenberg number investigated. We
Description
Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 2000.
Includes bibliographical references (leaves 275-280).
Subjects
Chemical Engineering.
MIT Department
Massachusetts Institute of Technology. Department of Chemical Engineering
Terms of Use
M.I.T. theses are protected by copyright. They may be viewed from this source for any purpose, but reproduction or distribution in any format is prohibited without written permission. See provided URL for inquiries about permission.
Persistent DSpace Link
