On the assumption of Saint-Venant's problem

In this paper we obtain uniquely the solution of Saint- Venant’s problem under theassumption of ∂^m/∂z^mσ_z=0(m≥2) for noncircular prismatic bars.     

The steady/pulsatile flow and macromolecular transport in T-bifurcation blood vessels

IntroductionThefluiddynamicsofbloodcyclesystemplayanimportantroleinthepathogenesisofatherosclerosis.ThephysiologicalanatomyfoundthattheatherosclerosisappearsoftenatthebifurcationorcurvedflowareatoallkindsofRefs.[1 ,2 ] .Theshearstressvariesgreatlyinthoseareaandinfluencesthemacromoleculartransportacrossthebloodwall[3,4].Thus ,theinvestigationoftheflowandmacromoleculartransportinthesecomplexbloodvesselaandthecorrelationbetweenthemareinterestingtoresearchers.Intheseaspect,Liepschstudiedtheflowi…     

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.     

Boundary layer development of pulsatile blood flow in a tapered vessel

Assuming that the tapered angle is small,the problems of developing flow under unsteady oscillatory condition are studied in this paper.The formula of velocity distribution is obtained.The analyses for the results show that the blood flow in a converging tapered vessel remains a developing flow throughout the length,and the effects of tapered angle on the developing flow are increased with the increment of the tapered angle.     

The effects of blood viscoelasticity on the pulse wave in arteries

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈａｔｔｈｅｂｌｏｏｄｈａｓｖｉｓｃｏｅｌａｓｔｉｃｐｒｏｐｅｒｔｉｅｓｉｓａｗｅｌｌ＿ｋｎｏｗｎｆａｃｔ.ＴｈｅｒｅｓｅａｒｃｈｅｓｆｏｒｔｈｅｂｌｏｏｄｖｉｓｃｏｅｌａｓｔｉｃｉｔｙｂｙＧ .Ｂ .Ｔｈｕｒｓｔｏｎ[1～4]ａｎｄＳ .Ｃｈｉｅｎ[5 ]ｓｈｏｗｔｈａｔｔｈｅｂｌｏｏｄｎｏｔｏｎｌｙａｐｐｅａｒｓｔｈｅｖｉｓｃｏｅｌａｓｔｉｃｉｔｙｉｎｖａｒｉｏｕｓｏｓｃｉｌｌａｔｏｒｙｂｌｏｏｄｆｌｏｗｓ,ｂｕｔａｌｓｏｈａｓｑｕｉｔｅｓｔｒｏｎｇｅｌａｓｔｉｃｉｔｙｉｎｓｏｍｅ…     

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.     

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 study of effects of pulsatile amplitude for transitional turbulent pulsatile flow in pipes with ring‐type constrictions

The effects of pulsatile amplitude on sinusoidal transitional turbulent flows through a rigid pipe in the vicinity of a sharp‐edged mechanical ring‐type constriction have been studied numerically. Pulsatile flows were studied for transitional turbulent flow with Reynolds number (Re) of the order of 10⁴, Womersley number (Nw) of the order of 50 with a corresponding Strouhal number (St) of the order of 0.04. The pulsatile flow considered is a sinusoidal flow with dimensionless amplitudes varying from 0.0 to 1.0. Transitional laminar and turbulent flow characteristics in an alternative manner within the pulsatile flow fields were observed and studied numerically. The flow characteristics were studied through the pulsatile contours of streamlines, vorticity, shear stress and isobars. It was observed that fluid accelerations tend to suppress the development of flow disturbances. All the instantaneous maximum values of turbulent kinetic energy, turbulent viscosity, turbulent shear stress are smaller during the acceleration phase when compared with those during deceleration period. Various parametric equations within a pulsatile cycle have also been formulated through numerical experimentations with different pulsatile amplitudes. In the vicinity of constrictions, the empirical relationships were obtained for the instantaneous flow rate (Q), the pressure gradient (dp/dz), the pressure loss (P_loss), the maximum velocity (V_max), the maximum vorticity (ζ_max), the maximum wall vorticity (ζ_w,max), the maximum shear stress (τ_max) and the maximum wall shear stress (τ_w,max). Elliptic relation was observed between flow rate and pressure gradient. Quadratic relations were observed between flow rate and the pressure loss, the maximum values of shear stress, wall shear stress, turbulent kinematic energy and the turbulent viscosity. Linear relationships exist between the instantaneous flow rate and the maximum values of vorticity, wall vorticity and velocity. The time‐average axial pressure gradient and the time average pressure loss across the constriction were observed to increase linearly with the pulsatile amplitude. Copyright © 1999 John Wiley & Sons, Ltd.     

ANALYSE OF PULSE WAVE PROPAGATION IN ARTERIES

Based upon the blood vessel of being regarded as the elasticity tube, and that the tissue restricts the blood vessel wall, the rule of pulse wave propagation in blood vessel was studied. The viscosity of blood, the elastic modulus of blood vessel, the radius of tube that influenced the pulse wave propagation were analyzed. Comparing the result that considered the viscosity of blood with another result that did not consider the viscosity of blood, we finally discover that the viscosity of blood that influences the pulse wave propagation can not be neglected; and with the accretion of the elastic modulus the speed of propagation augments and the press value of blood stream heightens; when diameter of blood vessel reduces, the press of blood stream also heightens and the speed of pulse wave also augments. These results will contribute to making use of the information of pulse wave to analyse and auxiliarily diagnose some causes of human disease.     

弯曲动脉的血流动力学数值分析

利用计算流体力学的理论和方法对弯曲动脉中的血流动力学进行数值分析，是研究心血管疾病流体动力学机理的一种行之有效的方法。本文将升主动脉、主动脉弓和降主动脉联系起来作为弯曲动脉几何模型，给出了血液流动的边界条件以及计算条件。根据生理脉动流条件，对狗的弯曲动脉几何模型内发展中的血液流动进行了有限元数值模拟，并利用可视化方法对血液流动的轴向速度、二次流、壁面切应力等计算结果进行了分析。研究结果表明，在弯管内侧壁处，同时存在主流方向和二次流方向的回流，此处容易形成涡流。弯管内侧壁比外侧壁的壁面切应力具有更强的脉动性。     

An unsteady analysis of non-Newtonian blood flow through tapered arteries with a stenosis

The problem of non-Newtonian and nonlinear blood flow through a stenosed artery is solved numerically where the non-Newtonian rheology of the flowing blood is characterised by the generalised Power-law model. An improved shape of the time-variant stenosis present in the tapered arterial lumen is given mathematically in order to update resemblance to the in vivo situation. The vascular wall deformability is taken to be elastic (moving wall), however a comparison has been made with nonlinear visco-elastic wall motion. Finite difference scheme has been used to solve the unsteady nonlinear Navier-Stokes equations in cylindrical coordinates system governing flow assuming axial symmetry under laminar flow condition so that the problem effectively becomes two-dimensional. The present analytical treatment bears the potential to calculate the rate of flow, the resistive impedance and the wall shear stress with minor significance of computational complexity by exploiting the appropriate physically realistic prescribed conditions. The model is also employed to study the effects of the taper angle, wall deformation, severity of the stenosis within its fixed length, steeper stenosis of the same severity, nonlinearity and non-Newtonian rheology of the flowing blood on the flow field. An extensive quantitative analysis is performed through numerical computations of the desired quantities having physiological relevance through their graphical representations so as to validate the applicability of the present model.     

Effect of the magnetic field on the pulsatile flow through a rigid tube

I.IntroductionThefunctionofanaPbliedmagneticfield\onhumanbodieshasbeenknownandservedasakindofthemedicaltreatmentforalongtime.ThemagneticstoneusedasamedicineisrdeordedinShennongMaterlbMedica,'which.isthefirstChinesemedicinebookandwaswrittenin200A.D..ItwasalsorecordedtyancientGreecethatholdingmagneticstQnesinhandsandfeetcouldrelievethepainandspasml'l.However,beinglimitedbythelevelQfproductionandmedicaltreatmentsatthattime,themagnetictreatmentwasmainlyusedforlocallydephlogisticating,relievin…     

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.     

The influence of low-frequency varying magnetic field on the pulsatile flow in a rigid round tube

ProjcctsupportedbytheNatiol1alNaturalScienceFoundatio11ofChinaThemagnatictherapy,asaneffccti\'emedicaltreatment,hasabroadapp1icatiQn.lnrecentyears,plentyofanimalandclinicalresearcheshavebeendoneandprovedthatmagneticfieldhassuchpositivecurativeefttctsasant…     

Performance modeling and analysis of blood flow in elastic arteries

Nomenclature　　τ　　wallshearstressγshearrateτy yieldstressηc Cassonviscosityktheconsistencyindexnnon_Newtonianindexτp shearstressofthepthelementωangularvelocityRvessel’sradiusCwavespeedM　　magneticparameter (Hartmannnumber)u,w　velocitycomponentinther_andz_directions,respectivelyP　　pressureα　　unsteadinessparameter k , R meanparametersTp relaxationtimeofthepthelementρ densityIntroductionTheimportancetoatherogenesisofarterialflowphenomenasuchasflowseparation ,recirculationands…     

A numerical study of pulsatile laminar flows in a pipe with a ring-type constriction

Numerical simulations have been carried out to study pulsatile laminar flows in a pipe with an axisymmetric ringtype constriction. Three types of pulsatile flows were investigated, namely a physiological flow, a pure sinusoidal flow and a non-zero mean velocity sinusoidal flow. The laminar flow governing equations were solved by the SIMPLE algorithm on a non-staggered grid and a modified Crank-Nicolson approximation was used to discretrize the momentum equations with respect to time. The maximum flow Reynolds numer (Re) is 100. The Womersley number (N_w) ranges from 0 to 50, with the corresponding Strouhal number (St) ranging from 0 to 3·98. The constriction opening ratio (d/D) and thickness ratio (h/D) are fixed at 0·5 and 0·1 respectively. Within the time period investigated, all these pulsatile flows include both forward and backward flows. The unsteady recirculation region and the recirculation points change in size and location with time. For N_w ≤ 1 and St≤ 1·56 x 10^?3 the three pulsatile flows have the same simple relation between the instantaneous flow rate and pressure loss (Δp) across the constriction and the pressure gradient in the axial direction (dp/dz) in the fully developed flow region. The phase angles between the flow rate and pressure loss and the pressure gradient are equal to zero. With increasing N_w and St, the phase angle between the flow rate and the dp/dz becomes larger and has its maximum value of 90° at N_w = 50 and St = 3·98. The three pulsatile flows also show different relations between the flow rate and the pressure gradient. The pure sinusoidal flow has the largest maximum pressure gradient and the non-zero mean velocity sinusoidal flow has the smallest. For larger N_w and St the fully developed velocity profiles in the fully developed flow region have a smaller velocity gradient along the radial direction in the central region. The maximum recirculation length increases for N_w ranging from 0 to 4·2, while this length becomes very small at N_w = 50 and St = 3·98. The deceleration tends to enlarge the recirculation region and this effect appears for N_w ≥ 3 and St ≥ 1·43×10^?2. Linear relations exist between the flow rate and the instantaneous maximum values of velocity, vorticity and shear stress.     

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.     

An experimental investigation of pulsatile flow through a smooth constriction

Measurements of flow disturbances in the downstream region of modeled stenoses in a rigid tube, with upstream pulsatile flow are reported. Experiments were conducted over physiologically relevant mean Reynolds numbers of 600; based on the tube diameter and the time-averaged value of upstream centerline velocity. Contoured constrictions with 25%, 50% and 75% area reductions were investigated and velocity data were obtained from ensemble averaging techniques (phase-locked waveform). Experimental data over extensive spatial regions of poststenotic fields were taken, employing a two-component laser Doppler velocimeter LDV. Constant time sampling techniques for performing data or frequency analyses were used to avoid velocity bias and to study the evolution of poststenotic flow disturbances. It is found that different types of flow disturbances exist downstream of the constriction. Data analysis methods with the aid of flow visualization allow accurate classification of the disturbances which are sensitive indicators of mild to moderate constrictions. Although the present study was motivated by a biological situation, sufficient data were reported in detail that they may also be used by investigators working in computational fluid dynamics.     

BLOOD FLOW AND MACROMOLECULAR TRANSPORT IN CURVED BLOOD VESSELS

A numerical analysis of the steady/pulsatile flow and macromolecular (such as LDL and Albumin) transport in curved blood vessels was carried out. The computational results predict that the vortex of the secondary flow is time-dependent in the aortic arch. The concentration of macromolecule concentrates at the region of sharp curve, and the wall concentration at the outer part is higher than that at the inner part. Atherosclerosis and thrombosis are prone to develop in such regions with sharp flow.     

Well‐balanced discontinuous Galerkin method and finite volume WENO scheme based on hydrostatic reconstruction for blood flow model in arteries

The blood flow model in arteries admits the steady state solutions, for which the flux gradient is nonzero, and is exactly balanced by the source term. In this paper, by means of hydrostatic reconstruction, we construct a high order discontinuous Galerkin method, which exactly preserves the dead‐man steady state, which is characterized by a discharge equal to zero (analogue to hydrostatic equilibrium). Moreover, the method maintains genuine high order of accuracy. Subsequently, we apply the key idea to finite volume weighted essentially non‐oscillatory schemes and obtain a well‐balanced finite volume weighted essentially non‐oscillatory scheme. Extensive numerical experiments are performed to verify the well‐balanced property, high order accuracy, as well as good resolution for smooth and discontinuous solutions.