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.     

Flow of viscoelastic fluids through a porous channel—I

The flow of viscoelastic fluids through a porous channel with one impermeable wall is computed. The flow is characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. Three solutions are developed: (i) an exact numerical solution, (ii) a perturbation solution for small R, the cross-flow Reynold's number and (iii) an asymptotic solution for large R. The results from exact numerical integration reveal that the solutions for a non-Newtonian fluid are possible only up to a critical value of the viscoelastic fluid parameter, which decreases with an increase in R. It is further demonstrated that the perturbation solution gives acceptable results only if the viscoelastic fluid parameter is also small. Two more related problems are considered: fluid dynamics of a long porous slider, and injection of fluid through one side of a long vertical porous channel. For both the problems, exact numerical and other solutions are derived and appropriate conclusions drawn.     

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.     

Three‐dimensional boundary element analysis of drop deformation in confined flow for Newtonian and viscoelastic systems

An adaptive (Lagrangian) boundary element approach is proposed for the general three‐dimensional drop deformation in confined flow. The adaptive method is stable as it includes remeshing capabilities of the deforming interface between drop and suspending fluid, and thus can handle large deformations. Both drop and surrounding fluid are viscous incompressible and can be Newtonian or viscoelastic. A boundary‐only formulation is implemented for fluids obeying the linear Jeffrey's constitutive equation. Similarly to the formulation for two‐dimensional Newtonian fluids (Khayat RE, Luciani A, Utracki LA. Boundary element analysis of planar drop deformation in confined flow. Part I. Newtonian fluids. Engineering Analysis of Boundary Elements 1997; 19 : 279), the method requires the solution of two simultaneous integral equations on the interface between the two fluids and the confining solid boundary. Although the problem is formulated for any confining geometry, the method is illustrated for a deforming drop as it is driven by the ambient flow inside a cylindrical tube. The accuracy of the method is assessed by comparison with the analytical solution for two‐phase radial spherical flow, leading to good agreement. The influence of mesh refinement is examined for a drop in simple shear flow. Copyright © 2000 John Wiley & Sons, Ltd.     

Extension of the Darcy�CForchheimer Law for Shear-Thinning Fluids and Validation via Pore-Scale Flow Simulations

Flow of non-Newtonian fluids through porous media at high Reynolds numbers is often encountered in chemical, pharmaceutical and food, as well as petroleum and groundwater engineering, and in many other industrial applications. Under the majority of operating conditions typically explored, the dependence of pressure drops on flow rate is non-linear and the development of models capable of describing accurately this dependence, in conjunction with non-trivial rheological behaviors, is of paramount importance. In this work, pore-scale single-phase flow simulations conducted on synthetic two-dimensional porous media are performed via computational fluid dynamics for both Newtonian and non-Newtonian fluids and the results are used for the extension and validation of the Darcy?CForchheimer law, herein proposed for shear-thinning fluid models of Cross, Ellis and Carreau. The inertial parameter ?? is demonstrated to be independent of the viscous properties of the fluids. The results of flow simulations show the superposition of two contributions to pressure drops: one, strictly related to the non-Newtonian properties of the fluid, dominates at low Reynolds numbers, while a quadratic one, arising at higher Reynolds numbers, is dependent on the porous medium properties. The use of pore-scale flow simulations on limited portions of the porous medium is here proposed for the determination of the macroscale-averaged parameters (permeability K, inertial coefficient ?? and shift factor ??), which are required for the estimation of pressure drops via the extended Darcy?CForchheimer law. The method can be applied for those fluids which would lead to critical conditions (high pressures for low permeability media and/or high flow rates) in laboratory tests.     

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.     

A note on the translation of a thin rod inside a cylinder

The problem of a thin rod moving longitudinally along the axis of symmetry of a cylindrical vessel is examined for Newtonian and non-Newtonian liquids. For non-Newtonian fluids, the inelastic power-law type solution predicts the experimental results particularly well. On account of wall effects, the induced pressure gradients are much greater for a Newtonian fluid than for a viscoelastic fluid. In fact, in the latter case, they may be considered negligible when the radius of the inner cylinder is small compared to the one of the outer cylinder.     

Flow of Newtonian and Non-Newtonian Fluids Through Packed Beds: An Experimental Study

The phenomena of flow reduction and flow enhancement was observed in case of viscoelastic and viscoinelastic fluids flowing through packed beds, respectively. In this study, the pressure drop-flow rate behaviors for the flow of Newtonian (water), non-Newtonian viscoinelastic (Carboxy methyl cellulose solution in water, CMC) and viscoelastic (Polyacrylamide solution in water, PAA) fluids have been experimentally studied and pressure drop behavior compared with existing models for viscoinelastic and viscoelastic fluids. Based on the observed data, an appropriate empirical correlation for pressure drop prediction in case of non-Newtonian fluid flowing through packed bed has been proposed. The correlation predicts the data well to within a reasonable accuracy.     

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.     

Frottements et pertes de pression des fluides non newtoniens dans des conduites non circulaires

The study of fluid flow in a duct requires characteristic parameters of the flow and dimensionless numbers to correlate and compare experimental results. For Newtonian fluids in simple configurations, the definition of the Reynolds number is quite standard, but for non-Newtonian fluid flows in ducts with arbitrary shape of cross section, the dependence of the apparent viscosity with the shear rate requires a generalization of this dimensionless number. This note proposes a general method valid for a large class of non-Newtonian fluids and for all duct shapes. An application is developed for a viscoelastic flow through a rectangular duct. Results obtained in the present investigation are in a good agreement with available correlations. To cite this article: M. Mahfoud et al., C. R. Mecanique 333 (2005).     

Constant heat flux solution for mixed convection boundary layer viscoelastic fluid

The mixed convection boundary layer of a viscoelastic fluid past a circular cylinder with constant heat flux is discussed. The boundary layer equations are an order higher than those for the Newtonian (viscous) fluid and the adherence boundary conditions are insufficient to determine the solution of these equations completely. The governing non-similar partial differential equations are transformed into dimensionless forms and then solved numerically using the Keller-box method by augmenting an extra boundary condition at infinity. Numerical results obtained in the form of velocity distributions and temperature profiles are presented for a range of values of the dimensionless viscoelastic fluid parameter. It is found that for some values of the viscoelastic parameter and some negative values of the mixed convection parameter (opposing flow) the momentum boundary layer separates from the cylinder. Heating the cylinder delays separation and can, if the cylinder is warm enough, suppress the separation completely. Similar to the case of a Newtonian fluid, cooling the cylinder brings the separation point nearer to the lower stagnation point.     

On the dispersion of a solute in Poiseuille flow of a third-grade fluid in a channel

Gill and Sankarasubramanian's analysis of the dispersion of Newtonian fluids in laminar flow between two parallel walls are extended to the flow of non-Newtonian viscoelastic fluid (known as third-grade fluid) using a generalized dispersion model which is valid for all times after the solute injection. The exact expression is obtained for longitudinal convective coefficient K₁(Γ), which shows the effect of the added viscosity coefficient Γ on the convective coefficient. It is seen that the value of the K₁(Γ) for Γ≠0 is always smaller than the corresponding value for a Newtonian fluid. Also, the effect of the added viscosity coefficient on the K₂(t,Γ) (which is a measure of the longitudinal dispersion coefficient of the solute) is explored numerically. Finally, the axial distribution of the average concentration C_m of the solute over the channel cross-section is determined at a fixed instant after the solute injection for several values of the added viscosity coefficient.     

Exact solutions for extended boundary value problems

Summary Extended definition of a stress tensor for a non-Newtonian fluid brings in higher degree derivatives with coefficients as powers of non-Newtonian parameter in the differential equations of motion. Yet, these differential equations need to be solved subject to the same boundary conditions as in the corresponding Newtonian flow problem. A technique is developed to obtain exact solutions for such an extended boundary value problem. Some flow problems forWalters liquidB are considered.     

Dynamics of stochastic non-Newtonian fluids driven by fractional Brownian motion with Hurst parameter H ∈ （1/4,1/2）

A two-dimensional (2D) stochastic incompressible non-Newtonian fluid driven by the genuine cylindrical fractional Brownian motion (FBM) is studied with the Hurst parameter $H \in \left( {\tfrac{1} {4},\tfrac{1} {2}} \right)$ under the Dirichlet boundary condition. The existence and regularity of the stochastic convolution corresponding to the stochastic non-Newtonian fluids are obtained by the estimate on the spectrum of the spatial differential operator and the identity of the infinite double series in the analytic number theory. The existence of the mild solution and the random attractor of a random dynamical system are then obtained for the stochastic non-Newtonian systems with $H \in \left( {\tfrac{1} {2},1} \right)$ without any additional restriction on the parameter H.     

Mass transfer from small capillaries with wall resistance in the laminar flow regime

The Graetz problem in heat transfer is extended to the analysis of mass transfer in circular ducts for the cases where wall resistance is included and where non-Newtonian fluids that obey Casson's equation are considered. The eigenvalues and fluid bulk coefficients are presented for the fluid between the extremes of Newtonian and slug flow. It is found that for fluids which are only slightly non-Newtonian, such as blood, which is closely approximated by Casson's equation, the mass transfer rate can be predicted by Newtonian fluid analysis without appreciable error. Some experimental results give support to the theory.     

Self-similar solutions in the problem of generalized vortex diffusion

The analysis of the group properties and the search for self-similar solutions in problems of mathematical physics and continuum mechanics have always been of interest, both theoretical and applied [1–3]. Self-similar solutions of parabolic problems that depend only on a variable of the type η = x/√t are classical fundamental solutions of the one-dimensional linear and nonlinear heat conduction equations describing numerous physical phenomena with initial discontinuities on the boundary [4]. In this study, the term “generalized vortex diffusion” is introduced in order to unify the different processes in mechanics modeled by these problems. Here, vortex layer diffusion and vortex filament diffusion in a Newtonian fluid [5] can serve as classical hydrodynamic examples. The cases of self-similarity with respect to the variable η are classified for fairly general kinematics of the processes, physical nonlinearities of the medium, and types of boundary conditions at the discontinuity points. The general initial and boundary value problem thus formulated is analyzed in detail for Newtonian and non-Newtonian power-law fluids and a medium similar in behavior to a rigid-ideally plastic body. New self-similar solutions for the shear stress are derived.     

The Falkner–Skan Wedge Flows of Power-Law Fluids Embedded in a Porous Medium

The non-Darcy flow characteristics of power-law non-Newtonian fluids past a wedge embedded in a porous medium have been studied. The governing equations are converted to a system of first-order ordinary differential equations by means of a local similarity transformation and have been solved numerically, for a number of parameter combinations of wedge angle parameter m, power-law index of the non-Newtonian fluids n, first-order resistance A and second-order resistance B, using a fourth-order Runge–Kutta integration scheme with the Newton–Raphson shooting method. Velocity and shear stress at the body surface are presented for a range of the above parameters. These results are also compared with the corresponding flow problems for a Newtonian fluid. Numerical results show that for the case of the constant wedge angle and material parameter A, the local skin friction coefficient is lower for a dilatant fluid as compared with the pseudo-plastic or Newtonian fluids.     

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     

A boundary integral method for two-dimensional (non)-Newtonian drops in slow viscous flow

A boundary integral method for the simulation of the time-dependent deformation of Newtonian or non-Newtonian drops suspended in a Newtonian fluid is developed. The boundary integral formulation for Stokes flow is used and the non-Newtonian stress is treated as a source term which yields an extra integral over the domain of the drop. The implementation of the boundary conditions is facilitated by rewriting the domain integral by means of the Gauss divergence theorem. To apply the divergence theorem smoothness assumptions are made concerning the non-Newtonian stress tensor. The correctness of these assumptions in actual simulations is checked with a numerical validation procedure. The method appears mathematically correct and the numerical algorithm is second order accurate. Besides this validation we present simulation results for a Newtonian drop and a drop consisting of an Oldroyd-B fluid. The results for Newtonian and non-Newtonian drops in two dimensions indicate that the steady state deformation is quite independent of the drop-fluid. The deformation process, however, appears to be strongly dependent on the drop-fluid. For the non-Newtonian drop a mechanical model is developed to describe the time-dependent deformation of the cylinder for small capillary numbers.     

Numerical study of electroosmotic micromixing of non-Newtonian fluids

Biofluids which exhibit non-Newtonian behavior are widely used in microfluidic devices which involve fluid mixing in microscales. In order to study the effects of shear depending viscosity of non-Newtonian fluids on characteristics of electroosmotic micromixing, a numerical investigation of flow of power-law fluid in a two-dimensional microchannel with nonuniform zeta potential distributions along the channel walls was carried out via finite volume scheme. The simulation results confirmed that the shear depending viscosity has a significant effect on the degree of mixing efficiency. It was shown as the fluid behavior index of power-law fluid, n, decreases, more homogeneous solution can be achieved at the microchannel outlet. Hence, electroosmotic micromixing was found more practical and efficient in microscale mixing of pseudoplastic fluids rather than those Newtonian and dilatant ones. Furthermore, it was found that increase in Reynolds number results in lower mixing efficiency while electroosmotic forces are kept constant.