On study of horizontal thin film flow of Sisko fluid due to surface tension gradient

The present paper is concerned with the steady thin film flow of the Sisko fluid on a horizontal moving plate, where the surface tension gradient is a driving mechanism. The analytic solution for the resulting nonlinear ordinary differential equation is obtained by the Adomian decomposition method (ADM). The physical quantities are derived including the pressure profile, the velocity profile, the maximum residue time, the stationary points, the volume flow rate, the average film velocity, the uniform film thickness, the shear stress, the surface tension profile, and the vorticity vector. It is found that the velocity of the Sisko fluid film decreases when the fluid behavior index and the Sisko fluid parameter increase, whereas it increases with an increase in the inverse capillary number. An increase in the inverse capillary number results in an increase in the surface tension which in turn results in an increase in the surface tension gradient on the Sisko fluid film. The locations of the stationary points are shifted towards the moving plate with the increase in the inverse capillary number, and vice versa locations for the stationary points are found with the increasing Sisko fluid parameter. Furthermore, shear thinning and shear thickening characteristics of the Sisko fluid are discussed. A comparison is made between the Sisko fluid film and the Newtonian fluid film.     

Steady flow and heat transfer of a magnetohydrodynamic Sisko fluid through porous medium in annular pipe

In this paper, the steady flow and heat transfer of a magnetohydrodynamic fluid is studied. The fluid is assumed to be electrically conducting in the presence of a uniform magnetic field and occupies the porous space in annular pipe. The governing nonlinear equations are modeled by introducing the modified Darcy's law obeying the Sisko model. The system is solved using the homotopy analysis method (HAM), which yields analytical solutions in the form of a rapidly convergent infinite series. Also, HAM is used to obtain analytical solutions of the problem for noninteger values of the power index. The resulting problem for velocity field is then numerically solved using an iterative method to show the accuracy of the analytic solutions. The obtained solutions for the velocity and temperature fields are graphically sketched and the salient features of these solutions are discussed for various values of the power index parameter. We also present a comparison between Sisko and Newtonian fluids. Copyright © 2011 John Wiley & Sons, Ltd.     

Thin film flow of non‐Newtonian MHD fluid on a vertically moving belt

Thin film flow of an Oldroyd 6‐constant fluid on a vertical moving belt is investigated analytically and numerically. The governing equations for the flow field are derived for a steady one‐dimensional flow. The effect of constant applied magnetic field is included and its influence on the flow field is studied. The nonlinear governing equations are solved analytically and the exact solution is obtained in an elegant way. Numerical solutions are also obtained using higher‐order Chebyshev spectral methods. The influence of various non‐Newtonian parameters, gravitational force and applied magnetic field is investigated. The results showing the effect of gravity, magnetic field and material constants α₁ and α₂ are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical and analytical solutions of an Oldroyd 8‐constant MHD fluid with nonlinear slip conditions

The flow of an Oldroyd 8‐constant non‐Newtonian MHD fluid is investigated analytically and numerically. The governing equations for the flow field are derived for a steady one‐dimensional flow. The effect of constant applied magnetic field is included and its influence on the flow field is studied. The nonlinear governing equation along with nonlinear boundary conditions is solved analytically and the solution is obtained in an elegant way. Numerical solutions are also obtained using higher order Chebyshev spectral methods. The influence of various non‐Newtonian parameters and applied magnetic field is investigated. Results showing the effect of various physical parameters of the flow are presented and investigated. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Stagnation flows for the Oldroyd-B fluid

Exact solutions to the plane and axi-symmetric stagnation flows of an Oldroyd-B fluid are reported. It is found that a steady flow is possible if the Weissenberg numberWi, defined by the product of the Maxwellian relaxation time and the shear rate at infinity, satisfies – 1/2 <Wi < 1/m, wherem = 1 in an axisym-metric flow andm = 2 in a plane flow. Furthermore, the fluid elasticity always decreases the boundary-layer thickness. An Oldroyd-B fluid with the parameters matched those of a typical Boger fluid behaves essentially like a Newtonian fluid in a stagnation flow.     

The micropolar fluid model for blood flow through a tapered artery with a stenosis

A micropolar model for axisymmetric blood flow through an axially nonsymmetreic but radially symmetric mild stenosis tapered artery is presented. To estimate the effect of the stenosis shape, a suitable geometry has been considered such that the axial shape of the stenosis can be changed easily just by varying a parameter （referred to as the shape parameter）. The model is also used to study the effect of the taper angle Ф. Flow parameters such as the velocity, the resistance to flow （the resistance impedance）, the wall shear stress distribution in the stenotic region and its magnitude at the maximum height of the stenosis （stenosis throat） have been computed for different values of the shape parameter n, the taper angle Ф, the coupling number N and the micropolar parameter m. It is shown that the resistance to flow decreases with increasing the shape parameter n and the micropolar parameter m while it increases with increasing the coupling number N. So, the magnitude of the resistance impedance is higher for a micropolar fluid than that for a Newtonian fluid model. Finally, the velocity profile, the wall shear stress distribution in the stenotic region and its magnitude at the maximum height of the stenosis are discussed for different values of the parameters involved on the problem.     

On stagnation point flow of Sisko fluid over a stretching sheet

The steady two-dimensional stagnation-point flow, represented by Sisko fluid constitutive model, over a stretching sheet is investigated theoretically. Using suitable similarity transformations, the governing boundary-layer equations are transformed into the self-similar non-linear ordinary differential equation. The transformed equation is then solved using a very efficient analytic technique namely the homotopy analysis method (HAM) and the HAM solutions are validated by the exact analytic solutions obtain in certain special cases. The influence of the power-law index (n), the material parameter (A) and the velocity ratio parameter (d/c) on the flow characteristics is studied and presented through several graphs. In addition, the local skin friction coefficient for several values of these parameters is tabulated and examined. The similarity solutions for both the Newtonian and the power-law fluids are presented as special cases of the analysis. The results obtained reveal that, in comparison with the Newtonian and the power-law fluids, the velocity profiles of the Sisko fluid are much faster (slower), for d/c<1 (d/c>1), respectively.     

Non‐Newtonian blood flow dynamics in a right internal carotid artery with a saccular aneurysm

Flow dynamics plays an important role in the pathogenesis and treatment of cerebral aneurysms. The temporal and spatial variations of wall shear stress in the aneurysm are hypothesized to be correlated with its growth and rupture. In addition, the assessment of the velocity field in the aneurysm dome and neck is important for the correct placement of endovascular coils. This work describes the flow dynamics in a patient‐specific model of carotid artery with a saccular aneurysm under Newtonian and non‐Newtonian fluid assumptions. The model was obtained from three‐dimensional rotational angiography image data and blood flow dynamics was studied under physiologically representative waveform of inflow. The three‐dimensional continuity and momentum equations for incompressible and unsteady laminar flow were solved with a commercial software using non‐structured fine grid with 283 115 tetrahedral elements. The intra‐aneurysmal flow shows complex vortex structure that change during one pulsatile cycle. The effect of the non‐Newtonian properties of blood on the wall shear stress was important only in the arterial regions with high velocity gradients, on the aneurysmal wall the predictions with the Newtonian and non‐Newtonian blood models were similar. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical simulation of MHD stagnation point flow and heat transfer of a micropolar fluid towards a heated shrinking sheet

A comprehensive study of magneto hydrodynamics two‐dimensional stagnation flow with heat transfer characteristics towards a heated shrinking sheet immersed in an electrically conducting incompressible micropolar fluid in the presence of a transverse magnetic field is analyzed numerically. The governing continuity, momentum, angular momentum and heat equations together with the associated boundary conditions are first reduced to a set of self similar nonlinear ordinary differential equations using a similarity transformation and are then solved by a method based on finite difference discretization. Some important features of the flow and heat transfer in terms of normal and streamwise velocities, microrotation and temperature distributions for different values of the governing parameters are analyzed, discussed and presented through tables and graphs. The results indicate that the reverse flow caused due to shrinking of the sheet can be stopped by applying a strong magnetic field. The magnetic field enhances the shear stresses and decreases the thermal boundary layer thickness. The heat loss per unit area from the sheet decreases with an increase in the shrinking parameter. Micropolar fluids exhibit reduction in shear stresses and heat transfer rate as compared with Newtonian fluids, which may be beneficial in the flow and thermal control of polymeric processing. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

MHD flow and heat transfer of micropolar fluid between two porous disks

A numerical study is carried out for the axisymmetric steady laminar incompressible flow of an electrically conducting micropolar fluid between two infinite parallel porous disks with the constant uniform injection through the surface of the disks. The fluid is subjected to an external transverse magnetic field. The governing nonlinear equations of motion are transformed into a dimensionless form through von Karman’s similarity transformation. An algorithm based on a finite difference scheme is used to solve the reduced coupled ordinary differential equations under associated boundary conditions. The effects of the Reynolds number, the magnetic parameter, the micropolar parameter, and the Prandtl number on the flow velocity and temperature distributions are discussed. The results agree well with those of the previously published work for special cases. The investigation predicts that the heat transfer rate at the surfaces of the disks increases with the increases in the Reynolds number, the magnetic parameter, and the Prandtl number. The shear stresses decrease with the increase in the injection while increase with the increase in the applied magnetic field. The shear stress factor is lower for micropolar fluids than for Newtonian fluids, which may be beneficial in the flow and thermal control in the polymeric processing.     

A spectral/hp least‐squares finite element analysis of the Carreau–Yasuda fluids

A least‐squares finite element model with spectral/hp approximations was developed for steady, two‐dimensional flows of non‐Newtonian fluids obeying the Carreau–Yasuda constitutive model. The finite element model consists of velocity, pressure, and stress fields as independent variables (hence, called a mixed model). Least‐squares models offer an alternative variational setting to the conventional weak‐form Galerkin models for the Navier–Stokes equations, and no compatibility conditions on the approximation spaces used for the velocity, pressure, and stress fields are necessary when the polynomial order (p) used is sufficiently high (say, p > 3, as determined numerically). Also, the use of the spectral/hp elements in conjunction with the least‐squares formulation with high p alleviates various forms of locking, which often appear in low‐order least‐squares finite element models for incompressible viscous fluids, and accurate results can be obtained with exponential convergence. To verify and validate, benchmark problems of Kovasznay flow, backward‐facing step flow, and lid‐driven square cavity flow are used. Then the effect of different parameters of the Carreau–Yasuda constitutive model on the flow characteristics is studied parametrically. Copyright © 2016 John Wiley & Sons, Ltd.     

Analytic solution for flow of Sisko fluid through a porous medium

This paper presents the analytic solution for flow of a magnetohydrodynamic (MHD) Sisko fluid through a porous medium. The non-linear flow problem in a porous medium is formulated by introducing the modified Darcy’s law for Sisko fluid to discuss the flow in a porous medium. The analytic solutions are obtained using homotopy analysis method (HAM). The obtained analytic solutions are explicitly expressed by the recurrence relations and can give results for all the appropriate values of material parameters of the examined fluid. Moreover, the well-known solutions for a Newtonian fluid in non-porous and porous medium are the limiting cases of our solutions.     

Further contributions on the flow past a stationary and confined cylinder: Creeping and slowly moving flow of Power law fluids

The steady flow of generalized Newtonian fluid around a stationary cylinder placed between two parallel plates was studied numerically. Finite volume method was applied to solve the momentum equations along with the continuity equation and the Power law rheological model within the laminar flow regime for a range of the Reynolds number Re and the Power law index n values. The values of the Reynolds number, based on physical and rheological properties, cylinder radius and bulk velocity, were varied between 0.0001≤Re≤10, while the Power law index values mapped the 0.50≤n≤1.50 range, allowing for the investigation of both shear-thinning and shear-thickening effects at the creeping as well as slowly moving fluid flow conditions. We report accurate results of a systematic study with a focus on the most important characteristics of fluid flow past circular cylinder. It is shown that for the creeping flow regime there exist finite sized redevelopment length, drag and loss coefficient. Last but not least, the present numerical results indicate that the shear-thinning viscous behaviour decreases the onset of flow separation.     

Slip effects on heat and mass transfer in MHD micropolar fluid flow over an inclined plate with thermal radiation and chemical reaction

The influence of partial slip, thermal radiation, chemical reaction and temperature‐dependent fluid properties on heat and mass transfer in hydro‐magnetic micropolar fluid flow over an inclined permeable plate with constant heat flux and non‐uniform heat source/sink is studied. The transverse magnetic field is assumed as a function of the distance from the origin. Also it is assumed that the fluid viscosity and the thermal conductivity vary as an inverse function and linear function of temperature, respectively. With the use of the similarity transformation, the governing system of non‐linear partial differential equations are transformed into non‐linear ordinary differential equations and are solved numerically using symbolic software MATHEMATICA 7.0 (Wolfram Research, Champaign, IL). The numerical values obtained for the velocity, microrotation, temperature, species concentration, skin friction coefficient and the Nusselt number are presented through graphs and tables for several sets of values of the parameters. The effects of various physical parameters on the flow and heat transfer characteristics are discussed.Copyright © 2011 John Wiley & Sons, Ltd.     

Peristaltic transport of a Herschel-Bulkley fluid in an inclined tube

Peristaltic flow of Herschel-Bulkley fluid in an inclined tube is analyzed. The velocity distribution, the stream function and the volume flow rate are obtained. Also, when the yield stress ratio τ→0, and when the inclination parameter α=0 and the fluid parameter n=1, the results agree with those of Jaffrin and Shapiro (Ann. Rev. Fluid Mech. 3 (1971) 13) for peristaltic transport of a Newtonian fluid in a horizontal tube. The effects of τ and n on the pressure drop and the mean flow are discussed through graphs. Furthermore, the results for the peristaltic transport of Bingham and power law fluids through a flexible tube are obtained and discussed. The results obtained for the flow characteristics reveal many interesting behaviors that warrant further study of the effects of Herschel-Bulkley fluid on the flow characteristics.     

Self-similar solution of the unsteady flow in the stagnation point region of a rotating sphere with a magnetic field

The unsteady flow and heat transfer of a viscous incompressible electrically conducting fluid in the forward stagnation point region of a rotating sphere in the presence of a magnetic field are investigated in this study. The unsteadiness in the flow field is caused by the velocity at the edge of the boundary layer and the angular velocity of the rotating sphere, both varying continuously with time. The system of ordinary differential equations governing the flow is solved numerically. For some particular cases, an analytical solution is also obtained. It is found that the surface shear stresses in x- and y-directions and the surface heat transfer increase with the acceleration, the magnetic and the rotation parameters whether the magnetic field is fixed relative to the fluid or body, except that the surface shear stress in x-direction and the surface heat transfer decrease with increasing the magnetic parameter when the magnetic field is fixed relative to the body. For a certain value of the acceleration parameter, the surface shear stress in the x-direction vanishes while the surface shear stress in the y-direction and the surface heat transfer remain finite. Also, below a certain value of the acceleration parameter, reverse flow occurs in the x-component of the velocity profile. Received on 18 May 1998     

Monolithic Newton‐multigrid solution techniques for incompressible nonlinear flow models

We present special Newton‐multigrid techniques for stationary incompressible nonlinear flow models discretized by the high order LBB‐stable Q₂P₁ element pair. We treat the resulting nonlinear and the corresponding linear discrete systems by a fully coupled monolithic approach to maintain high accuracy and robustness, particularly with respect to different rheological behaviors and also regarding different problem sizes and types of nonlinearity. Here, local pressure Schur complement techniques are presented as a generalization of the classical Vanka smoother. The discussed methodology is implemented for the well‐known flow around cylinder benchmark configuration for generalized Newtonian as well as non‐Newtonian flows including non‐isothermal, shear/pressure dependent and viscoelastic effects.Copyright © 2012 John Wiley & Sons, Ltd.