Particle deposition in a CT-scanned human lung airway

Effect of laryngeal jet on dry powder inhaler aerosol deposition: a numerical simulation

Although pulmonary drug delivery has been deeply investigated, the effect of the laryngeal jet on particle deposition during drug delivery with dry powder inhalers (DPI) has not been evaluated profoundly. In this study, the flow structure and particle deposition pattern of a DPI in two airway models, one with mouth-throat region including the larynx and the other one without it, are numerically investigated. The results revealed that the laryngeal jet has a considerable effect on particle deposition. The presence of laryngeal jet leads to 4-fold and 2-fold higher deposition efficiencies for inlet flow rates of 30 and 90 L/min, respectively.

Dosimetry Sensitivity in a Lower Dimensional Model of Patient-Specific Asthma Subjects

OBJECTIVE

Experimental uncertainty will impact in silico model calculations of aerosol delivery and deposition. Patient-specific dosimetry models are often parameterized based on medical imaging data, which contain inherent experimental variability.

METHODS

Here, we created and parameterized 1D models of three subject-specific asthmatic subjects and randomly assigned perturbations of up to 15 % on airway diameter, segmental volume, and defected volume. Sensitivity of imaging data experimental variability on dosimetry metrics were quantified.

RESULTS

Lobar particle delivery primarily depended on the distal segmental volumes; 15 % range of noise resulted in delivery to the upper right lobe to vary at most from 15.2 and 18.2 % for one of the severe subjects. Particle deposition was most sensitive to airway diameter; 95 % confidence intervals spanned from 8 to 10.6 % in the mild/moderate subject for 15 % variation on input metrics for 5 [Formula: see text] diameter particles. While these results provide possible ranges of dosimetry calculations for a specific subject, the perturbations were not sufficient to model the large observed inter-subject variability (8.9, 19, and 14.5 % deposition, subjects 1--3, respectively, 5 [Formula: see text] diameter particles).

CONCLUSION

This study highlights that in silico model predictions are robust in the presence of experimental uncertainty and that it continues to be necessary to perform subject-specific simulations, especially within the presence of heterogeneous airway disease.

SIGNIFICANCE

Sensitivity analysis provides confidence in calculating deposition in the airways of asthmatic subjects within the presence of experimental uncertainty.

A Numerical Simulation of the Airflow and Aerosol Particle Deposition in a Realistic Airway Model of a Healthy Adult

Determining the behavior of aerosol drug particles is of vital importance in the treatment of respiratory tract diseases. Despite the development of imaging techniques in the pulmonary region in recent years, current imaging techniques are insufficient to detect particle deposition. Computational fluid dynamics (CFD) methods can fill the gap in this field as they take into account the very different physical processes that occur during aerosol transport. This study aims to numerically investigate the airflow and the aerosol particle dynamics on a realistic human respiratory tract model during multiple breathing cycles. The simulations were conducted on the different breathing conditions for people under light, normal, and heavy physical activities, and the aerosol particles with different aerodynamic diameters (i.e., dp=2, 5, and 7 µm). The numerical results were validated by comparing extensively with experimental and numerical results. The results indicated that the airflow during inspiration and expiration was characteristically different from each other and changed with the inspiration flow rate. It was determined that small-sized particles followed the streamlines and moved towards the distal of the lung under low respiratory conditions. On the other hand, larger particles tended to deposit in higher generations due to the higher inertia. It was found that with the increase of inspiration flow rate the deposition of particles increased for all particles during multiple breaths. For light breathing conditions, low deposition efficiencies were obtained because the particles followed the streamlines and moved towards the distal part of the lung. The particle deposition efficiency under heavy breathing conditions was 28.2% for 2 µm, 33.05% for 5 µm, and 38.4% for 7 µm particles. The results showed that inertial impaction plays an active role in particle deposition.

Inter-species Variabilities of Droplet Transport, Size Change, and Deposition in Human and Rat Respiratory Systems: An In Silico Study

To speculate on human responses from animal studies, scale-up factors (body weight, lung volume, or lung surface area ratios) are currently used to extrapolate aerosol lung deposition from animal to human. However, those existing scale-up methods between animals and humans neglected two important inter-subject variability factors: (1) the effect of anatomical differences in respiratory systems from mouth/nose to peripheral lungs between human and rat, and (2) the effect of spatial distributions and temporal evolutions of temperature and relative humidity (RH) on droplet size change dynamics between the two species. To test the above-mentioned inter-species variability effects on droplet fates in pulmonary routes and generate correlations as a precise scale-up method for lung deposition estimation, this study simulated the transport of pure-water droplets in both human and Sprague-Dawley (SD) rat respiratory systems. Employing an experimentally validated Euler-Lagrange based Computational Fluid-Particle Dynamics (CFPD) model, simulations were performed for droplets with Stk/Fr between 8.36×10^-5 and 1.25×10^-2. Droplets were inhaled through human and rat nostrils with resting breathing conditions. Numerical results indicate that RH becomes uniformly distributed in rat airways sooner than in human airways, which significantly influences droplet size change dynamics and the resultant trajectories in pulmonary paths. Using the Stokes-Froude dimensionless number group (i.e., Stk/Fr) as the independent variable, the regional deposition fractions and evaporation fractions in both rat and human respiratory systems collapsed into unified correlations. The correlations can be used as a new rat-to-human scale-up method, estimating the lung depositions with consideration of anatomical differences. Furthermore, the necessity to employ realistic RH and temperature boundary conditions at airway walls was also confirmed for the accurate prediction of droplet size change using CFPD. Employing idealized boundary conditions leads the droplets to evaporate slower and deposit more than using realistic RH and temperature boundary conditions.

Investigation of airflow at different activity conditions in a realistic model of human upper respiratory tract

In the present study, the turbulent flows inside a realistic model of the upper respiratory tract were investigated numerically and experimentally. The airway model included the geometrical details of the oral cavity to the end of the trachea that was based on a series of CT-scan images. The topological data of the respiratory tract were used for generating the computational model as well as the 3D-printed model that was used in the experimental pressure drop measurement. Different airflow rates of 30, 45, and 60 L/min, which correspond to the light, semi-light, and heavy activity breathing conditions, were investigated numerically using turbulence and transition models, as well as experimentally. Simulation results for airflow properties, including velocity vectors, pressure drops, streamlines, eddy viscosity, and turbulent kinetic energy contours in the oral-trachea airway model, were presented. The simulated pressure drop was compared with the experimental data, and reasonable agreement was found. The obtained results showed that the maximum pressure drop occurs in the narrowest part of the larynx region. A comparison between the numerical results and experimental data showed that the transition (γ-Reθ) SST model predicts higher pressure losses, especially at higher breathing rates. Formations of the secondary flows in the oropharynx and trachea regions were also observed. In addition, the simulation results showed that in the trachea region, the secondary flow structures dissipated faster for the flow rate of 60 L/min compared to the lower breathing rates of 30 and 45 L/min.

Patient-specific modeling of aerosol delivery in healthy and asthmatic adults

The magnitude and regional heterogeneity of airway obstructions in severe asthmatics is likely linked to insufficient drug delivery, as evidenced by the inability to mitigate exacerbations with inhaled aerosol medications. To understand the correlation between morphometric features, airflow distribution, and inhaled dosimetry, we perform dynamic computational simulations in two healthy and four asthmatic subjects. Models incorporate computed tomography-based and patient-specific central airway geometries and hyperpolarized ³He MRI-measured segmental ventilation defect percentages (SVDPs), implemented as resistance boundary conditions. Particles [diameters (d_p) = 1, 3, and 5 μm] are simulated throughout inhalation, and we record their initial conditions, both spatially and temporally, with their fate in the lung. Predictions highlight that total central airway deposition is the same between the healthy subjects (26.6%, d_p = 3 μm) but variable among the asthmatic subjects (ranging from 5.9% to 59.3%, d_p = 3 μm). We found that by preferentially releasing the particles during times of fast or slow inhalation rates we enhance either central airway deposition percentages or peripheral particle delivery, respectively. These predictions highlight the potential to identify with simulations patients who may not receive adequate therapeutic dosages with inhaled aerosol medication and therefore identify patients who may benefit from alternative treatment strategies. Furthermore, by improving regional dose levels, we may be able to preferentially deliver drugs to the airways in need, reducing associated adverse side effects.NEW & NOTEWORTHY Although it is evident that exacerbation mitigation is unsuccessful in some asthmatics, it remains unclear whether or not these patients receive adequate dosages of inhaled therapeutics. By coupling MRI and computed tomography data with patient-specific computational models, our predictions highlight the large intersubject variability, specifically in severe asthma.

Regional flow and deposition variability in adult female lungs: A numerical simulation pilot study

BACKGROUND

Despite the promise of respiratory simulations improving diagnosis and treatment of pulmonary diseases, model predictions have yet to be translated into the clinical setting. Current state-of-the-art in silico models have not yet incorporated subject variability in their predictions of airflow distributions and extent of deposited particles. Until inter-subject variability is accounted for in lung modeling, it will remain impossible to translate model predictions into clinical practice.

METHODS

Airflow and particle trajectories (d_p=1,3,5μm) are calculated in three subject-specific female adults by performing physiologically-based simulations. The computation framework features the ability to track air and particles throughout the respiration cycle and in the entire lung. Airway resistances, air velocities, and local deposition sites are correlated to airway anatomical features.

FINDINGS

Smaller airway diameters are correlated to larger airway resistances and pressure gradients in one subject compared to the other two. Irregular shape of the airway and flow direction (e.g. inspiration or expiration) correspond with peak velocities and secondary flow motions. Largest subject variability in deposition between conducting and respiratory zones is seen for 1 μm diameter particles. Little difference in total deposition is found among subjects. Localized deposited particle concentration hotspots are linked to airway anatomy and flow motion.

INTERPRETATION

Simulation predictions provide a first look into the correlation of anatomical features with airflow characteristics and deposited particle concentrations. Global deposition percentages ranged (at most, by 20%) between subjects and variances in localized deposition hotspots are correlated to variances in flow characteristics.

The Effectiveness of An Averaged Airway Model in Predicting the Airflow and Particle Transport Through the Airway

Background: In this study, we proposed an averaged airway model design based on four healthy subjects and numerically evaluated its effectiveness for predicting the airflow and particle transport through an airway. Methods: Direct-averaged models of the conducting airways of four subjects were restored by averaging the three-dimensional (3D) skeletons of four healthy airways, which were calculated using an inverse 3D thinning algorithm. We simulated the airflow and particle transport in the individual and the averaged airway models using computational fluid dynamics. Results: The bifurcation geometry differs even among healthy subjects, but the averaged model retains the typical geometrical characteristics of the airways. The Reynolds number of the averaged model varied within the range found in the individual subject models, and the averaged model had similar inspiratory flow characteristics as the individual subject models. The deposition fractions at almost all individual lobes ranged within the variation observed in the subjects, however, the deposition fraction was higher in only one lobe. The deposition distribution at the main bifurcation point differed among the healthy subjects, but the characteristics of the averaged model fell within the variation observed in the individual subject models. On the contrary, the deposition fraction of the averaged model was higher than that of the average of the individual subject models and deviated from the range observed in the subject models. Conclusion: These results indicate that the direct-averaged model may be useful for predicting the individual airflow and particle transport on a macroscopic scale.

Detailed computational analysis of flow dynamics in an extended respiratory airway model

BACKGROUND

Understanding respiratory physiology can aid clinicians in diagnosing the cause of respiratory symptoms or shed light on drug delivery inhaler device optimisation. However, the sheer complexity of the human lung prohibits a full-scale study.

METHODS

In this study, a realistic respiratory airway model including large-to-small conducting airways was built. This airway model consists of subject-specific upper and lower airways, extending from nasal and oral openings to terminal bronchioles (up to the 15th generation). Based on the subject-specific airway model, topological information was extracted and a digital reference model that exhibits strong asymmetry and multi-fractal properties was provided. Inhalation flow rates 18 L/min and 50 L/min were adopted to understand inspiratory conditions subjecting to resting and light exercise inhalation modes. Regional airflow in terms of axial velocity and secondary flow vortices along the lung airway model was extracted.

FINDINGS

Obvious secondary flow currents were seen in the larynx-trachea segment and left main bronchus, while for the terminal conducting airway in the right lower lobe, the airflow tends to be much smoother with no secondary flow currents.

INTERPRETATION

This paper provides insights on respiratory physiology, especially in the lower lung airways, and will be potentially useful for diagnosis of lower airway diseases.

Targeting inhaled aerosol delivery to upper airways in children: Insight from computational fluid dynamics (CFD)

Despite the prevalence of inhalation therapy in the treatment of pediatric respiratory disorders, most prominently asthma, the fraction of inhaled drugs reaching the lungs for maximal efficacy remains adversely low. By and large drug delivery devices and their inhalation guidelines are typically derived from adult studies with child dosages adapted according to body weight. While it has long been recognized that physiological (e.g. airway sizes, breathing maneuvers) and physical transport (e.g. aerosol dynamics) characteristics are critical in governing deposition outcomes, such knowledge has yet to be extensively adapted to younger populations. Motivated by such shortcomings, the present work leverages in a first step in silico computational fluid dynamics (CFD) to explore opportunities for augmenting aerosol deposition in children based on respiratory physiological and physical transport determinants. Using an idealized, anatomically-faithful upper airway geometry, airflow and aerosol motion are simulated as a function of age, spanning a five year old to an adult. Breathing conditions mimic realistic age-specific inhalation maneuvers representative of Dry Powder Inhalers (DPI) and nebulizer inhalation. Our findings point to the existence of a single dimensionless curve governing deposition in the conductive airways via the dimensionless Stokes number (Stk). Most significantly, we uncover the existence of a distinct deposition peak irrespective of age. For the DPI simulations, this peak (∼ 80%) occurs at Stk ≈ 0.06 whereas for nebulizer simulations, the corresponding peak (∼ 45%) occurs in the range of Stk between 0.03-0.04. Such dimensionless findings hence translate to an optimal window of micron-sized aerosols that evolves with age and varies with inhalation device. The existence of such deposition optima advocates revisiting design guidelines for optimizing deposition outcomes in pediatric inhalation therapy.

Effect of upper airway on tracheobronchial fluid dynamics

The upper airways play a significant role in the tracheal flow dynamics. Despite many previous studies, however, the effect of the upper airways on the ventilation distribution in distal airways has remained a challenge. The aim of this study is to experimentally and computationally investigate the dynamic behaviour in the intratracheal flow induced by the upper respiratory tract and to assess its influence on the subsequent tributaries. Patient-specific images from 2 different modalities (magnetic resonance imaging of the upper airways and computed tomography of the lower airways) were segmented and combined. An experimental phantom of patient-specific airways (including the oral cavity, larynx, trachea, down to generations 6-8) was generated using 3D printing. The flow velocities in this phantom model were measured by the flow-sensitised phase contrast magnetic resonance imaging technique and compared with the computational fluid dynamics simulations. Both experimental and computational results show a good agreement in the time-averaged velocity fields as well as fluctuating velocity. The flows in the proximal trachea were complex and unsteady under both lower- and higher-flow rate conditions. Computational fluid dynamics simulations were also performed with an airways model without the upper airways. Although the flow near the carina remained unstable only when the inflow rate was high, the influence of the upper airways caused notable changes in distal flow distributions when the 2 airways models were compared with and without the upper airways. The results suggest that the influence of the upper airways should be included in the respiratory flow assessment as the upper airways extensively affect the flows in distal airways and consequent ventilation distribution in the lungs.

Particle Disposition in the Realistic Airway Tree Models of Subjects with Tracheal Bronchus and COPD

Dispositions of inhalable particles in the human respiratory tract trigger and exacerbate airway inflammatory diseases. However, the particle deposition (PD) in airway of subjects with tracheal bronchus (TB) and chronic obstructive pulmonary diseases (COPD) is unknown. We therefore propose to clarify the disrupted PD associated with TB and COPD using the computational fluid dynamics (CFD) simulation. Totally nine airway tree models are included. Six are extracted from CT images of different individuals (two with TB, two with COPD, and two healthy controls (HC)). The others are the artificially modified models (AMMs) generated by the virtual lesion. Specifically, they are constructed through artificially adding a tracheal bronchus or a stenosis on one HC model. The deposition efficiency (DE) and deposition fraction (DF) in these models are obtained by the Euler-Lagrange approach, analyzed, and compared across models, locations, and particle sizes (0.1-10.0 micrometers). It is found that the PD in models with TB and COPD has been disrupted by the geometrical changes and followed airflow alternations. DE of the tracheal bronchus is higher for TB models. For COPD, the stenosis location determines the effects on DE and DF. Higher DF at the trachea is observed in TB1, TB2, and COPD2 models. DE increases with the particle size, and DE of the terminal bronchi is higher than that of central regions. Combined with AMMs, the CFD simulation using realistic airway models demonstrates disruptions of DP. The methods and findings might help understand the etiology of pulmonary diseases and improve the efficacy of inhaled medicines.

Ultrafine particle transport and deposition in a large scale 17-generation lung model

To understand how to assess optimally the risks of inhaled particles on respiratory health, it is necessary to comprehend the uptake of ultrafine particulate matter by inhalation during the complex transport process through a non-dichotomously bifurcating network of conduit airways. It is evident that the highly toxic ultrafine particles damage the respiratory epithelium in the terminal bronchioles. The wide range of in silico available and the limited realistic model for the extrathoracic region of the lung have improved understanding of the ultrafine particle transport and deposition (TD) in the upper airways. However, comprehensive ultrafine particle TD data for the real and entire lung model are still unavailable in the literature. Therefore, this study is aimed to provide an understanding of the ultrafine particle TD in the terminal bronchioles for the development of future therapeutics. The Euler-Lagrange (E-L) approach and ANSYS fluent (17.2) solver were used to investigate ultrafine particle TD. The physical conditions of sleeping, resting, and light activity were considered in this modelling study. A comprehensive pressure-drop along five selected path lines in different lobes was calculated. The non-linear behaviour of pressure-drops is observed, which could aid the health risk assessment system for patients with respiratory diseases. Numerical results also showed that ultrafine particle-deposition efficiency (DE) in different lobes is different for various physical activities. Moreover, the numerical results showed hot spots in various locations among the different lobes for different flow rates, which could be helpful for targeted therapeutical aerosol transport to terminal bronchioles and the alveolar region.

Fluid flow and particle transport in mechanically ventilated airways. Part II: particle transport

The flow mechanisms that play a role on aerosol deposition were identified and presented in a companion paper (Timothy et al. in Med Biol Eng Comput. doi: 10.1007/s11517-015-1407-3 , 2015). In the current paper, the effects of invasive conventional mechanical ventilation waveforms and endotracheal tube (ETT) on the aerosol transport were investigated. In addition to the enhanced deposition seen at the carinas of the airway bifurcations, enhanced deposition was also seen in the right main bronchus due to impaction and turbulent dispersion resulting from the fluid structures created by jet caused by the ETT. The orientation of the ETT toward right bronchus resulted in a substantial deposition inside right lung compared to left lung. The deposition inside right lung was ~12-fold higher than left lung for all considered cases, except for the case of using pressure-controlled sinusoidal waveform where a reduction of this ratio by ~50 % was found. The total deposition during pressure constant, volume ramp, and ascending ramp waveforms was similar and ~1.44 times higher than deposition fraction when using pressure sinusoidal waveform. Varying respiratory waveform demonstrated a significant role on the deposition enhancement factors and give evidence of drug aerosol concentrations in key deposition sites, which may be significant for drugs with negative side effects in high concentrations. These observations are thought to be important for ventilation treatment strategy.

The role of coupled resistance-compliance in upper tracheobronchial airways under high frequency oscillatory ventilation

A large eddy simulation (LES) based computational fluid dynamics (CFD) study was conducted to investigate lung lobar ventilation and gas exchange under high frequency oscillatory ventilation conditions. Time-dependent pressure coupled with the airways resistance and compliance (R&C) were imposed as boundary conditions (BCs) in the upper tracheobronchial tree of patient-specific lung geometry. The flow distribution in the left and right lungs demonstrated significant variations compared to the case in which traditional BCs based on mass flow rate fractions was used and is in agreement with the in vivo data available in the literature. The gas transport due to the pendelluft mechanism was captured in the different lung lobes and units. The computed pendelluft elapsed time was consistent with available physiological data. In contrast to in vivo studies, our simulations were able to predict the volume associated with the pendelluft elapsed time at different frequencies. Significant differences in coaxial counter flow and flow structures were observed between different BCs. The consistency of the results with the physiological in vivo data indicates that computations with coupled R&C BCs provide a suitable alternative tool for understanding the gas transport, diagnosing lung pathway disease severity, and optimizing ventilation management techniques.

Flow transport and gas mixing during invasive high frequency oscillatory ventilation

A large Eddy simulation (LES) based computational fluid dynamics study was performed to investigate gas transport and mixing in patient specific human lung models during high frequency oscillatory ventilation. Different pressure-controlled waveforms (sinusoidal, exponential and square) and ventilator frequencies (15, 10 and 6Hz) were used (tidal volume=50mL). The waveforms were created by solving the equation of motion subjected to constant lung wall compliance and flow resistance. Simulations were conducted with and without endotracheal tube to understand the effect of invasive management device. Variation of pressure-controlled waveform and frequency exhibits significant differences on counter flow pattern, which could lead to a significant impact on the gas mixing efficiency. Pendelluft-like flow was present for the sinusoidal waveform at all frequencies but occurred only at early inspiration for the square waveform at highest frequency. The square waveform was most efficient for gas mixing, resulting in the least wall shear stress on the lung epithelium layer thereby reducing the risk of barotrauma to both airways and the alveoli for patients undergoing therapy.

Airway gas flow

Local characteristics of airflow and its global distribution in the lung are determined by interaction between resistance to flow through the airways and the compliance of the tissue, with tissue compliance dominating flow distribution in the healthy lung. Current understanding is that conceptualizing the airways of the lung as a system of smooth adjoined cylinders through which air traverses laminarly is insufficient for understanding flow and energy dissipation and is particularly poor for predicting physiologically realistic transport of particles by the airflow. With rapid advances in medical imaging, computer technologies, and computational techniques, computational fluid dynamics is now becoming a viable tool for providing detailed information on the mechanics of airflow in the human respiratory tract. Studies using such techniques have shown that the upper airway (specifically its development of a turbulent laryngeal jet in the trachea), airway geometry, branching and rotation angle, and the pattern of joining of successive bifurcations are important in determining airflow structures. It is now possible to compute airflow in physical domains that are anatomically accurate and subject specific, enabling comparisons among intersubjects, that among subjects of different ages, and that among different species.

A multiscale MDCT image-based breathing lung model with time-varying regional ventilation

A novel algorithm is presented that links local structural variables (regional ventilation and deforming central airways) to global function (total lung volume) in the lung over three imaged lung volumes, to derive a breathing lung model for computational fluid dynamics simulation. The algorithm constitutes the core of an integrative, image-based computational framework for subject-specific simulation of the breathing lung. For the first time, the algorithm is applied to three multi-detector row computed tomography (MDCT) volumetric lung images of the same individual. A key technique in linking global and local variables over multiple images is an in-house mass-preserving image registration method. Throughout breathing cycles, cubic interpolation is employed to ensure C₁ continuity in constructing time-varying regional ventilation at the whole lung level, flow rate fractions exiting the terminal airways, and airway deformation. The imaged exit airway flow rate fractions are derived from regional ventilation with the aid of a three-dimensional (3D) and one-dimensional (1D) coupled airway tree that connects the airways to the alveolar tissue. An in-house parallel large-eddy simulation (LES) technique is adopted to capture turbulent-transitional-laminar flows in both normal and deep breathing conditions. The results obtained by the proposed algorithm when using three lung volume images are compared with those using only one or two volume images. The three-volume-based lung model produces physiologically-consistent time-varying pressure and ventilation distribution. The one-volume-based lung model under-predicts pressure drop and yields un-physiological lobar ventilation. The two-volume-based model can account for airway deformation and non-uniform regional ventilation to some extent, but does not capture the non-linear features of the lung.

Flow oscillation - a measure to predict the surgery outcome for obstructed sleep apnea (OSA) subject

Obstructed sleep apnea (OSA) is a common disorder which may need surgery to widen the airway; however the success rate of surgery is limited. Here we report a finding that could be used to predict the outcome of the OSA surgery. We found that inspiratory flow oscillates due to flow separation near the larynx, and the resulting periodic signal (3-5Hz) is an intrinsic property of breathing. This flow oscillating signal may be the afferent stimulus to trigger respiratory events. It is found that the flow oscillation is attenuated for the OSA subjects. The computational fluid dynamics (CFD) simulation reveals that there exists flow separation near larynx and this separation is severely weakened in the OSA upper airway model. It is believed that the flow oscillating signal can serve as the measure to quantify the breathing quality of an OSA subject. This makes it possible to predict the surgery outcome of the OSA subject by applying CFD simulation.