黏弹流体挤出胀大的数值模拟研究进展

主要介绍了黏弹流体挤出胀大的数值模拟研究进展.给出了黏弹流体挤出胀大的 数学模型,回顾了近20多年以来挤出胀大的主要数值模拟研究工作,然后对主要模拟 方法的计算过程、方法特点和形成的结果进行了一定的总结.最后提出了作者对挤 出胀大研究的一些看法,包括目前研究中存在的问题和相关研究的发展趋势.     

用PSM模型模拟聚合物熔体轴对称收缩流动

采用模拟黏弹流体挤出胀大的方法，计算了IUPAC-LDPE熔体经过4：1轴对称收缩流道的流动．计算的相对涡强度、入口校正和献中的结果基本一致，给出的流场也显示出计算结果是合理的．表明该方法能够适用于用积分型PSM模型表征的黏弹流体在收缩流道内的流动模拟．     

积分型Maxwell流体挤出胀大的数值模拟

应用Luo提出的基于常规有限单元的应力计算方法和范毓润回避奇点的方法,模拟了积分型Maxwell黏弹流体的挤出胀大流动,计算得到了Weissenberg数达1．0下的合理结果。     

基于格子Boltzmann方法的挤出胀大的数值模拟

将单相格子Boltzmann方法(lattice Boltzmann method,LBM)引入到粘弹流体的瞬态挤出胀大的数值模拟中,建立了基于双分布函数的自由面粘弹性流动格子Boltzmann模型.分析得到的流道中流动速度分布和构型张量结果与理论解十分吻合.对粘弹流体瞬态挤出胀大过程进行了模拟,并分析了运动粘度比和剪切速率对挤出胀大率的影响,得到的胀大率结果与理论分析和其它模拟结果基本一致.表明给出的LBM可以捕捉挤出胀大的瞬态效应.     

基于格子Boltzmann两相模型的粘弹流体挤出胀大分析

挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题．本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象，采用颜色模型模拟出口处粘弹流体和空气的两相流动，通过重新标色获得两种流体的界面，并最终获得胀大的形状．Navier-Stokes方程和本构方程的求解采用双分布函数模型．将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较，发现格子Boltzmann两相模型结果与解析解和实验结果相吻合，相比于单相模型，收敛速度更快，解的稳定性更高．研究了流道尺寸对胀大率的影响，并对挤出胀大的内在机理进行了分析．     

牛顿流体挤出胀大模拟的网格依赖性

应用有限元法模拟了牛顿流体的挤出胀大，通过方差分析得到胀大结果具有网格依赖性．基于四边形等参单元，得到出口后首排单元沿轴向的长度和近壁处首行单元沿径向的高度显著影响牛顿流体的挤出胀大．     

聚合物熔体三维挤出胀大的数值模拟

采用有限元方法分析K-BKZ本构方程描述的聚合物熔体的三维挤出胀大.对于本构方程中偏应力张量的计算,首先给出质点的运动轨迹,分段求出局部的变形梯度张量,再求出整体的变形梯度、Cauchy-Green应变张量和 Finger应变张量,沿轨迹采用分段高斯积分计算应力.把应力作为方程的右端项,给出迭代方法,求解非线性方程组.并根据自由面处的边界条件,迭代得出出口处自由面的最终位置.对轴对称流道和矩形流道进行分析计算,并把结果与二维分析和实验结果进行了比较,显示方法是可行的.     

黏弹流体流动的数值模拟研究进展

综述了黏弹流体流动数值模拟的研究进展,突出介绍近十年来有限元法在黏弹流体流动数值模拟研究中取得的成果,通过动量方程的适当变形和本构方程离散权函数的合理选择,可以显著增强数值计算的稳定性。得到较高Weissenberg数下的解,同时文中对黏弹流体流动数值模拟中本构方程的应用、非等温情况和三维空间下的研究进行了介绍。     

液晶高分子-各向异性流体挤出拉伸分岔研究

应用共转导数型本构方程研究了液晶高分子纺丝挤出过程的拉伸黏度,应用计算机符号运算软件 Maple得出解析表达式,拉伸黏度与拉伸率之间关系(随剪切速率变化)表明存在分岔现象,得出拉伸黏度显著高于相应的剪切黏度,解释了液晶高分子熔体挤出时不发生挤出胀大的物理机制．     

高体积含量非线性黏弹复合材料有效性质

提出了一个细观力学模型,可用于预测高体积含量非线性黏弹复合材料有效性质.该模型基于广义割线模量法、双球法以及Laplace-Carson变换技术.所提出的模型对玻璃微珠填充高密度聚乙烯(GB/HDPE)复合材料的应力应变关系进行了预测,结果与文献实验结果吻合;计算结果还表明在高体积百分比下文中所提出的方法比基于MT方法预测的粘性效应明显减弱;最后还将所提方法与线黏弹框架下的均质化模型做了比较,结果表明GB/HDPE表现出明显的非线性,线黏弹本构无法描述应变率对其力学行为的影响.     

Numerical simulation of delayed die swell

In a recent paper, Joseph et al. showed that, for a number of viscoelastic fluids, one can observe the phenomenon of delayed die swell beyond a critical extrusion velocity, or beyond a critical value of the viscoelastic Mach number. Giesekus had also observed that delayed die swell is a critical phenomenon.In the present paper, we find a set of material and flow parameters under which it is possible to simulate delayed die swell. For the viscoelastic flow calculation, we use the finite element algorithm with sub-elements for the stresses and streamline upwinding in the discretized constitutive equations. For the free surface, we use an implicit technique which allows us to implement Newton's method for solving the non-linear system of equations. The fluid is Oldroyd-B which, in the present problem, is a singular perturbation of the Maxwell fluid. The results show very little sensitivity to the size of the retardation time. We also show delayed die swell for a Giesekus fluid.This paper is dedicated to Professor Hanswalter Giesekus on the occasion of his retirement as Editor of Rheologica Acta.     

A boundary element investigation of extrudate swell

Numerical studies of die swell have until now dealt perimarily with the Maxwell or Oldroyd-B viscoelastic models. However, these models exhibit features that often make them unsuitable for numerical work. Furthermore, they are not realistic representations of actual viscoelastic fluids. In this report a comparison is made between the behaviour of a variety of different viscoelastic models when applied to the die swell problem. A wide range of elongational and shear behaviour is exhibited by the models examined. Both types of behaviour are shown to be important in the die swell problem, and the observed swelling is related to these characteristics of the models.     

Prediction of extrudate swell in polymer melt extrusion using an Arbitrary Lagrangian Eulerian (ALE) based finite element method

Accurate prediction of extrudate (die) swell in polymer melt extrusion is important as this helps in appropriate die design for profile extrusion applications. Extrudate swell prediction has shown significant difficulties due to two key reasons. The first is the appropriate representation of the constitutive behavior of the polymer melt. The second is regarding the simulation of the free surface, which requires special techniques in the traditionally used Eulerian framework. In this paper we propose a method for simulation of extrudate swell using an Arbitrary Lagrangian Eulerian (ALE) technique based finite element formulation. The ALE technique provides advantages of both Lagrangian and Eulerian frameworks by allowing the computational mesh to move in an arbitrary manner, independent of the material motion. In the present method, a fractional-step ALE technique is employed in which the Lagrangian phase of material motion and convection arising out of mesh motion are decoupled. In the first step, the relevant flow and constitutive equations are solved in Lagrangian framework. The simpler representation of polymer constitutive equations in a Lagrangian framework avoids the difficulties associated with convective terms thereby resulting in a robust numerical formulation besides allowing for natural evolution of the free surface with the flow. In the second step, mesh is moved in ALE mode and the associated convection of the variables due to relative motion of the mesh is performed using a Godunov type scheme. While the mesh is fixed in space in the die region, the nodal points of the mesh on the extrudate free surface are allowed to move normal to flow direction with special rules to facilitate the simulation of swell. A differential exponential Phan Thien Tanner (PTT) model is used to represent the constitutive behavior of the melt. Using this method we simulate extrudate swell in planar and axisymmetric extrusion with abrupt contraction ahead of the die exit. This geometry allows the extrudate to have significant memory for shorter die lengths and acts as a good test for swell predictions. We demonstrate that our predictions of extrudate swell match well with reported experimental and numerical simulations.     

Das unterschiedliche Schwellverhalten eines extrudierten Schmelzestranges bei Variation des Relaxationsspektrums

The die swell of viscoelastic fluids after leaving a die depends not only on the geometry, temperature effects and the mass flow rate, but also on the specific properties of the media. Among these specific properties is the relaxation spectrum which determines the relaxation of stresses in the fluid. It is this relaxation that produces a die swell. This paper investigates theoretically the manner in which modifications of the relaxation spectrum affect the viscosity curve and the die swell behaviour. In particular it is shown which modifications of the spectrum have a large influence on the die swell and which are relatively unimportant. It is also shown how the time intervals in which the main die swell occurs dependend on the shear rates in the die and which parts of the relaxation spectrum are responsible for the size of the die swell.     

半解析有限元法分析非线性粘弹体在板材模具中的流动

聚合物熔体在模具中的均匀性决定板材成型的质量。熔体具有非线性记忆性的粘弹效应,同时流道十分复杂。为准确、简便地分析熔体流动情况,针对板材机头流道的特点,提出了流动计算新的方法——有限柱半解析法,即在流道的厚度方向采用富氏级数描述流动分布,把三维有限元方程简化成为二维,得到以流量和压力为未知量的有限元方程。同时采用K-BKZ积分型本构模型描述熔体的非线性记忆性,给出求解非线性方程组的迭代方法,取得较好的收敛性。对两种典型的板材模具进行了分析计算,得出模具出口处的流量分布,通过把结果与三维有限元分析的结果相比较,表明该方法可以准确地分析板材模具中熔体的流动。     

Polypropylene-polyethylene blends

Summary Die swell behaviour and morphology of melt blends of isotactic polypropylene (PP) and high density polyethylene for pure polymers and blends with 25, 50 and 75 weight % PP are described in the present study. A light interference contrast microscopy technique was used for the morphological characterization of melt blends and extrudate samples of the blends obtained with an Instron capillary rheometer. The results indicate that the domains from blends where the dispersed phase has higher viscosity than the continuous phase remain as continuous domains in the extrudate whereas domain destruction takes place when blends where the continuous phase has the higher viscosity are extruded.The die swell behaviour as well as the fiber forming properties of extrudates of melts having unstable domains extruded at high shear stresses resemble the behaviour of homopolymers, whereas samples with stable domains are significantly different, die swell increases with temperature at constant shear stress and stable fibers cannot be obtained after necking.With 10 figures and 1 table     

Quantitative evaluation of extrudate swell from viscoelastic properties of polystyrene

A theory of extrudate swell for short, intermediate or long dies is presented. In our experiment, we consider that the swelling phenomenon is mainly due to the recoverable elongational strain induced by the converging flow at the die entrance, as well as by recoverable shear strain originating within the die. From these concepts, an equation has been derived for the quantitative prediction of extrudate swell from the elastic material properties such as the entrance pressure drop, the relaxation modulus and the recoverable shear strain. Excellent agreement is found between predicted and measured values of extrudate swell obtained on commercial polystyrene melt, using capillaries of length-to-diameter ratios ranging from 1 to 20 and in a wide range of shear rates.     

Viscoelastic simulations of stick‐slip and die‐swell flows

Numerical solutions of viscoelastic flows are demonstrated for a time marching, semi‐implicit Taylor–Galerkin/pressure‐correction algorithm. Steady solutions are sought for free boundary problems involving combinations of die‐swell and stick‐slip conditions. Flows with and without drag flow are investigated comparatively, so that the influence of the additional component of the drag flow may be analysed effectively. The influence of die‐swell is considered that has application to various industrial processes, such as wire coating. Solutions for two‐dimensional axisymmetric flows with an Oldroyd‐B model are presented that compare favourably with the literature. The study advances our prior fixed domain formulation with this algorithm, into the realm of free‐surface viscoelastic flows. The work involves streamline‐upwind/Petrov–Galerkin weighting and velocity gradient recovery techniques that are applied upon the constitutive equation. Free surface solution reprojection and a new pressure‐drop/mass balance scheme are proposed. Copyright © 2001 John Wiley & Sons, Ltd.     

A mixed finite element scheme for viscoelastic flows with XPP model

A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.