Exact solution for compressive flow of viscoplastic fluids under perfect slip wall boundary conditions

Based on the perfect slip condition between rigid walls and fluids, the compressive flow of Herschel-Bulkley fluids and biviscous fluids was studied. The explicit expressions of stresses and fluid velocity were given. To move the rigid walls for a Herschel-Bulkley fluid with the yield stress (τ₀), the mean pressure applied onto the rigid wall should be larger than 2τ₀/. No yield surface exists in the interior of the fluids when flow occurs. For a biviscous fluid, a critical load was given. The fluid behaves like the Bingham fluid when the external applied load onto the wall is larger than the critical load, otherwise the fluid is Newtonian. Received: 10 June 1997 Accepted: 22 September 1997     

Flow of shear-thinning fluids in a concentric annulus

Distributions of mean axial velocity, axial and tangential turbulence intensities together with friction factor versus Reynolds number (f-Re) data are presented for three non-Newtonian liquids in fully developed laminar, transitional and turbulent flow in an annular geometry in the absence of centrebody rotation. Each of the non-Newtonian fluids was shear thinning and to some extent elastic and one was also thixotropic in character. For comparison purposes, measurements are also reported for a Newtonian fluid.In the case of the Newtonian fluid, a mixture of glucose syrup and water, the f-Re data in both laminar and turbulent flow follow the appropriate relationships for the annular geometry, with a clear demarcation at transition which is confirmed independently by a measured increase in the centre-channel axial turbulence intensities. The measured velocity profiles for laminar flow are in good agreement with those predicted theoretically, whilst the turbulent profiles obey the log-law relationship over much of the mid-channel region and tend to the u ⁺=y ⁺ relationship in the immediate vicinity of both walls.For the first non-Newtonian fluid, an aqueous solution of sodium carboxymethylcellulose (CMC), good agreement with theoretical predictions for a power-law fluid was observed in the f-Re data in the laminar regime with evidence of drag reduction in turbulent flow. Velocity profiles, determined in two planes, indicate minor circumferential asymmetry in laminar flow. Law-of-the-wall plots for fully turbulent flow indicate an upward shift in the data in the log-law region of the annulus consistent with the drag-reduction behaviour, as also observed in pipe-flow experiments for this fluid (Escudier et al. 1992). In the near-surface regions of both the outer and inner tubes the data again tend towards the u ⁺=y ⁺ relationship.Anomalous behaviour was observed in the f-Re curves for the second non-Newtonian fluid, 0.125% and 0.2% aqueous solutions of Xanthan gum, with data for both concentrations falling significantly below the appropriate f-Re relationship for a power-law fluid. The anomalies are attributed to the elastic character of Xanthan gum. In the near-surface region of the outer tube the velocity-profile data again tend towards the u ⁺=y ⁺ relationship but it proved impossible to obtain data in the near vicinity of the inner wall due to slight turbidity of the fluid.The third non-Newtonian fluid, a Laponite/CMC blend, again exhibits anomalous f-Re behaviour, attributed to the thixotropic nature of this fluid. Velocity profiles determined in two planes again indicate some circumferential asymmetry in the laminar regime. Law-of-the-wall plots for the transitional and turbulent profiles tend towards the u ⁺=y ⁺ relationship in both near-wall regions, again with an upward shift in the core of the annulus, consistent with drag reduction.In general terms, the experimental results are consistent with previous work for non-Newtonian fluid flow in circular pipes and with limited data for an annular geometry (Nouri et al. 1993), with regard to drag reduction, modified turbulence structure and scale effects.List of Symbols D _i centrebody diameter (m) - D _o outer pipe diameter (m) - (D _o -D _i ) hydraulic diameter (m) - f friction factor (2 · _A /U²) - n power-law exponent (-) - p fluid static pressure (Pa) - Q volumetric flow rate (m³/s) - r radial distance from pipe centreline (m) - R _i centrebody radius (m) - R _o outer pipe radius (m) - R _n refractive index (-) - Re Reynolds number U(D _o -D _i )/ _s - s geometric scaling factor (-) - u mean axial velocity (m/s) - u rms fluctuation in axial velocity (m/s) - u _c1 rms fluctuation in centreline axial velocity (m/s) - u non-dimensional value of u (u/u ) - u friction velocity (m/s) - w rms fluctuation in tangential velocity (m/s) - x axial distance along pipe (m) - y distance from pipe or centrebody wall (m) - y ⁺ non-dimensional value of y (u y/v _s ) - p/L pressure drop per unit length (N/m²/m) - shear rate (s^-1) - radius ratio R _i /R _o - _C constant in Cross model (s) - _CA constant in Carreau model (s) - _HB constant in Herschel-Bulkley model (s) - _n constant in power-law model (s) - _S constant in Sisko model (s) - dynamic viscosity (Pa · s) - _ref reference viscosity (1 Pa · s) - _s viscosity at wall at prevailing surface shear stress (Pa · s) - ₀ zero shear-rate viscosity (Pa · s) - infinite shear-rate viscosity (Pa · s) - v kinematic viscosity (/) (m²/s) - v _s kinematic viscosity at wall (m²/s) - non-dimensional radial location (R _o -r)/(R _o -R _i ) - fluid density (kg/m³) - shear stress (Pa) - _A weighted average wall shear stress (Pa) - _i shear stress on centrebody (Pa) - _o shear stress on outer wall (Pa) - _s surface shear stress (Pa) - _ref reference shear stress (1 Pa) - _y fluid yield stress (Pa) - ^* geometry function of Jones and Leung (1981) The work reported here represents part of programme of research which has received financial support from SERC (GR/F 87813), BP Exploration Company Ltd, Shell Research BV and AEA Petroleum Services. This support is gratefully acknowledged. Frequent meetings with Professor J. H. Whitelaw, Imperial College of Science, Technology and Medicine, Dr. C. F. Lockyear and Dr. D. Ryan, BP Research, Ms. B. Kampman, Shell Research BV, and Dr. W. J. Worraker, AEA Technology, were of considerable benefit to the research.     

LES of turbulent flow in a concentric annulus with rotating outer wall

Fully-developed turbulent flow in a concentric annulus, r₁/r₂ = 0.5, Re_h = 12,500, with the outer wall rotating at a range of rotation rates N = U_θ,wall/U_b from 0.5 up to 4 is studied by large-eddy simulations. The focus is on the effects of moderate to very high rotation rates on the mean flow, turbulence statistics and eddy structure. For N up to ∼2, an increase in the rotation rate dampens progressively the turbulence near the rotating outer wall, while affecting only mildly the inner-wall region. At higher rotation rates this trend is reversed: for N = 2.8 close to the inner wall turbulence is dramatically reduced while the outer wall region remains turbulent with discernible helical vortices as the dominant turbulent structure. The turbulence parameters and eddy structures differ significantly for N = 2 and 2.8. This switch is attributed to the centrifuged turbulence (generated near the inner wall) prevailing over the axial inertial force as well as over the counteracting laminarizing effects of the rotating outer wall. At still higher rotation, N = 4, the flow gets laminarized but with distinct spiralling vortices akin to the Taylor–Couette rolls found between the two counter-rotating cylinders without axial flow, which is the limiting case when N approaches to infinity. The ratio of the centrifugal to axial inertial forces, Ta/Re² ∝ N² (where Ta is the Taylor number) is considered as a possible criterion for defining the conditions for the above regime change.     

Analysis on creeping channel flows of compressible fluids subject to wall slip

Creeping channel flows of compressible fluids subject to wall slip are widely encountered in industries. This paper analyzes such flows driven by pressure in planar as well as circular channels. The analysis elucidates unsteady flows of Newtonian fluids subject to the Navier slip condition, followed by steady flows of viscoplastic fluids, in particular, Herschel–Bulkley fluids and their simplifications including power law and Newtonian fluids, that slip at wall with a constant coefficient or a coefficient inversely proportional to pressure. Under the lubrication assumption, analytical solutions are derived, validated, and discussed over a wide range of parameters. Analysis based on the derived solutions indicates that unsteadiness alters cross-section velocity profiles. It is demonstrated that compressibility of the fluids gives rise to a concave pressure distribution in the longitudinal direction, whereas wall slip with a slip coefficient that is inversely proportional to pressure leads to a convex pressure distribution. Energy dissipation resulting from slippage can be a significant portion in the total dissipation of such a flow. A distinctive feature of the flow is that, in case of the pressure-dependent slip coefficient, the slip velocity increases rapidly in the flow direction and the flow can evolve into a pure plug flow at the exit.     

Axial Couette–Poiseuille flow of Bingham fluids through concentric annuli

In this paper, the axial Couette–Poiseuille flow of Bingham fluids through concentric annuli is studied. Analytical solutions of different types of flow are derived. Compared to previous studies, we emphasize two new types of flow, which have been missed previously, are found in our results. Hence, there are eight different forms of the velocity profile depending on values of three dimensionless parameters, which are the Bingham, axial Couette numbers and the radius ratio. Distributions of these eight forms are specified in the parameter plane of axial Couette number vs. Bingham number for various radius ratios. These new flow regimes are analyzed from both a mathematical and physical perspective.     

Heat transfer to non-Newtonian fluids in laminar flow through concentric annuli with or without suction

Rheologica Acta - In this paper we have investigated the heat transfer aspects of the flow of a power law fluid in an annulus with porous walls. The case of no suction (solid walls annulus) is also...     

Purely tangential flow of a PTT-viscoelastic fluid within a concentric annulus

An analytical solution is presented for the steady state, purely tangential flow of a viscoelastic fluid obeying the Phan-Thien–Tanner (PTT) constitutive equation in a concentric annulus with relative rotation of the inner and outer cylinders. The influence on the velocity distribution within the annulus and on fRe of the Weissenberg number, aspect ratio and an elongational parameter are investigated. The results show that the differences between the radial location of the minimum velocity and of the critical angular velocity compared with their Newtonian counterparts increase as the fluid elasticity increases. The results also show that fRe decreases with increasing Weissenberg number, radius ratio and the elongational parameter in the case of inner-cylinder rotation. In contrast, fRe increases with increasing radius ratio when the outer cylinder is rotating while the inner cylinder is at rest.     

The flow of viscoplastic material between two concentric spheres

The motion of a viscoplastic medium between two concentric spheres is considered upon rotation of one sphere with constant angular velocity. This problem is solved by an heuristic iterative method. The boundary of the stagnation zones is found and its specific shape is shown. The flow characteristics versus the parameter of the medium are obtained. Voronezh State Engineering Academy, Voronezh 394017. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 1, pp. 133–139, January–February, 1999.     

Fountain flow of pseudoplastic and viscoplastic fluids

Numerical simulations have been undertaken for the benchmark problem of fountain flow present in injection-mold filling. The Finite Element Method (FEM) is used to provide numerical results for both cases of planar and axisymmetric domains under steady-state conditions. The Herschel–Bulkley model of viscoplasticity is used, which reduces with appropriate modifications to the Bingham, power-law and Newtonian models. The present results extend previous ones regarding the shape of the front, which is essential in correctly capturing the flow field. In particular the centreline front position is found as a function of the dimensionless power-law index (in the case of pseudoplasticity) and the dimensionless yield stress (in the case of viscoplasticity). The pressures from the simulations have been used to compute the excess pressure losses in the system (front pressure correction or exit correction). Both shear-thinning and shear-thickening lead to more extended front positions relative to the Newtonian values, which are 0.895 for the planar case and 0.835 for the axisymmetric one. Viscoplasticity leads also to more extended front positions as the dimensionless yield stress goes from zero (Newtonian behaviour) to higher values of the yield stress. In both cases of non-Newtonian behaviour, the front tends to follow the development of the fully developed Poiseuille velocity profile, which tends towards a plug-like profile at the extreme cases of non-Newtonianness. The front pressure (exit) correction increases monotonically with the decrease in the power-law index and the increase in the dimensionless yield stress.     

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.     

Unsteady flow in an annulus between two concentric rotating spheres

An analysis is presented for the unsteady laminar flow of an incompressible Newtonian fluid in an annulus between two concentric spheres rotating about a common axis of symmetry. A solution of the Navier-Stokes equations is obtained by employing an iterative technique. The solution is valid for small values of Reynolds numbers and acceleration parameters of the spheres. In applying the results of this analysis to a rotationally accelerating sphere, a virtual moment of inertia is introduced to account for the local inertia of the fluid.     

Unsteady flow in an annulus between two concentric rotating spheres

An analysis is presented for the unsteady laminar flow of an incompressible Newtonian fluid in an annulus between two concentric spheres rotating about a common axis of symmetry. A solution of the Navier-Stokes equations is obtained by employing an iterative technique. The solution is valid for small values of Reynolds numbers and acceleration parameters of the spheres. In applying the results of this analysis to a rotationally accelerating sphere, a virtual moment of intertia is introduced to account for the local inertia of the fluid.Nomenclature R _i radius of the inner sphere - R _o radius of the outer sphere - radial coordinate - r dimensionless radial coordinate, - meridional coordinate - azimuthal coordinate - time - t dimensionless time, - Re _i instantaneous Reynolds number of the inner sphere, _i R _k ² / - Re _o instantaneous Reynolds number of the outer sphere, _o R _o ² / - radial velocity component - V _r dimensionless radial velocity component, - meridional velocity component - V dimensionless meridional velocity component, - azimuthal velocity component - V dimensionless azimuthal velocity component, - viscous torque - T dimensionless viscous torque, - viscous torque at surface of inner sphere - T _i dimensionless viscous torque at surface of inner sphere, - viscous torque at surface of outer sphere - T _o dimensionless viscous torque at surface of outer sphere, - externally applied torque on inner sphere - T _p,i dimensionless applied torque on inner sphere, - moment of inertia of inner sphere - Z _i dimensionless moment of inertia of inner sphere, - virtual moment of inertia of inner sphere - Z _i,v dimensionless virtual moment of inertia of inner sphere, - virtual moment of inertia of outer sphere - _i instantaneous angular velocity of the inner sphere - _o instantaneous angular velocity of the outer sphere - density of fluid - viscosity of fluid - kinematic viscosity of fluid,/ - radius ratio,R _i/R _o - swirl function, - dimensionless swirl function, - stream function - dimensionless stream function, - _i acceleration parameter for the inner sphere, - _o acceleration parameter for the outer sphere, - shear stress - _r dimensionless shear stress,      

Entry flow behaviour of viscoelastic fluids in an annulus

Developing and fully developed velocity profiles in the entrance region of an abrupt 2-to-1 annular contraction were measured for a number of visco-elastic polymer solutions. Experimental results were obtained for Reynolds number and flow behaviour index in the range 9.8 ? Re ? 355 and 0.372 ? n ? 0.55 respectively. A momentum-energy integral technique was employed in the boundary layer analysis. The deviation from inelastic behaviour was indicated by the ratio of elastic to inertial forces, Ws/Re. Within the limits of confidence of the experimental results, good agreement with theoretical predictions was obtained and very little deviation from inelastic behaviour was observed for Ws/Re < 0.08. For the test fluids investigated, the entrance length was found to be longer than that predicted for the corresponding inelastic fluids of the same n.     

Axial Couette flow of an Oldroyd-B fluid in an annulus

This paper establishes the velocity field and the adequate shear stress corresponding to the motion of an Oldroyd-B fluid between two infinite coaxial circular cylinders by means of finite Hankel transforms. The flow of the fluid is produced by the inner cylinder which applies a time-dependent longitudinal shear stress to the fluid. The exact analytical solutions, presented in series form in terms of Bessel functions, satisfy all imposed initial and boundary conditions. The general solutions can be easily specialized to give similar solutions for Maxwell, second grade and Newtonian fluids performing the same motion. Finally, some characteristics of the motion as well as the influence of the material parameters on the behavior of the fluid motion are graphically illustrated.     

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.     

Heat transfer to non-Newtonian fluids in laminar flow through rectangular ducts

Heat transfer to non-newtonian fluids flowing laminarly through rectangular ducts is examined. The conservation equations of mass, momentum, and energy are solved numerically with the aid of a finite volume technique. The viscoelastic behavior of the fluid is represented by the Criminale-Ericksen-Filbey (CEF) constitutive equation. Secondary flows occur due to the elastic behavior of the fluid, and, consequently, heat transfer is strongly enhanced. It is observed that shear thinning yields negligible heat transfer enhancement effect, when compared with the secondary flow effect. Maximum heat transfer is shown to occur for some combinations of parameters. Thus, there are optimal combinations of aspect ratio and Reynolds numbers, which depend on the fluid's mechanical behavior. This result can be usefully explored in thermal designs of certain industrial processes.     

Heat transfer and pressure drop of a transversely finned concentric annulus with longitudinal flow

The characteristics of heat transfer and pressure drop have been experimentally studied for the fully developed concentric annular flow with transverse fins normal to the flow direction by the naphthalene sublimation technique. Correlations for calculating the heat transfer coefficient with different inner diametersD ₀ of the outer tube are presented. A characteristic Reynolds number has been proposed, by which the predominant role of the transverse fins can be evaluated. It has been indicated that the inner diameterD ₀ has much more effect on pressure drop than on heat transfer. The effect ofD ₀ on the overall performance is also compared under the same flow velocity or flow rate. It has been found that the effect of developing flow on heat transfer is significant and should be taken into account during experiment.

---

Wärmeübergang und Druckverlust in einem querberippten konzentrischen Ringkanal bei Längsströmung

Zusammenfassung Mit Hilfe der Naphthalin-Sublimationstechnik werden die Wärmeübergangs- und Druckverlustcharakteristiken in quer zur voll ausgebildeten Strömung bei berippten Ringkanälen experimentell ermittelt und Korrelationen zur Berechnung des Wärmeübergangskoeffizienten bei variablen InnendurchmesserD ₀ des umschließenden Rohres angegeben. Ferner wird eine charakteristische Reynolds-Zahl vorgeschlagen, über die sich der dominierende Einfluß der Querrippen erfassen läßt. Es zeigte sich, daß der InnendurchmesserD ₀ den Druckverlust wesentlich mehr beeinflußt als den Wärmeübergang. Auch wurde die Abhängigkeit des Gesamt-Übertragungsverhaltens vonD ₀ bei gleicher Strömungsgeschwindigkeit bzw. Volumenstromdichte ermittelt. Es zeigte sich, daß die Ausbildung des Strömungsprofils bei Einlaufströmung den Wärmeübergang signifikant beeinflußt und deshalb im Experiment zu berücksichtigen ist.

     

Simultaneously developing laminar flow in an isothermal micro-tube with slip flow models

In this study, a two-dimensional steady state simultaneously developing laminar flow inside a micro-tube is investigated numerically under slip flow conditions. The first and second-order slip flow models have been implemented for the case where the viscous dissipation and axial conduction are included and a constant wall temperature boundary condition is specified. The results are obtained for several combinations of the Knudsen number Kn, the coefficient β and the Brinkman number Br. The study reveals a significant impact of slip flow and temperature jump on the hydrodynamic and thermal fields. A comparison of the first and second-order slip flow shows a considerable variation of the hydrodynamic flow and a weak impact on the thermal field particularly when the flow is fully developed.     

Instability of a slip flow in a curved channel formed by two concentric cylindrical surfaces

Instability of a slip flow in a curved channel formed by two concentric cylindrical surfaces is investigated. Two cases are considered. In the first (Taylor–Couette flow) case the flow is driven by the rotation of the inner cylindrical surface; no azimuthal pressure gradient is applied. In the second case (Dean flow) both cylindrical surfaces are motionless, and the flow is driven by a constant azimuthal pressure gradient. The collocation method is used to find numerically the critical values of the Taylor and Dean numbers, which establish the instability criteria for these two cases. The dependencies of critical values of these numbers on the ratio between the radii of concave and convex walls and on the velocity slip coefficient are investigated.