Very slow flow of Bingham viscoplastic fluid around a circular cylinder

Numerical simulations have been used to study the flow of a Bingham viscoplastic fluid around a circular cylinder in an infinite medium with negligible inertia effects. Papanastasiou's regularisation technique has been adopted to approximate the model. The case corresponding to preponderant plasticity effects has been particularly studied and convergence of the solutions examined in detail. The flow kinematics and stresses have been determined. The rigid zones have been identified and characterised. At large Oldroyd numbers, when plasticity effects become preponderant, a viscoplastic boundary layer appears around the cylinder. The characteristics of this viscoplastic boundary layer are quantified. The results are compared with existing theoretical results, concerning particularly the predictions of the viscoplastic boundary layer theory and the plasticity theory.     

Unsteady shear flow of a viscoplastic material

A viscoplastic, or yield-stress, liquid occupies the space between two infinite parallel plates. Initially the whole system is at rest. The lower plate is suddenly jerked into motion with given speed or shear stress, while the upper plate is kept fixed. The flow consists of two regions; (1) a lower sheared region bounded above by the yield surface, (2) an upper unyielded region bounded below by the yield surface. The yield surface propagates to the upper plate as time proceeds. We first consider the equivalent one plate problem of flow over a jerked plate, and find similarity solutions and small time asymptotic solutions for prescribed shear and speed cases respectively. These solutions are used as initial solutions for the two plate case. The motion of the yield surface and the time taken for the entire material to yield are investigated. The problems considered here are two dimensional representations of some control devices, for example the light duty clutch, which consists of two corotating, coaxial discs separated by a layer of electrorheological material. In this application it is useful to know the time taken for all the material to yield.     

Analytical solutions of the boundary-value problem of nonstationary flow of viscoplastic medium between two plates

Summary Nonstationary flow of a viscoplastic medium between two parallel plates is considered for the case of a varying pressure gradient. The problem is reduced to the Stephan problem, with the condition on the boundary separating the flow domain from the quasi-rigid domain. Four multiparameter families of exact solutions are found. The first family describes the flow decelerations up to a full stop. The second family determines the development of the flow from the state of rest as the pressure gradient increases. The third family describes the development of the flow for the case where (1) the pressure gradient is constant and exceeds the threshold value related to the yield stress, (2) the upper plate does not move, and (3) the lower plate moves with a constant acceleration. Finally, the fourth family determines the flow retardation, when the pressure gradient is constant and is less than the threshold value. The decrease in the flow of the viscoplastic medium can be achieved for certain values of parameters by increasing the quasi-rigid domain, whereas the viscoplastic flow remains unchanged. Received 7 October 1998; accepted for publication 8 April 1999     

Numerical simulations of complex yield-stress fluid flows

Viscoplasticity is characterized by a yield stress, below which the materials will not deform and above which they will deform and flow according to different constitutive relations. Viscoplastic models include the Bingham plastic, the Herschel-Bulkley model and the Casson model. All of these ideal models are discontinuous. Analytical solutions exist for such models in simple flows. For general flow fields, it is necessary to develop numerical techniques to track down yielded/unyielded regions. This can be avoided by introducing into the models a regularization parameter, which facilitates the solution process and produces virtually the same results as the ideal models by the right choice of its value. This work reviews several benchmark problems of viscoplastic flows, such as entry and exit flows from dies, flows around a sphere and a bubble and squeeze flows. Examples are also given for typical processing flows of viscoplastic materials, where the extent and shape of the yielded/unyielded regions are clearly shown. The above-mentioned viscoplastic models leave undetermined the stress and elastic deformation in the solid region. Moreover, deviations have been reported between predictions with these models and experiments for flows around particles using Carbopol, one of the very often used and heretofore widely accepted as a simple “viscoplastic” fluid. These have been partially remedied in very recent studies using the elastoviscoplastic models proposed by Saramito.     

Interaction between two circular cylinders in slow flow of Bingham viscoplastic fluid

The interaction between two circular cylinders was studied in the slow flow of a Bingham viscoplastic fluid in an infinite medium without any inertia effects. The configuration studied is that in which the flow direction is parallel to the centre line of the cylinders. Finite-element numerical simulations were used with an approximation by Papanastasiou's regularisation method. The case of high yield stress effect was particularly examined. The convergence of the solutions was examined in detail. Changes in the rigid zones, kinematics and stresses were determined in relation to the degree of interaction, which is a function of the distance between the cylinders and the effect of yield stress. The results compared with the case of a single cylinder show that yield stress reduces interaction effects. The transition between configurations with interacting cylinders and configuration with isolated cylinders was examined as a function of the effect of yield stress. Correlations were proposed for the drag coefficient and the stability criterion when the cylinders are interacting.     

Computations with viscoplastic and viscoelastoplastic fluids

This numerical study focuses on regularised Bingham-type and viscoelastoplastic fluids, performing simulations for 4:1:4 contraction?Cexpansion flow with a hybrid finite element?Cfinite volume subcell scheme. The work explores the viscoplastic regime, via the Bingham?CPapanastasiou model, and extends this into the viscoelastoplastic regime through the Papanastasiou?COldroyd model. Our findings reveal the significant impact that elevation has in yield stress parameters, and in sharpening of the stress singularity from that of the Oldroyd/Newtonian models to the ideal Bingham form. Such aspects are covered in field response via vortex behaviour, pressure-drops, stress field structures and yielded?Cunyielded zones. With rising yield stress parameters, vortex trends reflect suppression in both upstream and downstream vortices. Viscoelastoplasticity, with its additional elasticity properties, tends to disturb upstream?Cdownstream vortex symmetry balance, with knock-on effects according to solvent-fraction and level of elasticity. Yield fronts are traced with increasing yield stress influences, revealing locations where relatively unyielded material aggregates. Analysis of pressure drop data reveals significant increases in the viscoplastic Bingham?CPapanastasiou case, O (12%) above the equivalent Newtonian fluid, that are reduced to 8% total contribution increase in the viscoelastoplastic Papanastasiou?COldroyd case. This may be argued to be a consequence of strengthening in first normal stress effects.     

Singular solutions in viscoplasticity under plane strain conditions

The paper deals with asymptotic behavior of viscoplastic solutions in the vicinity of maximum friction surfaces under plane strain conditions. The definition of maximum friction surfaces is that the friction stress is equal to the shear yield stress at sliding. The constitutive equations of the viscoplastic model adopted include a saturation stress. It is shown that it is possible to choose parameters of the viscoplastic model such that the regime of sliding is possible at maximum friction surfaces. In this case solutions are singular in the vicinity of such surfaces. Because of this feature of solutions, the viscoplastic model chosen possesses a smooth transition of qualitative behavior between rigid perfectly plastic and viscoplastic solutions, and this may prove to be advantageous for some applications.     

Flow of elasto-viscoplastic liquids through a planar expansion–contraction

The steady flow of incompressible elasto-viscoplastic liquids through a planar expansion–contraction is investigated. A novel constitutive model is employed to describe the mechanical behavior of the flowing liquids. Numerical solutions of the constitutive and conservation equations were obtained via a finite element method to investigate the role of elasticity, yield stress, and inertia. The fields of velocity, stress, elastic strain, and rate of strain were obtained for different combinations of the governing parameters. It was observed that these fields, as well as the shape and position of the yield surface, are all strong functions of elasticity, yield stress, and inertia. The trends observed agree well with previous numerical and visualization results available in the literature. The present work offers a detailed study on the effects of elasticity, presenting, in particular, the fields of elastic strain.     

Stagnant zones in viscoplastic media

The problem of the flow of a viscoplastic medium with linear viscosity, resulting from the motion of a rigid cylinder of arbitrary cross section has been considered by Oldroyd [1]. In [2] the author finds several exact particular solutions of this problem. Problems of a similar nature have also been examined in [3–4], In the following we consider the formation of stagnant zones in viscoplastic media.     

Numerical approximations for flow of viscoplastic fluids in a lid-driven cavity

The effect of inertia and rheology parameters on the flow of viscoplastic fluids inside a lid-driven cavity is investigated using a stabilized finite element approximation. The viscoplastic material behavior is described by the model introduced by de Souza Mendes and Dutra [30] – herein called SMD fluid – which is essentially a regularized viscosity function that involves only rheological properties of the material. The incompressible balance equations are coupled with the non-linear SMD model and are approximated by a multi-field Galerkin least-squares method in terms of extra-stress, pressure and velocity. The results obtained confirm the stability features of the multi-field formulation and the appropriateness of the rheological stress regularization introduced by the SMD fluid. The influence of inertia and rheological parameters on the morphology of the material yield surfaces is analyzed and discussed.     

The stability of two stratified non-newtonian liquids in couette flow

Consideration is given to the stability of the interface between two Oldroyd liquids with shear-dependent viscosities, flowing in distinct layers while undergoing plane Couette flow. Results are presented as regions of stability in the plane determined by the logarithms of the viscosity and depth ratios. The work of previous authors for two Newtonian, power-law and constant-viscosity Oldroyd liquids is revealingly presented in a similar fashion. It is found that the dependence of the viscosities on shear-rate can drastically affect the regions of interfacial stability in a way over and above that due to just a change in the effective viscosity ratio. It is also found that for the Oldroyd liquids this viscosity variation affects the stability when it is present in the less-viscous layer.     

Plane Couette flow of viscoplastic materials along a slippery vibrating wall

This study looks at the influence of slip at the wall on plane Couette flows of viscous and yield stress fluids with ultrasonic wall motion. These fluids are used in coating processes. A constant speed V at one wall creates the flow, and vibrations and slip take place at the other wall. Isothermal conditions and arbitrary (longitudinal or transverse) vibrations are considered, with negligible vibrational inertia.For the Bingham model, due to its nonlinearity, whatever the vibration direction and the wall slipperiness, significant decreases occur in the average stress as soon as moderate values of the dimensionless vibration velocity amplitude are involved. Such effects are associated with adherent or slippery walls, even with linear friction laws. They do not occur with linear viscous (Newtonian) models.Average stress reductions can reach nearly 100% for very high Oldroyd numbers, i.e. for stress values without vibration close to the yield limit. Slip velocity also decreases. The cost in terms of the power dissipated remains relatively less than in the Newtonian case, and may contribute to a change in the temperature field. Even when the flow without vibration is a pure slip one, large enough amplitude vibrations, either longitudinal or transverse, applied at the wall can reduce the average shear stress and slip velocity, giving rise to an average axial shear flow.Hence vibrations of moderate or high-velocity amplitude applied to adherent or slippery walls enhance plane Couette flow rates for viscoplastic materials. With moderate values of this amplitude, longitudinal vibrations may be 1.5–2 times more efficient than transverse vibrations with an equivalent cost. However, if for technological reasons transverse vibrations have to be preferred, they can also produce significant results. In any case, coating flows should benefit from an adequate application of ultrasound at the wall.     

Progress in numerical simulation of yield stress fluid flows

Numerical simulations of viscoplastic fluid flows have provided a better understanding of fundamental properties of yield stress fluids in many applications relevant to natural and engineering sciences. In the first part of this paper, we review the classical numerical methods for the solution of the non-smooth viscoplastic mathematical models, highlight their advantages and drawbacks, and discuss more recent numerical methods that show promises for fast algorithms and accurate solutions. In the second part, we present and analyze a variety of applications and extensions involving viscoplastic flow simulations: yield slip at the wall, heat transfer, thixotropy, granular materials, and combining elasticity, with multiple phases and shallow flow approximations. We illustrate from a physical viewpoint how fascinating the corresponding rich phenomena pointed out by these simulations are.     

The Rayleigh and Stokes problems with an incompressible non-Newtonian fluid

Summary The Rayleigh problem or impulsive motion of a flat plate has been solved using a perturbation scheme when the surrounding fluid is representable by the constitutive equations of Oldroyd or Coleman and Noll. The shear stress and normal stress at the wall were expressed analytically for this unsteady motion. Further, an exact solution of the equations was found for a special case of the constitutive equations.The motion of the fluid above a harmonically oscillating plate or the Stokes problem has been determined for a special non-Newtonian fluid. The penetration of the shear wave into the fluid, the energy dissipation, the velocity profiles and the shear and normal stresses at the wall were expressed and compared to an equivalent Newtonian fluid.Some of the features of these non-Newtonian fluids were examined in simple shearing flows, and techniques to calculate some of the material constants discussed.     

液晶高分子各向异性粘弹性流体本构方程理论

将液晶高分子－各向异性流体的本构方程,建立在Oldroyd随体导数观点基础上。推广上随机Oldroyd B流体模型,提出共转OldroydB流体模型,同时将微观结构的影响通过宏观参数表示出来,使在宏观理论中包含微观结构的贡献,即引入取向物质函数,非线性各向异性黏度函数和各向异性松弛时间及推迟时间等,表征取向运动对黏度和松弛及推迟现象的影响,在此基础上开展了一类新的液晶高分子－Oldroyd型本构方程理论,由该类型本构方程得出的物质函数,液晶高分子流体的第一、第二法向应力差与实验结果一致,解释了液晶高分子溶液的第一、第二法向应力差的特殊流变学行为。     

Extensional dynamics of viscoplastic filaments: II. Drips and bridges

A model for the dynamics of slender filaments of Herschel–Bulkley fluid is used to explore viscoplastic dripping under gravity and thinning under controlled extension (liquid bridges). The conditions required for fluid to yield are delineated, and the subsequent thinning and progression to pinch-off are tracked numerically. Calculations varying the dimensionless parameters of the problem are presented to illustrate the effect of surface tension, rheology, inertia (for dripping) and gravity. The theoretical solutions are compared with laboratory experiments using aqueous solutions of Carbopol and Kaolin suspensions. For drips and bridges, experiments with Carbopol are well matched by the theory, using a surface tension equal to that of water, even in situations when the fluid is not slender. Experiments with Kaolin do not compare well with theory for physically plausible values of the surface tension. Implications for rheometry and surface-tension inference are discussed.     

Transonic Euler computation in streamfunction co-ordinates

Numerical simulation by a finite element method is used to examine the problem of the rotating flow of a viscoelastic fluid in a cylindrical vessel agitated with a paddle impeller. The mathematical model consists of a viscoelastic constitutive equation of Oldroyd B type coupled to the hydrodynamic equations expressed in a rotating frame. This system is solved by using an unsteady approach for velocity, pressure and stress fields. For Reynolds numbers in the range 0.1–10, viscoelastic effects are taken into account up to a Deborah number De of 1.33 and viscoelasticity and inertia cross-effects are studied. Examining the velocity and stress fields as well as the power consumption, it is found that their evolutions are significantly different for low and moderate inertia. These results confirm the trends of experimental studies and show the specific contribution of elasticity without interference of the pseudoplastic character found in actual fluids.     

Turbulent Flow Around Fluid–Porous Interfaces Computed with a Diffusion-Jump Model for <Emphasis Type="Italic">k</Emphasis> and <Emphasis Type="Italic">ε</Emphasis> Transport Equations

Flow over vegetation and bottom of rivers can be characterized by some sort of porous structure of irregular surface through which a fluid permeates. Also, in engineering systems, one can have components that make use of a working fluid flowing over irregular layers of porous material. This article presents numerical solutions for such hybrid medium, considering here a channel partially filled with a flat porous layer saturated by a fluid flowing in turbulent regime. One unique set of transport equations is applied to both the regions. A diffusion-jump model for both the turbulent kinetic energy and its dissipation rate, across the interface, is presented and discussed upon. The discretization steps taken for numerically accommodating such model in the system of algebraic equations are presented. Numerical results show the effects of Reynolds number, porosity, and permeability on mean and turbulence fields. Results indicate that when negative values for the stress jump coefficient are applied, the peak of the turbulent kinetic energy distribution occurs at the macroscopic interface.     

Lower bound for LCF lifetime and its application to safe design of elastic viscoplastic structures

The theory of shakedown is applied to obtain an upper estimation of LCF lifetime of structures. A model of elastic viscoplastic material similar to the Perzyna one with isotropic strain hardening and isotropic damage is adopted. Assumptions: viscoplastic strain rate is proportional to the access of the yield function over zero; the rate of damage evolution is equal to a function of hardening and damage parameters with the coefficient of fluidity, as a factor of proportionality; damage process is coupled with the viscoplastic deformation process; the hardening parameter is equal to accumulated viscoplastic deformation. The yield surfaces form a family of self-similar surfaces with the diameter as the parameter. The shakedown condition of the Melan type is formulated relatively to the initial yield surface. Features of the stress path lead to an equation with min–max problem of the mathematical programming in the left side, which determines a safe value of the virtual residual stress. The equation provides an opportunity to compute the maximal value of the strain hardening parameter possible under the prescribed loading program. This value allows to obtain an upper estimate to safe work time of the structure, which results in a sufficient condition of the structure integrity during the prescribed time period. An example of the developed theory application to resolve various problems arising from designing of structures is considered.     

Extensional dynamics of viscoplastic filaments: I. Long-wave approximation and the Rayleigh instability

We derive an asymptotic reduced model for the extensional dynamics of long, slender, axisymmetric threads of incompressible Herschel–Bulkley fluids. The model describes the competition between viscoplasticity, gravity, surface tension and inertia, and is used to explore the viscoplastic Rayleigh instability. A finite-amplitude initial perturbation is required to yield the fluid and initiate capillary-induced thinning. The critical amplitude necessary for thinning depends on both the wavelength of the perturbation and on the yield stress. We also numerically examine the inertialess growth of the instability and the progression towards pinch-off. The final self-similar form of inertialess pinch-off is similar to that for a power-law fluid.