Displacement of a Newtonian fluid by a non-Newtonian fluid in a porous medium

This paper presents an analytical Buckley-Leverett-type solution for one-dimensibnal immiscible displacement of a Newtonian fluid by a non-Newtonian fluid in porous media. The non-Newtonian fluid viscosity is assumed to be a function of the flow potential gradient and the non-Newtonian phase saturation. To apply this method to field problems a practical procedure has been developed which is based on the analytical solution and is similar to the graphic technique of Welge. Our solution can be regarded as an extension of the Buckley-Leverett method to Non-Newtonian fluids. The analytical result reveals how the saturation profile and the displacement efficiency are controlled not only by the relative permeabilities, as in the Buckley-Leverett solution, but also by the inherent complexities of the non-Newtonian fluid. Two examples of the application of the solution are given. One application is the verification of a numerical model, which has been developed for simulation of flow of immiscible non-Newtonian and Newtonian fluids in porous media. Excellent agreement between the numerical and analytical results has been obtained using a power-law non-Newtonian fluid. Another application is to examine the effects of non-Newtonian behavior on immiscible displacement of a Newtonian fluid by a power-law non-Newtonian fluid.     

非牛顿流体力学研究的若干进展

简要介绍微极性流体,应力偶流体,非整数阶Maxwell流体和智能流体的基本概念,通过简单的例子阐述它们与牛顿流体的本质差别,并介绍非牛顿流体在以上热点研究领域的若干进展．     

基于舌/上颚微间隙下流体流动行为研究

为探讨口腔环境下流体的流动行为,采用数值方法与流变试验深入研究舌/上颚微间隙下流体流量的影响因素. 建立舌/上颚微间隙的简化模型及Reynolds方程,通过数值方法获取微间隙下流量变化;在DHR-2流变仪上研究非牛顿流体的黏度与剪切率的变化,探讨牛顿流体和非牛顿流体的流量影响. 结果表明:牛顿流体流量平方的倒数同载荷和黏度比值和时间均呈线性函数关系;所制备的非牛顿流体近似为幂律流体,其黏度随脂肪含量的增加而增大,而非牛顿流体流量率先高于后低于等效牛顿流体,其研究结果将为特定人群功能产品的研发提供技术支持.      

Free convection through conducting non-Newtonian fluids over a thin axisymmetric body

An analysis of free convection heat transfer in electrically conducting power law non-Newtonian fluid over a thin axisymmetric body of constant temperature is carried out. The uniform external magnetic field acts normally to the surface through the induced boundary layer. In view of the fact that most of the non-Newtonian fluids have large Prandtl number, the momentum equation is simplified. The equations of conservation of mass, momentum and energy which govern and describe the flow and heat transfer are solved numerically. The effect of the magnetic field on the velocity, temperature, the coefficient of friction and the Nusselt number are investigated. Numerical results are tabulated, presented graphically and discussed.     

Non Newtonian annular alloy solidification in mould

The annular solidification of an aluminium–silicon alloy in a graphite mould with a geometry consisting of horizontal concentric cylinders is studied numerically. The analysis incorporates the behavior of non-Newtonian, pseudoplastic (n?=?0.2), Newtonian (n?=?1), and dilatant (n?=?1.5) fluids. The fluid mechanics and heat transfer coupled with a transient model of convection diffusion are solved using the finite volume method and the SIMPLE algorithm. Solidification is described in terms of a liquid fraction of a phase change that varies linearly with temperature. The final results make it possible to infer that the fluid dynamics and heat transfer of solidification in an annular geometry are affected by the non-Newtonian nature of the fluid, speeding up the process when the fluid is pseudoplastic.     

Anomalous predictions of pressure drop and heat transfer in ducts of arbitrary cross-section with modified power law fluids

 The apparent viscosities of purely viscous non-Newtonian fluids are shear rate dependent. At low shear rates, many of such fluids exhibit Newtonian behaviour while at higher shear rates non-Newtonian, power law characteristics exist. Between these two ranges, the fluid's viscous properties are neither Newtonian or power law. Utilizing an apparent viscosity constitutive equation called the “Modified Power Law” which accounts for the above behavior, solutions have been obtained for forced convection flows. A shear rate similarity parameter is identified which specifies both the shear rate range for a given fluid and set of operating conditions and the appropriate solution for that range. The results of numerical solutions for the friction factor–Reynolds number product and for the Nusselt number as a function of a dimensionless shear rate parameter have been presented for forced fully developed laminer duct flows of different cross-sections with modified power law fluids. Experimental data is also presented showing the suitability of the “Modified Power Law” constitutive equation to represent the apparent viscosity of various polymer solutions. Received on 21 August 2000     

Soret and Dufour effect on double diffusion mixed convection from a vertical surface in a porous medium saturated with a non-Newtonian fluid

A non-similar boundary layer analysis is presented to study the flow, heat and mass transfer characteristics of non-Darcian mixed convection of a non-Newtonian fluid from a vertical isothermal plate embedded in a homogeneous porous medium with the effect of Soret and Dufour and in the presence of either surface injection or suction. The value of the mixed-convection parameter lies between 0 and 1. In addition, the power-law model is used for non-Newtonian fluids with exponent n < 1 for pseudoplastics n = 1 for Newtonian fluids and n > 1 for dilatant fluids. Furthermore, the coordinates and dependent variables are transformed to yield computationally efficient numerical solutions that are valid over the entire range of mixed convection, from the pure forced-convection limit to the pure free-convection limit, and the whole domain of non-Newtonian fluids, from pseudoplastics to dilatant fluids. The numerical solution of the problem is derived using a Runge–Kutta integration scheme with Newton–Raphson shooting technique. Distributions for velocity, temperature and concentration, as well as for the rate of wall heat and mass transfer, have been obtained and discussed for various physical parametric values.     

Hydrodynamic instability of extensional flow

The stability analysis for sheet stretching recently presented by Minoshima and White for Newtonian fluids is repeated for non-Newtonian fluids. For that purpose a constitutive law for nearly extensional flow is derived which, apart from a restriction to short memory of the fluid, is generally valid. Using this constitutive law the result found by Minoshima and White is generalized. Application of the general result to a Maxwell-type fluid shows that the elasticity of the fluid has a stabilizing influence.     

An analysis for forced convection heat transfer from external surfaces to non-Newtonian fluids

A theoretical analysis has been proposed for the forced convection heat transfer from external surfaces immersed in non-Newtonian fluids of the power-law model. The integral treatment previously introduced for Newtonian fluids has been successfully extended to the non-Newtonian fluids over a flat plate and a wedge of an arbitrary included angle. The integral momentum and energy equations are transformed into a pair of characteristic equations, which can readily be solved for the velocity shape factor and the boundary layer thickness ratio, once the exponents in the expressions for the power-law model, free stream velocity and wall temperature variation are specified. It has been also found that an asymptotic expression derived under the assumption of large Prandtl number, is valid practically for all power-law fluids, and hence, can be used for a speedy, and yet accurate estimation of the local heat transfer to non-Newtonian fluids.     

Studies on heat transfer to Newtonian and non-Newtonian fluids in agitated vessel

Heat transfer studies to Newtonian and non-Newtonian fluids are carried out in a stirred vessel fitted with anchor/turbine impeller and a coil for heating/cooling with an objective of determining experimentally the heat transfer coefficient of few industrially important systems namely castor oil and its methyl esters, soap solution, CMC and chalk slurries. The effect of impeller geometry, speed and aeration is investigated. Generalized Reynolds and Prandtl numbers are calculated using an apparent viscosity for non-Newtonian fluids. The data is correlated using a Sieder–Tate type equation. A trend of increase in heat transfer coefficient with RPM in presence and absence of solids has been observed. Relatively high values of Nusselt numbers are obtained for non-Newtonian fluids when aeration is coupled with agitation. The contribution of natural convection to heat transfer has been accounted for by incorporating the Grashof number. The correlations developed based on these studies are applied for design of commercial scale soponification reactor. Power per unit volume resulted in reliable design of a reactor.     

Turbulent flow of Newtonian and shear-thinning liquids through a sudden axisymmetric expansion

 Measurements are reported for the turbulent flow through a sudden expansion of a moderately elastic shear-thinning liquid and also for two Newtonian liquids. The differences in the mean velocity fields for the two fluid types are relatively small, including the length of the recirculation region which is essentially unaffected by the fluid rheology. Although turbulent kinetic energy levels for the non-Newtonian fluids are always lower than for the Newtonian fluids, no significant difference is found in the relative contributions to the turbulent kinetic energy of the axial, radial and tangential normal stresses. Since the vorticity thicknesses are much the same for all flows, viscoelasticity appears to be responsible for the reduced levels of turbulent kinetic energy for the non-Newtonian fluids. Received: 6 November 1998/Accepted: 27 January 1999     

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.     

Mixed convection boundary layer flow of non-Newtonian fluids along vertical wavy plates

Mixed convection boundary layer flows of non-Newtonian fluids over the wavy surfaces are studied by the coordinate transformation and the cubic spline collocation numerical method. The effects of the wavy geometry, the buoyancy parameter and the generalized Prandtl number for pseudoplastic fluids, Newtonian fluids and dilatant fluids on the skin-friction coefficient, local and mean Nusselt numbers have been graphically studied. Results show that both higher generalized Prandtl numbers and buoyancy parameters are seen to enhance the influence of wavy surfaces on the local Nusselt number, irrespective of whether the fluids are Newtonian fluids or non-Newtonian fluids. Moreover, the irregular surfaces have higher total heat flux than that of corresponding flats plate for any fluid.     

Interfacial instability of turbulent two-phase stratified flow: Pressure-driven flow and non-Newtonian layers

We derive a flat-interface model to describe the flow of two horizontal, stably stratified fluids, where the bottom layer exhibits non-Newtonian rheology. The model takes into account the yield stress and power-law nature of the bottom fluid. In the light of the large viscosity contrast assumed to exist across the fluid interface, and for large pressure drops in the streamwise direction, the possibility for the upper Newtonian layer to display fully developed turbulence must be considered, and is described in our model. We develop a linear-stability analysis to predict the conditions under which the flat-interface state becomes unstable, and pay particular attention to characterizing the influence of the non-Newtonian rheology on the instability. Increasing the yield stress (up to the point where unyielded regions form in the bottom layer) is destabilizing; increasing the flow index, while bringing a broader spectrum of modes into play, is stabilizing. In addition, a second mode of instability is found, which depends on conditions in the bottom layer. For shear-thinning fluids, this second mode becomes more unstable, and yet more bottom-layer modes can become unstable for a suitable reduction in the flow index. One further difference between the Newtonian and non-Newtonian cases is the development of unyielded regions in the bottom layer, as the linear wave on the interface grows in time. These unyielded regions form in the trough of the wave, and can be observed in the linear analysis for a suitable parameter choice.     

A correlation for yield stress fluid flow through packed beds

Relatively few correlations are available for non-Newtonian fluid flows through packed beds, even though such fluids are frequently used in industry. In this paper, a correlation is presented for yield stress fluid flow through packed beds. The correlation is developed by introducing the yield stress model in place of the Newtonian model used in deriving Erguns equation. The resulting model has three parameters that are functions of the geometry and roughness of the particle surfaces. Two of the parameters can be deduced in the limit as the yield stress becomes negligible and the model reduces to Erguns equation for Newtonian fluids. The third model parameter is determined from experimental data. The correlation relates a defined friction factor to the dimensionless Reynolds and Hedstrom numbers and can be used to predict pressure drop for flow of a yield stress fluid through a packed bed of spherical particles. Conditions for flow or no-flow are also determined in the correlation. Comparison of model calculations, between a Newtonian and a yield stress fluid for flow penetration into a packed bed of spheres, shows the yield stress fluid initially performs similar to the Newtonian fluid, at large Reynolds numbers. At lower Reynolds numbers the yield stress effect becomes important and the flow rate significantly decreases when compared to the Newtonian fluid.     

Thermal chaotic mixing of power-law fluids in a mixer with alternately rotating walls

We investigate the enhancement of both mixing and heat transfer in a two-rod mixer for highly viscous non-Newtonian fluids. The mixer is composed of two vertical circular rods in a cylindrical tank. Chaotic flows are obtained by imposing the temporal modulations of the rotational velocities of the walls. We study the effects of different stirring protocols, which lead to non-chaotic and chaotic flows, on the efficiency of both mixing and heat transfer for three different rheological fluid behaviors: shear-thinning, Newtonian and shear-thickening. For this purpose, we use statistical indicators that characterize the mean value of the fluid temperature and its homogenization. We find that chaotic mixing is suitable for shear-thickening fluids for which we observe a clear enhancement of the thermal mixing (heat extraction and homogenization). This is due to the increase in the apparent fluid viscosity in the vicinity of the rotating walls. This aspect confirms the relevance of chaotic mixing for highly viscous fluids.     

Stability of triple diffusive convection in a viscoelastic fluid-saturated porous layer

The triple diffusive convection in an Oldroyd-B fluid-saturated porous layer is investigated by performing linear and weakly nonlinear stability analyses. The condition for the onset of stationary and oscillatory is derived analytically. Contrary to the observed phenomenon in Newtonian fluids, the presence of viscoelasticity of the fluid is to degenerate the quasiperiodic bifurcation from the steady quiescent state. Under certain conditions, it is found that disconnected closed convex oscillatory neutral curves occur, indicating the requirement of three critical values of the thermal Darcy-Rayleigh number to identify the linear instability criteria instead of the usual single value, which is a novel result enunciated from the present study for an Oldroyd-B fluid saturating a porous medium. The similarities and differences of linear instability characteristics of Oldroyd-B, Maxwell, and Newtonian fluids are also highlighted. The stability of oscillatory finite amplitude convection is discussed by deriving a cubic Landau equation, and the convective heat and mass transfer are analyzed for different values of physical parameters.     

Studies on the axisymmetric sphere–sphere interaction problem in Newtonian and non-Newtonian fluids

In this research, experimental studies have been performed on the hydrodynamic interaction between two spheres by using particle image velocimetry and measuring the force between the spheres. To approach the system as a resistance problem, a servo-driving system was set-up by assembling a microstepping motor, a ball screw and a linear motion guide for the particle motion. Glycerin and a dilute solution of polyacrylamide in glycerin were used as Newtonian and non-Newtonian fluids, respectively. The polymer solution behaves like a Boger fluid when the concentration is 1000 ppm or less. The experimental results were compared with the asymptotic solution of Stokes equation. The result shows that fluid inertia and unsteadiness play important roles in the particle–particle interaction in the Newtonian fluid. This implies that the motion of two particles in suspension is not reversible even in the Newtonian fluid. In the non-Newtonian fluid, in addition to inertial effect, normal stress differences and viscoelasticity play important roles as expected. In dilute solutions weak shear thinning and the migration of polymer molecules in the inhomogeneous flow field also appear to affect the physics of the problem.     

Secondary flow in the neighborhood of a cylinder oscillating in a viscoelastic fluid

Measurements were made of the steady secondary flow which is caused by a cylinder executing small-amplitude oscillations perpendicular to the cylinder generators. Test fluids were a water-glycerine mixture and two non-Newtonian fluids: aqueous solutions of Separan AP-30 and PolyHall 295 polyacrylamide polymers. The fluids were contained in a cylinder whose center coincides with the rest position of the oscillating cylinder. Experimental results are consistent with a theoretical analysis for a three-constant Oldroyd fluid model in its simplified convected-Maxwell form. Both experiment and theory show that at an oscillation frequency of 40 Hz the secondary flow of the dilute polymer solutions is, essentially, in a sense opposite from that of a Newtonian fluid.     

On axisymmetric flow of Sisko fluid over a radially stretching sheet

This study presents an analysis of the axisymmetric flow of a non-Newtonian fluid over a radially stretching sheet. The momentum equations for two-dimensional flow are first modeled for Sisko fluid constitutive model, which is a combination of power-law and Newtonian fluids. The general momentum equations are then simplified by invoking the boundary layer analysis. Then a non-linear ordinary differential equation governing the axisymmetric boundary layer flow of Sisko fluid over a radially stretching sheet is obtained by introducing new suitable similarity transformations. The resulting non-linear ordinary differential equation is solved analytically via the homotopy analysis method (HAM). Closed form exact solution is then also obtained for the cases n=0 and 1. Analytical results are presented for the velocity profiles for some values of governing parameters such as power-law index, material parameter and stretching parameter. In addition, the local skin friction coefficient for several sets of the values of physical parameter is tabulated and analyzed. It is shown that the results presented in this study for the axisymmetric flow over a radially non-linear stretching sheet of Sisko fluid are quite general so that the corresponding results for the Newtonian fluid and the power-law fluid can be obtained as two limiting cases.