Search

We investigate the effect of pressure-dependent wall slip on the steady Newtonian annular Poiseuille flow employing Navier's slip law with a slip parameter that varies exponentially with pressure. The dimensionless governing equations and accompanying auxiliary conditions are solved analytically up to second order by implementing a regular perturbation scheme in terms of the small dimensionless pressure-dependence slip parameter. An explicit formula for the average pressure drop, required to maintain a constant volumetric flowrate, is also derived. This is suitably post-processed by applying a convergence acceleration technique to increase the accuracy of the original perturbation series. The effects of pressure-dependent wall slip are more pronounced when wall slip is weak. However, as the slip coefficient increases, these effects are moderated and eventually eliminated as the perfect slip case is approached. The results show that the average pressure drop remains practically constant until the Reynolds number becomes sufficiently large. It is worth noting that all phenomena associated with pressure-dependent wall slip are amplified as the annular gap is reduced., Comment: 19 pages, 9 figures

This study revisits the development of viscoplastic flow in pipes and channels, focusing on the flow of a Bingham plastic. Using finite element simulations and the Papanastasiou regularisation, results are obtained across a range of Reynolds and Bingham numbers. The novel contributions of this work include: (a) investigating a definition of the development length based on wall shear stress, a critical parameter in numerous applications; (b) considering alternative definitions of the Reynolds number in an effort to collapse the development length curves into a single master curve, independent of the Bingham number; (c) examining the patterns of yielded and unyielded regions within the flow domain; and (d) assessing the impact of the regularisation parameter on the accuracy of the results. The findings enhance the existing literature, providing a more comprehensive understanding of this classic flow problem.

The effect of wall slip on the apparent flow curves of viscoplastic materials obtained using torsional parallel plate rheometers is analysed by considering Herschel-Bulkley fluids and assuming that slip occurs above a critical wall shear stress, the slip yield stress {\tau}_c, taken to be lower than the yield stress, {\tau}_0. Thus, different flow regimes are encountered as the angular velocity of the experiment is increased. When the rim shear stress {\tau}_R is below the slip yield stress, the exerted torque is not sufficient to rotate the disk and the material remains still. When {\tau}_c<{\tau}_R<{\tau}_0 the material is still unyielded but exhibits wall slip and rotates as a solid at half the angular velocity of the rotating disk. Finally, when {\tau}_R>{\tau}_0, the material exhibits slip everywhere and yields only in the annulus r_0

An efficient numerical approach based on weighted average finite differences is used to solve the Newtonian plane Couette flow with wall slip, obeying a dynamic slip law that generalizes the Navier slip law with the inclusion of a relaxation term. Slip is exhibited only along the fixed plate, and the motion is triggered by the motion of the other plate. Three different cases are considered for the motion of the moving plate, i.e., constant speed, oscillating speed, and a single-period sinusoidal speed. The velocity and the volumetric flow rate are calculated in all cases and comparisons are made with the results of other methods and available results in the literature. The numerical outcomes confirm the damping with time and the lagging effects arising from the Navier and dynamic wall slip conditions and demonstrate the hysteretic behavior of the slip velocity in following the harmonic boundary motion., Comment: 21 pages, 15 figures

We provide an explicit analytical solution of the planar Poiseuille flow of a viscoplastic fluid governed by the constitutive equation proposed by De Kee and Turcotte (Chem. Eng. Commun. 6 (1980) 273-282). Formulae for the velocity and the flow rate are derived, making use of the Lambert W function. It is shown that a solution does not always exist because the flow curve is bounded from above and hence the rheological model can accommodate stresses only up to a certain limit. In fact, the flow curve reaches a peak at a critical shear rate, beyond which it exhibits a negative slope, giving rise to unstable solutions.

We propose a methodology for the rheological characterization of a semisolid metal slurry using experimental squeeze flow data. The slurry is modeled as a structural thixotropic viscoplastic material, obeying the regularized Herschel-Bulkley constitutive equation. All rheological parameters are assumed to vary with the structure parameter that is governed by a first-order kinetics accounting for the material structure breakdown and build-up. The squeeze flow is simulated using finite elements in a Lagrangian framework. The evolution of the sample height has been studied for wide ranges of the Bingham and Reynolds numbers, the power-law exponent as well as the kinetics parameters of the structure parameter. Systematic comparisons have been carried out with available experimental data on a semisolid aluminium alloy (A356), where the sample is compressed from its topside under a specified strain of 80% at a temperature of 582 oC while the bottom side remains fixed. Excellent agreement with the experimental data could be achieved provided that at the initial instances (up to 0.01s) of the experiment the applied load is much higher than the nominal experimental load and that the yield stress and the power-law exponent vary linearly with the structure parameter. The first assumption implies that a different model, such as an elastoviscoplastic one, needs to be employed during the initial stages of the experiment. As for the second one, the evolution of the sample height can be reproduced allowing the yield stress to vary from 0 (no structure) to a maximum nominal value (full structure) and the power-law exponent from 0.2 to 1.4, i.e., from the shear-thinning to the shear-thickening regime. These variations are consistent with the internal microstructure variation pattern known to be exhibited by semisolid slurries., Comment: 33 pages, 21 figures, 42 references

We study the flow of a thixotropic fluid around a cylinder. The rheology of the fluid is described by means of a structural viscoplastic model based on the Bingham constitutive equation, regularised using the Papanastasiou regularisation. The yield stress is assumed to vary linearly with the structural parameter, which varies from zero (completely broken structure) to one (fully developed skeleton structure), following a first-order rate equation which accounts for material structure break-down and build-up. The results were obtained numerically using the Finite Element Method. Simulations were performed for a moderate Reynolds number of 45, so that flow recirculation is observed behind the cylinder, but vortex shedding does not occur. The effects of the Bingham number and of the thixotropy parameters are studied. The results show that the viscous character of the flow can be controlled within certain limits through these parameters, despite the fact that the Reynolds number is fixed.

Numerical simulations of the flow in an extrusion damper are performed using a finite volume method. The damper is assumed to consist of a shaft, with or without a spherical bulge, oscillating axially in a containing cylinder filled with a viscoplastic material of Bingham type. The response of the damper to a forced sinusoidal displacement is studied. In the bulgeless case the configuration is the annular analogue of the well-known lid-driven cavity problem, but with a sinusoidal rather than constant lid velocity. Navier slip is applied to the shaft surface in order to bound the reaction force to finite values. Starting from a base case, several problem parameters are varied in turn in order to study the effects of viscoplasticity, slip, damper geometry and oscillation frequency to the damper response. The results show that, compared to Newtonian flow, viscoplasticity causes the damper force to be less sensitive to the shaft velocity; this is often a desirable damper property. The bulge increases the required force on the damper mainly by generating a pressure difference across itself; the latter is larger the smaller the gap between the bulge and the casing is. At high yield stresses or slip coefficients the amount of energy dissipation that occurs due to sliding friction at the shaft-fluid interface is seen to increase significantly. At low frequencies the flow is in quasi steady state, dominated by viscoplastic forces, while at higher frequencies the fluid kinetic energy storage and release also come into the energy balance, introducing hysteresis effects.

We extend our recent work on the creeping flow of a Bingham fluid in a lid-driven cavity, to the study of inertial effects, using a finite volume method and the Papanastasiou regularisation of the Bingham constitutive model [J. Rheology 31 (1987) 385-404]. The finite volume method used belongs to a very popular class of methods for solving Newtonian flow problems, which use the SIMPLE algorithm to solve the discretised set of equations, and have matured over the years. By regularising the Bingham constitutive equation it is easy to extend such a solver to Bingham flows since all that this requires is to modify the viscosity function. This is a tempting approach, since it requires minimum programming effort and makes available all the existing features of the mature finite volume solver. On the other hand, regularisation introduces a parameter which controls the error in addition to the grid spacing, and makes it difficult to locate the yield surfaces. Furthermore, the equations become stiffer and more difficult to solve, while the discontinuity at the yield surfaces causes large truncation errors. The present work attempts to investigate the strengths and weaknesses of such a method by applying it to the lid-driven cavity problem for a range of Bingham and Reynolds numbers (up to 100 and 5000 respectively). By employing techniques such as multigrid, local grid refinement, and an extrapolation procedure to reduce the effect of the regularisation parameter on the calculation of the yield surfaces (Liu et al. J. Non-Newtonian Fluid Mech. 102 (2002) 179-191), satisfactory results are obtained, although the weaknesses of the method become more noticeable as the Bingham number is increased.

We investigate the performance of the finite volume method in solving viscoplastic flows. The creeping square lid-driven cavity flow of a Bingham plastic is chosen as the test case and the constitutive equation is regularised as proposed by Papanastasiou [J. Rheol. 31 (1987) 385-404]. It is shown that the convergence rate of the standard SIMPLE pressure-correction algorithm, which is used to solve the algebraic equation system that is produced by the finite volume discretisation, severely deteriorates as the Bingham number increases, with a corresponding increase in the non-linearity of the equations. It is shown that using the SIMPLE algorithm in a multigrid context dramatically improves convergence, although the multigrid convergence rates are much worse than for Newtonian flows. The numerical results obtained for Bingham numbers as high as 1000 compare favourably with reported results of other methods.