A computed tomography imaging-based subject-specific whole-lung deposition model

The respiratory tract is an important route for beneficial drug aerosol or harmful particulate matter to enter the body. To assess the therapeutic response or disease risk, whole-lung deposition models have been developed, but were limited by compartment, symmetry or stochastic approaches. In this work, we proposed an imaging-based subject-specific whole-lung deposition model. The geometries of airways and lobes were segmented from computed tomography (CT) lung images at total lung capacity (TLC), and the regional air-volume changes were calculated by registering CT images at TLC and functional residual capacity (FRC). The geometries were used to create the structure of entire subject-specific conducting airways and acinar units. The air-volume changes were used to estimate the function of subject-specific ventilation distributions among acinar units and regulate flow rates in respiratory airway models. With the airway dimensions rescaled to a desired lung volume and the airflow field simulated by a computational fluid dynamics model, particle deposition fractions were calculated using deposition probability formulae adjusted with an enhancement factor to account for the effects of secondary flow and airway geometry in proximal airways. The proposed model was validated in silico against existing whole-lung deposition models, three-dimensional (3D) computational fluid and particle dynamics (CFPD) for an acinar unit, and 3D CFPD deep lung model comprising conducting and respiratory regions. The model was further validated in vivo against the lobar particle distribution and the coefficient of variation of particle distribution obtained from CT and single-photon emission computed tomography (SPECT) images, showing good agreement. Subject-specific airway structure increased the deposition fraction of 10.0-μm particles and 0.01-μm particles by approximately 10%. An enhancement factor increased the overall deposition fractions, especially for particle sizes between 0.1 and 1.0 μm.