纤维悬浮液搅拌流动的数值模拟

由于缺乏适当的本构方程，对纤维悬浮液流动的研究一直局限于纤维的牛顿流体悬浮液。本文采用ＭＵＣＭ模型对作者最近提出的纤维Ｏｌｄｒｏｙｄ－Ｂ流体悬浮液的本构方程作了改进，并对锚式桨搅拌槽的二维Ｏｌｄｒｏｙｄ－Ｂ流体和牛顿流体纤维悬浮液搅拌流动作了数值模拟。模拟的结果表明，本文所用的模型和方法能有效地抑制过大局部应力的影响并合理地处理流体的记忆效应。     

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.     

Der Einfluß der elastischen Eigenschaften des Trägermediums auf das Fließverhalten von Kugel- und Fasersuspensionen

The influence of the elastic properties of the suspending medium on the flow of viscoelastic glass bead and glass fibre suspensions through flat orifices was investigated. The results are discussed by contrasting the flow behaviour of the viscoelastic suspensions with that of corresponding suspensions in a Newtonian suspending medium. For the suspensions in a Newtonian oil linear relationships were always found between the pressure loss and the effective velocity gradient in orifice flow. Thus it can be concluded that in this case the influence of the filler on the flow behaviour is independent of the imposed strain. Increasing the filler content or using more anisotropic particles led to higher viscosities and thus to larger pressure losses. It is well known that viscoelastic polyisobutene solutions show strain rate dependent flow behaviour. Due to the increasing influence of elasticity with increasing strain rate, the apparent flow curves can be divided into characteristic regions with different slopes. The addition of filler to such solutions altered the shape of the flow curves and it was found that the onset of “flow hardening” occurred at lower imposed strains. In addition, characteristic changes in the hardening behaviour and flow stability were observed; these were most pronounced for the fibre suspensions, even at low concentrations. For the fibre suspensions, these phenomena could be related to the influence of the fibres on the undisturbed flow field near the orifice, leading, in general, to higher strain velocities between the fibres. On the other hand, enhanced extensional strains are induced at the ends of the fibres and, at the same time, shear flow occurs along the fibres. These two effects counteract each other with respect to the hardening behaviour of the polymer. Similar effects may also occur in the glass bead suspensions, although they would certainly be less pronounced.     

Resistance to uniaxial extensional flow of fibre suspensions

The tensile stress due to resistance to uniaxial extensional flow of fibre suspensions in Newtonian and non-Newtonian fluids has been measured using the filament stretching technique. It has been found that addition of fibres to a Newtonian fluid increases the extensional viscosity. The steady state results agree with Bachelors theory and the stress growth behaviour is qualitatively predicted by the theory of Dinh and Armstrong. Experimental results from this work have also shown that the behaviour of a fibre suspension in viscoelastic fluid is qualitatively described by Fans equation. The added fibres increase the extensional stress growth coefficient of the viscoelastic fluid at low strain but have marginal effect on the fluid after the onset of strain-hardening.     

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.     

Viscoelastic medium flows in convergent-divergent channels

The problem of viscoelastic fluid flow in a convergent-divergent channel is considered. A mathematical model and the results of numerical investigations are presented. A model of the differential type is used to describe the rheological properties of viscoelastic fluids. A comparative analysis of the flow parameters of generalized Newtonian and viscoelastic fluids in the channels considered is carried out on the basis of the results of numerical calculations.     

Rheology of a suspension of particles in viscoelastic fluids

In this paper we study the bulk stress of a suspension of rigid particles in viscoelastic fluids. We first apply the theoretical framework provided by Batchelor [J. Fluid Mech. 41 (1970) 545] to derive an analytical expression for the bulk stress of a suspension of rigid particles in a second-order fluid under the limit of dilute and creeping flow conditions. The application of the suspension balance model using this analytical expression leads to the prediction of the migration of particles towards the centerline of the channel in pressure-driven flows. This is in agreement with experimental observations. We next examine the effects of inertia (or flow Reynolds number) on the rheology of dilute suspensions in Oldroyd-B fluids by two-dimensional direct numerical simulations. Simulation results are verified by comparing them with the analytical expression in the creeping flow limit. It is seen that the particle contribution to the first normal stress difference is positive and increases with the elasticity of the fluid and the Reynolds number. The ratio of the first normal stress coefficient of the suspension and the suspending fluid decreases as the Reynolds number is increased. The effective viscosity of the suspension shows a shear-thinning behavior (in spite of a non-shear-thinning suspending fluid) which becomes more pronounced as the fluid elasticity increases.     

Flowing colloidal suspensions in non-Newtonian suspending fluids: Decoupling the composite birefringence

The rheology of dilute, colloidal suspensions in polymeric suspending fluids can be studied with simultaneous dichroism and birefringence measurements. The dichroism provides a direct measure of the particle dynamics, but the birefringence is a composite property with independent contributions from the suspended particles and the polymer molecules. For suspensions where the contribution from the particles is significant, the composite birefringence must be decoupled in order to analyze the dynamics of the polymeric suspending fluid. A method to perform the decoupling is derived and then demonstrated through transient shear flow experiments with dilute suspensions ofFeOOH particles in semi-dilute, xanthan gum suspending fluids. The birefringence of the xanthan gum suspending fluid is calculated from experimental measurements of the composite birefringence and the dichroism of the suspension. To gather information on particle/polymer interactions, the calculated birefringence is compared to the birefringence of xanthan gum solutions containing no suspended particles and the dirchoism is compared to that of a suspension in a Newtonian fluid.     

On co-rotational and other rate type constitutive equations

We derive expressions for the dilatational properties of suspensions of gas bubbles in incompressible fluids, using a cell model for the suspension. A cell, consisting of a gas bubble centered in a spherical shell of incompressible fluid, is subjected to a purely dilatational boundary motion and the resulting stress at the cell boundary is obtained. The same dilatational boundary motion is prescribed at the boundary of an “equivalent” cell composed of a one-phase, uniformly compressible fluid with unknown dilatational properties. By specifying that the stress at the boundary of the one-phase cell is equal to the stress at the boundary of the two-phase suspension cell, we obtain expressions for the unknown dilatational properties as a function of observable properties of the suspension. The dilatational viscosity of a suspension with a Newtonian continuous phase and the analogous properties for suspensions with non-Newtonian continuous phases are obtained as functions of the boundary motion, volume fraction of gas, and properties of the incompressible continuous phase. Results are presented for continuous phases which are Newtonian fluids, second-order fluids, and Goddard—Miller model fluids.     

Front‐tracking by the level‐set and the volume penalization methods in a two‐phase microfluidic network

Two‐phase immiscible fluids in a two‐dimensional micro‐channels network are considered. The incompressible Stokes equations are used to describe the Newtonian fluid flow, while the Oldroyd‐B rheological model is used to capture the viscoelastic behavior. In order to perform numerical simulations in a complex geometry like a micro‐channels network, the volume penalization method is implemented. To follow the interface between the two fluids, the level‐set method is used, and the dynamics of the contact line is modeled by Cox law. Numerical results show the ability of the method to simulate two‐phase flows and to follow properly the contact line between the two immiscible fluids. Finally, simulations with realistic parameters are performed to show the difference when a Newtonian fluid is pushed by a viscoelastic fluid instead of a Newtonian one. Copyright © 2015 John Wiley & Sons, Ltd.     

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

The main purpose of this study was to examine the viscous and elastic properties and capillary flow of fiber suspensions in Newtonian fluids as well as in polymer solutions. The fillers used were glass, carbon, nylon and vinylon fibers. Glycerin was used as a Newtonian suspending medium and HEC and Separan solutions as viscoelastic suspending media. The viscosity and the first normal-stress difference were measured using a coaxial cylindrical rotating viscometer and a parallel-plate rheogoniometer respectively. The influence of the concentration, aspect ratio, diameter and flexibility of the fibers on the viscous and elastic properties of the fiber suspensions was investigated. Empirical equations were obtained for the relative viscosity and first normal-stress difference for the fiber suspensions in glycerin. The capillary flow of these suspensions is discussed in part II.     

Dynamic simulation of non-spherical particulate suspensions

Particle-level simulation has been employed to investigate rheology and microstructure of non-spherical particulate suspensions in a simple shear flow. Non-spherical particles in Newtonian fluids are modeled as three-dimensional clusters of neutrally buoyant, non-Brownian spheres linked together by Hookean-type constraint force. Rotne–Prager correction to velocity disturbance has been employed to account for far-field hydrodynamic interactions. An isolated rod-like particle in simple shear flow exhibits a periodic orientation distribution, commonly referred to as Jeffery orbit. Lubrication-like repulsive potential between clusters have been included in simulation of rod-like suspensions at various aspect ratios over dilute to semi-dilute volume fractions. Shear viscosity evaluated by orientation distribution qualitatively agrees with one obtained by direct computation of shear stress.     

A note on the translation of a thin rod inside a cylinder

The problem of a thin rod moving longitudinally along the axis of symmetry of a cylindrical vessel is examined for Newtonian and non-Newtonian liquids. For non-Newtonian fluids, the inelastic power-law type solution predicts the experimental results particularly well. On account of wall effects, the induced pressure gradients are much greater for a Newtonian fluid than for a viscoelastic fluid. In fact, in the latter case, they may be considered negligible when the radius of the inner cylinder is small compared to the one of the outer cylinder.     

Dissipation in a dilute suspension of spheres in a second-order fluid

We examine the effective medium properties of a dilute suspension of spheres in a second-order fluid under linear shear. Since the second-order fluid is the first step toward the general viscoelastic fluid, the results obtained may provide a qualitative feel for the problem in which the suspending fluid obeys a more complicated (and realistic) constitutive relation.The dissipation in the medium is calculated by determining the rate of working by surface forces; this is compared to the dissipation in a homogeneous fluid to give the effective properties. The results show that the term linear in volume fraction increases the corresponding rheological coefficient, just as in the Newtonian case. It is to be noted that the second-order dissipation is zero for simple shear and other weak flows, whereas for strong flows the small correction may increase or decrease the overall dissipation.     

Experimental and numerical study of the rotation and the erosion of fillers suspended in viscoelastic fluids under simple shear flow

When a porous agglomerate immersed in a fluid is submitted to a shear flow, hydrodynamic stresses acting on its surface may cause a size reduction if they exceed the cohesive stress of the agglomerate. The aggregates forming the agglomerate are slowly removed from the agglomerate surface. Such a behaviour is known when the suspending fluid is Newtonian but unknown if the fluid is viscoelastic. By using rheo-optical tools, model fluids, carbon black agglomerates and particles of various shapes, we found that the particles had a rotational motion around the vorticity axis with a period which is independent on shape (flat particles not considered), but which is exponentially increasing with the elasticity of the medium expressed by the Weissenberg number (We). Spherical particles are always rotating for We up to 2.6 (largest investigated We in this study) but elongated particles stop rotating for We>0.9 while orienting along the flow direction. Erosion is strongly reduced by elasticity. Since finite element numerical simulation shows that elasticity increases the local stress around a particle, the origin of the erosion reduction is interpreted as an increase of cohesiveness of the porous agglomerate due to the infiltration of a viscoelastic fluid.     

Viscoelastic effects in double‐pipe single‐pass counterflow heat exchangers

The conducto‐convective heat loss from a viscoelastic liquid, in the core of a double‐pipe heat exchanger arrangement, to a cooler Newtonian fluid flowing in the outer annulus is investigated with direct numerical simulations. A numerical algorithm based on the finite difference method is implemented in time and space with the Giesekus constitutive model for the viscoelastic liquids. The flow of both the annulus and core‐fluids is considered to be Poiseuille flow, driven by respective pressure gradients. In general, the results show that a viscoelastic core‐fluid leads to slightly lower (albeit comparable) attainable temperatures in the core‐fluid stream as compared with a corresponding Newtonian fluid. Copyright © 2008 John Wiley & Sons, Ltd.     

Dispersion of linear gravity waves on a viscoelastic fluid in an horizontal canal

A characteristic equation is derived that describes the spatial decay of linear surface gravity waves on Maxwell fluids. Except at small frequencies, the derived dispersion relation is different from the temporal decay dispersion relation which is normally studied within fluid mechanics. The implications for waves on viscous Newtonian fluids using the two different dispersion relations is briefly discussed. The wave number is measured experimentally as function of the frequency in a horizontal canal. Seven Newtonian fluids and four viscoelastic liquids with constant viscosity have been used in the experiments. The spatial decay theory for Newtonian fluids fits well to the experimental data. The model and experiments are used to determine limits for the Maxwell fluid time numbers for the four viscoelastic liquids. As a result of low viscosity it was not possible within this study to obtain these time numbers from oscillatory experiments. Therefore, a comparison of surface gravity wave experiments with theory is applicable as a method to evaluate memory times of low viscosity viscoelastic fluids.     

Transient flows of a second grade fluid

Exact analytical solutions for a class of unsteady unidirectional flows of an incompressible second-order fluid are constructed. The flows are generated impulsively from rest by motion of a plate or two plates or by sudden application of a pressure gradient. Expressions for velocity, flux and skin friction are obtained for both large and small times. It is found that large and small times solutions are dependent on the coefficient of viscoelasticity. The solutions corresponding to Newtonian fluids can be easily obtained from those for fluids of second order by letting the viscoelastic parameter to be zero.