Numerical study on viscoelastic fluid flow past a rigid body

Viscoelastic non-Newtonian fluids can be achieved by adding a small amount of polymer additives to a Newtonian fluid. In this paper, numerical simulations are used to investigate the influence of such polymer additives on the behavior of flow past a circular cylinder. A numerical method is proposed that discretizes the non-linear viscoelastic system on a uniform Cartesian grid, with a penalization method to model the presence of the cylinder. The drag of the cylinder and the flow behavior under the effect of different Reynolds numbers ($\mathrm{Re}$), Weissenberg numbers (Wi) and polymer viscosity ratios (ε) are studied. Numerical results show that different flow characteristics are exhibited in different parameter zones. The polymer viscosity ratio plays an important role at low Weissenberg and Reynolds numbers, but as the Reynolds and Weissenberg numbers increase, the influence of ε weakens. The drag force of the cylinder is mostly affected by the Reynolds and Weissenberg numbers. At low Reynolds numbers, the drag of the cylinder and the flow fields are only affected by a large value of Wi when the elastic forces are strong. Non-trivial drag reduction occurs only when there is vortex shedding in the wake flow, whereas drag enhancement happens when the vortex shedding is inhibited.     

Hydromagnetic linear instability analysis of Giesekus fluids in plane Poiseuille flow

The effects of a fluid’s elasticity are investigated on the instability of plane Poiseuille flow on the presence of a transverse magnetic field. To determine the critical Reynolds number as a function of the Weissenberg number, a two-dimensional linear temporal stability analysis will be used assuming that the viscoelastic fluid obeys Giesekus model as its constitutive equation. Neglecting terms nonlinear in the perturbation quantities, an eigenvalue problem is obtained which is solved numerically by using the Chebyshev collocation method. Based on the results obtained in this work, fluid’s elasticity is predicted to have a stabilizing or destabilizing effect depending on the Weissenberg number being smaller or larger than one. Similarly, solvent viscosity and also the mobility factor are both found to have a stabilizing or destabilizing effect depending on their magnitude being smaller or larger than a critical value. In contrast, the effect of the magnetic field is predicted to be always stabilizing.     

On the non-orthogonal Stagnation flow of the Oldroyd-B fluid

This paper is concerned with a non-orthogonal stagnation flow of an Oldroyd-B fluid between two parallel plates. We reduce the problem to a set of ordinary differential equations (ODE's), which is then solved with finite differences using a parameter continuation method. Perturbation analyses are also carried out for small Reynolds numbers and small Weissenberg numbers respectively. The solution of the set of ODE's is discussed. It is known that for a Newtonian fluid, the stagnation point shifts from the potential flow case in the opposite direction of the tangential velocity. The effect of the fluid elasticity is to reduce this shift. It is also shown that the Oldroyd-B model has a limiting Weissenbeg number, depending on the angle of the injected flow.     

Computation of viscoelastic fluid flows using continuation methods

The numerical simulation of viscoelastic fluid flow becomes more difficult as a physical parameter, the Weissenberg number, increases. Specifically, at a Weissenberg number larger than a critical value, the iterative nonlinear solver fails to converge, a phenomenon known as the high Weissenberg number problem. In this work we describe the application and implementation of continuation methods to the nonlinear Johnson–Segalman model for steady-state viscoelastic flows. Simple, natural, and pseudo-arclength continuation approaches in Weissenberg number are investigated for a discontinuous Galerkin finite element discretization of the equations. Computations are performed for a benchmark contraction flow and, several aspects of the performance of the continuation methods including high Weissenberg number limits, are discussed.     

Sliding and squeezing flow of a viscoelastic fluid in a wedge

The paper reports an exact three-dimensional similarity solution for the Oldroyd fluid B. The flow involved is generated by closing as well as sliding the boundaries of a two-dimensional wedge. It is found that the squeezing motion is independent of the sliding motion, but not vice-versa. The squeezing load is shown to be a decreasing function of the Weissenberg number, while the frictional coefficient is only weakly-dependent on the Weissenberg number.     

A Projection Preconditioner for Solving the Implicit Immersed Boundary Equations

This paper presents a method for solving the linear semi-implicit immersed boundary equations which avoids the severe time step restriction presented by explicit-time methods. The Lagrangian variables are eliminated via a Schur complement to form a purely Eulerian saddle point system, which is preconditioned by a projection operator and then solved by a Krylov subspace method. From the viewpoint of projection methods, we derive an ideal preconditioner for the saddle point problem and compare the efficiency of a number of simpler preconditioners that approximate this perfect one. For low Reynolds number and high stiffness, one particular projection preconditioner yields an efficiency improvement of the explicit IB method by a factor around thirty. Substantial speed-ups over explicit-time method are achieved for Reynolds number below 100. This speedup increases as the Eulerian grid size and/or the Reynolds number are further reduced.     

Viscoelastic flow in a curved channel: A similarity solution for the Oldroyd-B fluid

A similarity solution is used to analyse the flow of the Oldroyd fluid B, which includes the Newtonian and Maxwell fluids, in a curved channel modelled by the narrow annular region between two circular concentric cylinders of large radius. The solution is exact, including inertial forces. It is found that the non-Netonian kinematics are very similar to the Newtonian ones, although some stress components can become very large. At high Reynolds number a boundary layer is developed at the inner cylinder. The structure of this boundary layer is asymptotically analysed for the Newtonian fluid. Non-Newtonian stress boundary layers are also developed at the inner cylinder at large Reynolds numbers.     

The flow of a viscoelastic fluid in a wedge

Summary It is shown that the kinematics of the flow of a general viscoelastic fluid in a wedge, one plate of which is being stretched at a rate proportional to the distance from the wedge apex, is Newtonian in character. Existence proof is given when non-Newtonian effects are slight. Furthermore, the stress field is multivalued at the wedge apex and the pressure field is logarithmically singular there. The strength of this singularity increases with the Weissenberg number.     

Fully-developed free-convective flow of micropolar and viscous fluids in a vertical channel

The problem of fully-developed laminar free-convection flow in a vertical channel is studied analytically with one region filled with micropolar fluid and the other region with a viscous fluid. Using the boundary and interface conditions proposed by previous investigators, analytical expressions for linear velocity, micro-rotation velocity and temperature have been obtained. Numerical results are presented graphically for the distribution of velocity, micro-rotation velocity and temperature fields for varying physical parameters such as the ratio of Grashof number to Reynolds number, viscosity ratio, width ratio, conductivity ratio and micropolar fluid material parameter. It is found that the effect of the micropolar fluid material parameter suppress the velocity whereas it enhances the micro-rotation velocity. The effect of the ratio of Grashof number to Reynolds number is found to enhance both the linear velocity and the micro-rotation velocity. The effects of the width ratio and the conductivity ratio are found to enhance the temperature distribution.     

Hydrodynamic and thermal slip flow boundary layers over a flat plate with constant heat flux boundary condition

In this paper the boundary layer flow over a flat plat with slip flow and constant heat flux surface condition is studied. Because the plate surface temperature varies along the x direction, the momentum and energy equations are coupled due to the presence of the temperature gradient along the plate surface. This coupling, which is due to the presence of the thermal jump term in Maxwell slip condition, renders the momentum and energy equations non-similar. As a preliminary study, this paper ignores this coupling due to thermal jump condition so that the self-similar nature of the equations is preserved. Even this fundamental problem for the case of a constant heat flux boundary condition has remained unexplored in the literature. It was therefore chosen for study in this paper. For the hydrodynamic boundary layer, velocity and shear stress distributions are presented for a range of values of the parameter characterizing the slip flow. This slip parameter is a function of the local Reynolds number, the local Knudsen number, and the tangential momentum accommodation coefficient representing the fraction of the molecules reflected diffusively at the surface. As the slip parameter increases, the slip velocity increases and the wall shear stress decreases. These results confirm the conclusions reached in other recent studies. The energy equation is solved to determine the temperature distribution in the thermal boundary layer for a range of values for both the slip parameter as well as the fluid Prandtl number. The increase in Prandtl number and/or the slip parameter reduces the dimensionless surface temperature. The actual surface temperature at any location of x is a function of the local Knudsen number, the local Reynolds number, the momentum accommodation coefficient, Prandtl number, other flow properties, and the applied heat flux.     

Steady magnetohydrodynamic flow in a diverging channel with suction or blowing

An analysis is made of steady two-dimensional divergent flow of an electrically conducting incompressible viscous fluid in a channel formed by two non-parallel walls, the flow being caused by a source of fluid volume at the intersection of the walls. The fluid is permeated by a magnetic field produced by an electric current along the line of intersection of the channel walls. The walls are porous and subjected to either suction (k > 0) or blowing (k < 0) of equal magnitude on both the walls. It is found that when the Reynolds number for the flow is large and the magnetic Reynolds number is very small, boundary layers are formed on the channel walls such that a sufficient condition for the existence of a unique boundary layer solution (without separation) in the case of suction is N > 2, N being the magnetic parameter. When k = 0, boundary layer exists without separation only when N > 2. Further, it is found that the necessary and sufficient condition for the existence of a unique solution for boundary layer flow (without separation) even in the presence of blowing (k < 0) is N > 2. For given value of k, velocity at a point increases with increase in N. It is also shown that when N > 2, blowing makes the boundary layer thinner. A similarity solution for steady temperature distribution in the divergent flow is also presented when the channel walls are held at variable temperature. It is found that for fixed value of wall suction, temperature at a point decreases with increase in N. It is further shown that when N > 2, steady distribution of temperature exists even in the case of blowing at the walls.     

Bifurcation of multimode flows of a viscous fluid in a plane diverging channel

The pattern of steady multimode flow of a viscous incompressible fluid in a plane diverging channel is constructed and investigated. It is shown that odd-mode flows have velocity profiles that are symmetrical about the axis of the channel and from one to three different flows with a fixed number of modes exist. The even-mode flows are asymmetric and exist as pairs. The existence of a denumerable set of finite ranges adjoining one another, in which a single-type of complex bifurcation of the flow occurs, is established in the case of an unbounded range of values of the Reynolds number. As the Reynolds number increases, transitions to flows with an increasing number of modes, containing domains of forward and backward flows, occur successively. Flow patterns with a smaller number of modes do not occur. An increase in the number of an range corresponding to an increase in the Reynolds number leads to an unlimited increase in the length of the range and the number of modes of permissible flows.     

The initial value problem for creeping flow of the upper convected Maxwell fluid at high Weissenberg number

We consider the equations for time dependent creeping flow of an upper convected Maxwell fluid. For finite Weissenberg number, these equations can be reformulated as a coupled system of a hyperbolic equation for the stresses and an elliptic equation for the velocity. In the high Weissenberg number limit, however, the elliptic equation becomes degenerate. As a consequence, the initial value problem is no longer uniquely solvable if we just naively let the Weissenberg number go to infinity in the equations. In this paper, we make an a priori assumption on the stresses, which is motivated by the behavior in shear flow. We formulate a systematic perturbation procedure to solve the resulting initial value problem. Copyright © 2014 JohnWiley & Sons, Ltd.     

On the flow of a non-Newtonian fluid between rotating,coaxial discs

In recent times, workers in numerical methods for non-Newtonian flows have been concerned with the high Weissenberg number problem. This problem which may be caused by a number of different things, manifests itself in the failure of numerical methods at some finite and often small Weissenberg number.This paper is concerned with the torsional flow of an Oldroyd-B fluid. The kinematics is restricted to that commonly referred to as Von Kármán kinematics. These restrictions allow the reduction of the problem to a set of ordinary differential equations. The problem is then solved with finite differences using well-known branch following and jumping techniques.The solution of this set of equations is discussed, and it is found that the solutions either lose stability at subcritical bifurcation points, or fold back on themselves at limit points. Either of these will cause a high Weissenberg number problem.Comparisons are also made with the known solutions to the Newtonian problem by considering small values of the Weissenberg number.     

竖直毛细管中有限长液柱的粘性流体运动

本文考虑了竖直毛细管中具有两个自由面的有限长液柱的粘性流体运动.假设流体是牛顿的对边界条件进行线性化后,得到了小雷诺数情形下速度、压力和自由面形状的级数形式的分析表达式,对水和血液在多种液柱长度下求出了数值结果.分析表明在上下弯月面处有强回流.最大回流速度可达主流平均速度57%左右此外,本文还研究了惯性效应.采用时间相关的有限差分法求出了R_e=24.5时非线性方程的数值解.将此数值解和小雷诺数时的分析解进行比较表明,当R_e≤24.5时惯性效应不大.     

Peristatcally induced MHD slip flow in a porous medium due to a surface acoustic wavy wall

The influences of Hall current and slip condition on the MHD flow induced by sinusoidal peristaltic wavy wall in two dimensional viscous fluid through a porous medium for moderately large Reynolds number is considered on the basis of boundary layer theory in the case where the thickness of the boundary layer is larger than the amplitude of the wavy wall. Solutions are obtained in terms of a series expansion with respect to small amplitude by a regular perturbation method. Graphs of velocity components, both for the outer and inner flows for various values of the Reynolds number, slip parameter, Hall and magnetic parameters are drawn. The inner and outer solutions are matched by the matching process. An interesting application of the present results to mechanical engineering may be the possibility of the fluid transportation without an external pressure.     

Construction of circle bifurcations of a two-dimensional spatially periodic flow

The study by Yudovich [V.I. Yudovich, Example of the generation of a secondary stationary or periodic flow when there is loss of stability of the laminar flow of a viscous incompressible fluid, J. Math. Mech. 29 (1965) 587-603] on spatially periodic flows forced by a single Fourier mode proved the existence of two-dimensional spectral spaces and each space gives rise to a bifurcating steady-state solution. The investigation discussed herein provides a structure of secondary steady-state flows. It is constructed explicitly by an expansion that when the Reynolds number increases across each of its critical values, a unique steady-state solution bifurcates from the basic flow along each normal vector of the two-dimensional spectral space. Thus, at a single Reynolds number supercritical value, the bifurcating steady-state solutions arising from the basic solution form a circle.     

Effect of variable density on hydromagnetic mixed convection flow of a non-Newtonian fluid past a moving vertical plate

The effect of temperature-dependent density on MHD mixed convection flow of power-law fluid past a moving semi-infinite vertical plate for high temperature differences between the plate and the ambient fluid is studied. The fluid density is assumed to decrease exponentially with temperature. The usual Boussinesq approximations are not considered due to the large temperature differences. The surface temperature of the moving plate was assumed to vary according to a power-law form, that is, T_w(x) = T_∞ + Ax^γ. The fluid is permeated by a uniform magnetic field imposed perpendicularly to the plate on the assumption of small magnetic Reynolds number. A numerical shooting algorithm for two unknown initial conditions with fourth-order Runge–Kutta integration scheme has been used to solve the coupled non-linear boundary value problem. The effects of various parameters on the velocity and temperature profiles as well as the local skin-friction coefficient and the local Nusselt number are presented graphically and in the tabular form. The results show that application of Boussinesq approximations in a non-Newtonian fluid subjected to high temperature differences gives a significant error in the values of the skin-friction coefficient and the application of an external magnetic field reduces this error markedly in the case of shear-thickening fluid.     

Hydromagnetic flow near an accelerated plate with suction

The two-dimensional unsteady flow of a conducting viscous incompressible fluid past, an infinite flat plate with uniform suction, is considered in the presence of a uniform magnetic field. For a constant time, it is shown that for a given Hartmann numbera, as the cross Reynolds number β (corresponding to the suction velocity of the plate) increases, the velocity at any point of the fluid decreases and the skin friction at the plate increases. The results also hold good for a given β, asa increases if the magnetic lines of force are fixed relative to the fluid and are just opposite for the magnetic lines of force fixed relative to the plate.     

Effect of aligned magnetic field on the steady viscous flow past a circular cylinder

The steady viscous incompressible and slightly conducting fluid flow around a circular cylinder with an aligned magnetic field is simulated for the range of Reynolds numbers 100 ? Re ? 500 using the Hartmann number, M. The multigrid method with defect correction technique is used to achieve the second order accurate solution of complete non-linear Navier–Stokes equations. The magnetic Reynolds number is assumed to be small. It is observed that volume of the separation bubble decreases and drag coefficient increases as M is increased. We noticed that the upstream base pressure increases slightly with increase of M whereas downstream base pressure decreases with increase of M. The effect of the magnetic field on the flow is discussed with contours of streamlines, vorticity, plots of surface pressure and surface vorticity.