An intergral method for the analysis of flow in arterial stenoses

An approximate solution is presented to the problem of incompressible flow through an axisymmetric constriction. The geometry is intended to simulate an arterial stenosis, and the solution is applicable to both mild and severe stenoses for Reynolds numbers below transition. Theoretical results obtained for specific geometries are given for the velocity distribution, pressure drop, wall shearing stress, and separation phenomena. These results reveal the significant alterations in flow caused by a stenosis. Experiments using model stenoses are described and compared with the theoretical results. Theoretical predictions of pressure drop and separation characteristics are in reasonably good agreement with the experimental observations.     

Physiological flow simulation in residual human stenoses after coronary angioplasty

To evaluate the local hemodynamic implications of coronary artery balloon angioplasty, computational fluid dynamics (CFD) was applied in a group of patients previously reported by [Wilson et al. (1988), 77, pp. 873-885] with representative stenosis geometry post-angioplasty and with measured values of coronary flow reserve returning to a normal range (3.6 +/- 0.3). During undisturbed flow in the absence of diagnostic catheter sensors within the lesions, the computed mean pressure drop delta p was only about 1 mmHg at basal flow, and increased moderately to about 8 mmHg for hyperemic flow. Corresponding elevated levels of mean wall shear stress in the midthroat region of the residual stenoses, which are common after angioplasty procedures, increased from about 60 to 290 dynes/cm2 during hyperemia. The computations (Ree approximately equal to 100-400; alpha e = 2.25) indicated that the pulsatile flow field was principally quasi-steady during the cardiac cycle, but there was phase lag in the pressure drop-mean velocity (delta p - u) relation. Time-averaged pressure drop values, delta p, were about 20 percent higher than calculated pressure drop values, delta ps, for steady flow, similar to previous in vitro measurements by Cho et al. (1983). In the throat region, viscous effects were confined to the near-wall region, and entrance effects were evident during the cardiac cycle. Proximal to the lesion, velocity profiles deviated from parabolic shape at lower velocities during the cardiac cycle. The flow field was very complex in the oscillatory separated flow reattachment region in the distal vessel where pressure recovery occurred. These results may also serve as a useful reference against catheter-measured pressure drops and velocity ratios (hemodynamic endpoints) and arteriographic (anatomic) endpoints post-angioplasty. Some comparisons to previous studies of flow through stenoses models are also shown for perspective purposes.     

Effects of surface irregularities on flow resistance in differently shaped arterial stenoses

The combined influence of an asymmetric shape and surface irregularities has been explored in a computational study of flow through arterial stenoses with 48% areal occlusion. Contrary to the conclusion of an earlier investigation, namely that the resistance to laminar flow through a stenosed artery is being reduced in the presence of surface irregularities, the present predictions demonstrate that the flow resistance is practically unaffected by surface irregularities at low Reynolds numbers, whereas an excess pressure drop up to 10% above that for a smooth stenosis is observed for higher Reynolds numbers. For a given areal occlusion, the flow resistance is reduced with increasing degree of stenosis asymmetry and this effect may more than outweigh the influence of surface irregularities. This effect is moreover prevailing throughout the entire range of Reynolds numbers considered.     

Effects of elastic property of the wall on flow characteristics through arterial stenoses

Hemodynamic characteristics of blood flow through arterial stenoses are numerically investigated in this work. The blood is assumed as a Newtonian fluid and the pulsatile nature of flow is modeled by using measured values of the flowrate and pressure for the canine femoral artery. An isotropic elastic and incompressible material is assumed for the wall at each axial section, but a non-uniform distribution of the shear modulus in axial direction is used to model the high stiffness of the wall at the stenosis location. Full Navier equations for a thick wall are used as the governing equations for the wall displacements. A continuous grid extending over the flow field and the wall is considered and governing equations are transformed for use in the computational domain. Discretized forms of the transformed wall and flow equations, which are coupled through the boundary conditions at their interface, are obtained by control volume method and simultaneously solved using the well-known SIMPLER algorithm. To study the effects of wall deformability, solutions are obtained for both rigid and elastic walls. The results indicate that deformability of the wall causes an increase in the time average of pressure drop, but a decrease in the maximum wall shear stress. Displacement and stress distributions in the wall are presented.     

Computer simulation of intraventricular flow and pressure gradients during diastole

A two-dimensional axisymmetric computer model is developed for the simulation of the filling flow in the left ventricle (LV). The computed results show that vortices are formed during the acceleration phases of the filling waves. During the deceleration phases these are amplified and convected into the ventricle. The ratio of the maximal blood velocity at the mitral valve (peak E velocity) to the flow wave propagation velocity (WPV) of the filling wave is larger than 1. This hemodynamic behavior is also observed in experiments in vitro (Steen and Steen, 1994, Cardiovasc. Res., 28, pp. 1821-1827) and in measurements in vivo with color M-mode Doppler echocardiography (Stugaard et al., 1994, J. Am. Coll. Cardiol., 24, 663-670). Computed intraventricular pressure profiles are similar to observed profiles in a dog heart (Courtois et al., 1988, Circulation, 78, pp. 661-671). The long-term goal of the computer model is to study the predictive value of noninvasive parameters (e.g., velocities measured with Doppler echocardiography) on invasive parameters (e.g., pressures, stiffness of cardiac wall, time constant of relaxation). Here, we show that higher LV stiffness results in a smaller WPV for a given peak E velocity. This result may indicate an inverse relationship between WPV and LV stiffness, suggesting that WPV may be an important noninvasive index to assess LV diastolic stiffness, LV diastolic pressure and thus atrial pressure (preload).     

Disorder distal to modeled stenoses in steady and pulsatile flow

Measurements of the velocity and energy spectra were made in the distal region of modeled stenoses in a rigid tube with both steady and pulsatile water flows. Reynolds numbers of 318–2540 and a pulsatile flow frequency parameter of 15 were employed. The effects of the degree of stenosis, the stenosis geometry and the presence or absence of the downstream confining wall on the development of flow disturbances were investigated. Visualization of the distal flow patterns in stenotic and free jets illustrated the existence of complex fields which included vortex shedding, highly turbulent regions, and recirculation zones. Significant flow disorder was created by a mild stenosis in pulsatile, but not in steady, flow. Nondimensionalization employing the stenosis diameter and flow velocity in the throat of the constriction correlates the vortex shedding frequency and energy spectra within a limited postestenotic region.     

Numerical simulation of two-phase non-Newtonian blood flow with fluid-structure interaction in aortic dissection

The behavior of blood cells and vessel compliance significantly influence hemodynamic parameters, which are closely related to the development of aortic dissection. Here the two-phase non-Newtonian model and the fluid-structure interaction (FSI) method are coupled to simulate blood flow in a patient-specific dissected aorta. Moreover, three-element Windkessel model is applied to reproduce physiological pressure waves. Important hemodynamic indicators, such as the spatial distribution of red blood cells (RBCs) and vessel wall displacement, which greatly influence the hemodynamic characteristics are analyzed. Results show that the proximal false lumen near the entry tear appears to be a vortex zone with a relatively lower volume fraction of RBCs, a low time-averaged wall shear stress (TAWSS) and a high oscillatory shear index (OSI), providing a suitable physical environment for the formation of atherosclerosis. The highest TAWSS is located in the narrow area of the distal true lumen which might cause further dilation. TAWSS distributions in the FSI model and the rigid wall model show similar trend, while there is a significant difference for the OSI distributions. We suggest that an integrated model is essential to simulate blood flow in a more realistic physiological environment with the ultimate aim of guiding clinical treatment.     

Computer simulation of blood flow in the human arm

This paper considers a finite element method to characterize blood flow in the human arm arteries. A set of different pressure waveforms, which represent normal and diseased heart pulses, is used for the proximal boundary conditions, and a modified Windkessel model is used for the distal arterial boundary conditions. A comparison of the distal pressure and flow waveforms, for each different proximal pressure, is made to determine whether such waveforms are significantly altered from normal waveforms. The results show that the distal pressure and/or flow waveforms in certain cases are sufficiently different to be possibly used as a diagnostic indicator of an abnormal heart condition. Also considered is the effect of stenosis, change of compliance, and dilatation of the distal beds on the pressure and flow waveforms. A stenosis which has an area reduction of greater than approximately 75% is found to significantly alter both the distal pressure and flow waveforms. Changes in arterial compliance, however, do not strongly influence the waveforms. Dilatation of distal vascular beds is simulated by reducing the lumped resistance of these beds, and this reduction increases mean flow and decreases mean distal pressure, but has little effect on the basic shape of either the pressure or flow waveform.     

Finite element simulation of pulsatile flow through arterial stenosis.

The problem of blood flow through a stenosis is solved using the incompressible Navier-Stokes equations in a rigid circular tube presenting a partial occlusion. Calculations are based on a Galerkin finite element method. The time marching scheme employs a predictor-corrector technique using a variable time step. Results are obtained for steady and physiological pulsatile flows. Computational experiments analyse the effect of varying the degree of stenosis, the stricture length, the Reynolds number and Womersley number. The method gives results which agree well with previous computations for steady flows and experimental findings for steady and pulsatile flows.     

Steady flow of a viscous fluid through a network of tubes with applications to the human arterial system

The paper presents a finite-element model for the analysis of steady flow of a viscous fluid through a connected system of elastic tubes with the aim of simulating the conditions of blood flow through the human arterial system. The governing equations of the model are non-linear in character and are solved through an iterative computational procedure. This model is capable of incorporating the effects of stenosis on flow and pressure. Typical results are presented and discussed. Quantitative results have been obtained on blood flow through a model of the human arterial system corresponding to the sets of prescribed conditions at the terminations. Also computational results on the effect of stenosis in typical arteries of the system are presented.     

Unsteady and three-dimensional simulation of blood flow in the human aortic arch

A three-dimensional and pulsatile blood flow in a human aortic arch and its three major branches has been studied numerically for a peak Reynolds number of 2500 and a frequency (or Womersley) parameter of 10. The simulation geometry was derived from the three-dimensional reconstruction of a series of two-dimensional slices obtained in vivo using CAT scan imaging on a human aorta. The numerical simulations were obtained using a projection method, and a finite-volume formulation of the Navier-Stokes equations was used on a system of overset grids. Our results demonstrate that the primary flow velocity is skewed towards the inner aortic wall in the ascending aorta, but this skewness shifts to the outer wall in the descending thoracic aorta. Within the arch branches, the flow velocities were skewed to the distal walls with flow reversal along the proximal walls. Extensive secondary flow motion was observed in the aorta, and the structure of these secondary flows was influenced considerably by the presence of the branches. Within the aorta, wall shear stresses were highly dynamic, but were generally high along the outer wall in the vicinity of the branches and low along the inner wall, particularly in the descending thoracic aorta. Within the branches, the shear stresses were considerably higher along the distal walls than along the proximal walls. Wall pressure was low along the inner aortic wall and high around the branches and along the outer wall in the ascending thoracic aorta. Comparison of our numerical results with the localization of early atherosclerotic lesions broadly suggests preferential development of these lesions in regions of extrema (either maxima or minima) in wall shear stress and pressure.     

Application of large-eddy simulation to the study of pulsatile flow in a modeled arterial stenosis

The technique of large-eddy simulation (LES) has been applied to the study of pulsatile flow through a modeled arterial stenosis. A simple stenosis model has been used that consists of a one-sided 50 percent semicircular constriction in a planar channel. The inlet volume flux is varied sinusoidally in time in a manner similar to the laminar flow simulations of Tutty (1992). LES is used to compute flow at a peak Reynolds number of 2000 and a Strouhal number of 0.024. At this Reynolds number, the flow downstream of the stenosis transitions to turbulence and exhibits all the classic features of post-stenotic flow as described by Khalifa and Giddens (1981) and Lieber and Giddens (1990). These include the periodic shedding of shear layer vortices and transition to turbulence downstream of the stenosis. Computed frequency spectra indicate that the vortex shedding occurs at a distinct high frequency, and the potential implication of this for noninvasive diagnosis of arterial stenoses is discussed. A variety of statistics have been also extracted and a number of other physical features of the flow are described in order to demonstrate the usefulness of LES for the study of post-stenotic flows.     

Influence of arterial wall compliance on the pressure drop across coronary artery stenoses under hyperemic flow condition

Hemodynamic endpoints such as flow and pressure drop are often measured during angioplasty procedures to determine the functional severity of a coronary artery stenosis. There is a lack of knowledge regarding the influence of compliance of the arterial wall-stenosis on the pressure drop under hyperemic flows across coronary lesions. This study evaluates the influence in flow and pressure drop caused by variation in arterial-stenosis compliance for a wide range of stenosis severities. The flow and pressure drop were evaluated for three different severities of stenosis and tested for limiting scenarios of compliant models. The Mooney-Rivlin model defined the non-linear material properties of the arterial wall and the plaque regions. The non-Newtonian Carreau model was used to model the blood flow viscosity. The fluid (blood)-structure (arterial wall) interaction equations were solved numerically using the finite element method. Irrespective of the stenosis severity, the compliant models produced a lower pressure drop than the rigid artery due to compliance of the plaque region. A wide variation in the pressure drop was observed between different compliant models for significant (90% area occlusion) stenosis with 41.0, 32.1, and 29.8 mmHg for the rigid artery, compliant artery with calcified plaque, and compliant artery with smooth muscle cell proliferation, respectively. When compared with the rigid artery for significant stenosis the pressure drop decreased by 27.7% and 37.6% for the calcified plaque and for the smooth muscle cell proliferation case, respectively. These significant variations in pressure drop for the higher stenosis may lead to misinterpretation and misdiagnosis of the stenosis severity.     

Combined MR imaging and numerical simulation of flow in realistic arterial bypass graft models

We report methods for (a) transforming a three-dimensional geometry acquired by magnetic resonance angiography (MRA) in vivo, or by imaging a model cast, into a computational surface representation, (b) use of this to construct a three dimensional numerical grid for computational fluid dynamic (CFD) studies, and (c) use of the surface representation to produce a stereo-lithographic replica of the real detailed geometry, at a scale convenient for detailed magnetic resonance imaging (MRI) flow studies. This is applied to assess the local flow field in realistic geometry arterial bypass grafts. Results from a parallel numerical simulation and MRI measurement of flow in an aorto-coronary bypass graft with various inlet flow conditions demonstrate the strong influence of the graft inlet waveform on the perianastomotic flow field. A sinusoidal and a multi harmonic coronary flow waveform both with a mean Reynolds number (Re) of 100 and a Womersley parameter of 2.7 were applied at the graft inlet. A weak axial flow separation region just distal to the toe was found in sinusoidal flow near end deceleration (Re = 25). At the same location and approximately the same point in the cycle (Re = 30) but in coronary flow, the axial flow separation was stronger and more spatially pronounced. No axial flow separation occurred in steady flow for Re = 100. Numerical predictions indicate a region in the vicinity of the suture line (where there is a local narrowing of the graft) with a wall shear magnitude in excess of five times that associated with fully developed flow at the graft inlet.     

Multiphysics simulation of blood flow and LDL transport in a porohyperelastic arterial wall model

Atherosclerosis localizes at a bend andor bifurcation of an artery, and low density lipoproteins (LDL) accumulate in the intima. Hemodynamic factors are known to affect this localization and LDL accumulation, but the details of the process remain unknown. It is thought that the LDL concentration will be affected by the filtration flow, and that the velocity of this flow will be affected by deformation of the arterial wall. Thus, a coupled model of a blood flow and a deformable arterial wall with filtration flow would be invaluable for simulation of the flow field and concentration field in sequence. However, this type of highly coupled interaction analysis has not yet been attempted. Therefore, we performed a coupled analysis of an artery with multiple bends in sequence. First, based on the theory of porous media, we modeled a deformable arterial wall using a porohyperelastic model (PHEM) that was able to express both the filtration flow and the viscoelastic behavior of the living tissue, and simulated a blood flow field in the arterial lumen, a filtration flow field and a displacement field in the arterial wall using a fluid-structure interaction (FSI) program code by the finite element method (FEM). Next, based on the obtained results, we further simulated LDL transport using a mass transfer analysis code by the FEM. We analyzed the PHEM in comparison with a rigid model. For the blood flow, stagnation was observed downward of the bends. The direction of the filtration flow was only from the lumen to the wall for the rigid model, while filtration flows from both the wall to the lumen and the lumen to the wall were observed for the PHEM. The LDL concentration was high at the lumenwall interface for both the PHEM and rigid model, and reached its maximum value at the stagnation area. For the PHEM, the maximum LDL concentration in the wall in the radial direction was observed at the position of 3% wall thickness from the lumenwall interface, while for the rigid model, it was observed just at the lumenwall interface. In addition, the peak LDL accumulation area of the PHEM moved about according to the pulsatile flow. These results demonstrate that the blood flow, arterial wall deformation, and filtration flow all affect the LDL concentration, and that LDL accumulation is due to stagnation and the presence of filtration flow. Thus, FSI analysis is indispensable.