同心旋转圆柱间粘弹性流的非线性动力学模型

基于小间隙假设,将速度场和应力场有杉Ｆｏｕｒｉｅｒ展开式截断近似,推导了同心旋转圆柱间Ｏｌｄｒｏｙｄ－Ｂ型流体的六维动力系统,探讨了高分子添加对滑动轴承间油膜稳定性的影响。结果表明,少量的高分子添加剂具有推迟流体层流失稳的作用。     

非牛顿流体力学研究的若干进展

简要介绍微极性流体,应力偶流体,非整数阶Maxwell流体和智能流体的基本概念,通过简单的例子阐述它们与牛顿流体的本质差别,并介绍非牛顿流体在以上热点研究领域的若干进展．     

狭窄血管内偶应力流体流动分析

本文考察了血管狭窄对血液流动的影响,血液以偶应力流体表示,并在求解过程中采用了在管壁上流体质点无相对涡量的边界条件,结果表明,和Young的经典工作相比流动阻抗和壁切应力大于同样程度狭窄下牛顿流体的相应值,偶应力流体对狭窄的敏感性大于牛顿流体;在狭窄发展过程中,偶应力流体的流量要小于牛顿流体的流量,和牛顿流体相比,这些结果更符合生理实际。     

弹性湍流及其对微混合效率影响研究的若干进展

介绍和评述弹性湍流的产生及其对于微混合效率的影响等问题上的若干研究进展. 弹性 湍流和惯性湍流具有类似的流场特征,但引发机理有所不同. 惯性湍流产生的原因 是惯性引起的Reynolds应力,而弹性湍流则是由弹性应力所引起的. 鉴于在微流动 中,惯性力可忽略不计,因此牛顿流体的混合变得十分困难. 此时可在流体中加入 微量高分子聚合物以生成黏弹性流体. 由黏弹性流体所引发的弹性湍流在提高流体 微混合的效率上可发挥重要作用.     

基于舌/上颚微间隙下流体流动行为研究

为探讨口腔环境下流体的流动行为，采用数值方法与流变试验深入研究舌/上颚微间隙下流体流量的影响因素. 建立舌/上颚微间隙的简化模型及Reynolds方程，通过数值方法获取微间隙下流量变化；在DHR-2流变仪上研究非牛顿流体的黏度与剪切率的变化，探讨牛顿流体和非牛顿流体的流量影响. 结果表明:牛顿流体流量平方的倒数同载荷和黏度比值和时间均呈线性函数关系；所制备的非牛顿流体近似为幂律流体，其黏度随脂肪含量的增加而增大，而非牛顿流体流量率先高于后低于等效牛顿流体，其研究结果将为特定人群功能产品的研发提供技术支持.      

基于分形导数对非牛顿流体层流的数值研究

非牛顿流体具有复杂的流变特性,揭示该流变特性可以更加合理地指导非牛顿流体在工农业生产中的应用.经典的非牛顿流体本构模型往往形式复杂,仅能应用于某些特定的情况.分数阶导数模型具有参数少和形式简单的特点,己成功地应用于描述非牛顿流体的运动.Hausdorff分形导数作为一个备选的建模方法,相比分数阶导数具有更简单的形式以及更高的计算效率.本文基于Hausdorff分形导数改进现有牛顿黏性模型,提出分形黏壶模型.通过研究分形黏壶在常应变率下表观黏度的变化情况,以及在加、卸载条件下的蠕变及恢复特性,发现分形黏壶模型适合于描述具有黏弹性的非牛顿流体(本文称之为分形流体).结合连续性方程及运动微分方程,推导出分形流体在平行板间层流的基本方程.按是否拖动上板和是否存在水平的压力梯度分为3种工况,分别用数值方法计算这3种工况下流速在板间的分布及其随时间变化的情况.通过分析不同工况下的流速分布,发现水平的压力梯度会改变流速随时间变化的形状,且会推迟流速到达稳定的时间.在水平压力梯度不存在的情况下,不同阶数的分形流体具有相同的流速分布或是演变过程.另外,在水平压力梯度存在的情况下,上板速度不影响不同阶数分形流体间稳定速度的差值.     

填隙二阶流体下圆球平行于壁平移的粘性阻力

离散元法是分析散体力学行为的数值方法。存在填隙流体时，颗粒之间或颗粒与壁之间产生的法向挤压力和切向阻力、阻力矩，是湿颗粒离散元法的理论基础。二阶流体是以微小偏离牛顿流体本构而考虑时间影响的一种流体。它具有常粘度，并且第一和第二法向应力差正比于剪切率的平方。根据Reynolds润滑理论，采用小参数法，导出了存在填隙二阶流体时，圆球沿平行于平壁缓慢移动时流体的速度场和压力方程，进而求出切向阻力和阻力矩的解析解。有趣的是在推导时所得的速度场和压力方程形式比牛顿流体要复杂得多，但最终结果表明圆球沿平行于平壁移动时因填隙二阶流体引起的切向阻力和阻力矩与牛顿流体时的结果相同。     

存在填隙幂律流体时圆球间切向作用的近似解

在Reynolds润滑理论的基础上，导出了存在填隙幂律流体时，两刚性圆球有相对切向运动时流体压力的近似方程，并进一步求得了圆球所受阻力矩的近似积分表达式，给出了问题的数值解。结果发现幂流体的幂指数、颗粒间最小间隙以及颗粒的大小都是流体压力和颗粒间切向阻力的重要因素。与Godman等人的牛顿流体下圆球平行于平面缓慢运动时圆球所受阻力以及阻力矩的渐近解作的比较表明，本文的数值解优于渐近解。     

基于虚拟区域法的黏弹性流体中微生物游动的数值模拟和应用

微生物是自然生态系统的重要组成部分,掌握微生物在复杂流体中的运动特性可以为微型器件的设计制造提供理论指导.壁面效应是微生物游动研究中的重要问题之一,已有研究表明微生物在壁面附近存在复杂的行为特征.然而已有研究大多集中于微生物在牛顿流体中的游动模拟,仅有少数涉及黏弹性流体等非牛顿流体.本文采用直接力虚拟区域法与乔列斯基分解相结合的数值方法,引入Squirmer微生物游动模型,研究了微生物在黏弹性流体中的游动问题.首先给出求解黏弹性流体本构方程的数值格式;并将该方法应用于研究微生物游动中的壁面效应.研究结果表明,游动方向是影响微生物颗粒壁面效应的重要因素.流体弹性应力会对微生物产生一个反向转矩,影响微生物的游动方向,从而阻碍微生物逃离壁面.微生物颗粒在黏弹性流体中与壁面作用时间较长,几乎达到牛顿流体的两倍以上.     

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

基于Oldroyd-B型粘弹性流体模型，采用同心旋转圆柱间非线性动力系统分析了流体的弹性对轴对称Taylor涡稳定性的影响．分析结果表明，对于弱弹性流体，Taylor涡出现时，系统存在超临界分岔；而对于强弹性流体则出现亚临界分岔．在小间隙大扰动条件下，采用有限差分法分析了非线性效应对系统稳定性的影响．数值计算结果表明，随着流动速度的增加，润滑油膜的失稳结构与流体的弹性有关，对于弱弹性流，流体以同宿轨道分岔失稳；强弹性流则出现倍周期分岔，直至发生混沌，流场最终发展为湍流．     

Experimental study of viscoelastic effectives in a cylinder-plane lubricated contact

This experimental study is related to the flow between a stationary small diameter needle and a moving drum of much larger diameter, seen as a moving plane. Our visualisation experimental set-up adds new results regarding the flow behaviour of a polymer solution. This flow is compared to that of a Newtonian fluid used under the same kinematic conditions. Important differences between Newtonian and polymeric flows concern mainly streamlines, stagnation point positions, and reverse flow. A second experimental set-up enabled us to obtain normal stress profiles along the needle wall and the flow rate through the gap. Reslts related to Newtonian flows are compared to the well-known analytical solutions for a bidimensional laminar flow.     

Electrokinetic Mixing and Displacement of Charged Droplets in Hydrogels

Mixing in droplets is an essential task in a variety of microfluidic systems. Inspired by electrokinetic mixing, electric field-induced hydrodynamic flow inside a charged droplet embedded in an unbounded polyelectrolyte hydrogel is investigated theoretically. In this study, the polyelectrolyte hydrogel is modeled as a soft, and electrically charged porous solid saturated with a salted Newtonian fluid, and the droplet is considered an incompressible Newtonian fluid. The droplet-hydrogel interface is modeled as a surface, which is located at the plane of shear, with the electrostatic potential \(\zeta \) . The fluid inside the droplet attains a finite velocity owing to hydrodynamic coupling with the electroosmotic flow arising from the droplet and polymer charge. The fluid velocity inside the droplet is linearly proportional to the electroosmotic flow velocity in the charged gel and the electroosmotic flow velocity beyond the electrical double layer of a charged interface. It is found that the polymer boundary condition at the droplet surface and the viscosities of the fluids inside and outside the droplet significantly modulate the interior fluid flow. The ionic strength and the permeability of the polymer network impact the flow differently depending on whether the flow arises from the droplet or polymer charge. Finally, the displacement of a charged droplet embedded in a gel under the influence of an external electric field is undertaken. This work is motivated by experimental attempts, which can register sub-nanometer-scale inclusion displacements in hydrogels, to advance electrical microrheology as a diagnostic tool for probing inclusion-hydrogel interfaces. In the absence of polymer charge, a close connection is found between the electrical response of a charged droplet when it is immobilized in an uncharged incompressible gel and when it is dispersed in a Newtonian electrolyte.     

Selective withdrawal of polymer solutions: Experiments

Selective withdrawal refers to the process of drawing one or both components of stratified fluids through a tube placed near their interface. This paper reports an experimental study of selective withdrawal of viscous and viscoelastic liquids under air. The key mechanism of interest is how the viscoelasticity in the bulk liquid affects the evolution of the free surface. This is investigated by comparing the interfacial behavior between a Newtonian silicone oil and two dilute polymer solutions. While the surface undergoes smooth and gradual deformation for Newtonian liquids, for the polymer solutions there is a critical transition where the surface forms a cusp from which an air jet emanates toward the suction tube. This transition shows a hysteresis when the flow rate or location of the tube is varied. In the subcritical state, the surface of polymer solutions deform much more than its Newtonian counterpart but the deformation is more localized. The interfacial behavior of the polymer solutions can be attributed to the large polymer stress that develops under the surface because of predominantly extensional deformation.     

Laminar,transitional and turbulent annular flow of drag-reducing polymer solutions

Mean and rms axial velocity-profile data obtained using laser Doppler anemometry are presented together with pressure-drop data for the flow through a concentric annulus (radius ratio κ = 0.506) of a Newtonian (a glycerine–water mixture) and non-Newtonian fluids—a semi-rigid shear-thinning polymer (a xanthan gum) and a polymer known to exhibit a yield stress (carbopol). A wider range of Reynolds numbers for the transitional flow regime is observed for the more shear-thinning fluids. In marked contrast to the Newtonian fluid, the higher shear stress on the inner wall compared to the outer wall does not lead to earlier transition for the non-Newtonian fluids where more turbulent activity is observed in the outer wall region. The mean axial velocity profiles show a slight shift (～5%) of the location of the maximum velocity towards the outer pipe wall within the transitional regime only for the Newtonian fluid.     

The flow behavior of fiber suspensions in Newtonian fluids and polymer solutions.

This is the second part of a study examining the mechanical properties and capillary flow of fiber suspensions in Newtonian fluids and in polymer solutions. In part I results for the viscous and elastic properties of the fiber suspensions were presented and it was shown that the fiber suspensions exhibited normal stresses in Newtonian as well as in viscoelastic suspending media. It was thus expected that circulating secondary flows would occur near the entrance to a capillary. Four types of fillers (glass, carbon, nylon and vinylon fibers) suspended in glycerin, HEC solutions and Separan solutions were investigated. The entrance flow patterns were visualized and the pressure fluctuations measured. The visualization enabled the eddies occurring in the fiber suspensions in Newtonian fluids to be analysed and classified into two tpyes. The results from the flow visualization were correlated with the pressure fluctuations. Empirical equations for the tube length correction factor due to the elasticity were obtained.     

A convenient way of interpreting steady shear rheo-optical data of semi-dilute polymer solutions

Rheo-mechanical and rheo-optical investigations were carried out with the aim of determining the influence of deformation and orientation or disentangling of polymer coils on the flow behavior in the non-Newtonian region of the flow curve, for a moderately concentrated network solution. To avoid the influence of polydispersity this was done on a series of narrowly distributed polystyrene standards (dissolved in toluene). By using steady state shear flow measurements it was possible to detect qualitatively a reduction in the entanglement density within the non-Newtonian flow region. Birefringence experiments were able to show that deformation of the polymer coils also occurs in the Newtonian flow region, which has no effect on the flow behavior in this range, whereas in the non-Newtonian flow region the increase in deformation is lower than in the Newtonian range. The flow birefringence and its orientation can be described over the whole range of the flow curve with a newly developed equation system (Eq. 8 and 14) derived from the stress states of a sheared solution using the stress-optical rule. Starting from these equations, it could be shown, that in the Newtonian flow region a mastercurve in form of a reduced birefringence Δn′/η₀=f(γ˙) and a reduced orientation φ= f(γ˙/γ˙ _crit) can be plotted, independent from concentration and molar mass. A comparison of the experimentally determined orientation angle and birefringence curve form with theoretical deformations and orientations of polymer coils in a solution state, without intermolecular interactions, was able to demonstrate that the flow behavior of a moderately concentrated network solution is determined decisively (approximately to 85%) by the disentanglement. Received: 8 May 2000 Accepted: 12 September 2000     

Turbulent expansion flow of low molecular weight shear-thinning solutions

A Laser-Doppler anemometer and a pressure transducer were used to carry out detailed measurements of the mean and root mean square of the velocity and wall-pressure in an axisymmetric sudden expansion flow, with 0.4 and 0.5% by weight shear-thinning aqueous solutions of a low molecular weight polymer (6,000), after appropriate rheological characterisation. In spite of their very low molecular weight, these solutions still exhibited elongational elastic effects through drag reduction of up to 35% relative to Newtonian turbulent pipe flow, as shown by Pereira and Pinho (1994). The results showed small variations of the recirculation bubble length with polymer concentration and Reynolds number and reductions of the normal Reynolds stresses of up to 30%, especially in the tangential and radial directions. The reduction in normal Reynolds stresses within the shear layer is an elongational elastic effect, but this elasticity needs to be considerably more intense, such as with high molecular weight polymers, in order to strongly affect the mean flow characteristics. The observed mean flow patterns with these low molecular weight polymer solutions were indeed similar to those exhibited by Newtonian and inelastic fluids.     

Surface tension driven oscillations of a bubble in a viscoelastic liquid

Small amplitude surface tension driven oscillations of a spherical bubble in a dilute polymer solution are considered. The rheological properties of the liquid are modelled by using a 3-constant constitutive equation of the Oldroyd type. The Laplace transform of the solution of the initial value problem is inverted numerically. As in the Newtonian fluid case, both a discrete and a continuous spectrum occurs. In addition to the non-dimensional parameters in the corresponding problem for a Newtonian fluid, the results depend on two other parameters: the ratio of the relaxation time of the polymer solution and the time scale of the flow (the Deborah number) and the product of the polymer concentration and the intrinsic viscosity. For small bubbles in an aqueous solution having a small relaxation time, significant additional damping is found even for dilute solutions.     

Models for the deformation of a single ellipsoidal drop: a review

Dilute polymer blends and immiscible liquid emulsions are characterized by a globular morphology. The dynamics of a single drop subjected to an imposed flow field has been considered to be a valuable model system to get information on dilute blends. This problem has been studied either theoretically by developing exact theories for small drop deformations or by developing simplified models often based on phenomenological assumptions. In this paper, a critical overview of the available models for the dynamics of a single drop is presented, discussing four different systems, namely the Newtonian system, where a single Newtonian drop is immersed in an infinite Newtonian matrix; the non-Newtonian system, where at least one of the components, the drop fluid or the matrix one, is non-Newtonian; the confined Newtonian system, where the matrix is confined and wall effects alter the drop dynamics; and the confined non-Newtonian system.