Numerical investigation of the thermally and flow induced crystallization behavior of semi-crystalline polymers by using finite element-finite difference method

Abstract The thermally and flow induced crystallization behavior of semi-crystalline polymer in processing can significantly influence the quality of final products. The investigation of its mechanism has both scientific and industrial interest. A mathematical model in three dimensions for thermally and flow induced crystallization of polymer melts obeying differential Phan-Thien and Tanner (PTT) constitutive model has been developed and solved by using the finite element–finite difference method. A penalty method is introduced to solve the nonlinear governing equations with a decoupled algorithm. The corresponding finite element–finite difference model is derived by using the discrete elastic viscous split stress algorithm incorporating the streamline upwind scheme. A modified Schneider's approach is employed to discriminate the relative roles of the thermal state and the flow state on the crystallization phenomenon. The thermally and flow induced crystallization characteristics of polypropylene is investigated based on the proposed mathematical model and numerical scheme. The half crystallization time of polypropylene in a cooled couette flow configuration obtained by simulation are compared with Koscher's experimental results, which show that they agree well with each other. Two reasons to cause crystallization of polypropylene in pipe extrusion process including the thermal state and the flow state are investigated. Both the crystalline distribution and crystalline size of polypropylene are obtained by using the finite element–finite difference simulation of three-dimensional thermally and flow induced crystallization. The effects of processing conditions including the volume flow rate and the temperature boundary on the crystallization kinetics process are further discussed.

[1]  Hongbin Zhang,et al.  Experimental study of flow-induced crystallization in the blends of isotactic polypropylene and poly(ethylene-co-octene) , 2008 .

[2]  Antonios K. Doufas,et al.  A continuum model for flow-induced crystallization of polymer melts , 1999 .

[3]  C. Y. Wang,et al.  Prediction of isothermal crystallization parameters in monomer cast nylon 6 , 2008, Comput. Chem. Eng..

[4]  Three‐dimensional simulation of planar contraction viscoelastic flow by penalty finite element method , 2009 .

[5]  Jie Ouyang,et al.  Computer modeling of isothermal crystallization in short fiber reinforced composites , 2011, Comput. Chem. Eng..

[6]  G. G. Peters,et al.  Development and Validation of a Recoverable Strain‐Based Model for Flow‐Induced Crystallization of Polymers , 2001 .

[7]  N. Hudson,et al.  Polymer complexation effects in extensional flows , 1988 .

[8]  V. Brucato,et al.  On the Simulation of Thermoplastic Injection Moulding Process , 1997 .

[9]  R. Tanner,et al.  A new constitutive equation derived from network theory , 1977 .

[10]  Clara Silvestre,et al.  Non-isothermal crystallization of polymers , 1999 .

[11]  R. Tanner,et al.  A comparison of some models for describing polymer crystallization at low deformation rates , 2005 .

[12]  R. Tanner.  On the flow of crystallizing polymers: I. Linear regime , 2003 .

[13]  Caliane Bastos Borba Costa,et al.  Considerations on the crystallization modeling: Population balance solution , 2007, Comput. Chem. Eng..

[14]  Anthony J. McHugh,et al.  A continuum model for the dynamics of flow‐induced crystallization , 1996 .

[15]  R. Tanner.  A suspension model for low shear rate polymer solidification , 2002 .

[16]  S. Coppola,et al.  Effects of molecular weight distribution on the flow-enhanced crystallization of poly(1-butene) , 2008 .

[17]  R. Zheng,et al.  A model for post-flow induced crystallization: General equations and predictions , 2004 .

[18]  Hrvoje Jasak,et al.  Viscoelastic fluid analysis in internal and in free surface flows using the software OpenFOAM , 2010, Comput. Chem. Eng..

[19]  L. Huiping,et al.  Research on key technologies of FEM simulation of temperature field in the process of quenching , 2007 .

[20]  R. Fulchiron,et al.  Influence of shear on polypropylene crystallization: morphology development and kinetics , 2002 .

[21]  A. Doufas.  Analysis of the rheotens experiment with viscoelastic constitutive equations for probing extensional rheology of polymer melts , 2006 .

[22]  M. Fortin,et al.  A new mixed finite element method for computing viscoelastic flows , 1995 .

[23]  H. Janeschitz-Kriegl,et al.  New extensional rheometer for creep flow at high tensile stress. Part II. Flow induced nucleation for the crystallization of iPP , 2004 .

[24]  P. Chantrenne,et al.  Numerical modeling of nonisothermal polymer crystallization kinetics: Flow and thermal effects , 2010 .

[25]  K. Funatsu,et al.  Numerical simulation of three‐dimensional viscoelastic flow within dies , 1999 .

[26]  W. Schneider,et al.  Non – Isothermal Crystallization Crystallization of Polymers , 1988 .

[27]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[28]  M. Boutaous,et al.  Thermally and flow induced crystallization of polymers at low shear rate , 2010 .

[29]  J. Hoffman,et al.  Kinetic of crystallization from the melt and chain folding in polyethylene fractions revisited: theory and experiment , 1997 .

[30]  M. Crochet,et al.  A new mixed finite element for calculating viscoelastic flow , 1987 .

[31]  H. Janeschitz-Kriegl,et al.  The Role of Long Molecules and Nucleating Agents in Shear Induced Crystallization of Isotactic Polypropylenes** , 1997 .

[32]  Roland Keunings,et al.  Parallel finite element algorithms applied to computational rheology , 1995 .