A two-fluid model for pulsatile flow in catheterized blood vessels

The pulsatile flow of blood through a catheterized artery is analyzed, assuming the blood as a two-fluid model with the suspension of all the erythrocytes in the core region as a Casson fluid and the peripheral region of plasma as a Newtonian fluid. The resulting non-linear implicit system of partial differential equations is solved using perturbation method. The expressions for shear stress, velocity, flow rate, wall shear stress and longitudinal impedance are obtained. The variations of these flow quantities with yield stress, catheter radius ratio, amplitude, pulsatile Reynolds number ratio and peripheral layer thickness are discussed. It is observed that the velocity distribution and flow rate decrease, while, the wall shear, width of the plug flow region and longitudinal impedance increase when the yield stress increases. It is also found that the velocity increases, but, the longitudinal impedance decreases when the thickness of the peripheral layer increases. The wall shear stress decreases non-linearly, while, the longitudinal impedance increases non-linearly when the catheter radius ratio increases. The estimates of the increase in the longitudinal impedance are considerably lower for the present two-fluid model than those of the single-fluid model.     

Oxygen transfer in human carotid artery bifurcation

Arterial bifurcations are places where blood flow may be disturbed and slow recirculation flow may occur. To reveal the correlation between local oxygen transfer and atherogenesis, a finite element method was employed to simulate the blood flow and the oxygen transfer in the human carotid artery bifurcation. Under steady-state flow conditions, the numerical simulation demonstrated a variation in local oxygen transfer at the bifurcation, showing that the convective condition in the disturbed flow region may produce uneven local oxygen transfer at the blood/wall interface. The disturbed blood flow with formation of slow eddies in the carotid sinus resulted in a depression in oxygen supply to the arterial wall at the entry of the sinus, which in turn may lead to an atherogenic response of the arterial wall, and contribute to the development of atherosclerotic stenosis there. The project supported by the National Natural Science Research Council of China (10632010, 10572017, 30670517).     

A two-fluid model of non-Newtonian blood flow induced by peristaltic waves

The problem of blood flow induced by peristaltic waves in a uniform small diameter tube has been investigated. Blood has been represented by a two-fluid model consisting of a core region of suspension of all the erythrocytes, assumed to be a Casson fluid, and a peripheral layer of plasma as a Newtonian fluid. The expressions for dimensionless pressure drop and friction force have been obtained. The results obtained in the analysis have been evaluated numerically and discussed briefly. The significance of the present model over the existing models has been pointed out by comparing the results with other theories both analytically and numerically.     

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.     

冠状动脉对称双路移植管中脉动血流的数值模拟

冠状动脉旁路移植管搭桥术后,常会产生血管再狭窄,导致手术失败,。这与移植管的几何结构及血流动力学是密切相关的。作为改进措施,作者提出了采用对称双路搭桥的设想。本文利用有限元分析方法,对冠状动脉搭桥术中对称双路移植管内的生理流动进行了数值仿真。计算了缝合区附近的流场、壁面切应力、压力等血流动力学因素在心动周期内的时空分布情况。计算结果表明,对称双路搭桥具有较好的血流动力学,可以改善血管流场状况和减轻再狭窄发生。这对临床手术计划是很有帮助和指导意义的。     

颈动脉分支的血流动力学数值模拟

采用有限元法数值模拟颈动脉分支的血流动力学。根据在体测量的实际尺寸来构造颈动脉分支的几何模型,以保持模型的解剖精确度;利用在体测量的颈内动脉和颈外动脉流量波形以及主颈动脉的压力波形来确定数值计算的边界条件,以保持数值计算的生理真实性。关注的重点是颈动脉窦内的局部血流形态、二次流和壁面剪应力。在心脏收缩的减速期和舒张期的某些时刻,颈动脉窦中部外侧壁面附近产生了流动分离,形成了一个低速回流区。该流动分离是瞬态的,导致了壁面剪应力的振荡,其振荡范围在-2～6dyn／cm^2之间。同时,颈动脉窦中部横截面内的二次流存在于整个心动周期,最大的二次流速度为同时刻轴向速度平均值的1／3左右。     

窦部对称狭窄对颈动脉内流场影响的数值研究

以颈动脉分岔血管为例，采用数值方法研究了窦部环缩狭窄之后的流场分布情况，并和正 常血管情况下的流场分布进行了比较. 结果表明，采用环缩方式给颈动脉分岔血管施加对称 的狭窄改变了颈动脉窦内流场，特别是壁面剪应力的分布规律. 低剪应力区出现在狭窄段之 后的窦内，并且沿整个周向均匀分布. 根据低剪应力和动脉粥样硬化的关系，指出： 若人为地给颈动脉窦内施加对称狭窄，则脂质沉积将在狭窄下游的窦内沿周向轴对称 发展. 为了更真实地反映颈动脉窦内的狭窄，建议根据动脉血管中的实际狭窄情况，采用非 对称的狭窄分布模式.     

Multidimensional modeling of the stenosed carotid artery: A novel CAD approach accompanied by an extensive lumped model

This study describes a multidimensional 3D/lumped parameter(LP) model which contains appropriate inflow/outflow boundary conditions in order to model the entire human arterial trees. A new extensive LP model of the entire arterial network(48 arteries) was developed including the effect of vessel diameter tapering and the parameterization of resistance, conductor and inductor variables. A computer aided-design(CAD) algorithm was proposed to effciently handle the coupling of two or more 3D models with the LP model, and substantially lessen the coupling processing time. Realistic boundary conditions and Navier–Stokes equations in healthy and stenosed models of carotid artery bifurcation(CAB) were used to investigate the unsteady Newtonian blood flow velocity distribution in the internal carotid artery(ICA). The present simulation results agree well with previous experimental and numerical studies. The outcomes of a pure LP model and those of the coupled 3D healthy model were found to be nearly the same in both cases. Concerning the various analyzed 3D zones, the stenosis growth in the ICA was not found as a crucial factor in determining the absorbing boundary conditions.This paper demonstrates the advantages of coupling local and systemic models to comprehend physiological diseases of the cardiovascular system.     

The numerical simulation of turbulent flows of Newtonian and a sort of non-Newtonian fluid in 180-degree square sectioned bend

Using k-εmodel of turbulence and measured wall functions.turbulent flows ofNewtonian(pure water)and a sort of non-Newtonian fluid(dilute,drag-reduction solutionof polymer in a180-degree curved bend were simulated numerically.The calculated resultsagreed well with the measured velocity profiles.On the basis of calculation andmeasurement,appropriateness of turbulence model to complicated flow in which the large-scale vortex exists was analyzed and discussed.     

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.     

Multiphase non‐Newtonian effects on pulsatile hemodynamics in a coronary artery

Hemodynamic stresses are involved in the development and progression of vascular diseases. This study investigates the influence of mechanical factors on the hemodynamics of the curved coronary artery in an attempt to identify critical factors of non‐Newtonian models. Multiphase non‐Newtonian fluid simulations of pulsatile flow were performed and compared with the standard Newtonian fluid models. Different inlet hematocrit levels were used with the simulations to analyze the relationship that hematocrit levels have with red blood cell (RBC) viscosity, shear stress, velocity, and secondary flow. Our results demonstrated that high hematocrit levels induce secondary flow on the inside curvature of the vessel. In addition, RBC viscosity and wall shear stress (WSS) vary as a function of hematocrit level. Low WSS was found to be associated with areas of high hematocrit. These results describe how RBCs interact with the curvature of artery walls. It is concluded that although all models have a good approximation in blood behavior, the multiphase non‐Newtonian viscosity model is optimal to demonstrate effects of changes in hematocrit. They provide a better stimulation of realistic blood flow analysis. Copyright © 2008 John Wiley & Sons, Ltd.     

Unified analysis of liquefaction and the ground flow phenomenon

Loose saturated sand behaves as a solid before liquefaction but as a fluid when the excess pore water pressure equals the initial confining stress, after which it recovers its strength. A simple constitutive equation for loose saturated sand was developed to express the phase transformation between a solid and fluid during liquefaction and the ground flow phenomenon. This constitutive equation was used for a shaking table test, and its applicability was investigated by comparing numerical and experimental results Published in Prikladnaya Mekhanika, Vol. 43, No. 8, pp. 129–144, August 2007. An erratum to this article is available at .     

Effects of renal artery stenosis on realistic model of abdominal aorta and renal arteries incorporating fluid-structure interaction and pulsatile non-Newtonian blood flow

The effects of the renal artery stenosis(RAS) on the blood flow and vesselwalls are investigated.The pulsatile blood flow through an anatomically realistic model ofthe abdominal aorta and renal arteries reconstructed from CT-scan images is simulated,which incorporates the fluid-structure interaction(FSI).In addition to the investigationof the RAS effects on the wall shear stress and the displacement of the vessel wall,it isdetermined that the RAS leads to decrease in the renal mass flow.This may cause theactivation of the renin-angiotension system and results in severe hypertension.     

Numerical simulation of the non-Newtonian blood flow through a mechanical aortic valve

This work focuses on the comparison between Newtonian and non-Newtonian blood flows through a bileaflet mechanical heart valve in the aortic root. The blood, in fact, is a concentrated suspension of cells, mainly red blood cells, in a Newtonian matrix, the plasma, and consequently its overall behavior is that of a non-Newtonian fluid owing to the action of the cells’ membrane on the fluid part. The common practice, however, assumes the blood in large vessels as a Newtonian fluid since the shear rate is generally high and the effective viscosity becomes independent of the former. In this paper, we show that this is not always the case even in the aorta, the largest artery of the systemic circulation, owing to the pulsatile and transitional nature of the flow. Unexpectedly, for most of the pulsating cycle and in a large part of the fluid volume, the shear rate is smaller than the threshold level for the blood to display a constant effective viscosity and its shear thinning character might affect the system dynamics. A direct inspection of the various flow features has shown that the valve dynamics, the transvalvular pressure drop and the large-scale features of the flow are very similar for the Newtonian and non-Newtonian fluid models. On the other hand, the mechanical damage of the red blood cells (hemolysis), induced by the altered stress values in the flow, is larger for the non-Newtonian fluid model than for the Newtonian one.     

动脉分岔血管内膜增生过程的数值模拟

内膜增生从发生到阻塞血管是一个复杂的变化过程,在这个过程中,内膜的增生、血管腔体形状的改变和血流动力学之间是相互影响的。为了研究这些变化,本文提出一种单元填充方法数值模拟了三维颈动脉分岔血管在低切应力作用下血管内膜增生的过程。该方法既可以克服节点移动方法所不可避免的内膜增生的不连续性,也可以避免网格重划分的困难。结果发现,如果单纯以切应力阈值作为内膜增生的判据,低切应力的作用将无法导致血管完全阻塞,但内膜增生和血流动力学之间的相互影响是可以通过数值方法进行模拟的。在本数值模拟中,内膜增生的过程分为"增厚"(先)和"扩展"(后)两个阶段,最大狭窄率为34.4%,发生在距血管分岔5mm处动脉窦的外侧壁面。其发生位置和形状与临床观察吻合。     

Basic equations of turbulent flow for variable density and variable viscosity Newtonian fluid in open channel

I.Intr0ducti0nManyscholarshaveresearchedtheturbulenceproblemsofnon-compressiblefluidll'2l,e.g.inRef[1],DouGu0renhaspresentedtheturbulentmotionaldifferentialequationsofconstantdensityandconstantvisc0sityfluidin0penchannel;inRef[3],XiaoTianduohasstudiedthet…     

Numerical simulation of mass transfer to micropolar fluid flow past a stenosed artery

A mathematical model of unsteady non‐Newtonian blood flow together with the mass transfer through constricted arteries has been developed. The mass transport refers to the movement of atherogenic molecules, i.e. blood‐borne components, such as low‐density lipoproteins from flowing blood into the arterial walls or vice versa. The flowing blood is represented as the suspension of all erythrocytes assumed to be Eringen's micropolar fluid and the arterial wall is considered to be rigid having cosine‐shaped stenosis in its lumen. The mass transfer to blood is controlled by the convection–diffusion equation. The governing equations of motion accompanied by the appropriate choice of the boundary conditions are solved numerically by Marker and Cell method and the results obtained are checked for numerical stability with the desired degree of accuracy. The quantitative analysis carried out finally includes the respective profiles of the flow‐field and the mass concentration along with their distributions over the entire arterial segment as well. The key factors, such as the wall shear stress and Sherwood number, are also examined for further quantitative insight into the flow and the mass transport phenomena through arterial stenosis. The present results show consistency with several existing results in the literature which substantiate sufficiently to validate the applicability of the model under consideration. Copyright © 2010 John Wiley & Sons, Ltd.     

峰后岩石非Darcy渗流的分岔行为研究

煤矿采动围岩大多处于峰后应力状态或破碎状态，其渗流一般不符合Darcy定律，为非Darcy渗流系统．峰后岩石非Darcy渗流系统的失稳和分岔是煤矿突水和煤与瓦斯突出动力灾害发生的根源．文中用谱截断方法建立了Ahmed-Sunada型非Darcy渗流系统的降阶动力学方程，再由变量代换得到以无量纲变量表示的平衡态附近的演化方程，分析了系统的分岔条件，给出了系统的各种吸引子图案，并结合采矿工程实际，用非线性数学的观点揭示了煤矿突水和煤与瓦斯突出的机理．研究表明：当非Darcy渗流系统渗流特性和边界压力的初始值满足一定条件时，系统由平衡转向不稳定，即存在跨临界Hopf分岔和切分岔，并且，系统的动力学响应不随渗透特性连续变化，即该系统存在突变性．     

Tomographic particle image velocimetry investigation of the flow in a modeled human carotid artery bifurcation

Hemodynamic forces within the human carotid artery are well known to play a key role in the initiation and progression of vascular diseases such as atherosclerosis. The degree and extent of the disease largely depends on the prevailing three-dimensional flow structure and wall shear stress (WSS) distribution. This work presents tomographic PIV (Tomo-PIV) measurements of the flow structure and WSS in a physiologically accurate model of the human carotid artery bifurcation. The vascular geometry is reconstructed from patient-specific data and reproduced in a transparent flow phantom to demonstrate the feasibility of Tomo-PIV in a complex three-dimensional geometry. Tomographic reconstruction is performed with the multiplicative line-of-sight (MLOS) estimation and simultaneous multiplicative algebraic reconstruction (SMART) technique. The implemented methodology is validated by comparing the results with Stereo-PIV measurements in the same facility. Using a steady flow assumption, the measurement error and RMS uncertainty are directly inferred from the measured velocity field. It is shown that the measurement uncertainty increases for increasing light sheet thickness and increasing velocity gradients, which are largest near the vessel walls. For a typical volume depth of 6 mm (or 256 pixel), the analysis indicates that the velocity derived from 3D cross-correlation can be measured within ±2% of the maximum velocity (or ±0.2 pixel) near the center of the vessel and within ±5% (±0.6 pixel) near the vessel wall. The technique is then applied to acquire 3D-3C velocity field data at multiple axial locations within the carotid artery model, which are combined to yield the flow field and WSS in a volume of approximately 26 mm × 27 mm × 60 mm. Shear stress is computed from the velocity gradient tensor and a method for inferring the WSS distribution on the vessel wall is presented. The results indicate the presence of a complex and three-dimensional flow structure, with regions of flow separation and strong velocity gradients. The WSS distribution is markedly asymmetric confirming a complex swirling flow structure within the vessel. A comparison of the measured WSS with Stereo-PIV data returns an acceptable agreement with some differences in stress magnitude.     

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube

The effect of yield stress on the flow characteristics of a Casson fluid in a homogeneous porous medium bounded by a circular tube is investigated by employing the Brinkman model to account for the Darcy resistance offered by the porous medium. The non-linear coupled implicit system of differential equations governing the flow is first transformed into suitable integral equations and are solved numerically. Analytical solution is obtained for a Newtonian fluid in the case of constant permeability, and the numerical solution is verified with that of the analytic solution. The effect of yield stress of the fluid and permeability of the porous medium on shear stress and velocity distributions, plug flow radius and flow rate are examined. The minimum pressure gradient required to start the flow is found to be independent of the permeability of the porous medium and is equal to the yield stress of the fluid.