Local similarity solution for the flow of a “second-grade” viscoelastic fluid above a moving plate

The boundary layer flow of a viscoelastic fluid of the second-grade type over a rigid continuous plate moving through an otherwise quiescent fluid with constant velocity U is studied. Assuming the flow to be laminar and two-dimensional, local similarity solution is found with fluid's elasticity and plate's withdrawal speed as the main variables. Results are presented for velocity profiles, boundary layer thickness, wall skin friction coefficient and fluid entrainment in terms of the local Deborah number. A marked formation of boundary layer is predicted, even at low Reynolds numbers, provided the Deborah number is sufficiently large. The boundary layer thickness and the wall skin friction coefficient are found to scale with fluid's elasticity—both decreasing the higher the fluid's elasticity. The amount of fluid entrained is also predicted to decrease whenever a fluid exhibits elastic behavior.     

Melting heat transfer in an axisymmetric stagnation-point flow of the Jeffrey fluid

This investigation explores the characteristics of melting heat transfer in a boundary layer flow of the Jeffrey fluid near the stagnation point on a stretching sheet subject to an applied magnetic field. The governing boundary layer equations are transformed to ordinary differential equations by similarity transformations. Resulting nonlinear problems are solved analytically by the homotopy analysis method. It is noticed that an increase in the melting parameter decreases the dimensionless velocity and temperature, while an increase in the Deborah number increases the velocity and momentum boundary layer thickness.     

MHD flow and heat transfer for the upper-convected Maxwell fluid over a stretching sheet

In the present work, the effect of MHD flow and heat transfer within a boundary layer flow on an upper-convected Maxwell (UCM) fluid over a stretching sheet is examined. The governing boundary layer equations of motion and heat transfer are non-dimensionalized using suitable similarity variables and the resulting transformed, ordinary differential equations are then solved numerically by shooting technique with fourth order Runge–Kutta method. For a UCM fluid, a thinning of the boundary layer and a drop in wall skin friction coefficient is predicted to occur for higher the elastic number. The objective of the present work is to investigate the effect of Maxwell parameter β, magnetic parameter Mn and Prandtl number Pr on the temperature field above the sheet.     

A numerical study for three-dimensional viscoelastic flow inspired by non-linear radiative heat flux

Present article examines the three-dimensional flow of upper-convected Maxwell (UCM) fluid over a radiative bi-directional stretching surface. Novel non-linear Rosseland formula for thermal radiation is utilized in the formulation of energy equation. The conventional transformations lead to a strongly non-linear differential system which is treated numerically through Runge–Kutta integration procedure together with the shooting approach. We found that heat transfer rate from the sheet has inverse as well as non-linear relationship with wall to ambient temperature ratio. Moreover an increase in viscoelastic fluid parameter (Deborah number) corresponds to a decrease in the fluid velocity and the boundary layer thickness.     

Melting heat transfer in the stagnation‐point flow of an upper‐convected Maxwell (UCM) fluid past a stretching sheet

The steady laminar boundary layer flow and heat transfer past a stretching sheet arre considered. Upper‐convected Maxwell (UCM) fluid is treated as a rheological model. The resulting nonlinear differential system is solved by homotopy analysis method (HAM). The influence of melting parameter (M), Prandtl number (Pr), Deborah number (β) and stretching ratio (A = a/c) on the velocity and temperature profiles is thoroughly examined. It is noticed that fields are effected appreciably with the variation of parameters. Furthermore, it is seen that the local Nusselt number is a decreasing function of melting parameter. Copyright © 2011 John Wiley & Sons, Ltd.     

The behavior of viscoelastic squeeze films subject to normal oscillations, including the effect of fluid inertia

The problem of the squeeze film flow of a viscoelastic fluid between parallel, circular disks is analyzed. The upper disk is subject to small, axial oscillations. Lodge's “rubber-like liquid” is used as the viscoelastic fluid model, and fluid inertia forces are included. An exact solution to the equations of motion is obtained involving in-phase and out-of-phase components of velocity field and load, with respect to the plate velocity. Peculiar resonance phenomena in the load amplitude are exhibited at high Deborah number. At certain combinations of Reynolds number and Deborah number, the in-phase and/or out-of-phase velocity field components may attain an unusual circulating type of motion in which the flow reverses direction across the film. In the low Deborah number limit, and in the low Reynolds number limit, the results of this study reduce to those obtained by other workers.     

Viscoelastic flow between eccentric rotating cylinders: unstructured control volume method

This paper reports a convergent numerical algorithm for the Upper-Convected Maxwell (UCM) fluid between two eccentric cylinders at various eccentricity ratios (?); the outer cylinder is stationary, and the inner one rotating. The problem is solved by an unstructured control volume method (UCV), which is designed for a general viscoelastic flow problem with an arbitrary computational domain. A self-consistent false diffusion technique and an iteration scheme are used in combination to solve the problem. The computations of the UCM fluid using the numerical algorithm are carried out to a higher value of the Deborah number (De) at each eccentricity tested than hitherto possible with previous numerical simulations. The solutions are compared with previous numerical results, confirming the effectiveness of the UCV method as a general technique for solving viscoelastic flow problems.     

PIV measurements of slow flow of a viscoelastic fluid within a porous medium

Particle image velocimetry (PIV) and pressure loss measurements were used to investigate slow flow through a square array of cylinders having a solid fraction of 10%. The test fluids were a Newtonian fluid and a Boger fluid, both of high viscosity such that the Reynolds number did not exceed 0.1. The pressure loss data reveal that the onset of elastic effects occurred at a Deborah number around 0.5 and that flow resistance was up to several times Newtonian values at Deborah numbers up to 3. PIV showed that the transverse velocity profiles for the Newtonian and non-Newtonian fluid were the same at Deborah numbers below onset. Above onset, the profiles became skewed, increasingly so as the Deborah number increased. In the wake regions between cylinders in a column, periodic flow structures formed in the spanwise direction. The structures were staggered from column to column, consistent with the skewing and were offset. These flow patterns are the result of an apparent elastic instability.     

Stress boundary layers in the viscoelastic flow past a cylinder in a channel： limiting solutions

The flow past a cylinder in a channel with the aspect ratio of 2:1 for the upper convected Maxwell (UCM) fluid and the Oldroyd-B fluid with the viscosity ratio of 0.59 is studied by using the Galerkin/Least-square finite element method and a p-adaptive refinement algorithm. A posteriori error estimation indicates that the stress-gradient error dominates the total error. As the Deborah number, D_e, approaches 0.8 for the UCM fluid and 0.9 for the Oldroyd-B fluid, strong stress boundary layers near the rear stagnation point are forming, which are characterized by jumps of the stress-profiles on the cylinder wall and plane of symmetry, huge stress gradients and rapid decay of the gradients across narrow thicknesses. The origin of the huge stress-gradients can be traced to the purely elongational flow behind the rear stagnation point, where the position at which the elongation rate is of 1/2D_e approaches the rear stagnation point as the Deborah number approaches the critical values. These observations imply that the cylinder problem for the UCM and Oldroyd-B fluids may have physical limiting Deborah numbers of 0.8 and 0.9, respectively.The project supported by the National Natural Science Foundation of China (50335010 and 20274041) and the MOLDFLOW Comp. Australia.     

Transport Phenomenon in a Third-Grade Fluid Over an Oscillating Surface

The heat and mass transfer effects on the flow of a conducting third-grade fluid over an oscillating vertical porous plate with chemical reactions are considered. Highly nonlinear governing equations of the third-grade fluid are solved analytically by using a multi-parameter perturbation technique and compared with the numerical results obtained by the parallel shooting method. The fluid flow velocity, temperature, and concentration are analyzed as functions of the Hartmann number, suction parameter, Prandtl and Schmidt numbers, and chemical reaction parameter.     

The effect of an additional boundary condition on the plane creeping flow of a second-order fluid

A similarity solution is obtained for the two-dimensional creeping flow of a second-order fluid with non-parallel porous walls. The resulting ordinary differential equation is fifth order. Thus, an additional velocity boundary condition is needed, the other four being due to the usual no-slip conditions. Having chosen to prescribe the rate of shear at the wall, the problem is solved by a standard numerical routine. A singular perturbation analysis is developed for small values of the Deborah number. A type of boundary layer forms for which the viscous Newtonian case is the outer solution.     

Stagnation-point flow of upper-convected Maxwell fluids

Two-dimensional stagnation-point flow of viscoelastic fluids is studied theoretically assuming that the fluid obeys the upper-convected Maxwell (UCM) model. Boundary-layer theory is used to simplify the equations of motion which are further reduced to a single non-linear third-order ODE using the concept of stream function coupled with the technique of the similarity solution. The equation so obtained was solved using Chebyshev pseudo-spectral collocation-point method. Based on the results obtained in the present work, it is concluded that the well-established but controversial prediction that in stagnation-point flows of viscoelastic fluids the velocity inside the boundary layer may exceed that outside the layer may just be an artifact of the rheological model used in previous studies (namely, the second-grade model). No such peculiarity is predicted to exist for the Maxwell model. For a UCM fluid, a thickening of the boundary layer and a drop in wall skin friction coefficient is predicted to occur the higher the elasticity number. These predictions are in direct contradiction with those reported in the literature for a second-grade fluid.     

Visco-elastic fluid flow past an infinite vertical porous plate in the presence of first-order chemical reaction

An analysis has been developed to study the unsteady free convection flow of an incompressible visco-elastic fluid on a continuously moving vertical porous plate in the presence of a first-order chemical reaction. The governing equations are solved numerically using an implicit finite difference technique. The obtained numerical solutions are compared with the analytical solutions. The velocity profiles are presented. A parametric analysis is performed to illustrate the influences of the visco-elastic parameter, the dimensionless chemical reaction parameter, and the plate moving velocity on the steady state velocity profiles, the time dependent friction coefficient, the Nusselt number, and the Sherwood number.     

Effects of thermophoresis and radiation on laminar flow along a semi-infinite vertical plate

The present contribution deals with the thermophoresis particle deposition and thermal radiation effects on the flow, heat and mass transfer characteristics in a viscous fluid over a semi-infinite vertical porous plate. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed coupled nonlinear ordinary differential equations are solved numerically by means of the fourth-order Runge–Kutta method with a shooting technique. The effects of different parameters on the dimensionless velocity, temperature, and concentration profiles are shown graphically. In addition, results for the local skin-friction coefficient, the local Nusselt number, and the local Sherwood number are tabulated and discussed.     

Perturbation theory for viscoelastic fluids between eccentric rotating cylinders

The circumferential and radial profiles of velocity, pressure and stress are derived for the flow of model viscoelastic liquids between two slightly eccentric cylinders with the inner one rotating. Singular perturbation methods are used to derive expansions valid for small gaps between the cylinders, but for all Deborah numbers. Results for Newtonian, second-order, Criminale-Ericksen-Filbey, upper-convected Maxwell, and White-Metzner constitutive equation separate the effects of elasticity, memory, and shear thinning on the development of the large stress gradients that hinder numerical solutions with these models in more complicated geometries. The effect of the constitutive equation on the critical Deborah number for flow separation is presented.     

Predictions of the excess pressure drop of Boger fluids through a 2:1:2 contraction–expansion geometry using non-equilibrium molecular dynamics

Non-equilibrium molecular dynamics are used to generate the flow of polymer solutions, specifically of Boger fluids, through a planar 2:1:2 contraction–expansion geometry. The solvent molecules are represented by Lennard–Jones particles, while linear molecules are described by spring-monomers with a finite extensible non-linear elastic spring potential. The equations for Poiseuille flow are solved using a multiple time-scale algorithm extended to non-equilibrium situations. Simulations are performed at constant temperature using Nose–Hoover dynamics. At simulation conditions, changes in concentration show no significant effect on molecular conformation, velocity profiles, and stress fields, while variations in the Deborah number have a strong influence on fluid response. Increasing the magnitude of the Deborah number (De), larger deformation rates are developed in the flow region. For a Deborah number of one, the non-dimensional pressure drop presents values lower than the correspondent Newtonian case. However, for large Deborah numbers, the pressure drop increases above the Newtonian reference. An effective excess pressure drop above the Newtonian value is predicted for Boger fluids along this geometry.     

Boundary layer flow of Maxwell fluid due to torsional motion of cylinder: modeling and simulation

This paper investigates the boundary layer flow of the Maxwell fluid around a stretchable horizontal rotating cylinder under the influence of a transverse magnetic field. The constitutive flow equations for the current physical problem are modeled and analyzed for the first time in the literature. The torsional motion of the cylinder is considered with the constant azimuthal velocity E. The partial differential equations (PDEs) governing the torsional motion of the Maxwell fluid together with energy transport are simplified with the boundary layer concept. The current analysis is valid only for a certain range of the positive Reynolds numbers. However, for very large Reynolds numbers, the flow becomes turbulent. Thus, the governing similarity equations are simplified through suitable transformations for the analysis of the large Reynolds numbers. The numerical simulations for the flow, heat, and mass transport phenomena are carried out in view of the bvp4c scheme in MATLAB. The outcomes reveal that the velocity decays exponentially faster and reduces for higher values of the Reynolds numbers and the flow penetrates shallower into the free stream fluid. It is also noted that the phenomenon of stress relaxation, described by the Deborah number, causes to decline the flow fields and enhance the thermal and solutal energy transport during the fluid motion. The penetration depth decreases for the transport of heat and mass in the fluid with the higher Reynolds numbers. An excellent validation of the numerical results is assured through tabular data with the existing literature.     

Purely elastic flow asymmetries in flow-focusing devices

The flow of a viscoelastic fluid through a microfluidic flow-focusing device is investigated numerically with a finite-volume code using the upper-convected Maxwell (UCM) and Phan-Thien–Tanner (PTT) models. The conceived device is shaped much like a conventional planar “cross-slot” except for comprising three inlets and one exit arm. Strong viscoelastic effects are observed as a consequence of the high deformation rates. In fact, purely elastic instabilities that are entirely absent in the corresponding Newtonian fluid flow are seen to occur as the Deborah number (De) is increased above a critical threshold. From two-dimensional numerical simulations we are able to distinguish two types of instability, one in which the flow becomes asymmetric but remains steady, and a subsequent instability at higher De in which the flow becomes unsteady, oscillating in time. For the UCM model, the effects of the geometric parameters of the device (e.g. the relative width of the entrance branches, WR) and of the ratio of inlet average velocities (VR) on the onset of asymmetry are systematically examined. We observe that for high velocity ratios, the critical Deborah number is independent of VR (e.g. De_c ≈ 0.33 for WR = 1), but depends non-monotonically on the relative width of the entrance branches. Using the PTT model we are able to demonstrate that the extensional viscosity and the corresponding very large stresses are decisive for the onset of the steady-flow asymmetry.     

Chemical reaction effects on unsteady MHD flow past semi-infinite vertical porous plate with viscous dissipation

The chemical reaction effect on an unsteady magnetohydrodynamic (MHD) flow past a semi-infinite vertical porous plate with viscous dissipation is analyzed. The governing equations of motion, energy, and species are transformed into ordinary differential equations (ODEs) using the time dependent similarity parameter. The resultant ODEs are then solved numerically by a finite element method. The effects of various parameters on the velocity, temperature, and concentration profiles are presented graphically, and the values of the skin-friction, Nusselt number, and Sherwood number for various values of physical parameters are presented through tables.