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.     

Flow of Robertson-Stiff fluids through an eccentric annulus

I.IntroductionInthepetroleumindustry,itisusuallythecasethatthedrillingstring(orcasingpipe)isnotlocatedinsidethecenteroftheflowgeometryduringdrillingandcompletillgpl'ocess,this,inturll,willaltertile'11owingbehaviorofdrillingmudandcementslurl.ieswhichtlowinginaneccentricallllulus.OilaccoUlltofviscousofnon-Newtoniantluid,flowillaneccentricannulusdifTeresmuchwiththatinacollcentricannulus.Mailyinvestigatorshaveconductedresearchworksonthisstlbject.Earlyin1935,TaoandDollovallrealizedthattheimpel.t…     

Flow instabilities during annular displacement of one non-Newtonian fluid by another

This paper describes an experimental setup for axial laminar flow of liquids in the annulus between two eccentered cylinders. The design uses a conductivity method for measuring peak axial velocities around the annulus, and for the determination of displacement efficiency when displacing one fluid by another (displacement efficiency being defined as the ratio of volume of displaced fluid removed from the annulus, to the volume of the annulus, after a given number of annular volumes have been pumped). In an eccentric annulus, lower axial velocity in the narrow side produces “channeling” of the displacing fluid in the wide side and reduces the displacement efficiency. A positive density contrast between the two fluids can increase the efficiency by promoting azimuthal flow of the (denser) displacing fluid towards the narrow side. In this paper we report that gravity driven azimuthal flow is prone to severe instabilities which accelerate the displacement process but may leave behind an immobile strip of the displaced fluid in the narrow side.     

Numerical study of the combined effects of inertia,slip, and compressibility in extrusion of yield stress fluids

The axisymmetric extrudate swell flow of a compressible Herschel–Bulkley fluid with wall slip is solved numerically. The Papanastasiou-regularized version of the constitutive equation is employed, together with a linear equation of state relating the density of the fluid to the pressure. Wall slip is assumed to obey Navier’s slip law. The combined effects of yield stress, inertia, slip, and compressibility on the extrudate shape and the extrudate swell ratio are analyzed for representative values of the power-law exponent. When the Reynolds number is zero or low, swelling is reduced with the yield stress and eventually the extrudate contracts so that the extrudate swell ratio reaches a minimum beyond which it starts increasing asymptotically to unity. Slip suppresses both swelling and contraction in this regime. For moderate Reynolds numbers, the extrudate may exhibit necking and the extrudate swell ratio initially increases with yield stress reaching a maximum; then, it decreases till a minimum corresponding to contraction, and finally, it converges asymptotically to unity. In this regime, slip tends to eliminate necking and may initially cause further swelling of the extrudate, which is suppressed if slip becomes stronger. Compressibility was found to slightly increase swelling, this effect being more pronounced for moderate yield stress values and wall slip.     

Maxwell流体环管不定常流动解析解

本文对Maxwell非牛顿流体在环管内的不定常流动进行分析并得到了解析结果。依此,详尽分析了速度剖面、剪切应力随时间、非牛顿系数和环管形状变化的特征,为管道工程设计提供了必要的理论依据。     

Nonhomogeneous Incompressible Herschel–Bulkley Fluid Flows Between Two Eccentric Cylinders

The equations for the nonhomogeneous incompressible Herschel–Bulkley fluid are considered and existence of a weak solution is proved for a boundary-value problem which describes three-dimensional flows between two eccentric cylinders when in each two-dimensional cross-section annulus the flow characteristics are the same. The rheology of such a fluid is defined by a yield stress τ* and a discontinuous stress-strain law. A fluid volume stiffens if its local stresses do not exceed τ*, and a fluid behaves like a nonlinear fluid otherwise. The flow equations are formulated in the stress–velocity–density–pressure setting. Our approach is different from that of Duvaut–Lions developed for the classical Bingham viscoplastic fluids. We do not apply the variational inequality but make use of an approximation of the generalized Bingham fluid by a non-Newtonian fluid with a continuous constitutive law.     

Influence of inertia,topography and gravity on transient axisymmetric thin‐film flow

This study examines theoretically the development of early transients for axisymmetric flow of a thin film over a stationary cylindrical substrate of arbitrary shape. The fluid is assumed to emerge from an annular tube as it is driven by a pressure gradient maintained inside the annulus, and/or by gravity in the axial direction. The interplay between inertia, annulus aspect ratio, substrate topography and gravity is particularly emphasized. Initial conditions are found to have a drastic effect on the ensuing flow. The flow is governed by the thin‐film equations of the ‘boundary‐layer’ type, which are solved by expanding the flow field in terms of orthonormal modes in the radial direction. The formulation is validated upon comparison with the similarity solution of Watson (J. Fluid Mech 1964; 20 :481) leading to an excellent agreement when only 2–3 modes are included. The wave and flow structure are examined for high and low inertia. It is found that low‐inertia fluids tend to accumulate near the annulus exit, exhibiting a standing wave that grows with time. This behaviour clearly illustrates the difficulty faced with coating high‐viscosity fluids. The annulus aspect is found to be influential only when inertia is significant; there is less flow resistance for a film over a cylinder of smaller diameter. For high inertia, the free surface evolves similarly to two‐dimensional flow. The substrate topography is found to have a significant effect on transient behaviour, but this effect depends strongly on inertia. It is observed that the flow of a high‐inertia fluid over a step‐down exhibits the formation of a secondary wave that moves upstream of the primary wave. Gravity is found to help the film (coating) flow by halting or prohibiting the wave growth. The initial film profile and velocity distribution dictate whether the fluid will flow downstream or accumulate near the annulus exit. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Model-free inversion of squeeze-flow rheometer data

A method for extracting rheological data from squeeze-flow tests is proposed. The analysis, based upon the lubrication approximation for a generalized Newtonian fluid, differentiates experimental data in order to obtain an estimate of the wall shear rate (as in the Weissenberg–Rabinowitsch correction for the capillary rheometer) and of the wall shear stress. Two examples are discussed. The first is based on an approximate expression for the force required to squeeze a Herschel–Bulkley fluid. The second example concerns a power-law fluid with partial slip at the plates (but non-zero wall shear stress). Second derivatives of the experimental data are required: the interpretation of noisy results is therefore likely to be difficult.     

Estimation of the parameters of Herschel-Bulkley fluid under wall slip using a combination of capillary and squeeze flow viscometers

The determination of the parameters of viscoplastic fluids subject to wall slip is a special challenge and accurate results are generally obtained only when a number of viscometers are utilized concomitantly. Here the characterization of the parameters of the Herschel-Bulkley fluid and its non-linear wall slip behavior is formulated as an inverse problem which utilizes the data emanating from capillary and squeeze flow rheometers. A finite element method of the squeeze flow problem is employed in conjunction with the analytical solution of the capillary data collected following Mooneys procedure, which uses dies with differing surface to volume ratios. The uniqueness of the solution is recognized as a major problem which limits the accuracy of the solution, suggesting that the search methodology should be carefully selected.     

Pressure-driven and electroosmotic non-Newtonian flows through microporous media via lattice Boltzmann method

The lattice Boltzmann method is developed to simulate the pressure-driven flow and electroosmotic flow of non-Newtonian fluids in porous media based on the representative elementary volume scale. The flow through porous media was simulated by including the porosity into the equilibrium distribution function and adding a non-Newtonian force term to the evolution equation. The non-Newtonian behavior is considered based on the Herschel–Bulkley model. The velocity results for pressure-driven non-Newtonian flow agree well with the analytical solutions. For the electroosmotic flow, the influences of porosity, solid particle diameter, power law exponent, yield stress and electric parameters are investigated. The results demonstrate that the present lattice Boltzmann model is capable of modeling non-Newtonian flow through porous media.     

Exact solutions for generalized Burgers’ fluid in an annular pipe

This paper deals with some unsteady unidirectional transient flows of generalized Burgers’ fluid in an annular pipe. Exact solutions of some unsteady flows of generalized Burgers’ fluid in an annular pipe are obtained by using Hankel transform and Laplace transform. The following two problems have been studied: (1) Poiseuille flow due to a constant pressure gradient; (2) axial Couette flow in a annulus. The well known solutions for Navier-Stokes fluid, as well as those corresponding to a Maxwell fluid, a second grade fluid and an Oldroyd-B fluid appear as limiting cases of our solutions.     

Conjugate heat transfer within a concentric annulus filled with micropolar fluid

This paper investigates numerically the conjugate heat transfer in an annulus between two concentric cylinders. The annulus contains micropolar fluid and is heated isothermally from its inner wall. The effect of Rayleigh number, thickness of inner wall, inner wall-fluid thermal conductivity ratio, and material parameters of micropolar fluid on heat transfer rate within the annulus has been investigated. The study has shown that for low Rayleigh number regimes and for thermal conductivity of the inner wall greater than that of the fluid, the increase of inner wall thickness increases the heat transfer rate through the annulus and vice versa. While for convection dominating regimes Ra ≥ 10⁴ the increase of inner wall thickness decreases the heat transfer rate. Moreover, the study has shown that for fixed geometrical and flow parameters the heat transfer decreases in case of micropolar fluids in comparison with that of Newtonian fluids.     

Stokes flow of micro-polar fluids by peristaltic pumping through tube with slip boundary condition

This paper studies the Stokes flow of micro-polar fluids by peristaltic pumping through the cylindrical tube under the effect of the slip boundary condition. The motion of the wall is governed by the sinusoidal wave equation. The analytical and numerical solutions for the axial velocity, the micro-polar vector, the stream function, the pressure gradient, the friction force, and the mechanical efficiency are obtained by using the lubrication theory under the low Reynolds number and long wavelength approximations. The impacts of the emerging parameters, such as the coupling number, the micro-polar parameter, the slip parameter on pumping characteristics, the friction force, the velocity profile, the mechanical efficiency, and the trapping phenomenon are depicted graphically. The numerical results infer that large pressure is required for peristaltic pumping when the coupling number is large, while opposite behaviors are found for the micro-polar parameter and the slip parameter. The size of the trapped bolus reduces with the increase in the coupling number and the micro-polar parameter, whereas it blows up with the increase in the slip parameter.     

Rheology of pulp suspensions using ultrasonic Doppler velocimetry

Conventional rheometry coupled with local velocity measurements (ultrasonic Doppler velocimetry) are used to study the flow behaviour of various commercial pulp fibre suspensions at fibre mass concentrations ranging from 1 to 5 wt.%. Experimental data obtained using a stress-controlled rheometer by implementing a vane in large cup geometry exhibits apparent yield stress values which are lower than those predicted before mainly due to existence of apparent slip. Pulp suspensions exhibit shear-thinning behaviour up to a high shear rate value after which Newtonian behaviour prevails. Local velocity measurements prove the existence of significant wall slippage at the vane surface. The velocimetry technique is also used to study the influence of pH and lignin content on the flow behaviour of pulp suspensions. The Herschel–Bulkley constitutive equation is used to fit the local steady-state velocity profiles and to predict the steady-state flow curves obtained by conventional rheometry. Consistency between the various sets of data is found for all suspensions studied, including apparent yield stress, apparent wall slip and complete flow curves.     

环形截面螺旋管道内二次流动特性的研究

从曲线柱坐标系下的N-S方程出发,以曲率和挠率为小参数,采用摄动法求解了环形截面螺旋管道内的黏性流动,给出了完全二阶摄动解,结果表明：当挠率为零时,二次流表现为上下对称的四个涡;当挠率不为零,涡的对称性遭到破坏,二次涡的强度和个数受De数和环形截面内外径之比δ的影响,轴向速度最大值在De数较小时靠近管道的内侧,随着De数的增加,其最大值向外侧移动。     

Transport of reactive solutes in unsteady annular flow subject to wall reactions

The present paper concerns with the dispersion process in steady and oscillatory flows through an annular pipe in presence of reversible and irreversible reactions at the wall. Method of homogenization, a multiple-scale method of averaging, is adopted for deriving the effective transport equations. The main objective is to look into the effect of aspect ratio of the annular pipe on the dispersion coefficient due to the combined effect of axial convection and radial diffusion in steady and oscillatory flows along the annulus, subject to the kinetic reversible phase exchange and irreversible absorption at the outer wall. Results demonstrate that upto a certain critical value of aspect ratio, dispersion coefficient increases with increase of aspect ratio when the wall is retentive, though the wall inertness may lead to decrease of dispersion coefficient with increase of aspect ratio. The results would be useful to the medical practitioners working in the domain of catheterized artery.     

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.     

Using the Euler–Lagrange variational principle to obtain flow relations for generalized Newtonian fluids

The Euler–Lagrange variational principle is used to obtain analytical and numerical flow relations in cylindrical tubes. The method is based on minimizing the total stress in theflow duct using the fluid constitutive relation between stress and rate of strain. Newtonian and non-Newtonian fluid models, which include power law, Bingham, Herschel–Bulkley, Carreau, and Cross, are used for demonstration.     

Peristaltic motion of third grade fluid in curved channel

＂ Analysis is performed to study the slip effects on the peristaltic flow of non-Newtonian fluid in a curved channel with wall properties. The resulting nonlinear partial differential equations are transformed to a single ordinary differential equation in a stream function by using the assumptions of long wavelength and low Reynolds number. This differential equation is solved numerically by employing the built-in routine for solving nonlinear boundary value problems （BVPs） through the software Mathematica. In addition, the analytic solutions for small Deborah number are computed with a regular perturbation technique. It is noticed that the symmetry of bolus is destroyed in a curved channel. An intensification in the slip effect results in a larger magnitude of axial velocity. Further, the size and circulation of the trapped boluses increase with an increase in the slip parameter. Different from the case of planar channel, the axial velocity profiles are tilted towards the lower part of the channel. A comparative study between analytic and numerical solutions shows excellent agreement.