聚合物熔体三维等温流动的罚有限元分析

Penalty Finite Element Analysis of Three Dimensional Isothermal Polymer Melt Flow

本文针对幂律流体和非线性粘弹性PTT流体,采用罚有限元法,形成了求解聚合物熔体三维等温流动速度场和粘弹性应力场的有限元模型。对非线性粘弹性PTT流体,为了降低模拟计算对计算机硬件的要求,并使模拟计算更加稳定,采用了去耦算法,包括拟体力方法和动量方程的椭圆类方程转化方法等。文中还给出了总体有限元方程组的形成和求解过程。

In this paper, the penalty finite element method is adopted to formulate the finitr element equations for computation of the velocity and viscoelastic stress distribution s three dimensional isothermal polymer melt flow under the constitutive relations power-law and non-linear viscoelastic Phan Thien-Tanner(PTT). For which a couple of decoupled numerical methods, such as pseudo-body force and transformation of momentum equation into the elliptic equation, are used to decrease the capacity of computer memory required for numerical simulation and raise the convergence stability. Finally, the formation and solution of the general element equations were stated briefly.