Numerical computation of pulsatile flow through a locally constricted channel

This paper deals with the numerical solution of a pulsatile laminar flow through a locally constricted channel. A finite difference technique has been employed to solve the governing equations. The effects of the flow parameters such as Reynolds number, flow pulsation in terms of Strouhal number, constriction height and length on the flow behaviour have been studied. It is found that the peak value of the wall shear stress has significantly changed with the variation of Reynolds numbers and constriction heights. It is also noted that the Strouhal number and constriction length have little effect on the peak value of the wall shear stress. The flow computation reveals that the peak value of the wall shear stress at maximum flow rate time in pulsatile flow situation is much larger than that due to steady flow. The constriction and the flow pulsation produce flow disturbances at the vicinity of the constriction of the channel in the downstream direction.     

Jeffrey流体流过有限长圆柱管时磁流体动力学的蠕动流

研究Jeffrey流体流过有限长管道时的蠕动流.在外加磁场作用时,流体呈导电性.分析是在长波长和低Reynolds数近似假设下完成.得到了压力梯度、体积流量、平均体积流量和局部壁面剪应力的表达式.研究了松弛时间、延迟时间和Hartman数,对压力、局部壁面剪应力以及蠕动泵机械效率的影响.还研究了回流现象,调查了沿管道壁波数非整倍数时的传播情况,研究有限长管道传播的内在特性.     

Mathematical modeling of pulsatile flow of non-Newtonian fluid in stenosed arteries

The pulsatile flow of blood through mild stenosed artery is studied. The effects of pulsatility, stenosis and non-Newtonian behavior of blood, treating the blood as Herschel–Bulkley fluid, are simultaneously considered. A perturbation method is used to analyze the flow. The expressions for the shear stress, velocity, flow rate, wall shear stress, longitudinal impedance and the plug core radius have been obtained. The variations of these flow quantities with different parameters of the fluid have been analyzed. It is found that, the plug core radius, pressure drop and wall shear stress increase with the increase of yield stress or the stenosis height. The velocity and the wall shear stress increase considerably with the increase in the amplitude of the pressure drop. It is clear that for a given value of stenosis height and for the increasing values of the stenosis shape parameter from 3 to 6, there is a sharp increase in the impedance of the flow and also the plots are skewed to the right-hand side. It is observed that the estimates of the increase in the longitudinal impedance increase with the increase of the axial distance or with the increase of the stenosis height. The present study also brings out the effects of asymmetric of the stenosis on the flow quantities.     

Effect of constriction height on flow separation in a two-dimensional channel

We studied numerically the effect of the constriction height on viscous flow separation past a two-dimensional channel with locally symmetric constrictions. A numerically stable scheme in primitive variables (velocity and pressure) for the solution of two-dimensional incompressible time-dependent Navier–Stokes equations is employed using finite-difference approximation in staggered grid. The wall shear stresses at different heights of the constriction are computed and presented graphically. It is noticed that the maximum stress and the length of the recirculating region associated with two shear layers of the constriction increase with the increase of the area reduction of the constriction. The critical Reynolds number for symmetry breaking bifurcation for the 50%, 60% and 70% area reduction are obtained numerically. The flow field separates after the symmetry breaking bifurcation and the symmetry of the flow depends on the Reynolds number and the height of the constriction.     

Non-Newtonian characteristics of peristaltic flow of blood in micro-vessels

Of concern in the paper is a generalized theoretical study of the non-Newtonian characteristics of peristaltic flow of blood through micro-vessels, e.g. arterioles. The vessel is considered to be of variable cross-section and blood to be a Herschel–Bulkley type of fluid. The progressive wave front of the peristaltic flow is supposed sinusoidal/straight section dominated (SSD) (expansion/contraction type); Reynolds number is considered to be small with reference to blood flow in the micro-circulatory system. The equations that govern the non-Newtonian peristaltic flow of blood are considered to be non-linear. The objective of the study has been to examine the effect of amplitude ratio, mean pressure gradient, yield stress and the power law index on the velocity distribution, wall shear stress, streamline pattern and trapping. It is observed that the numerical estimates for the aforesaid quantities in the case of peristaltic transport of blood in a channel are much different from those for flow in an axisymmetric vessel of circular cross-section. The study further shows that peristaltic pumping, flow velocity and wall shear stress are significantly altered due to the non-uniformity of the cross-sectional radius of blood vessels of the micro-circulatory system. Moreover, the magnitude of the amplitude ratio and the value of the fluid index are important parameters that affect the flow behaviour. Novel features of SSD wave propagation that affect the flow behaviour of blood have also been discussed.     

Nonlinear mathematical analysis for blood flow in a constricted artery under periodic body acceleration

This study analyses the pulsatile flow of blood through mild stenosed narrow arteries, treating the blood in the core region as a Casson fluid and the plasma in the peripheral layer as a Newtonian fluid. Perturbation method is employed to solve the resulting coupled implicit system of non-linear partial differential equations. The expressions for shear stress, velocity, wall shear stress, plug core radius, flow rate and longitudinal impedance to flow are obtained. The effects of pulsatility, stenosis depth, peripheral layer thickness, body acceleration and non-Newtonian behavior of blood on these flow quantities are discussed. It is noted that the plug core radius, wall shear stress and longitudinal impedance to flow increase as the yield stress and stenosis depth increase and they decrease with the increase of the body acceleration, pressure gradient, width of the peripheral layer thickness. It is observed that the plug flow velocity and flow rate increase with the increase of the pulsatile Reynolds number, body acceleration, pressure gradient and the width of the peripheral layer thickness and the reverse behavior is found when the yield stress, stenosis depth and lead angle increase. It is also recorded that the wall shear stress and longitudinal impedance to flow are considerably lower for the two-fluid Casson model than that of the single-fluid Casson model. It is found that the presence of body acceleration and peripheral layer influences the mean flow rate and mean velocity by increasing their magnitude significantly in the arteries.     

Analysis of blood flow through a modelled artery with an aneurysm

The intention of the present work is to carry out a systematic analysis of flow features in a tube, modelled as artery, having a local aneurysm in presence of haematocrit. The arterial model is treated to be axi-symmetric and rigid. The blood, flowing through the modelled artery, is treated to be Newtonian and non-homogeneous. For a thorough quantitative analysis of the flow characteristics such as wall pressure, flow velocity, wall shear stress, the unsteady incompressible Navier-Stokes equations in cylindrical polar co-ordinates under the laminar flow conditions are solved by using the finite-difference method. Finally, the numerical illustrations presented in this paper provide an effective measure to estimate the combined influence of haematocrit and aneurysm on flow characteristics. It is found that the magnitude of wall shear stress and also the length of separation increase with increasing values of the haematocrit parameter. The length of flow separation increases but the peak value of wall shear stress decreases with the increasing length of aneurysm. The peak value of wall shear stress as well as the length of separation increases with the increasing height of the aneurysm.     

Three-dimensional rotating flow induced by a shrinking sheet for suction

Series solution of magnetohydrodynamic (MHD) and rotating flow over a porous shrinking sheet is obtained by a homotopy analysis method (HAM). The viscous fluid is electrically conducting in the presence of a uniform applied magnetic field and the induced magnetic field is neglected for small magnetic Reynolds number. Similarity solutions of coupled non-linear ordinary differential equations resulting from the momentum equation are obtained. Convergence of the obtained solutions is ensured by the proper choice of auxiliary parameter. Graphs are sketched and discussed for various emerging parameters on the velocity field. The variations of the wall shear stress f″(0) and ?g′(0) are also tabulated and analyzed.     

Application of Adomian's approximation to blood flow through arteries in the presence of a magnetic field

The present investigation deals with the application of Adomian's decomposition method to blood flow through a constricted artery in the presence of an external transverse magnetic field which is applied uniformly. The blood flowing through the tube is assumed to be Newtonian in character. The expressions for the two-term approximation to the solution of stream function, axial velocity component and wall shear stress are obtained in this analysis. The numerical solutions of the wall shear stress for different values of Reynold number and Hartmann number are shown graphically. The solution of this theoretical result for a particular Hartmann number is compared with the integral method solution of Morgan and Young [17].     

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.     

Turbulent Poiseuille flow with near-critical wall transpiration

An incompressible, pressure-driven, fully developed turbulent flow between two parallel walls, with an extra constant transverse velocity component, is considered. A closure condition is formulated, which relates the shear stress with the first and the second derivatives of the longitudinal mean velocity. The closure condition is derived without invoking any special hypotheses on the nature of turbulent motion, only taking advantage of the fact that the flow depends on a finite number of governing parameters. By virtue of the closure condition, the momentum equation is reduced to the boundary-value problem for a second-order differential equation, which is solved by the method of matched asymptotic expansions at high values of the logarithm of the Reynolds number based on the friction velocity. The case of near-critical transpiration, when the shear stress at the injection wall vanishes, is considered. It is shown that the maximum point on the mean velocity profile lies in a thin sublayer near the suction wall in this case. A formula for the position of the maximum point as a function of the transpiration factor is obtained. The mean velocity profiles near the suction wall are calculated. A friction law for Poiseuille flow with near-critical transpiration is found, which makes it possible to describe the relation between the shear stress at the wall, the Reynolds number, and the transpiration velocity by a single function of one variable. Direct numerical simulation of the flow for some transpiration factors is performed. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

T型分叉血管的定常/脉动流动和大分子传质

采用计算流体动力学方法,数值求解了T型分叉流动的定常/脉动流场和低密度脂蛋白(LDL)以及血清白蛋白(Albumin)的浓度分布。计算了雷诺数、主管和支管的流量比等参数对流场和大分子传质的影响,计算结果表明,流体动力学因素影响大分子的分布和跨壁渗透,在动脉硬化的发生和发展过程中起着重要的作用。在流动发生分离处,即支管入口外侧壁面剪应力变化最剧烈,这儿LDL和Albumin的壁面浓度变化也是最剧烈,是动脉硬化危险区。     

Asymptotic study of turbulent Poiseuille flow with wall transpiration

An incompressible, pressure–driven, fully developed turbulent flow between two parallel walls, with an extra constant transverse velocity component, is considered. A closure condition is formulated, which relates the shear stress to the first and second derivatives of the longitudinal mean velocity. The closure condition is derived without invoking any special hypotheses on the nature of turbulent motion, only taking advantage of the fact that the flow depends on a finite number of governing parameters. By virtue of the closure condition, the momentum equation is reduced to the boundary–value problem for a second–order differential equation, which is solved by the method of matched asymptotic expansions at high values of the logarithm of the Reynolds number based on the friction velocity. A limiting transpiration velocity is obtained, such that the shear stress at the injection wall vanishes, while the maximum point on the velocity profile approaches the suction wall. In this case, a sublayer near the suction wall appears where the mean velocity is proportional to the square root of the distance from the wall. A friction law for Poiseuille flow with transpiration is found, which makes it possible to describe the relation between the wall shear stress, the Reynolds number, and the transpiration velocity by a function of one variable. A velocity defect law, which generalizes the classical law for the core region in a channel with impermeable walls to the case of transpiration, is also established. In similarity variables, the mean velocity profiles across the whole channel width outside viscous sublayers can be described by a one–parameter family of curves. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of Mass Loading on Particle-Laden Turbulent Pipe Flow

A CFD code in the framework of OpenFOAM was validated for simulations of particle-laden pipe and channel flows at low to intermediate mass loadings. The code is based on an Eulerian two-fluid approach with Reynolds-averaged conservation equations, including turbulence modeling and four-way coupling. Pipe flow simulations of particles in air against gravity were conducted at Reynolds numbers up to 50000. The particle mass loading was varied and its effect on the mean velocities and turbulent fluctuations of the two phases was studied. Special attention was paid to the influence of mass loading on the centerline velocity and the wall shear velocity of the fluid phase for various flow parameters and particle properties. Empirical correlations were established between these two quantities and the flow Reynolds number, particle Reynolds number, Stokes number and particle to fluid density ratio for a range of particle mass loadings. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

变截面圆形管道粘性流动的伽辽金-摄动杂交解

采用伽辽金-摄动杂交法来研究壁面是正弦形状的变截面圆形管道的粘性流动,从而避免了摄动小参数的局限性和单纯伽辽金法基函数选取的任意性的困难．讨论了边界和雷诺数对流动的影响,获得流动分离点和附着点的位置,还分析了壁面剪应力和摩擦系数沿轴向的变化情况．在小参数的情况下,计算所获得的结果与摄动解吻合良好．     

The influence of heat transfer on peristaltic transport of a Jeffrey fluid in a vertical porous stratum

The peristaltic flow of a Jeffrey fluid in a vertical porous stratum with heat transfer is studied under long wavelength and low Reynolds number assumptions. The nonlinear governing equations are solved using perturbation technique. The expressions for velocity, temperature and the pressure rise per one wave length are determined. The effects of different parameters on the velocity, the temperature and the pumping characteristics are discussed. It is observed that the effects of the Jeffrey number λ₁, the Grashof number Gr, the perturbation parameter N = EcPr, and the peristaltic wall deformation parameter ϕ are the strongest on the trapping bolus phenomenon. The results obtained for the flow and heat transfer characteristics reveal many interesting behaviors that warrant further study on the non-Newtonian fluid phenomena, especially the shear-thinning phenomena. Shear-thinning reduces the wall shear stress.     

Two-fluid Casson model for pulsatile blood flow through stenosed arteries: A theoretical model

Pulsatile flow of blood through mild stenosed narrow arteries is analyzed by treating the blood in the core region as a Casson fluid and the plasma in the peripheral layer as a Newtonian fluid. Perturbation method is used to solve the coupled implicit system of non-linear differential equations. The expressions for velocity, wall shear stress, plug core radius, flow rate and resistance to flow are obtained. The effects of pulsatility, stenosis, peripheral layer and non-Newtonian behavior of blood on these flow quantities are discussed. It is found that the pressure drop, plug core radius, wall shear stress and resistance to flow increase with the increase of the yield stress or stenosis size while all other parameters held constant. The percentage of increase in the resistance to flow over the uniform diameter tube is considerably very low for the present two-fluid model compared with those of the single-fluid model.     

Pulsatile flow of particle-fluid suspension model of blood under periodic body acceleration

A particle-fluid suspension model is applied to the problem of pulsatile blood flow through a circular tube under the influence of body acceleration. With the help of finite Hankel and Laplace transforms, analytic expressions for axial velocity for both fluid and particle phase, fluid acceleration, wall shear stress and instantaneous flow rate have been obtained. It is observed that the solutions can be used for all feasible values of pulsatile and body acceleration Reynolds numbers R_p and R_b. Using physiological data, the following qualitative and quantitative results have been obtained. The amplitude Q_b of instantaneous flow rate due to body acceleration decreases as the tube radius decreases. The effect of the volume fraction of particle C on Q_b is to increase it with increase of C in arteriole and to decrease Q_b as C increases in coronary and femoral arteries. The maximum of the axial velocity and fluid acceleration shifts from the axis of the tube to the vicinity of the tube wall as the tube diameter increases. The effect of C on the velocity and acceleration are nonuniform. The wall shear amplitude t_b\tau_b due to body acceleration increases as the tube diameter decreases from femoral to coronary and a further decrease in the tube diameter leads to a decrease in t_b\tau_b. The effects of C on t_b\tau_b are again nonuniform.     

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.     

微极液体在两个多孔圆盘间的MHD流动及其热传导

两个平行的无限大多孔圆盘,圆盘表面有均匀注入时,数值地研究圆盘间不可压缩导电微极流体,在横向外加磁场作用下的轴对称稳定层流.运用von Krmn的相似变换,将非线性运动的控制方程转化为无量纲形式.使用基于有限差分格式的算法,在相应的边界条件下,求解简化后耦合的常微分方程组.讨论Reynolds数、磁场参数、微极参数和Prandtl数,对流动速度和温度分布的影响.在特殊情况下,所得结果与已有文献的工作有着很好的一致性.研究表明,圆盘表面的传热率随着Rynolds数、磁场参数和Prandtl数的增加而增加;剪切应力随着注入的增加而减少,但它随着外部磁场的加强而增加.和Newton流体相比较,微极流体的剪切应力因素较弱,有利于聚合体加工过程中流动和温度的控制.