Analytical and numerical solutions of a generalized hyperbolic non-Newtonian fluid flow

Generalized hyperbolic non-Newtonian fluid model first proposed by Al-Zahrani [1] is considered. The model was successfully applied to some drilling fluids with better performance in relating shear stress and velocity gradient compared to power-law and Hershel-Bulkley model. Special flow geometries namely pipe flow, parallel plate flow and flow between two rotating cylinders are treated. For the first two cases, analytical solutions of velocity profiles in the form of integrals are presented. For the flow between two rotating cylinders, the differential equation is solved by Runge-Kutta method combined with shooting. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Flow field of a third-grade non-Newtonian fluid in the annulus of rotating concentric cylinders in the presence of magnetic field

In this paper, the flow field of a third-grade non-Newtonian fluid in the annulus of rotating concentric cylinders has been investigated in the presence of magnetic field. For this purpose, the constitutive equation of such a fluid flow was simplified, and the existence of the solution to the governing equation was established using Schauder's fixed point theorem. Using the finite difference method, the numerical solution of the non-dimensionalized form of the established governing equation was obtained. The effect of sundry parameters such as the rotating speed of the cylinders, the physical properties of fluid, and magnetic field intensity on the fluid velocity field was studied as well.     

Rheological effect on thermocapillary flow of a liquid film jet painted on a moving boundary

In the present paper, a liquid (or melt) film of relatively high temperature ejected from a vessel and painted on the moving solid film is analyzed by using the second-order fluid model of the non-Newtonian fluid. The thermocapillary flow driven by the temperature gradient on the free surface of a Newtonian liquid film was discussed before. The effect of rheological fluid on thermocapillary flow is considered in the present paper. The analysis is based on the approximations of lubrication theory and perturbation theory. The equation of liquid height and the process of thermal hydrodynamics of the non-Newtonian liquid film are obtained, and the case of weak effect of the rheological fluid is solved in detail.     

二维分岔槽道内非牛顿流体流动的有限元分析

本文介绍二维分岔槽道内非牛顿流体流动的有限元分析.采用Galerkin法及混合有限元法,流体看作不可压缩的非牛顿流体,满足Oldyord微分型本构方程.由有限元法形成的非线性代数方程组用连续微分法求解.结果表明有限元法适于分析复杂流场中非牛顿流体的流动.     

非牛顿流体内流传热研究

在本文中,研究了注入轴对称模腔非牛顿流体非定常流动.本文的第二部份研究了上随体Maxwell流体管内热流动.对于注入模腔流动.其本构方程采用幂律流体模型方程.为了避免在表现粘度中温度关系引起的非线性.引进了一特征粘度的概念.描述本力学过程的基本方程是,本构方程、定常状态的运动方程、非定常能量方程及连续方程.该方程组在空间是二维问题,在数学上是三维问题.采用分裂差分格式求得本方程组的数值解答.分裂法曾成功应用于求解牛顿流体问题.在本文中,首次将分裂法成功地应用解决非牛顿流体流动问题.对于圆管内热流,给出了差分格式,使基本方程组化为一个三对角方程组.其结果,给出了不同时刻的模腔内二维温度分布.     

Direct numerical simulation of turbulent non-Newtonian flow using a spectral element method

A spectral element—Fourier method (SEM) for Direct Numerical Simulation (DNS) of the turbulent flow of non-Newtonian fluids is described and the particular requirements for non-Newtonian rheology are discussed. The method is implemented in parallel using the MPI message passing kernel, and execution times scale somewhat less than linearly with the number of CPUs, however this is more than compensated by the improved simulation turn around times. The method is applied to the case of turbulent pipe flow, where simulation results for a shear-thinning (power law) fluid are compared to those of a yield stress (Herschel–Bulkley) fluid at the same generalised Reynolds number. It is seen that the yield stress significantly dampens turbulence intensities in the core of the flow where the quasi-laminar flow region there co-exists with a transitional wall zone. An additional simulation of the flow of blood in a channel is undertaken using a Carreau–Yasuda rheology model, and results compared to those of the one-equation Spalart-Allmaras RANS (Reynolds-Averaged Navier–Stokes) model. Agreement between the mean flow velocity profile predictions is seen to be good. Use of a DNS technique to study turbulence in non-Newtonian fluids shows great promise in understanding transition and turbulence in shear thinning, non-Newtonian flows.     

Flow of polymer melts in a combined extruder

A viscoelastic incompressible fluid, modeled by the White-Metzner equation, flows in a combined extruder. Approximate expressions are obtained for the secondary flows in the disk zone of the machine whose presence has been confirmed experimentally. By means of these expressions it is possible to determine the total shear deformation of an element of the melt as it passes through the extruder.Kiev Polytechnic Institute. L'vov Polytechnic Institute. Translated from Mekhanika Polimerov, No. 2, pp. 351–358, March–April, 1972.     

Effect of Deborah number and phase difference on peristaltic transport of a third-order fluid in an asymmetric channel

The effect of a third-order fluid on the peristaltic transport in an asymmetric channel is studied. The wavelength of the peristaltic waves is assumed to be large compared to the varying channel width, whereas the wave amplitudes need not be small compared to the varying channel width. The channel asymmetry is produced by choosing the peristaltic wave train on the walls to have different amplitudes and phase. The flow is investigated in a wave frame of reference moving with velocity of the wave. The effects of Deborah number, phase difference, varying channel width and wave amplitudes on the pumping characteristics, streamline pattern and trapping phenomena are investigated. It is observed that the trapping regions increase as the channel becomes more and more symmetric and the trapped bolus volume decreases for increasing Deborah number, phase difference and varying channel width whereas it increases for increasing flow rate and wave amplitudes. Furthermore, the obtained results could also have applications to a range of peristaltic flows for a variety of non-Newtonian fluids such as aqueous solutions of high-molecular weight polyethylene oxide and polyacrylamide.     

Pulsatile flow of Herschel–Bulkley fluid through catheterized arteries – A mathematical model

The pulsatile flow of blood through catheterized artery has been studied in this paper by modeling blood as Herschel–Bulkley fluid and the catheter and artery as rigid coaxial circular cylinders. The Herschel–Bulkley fluid has two parameters, the yield stress θ and the power index n. Perturbation method is used to solve the resulting quasi-steady nonlinear coupled implicit system of differential equations. The effects of catheterization and non-Newtonian nature of blood on yield plane locations, velocity, flow rate, wall shear stress and longitudinal impedance of the artery are discussed. The existence of two yield plane locations is investigated and their dependence on yield stress θ, amplitude A, and time t are analyzed. The width of the plug core region increases with increasing value of yield stress at any time. The velocity and flow rate decrease, whereas wall shear stress and longitudinal impedance increase for increasing value of yield stress with other parameters held fixed. On the other hand, the velocity, flow rate and wall shear stress decrease but resistance to flow increases as the catheter radius ratio (ratio of catheter radius to vessel radius) increases with other parameters fixed. The results for power law fluid, Newtonian fluid and Bingham fluid are obtained as special cases from this model.     

Laminar non-Newtonian fluid flow in a porous annulus

The effect of suction and injection through the porous walls of an annulus on a two-dimensional steady laminar flow of a non-Newtonian fluid is investigated by solving directly the Navier-Stokes equations in cylindrical coordinates. A quasi-linearization method is used to solve the flow problem for both small and moderately large cross flow through the porous walls. With very approximate starting values of functions, only a few iterations are needed to obtain results with very high accuracy. For comparison, a perturbation method is also used to solve the same flow problem. For small values of cross flow through the porous wall, results obtained by both methods agree very well; while for the moderately large values of cross flow, the perturbation technique is impractical and may give an inaccurate result. The effects of the non-Newtonian fluid on the axial velocity, the pressure drop, and skin friction are also examined.     

Motion of a non-Newtonian liquid in the annular channel of an insulation-coating tool

The axisymmetric flow of a non-Newtonian liquid with an exponential rheological law between two coaxial cylinders is considered for the case when the outer cylinder is fixed and the inner cylinder moves freely together with the liquid. At the outer cylinder the liquid moves either with or without slip. The individual theoretical conclusions are confirmed by the experimental results.Mekhanika Polimerov, Vol. 3, No. 1, pp. 151–155, 1967     

The mixing of a viscous fluid in a layer between rotating eccentric cylinders

The motion of an incompressible viscous fluid in a thin layer between two circular cylinders, inserted into one another, with parallel axes is investigated. The cylinders rotate relative to one another about an axis parallel to the axes of the cylinders. The stream function of the unsteady plane-parallel flow that occurs is found by solving the boundary-value problem for the equations of hydrodynamic lubrication theory. The motion of the fluid particles is found from the solution of a non-autonomous time-periodic Hamiltonian system with a Hamiltonian equal to the stream function. The positions of fluid particles over time intervals that are a multiple of the period of rotation (Poincaré points) are calculated. The set of points is investigated using a Poincaré mapping on the phase flow. The observed transition to chaotic motion is related to the mixing of the fluid particles and is investigated both numerically and using a mapping, calculated with an accuracy up to the third power of the small eccentricity. The optimum mode of motion is observed when the area of the mixing (chaos) region reaches its highest value.     

Steady flow in the Willis circle using a Quemada model

The paper presents new results on the dynamics of blood flow in the Willis circle, by taking into account the non-Newtonian nature of the blood through a Quemada model. In previous papers we performed a steady analysis of the blood flow in the Willis polygon, by considering the blood as a Newtonian fluid. In hemodynamics this is the first step, while taking into account the non-Newtonian behaviour of the fluid is called a second order approach. This is in fact the aim of our paper, to compare the results (blood flow rates and pressures) obtained in a non-Newtonian approach with those furnished by a Newtonian modeling of the blood. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Slow rotation of a sphere in a non-Newtonian fluid with or without suction or injection

The flow generated by the rotation of a sphere in an infinitely extending fluid has recently been studied by Goldshtik. The corresponding problem for non-Newtonian Reiner-Rivlin fluids has been studied by Datta. Bhatnagar and Rajeswari have studied the secondary flow between two concentric spheres rotating about an axis in the non-Newtonian fluids. This last investigation was further generalised by Rajeswari to include the effects of small radial suction or injection. In Part A of the present investigation, we have studied the secondary flow generated by the slow rotation of a single sphere in non-Newtonian fluid obeying the Rivlin-Ericksen constitutive equation. In Part B, the effects of small suction or injection have been studied which is applied in an arbitrary direction at the surface of the sphere. In the absence of suction or injection, the secondary flow for small values of the visco-elastic parameter is similar to that of Newtonian fluids with inclusion of inertia terms in the Oseen approximation. If this parameter exceeds K_c = 18R/219, whereR is the Reynolds number, the breaking of the flow field takes place into two domains, in one of which the stream lines form closed loops. For still higher values of this parameter, the complete reversal of the sense of the flow takes place. When suction or injection is included, the breaking of the flow persists under certain condition investigated in this paper. When this condition is broken, the breaking of the flow is obliterated.     

Non-linear peristaltic transport of a second-order fluid through a porous medium

The present paper investigates phenomena brought about into the classic peristaltic mechanism by inclusion of non-Newtonian effects through a porous space in a channel. The peristaltic motion of a second-order fluid through a porous medium was studied for the case of a planar channel with harmonically undulating extensible walls. The system of the governing nonlinear PDE is solved by using the perturbation method to second-order in dimensionless wavenumber. The analytic solution has been obtained in the form of a stream function from which the axial pressure gradient has been derived. The flow is investigated in a wave frame of reference moving with velocity of the wave. Numerical calculations are carried out for the pressure rise and frictional force. The features of the flow characteristics are analyzed by plotting graphs and discussed in detail.     

生物磁粘弹性流体的流动:应用动脉电磁过热评估血液的流动,癌症治疗进程

利用生物磁流体动力学(BFD)原理,在生物磁流体经由遭受磁场作用的多孔介质时,研究其流动的基本理论.所研究流体的磁化强度随温度而变化.流体被认为是非Newton流体,其流动由二阶梯度流体方程所控制,并考虑了流体的粘弹性效应.假设管道壁是能够伸展的,管壁表面的速度与到坐标原点的纵向距离成正比.首先将问题简化为包括7个参数的、耦合的非线性微分方程组的求解.将血液看作生物磁流体,并用上述方法分析,目的是计算某些血液的流动参数,并配以适当的数值方法,导数用差分格式近似.计算结果用图形给出,从而在磁场作用下,得到过热状态中关系血液的、血流动力学流动的理论预测.结果清楚地表明,在电磁过热治疗进程期间,磁偶极子对动脉中血液流动特征的影响起着重大作用.该研究引起了临床医学的关注,其结果有益于癌症病人采用电磁过热的治疗.     

A non-Newtonian fluid flow model for blood flow through a catheterized artery—Steady flow

In this paper the effects catheterization and non-Newtonian nature of blood in small arteries of diameter less than 100 μm, on velocity, flow resistance and wall shear stress are analyzed mathematically by modeling blood as a Herschel–Bulkley fluid with parameters n and θ and the artery and catheter by coaxial rigid circular cylinders. The influence of the catheter radius and the yield stress of the fluid on the yield plane locations, velocity distributions, flow rate, wall shear stress and frictional resistance are investigated assuming the flow to be steady. It is shown that the velocity decreases as the yield stress increases for given values of other parameters. The frictional resistance as well as the wall shear stress increases with increasing yield stress, whereas the frictional resistance increases and the wall shear stress decreases with increasing catheter radius ratio k (catheter radius to vessel radius). For the range of catheter radius ratio 0.3–0.6, in smaller arteries where blood is modeled by Herschel–Bulkley fluid with yield stress θ = 0.1, the resistance increases by a factor 3.98–21.12 for n = 0.95 and by a factor 4.35–25.09 for n = 1.05. When θ = 0.3, these factors are 7.47–124.6 when n = 0.95 and 8.97–247.76 when n = 1.05.     

Modelling of Newtonian and Non-Newtonian Incompressible Fluid Flow

This paper deals with numerical simulation of steady laminar incompressible Newtonian and non-Newtonian fluid flow. The viscosity of non-Newtonian fluid is modelled using the fluid power law. The mathematical model is described by the non-linear conservative system of the incompressible Navier-Stokes equations that is solved using the explicit two-step MacCormack scheme with artificial viscosity together with the pseudo-compressibility method. The spatial discretization is performed by the cell-centred finite volume method on a structured quadrilateral grid. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exact solutions for the rotational flow of a generalized Maxwell fluid between two circular cylinders

Unsteady flow of an incompressible generalized Maxwell fluid between two coaxial circular cylinders is studied by means of the Laplace and finite Hankel transforms. The motion of the fluid is produced by the rotation of cylinders around their common axis. The solutions that have been obtained, written in integral and series form in terms of the generalized G_a,b,c(·, t)-functions, are presented as a sum of the Newtonian solutions and the corresponding non-Newtonian contributions. They satisfy all imposed initial and boundary conditions and for λ → 0 reduce to the solutions corresponding to the Newtonian fluids performing the same solution. Furthermore, the corresponding solutions for ordinary Maxwell fluids are also obtained for β = 1. Finally, in order to reveal some relevant physical aspects of the obtained results, the diagrams of the velocity field ω(r, t) have been depicted against r and t for different values of the material and fractional parameters.     

Investigation of the Free-Surface Flow in Partially Filled Extruders

In this work we present the free-surface flow of high viscous Newtonian liquids in a simplified extruder model, namely the flat-plate model. Due to high density and viscosity ratios, a free-surface approach based on the volume of fluid method is used. The material distribution for different degree of fillings and the resulting power consumption is investigated. In addition, the obtained power characteristic is compared against analytical consideration and numerical simulation of a fully filled case. The investigation shows that the power characteristic of a partially filled extruder can be estimated by a fully filled computation. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)