Anbiao Chen

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 Schneiders 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 Koschers 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.

Numerical investigation of the crystallization and orientation behavior in polymer processing with a two-phase model

Abstract The crystallization and orientation behavior of a polymeric material can significantly influence the performance of products in practical processing. In this study, the variations in morphology that occur during solidification in polymer processing are mathematically modeled using a two-phase model. The amorphous phase is approximated as a finite extensible nonlinear elastic dumbbell with a Peterlin closure approximation (FENE-P) fluid, and the semi-crystalline phase is modeled as rigid rods oriented within the flow field. The crystallization and orientation behavior are numerically investigated using the penalty finite element–finite difference method with a decoupled algorithm. The evolution of the crystallization process is described by Schneiders equation, which differentiates between the effects of thermal and flow states. The hybrid closure approximation is adopted for the calculation of the three-dimensional orientation tensor. The discrete elastic viscous split stress (DEVSS) algorithm, which incorporates the streamline upwind scheme, is introduced to improve calculation stability. The variations in morphology during polymer processing are successfully predicted using the proposed mathematical model and numerical method. The influence of processing conditions on the crystallization and orientation behavior is further discussed.

Prediction for the mechanical property of short fiber-reinforced polymer composites through process modeling method

The fiber-reinforced polymer composites are important alternative for conventional structural materials because of their excellent comprehensive performance and weight reduction. The mechanical properties of such composite materials are mainly determined by the fiber orientation induced through practical manufacturing process. In the study, a through process modeling (TPM) method coupling the microstructure evolution and the mechanical properties of fiber-reinforced composites in practical processing is presented. The numerical methodology based on the finite volume method is performed to investigate three-dimensional forming process in the injection molding of fiber-reinforced composites. The evolution of fiber orientation distribution is successfully predicted by using a reduced strain closure model. The corresponding finite volume model for TPM is detailedly derived and the pressure implicit with splitting of operators (PISO) algorithm is employed to improve computational stability. The flow-induced multilayer structure is successfully predicted according to essential flow characteristics and the fiber orientation distribution. The mechanical properties of such anisotropy composites is further calculated based on the stiffness analysis and the Tandon–Weng model. The improvement of mechanical properties in each direction of the injection molded product are evaluated by using the established mathematical model and numerical algorithm. The influences of the geometric structure of injection mold cavity, the fiber volume fractions, and the fiber aspect ratios on the mechanical properties of composite products are further discussed. The mathematical model and numerical method proposed in the study can be successfully adopted to investigate the structural response of composites in practical manufacturing process that will be helpful for optimum processing design.

Numerical investigation of three-dimensional fiber suspension flow by using finite volume method

Fiber suspension flow is common in many industrial processes like papermaking and fiber-reinforcing polymer-based material forming. The investigation of the mechanism of fiber suspension flow is of significant importance, since the orientation distribution of fibers directly influences the mechanical and physical properties of the final products. A numerical methodology based on the finite volume method is presented in the study to simulate three-dimensional fiber suspension flow within complex flow field. The evolution of fiber orientation is described using different formulations including FT model and RSC model. The pressure implicit with splitting of operators algorithm is adopted to avoid oscillations in the calculation. A laminate structure of fiber orientation including the shell layer, the transition layer and the core layer along radial direction within a center-gated disk flow channel is predicted through a three-dimensional simulation, which agrees well with Mazahir’s experimental results. The evolution of fiber orientation during the filling process within the complex flow field is further discussed. The mathematical model and numerical method proposed in the study can be successfully adopted to predict fiber suspension flow patterns and hence to reveal the fiber orientation mechanism.

Finite element simulation of three-dimensional viscoelastic planar contraction flow with multi-mode FENE-P constitutive model

AbstractViscoelasticity is a characteristic of many complex fluids like polymer melts, petroleum, blood, etc. The investigation of viscoelastic flow mechanism has practical significance in both scientific and engineering field. Owing to strongly nonlinear, numerical method becomes a practical way to solve viscoelastic flow problem. In the study, the mathematical model of three-dimensional flow of viscoelastic fluids is established. The planar contraction flow as a benchmark problem for the numerical investigation of viscoelastic flow is solved by using the penalty finite element method with a decoupled algorithm. The multi-mode finitely extensible nonlinear elastic dumbbell with a Peterlin closure approximation (FENE-P) constitutive model is used to describe the viscoelastic rheological properties. The discrete elastic viscous split stress formulation in cooperating with the inconsistent streamline upwind scheme is employed to improve the computation stability. The numerical methods proposed in the study can be well used to predict complex flow patterns of viscoelastic fluids.