Pulsatile flow of blood through a stenosed tube: effect of periodic body acceleration and a magnetic field

A proper understanding of the interactions of body acceleration and a magnetic field with blood flow could be useful in the diagnosis and treatment of some health problems. In the work reported in this paper we studied the pulsatile flow of blood through stenosed arteries, including the effects of body acceleration and a magnetic field. Blood is regarded as an electrically conducting, incompressible, couple-stress fluid in the presence of a magnetic field along the radius of the tube. The effects of the body acceleration and the magnetic field on the axial velocity, flow rate, and fluid acceleration were obtained analytically by use of the Hankel transform and the Laplace transform. Velocity variations under different conditions are shown graphically. The results have been compared with those from other theoretical models, and are in good agreement. Finally, our mathematical model gives a simple velocity expression for blood flow so it will help not only in the field of physiological fluid dynamics but will also help medical practitioners with elementary knowledge of mathematics.     

Inflation of an artery leading to aneurysm formation and rupture

Formation and rupture of aneurysms due to the inflation of an artery with collagen fibers distributed in two preferred directions, subjected to internal pressure and axial stretch are examined within the framework of nonlinear elasticity. A two layer tube model with a fiber-reinforced composite based incompressible anisotropic hyperelastic constitutive material is employed to model the stress-strain behavior of the artery wall with distributed collagen fibers. The artery wall takes up a uniform inflation deformation, and there are no aneurysms in the artery under the normal condition. But an aneurysm may be formed in arteries when the stiffness of the fibers is decreased to a certain value or the direction of the fibers is changed to a certain degree towards the circumferential direction. The aneurysm may expand to much large extent and become complex in shape. One portion of the aneurysm becomes highly distended as a bubble while the rest remains lightly inflated. The rupture of the aneurysm is discussed along with the distribution of stresses. Critical pressures and the rupture pressures are given for different collagen fiber orientations or stiffness. Furthermore, the stability of the solutions is discussed to explain the formation of aneurysm.     

Dynamic effects on the formation and rupture of aneurysms

Dynamic analysis of an axially stretched arterial wall with collagen fibers distributed in two preferred directions under a suddenly applied constant internal pressure along with the possibility of the formation and rupture of aneurysm are examined within the framework of nonlinear dynamics. A two layer tube model with the fiber-reinforced composite-based incompressible anisotropic hyper-elastic material is employed to model the mechanical behavior of the arterial wall. The maximum amplitudes and the phase diagrams are given by numerical computation of the differential relation. It is shown that the arterial wall undergoes nonlinear periodic oscillation and no aneurysms are formed under the normal condition. However, an aneurysm may be formed under such abnormal conditions as the stiffness of the fibers is deduced or the direction of the fibers is oriented towards the axial direction. Furthermore, the possibility for the rupture of aneurysm is discussed with the distribution of stresses.     

Wave Propagation in a Viscous Fluid Contained in an Orthotropic Elastic Tube

The problem of pressure wave propagation through a viscous fluid contained in an orthotropic elastic tube is considered in connection with arterial blood flow. Solutions to the fluid flow and elasticity equations are obtained for the presence of a reflected wave. Numerical results are presented for both isotropic and orthotropic elastic tubes. In particular, the pressure pulse, flow rate, axial fluid velocity, and wall displacements are plotted vs. time at various stations along the ascending aorta of man. The results indicate an increase in the peak value of the pressure pulse and a decrease in the flow rate as the pulse propagates away from the heart. Finally, the velocity of wave propagation depends mainly on the tangential modulus of elasticity of the arterial wall, and anisotropy of the wall accounts in part for the reduction of longitudinal movements and an increase in the hydraulic resistance.     

Axial prestretch and circumferential distensibility in biomechanics of abdominal aorta

Elastic arteries are significantly prestretched in an axial direction. This property minimises axial deformations during pressure cycle. Ageing-induced changes in arterial biomechanics, among others, are manifested via a marked decrease in the prestretch. Although this fact is well known, little attention has been paid to the effect of decreased prestretch on mechanical response. Our study presents the results of an analytical simulation of the inflation–extension behaviour of the human abdominal aorta treated as nonlinear, anisotropic, prestrained thin-walled as well as thick-walled tube with closed ends. The constitutive parameters and geometries for 17 aortas adopted from the literature were supplemented with initial axial prestretches obtained from the statistics of 365 autopsy measurements. For each aorta, the inflation–extension response was calculated three times, with the expected value of the initial prestretch and with the upper and lower confidence limit of the initial prestretch derived from the statistics. This approach enabled age-related trends to be evaluated bearing in mind the uncertainty in the prestretch. Despite significantly decreased longitudinal prestretch with age, the biomechanical response of human abdominal aorta changes substantially depending on the initial axial stretch was used. In particular, substituting the upper limit of initial prestretch gave mechanical responses which can be characterised by (1) low variation in axial stretch and (2) high circumferential distensibility during pressurisation, in contrast to the responses obtained for their weakly prestretched counterparts. The simulation also suggested the significant effect of the axial prestretch on the variation of axial stress in the pressure cycle. Finally, the obtained results are in accordance with the hypothesis that circumferential-to-axial stiffness ratio is the quantity relatively constant within this cycle.     

A stress-strain relation for a rat abdominal aorta

Assuming the arterial wall is homogeneous, incompressible, isotropic and elastic, a stress-strain relation has been presented for a rat's abdominal aorta. As an illustrating example, the problem of simultaneous inflation and the axial stretch of a cylindrical artery under physiological loading has been solved and then the material coefficients are determined by comparing theoretical results with the existing experiments. The result indicates that the maximum deviation between the theory and experiment for various pressure levels is 3.7% which seems to be a good approximation of theory to the experiments. The variation of circumferential stress and the incremental pressure modulus with inner pressure are also depicted in the work.     

Two-layered blood flow in stenosed tubes for different diseases

A two-fluid model for blood flow through a stenosed tube has been developed. The model consists of a core (suspension of RBCs) and peripheral plasma layer. The core is assumed to be represented by a polar fluid and the plasma layer by a Newtonian fluid. The flow is assumed to be steady and laminar, and the fluids incompressible. The flow variables are computed for normal blood and for the cases of polycythemia, plasma cell dyscrasias and for Hb SS diseases. Resistance to flow has been computed for different stenosis length and for different stenosis height. Shear stress distribution along the axial distance has been computed for different stenosis height. The impact of size effects (particle size to tube diameter) on blood diseases is discussed.     

Pulse waves in prestressed arteries

In order to better understand the effect of initial stress in blood flow in arteries, a theoretical analysis of wave propagation in an initially inflated and axially stretched cylindrical thick shell is investigated. For simplicity in the mathematical analysis, the blood is assumed to be an incompressible inviscid fluid while the arterial wall is taken to be an isotropic, homogeneous and incompressible elastic material. Employing the theory of small deformations superimposed on a large initial field the governing differential equations of perturbed solid motions are obtained in cylindrical polar coordinates. Considering the difficulty in obtaining a closed form solution for the field equations, an approximate power series method is utilized. The dispersion relations for the most general case of this approximation and for the thin tube case are thoroughly discussed. The speeds of waves propagating in an unstressed tube are obtained as a special case of our general treatment. It is observed that the speeds of both waves increase with increasing inner pressure and axial stretch.     

Contribution of elastin and collagen to the inflation response of the pig thoracic aorta: assessing elastin's role in mechanical homeostasis

This study was undertaken to understand elastin's role in the mechanical homeostasis of the arterial wall. The mechanical properties of elastin vary along the aorta, and we hypothesized this maintained a uniform mechanical environment for the elastin, despite regional variation in loading. Elastin's physiological loading was determined by comparing the inflation response of intact and autoclave purified elastin aortas from the proximal and distal thoracic aorta. Elastin's stretch and stress depend on collagen recruitment. Collagen recruitment started in the proximal aorta at systolic pressures (13.3 to 14.6 kPa) and in the distal at sub-diastolic pressures (9.3 to 10.6 kPa). In the proximal aorta collagen did not contribute significantly to the stress or stiffness, indicating that elastin determined the vessel properties. In the distal aorta, the circumferential incremental modulus was 70% higher than in the proximal aorta, half of which (37%) was due to a stiffening of the elastin. Compared to the elastin tissue in the proximal aorta, the distal elastin suffered higher physiological circumferential stretch (29%, P=0.03), circumferential stress (39%, P=0.02), and circumferential stiffness (37%, P=0.006). Elastin's physiological axial stresses were also higher (67%, P=0.003). These findings do not support the hypothesis that the loading on elastin is constant along the aorta as we expected from homeostasis.     

Developing steady laminar flow through uniform straight tubes with varying wall cross curvature

Numerical calculations are used to determine not only the wall shear stress but also the entry length in a laminar steady flow of an incompressible Newtonian fluid. The fluid is conveyed through rigid straight tubes with axially uniform cross sections, which mimic collapsed vessels. For each tube configuration, the "Navier-Stokes" equations are solved using the finite element method. The numerical tests are performed with the same value of the volume flow-rate whatever the tube configuration for three "Reynolds numbers". The wall shear stress is computed and determined along the axis of the tube, then the entry length is estimated by introducing two indexes by using: (i) the axial fluid velocity, and (ii) the wall shear stress. The results are analysed in order to exhibit the mechanical environment of cultured endothelial cells in the flow chamber for which the test conditions will be well-defined. For example, in a tube configuration where the opposite walls are in contact for which the inner perimeter and the area of the cross section are respectively given by 45 mm and 37.02 mm(2), the computed entry lengths with the criteria defined by (i) and (ii) are equals to about 118 and 126 mm, respectively for R(e0) = 500.     

Wave Propagation through a Viscous Incompressible Fluid Contained in an Initially Stressed Elastic Tube

To have a better understanding of the flow of blood in arteries a theoretical analysis of the pressure wave propagation through a viscous incompressible fluid contained in an initially stressed tube is considered. The fluid is assumed to be Newtonian. The tube is taken to be elastic and isotropic. The analysis is restricted to tubes with thin walls and to waves whose wavelengths are very large compared with the radius of the tube. It is further assumed that the amplitude of the pressure disturbance is sufficiently small so that nonlinear terms of the inertia of the fluid are negligible compared with linear ones. Both circumferential and longitudinal initial stresses are considered; however, their origins are not specified. Initial stresses enter equations as independent parameters. A frequency equation, which is quadratic in the square of the propagation velocity is obtained. Two out of four roots of this equation give the velocity of propagation of two distinct outgoing waves. The remaining two roots represent incoming waves corresponding to the first two waves. One of the waves propagates more slowly than the other. As the circumferential and/or longitudinal stress of the wall increases, the velocity of propagation and transmission per wavelength of the slower wave decreases. The response of the fast wave to a change in the initial stress is on the opposite direction.     

A piece-wise non-linear elastic stress expression of human and pig coronary arteries tested in vitro.

Eight human and nineteen pig unembalmed proximal left anterior descending and circumflex coronary arteries were subjected to linear volume changes (2 s ramp time) at three fixed axial extensions while immersed in a physiological saline bath at body temperature. Measured parameters included: lumen pressure, outside diameter, axial force, and axial extension. The deformations were measured using a video dimensional analyzer. The arteries were inflated to pressures well above the physiological range at each axial extension. A latex inner tube was placed inside of each specimen to prevent leakage, and its effects upon the measured stresses were corrected analytically. With this method, the average circumferential and axial stresses could be computed directly from the experimental data. In both directions the average stresses measured displayed two distinct regions: stresses occurring for small diameter changes (physiological pressures) and stresses occurring for large diameter changes (high pressures). The resulting average small strain and large strain stress components were curve-fit separately and, when reassembled, provided a piece-wise model of the stress response of coronary arteries over a wide range of inflation pressures and axial extensions.     

Influence of cardiac flow rate on turbulent shear stress from a prosthetic heart valve

Elevated turbulent shear stresses associated with sufficient exposure times are potentially damaging to blood constituents. Since these conditions can be induced by mechanical heart valves, the objectives of this study were to locate the maximum turbulent shear stress in both space and time and to determine how the maximum turbulent shear stress depends on the cardiac flow rate in a pulsatile flow downstream of a tilting disk valve. Two-component, simultaneous, correlated laser velocimeter measurements were recorded at four different axial locations and three different flow rates in a straight tube model of the aorta. All velocity data were ensemble averaged within a 15 ms time window located at approximately peak systolic flow over more than 300 cycles. Shear stresses as high as 992 dynes/cm2 were found 0.92 tube diameters downstream of the monostrut, disk valve. The maximum turbulent shear stress was found to scale with flow rate to the 0.72 power. A repeatable starting vortex was shed from the disk at the beginning of each cycle.     

Towards in vivo aorta material identification and stress estimation

This paper addresses the problem of constructing a mechanical model for the abdominal aorta and calibrating its parameters to in vivo measurable data. The aorta is modeled as a pseudoelastic, thick-walled, orthotropic, residually stressed cylindrical tube, subjected to an internal pressure. The model parameters are determined by stating a minimization problem for the model pressure and computing the optimal solution by a minimization algorithm. The data used in this study is in vivo pressure–diameter data for the abdominal aorta of a 24-year-old man. The results show that the axial, circumferential and radial stresses have magnitudes in the span 0 to 180 kPa. Furthermore, the results show that it is possible to determine model parameters directly from in vivo measurable data. In particular, the parameters describing the residual stress distribution can be obtained without interventional procedures.     

Effects of periodic accelerations on blood flow in arteries

This paper discusses the flow of blood in large artries under the influence of linear periodic acceleration. The governing equations and boundary conditions are established and analytical solutions for the velocity, fluid acceleration, bulk flow and shear stress are obtained. The results for these physical quantitites are computed for the case of an artery the size of a normal human aorta. It is found that the flow field variables are directly proportional to the external accelerating force. The behaviour of the velocity profile along the radial distance at different stages of times at fixed applied acceleration is also shown.     

Fluid-Structure Interaction Analysis of Pulsatile Flow within a Layered and Stenotic Aorta

In this paper, the hemodynamic characteristics of blood flow and stress distribution in a layered and stenotic aorta are investigated. By introducing symmetrical and unsymmetrical stenosis, the influence of stenosis morphology and stenotic ratio on the coupled dynamic responses of aorta is clarified. In the analysis, the in-vivo pulsatile waveforms and fully fluid–structure interaction (FSI) between the layered elastic aorta and the blood are considered. The results show that the fluid domain is abnormal in the stenotic aorta, and the whirlpool forms at the obstructed and downstream unobstructed regions. The maximum wall shear stresses appear at the throat of the stenosis. Downstream region appears low and oscillated shear stresses. In addition, along with the increase of the stenotic ratio, the amplitude of the maximum shear stress will be intensively increased and localized, and the sensitivity is also increased. In the aorta with unsymmetrical stenosis, the Von Mises stresses reach the peak value at the side with the surface protuberance, but they are reduced at the side with no protuberance. The sign variation of the layer interface shear stresses near the throat indicates the variation of the shear direction which increases the opportunity of shear damage at the transition plane. Moreover, the shear stress levels at the fluid-solid and intima-media interfaces are higher than that at the media-adventitia interface. The unsymmetrical stenosis causes higher stresses at the side with the surface protuberance than symmetrical one, but lower at the side with no protuberance. These results provide an insight in the influence of the stenosis, as well as its morphology, on the pathogenesis and pathological evolution of some diseases, such as arteriosclerosis and aortic dissection.     

Dynamics of flow in a mechanical heart valve: the role of leaflet inertia and leaflet compliance

In this work, we examine the dynamics of fluid flow in a mechanical heart valve when the solid inertia and leaflet compliance are important. The fluid is incompressible and Newtonian, and the leaflet is an incompressible neo-Hookean material. In the case of an inertialess leaflet, we find that the maximum valve opening angle and the time that the valve remains closed increase as the shear modulus of the leaflet decreases. More importantly, the regurgitant volume decreases with decreasing shear modulus. When we examined the forces exerted on the leaflet, we found that the downward motion of the leaflet is initiated by a vertical force exerted on its right side and, later on, by a vertical force exerted on the top side of the leaflet. In the case of solid inertia, we find that the maximum valve opening angle and the regurgitant volume are larger than in the case of an inertialess leaflet. These results highlight the importance of solid compliance in the dynamics of blood flow in a mechanical heart valve. More importantly, they indicate that mechanical heart valves with compliant leaflets may have smaller regurgitant volumes and smaller shear stresses than the ones with rigid leaflets.     

Axial nonuniformity of geometric and mechanical properties of mouse aorta is increased during postnatal growth

The hemodynamic conditions of aorta are relatively uniform prenatally and become more heterogeneous postnatally. Our objective was to quantify the heterogeneity of geometry and mechanical properties during growth and development. To accomplish this objective, we obtained a systematic set of data on the geometry and mechanical properties along the length of mouse aorta during postnatal development. C57BL/6 mice of ages 1-33 days were studied. The ascending aorta was cannulated in situ and preconditioned with several cyclic changes in pressure. We investigated the axial variations of geometry (diameter and length) and mechanical properties (stress-stain relation, elastic modulus and compliance) of the mouse aorta from the aortic valve to the common iliac. Our results show that the arterial blood pressure of mice increased from approximately 30 to 80 mmHg during the first 2 wk of life. The stretch ratio, diameter, wall (intima-media) thickness, and total lumen volume of mouse aorta increased with age. The aorta was transformed from a cylindrical tube at birth to a tapered structure during growth. Furthermore, we found the mechanical properties were fairly uniform along the length of the aorta at birth and become more nonuniform with age. We conclude that the rapid change of blood pressure and blood flow after birth alter the geometric and mechanical properties differentially along the length of the aorta. Hence, the axial nonuniformity of the aorta increases as the organ becomes more specialized during growth and development.     

Natural Convection Flow near a Vertical Plate that Applies a Shear Stress to a Viscous Fluid

The unsteady natural convection flow of an incompressible viscous fluid near a vertical plate that applies an arbitrary shear stress to the fluid is studied using the Laplace transform technique. The fluid flow is due to both the shear and the heating of the plate. Closed-form expressions for velocity and temperature are established under the usual Boussinesq approximation. For illustration purposes, two special cases are considered and the influence of pertinent parameters on the fluid motion is graphically underlined. The required time to reach the steady state in the case of oscillating shear stresses on the boundary is also determined.     

Low pressure condition of a lipid core in an eccentrically developed carotid atheromatous plaque: a static finite element analysis

Plaque ruptures in atherosclerotic carotid arteries cause cerebral strokes. Accumulation of lipoproteins in the deep intimal layer forms a lipid core (LC), whose progression may be enhanced by mechanical conditions on the arterial wall. In this study, we investigated the pressure conditions of a liquid LC through numerical simulations of a sliced segment finite element (FE) model and a three-dimensional (3D) symmetric FE model. A model of an LC filled with nearly incompressible fluid was compared with incompressible and soft neo-Hookean LC models in a static FE analysis. Material constants for a nonlinear hyperelastic model of the arterial wall were identified based on an inflation test using a tube specimen. The results from the FE analysis of a sliced segment model show an LC fluid pressure as low as 1.9 kPa at a blood pressure of 16 kPa. A neo-Hookean LC model with a Young’s modulus of 0.06 kPa produced an almost uniform pressure in the LC within an error of 1.3 %. The 3D model predicted a similar level of LC pressure. Such low fluid pressure in the LC region may enhance the infiltration of lipoproteins and other substances from the lumen and facilitate transport through microvessels from the adventitia to the LC.