The solution of large viscoelastic flow problems using parallel iterative techniques
Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 2000.
| Main Author: | |
|---|---|
| Other Authors: | |
| Format: | Thesis |
| Language: | eng |
| Published: |
Massachusetts Institute of Technology
2005
|
| Subjects: | |
| Online Access: | http://hdl.handle.net/1721.1/8827 |
| _version_ | 1826193307545370624 |
|---|---|
| author | Caola, Anthony E |
| author2 | Robert A. Brown and Robert C. Armstrong. |
| author_facet | Robert A. Brown and Robert C. Armstrong. Caola, Anthony E |
| author_sort | Caola, Anthony E |
| collection | MIT |
| description | Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 2000. |
| first_indexed | 2024-09-23T09:36:55Z |
| format | Thesis |
| id | mit-1721.1/8827 |
| institution | Massachusetts Institute of Technology |
| language | eng |
| last_indexed | 2024-09-23T09:36:55Z |
| publishDate | 2005 |
| publisher | Massachusetts Institute of Technology |
| record_format | dspace |
| spelling | mit-1721.1/88272019-04-12T13:46:58Z The solution of large viscoelastic flow problems using parallel iterative techniques Caola, Anthony E Robert A. Brown and Robert C. Armstrong. Massachusetts Institute of Technology. Dept. of Chemical Engineering. Massachusetts Institute of Technology. Dept. of Chemical Engineering. Chemical Engineering. Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, 2000. Includes bibliographical references (p. 306-325). Efficient parallel computation of complex viscoelastic flows is essential to bring modem computer power to bear on fluid calculations where complicated constitutive descriptions are required. We have developed an efficient, parallel time integration scheme for calculating viscoelastic flows. The method is a synthesis of an operator splitting and time integration method that decouples the calculation of the polymer portion of the stress by solution of a hyperbolic constitutive equation from the temporal updating of the velocity and pressure fields by solution of a generalized Stokes problem. Finite element discretizations using the DEVSS-G 1 method are used; here a direct interpolation of the components of the velocity gradient is introduced. The key to the parallelization is the incorporation of iterative matrix solution for each portion of the split calculation. The linear system that results from the generalized Stokes problem is asymmetric, indefinite, and block singular. At each timestep, this system is solved using an algorithm which combines a parallel preconditioner with a Krylov subspace method. The parallel preconditioner, called the Block Complement and Additive Levels Method (BCALM) preconditioner, is based on treating pressure unknowns separately from the variables velocities and gradients. A pressure preconditioner is constructed from a factor of the Schur complement of the pressures using a Jacobi iteration. The viscous operator is treated using the additive Schwarz method. The number of Krylov iterations required per implicit step with increasing numbers of processors is stabilized using a two-level additive Schwarz method. The coarse grid is generated using addition as a restriction operator and insertion for extension. The resulting iterative method is demonstrated to have high parallel efficiency, subject to effective domain decomposition. Robustness and good performance result from using geometric partitioning of the unknowns by the nested bisection routines supplied in CHAC02. The software is implemented using MPI and the PETSc toolkit, so as to be readily portable to a variety of parallel computers. This solver is developed in the context of the natural convection problem, where in the limit of high Grashof number, the linear system that arises during solution for steady state using Newton's method is stiff, assymmetric and indefinite. With this model problem, the preconditioner's robustness is tested and its efficiency is measured. The solver's generalizability for viscous flow problems is also discussed, as the only assumption about problem physics is that the continuity equation exists. For a given velocity field, the discretized equations from a typical differential constitutive model, e.g. the Oldroyd-B model, yields an asymmetric, nonsingular system of linear equations when nonlinearities are updated explicitly. Krylov iterations preconditioned with inexact solves in a two-level additive Schwarz framework are used to solve this system at each time step. Added efficiency is introduced by using a discontinuous Galerkin spatial discretization, which allows element-by-element condensation of the linear systems and higher parallel efficiency. The effectiveness of the algorithm is demonstrated by the calculation of transient, two dimensional flow of an Oldroyd-B fluid past an isolated cylinder confined symmetrically between parallel plates and by calculation of the stability of the steady-state motion to small-amplitude, three-dimensional disturbances The resulting method is far superior to the use of serial, frontal methods. by Anthony Edward Caola. Ph.D. 2005-08-23T15:39:14Z 2005-08-23T15:39:14Z 2000 2000 Thesis http://hdl.handle.net/1721.1/8827 48385756 eng 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. http://dspace.mit.edu/handle/1721.1/7582 325 p. 20421494 bytes 20421253 bytes application/pdf application/pdf application/pdf Massachusetts Institute of Technology |
| spellingShingle | Chemical Engineering. Caola, Anthony E The solution of large viscoelastic flow problems using parallel iterative techniques |
| title | The solution of large viscoelastic flow problems using parallel iterative techniques |
| title_full | The solution of large viscoelastic flow problems using parallel iterative techniques |
| title_fullStr | The solution of large viscoelastic flow problems using parallel iterative techniques |
| title_full_unstemmed | The solution of large viscoelastic flow problems using parallel iterative techniques |
| title_short | The solution of large viscoelastic flow problems using parallel iterative techniques |
| title_sort | solution of large viscoelastic flow problems using parallel iterative techniques |
| topic | Chemical Engineering. |
| url | http://hdl.handle.net/1721.1/8827 |
| work_keys_str_mv | AT caolaanthonye thesolutionoflargeviscoelasticflowproblemsusingparalleliterativetechniques AT caolaanthonye solutionoflargeviscoelasticflowproblemsusingparalleliterativetechniques |