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.     

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.     

Flow of non‐Newtonian polymeric solutions through fibrous media

The equations of motion (continuity and momentum) describing the steady flow of incompressible power law liquids in a model porous medium consisting of an assemblage of long cylinders have been solved numerically using the finite difference method. The field equations as well as the pertinent boundary conditions have been re‐cast in terms of the stream function and vorticity. The inter‐cylinder interactions have been simulated using a simple “concentric cylinders” cell model. Extensive information on the detailed structure of the flow field in terms of the surface vorticity distribution, streamlines, and viscosity distribution on the surface of the solid cylinder as well as on the values of the pressure and friction drag coefficients under wide ranges of physical (0.4 ≤ ϵ ≤ 0.95; 1 ≥ n ≥ 0.4) and kinematic (0.01 ≤ Re ≤ 10) conditions have been obtained. The numerical results presented herein have been validated using the experimental results for the flow of Newtonian and power law fluids available in the literature; the match between the present predictions and the experiments was found to be satisfactory. © 2000 John Wiley & Sons, Inc. J Appl Polym Sci 76: 1171–1185, 2000     

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.     

Experimental study of flow regimes in three‐dimensional confined impinging jets reactor

Dynamic behaviors in a three‐dimensional confined impinging jets reactor (CIJR) were experimentally studied by a flow visualization technique at 100 ≤ Re ≤ 2000 and 2 ≤ D/d ≤ 12 (where D is the reactor diameter and d is the nozzle diameter). The effects of inlet Reynolds numbers (Re) and geometry configurations of the CIJR on the flow regimes have been investigated by a particle image velocimetry and a high‐speed camera. Results show that with the increasing Re, a segregated flow regime, a radial deflective oscillation, an axial oscillation and a vortex shedding regime emerge in turns in CIJR. A map of parameter space formed by the inlet Reynolds number (Re) and the normalized reactor diameter (D/d) has been presented. The effects of jet instability and confined boundary of the chamber on the flow regimes and their transition are also investigated and discussed. © 2014 American Institute of Chemical Engineers AIChE J, 60: 3033–3045, 2014     

Viscous flow through particle assemblages at intermediate reynolds numbers — a cell model for transport in bubble swarms

The fluid mechanical behaviour of a bubble swarm was simulated using a cell model. The Navier-Stokes equations were solved numerically for the liquid flow in a uniform assemblage of circulating, spherical bubbles. Ranges of parameters studied included, Reynolds numbers, 0–1000 and porosities, 0.4–1. The numerical calculations show the effects of variations in Reynolds numbers and porosity on: surface vorticity and pressure distributions and form and friction drag coefficients. For all Reynolds numbers investigated a standing vortex ring was absent Predicted drag coefficients and Sherwood or Nusselt numbers agree with limiting analytical solutions for low and high Reynolds numbers. The theoretical results show good agreement with experimental data for porosity as a function of superficial gas velocity. Predicted and measured Sherwood and Nusselt Numbers were in substantial disagreement, making detailed comparison unwarranted The calculations should also be valid for dispersions of uniform, circulating, spherical droplets for the special case where the droplet viscosity is much less than the viscosity of the continuous fluid     

Pressure drop in a split‐and‐recombine caterpillar micromixer in case of newtonian and non‐newtonian fluids

Microreactors are very promising tools for the design of future chemical processes. For example, emulsions of very narrow size distribution are obtained at much lower energy consumption than the one spent with usual processes. Micromixers play thereby an eminent role. The goal of this study is to better understand the hydrodynamic properties of a split‐and‐recombine Caterpillar micromixer (CPMM) specially with regard to handling viscoelastic fluids, a topic hardly addressed so far in the context of micromixers in general, although industrial fluids like detergent, cosmetic, or food emulsions are non‐Newtonian. Friction factor was measured in a CPMM for both Newtonian and non‐Newtonian fluids. For Newtonian fluids, the friction factor in the laminar regime is f/2 = 24/Re. The laminar regime exists up to Reynolds numbers of 15. For shear‐thinning fluids like Carbopol 940 or viscoelastic fluids like Poly Acryl Amide (PAAm) aqueous solutions, the friction factor scales identically within statistical errors up to a generalized Reynolds number of 10 and 0.01, respectively. Above that limit, there is an excess pressure drop for the viscoelastic PAAm solution. This excess pressure drop multiplies the friction factor by more than a decade over a decade of Reynolds numbers. The origin of this excess pressure drop is the high elongational flow present in the Caterpillar static mixer applied to a highly viscoelastic fluid. This result can be extended to almost all static mixers, because their flows are generally highly elongational. © 2013 American Institute of Chemical Engineers AIChE J, 59: 2679–2685, 2013     

Steady incompressible fluid flow over a bundle of cylinders at moderate reynolds numbers

The complete Navier-Stokes equations describing the steady flow of incompressible Newtonian fluids normal to an array of long cylinders have been solved numerically using the finite difference method in terms of the stream function and vorticity variables. The inter-cylinder interactions have been mimicked using the well known free surface cell model. Extensive information on the detailed structure of the flow field in terms of the surface vorticity distribution, stagnation pressure, stream line and iso-vorticity line plots, as well as on the values of the integral quantities, such as pressure, friction and total drag coefficients, have been obtained under wide ranges of conditions as follows: 0.3 ≤ ? 0.99 and 0.01 ≤ Re ≤ 100. The numerical results presented herein have been validated using the appropriate theoretical and experimental results available in the literature; the match between the present predictions and the scant experimental results is good.     

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.     

Slow non-newtonian flow through packed beds: Effect of zero shear viscosity

The slow non-Newtonian (inelastic) flow through packed beds of mono-size spherical particles has been simulated by solving the equations of motion numerically. The inter-particle interactions have been modelled by using a simple cell model. Theoretical estimates of pressure, friction and total drag coefficients as function of the pertinent physical (l≥n≥ 0.2; 0.3 ≤ e ≤ 0.5) and kinematic parameters (0.01 ≤ '≤ 100) for a fixed value of Reynolds number {Re = 0.001) have been obtained. The theoretical predictions reported herein have been validated using the suitable experimental results available in the literature, and the importance of including the zero shear viscosity in analyses for the creeping flow problems is convincingly demonstrated.     

Low Reynolds number flow of power-law fluids over two square cylinders in tandem

The governing partial differential equations have been solved numerically for the 2-D and steady powerlaw fluid flow over two square cylinders in tandem arrangement. Extensive numerical results spanning wide ranges of the governing parameters as Reynolds number (0.1≤Re≤40), power-law index (0.2≤n≤1) and inter-cylinder spacing (2≤L/d≤6) are presented herein; limited results for L/d=24 are also obtained to approach the single cylinder behavior. The detailed flow visualization is done by means of the streamline and vorticity contours in the vicinity of two cylinders. The global characteristics are analyzed in terms of the surface pressure distribution and pressure drag coefficient. The drag coefficient shows the classical inverse dependence on the Reynolds number irrespective of the value of the powerlaw index; the drag on the upstream cylinder is always greater than that for the downstream cylinder.     

Power Requirement when Mixing a Shear‐Thickening Fluid with a Helical Ribbon Impeller Type

The optimal design of close clearance impellers requires the knowledge of the power demand of the mixing equipment. In non‐Newtonian mixing, this can be readily obtained using the Metzner and Otto concept [1]. In this work, this concept and the determination of the K_s value for an atypical helical agitator (PARAVISC system from Ekato firm) have been revised in the case of shear‐thinning fluids and a shear‐thickening fluid. For poor shear‐thinning fluids, it has been shown that for our mixing system the K_s value does not vary strongly with the flow behavior index, and may be regarded as a constant for the mixing purpose design. By contrast, for the shear‐thickening fluid, power consumption measurements indicate that the relationship between the K_s values and the flow behavior index is much more complex due to a partial solidification of the product around the impeller.     

Viscoelastic Liquids in Stirred Vessels – Part I: Power Consumption in Unaerated Vessels

Different shear‐thinning and elastic fluids (STE fluids) have been stirred under unaerated conditions, in vessels equipped with Rushton disc turbines. Their power consumption has been evaluated over a wide range of stirring rates and their Metzner‐Otto constant (k_s) has been measured. A correlation has then been proposed to predict k_s values for a Rushton turbine operating in non‐Newtonian solutions. Power curves of STE fluids have been drawn and compared with reference curves (Newtonian, shear‐thinning inelastic and elastic with constant shear viscosity fluids). The STE fluids have thus been divided into two categories. The STE fluids of the first category (STE I fluids), which are concentrated viscous solutions of polymers (guar, CMC) reducing the power consumption at the beginning of the transitional region and connecting with the Newtonian reference at higher Reynolds numbers. In contrast, STE solutions of the second category (STE II fluids), which are solutions of drag reducing polymers (PAA), are less viscous and more elastic. They reduce the power consumption at the end of the transitional region and do not connect with the Newtonian reference, at least until Re = 6000. A general correlation has finally been proposed to model the power curve of STE fluids stirred by a Rushton turbine from the laminar to the turbulent regions, as a function of their elasticity.     

Pressure drop for the slow flow of dilatant fluids through a fixed bed of spherical particles

The equations of motion in conjunction with the free surface cell model have been solved numerically for the creeping flow of dilatant fluids through a fixed bed of spheres. The effect of the flow behaviour index (1 ≤ n ≤ 2) on the individual as well as the total drag coefficients is elucidated. Finally, the theoretical predictions have been validated using the scant experimental results available in the literature.     

Convective heat transfer for power law fluids in packed and fluidised beds of spheres

The forced convection heat transfer characteristics for an incompressible and steady flow of power law liquids in fixed and extended beds of spherical particles has been studied numerically. The sphere-sphere hydrodynamic interactions have been accounted for by using a simple cell model. Within the framework of such a cell model, the momentum and energy equations have been solved using a finite difference method to obtain the velocity and temperature fields. Extensive numerical estimates of the local and average Nusselt numbers as functions of the physical, rheological and kinematic variables have been presented and discussed for the two commonly employed thermal boundary conditions. In broad terms, the Nusselt number for power law fluids (both shear-thinning and shear-thickening conditions) normalized with respect to the corresponding value for a Newtonian fluid shows weak additional dependence on the power law flow behaviour index. The shear-thinning behaviour is seen to promote heat transfer and as expected the shear-thickening behaviour impedes heat transfer in fixed and fluidised beds. All in all, the present results encompass wide ranges of conditions as follows: Reynolds number: 1-500; Peclet number: 1-500; bed voidage: 0.4-0.8 and the flow behaviour index: 0.5-1.8 thereby covering extremely shear-thinning and shear-thickening types of fluid behaviours. The paper is concluded by presenting detailed comparisons with the limited analytical and/or experimental results available for liquid-solid mass transfer in such systems.     

Numerical studies of viscous stagnation point flows

A rapid finite-difference procedure for computing stagnation point flows at a given Reynolds number in terms of stream function and vorticity distributions is described. The method is illustrated for spherical and cylindrical curvature for Reynolds number ranges of 1 ≤ Re ≤ 400 and 4 ≤ Re ≤ 100 respectively. Predicted radial and angular velocity distributions in the frontal region are compared with results from numerical solutions of the full elliptic equations of motion. The use of a modified potential flow outer boundary condition reduces the region of numerical computation.     

Effect of Non‐Newtonian Characteristics on Convective Liquid‐Solid Heat Transfer in Packed and Fluidised Beds of Spherical Particles

In this work, the field (continuity, momentum and thermal energy) equations togetherwith a cell model have been solved numerically to elucidate the influence of non‐Newtonian (Power law rheology) liquid characteristics on liquid‐solid heat transfer in packed and fluidised beds of spherical particles. The results presented herein relate to wide ranges of conditions of bed voidage, power‐law index and Peclet number but are limited to low Reynolds number (≤1) flow conditions. Within the range of conditions, the effect of power‐law index is found to be small and this is also consistent with the available experimental results on liquid‐solid mass transfer in these systems.     

Power‐Law Fluid Flow Over a Sphere: Average Shear Rate and Drag Coefficient

Based on the consideration of the rate of mechanical energy dissipation, an expression for the average shear rate for a sphere falling in a power‐law fluid in the creeping flow regime has been deduced. The average shear rate in a power‐law fluid (n<1) appears to be higher than that in an equivalent Newtonian fluid. This in turn has been combined with the numerical predictions of drag coefficient (up to Reynolds number of 100) of a sphere to develop a generalized drag correlation for power‐law liquids encompassing both n > 1 and n < 1 which appears to apply up to much higher values of the Reynolds number. The available experimental data have been used to demonstrate the reliability and accuracy of the new correlation for shearthinning liquids. Also, in the limit of n = 1, this expression reproduces the standard drag curve with a very high accuracy.     

FLOW OF NON-NEWTONIAN POWER LAW LIQUIDS THROUGH PACKED AND FLUIDIZED BEDS

The creeping flow of power law liquids through assemblages of spherical particles has been studied theoretically. The inter-particle interactions are modelled via the zero vorticity cell model. The governing equations have been solved numerically to obtain the theoretical estimates of the drag force experienced by an assemblage placed in a streaming power law fluid. The results reported herein encompass wide ranges of fluid behaviour (values of power law index) and bed voidage thereby covering packed and fluidized bed conditions. Detailed comparisons with experimental data suggest that this theory can be used to predict pressure drop for power law fluid flow in packed beds as well as velocity-bed expansion characteristics for a fluidized bed. By analogy with the Newtonian case, intuitively, one would expect these results to be applicable to hindered settling in power law fluids, and indeed this is borne out by the limited amount of data covering the range 1≥n≥0.8 available in the literature.     

Mixed convection flow of non‐Newtonian fluid from a slotted vertical surface with uniform surface heat flux

In the present paper, the combined convection flow of an Ostwald–de Waele type power‐law non‐Newtonian fluid past a vertical slotted surface has been investigated numerically. The boundary condition of uniform surface heat flux is considered. The equations governing the flow and the heat transfer are reduced to local non‐similarity form. The transformed boundary layer equations are solved numerically using implicit finite difference method. Solutions for the heat transfer rate obtained for the rigid surface compare well with those documented in the published literature. From the present analysis, it is observed that, an increase in χ leads to increase in skin friction as well as reduction in heat transfer at the surface. As the power‐law index n increases, the friction factor as well as heat transfer increase.