Two-fluid non-linear model for flow in catheterized blood vessels

Blood flow through a catheterized artery is analyzed, assuming the flow is steady and blood is treated as a two-fluid model with the suspension of all the erythrocytes in the core region as a Casson fluid and the plasma in the peripheral region as a Newtonian fluid. The expressions for velocity, flow rate, wall shear stress and frictional resistance are obtained. The variations of these flow quantities with yield stress, catheter radius ratio and peripheral layer thickness are discussed. It is noticed that the velocity and flow rate decrease while the wall shear stress and resistance to flow increase when the yield stress or the catheter radius ratio increases while all the other parameters were held fixed. It is found that the velocity and flow rate increase while the wall shear stress and frictional resistance decrease with the increase of the peripheral layer thickness. The estimates of the increase in the frictional resistance are significantly very small for the present two-fluid model than those of the single-fluid Casson model.     

Two-phase non-linear model for blood flow in asymmetric and axisymmetric stenosed arteries

The pulsatile flow of a two-phase model for blood flow through axisymmetric and asymmetric stenosed narrow arteries is analyzed, treating blood as a two-phase model with the suspension of all the erythrocytes in the core region as the Herschel-Bulkley material and plasma in the peripheral layer as the Newtonian fluid. The perturbation method is applied to solve the resulting non-linear implicit system of partial differential equations. The expressions for various flow quantities are obtained. It is found that the pressure drop, plug core radius, wall shear stress increase as the yield stress or stenosis height increases. It is noted that the velocity increases, longitudinal impedance decreases as the amplitude increases. For asymmetric stenosis, the wall shear stress increases non-linearly with the increase of the axial distance. The estimates of the increase in longitudinal impedance to flow of the two-phase Herschel-Bulkley material are significantly lower than those of the single-phase Herschel-Bulkley material. The results show the advantages of two-phase flow over single-phase flow in small diameter arteries with stenosis.     

Effect of couple stresses on the flow in a constricted annulus

The flow of an incompressible couple stress fluid in an annulus with local constriction at the outer wall is considered. This configuration is intended as a simple model for studying blood flow in a stenosed artery when a catheter is inserted into it. The effects couple stress fluid parameters α and σ, height of the constriction (ε), and ratio of radii (k) on the impedance and wall shear stresses are studied graphically. Graphical results show that the resistance to the flow as well as the wall shear stress increases as the ratio of the radii increases and decreases as the couple stress fluid parameters increases.     

Blood flow in stenosed arteries with radially variable viscosity, peripheral plasma layer thickness and magnetic field

A novel approach of combined mathematical and computational models has been developed to investigate the oscillatory two-layered flow of blood through arterial stenosis in the presence of a transverse uniform magnetic field applied. Blood in the core region and plasma fluid in the peripheral layer region are assumed to obey the law of Newtonian fluid. An analytical solution is obtained for velocity profile and volumetric flow rate in the peripheral plasma region and also wall shear stress. Finite difference method is employed to solve the momentum equation for the core region. The numerical solutions for velocity, flow rate and flow resistance are computed. The effects of various parameters associated with the present flow problem such as radially variable viscosity, hematocrit, plasma layer thickness, magnetic field and pulsatile Reynolds number on the physiologically important flow characteristics namely velocity distribution, flow rate, wall shear stress and resistance to flow have been investigated. It is observed that the velocity increases with the increase of plasma layer thickness. An increase or a decrease in the velocity and wall shear stress against the increase in the value of magnetic parameter (Hartmann number) and hematocrit is dependent on the value of t. An increase in magnetic field leads to an increase in the flow resistance and it decreases with the increase in the plasma layer thickness and pulsatile Reynolds number. The information concerning the phase lag between the flow characteristics and how it is affected by the hematocrit, plasma layer thickness and Hartmann number has, for the first time, been added to the literature.     

Mathematical modelling of couple stresses on fluid flow in constricted tapered artery in presence of slip velocity-effects of catheter

This paper explores the mathematical model for couple stress fluid flow through an annular region. The above model is used for studying the blood flow be-tween the clogged （stenotic） artery and the catheter. The asymmetric nature of the stenosis is considered. The closed form expressions for the physiological parameters such as impedance and shear stress at the wall are obtained. The effects of various geomet-ric parameters and the parameters arising out of the fluid considered are discussed by considering the slip velocity and tapering angle. The study of the above model is very significant as it has direct applications in the treatment of cardiovascular diseases.     

A Two-Layer Mathematical Model of Blood Flow in Porous Constricted Blood Vessels

In this paper, we discussed a mathematical model for two-layered non-Newtonian blood flow through porous constricted blood vessels. The core region of blood flow contains the suspension of erythrocytes as non-Newtonian Casson fluid and the peripheral region contains the plasma flow as Newtonian fluid. The wall of porous constricted blood vessel configured as thin transition Brinkman layer over layered by Darcy region. The boundary of fluid layer is defined as stress jump condition of Ocha-Tapiya and Beavers–Joseph. In this paper, we obtained an analytic expression for velocity, flow rate, wall shear stress. The effect of permeability, plasma layer thickness, yield stress and shape of the constriction on velocity in core & peripheral region, wall shear stress and flow rate is discussed graphically. This is found throughout the discussion that permeability and plasma layer thickness have accountable effect on various flow parameters which gives an important observation for diseased blood vessels.     

Effects of the side walls on unsteady flow of a second grade fluid over a plane wall

The effects of the side walls on unsteady flow of a second grade fluid over a plan wall are considered. The solution of the governing equation for velocity is obtained by the sine transform method. This gives a correct result for the shear stress at the bottom wall. The shear stress at the bottom wall is minimum at the middle of the plate and it increases near the side walls. It is shown that the mean thickness of the layer of the liquid over the plate increases with time and the ratio of the mean thickness to the distance between the side walls becomes ultimately 0.2714.     

Pulsatile flow of Herschel-Bulkley fluid through stenosed arteries—A mathematical model

In this paper, the pulsatile flow of blood through stenosed artery is studied. The effects of pulsatility, stenosis and non-Newtonian behavior of blood, assuming the blood to be represented by Herschel-Bulkley fluid, are simultaneously considered. A perturbation method is used to analyze the flow assuming the thickness of plug core region to be non-uniform changing with axial distance. An expression for the variation of plug core radius with time and axial distance is obtained. The variation of pressure gradient with steady flow rate is given. Also the variation of wall shear stress distribution as well as resistance to flow with axial distance for different values of time and for different values of yield stress is given and the results analyzed.     

A fluid flow model in the lacunar-canalicular system under the pressure gradient and electrical field driven loads

The lacunar-canalicular system (LCS) is acknowledged to directly participate in bone tissue remodeling. The fluid flow in the LCS is synergic driven by the pressure gradient and electric field loads due to the electro-mechanical properties of bone. In this paper, an idealized annulus Maxwell fluid flow model in bone canaliculus is established, and the analytical solutions of the fluid velocity, the fluid shear stress, and the fluid flow rate are obtained. The results of the fluid flow under pressure gradient driven (PGD), electric field driven (EFD), and pressure-electricity synergic driven (P-ESD) patterns are compared and discussed. The effects of the diameter of canaliculi and osteocyte processes are evaluated. The results show that the P-ESD pattern can combine the regulatory advantages of single PGD and EFD patterns, and the osteocyte process surface can feel a relatively uniform shear stress distribution. As the bone canalicular inner radius increases, the produced shear stress under the PGD or P-ESD pattern increases slightly but changes little under the EFD pattern. The increase in the viscosity makes the flow slow down but does not affect the fluid shear stress (FSS) on the canalicular inner wall and osteocyte process surface. The increase in the high-valent ions does not affect the flow velocity and the flow rate, but the FSS on the canalicular inner wall and osteocyte process surface increases linearly. In this study, the results show that the shear stress sensed by the osteocyte process under the P-ESD pattern can be regulated by changing the pressure gradient and the intensity of electric field, as well as the parameters of the annulus fluid and the canaliculus size, which is helpful for the osteocyte mechanical responses. The established model provides a basis for the study of the mechanisms of electro-mechanical signals stimulating bone tissue (cells) growth.

     

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.     

Wave propagation in an elastic tube: a numerical study

In this paper, a fluid–wall interaction model, called the elastic tube model, is introduced to investigate wave propagation in an elastic tube and the effects of different parameters. The unsteady flow was assumed to be laminar, Newtonian and incompressible, and the vessel wall to be linear-elastic, isotropic and incompressible. A fluid–wall interaction scheme is constructed using a finite element method. The results demonstrate that the elastic tube plays an important role in wave propagation. It is shown that there is a time delay between the velocity waveforms at two different locations and that the peak velocity increases while the low velocity decreases in the elastic tube model, contrary to the rigid tube model where velocity waveforms overlap each other. Compared with the elastic tube model, the increase of the wall thickness makes wave propagation faster and the time delay cannot be observed clearly, however, the velocity amplitude is reduced slightly due to the decrease of the internal radius. The fluid–wall interaction model simulates wave propagation successfully and can be extended to study other mechanical properties considering complicated geometrical and material factors.     

一种确定均匀动脉壁面切应力的非线性方法

从Ling和Atabek提出的``局部流'理论出发，提出一种利用测量血液黏度、管轴上 的血流速度、压力和管径波形计算均匀动脉管壁切应力的非线性方法. 将这种方法与柳兆荣 等提出的利用测量血液黏度、管轴上的血流速度和平均管径计算切应力的线性方法比较，结 果表明，当管壁脉动幅度较小时，两种方法计算的压力梯度、流速剖面和管壁切应力差别较 小；而当管壁脉动幅度增大时，两种方法计算的压力梯度、流速剖面和管壁切应力差别增大. 对于小幅脉动均匀动脉，用线性方法计算管壁切应力有较高的精度；而对于大变形 均匀动脉，则需要考虑非线性因素对管壁切应力的影响. 由于作为输入量的血液黏度、轴心 血流速度、压力波形和管径波形可在活体上通过无损伤或微损伤的检测方法得到， 所提出的计算切应力的方法为在体或离体研究切应力与动脉重建的关系提供了方法学基础.     

动脉瘤内流场以及瘤体尺寸的影响的数值研究

采用计算流体动力学(CFD)数值模拟的方法,在周期性脉动速度入流条件下,建立刚性动脉瘤模型并研究了动脉瘤模型中流场的特征(速度、压力、壁面剪切应力)。得到了脉动入流一个周期内流场特征的变化规律,发现动脉瘤的后端有相当高的压力和壁面剪切应力,而且高压力和壁面剪切应力分布的位置几乎是固定的。探讨了不同动脉瘤尺寸对内部流场的影响,动脉瘤的直径与瘤体长度之比越大,瘤壁承受的剪切应力就越大,动脉瘤破裂的危险性就越高。     

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.     

Flow pattern transition,pressure gradient,hold-up predictions in gas/non-Newtonian power-law fluid stratified flow

This work focuses on gas/non-Newtonian power-law fluid stratified pipe flow. Two different theoretical approaches to obtain pressure gradient and hold-up predictions are presented: the steady fully developed two-fluid model and the pre-integrated model. The theoretical predictions are compared with experimental data available for horizontal and for slightly downward inclined air/shear thinning fluid stratified flow taken from literature. The predictions of the pre-integrated model are validated showing a good agreement when compared with experimental data. The criteria for the transition from the stratified flow pattern are applied to gas/non-Newtonian stratified flow. The neutral stability analysis (smooth/wavy stratified flow) and the well-posedness (existence region of stratified flow) of governing equations are carry out. The predicted transition boundaries are obtained using the steady fully developed two-fluid model and the pre-integrated model, where the shape factors and their derivatives are accounted for. A comparison between the predicted boundaries and experimental flow pattern maps is presented and shows a good agreement. A comment on the shear stress modeling by the pre-integrated model is provided.     

一种骨小管中液体流动产生的流量及切应力模型

骨组织受力变形后其内部液体就会流动,同时在其微观结构——骨单元壁中扩散,并进一步产生一系列与骨液流动相关的物理效应,如流体剪切应力、流动电位等,这些物理效应被细胞感知并做出破骨或成骨等反应,来使骨适应外部载荷环境.鉴于骨组织产生的内部液体流动很难实验测定,理论模拟是目前的主要研究手段.基于骨单元的多孔弹性性质建立了骨小管内部液体的流动模型,该模型将骨单元所受的外部载荷与骨小管内部液体的压力、流速、流量和切应力联系起来,并进一步可以研究其力传导与力电传导机制.骨小管模型的建立分别基于中空和考虑哈弗液体的骨单元模型,并考虑了骨单元外壁的弹性约束和刚性位移约束两种边界条件.最终得到骨单元在外部轴向载荷作用下,骨小管内部液体的流量及流体切应力的解析解.结果表明:骨小管中的液体流量与流体切应力都正比于应变载荷幅值和频率,并由载荷的应变率决定.因此应变率可以作为控制流量和流体切应力的一种生理载荷因素.流量随着骨小管半径的增大而非线性增大,而流体切应力则随着骨小管半径的增大而线性增大.此外,在相同的载荷下,含哈弗液体的骨单元的模型中,骨小管中液体的流量和切应力均大于中空骨单元模型.     

磁场和Hall电流对狭窄动脉中血液流动的影响

A micropolar model for blood simulating magnetohydrodynamic flow through a horizontally nonsymmetric but vertically symmetric artery with a mild stenosis is presented. To estimate the effect of the stenosis shape, a suitable geometry has been considered such that the horizontal shape of the stenosis can easily be changed just by varying a parameter referred to as the shape parameter. Flow parameters, such as 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 shape parameters, the Hartmann number and the Hall parameter. This shows that the resistance to flow decreases with the increasing values of the parameter determining the stenosis shape and the Hail parameter, while it increases with the increasing Hartmann number. The wall shear stress and the shearing stress on the wall at the maximum height of the stenosis possess an inverse characteristic to the resistance to flow with respect to any given value of the Hartmann number and the Hall parameter. Finally, the effect of the Hartmann number and the Hall parameter on the horizontal velocity is examined.     

Rheology of shear thickening suspensions and the effects of wall slip in torsional flow

The rheological characterisation of concentrated shear thickening materials suspensions is challenging, as complicated and occasionally discontinuous rheograms are produced. Wall slip is often apparent and when combined with a shear thickening fluid the usual means of calculating rim shear stress in torsional flow is inaccurate due to a more complex flow field. As the flow is no longer “controlled”, a rheological model must be assumed and the wall boundary conditions are redefined to allow for slip. A technique is described where, by examining the angular velocity response in very low torque experiments, it is possible to indirectly measure the wall slip velocity. The suspension is then tested at higher applied torques and different rheometer gaps. The results are integrated numerically to produce shear stress and shear rate values. This enables the measurement of true suspension bulk flow properties and wall slip velocity, with simple rheological models describing the observed complex rheograms.     

Investigation of third-grade non-Newtonian blood flow in arteries under periodic body acceleration using multi-step differential transformation method

In this paper, a non-Newtonian third-grade blood in coronary and femoral arteries is simulated analytically and numerically. The blood is considered as the thirdgrade non-Newtonian fluid under the periodic body acceleration motion and the pulsatile pressure gradient. The hybrid multi-step differential transformation method (Hybrid-MsDTM) and the Crank-Nicholson method (CNM) are used to solve the partial differential equation (PDE), and a good agreement between them is observed in the results. The effects of the some physical parameters such as the amplitude, the lead angle, and the body acceleration frequency on the velocity and shear stress profiles are considered. The results show that increasing the amplitude, A_g, and reducing the lead angle of body acceleration, φ, make higher velocity profiles on the center line of both arteries. Also, the maximum wall shear stress increases when A_g increases.     

局部狭窄血管中血液振荡流的速度和切应力

本文建立一种分析局部缓慢狭窄血管中血液振荡流的数学模型，给出了血液的轴向流速，径向流速和切应力的包含压力梯度项的解析表达式，并讨论了血管内由局部狭窄引起的压力梯度沿轴向变化的规律。文章以局部余弦狭窄为例进行数值计算，详细讨论上游均匀管段压力梯度的定常部分和不同次谐波对狭窄管段内流速和切应力的影响。数值结果表明，与均匀管情况相比，在狭窄段内，血液振荡流轴向流速无论平均值还是脉动幅值均明显增大，且径向流速不再为零。但径向流速仍远小于轴向流速。同时，切应力也不再仅由轴向流速梯度提供，径向流速梯度也将产生切应力，但是在计算管壁切向上的切应力时，径向流速梯度的贡献仍相当大。与均匀管管壁切应力沿流运方向保持恒定不同。狭窄管管壁切应力（平均值和脉动值）将随着狭窄高度的增大而增大，在狭窄最大高度处达到最大，因而沿流动方向产生了较大的切应力梯度。