2026-09-18T22:30:43Zhttps://dspace.mit.edu/server/oai/request

oai:dspace.mit.edu:1721.1/292932022-01-13T07:54:19Zcom_1721.1_7582com_1721.1_7581col_1721.1_131022

Robert C. Armstrong. Suen, Jason Ka-Chun, 1973- Massachusetts Institute of Technology. Dept. of Chemical Engineering. Massachusetts Institute of Technology. Department of Chemical Engineering 2005-10-14T19:41:29Z 2005-10-14T19:41:29Z 2002 2003 http://hdl.handle.net/1721.1/29293 52297506 Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Chemical Engineering, February 2003. Includes bibliographical references (v. 2, leaves 351-363). Knowing and understanding the dynamics and molecular configurations of polymer molecules is important for efficient process design and novel product development. Much research has been focused on combining molecular simulations and traditional fluid mechanics computations to simulate the behavior of polymeric liquids in a fabrication process. These simulation approaches require solution of the coupled set of the equations of change, the governing equations from kinetic theory, and the flux expressions that map molecular configurations to macroscopic quantities. Most complex flow simulations so far make use of a mixed finite element method to calculate the velocity field, with stress tensor evaluated by using a stochastic simulation method. This so-called CONNFFESSIT approach suffers from both a large memory requirement and stochastic noise. This thesis focuses on the development and application of a fully deterministic numerical approach for computing viscoelastic flows with constitutive descriptions based directly on diffusion equations from kinetic theory. The numerical approach is based on an operator splitting time integration method that decouples the calculation of microstructure by solution of a hyperbolic diffusion equation from the velocity and pressure field evolution, which is obtained by solution of a generalized Stokes problem. The generalized Stokes problem is written in the DEVSS-G formulation, where a direct interpolation of the components of the velocity gradient tensor is introduced. The efficiency and robustness of this numerical method is demonstrated through calculating the viscoelastic flows of a modified Doi model for liquid crystalline polymer and a number of reptation models for polymer melts in different flow geometries. (cont.) Simulations of the original Doi model with the Maier-Saupe potential in a pressure-driven channel flow by Nayak showed disclination formation due to the shear-rate-dependent frequencies of the tumbling/wagging states of the Doi model in a shear or mixed shear flow. The lack of an instrinsic length scale in the model leads to an infinitesimal structure refinement that eventually causes numerical instabilities. In this thesis, the effect of concentration variation is incorporated to develop a modified Doi model for introducing an intrinsic length scale through translational diffusion. This changes the mathematical characteristics of the spatial variation of the underlying diffusion equation from that of a hyperbolic equation to that of an elliptic equation. The resulting elliptic diffusion equation is then solved by using a local discontinuous Galerkin method, where an auxillary variable is introduced to rewrite the elliptic diffusion equation into a pair of formal, hyperbolic equations, which in turn is solved by the standard discontinuous Galerkin method. Unlike the original Doi model, a steady state is reached for a variety of De. Although there is structure variation across the channel width, the director profiles point uniformly along the flow direction. The lack of disclination formation may be rectified by introducing Frank elasticity into the modified Doi model ... by Jason Ka-Chun Suen. Ph.D. 2 v. (363 leaves) 13887149 bytes 13886950 bytes application/pdf application/pdf eng Massachusetts Institute of Technology 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 Chemical Engineering. Multiscale simulation of viscoelastic flows : applications to kinetic theory models of polymer melts and liquid crystalline polymers Thesis application/pdf false Publication unknown unknown <Publication xmlns="https://www.openaire.eu/cerif-profile/1.1/" id="1fa88568-2eac-4995-8d90-4b20c6c98775"> <Type xmlns="https://www.openaire.eu/cerif-profile/vocab/COAR_Publication_Types">http://purl.org/coar/resource_type/c_1843</Type> <Language>eng</Language> <Title>Multiscale simulation of viscoelastic flows : applications to kinetic theory models of polymer melts and liquid crystalline polymers</Title> <PublishedIn> <Publication> </Publication> </PublishedIn> <PublicationDate>2003</PublicationDate> <Authors> <Author> <DisplayName>Suen, Jason Ka-Chun, 1973-</DisplayName> <Affiliation> <OrgUnit> </OrgUnit> </Affiliation> </Author> </Authors> <Editors> </Editors> <Publishers> <Publisher> <DisplayName>Massachusetts Institute of Technology</DisplayName> <OrgUnit /> </Publisher> </Publishers> <License>http://dspace.mit.edu/handle/1721.1/7582</License> <Keyword>Chemical Engineering.</Keyword> <Abstract>Knowing and understanding the dynamics and molecular configurations of polymer molecules is important for efficient process design and novel product development. Much research has been focused on combining molecular simulations and traditional fluid mechanics computations to simulate the behavior of polymeric liquids in a fabrication process. These simulation approaches require solution of the coupled set of the equations of change, the governing equations from kinetic theory, and the flux expressions that map molecular configurations to macroscopic quantities. Most complex flow simulations so far make use of a mixed finite element method to calculate the velocity field, with stress tensor evaluated by using a stochastic simulation method. This so-called CONNFFESSIT approach suffers from both a large memory requirement and stochastic noise. This thesis focuses on the development and application of a fully deterministic numerical approach for computing viscoelastic flows with constitutive descriptions based directly on diffusion equations from kinetic theory. The numerical approach is based on an operator splitting time integration method that decouples the calculation of microstructure by solution of a hyperbolic diffusion equation from the velocity and pressure field evolution, which is obtained by solution of a generalized Stokes problem. The generalized Stokes problem is written in the DEVSS-G formulation, where a direct interpolation of the components of the velocity gradient tensor is introduced. The efficiency and robustness of this numerical method is demonstrated through calculating the viscoelastic flows of a modified Doi model for liquid crystalline polymer and a number of reptation models for polymer melts in different flow geometries.</Abstract> <Abstract>(cont.) Simulations of the original Doi model with the Maier-Saupe potential in a pressure-driven channel flow by Nayak showed disclination formation due to the shear-rate-dependent frequencies of the tumbling/wagging states of the Doi model in a shear or mixed shear flow. The lack of an instrinsic length scale in the model leads to an infinitesimal structure refinement that eventually causes numerical instabilities. In this thesis, the effect of concentration variation is incorporated to develop a modified Doi model for introducing an intrinsic length scale through translational diffusion. This changes the mathematical characteristics of the spatial variation of the underlying diffusion equation from that of a hyperbolic equation to that of an elliptic equation. The resulting elliptic diffusion equation is then solved by using a local discontinuous Galerkin method, where an auxillary variable is introduced to rewrite the elliptic diffusion equation into a pair of formal, hyperbolic equations, which in turn is solved by the standard discontinuous Galerkin method. Unlike the original Doi model, a steady state is reached for a variety of De. Although there is structure variation across the channel width, the director profiles point uniformly along the flow direction. The lack of disclination formation may be rectified by introducing Frank elasticity into the modified Doi model ...</Abstract> <Access xmlns="http://purl.org/coar/access_right" > </Access> </Publication>