Steady Flow of Power Law Fluids across a Circular Cylinder

The momentum equations describing the steady cross‐flow of power law fluids past an unconfined circular cylinder have been solved numerically using a semi‐implicit finite volume method. The numerical results highlighting the roles of Reynolds number and power law index on the global and detailed flow characteristics have been presented over wide ranges of conditions as 5 ≤ Re ≤ 40 and 0.6 ≤ n ≤ 2. The shear‐thinning behaviour (n < 1) of the fluid decreases the size of recirculation zone and also delays the separation; on the other hand, the shear‐thickening fluids (n > 1) show the opposite behaviour. Furthermore, while the wake size shows non‐monotonous variation with the power law index, but it does not seem to influence the values of drag coefficient. The stagnation pressure coefficient and drag coefficient also show a complex dependence on the power law index and Reynolds number. In addition, the pressure coefficient, vorticity and viscosity distributions on the surface of the cylinder have also been presented to gain further physical insights into the detailed flow kinematics.     

Drag on chains and agglomerates of spheres in viscous newtonian and power law fluids

New experimental data on the free settling velocity of straight chains (up to twenty spheres) and planar clusters of touching spheres in Newionian and power law media are reported. The results embrace the following ranges of conditions: 0.65 ≤ n ≤ 1; Re < - 2.5 and 1.22 < m < 48.87 Pa·sⁿ. The straight chain drag measurements are in line with theoretical predictions for Newtonian fluids. The present results in power law fluids seem to suggest that it is possible to express the drag on a straight chain of spheres in terms of that on a single sphere of equal volume. Limited results with planar clusters are satisfactorily correlated using a volume equivalent sphere diameter.     

Sedimentation of a sphere along the axis of a long square duet filled with non-newtonian liquids

Based on extensive experimental results, it is shown that the retardation effect caused by the confining walls on the free settling velocity of a sphere is smaller with square walls than that with cylindrical boundaries. This is true for both Newtonian and power law fluids, provided the particle Reynolds number is small (< about 5). The values of the wall factor for Newtonian liquids are in excellent agreement with theory (up to R / L ≤ 0.1) while those for power law fluids have been correlated empirically via a linear relationship. The results reported here encompass the following ranges of conditions: 1 ≥ n ≥ 0.7; Re < 15 and 0.024 < R/L < 0.238.     

Settling velocity of cubes in Newtonian and power law liquids

New extensive data on the free settling velocity of thirty cubes of various densities and sizes falling in scores of Newtonian and Power law liquids are reported herein to supplement the existing data, for there is very little prior data on cubes in power law liquids. The new data embrace the range of conditions as follows: sphericity of 0.805; power law index, 0.61 to 1 and consistency index, 0.0078-15.31 Pa sⁿ; Reynolds number, 0.0013 to 860. The new results are shown to be consistent with an existing drag correlation which has been tested extensively using the literature data for spherical and non-spherical particles falling in Newtonian and power law liquids with acceptable levels of accuracy.     

Steady Two‐Dimensional Non‐Newtonian Flow Past an Array of Long Circular Cylinders up to Reynolds number 500: A Numerical Study

The free surface and zero vorticity cell models have been combined with the equations of motion to investigate numerically the steady flow of incompressible power‐law (shear‐thinning and shear‐thickening) fluids across banks of long cylinders. The equations of motion in the stream function/vorticity formulation have been solved numerically using a second order accurate finite difference method to obtain extensive information on the behaviour of the drag coefficient, surface vorticity distribution, streamlines and iso‐vorticity patterns, for high Reynolds numbers (Re = 50 500) and using a wide range of power‐law index (0.3 ≤ n ≤ 2.0), and porosity (0.4 ≤ e ≤0.9) values. The behaviour of the aforementioned parameters at low Reynolds numbers has also been investigated and validated using theoretical and numerical work from the literature. The results reported here enable extension of the limits of creeping flow behaviour up to Re = 50 for fluids with highly shear‐thickening characteristics under low porosity conditions.     

A numerical study of the accelerating motion of a dense rigid sphere in non-newtonian power law fluids

The equations of motion of an accelerating sphere falling through non-Newtonian fluids with power law index n in the range 0.2 ≤ n ≤ 1.8 were integrated numerically using the assumption that the drag on the sphere was a function of both power law index and terminal Reynolds number, Re_t For 10^?2 ≤ Re_t ≤ 10³ both dimensionless time and distance travelled by the sphere under transient conditions showed a much stronger dependence on the flow behaviour index, n, for shear-thinning than for shear-thickening fluids. The form of this dependence is investigated here. Furthermore, results in four typical shear-thinning fluids suggested a strong correlation between the distance and time travelled by the sphere under transient conditions and the value of the fluid consistency index. The analysis reported herein is, however, restricted to dense spheres falling in less dense fluids, when additional effects arising from the Basset forces can be neelected.     

Mass transfer from a contaminated fluid sphere

The unsteady mass transfer from a contaminated fluid sphere moving in an unbounded fluid is examined numerically for unsteady‐state transfer. The effect of the interface contamination and the flow regime on the concentration profiles, inside and outside a fluid sphere, is investigated for different ranges of Reynolds number (0 < Re < 200) and Peclet number (0 < Pe < 10⁵), viscosity ratio between the dispersed phase and the continuous phase (0 < κ < 10), and the stagnant‐cap angle (0° < θ_cap < 180°). It was found that the stagnant‐cap angle significantly influences the mass transfer from the sphere to a surrounding medium. For all Peclet and Reynolds numbers and κ, the contamination reduces the mass transfer flux. The average Sherwood number increases with an increase of stagnant‐cap angle and reaches a maximum equal to the average one for a clean fluid sphere at low viscosity ratio and large Peclet numbers. A predictive equation for the Sherwood number is derived from these numerical results. © 2010 American Institute of Chemical Engineers AIChE J, 2011     

Drag on Freely Falling Cones in Newtonian and in Power Law Fluids

New extensive data on the terminal falling velocities of conical shaped bodies in scores of Newtonian and power law fluids are reported. Altogether, 11 Newtonian and 11 non‐Newtonian test liquids together with 33 cones made from four different materials and 14 spheres of three different materials have been used to gather 486 individual data points covering wide ranges of conditions as follows: Reynolds number 0.0019 to 507; power law flow behaviour index 0.4 to 1, the value of sphericity 0.59 to 0.79; and the cone‐to‐fall tube diameter ratios up to 0.264 to assess the extent of wall effects. A simple expression is developed to estimate the terminal falling velocity of a cone from a knowledge of its dimensions, and the terminal velocity of an equivalent sphere. A generalized drag equation applicable to both Newtonian and power law liquids is also presented.     

Particle trajectories and terminal velocities in vertically oscillating fluids

The trajectories and terminal velocities of particles in vertically oscillating fluids have been studied by obtaining analytic and numerical solutions to the nonlinear Langevin equation representing a superposition of forces arising from particle acceleration, displaced fluid acceleration, buoyancy and particle-fluid drag as represented by an n^th-power drag law. In vibrating fluids the directional particle velocities are found to be lower than in stationary fluids for drag exponents n > 1 and correspondingly higher for 0 < n < 1. The engineering significance of the results is discussed in relation to hold-up, separation and transport phenomena in vertically oscillating multiparticle-fluid systems.     

Drag on a single fluid sphere translating in power-law liquids at moderate Reynolds numbers

A numerical investigation has been carried out to obtain the steady state drag coefficients and flow patterns of a single Newtonian fluid sphere sedimenting in power-law liquids. A finite difference method based simplified marker and cell (SMAC) algorithm has been implemented on a staggered grid arrangement to solve the continuity and momentum equations. For both phases, the convective terms have been discretized using the quadratic upstream interpolation for convective kinematics (QUICK) scheme, and diffusive and non-Newtonian terms with central differencing scheme. An exponential transformation has been applied in the radial direction for the continuous phase computational domain. In order to ensure the accuracy of the solver, extensive validation has been carried out by comparing the present results with the existing literature values for a few limiting cases. Further, in this study the effects of the Reynolds number (Re_o), internal to external fluid characteristic viscosity ratio (k) and power-law index (n_o) on the continuous phase flow field, pressure drag (C_dp), friction drag (C_df) and total drag (C_D) coefficients have been analyzed over the range of parameters: 5?Re_o?500, 0.1?k?50 and 0.6?n_o?1.6. Based on numerical results obtained in this work, a simple correlation has been proposed for the total drag coefficient, which can be used to predict the rate of sedimentation of a fluid sphere in power-law liquids.     

Hydrodynamic drag forces on two porous spheres moving along their centerline

This paper numerically evaluates the hydrodynamic drag force exerted on two highly porous spheres moving steadily along their centerline (sphere #1 and sphere #2) through a quiescent Newtonian fluid over a Reynolds number ranging from 0.1 to 40. At creeping flow limit, the drag forces exerted on both spheres were identical. At higher Reynolds numbers the drag force on sphere #1 was higher than sphere #2, revealing the shading effects produced by sphere #1 on sphere #2. At dimensionless diameter (β, =d_f/2k^0.5, d_f and k are floc diameter and interior permeability, respectively) >20, the spheres can be regarded nonporous. At β<20, the drag forces dropped. At β<2, the drag forces approached “no-spheres” limit. An increased size ratio of two spheres (d_f1/d_f2) would increase the drag force on sphere #1 and reduce that on sphere #2. At increasing β for both spheres, the drag force on sphere #2 was increased because of the more difficult advective flow through its interior, and at the same time the drag was reduced owing to the stronger wake flow produced by the denser sphere #1. The competition between these two effects leads to complicated dependence of drag force on sphere #2 on β value. These effects were minimal when β became low. Two identical spheres could move steadily along their centerline. At higher Reynolds number, the two spheres would move closer because of the incorporation of inertia force. For spheres of different diameters, the sphere # 2 would move faster than sphere #1 regardless of their size ratio and β value. This occurrence yielded efficient coagulation when two porous spheres were moving in-line.     

Terminal velocity of porous spheres

Terminal velocity of porous spheres was experimentally measured for a Reynolds number range of 0.2 to 120 for a normalized sphere radius, β = R/R of 15.6 to 33, where R and k are the sphere radius and permeability, respectively. The drag coefficient for 15 < β < 33 was found to be C_D = 24Ω/Re [1 + 0.1315 Re^{(0.82 - 0.05w)}] for 0.1 < Re ≤ 7 and C_D = 24Ω/Re [1 + 0.0853 Re^{(1.093 - 0.105w)}] for 7 < Re < 120 with w = log₁₀Re where Re is the sphere Reynolds number and Ω=2β² [1 - (tanh β/β)] / 2β² + 3[1 - tanh β/β)] At high Reynolds numbers, it was found that the porous sphere terminal velocity was less affected by the container walls than for the case of an impermeable sphere. However, at very low Reynolds numbers, the wall effects were found to be similar for both the permeable and the impermeable spheres.     

The Effects of Power Law Fluid Rheology on Coil Heat Transfer for Agitated Vessel in the Transitional Flow Regime

A CFD model of heat transfer from power‐law fluids to helical cooling coils in the transitional flow regime of a baffled tank mixed with a pitched blade turbine was developed with Fluent^TM. The model captured local temperature and velocity gradients. Simulations were run, varying Re, Pr, K and n. The results indicate that a Sieder‐Tate type correlation, with the exponent on and the coefficient in front of the Reynolds number being a function of n, is recommended for estimating ho. Also, a new two coil bank design was found to be more efficient when 450 < Re < 650.     

Shear and skin friction on particles in power-law fluids agitated by flat-blade and fluid foil impellers

The shear rates that exert angular deformation on spherical particles have been measured. The particles are mimiced by a spherical probe. The probe has been immersed in various impeller-agitated power law fluids. The fluids are aqueous dispersions of polymers, e.g. CMC, xanthan gum and starch. The probe has been positioned in various points of a stirred vessel and at various angles. Angle-averaged shear rate distributions were produced. The distributions obtained are characteristic for the specific impeller flow patterns. The flow patterns have been identified by computational fluid dynamics (CFD). Two types of impellers representative for the flat and the fluid-foil blade design, i.e., a Rushton flat-blade turbine (RT) and a Narcissus impeller (NS) are studied. The effects of rheological properties and blade design on the ‘shear-rate-on-particles’ distribution are examined. The local shear field non-uniformity has been uncovered and compared in terms of the CFD-generated time-averaged velocity and deformation rate profiles. The ‘shear-rate-on-particles’ distribution apart from the impeller is found to follow qualitatively the time-averaged inner flow shear rate distribution. Referring to impeller speed 5-12.5 Hz, the dimensionless wall shear rate varied between 200 and 1000. In power law fluids, the shear rate on particles decreased up to 50%. The fluid-foil NS-generated shear field was found comparable to the shear field induced by conventional flat-blade turbines and appeared in cases less sensitive to polymer presence. The shear rate produced by the fluid-foil impeller in the highly shear-thinning model solution (n∼0.4) exceeded the flat-blade RT-imposed shear rate. The analysis has been extended to skin friction drag on particles. It is shown that, while exerting an undoubtedly greater angular deformation in water-like fluids, in polymer presence the conventional flat-blade turbine introduces a flow geometry that imposes particle drag that is close or in some cases even less than the one generated by the fluid-foil impeller. The fact implies a weak shape effect of radial turbines on shear-sensitive particles or particle dispersions in power law liquids.     

Settling of cylinders in power law liquids

New measurements on drag coefficient and wall effects on the free settling motion of cylinders in two shearthinning polymer solutions are reported. The ranges of conditions covered in this study are: 1 < Re_∞∞ < 40; 0.25 < (LID) < 2; 0.079 < dID < 0.4, and n = 0.65 and 0.74. Both wall correction factor and drag coefficient results are in line with the Newtonian behavior provided a modified Reynolds number is used to represent the results.     

Steady flow of power-law fluids across an unconfined elliptical cylinder

The momentum transfer characteristics of the power-law fluid flow past an unconfined elliptic cylinder is investigated numerically by solving continuity and momentum equations using FLUENT (version 6.2) in the two-dimensional steady cross-flow regime. The influence of the power-law index (0.2?n?1.8), Reynolds number (0.01?Re?40) and the aspect ratio of the elliptic cylinder (0.2?E?5) on the local and global flow characteristics has been studied. In addition, flow patterns showing streamline and vorticity profiles, and the pressure distribution on the surface of the cylinder have also been presented to provide further physical insights into the detailed flow kinematics. For shear-thinning (n<1) behaviour and the aspect ratio E>1, flow separation is somewhat delayed and the resulting wake is also shorter; on the other hand, for shear-thickening (n>1) fluid behaviour and for E<1, the opposite behaviour is obtained. The pressure coefficient and drag coefficient show a complex dependence on the Reynolds number and power-law index. The decrease in the degree of shear-thinning behaviour increases the drag coefficient, especially at low Reynolds numbers. While the aspect ratio of the cylinder exerts significant influence on the detailed flow characteristics, the total drag coefficient is only weakly dependent on the aspect ratio in shear-thickening fluids. The effect of the flow behaviour index, however, diminishes gradually with the increasing Reynolds number. The numerical results have also been presented in terms of closure relations for easy use in a new application.     

Hydrodynamic force between two hard spheres tangentially translating in a power-law fluid

The hydrodynamic interaction between two hard spheres tangentially translating in a power-law fluid is investigated. By considering the gap between the two spheres being sufficiently small such that the Reynolds’ lubrication theory applies, an analytical equation to the pressure in the gap is obtained using truncated Fourier series. To a good approximation, the pressure equation can be further simplified. The simplified approximate equation over-predicts the pressure for shear thickening fluid (n>1) but under-predicts the pressure for shear-thinning fluid (n<1). However, the errors in the predicted tangential force and moment are relatively small. In particular, for a Newtonian fluid, the accurate solution and the simplified approximate solution degenerate to the asymptotic solution of Goldman et al. [1967. Slow viscous motion of a sphere parallel to a plane wall-motion through a quiescent fluid. Chemical Engineering Science 22, 637-651.] and O’Neill and Stewartson [1967. On the slow motion of a sphere parallel to a nearby plane wall. Journal of Fluid Mechanics 27, 705-724.]. Both solutions predict that for shear thickening fluid (n>1), the hydrodynamic force converged in the inner region of the gap between the two spheres and the contribution from the outer region is sufficiently small. For shear thinning fluid (n<1), the contribution from the outer region is also significant.     

Effect of Drag‐Reducing Polymer on the Rate of Cementation of Copper Ion on Zinc Pellets

The effect of Polyox WSR‐301 drag‐reducing polymer on the rate of cementation of copper from copper sulfate solution over packed beds of zinc pellets was investigated. The mass transfer coefficient was found to increase with increasing superficial liquid velocity for solution free of polymer. The mass transfer data were correlated in the absence of drag‐reducing polymer, using the following equation J_d = 12.55 Re^–0.5 for the conditions 10 < Re < 1970 and 1265 < Sc < 1393. In the presence of polymer, starting from a Reynolds number (Re) 550, the rate of mass transfer was found to decrease by an amount ranging from 7.5 to 51 % depending on Re and polymer concentrations. The percentage decrease in the rate of mass transfer increased with increasing Re, passed through a maximum at Re = 1400 and then decreased rapidly with further increase in Re.     

Power Requirement for Mixing Shear‐Thinning Fluids with a Planetary Mixer

The optimal control of processes dealing with non‐Newtonian liquids requires the knowledge and control of the power demand of the mixing equipment. In this context, an extension of the Metzner and Otto concept to planetary mixers is proposed to adapt this concept to planetary mixers. The theoretical part of this work defines modified expressions of Reynolds and power numbers. These definitions introduce a characteristic velocity u_ch that is used to define the parameter K_s. A planetary mixer is employed to experimentally ascertain this guideline. Power consumption measurements carried out by mixing shear‐thinning fluids permit to determine the K_s factor. This factor varies only slightly with the flow behavior index and may be regarded as a defined constant for this geometry. Finally, experiments with an additional shear‐thickening fluid confirm the validity of this approach.     

Effects of Viscous Dissipation on Heat Transfer between an Array of Long Circular Cylinders and Power Law Fluids

The free surface model has been combined with the equations of motion and of thermal energy to investigate the role of viscous dissipation on heat transfer between banks of long cylinders and power law (shear‐thinning and shear‐thickening) fluids. The equations of motion cast in the stream function/vorticity formulation have been solved numerically using a second‐order accurate finite difference method to obtain extensive information on the behaviour of local and surface‐averaged Nusselt numbers over a range of Reynolds numbers 1 – 500, for a wide range of power law indices (0.4 ≤ n ≤ 2.0), Brinkman numbers (0 ≤ Br ≤ 5) and Prandtl numbers (Pr = 1, 1000) at two representative solid volume fractions corresponding to the porosities of e = 0.4 and 0.9. Two different thermal boundary conditions are considered at the cylinder surface: constant temperature (CT) and constant heat flux (CHF). The results presented herein provide a fundamental knowledge about the influence of viscous dissipation on the heat transfer characteristics. The results reported herein further show that the effect of Brinkman number on heat transfer is strongly conditioned by the thermal boundary condition, Prandtl number and the power law index.