同心旋转圆柱间弱弹性流的非线性稳定性分析

基于同心旋转圆柱间Ｏｌｄｒｏｙｄ－Ｂ型流体的六维动力系统,探讨了小间隙大扰动条件下高分子添加剂对滑动轴承间油膜非线性稳定性的影响。结果表明,弱弹性流体的失稳结构与牛顿流体相似,随着转速的增加,流体以同宿轨道分岔失稳,与牛顿流体相比,少量的高分子添加剂具有推迟流体层流的稳的作用。     

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

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

单晶铁电体的宏细观本构关系描述

以高技术电子材料－铁电单晶体为研究对象，在对铁电体中电畸内的细观力学耦合场分析的基础上，采用Ｍｏｒｉ－Ｔａｎａｋａ方法以自洽的方式导出了材料构元的Ｈｅｌｍｈｏｌｔｚ自由能及Ｇｉｂｂｓ自由能函数的解析表达式。并分析在广义应力空间和广义应变空间中，按Ｈｉｌｌ－Ｒｉｃｅ内变量本构理论框架，导出了铁电体畸变屈服面方程，增量型本构关系及内变量的演化方程。文末给出了对ＢａＴｉＯ３单晶材料力电行为的一维数值模拟     

Sn—Pb共晶钎料合金的Bodner—Partom本构方程

采用Ｂｏｄｎｅｒ－Ｐａｒｔｏｍ粘塑性本构理论研究了Ｓｎ－Ｐｂ共晶合金的本构方程。在应变速率１０－５—１０－２Ｓ－１、温度为－５５—１２５℃的外部变量范围内进行了单轴拉伸和稳态蠕变数值模拟，结果与实验吻合良好。同时讨论了该本构方程用于ＳＭＴ焊点热循环寿命预测的实际意义。     

热粘塑性体的积分—微分型本构关系

应用ИЛЪЮШИН关于应力是五维偏应变空间变形历史的泛函的概念和Ｖａｌａｎｉｓ有关内时理论的描述，本文提出，对热粘塑性体，应力可设为应变、应变率和温度历史泛；并应用Ｍｉｌｌｅｒ和其它一些作者有关内变量演化方程的描述，由此建立了热粘塑性体的积分－微分型本构方程这一积分－微分型本构关系大体和Ｍｉｌｌｅｒ微分型模型等价。对１０２０钢的单轴本构响应进行了数值模拟，和Ｔａｎａｋａ与Ｍｉｌｌｅｒ的分析及一些实     

热粘塑性体的积分－微分型本构关系

应用　　　关于应力是五维偏应变空间变形历史的泛函的概念和Ｖａｌａｎｉｓ有关内＊时理论的描述，本文提出，对热粘塑性体，应力可设为应变、应变率和温度历史的泛函；并应用Ｍｉｌｌｅｒ和其它一些作者有关内变量演化方程的描述，由此建立了热粘塑性体的积分－微分到本构方程．这一积分－微分型本构关系大体和Ｍｉｌｌｅｒ微分型模型等价．对１０２０钢的单轴本构响应进行了数值模拟，和Ｔａｎａｋａ与Ｍｉｌｌｅｒ的分析及一些实验结果符合较好．     

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

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

哑铃式聚合物分子模型流变学数值研究

本文介绍了求解哑铃式分子模型位形空间分布函数扩散方程的数值方法,以及用这种方法计算的若干分子模型的流变性质。在通常情况下,将这一方法与求解流动守恒方程的边界元法相结合,便有可能用一个得不到本构方程的分子模型去代替连续介质力学本构方程,来模拟聚合物流体的复杂流动。本文还讨论了这一方法某些令人感兴趣的特点。      

宾汉流体物性参数对湍流强度影响的研究

考虑宾汉流体本构关系的特点，建立了适用于宾汉流体的修正动量方程和Ｋτ－ετ模型方程。采用压力耦合半隐式ＳＩＭＰＬＥ算法，通过流场的数值计算研究，指出了宾汉流体本身的物性参数（如屈服应力、塑性粘度）对其湍流强度有着重要影响。     

粘弹流体轴对称流动的格子Boltzmann方法

基于BGK碰撞模型，通过在迁移方程中引入作用力项，建立了粘弹流体的轴对称格子Boltzmann模型．通过Chapman-Enskog展开，获得了准确的柱坐标下轴对称宏观流动方程．采用双分布函数对运动方程和本构方程进行迭代求解，模拟分析了粘弹流体管道流动，获得了流场中的速度和构型张量的分布，通过与解析解进行比较，验证了模型的准确性．研究了作为粘弹流体流动基准问题的收敛流动，对涡旋位置进行了定量分析，将回转长度的计算结果与有限体积法进行了比较，两种数值结果十分吻合．研究结果表明，模型能够准确表征粘弹流体的轴对称流动，具有较广阔的应用前景．     

Numerical analysis of complex flow of fibre suspensions in a viscoelastic fluid

Research efforts made so far to simulate fibre suspension flows are limited to fibre suspensions in Newtonian fluids. Though short fibre composites are mostly made of polymers, lack of suitable constitutive equations for fibre suspensions in viscoelastic fluids frustrates attempts to simulate flows of these suspensions. A preliminary work done by the author led to a constitutive equation for semiconcentrated fibre suspensions in the Oldroyd-B fluid. This paper describes the mathematical formulation of the flow problems for such a suspension and the numerical procedure to solve them. Some numerical results of flow past a sphere in a tube are also presented. This paper is supported by the University of Melbourne, the National Natural Science Foundation of China and Zhejiang Province.     

Rheology of semiconcentrated fibre suspension in the Oldroyd-B fluid

The constitutive equation for a semiconcentrated fibre suspension in the Oldroyd-B fluid has been derived from a statistical model of such a suspension by employing the molecular theory for polymeric liquids. To circumvent theoretical difficulties in viscoelastic fluid mechanics, several simplified models are used to account for the interactions of fibres and polymer molecules. Some of material functions are calculated in terms of the constitutive equation. The project Supported by the University of Melbourne of Australia, the National Natural Science Foundation of China and Zhejiang Province     

非稳态分数阶Oldroyd-B流体绕楔形体流动研究

研究了非稳态分数阶Oldroyd-B流体在多孔介质中通过楔形拉伸板的驻点流动问题。基于分数阶Oldroyd-B流体的本构模型建立了动量方程,并在其中引入了浮升力和驻点流动特征。此外,考虑了具有热松弛延迟时间的修正的分数阶Fourier定律,并将其应用于能量方程和对流换热边界条件。接着,采用与L1算法相结合的有限差分法求解控制偏微分方程。最后,分析了相关物理参数对流动的影响。结果表明,随着楔角参数的增加,流体受到的浮升力增大,导致速度加快;达西数越大,介质的孔隙度变大,流体的流动越快;此外,温度分布先略有上升后明显下降,这表明Oldroyd-B流体具有热延迟特性。     

Stokes＇ first problem for a viscoelastic fluid with the generalized Oldroyd-B model

The flow near a wall suddenly set in motion for a viscoelastic fluid with the generalized Oldroyd-B model is studied. The fractional calculus approach is used in the constitutive relationship of fluid model. Exact analytical solutions of velocity and stress are obtained by using the discrete Laplace transform of the sequential fractional derivative and the Fox H-function. The obtained results indicate that some well known solutions for the Newtonian fluid, the generalized second grade fluid as well as the ordinary Oldroyd-B fluid, as limiting cases, are included in our solutions. The project supported by the National Natural Science Foundation of China (10272067), the Doctoral Program Foundation of the Education Ministry of China (20030422046), the Natural Science Foundation of Shandong Province, China (Y2006A14) and the Research Foundation of Shandong University at Weihai. The English text was polished by Keren Wang.     

Three dimensional finite element analysis of the flow of polymer melts

The finite element simulation of a selection of two- and three-dimensional flow problems is presented, based upon the use of four different constitutive models for polymer melts (Oldroyd-B, Rolie-Poly, Pom-Pom and XPP). The mathematical and computational models are first introduced, before their application to a range of visco-elastic flows is described. Results demonstrate that the finite element models used here are able to re-produce predictions made by other published numerical simulations and, significantly, by carefully conducted physical experiments using a commercial-grade polystyrene melt in a three-dimensional contraction geometry. The paper also presents a systematic comparison and evaluation of the differences between two- and three-dimensional simulations of two different flow regimes: flow of an Oldroyd-B fluid around a cylinder and flow of a Rolie-Poly fluid into the contraction geometry. This comparison allows new observations to be made concerning the relatively poor quality of two-dimensional simulations for flows in even quite deep channels.     

Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method

An efficient immersed boundary-lattice Boltzmann method (IB-LBM) is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows. Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor. A boundary thickening-based direct forcing IB method is adopted to solve the particle–fluid interactions with high accuracy for non-slip boundary conditions. A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computational parameters. An asynchronous calculation strategy is adopted to further improve the computing efficiency. The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods. The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed. Subsequently, 504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh–Taylor instabilities was observed. At last, simulations of a dense flow involving 11001 particles, the largest number of particles to date, were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid. The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time. These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.     

Transient viscoelastic helical flows in pipes of circular and annular cross-section

Start-up helical flows for Oldroyd-B and upper-convected Maxwell fluids are studied in straight pipes of circular and annular cross-section. The differential form of the constitutive equation leads to partial differential equations which are second-order in space and time. Apart from the condition that the fluid is initially at rest another initial condition is required to complete the solution process. By comparing results derived from the integral form of the constitutive equation we show that an appropriate initial condition may be found. Numerical results for start-up rotational flow in pipes of annular cross-section are presented.     

Influence of viscoelasticity on drop deformation and orientation in shear flow: Part 1. Stationary states

The influence of matrix and droplet viscoelasticity on the steady deformation and orientation of a single droplet subjected to simple shear is investigated microscopically. Experimental data are obtained in the velocity–vorticity and velocity–velocity gradient plane. A constant viscosity Boger fluid is used, as well as a shear-thinning viscoelastic fluid. These materials are described by means of an Oldroyd-B, Giesekus, Ellis, or multi-mode Giesekus constitutive equation. The drop-to-matrix viscosity ratio is 1.5. The numerical simulations in 3D are performed with a volume-of-fluid algorithm and focus on capillary numbers 0.15 and 0.35. In the case of a viscoelastic matrix, viscoelastic stress fields, computed at varying Deborah numbers, show maxima slightly above the drop tip at the back and below the tip at the front. At both capillary numbers, the simulations with the Oldroyd-B constitutive equation predict the experimentally observed phenomena that matrix viscoelasticity significantly suppresses droplet deformation and promotes droplet orientation. These two effects saturate experimentally at high Deborah numbers. Experimentally, the high Deborah numbers are achieved by decreasing the droplet radius with other parameters unchanged. At the higher capillary and Deborah numbers, the use of the Giesekus model with a small amount of shear-thinning dampens the stationary state deformation slightly and increases the angle of orientation. Droplet viscoelasticity on the other hand hardly affects the steady droplet deformation and orientation, both experimentally and numerically, even at moderate to high capillary and Deborah numbers.     

Direct simulation for concentrated fibre suspensions in transient and steady state shear flows

A numerical method to simulate fibre suspensions in transient and steady state shear flows for concentrated solutions is reported, which takes into account short-range hydrodynamic interaction via lubrication forces, contact forces and hydrodynamic forces. Fibres are assumed to have varying lengths as observed in industrial composites and this can be described with a fibre distribution length function. Stresses in the composite have been calculated from a constitutive equation or from an approximate expression for the effective stress of fibre suspensions. This simulation will assist in modelling of short fibre filled industrial composites. Our first results shown here qualitatively agree with experiments.