A nonlinear axisymmetric model with fluid-wall interactions for steady viscous flow in stenotic elastic tubes.

Arteries with high-grade stenoses may compress under physiologic conditions due to negative transmural pressure caused by high-velocity flow passing through the stenoses. To quantify the compressive conditions near the stenosis, a nonlinear axisymmetric model with fluid-wall interactions is introduced to simulate the viscous flow in a compliant stenotic tube. The nonlinear elastic properties of the tube (tube law) are measured experimentally and used in the model. The model is solved using ADINA (Automatic Dynamic Incremental Nonlinear Analysis), which is a finite element package capable of solving problems with fluid-structure interactions. Our results indicate that severe stenoses cause critical flow conditions such as negative pressure and high and low shear stresses, which may be related to artery compression, plaque cap rupture, platelet activation, and thrombus formation. The pressure filed near a stenosis has a complex pattern not seen in one-dimensional models. Negative transmural pressure as low as -24 mmHg for a 78 percent stenosis by diameter is observed at the throat of the stenosis for a downstream pressure of 30 mmHg. Maximum shear stress as a high as 1860 dyn/cm2 occurs at the throat of the stenoses, while low shear stress with reversed direction is observed right distal to the stenosis. Compressive stresses are observed inside the tube wall. The maximal principal stress and hoop stress in the 78 percent stenosis are 80 percent higher than that from the 50 percent stenosis used in our simulation. Flow rates under different pressure drop conditions are calculated and compared with experimental measurements and reasonable agreement is found for the prebuckling stage.     

Numerical modeling of stress in stenotic arteries with microcalcifications: a parameter sensitivity study

As a follow-up to the work presented in Wenk et al. (2010, "Numerical Modeling of Stress in Stenotic Arteries With Microcalcifications: A Micromechanical Approximation," ASME J. Biomech. Eng., 132, p. 091011), a formal sensitivity study was conducted in which several model parameters were varied. The previous work only simulated a few combinations of the parameters. In the present study, the fibrous cap thickness, longitudinal position of the region of microcalcifications, and volume fraction of microcalcifications were varied over a broader range of values. The goal of the present work is to investigate the effects of localized regions of microcalcifications on the stress field of atherosclerotic plaque caps in a section of carotid artery. More specifically, the variations in the magnitude and location of the maximum circumferential stress were assessed for a range of parameters using a global sensitivity analysis method known as Sobol' indices. The stress was calculated by performing finite element simulations of three-dimensional fluid-structure interaction models, while the sensitivity indices were computed using a Monte Carlo scheme. The results indicate that cap thickness plays a significant role in the variation in the magnitude of the maximum circumferential stress, with the sensitivity to volume fraction increasing when the region of microcalcification is located at the shoulder. However, the volume fraction played a larger role in the variation in the location of the maximum circumferential stress. This matches the finding of the previous study (Wenk et al., 2010, "Numerical Modeling of Stress in Stenotic Arteries With Microcalcifications: A Micromechanical Approximation," ASME J. Biomech. Eng., 132, p. 091011), which indicates that the maximum circumferential stress always shifts to the region of microcalcification.     

Pulsatile non-Newtonian flow characteristics in a three-dimensional human carotid bifurcation model.

Numerical analysis of flow phenomena and wall shear stresses in the human carotid artery bifurcation has been carried out using a three-dimensional geometrical model. The primary aim of this study is the detailed discussion of non-Newtonian flow velocity and wall shear stress during the pulse cycle. A comparison of non-Newtonian and Newtonian results is also presented. The applied non-Newtonian behavior of blood is based on measured dynamic viscosity. In the foreground of discussion are the flow characteristics in the carotid sinus. The investigation shows complex flow patterns especially in the carotid sinus where flow separation occurs at the outer wall throughout the systolic deceleration phase. The changing sign of the velocity near the outer sinus wall results in oscillating shear stress during the pulse cycle. At the outer wall of the sinus at maximum diameter level the shear stress ranges from -1.92 N/m2 to 1.22 N/m2 with a time-averaged value of 0.04 N/m2. At the inner wall of the sinus at maximum diameter level the shear stress range is from 1.16 N/m2 to 4.18 N/m2 with a mean of 1.97 N/m2. The comparison of non-Newtonian and Newtonian results indicates unchanged flow phenomena and rather minor differences in the basic flow characteristics.     

Effects of wall shear stress and fluid recirculation on the localization of circulating monocytes in a three-dimensional flow model

There is a correlation between the location of early atherosclerotic lesions and the hemodynamic characteristics at those sites. Circulating monocytes are key cells in the pathogenesis of atherosclerotic plaques and localize at sites of atherogenesis. The hypothesis that the distribution of monocyte adhesion to the vascular wall is determined in part by hemodynamic factors was addressed by studying monocyte adhesion in an in vitro flow model in the absence of any biological activity in the model wall.

Suspensions of U937 cells were perfused (Re = 200) through an axisymmetric silicone flow model with a stenosis followed by a reverse step. The model provided spatially varying wall shear stress, flow separation and reattachment, and a three-dimensional flow pattern. The cell rolling velocity and adhesion rates were determined by analysis of videomicrographs. Wall shear stress was obtained by numerical solution of the equations of fluid motion. Cell adhesion patterns were also studied in the presence of chemotactic peptide gradients.

The cell rolling velocity varied linearly with wall shear stress. The adhesion rate tended to decrease with increasing local wall shear stress, but was also affected by the radial component of velocity and the dynamics of the recirculation region and flow reattachment. Adhesion was increased in the vicinity of chemotactic peptide sources downstream of the expansion site. Results with human monocytes were qualitatively similar to the U937 experiments.

Differences in the adhesion rates of U937 cells occurring solely as a function of the fluid dynamic properties of the flow field were clearly demonstrated in the absence of any biological activity in the model wall.     

Effect of a lipid pool on stress/strain distributions in stenotic arteries: 3-D fluid-structure interactions (FSI) models

Nonlinear 3-D models with fluid-structure interactions (FSI) based on in vitro experiments are introduced and solved by ADINA to perform flow and stress/strain analysis for stenotic arteries with lipid cores. Navier-Stokes equations are used as the governing equations for the fluid. Hyperelastic Mooney-Rivlin models are used for both the arteries and lipid cores. Our results indicate that critical plaque stress/strain conditions are affected considerably by stenosis severity, eccentricity, lipid pool size, shape and position, plaque cap thickness, axial stretch, pressure, and fluid-structure interactions, and may be used for possible plaque rupture predictions.     

Nasal vestibule wall elasticity: interactions with a nasal dilator strip.

We studied the effect of an adhesive external nasal dilator strip (ENDS) on external nasal geometry in 20 healthy Caucasian adults (10 men, 10 women; age 21-45 yr). The recoil force exerted by ENDS was estimated by bending the device (n = 10) with known weights. In the horizontal direction, a small/medium-sized ENDS in situ exerted a unilateral recoil force of 21.4-22.6 g. Application of ENDS resulted in a displacement of the lateral nasal vestibule walls that had both anterosuperior and horizontal components and that was maintained over an 8-h period. The resultant unilateral nasal vestibule wall displacement at the tip of the device was at 47.6 +/- 2.0 degrees to the horizontal (as related to the plane of the device when in situ) and had a magnitude of 3.5 +/- 0.1 mm. ENDS increased external nasal cross-sectional area by 23.0-65.3 mm2. Nasal vestibule wall compliance was estimated at 0.05-0.16 mm/g. Thus ENDS applies a relatively constant abducting force irrespective of nasal width. Variable responsiveness to ENDS may be related to differences in elastic properties of the nasal vestibule wall.     

Steady flow velocity measurements in a pulmonary artery model with varying degrees of pulmonic stenosis

Velocity and flow visualization studies were conducted in an adult size pulmonary artery model with varying degrees of valvular stenosis, using a two dimensional laser Doppler anemometer system. Velocity measurements in the main, left and right branches of the pulmonary artery revealed that as the degree of pulmonic stenosis increased, the jet type flow created by the valve hit the distal wall of the LPA farther downstream from the junction of the bifurcation. This in turn led to higher levels of turbulent and disturbed flow, and larger secondary flow motion in the LPA compared to the RPA. The high levels of turbulence measured in the main and left pulmonary arteries with the stenotic valves, could lead to the clinically observed phenomenon of post stenotic dilatation in the MPA extending into the LPA.     

Local hemodynamics affect monocytic cell adhesion to a three-dimensional flow model coated with E-selectin.

Monocyte adhesion to the endothelium depends on concentrations of receptors/ligands, local concentrations of chemoattractants, monocyte transport to the endothelial surface and hemodynamic forces. Monocyte adhesion to the inert surface of a three-dimensional perfusion model was shown to correlate inversely with wall shear stress, but was also affected by flow patterns which influenced the near-wall cell availability. We hypothesized that (a) under the same flow conditions, insolubilized E-selectin on the model's surface may mediate adhesive interactions at higher wall shear stresses, compared to an uncoated model, and (b) pulsatile flow may modify the adhesion profile obtained under steady flow. An axisymmetric flow model with a stenosis and a sudden expansion produced a range of wall shear stresses and a separated flow with recirculation and reattachment. Pre-activated U937 cells were perfused through the model under either steady (Re = 100, 140) or pulsatile (Remean = 107) flow. The velocity field was characterized through computational fluid dynamics and validated by inert particle tracking. Surface E-selectin greatly increased cell adhesion in all regions at Re = 100 and 140, compared to an uncoated model under the same flow conditions. In regions where the cells near the wall were abundant (taper and stenosis), adhesion to E-selectin correlated with the reciprocal of local wall shear stress when flow was steady. Pulsatile flow distributed the adherent cells more evenly throughout the coated model. Hence, characterizing both the local hemodynamics and the biological activity on the vessel wall is important in leukocyte adhesion.     

Steady flow through collapsible tubes: measurements of flow and geometry.

Compliant tubes attain a complex three-dimensional geometry when the external pressure exceeds the internal pressure and the tube is partially collapsed. A new technique for remote measurement of dynamic surfaces was applied to classical experiments with collapsible tubes. This work presents measurements of the three-dimensional structure of the tube as well as pressure and flow measurements during static loading and during steady-state fluid flow. Results are shown for two tubes of the same material and internal diameter but with different wall thicknesses. The measured tube laws compare well with previously published data and suggest the possible existence of a similarity tube law. The steady flow measurements did not compare well with the one-dimensional theoretical predictions.     

A model of microbial growth in a plug flow reactor with wall attachment

A mathematical model of microbial growth for limiting nutrient in a plug flow reactor which accounts for the colonization of the reactor wall surface by the microbes is formulated and studied analytically and numerically. It can be viewed as a model of the large intestine or of the fouling of a commercial bio-reactor or pipe flow. Two steady state regimes are identified, namely, the complete washout of the microbes from the reactor and the successful colonization of both the wall and bulk fluid by the microbes. Only one steady state is stable for any particular set of parameter values. Sharp and explicit conditions are given for the stability of each, and for the long term persistence of the bacteria in the reactor.     

Steady flow through a double converging-diverging tube model for mild coronary stenoses

Velocity profiles and the pressure drop across two mild (62 percent) coronary stenoses in series have been investigated numerically and experimentally in a perspex-tube model. The mean flow rate was varied to correspond to a Reynolds number range of 50-400. The pressure drop across two identical (62 percent) stenoses show that for low Reynolds numbers the total effect of two stenoses equals that of two single stenoses. A reduction of 10 percent is found for the higher Reynolds numbers investigated. Numerical and experimental results obtained for the velocity profiles agree very well. The effect of varying the converging angle of a single mild (62 percent) coronary stenosis on the fluid flow has been determined numerically using a finite element method. Pressure-flow relation, especially with respect to relative short stenoses, is discussed.     

The influence of the non-Newtonian properties of blood on the flow in large arteries: steady flow in a carotid bifurcation model.

Laser Doppler anemometry experiments and finite element simulations of steady flow in a three dimensional model of the carotid bifurcation were performed to investigate the influence of non-Newtonian properties of blood on the velocity distribution. The axial velocity distribution was measured for two fluids: a non-Newtonian blood analog fluid and a Newtonian reference fluid. Striking differences between the measured flow fields were found. The axial velocity field of the non-Newtonian fluid was flattened, had lower velocity gradients at the divider wall, and higher velocity gradients at the non-divider wall. The flow separation, as found with the Newtonian fluid, was absent. In the computations, the shear thinning behavior of the analog blood fluid was incorporated through the Carreau-Yasuda model. The viscoelastic properties of the fluid were not included. A comparison between the experimental and numerical results showed good agreement, both for the Newtonian and the non-Newtonian fluid. Since only shear thinning was included, this seems to be the dominant non-Newtonian property of the blood analog fluid under steady flow conditions.     

Numerical investigation of the three-dimensional flow in a human lung model

The flow field at inspiration and expiration in the upper human airways consisting of the trachea down to the sixth generation of the bronchial tree is numerically simulated. The three-dimensional steady flow at a hydraulic diameter-based Reynolds number Re(D)=1250 is computed via a lattice-Boltzmann method (LBM). The simulation is validated by the experimental data based on particle-image velocimetry (PIV) measurements. The good agreement between numerical and experimental results is evidenced by comparing velocity contours and distributions in a defined reference plane. The results show the LBM to be an accurate tool to numerically predict flow structures in the human lung. Using an automatic Cartesian grid generator, the overall process time from meshing to a steady-state solution is <12h. Moreover, the numerical simulation allows a closer analysis of the secondary flow structures than in the experimental investigation. The three-dimensional streamline patterns reveal some insight on the air exchange mechanism at inspiration and expiration. At inspiration, the slower near-wall tracheal flow enters through the right principal bronchus into the right upper lobar bronchus. The bulk mass flux in the trachea is nearly evenly distributed over the left upper, center and lower lobar bronchi and the right center and lower bronchi. At expiration, the air from the right upper lobar bronchus enters the right center of the trachea and displaces the airflow from the lower and center right bronchi such that the tracheal positions of the streamlines at inspiration and expiration are switched. The flow in the left bronchi does not show this kind of switching. The findings emphasize the impact of the asymmetry of the lung geometry on the respiratory air exchange mechanism.     

Influence of model boundary conditions on blood flow patterns in a patient specific stenotic right coronary artery

Background

In literature, the effect of the inflow boundary condition was investigated by examining the impact of the waveform and the shape of the spatial profile of the inlet velocity on the cardiac hemodynamics. However, not much work has been reported on comparing the effect of the different combinations of the inlet/outlet boundary conditions on the quantification of the pressure field and flow distribution patterns in stenotic right coronary arteries.

Method

Non-Newtonian models were used to simulate blood flow in a patient-specific stenotic right coronary artery and investigate the influence of different boundary conditions on the phasic variation and the spatial distribution patterns of blood flow. The 3D geometry of a diseased artery segment was reconstructed from a series of IVUS slices. Five different combinations of the inlet and the outlet boundary conditions were tested and compared.

Results

The temporal distribution patterns and the magnitudes of the velocity, the wall shear stress (WSS), the pressure, the pressure drop (PD), and the spatial gradient of wall pressure (WPG) were different when boundary conditions were imposed using different pressure/velocity combinations at inlet/outlet. The maximum velocity magnitude in a cardiac cycle at the center of the inlet from models with imposed inlet pressure conditions was about 29% lower than that from models using fully developed inlet velocity data. Due to the fact that models with imposed pressure conditions led to blunt velocity profile, the maximum wall shear stress at inlet in a cardiac cycle from models with imposed inlet pressure conditions was about 29% higher than that from models with imposed inlet velocity boundary conditions. When the inlet boundary was imposed by a velocity waveform, the models with different outlet boundary conditions resulted in different temporal distribution patterns and magnitudes of the phasic variation of pressure. On the other hand, the type of different boundary conditions imposed at the inlet and the outlet did not have significant effect on the spatial distribution patterns of the PD, the WPG and the WSS on the lumen surface, regarding the locations of the maximum and the minimum of each quantity.

Conclusions

The observations from this study indicated that the ways how pressure and velocity boundary conditions are imposed in computational models have considerable impact on flow velocity and shear stress predictions. Accuracy of in vivo measurements of blood pressure and velocity is of great importance for reliable model predictions.

     

Steady flow development past valve prostheses in a model human aorta—I. Centrally occluding valves

In this paper, laser-Doppler anemometry measurement of steady flow development in a model human aorta has been reported. Studies were made with uniform entry flow at the root of the aorta and our measurements showed the establishment of a pair of Dean vortices in the mid-arch region. Subsequently, the nature of flow development past centrally occluding caged ball valves in the model aorta was investigated. Our studies showed that in the ascending aorta, an asymmetric velocity profile is obtained with larger velocity gradients towards the inner wall of tertiary curvature (anatomically the left lateral wall) with centrally occluding valves. The peripheral flow past these valves prevented the development of Dean vortices in the mid-arch region. The caged ball valves at the root of the aorta had no discernible effect on the velocity profiles in the brachio-cephalic artery.     

Steady flow development past valve prostheses in a model human aorta—II. Tilting disc valves

The flow development in the model human aorta with uniform entry as well as with centrally occuluding valves mounted at the root of the aorta was described in Part I of this two-paper sequence. Part II deals with the flow development in the model aorta with tilting disc valves mounted at the root of the aorta. Bjork-Shiley and Hall-Kaster tilting disc valves were mounted in three different orientations with respect to the root of the aorta. The velocity profiles and turbulent stresses were measured with laser-Doppler anemometry. Our results under steady flow conditions in the model human aorta show quantitatively that the flow development in the ascending aorta as well as in the brachio-cephalic artery are strongly dependent on the orientation of the tilting disc valves. With the valves tilting towards the outer wall of curvature, our results suggest a tendency for flow separation at the flow divider region of the brachio-cephalic artery.