Research Papers

Numerical Simulation of Particle Transport and Deposition in the Pulmonary Vasculature

[+] Author and Article Information
Salman Sohrabi

Department of Mechanical
Engineering & Mechanics,
Lehigh University,
Bethlehem, PA 18015
e-mail: Sas713@Lehigh.edu

Junda Zheng

Department of Mechanical
Engineering and Mechanics,
Lehigh University,
Bethlehem, PA 18015
e-mail: Juz212@Lehigh.edu

Ender A. Finol

Department of Biomedical Engineering,
The University of Texas at
San Antonio, TX 78249
e-mail: Ender.finol@utsa.edu

Yaling Liu

Department of Mechanical
Engineering & Mechanics,
Bioengineering Program Lehigh University,
Bethlehem, PA 18015
e-mail: Yal310@lehigh.edu

1Corresponding author.

Manuscript received April 21, 2014; final manuscript received October 6, 2014; accepted manuscript posted October 15, 2014; published online November 7, 2014. Assoc. Editor: Naomi Chesler.

J Biomech Eng 136(12), 121010 (Dec 01, 2014) (11 pages) Paper No: BIO-14-1173; doi: 10.1115/1.4028800 History: Received April 21, 2014; Revised October 06, 2014; Accepted October 15, 2014

To quantify the transport and adhesion of drug particles in a complex vascular environment, computational fluid particle dynamics (CFPD) simulations of blood flow and drug particulate were conducted in three different geometries representing the human lung vasculature for steady and pulsatile flow conditions. A fully developed flow profile was assumed as the inlet velocity, and a lumped mathematical model was used for the calculation of the outlet pressure boundary condition. A receptor–ligand model was used to simulate the particle binding probability. The results indicate that bigger particles have lower deposition fraction due to less chance of successful binding. Realistic unsteady flow significantly accelerates the binding activity over a wide range of particle sizes and also improves the particle deposition fraction in bifurcation regions when comparing with steady flow condition. Furthermore, surface imperfections and geometrical complexity coupled with the pulsatility effect can enhance fluid mixing and accordingly particle binding efficiency. The particle binding density at bifurcation regions increases with generation order and drug carriers are washed away faster in steady flow. Thus, when studying drug delivery mechanism in vitro and in vivo, it is important to take into account blood flow pulsatility in realistic geometry. Moreover, tissues close to bifurcations are more susceptible to deterioration due to higher uptake.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Liu, Y., Tan, J., Thomas, A., Ou-Yang, D., and Muzykantov, V. R., 2012, “The Shape of Things to Come: Importance of Design in Nanotechnology for Drug Delivery,” Ther. Delivery, 3(2), pp. 181–194. [CrossRef]
Saul, J. M., Annapragada, A., Natarajan, J. V., and Bellamkonda, R. V., 2003, “Controlled Targeting of Liposomal Doxorubicin via the Folate Receptor in vitro,” J. Controlled Release, 92(1–2), pp. 49–67. [CrossRef]
Wang, S., Zhou, Y., Tan, J., Xu, J., Yang, J., and Liu, Y., 2014, “Computational Modeling of Magnetic Nanoparticle Targeting to Stent Surface Under High Gradient Field,” Comput. Mech., 53(3), pp. 403–412. [CrossRef] [PubMed]
Kona, S., Dong, J. F., Liu, Y., and Nguyen, K. T., 2012, “Biodegradable Nanoparticles Mimicking Platelet Binding Facilitate Targeted Drug Delivery to Endothelial Cells under Shear Flow,” Int. J. Pharm., 2(423), pp. 516–524. [CrossRef]
Tan, J., Shah, S., Thomas, A., Ou-Yang, H. D., and Liu, Y., 2013, “The Influence of Size, Shape and Vessel Geometry on Nanoparticle Distribution,” Microfluid. Nanofluid., 14(1–2), pp. 77–87. [CrossRef]
Tan, J., Wang, S., Yang, J., and Liu, Y., 2013, “Coupled Particulate and Continuum Model for Nanoparticle Targeted Delivery,” Comput. Struct., 122(6), pp. 128–134. [CrossRef] [PubMed]
Asgharian, B., Hofmann, W., and Bergmann, R., 2001, “Particle Deposition in a Multiple-Path Model of the Human Lung,” Aerosol Sci. Technol., 34(4), pp. 332–339. [CrossRef]
Lai, T., Morsi, Y., and Mazumdar, J., 2002, “Modelling and Simulation of Particle Deposition in the Human Lung,” Profiles in Industrial Research Knowledge and Innovation, pp. 313–320.
Johnston, B. M., Johnston, P. R., Corney , S., and Kilpatrick , D., 2006, “Non-Newtonian Blood Flow in Human Right Coronary Arteries: Transient Simulations,” J. Biomech., 39(6), pp. 1116–1128. [CrossRef] [PubMed]
Tan, J., Thomas, A., and Liu, Y., 2012, “Influence of Red Blood Cells on Nanoparticle Targeted Delivery in Microcirculation,” Soft Matter, 8(6), pp. 1934–1946. [CrossRef]
Liu, Y., Shah, S., and Tan, J., 2012, “Computational Modeling of Nanoparticle Targeted Drug Delivery,” Rev. Nanosci. Nanotechnol., 1(1), pp. 66–83. [CrossRef]
Nowak, N., Kakade, P. P., and Annapragada, A. V., 2003, “Computational Fluid Dynamics Simulation of Airflow and Aerosol Deposition in Human Lungs,” Ann. Biomed. Eng., 31(4), pp. 374–390. [CrossRef] [PubMed]
Lin, C.-L., and HoffmanE. A., 2004, “Comparison of Airflows in Weibel-based and CT-based Human Lung Geometries,” APS Division of Fluid Dynamics Meeting Abstracts, 1, pp. 21–23.
Lin, C-L., Tawhai, M. H., McLennan, G., and Hoffman, E. A., 2006, “Multiscale Simulation of Air Flow in the CT-Based Lung Model,” J. Biomech., 39(4), p. S265. [CrossRef]
Sanghun, C.,Hoffman, E. A., Tawhai, M., Castro, M., and Lin, C.-L., 2012, “A Numerical Study Of Airway Resistance And Particle Deposition In Normal And Asthmatic Lungs,” in A66, Modeling, Mechanics and Gas Exchange, American Thoracic Society. Abstract No. A2070.
Hasegawa, M., Burnette, N., Yin, Y., Lin, C.-L., and Hoffman, E. A., 2011, “Ct-Based Characterization Of Male And Female Airway Tree Branching Patterns: A Revisit To A Half-Century-Old Post Mortem Study With New Non-Invasive Measures,” Am. J. Respir. Crit. Care Med., 183(3), pp. A5225–A5227.
Marsden, A. L., Vignon-Clementel, I. E., Chan, F. P., Feinstein, J. A., and Taylor, C. A., 2007, “Effects of Exercise and Respiration on Hemodynamic Efficiency in CFD Simulations of the Total Cavopulmonary Connection,” Ann. Biomed. Eng., 35(2), pp. 250–263. [CrossRef] [PubMed]
Oakes, J. M., Marsden, L. A., Grandmont, C., Shadden, C. S., Darquenne, C., and Vignon-Clementel, E., 2014, “Airflow and Particle Deposition Simulations in Health and Emphysema: From In Vivo to In Silico Animal Experiments,” Ann. Biomed. Eng., 42(4), pp. 899–914. [CrossRef] [PubMed]
Schroeter, J. D., Kimbell, J. S., Asgharian, B., Tewksbury, E. W., and Singal, M., 2012, “Computational Fluid Dynamics Simulations of Submicrometer and Micrometer Particle Deposition in the Nasal Passages of a Sprague-Dawley Rat,” J. Aerosol Sci., 43(1), pp. 31–44. [CrossRef]
Chern, M.-J., Wu, M.-T., and Her, S.-W., 2012, “Numerical Study for Blood Flow in Pulmonary Arteries After Repair of Tetralogy of Fallot,” Comput. Math. Methods Med., 11(20), p. 198108. [CrossRef]
Chern, M.-J., Kanna, P. R., Lu, Y-J., Cheng, I-C., and Chang, S-C., 2010, “A CFD Study of the Interaction of Oscillatory Flows With a Pair of Side-By-Side Cylinders,” J. Fluids Struct., 26(4), pp. 626–643. [CrossRef]
Zhang, Z., C. Kleinstreuer, C., Donohue, J. F., and Kim, C. S., 2005, “Comparison of Micro- and Nano-Size Particle Depositions in a Human Upper Airway Model,” J. Aerosol Sci., 36(2), pp. 211–233. [CrossRef]
Kleinstreuer, C., and Zhang, Z., 2003, “Laminar-To-Turbulent Fluid-Particle Flows in a Human Airway Model,” Int. J. Multiphase Flow, 29(2), pp. 271–289. [CrossRef]
Worth Longest, P., and Kleinstreuer, C., 2003, “Comparison of Blood Particle Deposition Models for Non-Parallel Flow Domains,” J. Biomech., 36(3), pp. 421–430. [CrossRef] [PubMed]
Decuzzi, P., and Ferrari, M., 2006, “The Adhesive Strength of Non-Spherical Particles Mediated by Specific Interactions,” Biomaterials, 27(30), pp. 5307–5314. [CrossRef] [PubMed]
Kheyfets, V. O., O'Dell, W., Smith, T., Reilly, J. J., and Finol, E. A., 2013, “Considerations for Numerical Modeling of the Pulmonary Circulation—A Review With a Focus on Pulmonary Hypertension,” ASME J. Biomech. Eng., 135(6), p. 061011. [CrossRef]
Goldman, A., Cox, R., and Brenner, H., 1967, “Slow Viscous Motion of a Sphere Parallel to a Plane Wall—II Couette Flow,” Chem. Eng. Sci., 22(4), pp. 653–660. [CrossRef]
Haun, J. B., and Hammer, D. A., 2008, “Quantifying Nanoparticle Adhesion Mediated by Specific Molecular Interactions,” Langmuir, 24(16), pp. 8821–8832. [CrossRef] [PubMed]
Smith, L. A., Aranda-Espinoza, H., Haun, J. B., and Hammer, D. A., 2007, “Interplay Between Shear Stress and Adhesion on Neutrophil Locomotion,” Biophys. J., 92(2), pp. 632–640. [CrossRef] [PubMed]
Henk, C. B., Schlechta, B., Grampp, S., Gomischek, G., Klepetko, W., and Mostbeck, G. H., 1998, “Pulmonary and Aortic Blood Flow Measurements in Normal Subjects and Patients After Single Lung Transplantation at 0.5 T Using Velocity Encoded Cine MRI,” Chest J., 114(3), pp. 771–779. [CrossRef]
Strahler, A. N., 1952, “Hypsometric (Area-Altitude) Analysis of Erosional Topography,” Geol. Soc. Am. Bull., 63(11), pp. 1117–1142. [CrossRef]
Rider, C. A., Tang, N. J., and Fung, B. Y-C., 1993, “Morphometry of Pig Coronary Arterial Trees,” Am. J. Physiol., 265(1), pp. H350–H365. [CrossRef] [PubMed]
Huang, W., Yen, R. T., McLaurine, M., and Bledsoe, G., 1996, “Morphometry of the Human Pulmonary Vasculature,” J. Appl. Physiol., 81(5), pp. 2123–2133. [CrossRef] [PubMed]
Spilker, R. L., Feinstein, J. A., Parker, D. W., Reddy, V. M., and Taylor, C. A., 2007, “Morphometry-Based Impedance Boundary Conditions for Patient-Specific Modeling of Blood Flow in Pulmonary Arteries,” Ann. Biomed. Eng, 35(4), pp. 546–559. [CrossRef] [PubMed]
Mach, A. J., Yi,Xu., Sollier, E., Amini, H., Go, D. E., and Di Carlo, D., Mechanics of Particle Trapping and Maintenance in Micro-Scale Fluid Vortices, University of California, Los Angeles, CA.
Olufsen, M. S., and Nadim, A., 2004, “On Deriving Lumped Models for Blood Flow and Pressure in the Systemic Arteries,” Math. Biosci. Eng., 1(1), pp. 61–80. [CrossRef] [PubMed]
Savithiri, S., Pattamatta, A., and Das, S. K., 2011, “Scaling Analysis for the Anvestigation of Slip Mechanisms in Nanofluids,” Nanoscale Res. Lett., 6(1), pp. 1–15. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Artificial ideal Weibel model and (b) lung vasculature geometry reconstructed from a CT image dataset; a four-generation original (c) and oversmoothed (d) subunit of pulmonary vasculature

Grahic Jump Location
Fig. 2

Illustration of a spherical particle in contact with vascular wall. δeq is the separation distance between the spheroidal particle and the substrate; h0 is the maximum distance of the spheroid from the substrate at which a specific bond can occur; d is the diameter of the particle.

Grahic Jump Location
Fig. 3

Instantaneous inlet velocity profile estimated from experimentally measured data [30] for left and right pulmonary arteries (LPA and RPA). 0.1 m/s is approximately the average of LPA flow which occurs at 0.1 s and 0.45 s during cardiac cycle.

Grahic Jump Location
Fig. 4

Velocity contours and selected streamlines in a planar slice of original geometry in (a) steady flow case with 0.1 m/s inlet velocity and in (b), and (c) pulsatile flow case at t = 0.1 s and t = 0.45 s with 0.1 m/s instantaneous inlet velocity (refer to Fig. 6), respectively

Grahic Jump Location
Fig. 5

Demonstration of shear stress on the vessel wall in (a) steady flow case with 0.1 m/s inlet velocity and pulsatile flow case at t = 0.1 s, (b) and t = 0.45 s, and (c) with 0.1 m/s instantaneous inlet velocity

Grahic Jump Location
Fig. 6

Particle deposition fraction at two different ranges of particle sizes under various steady flow rates (a) 1 to 100 nm and (b) 1 nm to 3 μm

Grahic Jump Location
Fig. 7

Study of streamline disturbance effect on particle deposition fraction using original, oversmoothed and artificial geometry with the same inlet flow rate

Grahic Jump Location
Fig. 8

Particle deposition fraction for particles of different sizes under steady and pulsatile blood flow. It is demonstrated at two different particle diameter ranges (a) 1 to 100 nm and (b) 1 nm to 3 μm.

Grahic Jump Location
Fig. 9

Effect of vessel surface smoothness on particle deposition fraction under pulsatile flow condition

Grahic Jump Location
Fig. 10

Illustrations of deposition profile versus time for various simulation parameters. (a) Accumulative profiles of particle binding under different flow conditions and geometries (b) instantaneous deposition profile of cardiac flow in geometries with different surface smoothness.

Grahic Jump Location
Fig. 11

3D demonstration of deposition pattern when releasing two million particles with diameters ranging from 1 nm to 0.5 μm under pulsatile inlet flow in (a) original lung vascular geometry and (b) artificial Weibel geometry

Grahic Jump Location
Fig. 12

Normalized deposited number of particles with different diameter ranges deposited on various regions of the original geometry

Grahic Jump Location
Fig. 13

3D representation for trapped location of injected particles at different times (a) 1 s, (b) 2 s, and (c) 8 s after starting the simulation with cardiac inlet flow

Grahic Jump Location
Fig. 14

Normalized deposited number on a bifurcating region of vasculature coated with different anti-ICAM densities for different particle sizes

Grahic Jump Location
Fig. 15

Mesh models over a cross section of a lung geometry. Mesh size is large in the middle and small near boundary. The minimum mesh size limits are 0.2, 0.15, 0.125, and 0.1 mm for (a), (b), (c), and (d) cases, respectively.

Grahic Jump Location
Fig. 16

Total number of deposited particle versus number of elements in different mesh models

Grahic Jump Location
Fig. 17

Mesh analysis using elements of four different sizes for (a) steady inlet flow and (b) unsteady inlet flow

Grahic Jump Location
Fig. 18

Mesh model of lung vasculature has been split into different regions. Gen: generation; Bif: bifurcation.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In