appears at 0.15 s in systolic phase and outlet peak pressure (the blue dash line) is shown at 0.4 s in diastolic phase.

100 at each four-time point (0 s right before acceleration, 0.15 s for inlet peak velocity in systole, 0.4 s for outlet peak pressure in mid-cycle, and 0.7 s in diastole) during the cardiac cycle.

2.3. Numerical Methods

The fully coupled solid and fluidmodels were solved interactively with 0.01 s of the time step and 1 s of the total time (i.e., one cardiac cycle) based on a heart rate of 60 beats per minute using a commercial program ANSYS for structural mechanics analysis and ANSYS-CFX for fluid mechanics analysis. The fluid pressure at the fluid- solid interface can be applied as a load for the structural analysis, and the resulting displacement, velocity, or acceleration obtained in the structural analysis can be passed back for fluid analysis. As a solver scheme, a high resolution for advection and second order backward Euler for the transientscheme were utilized, and the convergence criterion was set as 10−4. Iteration between the structural analysis solution and fluid mechanics solution continues until overall equilibrium is reached.

3. Results and Discussion

3.1. Comparison of Velocity Profiles of Newtonian and Non-Newtonian Models

In reality, the blood has been known to be a non-Newtonian fluid where the blood viscosity changes with the hematocrit and macroglobulin concentration [18] , but many studies have assumed the blood as a Newtonian fluid under the high shear rate (i.e., shear rate > 100 s−1) [17] . As mentioned previously, I assumed the blood as a Newtonian fluid for the FSI analysis. In order to justify this assumption, it should be guaranteed that physical quantities such as velocity profiles, shear rates, or wall shear stress between the Newtonian and non-Newtonian models should be similar. For this, I employed various blood viscosity models for examining their effect on the hemodynamics. Table 2 shows several blood models for the effective viscosity. The non-Newtonian models include power law, Carreau, Casson, and Walburn-Schneck. As a Newtonian model, plasma and blood with the constant viscosity were utilized. For the comparison between viscosity models, the rigid RCA wall was used.

Figure 3(a) shows five major areas of the stenotic RCA to extract velocity profiles: in the upstream (L1), in the pre-stenotic region (L2), at the throat of the stenosis (L3), in the post-stenotic regions (L4), and in the downstream (L5). The diameter of each location is 3.47 mm (L1), 3.16 mm (L2), 1.43 mm (L3), 2.61 mm (L4), and 2.50 mm (L5), respectively. Figure 3(b) shows velocity profiles at three time points (0.15 s, 0.4 s, and 0.7 s) for each location. Overall, the velocity magnitude at the throat of the stenosis (L3) is highest for all the time points, regardless of blood models. In addition, at the peak flow in systolic phase (0.15 s), the maximum velocity values (around 4 m/s) are approximately two times as high as those (around 2 m/s) of other time points (0.4 s and 0.7 s). Noticeably, at the post-stenotic region (L4) and the downstream (L5), I found that the pattern of the velocity profiles were different from other locations (L1, L2, and L3). The velocity profiles are highly skewed due to the curved shape of the RCA and the resulting higher inertial force of the fluid.

Table 2. The Newtonian and non-Newtonian blood models for the effective viscosity.

Figure 3. Comparison of Newtonian and non-Newtonian models: (a) Five different locations to extract velocity profiles: in the upstream (L1), in the pre-stenotic region (L2), at the throat of the stenosis (L3), in the post-stenotic regions (L4), and in the downstream (L5). (b) Velocity profiles with three time points for five locations.

For the comparison of the viscosity models, the Casson model predictsslightly higher velocity magnitudes than others, especially for L1, L2, and L4, whereas the plasma model predicts relatively lower values at the same locations. Among non-Newtonian models, the Walburn-Schneck model tends to underpredict the velocity magnitudes, particularly for L1, L2, and L4 compared to others. Overall, only subtle difference in velocity mag- nitudes among the viscosity models―except for the plasma model and Walburn-Schneck―was observed for all locations. It should be noted that the Newtonian blood model in this RCA model presents relatively similar velocity patterns and values for all times and locations to other non-Newtonian models. For example, for the comparison of the blood Newtonian model to non-Newtonian models except for the Walburn-Schneck model, the maximum difference is found to be within 4%. In light of the result, it is believed that applying the Newtonian blood model to FSI analysis (i.e., flexible wall) does not have a significant impact on the hemodynamic results, compared to non-Newtonian models.

3.2. Total Displacement and von Mises Stress of Flexible Wall

Figure 4 shows a total displacement of the flexible RCA wall. The extracted total displacement is a quantity of the movement at the wall interface at each time relative to the initial state. Overall, the total displacement was

Figure 4. Total displacement of the flexible RCA wall. The center of the stenotic area at the initial state is expressed as a black cross.

found to be highest at the stenotic throat. The maximum displacements for 0.15 s, 0.4 s, and 0.7 s were shown to be about 1.67 mm, 2.14 mm, and 1.46 mm, respectively. These values relatively agree with data reported in the other study by showing the same order of the magnitude [17] .

The maximal von Mises stress was found to be 92.06 kPa at the downstream away from the throat. This is less than 10% of the breaking strength of the RCA wall, which has been known to be approximately 1 MPa [24] . However, the von Mises stress at the throat yields a very low value.

3.3. Comparison of Pressure and Velocity Profiles between the Rigid and Flexible RCA Wall

Figure 5 shows the comparison of temporal averaged velocity (a) and pressure profiles (b) at the cross-section of the stenotic throat (L3). It is presented that velocity and pressure patterns between the two models are similar but there exists quantitative difference. First, the averaged velocity magnitudes for the flexible wall are smaller than the rigid wall for most of the time points. In particular, the degree of the difference in velocity is found to be largest (15%) at the inlet peak flow (0.15 s). This is because the high fluid pressure force (due to the peak flow in the flexible RCA) is transferred to the wall and subsequently causes wall movement and deformation. Correspondingly, there are more spaces where fluid can move and thus, the overall velocity magnitude might be lower. On the other hand, there is very small difference (less than 5%) in velocity at the 0.43 s. This is thought to be because the low inlet velocity conditions lead to the small fluid pressure exerting on the wall, which does not much affect the wall movement or deformation. Consequently, this might lead to no significant difference between the two models at this time point. In contrast to the case of the velocity, the pressure values for the flexible wall are slightly larger, but the difference in pressure between the two models is on average smaller.

Figure 5. Comparison of averaged velocity and pressure profiles between the rigid and flexible RCA wall at the cross- section of the stenotic throat (L3): (a) velocity and (b) pressure profiles, respectively.

3.4. Comparison of the Particle Behavior between Rigid and Flexible RCA Wall

Figure 6 shows the comparison of the particle behavior within the rigid and flexible RCAs using four different injecting time points (tin = 0 s, 0.15 s, 0.4 s and 0.7 s). The red color indicates trajectories of deposited particles and the green indicates trajectories of particles exiting from the RCA. This qualitative result clearly shows that the number of particles exiting (green) is higher than deposited particles (red) within the RCA.

Such a trend was shown in quantitative results (Figure 7) as well. Figure 7(a) shows one of the sample images of the deposited particles in stenotic (yellow oval) and non-stenotic area. Figure 7(b) shows quantitative results of particle deposition in the RCAs. As in the case of the qualitative results in Figure 6, it was clearly observed that the overall deposited particles were less than 50% of all injected particles, regardless of rigid and flexible wall. In other words, more particles tend to escape the RCAs than deposit. The ratio of total particle deposition (blue color) in the rigid RCA wall is 43%, 37%, 25%, and 23% for 0 s, 0.15 s, 0.4 s, and 0.7 s, respectively. In the case of the flexible RCA wall, the total particle deposition (red color) of 41%, 29%, 24%, and 25% for the same time points above was shown. The result clearly shows that the number of particles deposited in the rigid wall is generally higher than the flexible wall, except for 0.7 s. In particular, at 0.15 s (inlet peak flow), there was the most notable difference in particle deposition between the rigid and flexible wall showing 8% discrepancy. On the other hand, the difference was not more than 2% for other cases. When considering the number of particles deposited only in stenotic areas, the difference between rigid and flexible models became substantially lower. Such a small difference might be associated with the fact that dynamic vessel motion on flow behaviors is less significant than other factors such as arterial geometry or pulsatile flow condition [25] . Similarly, the flexible wall might not affect the particle behaviors that much either. This result implies that for the prediction of drug transport within the RCAs, rigid wall model as well as flexible wall model can be one of the options.

Interestingly, for the time period up to 0.15 s as shown in Figure 7(b), the number of particles deposited in the stenotic area was larger than the non-stenotic area for both rigid and flexible wall. On the contrary, from the time point of 0.4 s, the trend became inverse. The most notable feature is that the injecting time of 0 s was most effective for the particle deposition for both two models. During this injecting time, the particle deposition in the stenotic area is found to be highest presenting approximately 24% for both the rigid wall and flexible wall. This is presumed due to the fact that, from the initial state (0 s) to peak flow (0.15 s), fluid flow tends to be accelerated, whereas fluid decelerates after then as shown in the inlet velocity condition of Figure 2. During this time period, the maximum particle traveling time is shown to be slightly less than 0.2 s. Thus, the injected particles at 0 s might be strongly affected by higher inertia force due to high velocity magnitudes and thus deposit more, compared to other time points.

Figure 6. Comparison of particle deposition rate of the rigid and flexible RCA wall. The red color indicates trajectories of deposited particles and the green is trajectories of particles exiting from the RCA.

Figure 7. A sample image (a) and quantitative results (b) of particle deposition of the rigid and flexible wall: stenosis, non- stenosis, and total indicate particles deposited in the stenotic area, non-stenotic area, and both in the stenotic and non-stenotic area.

3.5. Perspectives

In this study, the numerical simulations of the curved stenotic RCA were quantitatively performed in terms of velocity profiles, wall displacement, maximal von Mises stress, and particle deposition in the rigid and flexible RCAs. I employed the real pulse wave for the inlet and outlet condition to make the simulation closer to the actual physiological phenomena. In addition, I justified the utilization of the Newtonian blood viscosity model by comparing velocity profiles in five major areas of the RCAs with non-Newtonian viscosity models.

Since the characteristics of blood vessel are highly complicated, modeling a 3D realistic blood vessel is essential for accurate surgical prediction. Complex flow such as turbulent flow cannot be precisely described in two- dimension (2D). In contrast, our 3D model and simulation of deformable blood vessel and drug transport enable to describe more realistic phenomena, and thereby can give exact information and feedback to the patient. In addition, it can reduce the amount of time for model validation.

Nevertheless, several limitations still exist that need to be addressed. First, only a specific particle size (3.8 µm, microparticle composed of aggregates of nanoparticles) was adopted as reported [9] . However, in reality the size of the drug varies, ranging from submicron to micron, and thus different-sized particles should be employed to investigate the size effect. In particular, submicron particles tend to moverandomly in a fluid, Brownian motion, so the related force (i.e., Brownian force) should be taken into account for the particle transport in this case. In addition, I do not consider surface property or signaling factor for the attachment of the drugs. Finally, the simulation was carried out based on the patient-specific size and shape of the stenotic RCA. Thus,varioussized and shaped RCA models of different patients should be tested to guarantee the reliability of the efficacy of the drug. These addressed factors will be the primary goals in our later research.

4. Conclusions

I have developed a computational model with FSI for blood flow and drug delivery of the curved RCA with stenosis. The key points in this research are as follows:

1) The Newtonian blood viscosity model used in this study shows similar flow patterns and velocity magnitudes to non-Newtonian models.

2) The maximal total displacement and von Mises stress of the RCA wall are 2.14 mm and 92.06 kPa, which are both on the same order of the magnitude reported in the other study.

3) There is no notable difference in the particle sedimentation to the stenotic areas between the rigid and flexible RCA models.

4) The effective time point for the drug injection is found to be between 0 s and 0.15 s (i.e., fluid acceleration region).

I believe that this information will significantly contribute to the design of a drug delivery system for the treatment of the stenotic arteries by targeting drugs selectively to the stenotic sites.


This work was done in the mechanical engineering computer laboratory (MECL) at Purdue University.

Cite this paper
Park, S. (2016) Fluid-Structure Interaction Analysis for Drug Transport in a Curved Stenotic Right Coronary Artery. Journal of Biosciences and Medicines, 4, 105-115. doi: 10.4236/jbm.2016.45011.

[1]   Bonow, R.O., et al. (2002) World Heart Day 2002—The International Burden of Cardiovascular Disease: Responding to the Emerging Global Epidemic. Circulation, 106, 1602-1605.

[2]   Valencia, A. and Villanueva, M. (2006) Unsteady Flow and Mass Transfer in Models of Stenotic Arteries Considering Fluid-Structure Interaction. International Communications in Heat and Mass Transfer, 33, 966-975.

[3]   Ambrose, J.A., et al. (1988) Angiographic Progression of Coronary Artery Disease and the Development of Myocardial Infarction. Journal of the American College of Cardiology, 12, 56-62.

[4]   Kirpalani, A., et al. (1999) Velocity and Wall Shear Stress Patterns in the Human Right Coronary Artery. Journal of Biomechanical Engineering-Transactions of the Asme, 121, 370-375.

[5]   Bathe, M. and Kamm, R.D. (1999) A Fluid-Structure Interaction Finite Element Analysis of Pulsatile Blood Flow through a Compliant Stenotic Artery. Journal of Biomechanical Engineering-Transactions of the Asme, 121, 361-369.

[6]   Tang, D.L., et al. (2002) Simulating Cyclic Artery Compression Using a 3D Unsteady Model with Fluid-Structure Interactions. Computers & Structures, 80, 1651-1665.

[7]   Giannoglou, G.D., et al. (2005) Wall Pressure Gradient in Normal Left Coronary Artery Tree. Medical Engineering & Physics, 27, 455-464.

[8]   Giannoglou, G.D., et al., (2002) Haemodynamic Factors and the Important Role of Local Low Static Pressure in Coronary Wall Thickening. International Journal of Cardiology, 86, 27-40.

[9]   Korin, N., et al. (2012) Shear-Activated Nanotherapeutics for Drug Targeting to Obstructed Blood Vessels. Science, 337, 738-742.

[10]   Vijayaratnam, P.R., et al. (2015) The Impact of Blood Rheology on Drug Transport in Stented Arteries: Steady Simulations. PLoS One, 10, e0128178.

[11]   Stangeby, D.K. and Ethier, C.R. (2002) Computational Analysis of Coupled Blood-Wall Arterial LDL Transport. Journal of Biomechanical Engineering, 124, 1-8.

[12]   Liu, G.Y., et al. (2014) Numerical Simulation of Flow in Curved Coronary Arteries with Progressive Amounts of Stenosis Using Fluid-Structure Interaction Modelling. Journal of Medical Imaging and Health Informatics, 4, 605-611.

[13]   Park, S.M., et al. (2010) In Vitro Hemodynamic Study on the Stenotic Right Coronary Artery Using Experimental and Numerical Analysis. Journal of Mechanics in Medicine and Biology, 10, 695-712.

[14]   Gradus-Pizlo, I., et al. (2003) Left Anterior Descending Coronary Artery Wall Thickness Measured by high-Frequency Transthoracic and Epicardial Echocardiography Includes Adventitia. American Journal of Cardiology, 91, 27-32.

[15]   Perry, R., et al. (2013) Coronary Artery Wall Thickness of the Left Anterior Descending Artery Using High Resolution Transthoracic Echocardiography—Normal Range of Values. Echocardiography—A Journal of Cardiovascular Ultrasound and Allied Techniques, 30, 759-764.

[16]   Koshiba, N., et al. (2007) Multiphysics Simulation of Blood Flow and LDL Transport in a Porohyperelastic Arterial Wall Model. Journal of Biomechanical Engineering-Transactions of the ASME, 129, 374-385.

[17]   Torii, R., et al. (2009) Fluid-Structure Interaction Analysis of a Patient-Specific Right Coronary Artery with Physiological Velocity and Pressure Waveforms. Communications in Numerical Methods in Engineering, 25, 565-580.

[18]   Mannik, M. (1974) Blood Viscosity in Waldenstrom’s Macroglobulinemia. Blood, 44, 87-98.

[19]   Pedersen, L., Nielsen, E.B., Christensen, M.K., Buchwald, M. and Nybo, M. (2014) Measurement of Plasma Viscosity by Free Oscillation Rheometry: Imprecision, Sample Stability and Establishment of a New Reference Range. Annals of Clinical Biochemistry, 51, 495-498.

[20]   Ballyk, P.D., Steinman, D.A. and Ethier, C.R. (1994) Simulation of Non-Newtonian Blood Flow in an End-to-Side Anastomosis. Biorheology, 31, 565-586.

[21]   Cho, Y.I. and Kensey, K.R. (1991) Effects of the Non-Newtonian Viscosity of Blood on Flows in a Diseased Arterial Vessel. Part 1: Steady Flows. Biorheology, 28, 241-262.

[22]   Razavi, A., Shirani, E. and Sadeghi, M.R. (2011) Numerical Simulation of Blood Pulsatile Flow in a Stenosed Carotid Artery Using Different Rheological Models. Journal of Biomechanics, 44, 2021-2030.

[23]   Walburn, F.J. and Schneck, D.J. (1976) A Constitutive Equation for Whole Human Blood. Biorheology, 13, 201-210.

[24]   Claes, E., et al. (2010) Mechanical Properties of Human Coronary Arteries. 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology, Buenos Aires, 31 August-4 September 2010, 3792-3795.

[25]   Myers, J.G., Moore, J.A., Ojha, M., Johnston, K.W. and Ethier, C.R. (2001) Factors Influencing Blood Flow Patterns in the Human Right Coronary Artery. Annals of Biomedical Engineering, 29, 109-120.