化工学报 ›› 2016, Vol. 67 ›› Issue (5): 2144-2151.DOI: 10.11949/j.issn.0438-1157.20151472

• 材料化学工程与纳米技术 • 上一篇    下一篇

剪切流场中聚乙烯结晶过程的建模与模拟

阮春蕾1,2, 刘春太2   

  1. 1 河南科技大学数学与统计学院, 河南 洛阳 471023;
    2 橡塑模具国家工程研究中心, 河南 郑州 450001
  • 收稿日期:2015-09-21 修回日期:2016-01-11 出版日期:2016-05-05 发布日期:2016-05-05
  • 通讯作者: 阮春蕾
  • 基金资助:

    国家自然科学基金项目(11402078,51375148)。

Numerical simulation of morphology and kinetics of polyethylene in shear flow

RUAN Chunlei1,2, LIU Chuntai2   

  1. 1 School of Mathematics and Statistics, Henan University of Science and Technology, Luoyang 471023, Henan, China;
    2 National Engineering Research Center for Advanced Polymer Processing Technology, Zhengzhou 450001, Henan, China
  • Received:2015-09-21 Revised:2016-01-11 Online:2016-05-05 Published:2016-05-05
  • Supported by:

    supported by the National Natural Science Foundation in China (11402078, 51375148).

摘要:

基于Eder模型推导了剪切流场中球晶、串晶形态演化的数学模型,将第一法向应力差作为串晶成核的驱动,并引入两相悬浮模型描述体系,认为其由无定形相和半结晶相组成,分别用FENE-P模型和刚性哑铃模型描述。基于上述数学模型,分别构造了捕捉球晶、串晶生长的Monte Carlo法与体系控制方程求解的有限差分法,成功模拟了二维剪切流场中聚乙烯的结晶过程,给出了球晶、串晶的形态演化,分析了剪切时间、剪切速率对串晶成核密度、结晶速率、流体黏度等的影响。数值结果表明:所构造的Monte Carlo法合理有效,不仅成功捕捉了晶体的生长与碰撞,而且较为准确地预测了结晶速率。此外,提高剪切时间或剪切速率,将增加串晶成核密度、提高结晶速率、使流体黏度突增的时间点提前。

关键词: 流动诱导结晶, 形态学, 动力学, Monte Carlo模拟, 结晶

Abstract:

The mathematical model of morphology evolution of spherulites and shish-kebabs in the shear flow is deduced based on the Eder model. The model considers that the nucleation and growth of spherulites are determined by the static temperature, while the nucleation and growth of shish-kebabs are determined by the flow, which is depended on the first normal stress difference of the system and the shear rate respectively. In order to calculate the nucleation density of the shish-kebabs, the two-phase suspension model of Zheng is introduced. The model treats the stress as the combination of the amorphous phase and semi-crystalline phase. The amorphous phase is described by FENE-P model while the semi-crystalline phase is depicted by a rigid dumbbell model. Based on the mathematical model, the Monte Carlo method and the finite difference method are constructed, respectively. The former is to capture the crystal growth while the latter is to calculate the equation of the system. By using these methods, the 2D crystallization in the shear flow is simulated. The evolution of spherulites and shish-kebabs is given. Also, the effects of shear time and shear rate on the nucleation density of shish-kebabs, crystallization rate, the viscosity of the fluid and the system stress are discussed. Numerical results show that the Monte Carlo method is valid which not only captures the morphology evolution of crystals successfully, but also predicts the crystallization rate well. In addition, the increase of the shear time or the shear rate will increase the nucleation density of shish-kebabs and the crystallization rate.

Key words: flow induced crystallization, morphology, kinetics, Monte Carlo simulation, crystallization

中图分类号: