Quantification of particle deposition in asymmetrical tracheobronchial model geometry

Particle Deposition in Large-Scale Human Tracheobronchial Airways Predicted by Single-Path Modelling

The health effects of particles are directly related to their deposition patterns (deposition site and amount) in human airways. However, estimating the particle trajectory in a large-scale human lung airway model is still a challenge. In this work, a truncated single-path, large-scale human airway model (G3-G10) with a stochastically coupled boundary method were employed to investigate the particle trajectory and the roles of their deposition mechanisms. The deposition patterns of particles with diameters (dp) of 1-10 μm are investigated under various inlet Reynolds numbers (Re = 100-2000). Inertial impaction, gravitational sedimentation, and combined mechanism were considered. With the increasing airway generations, the deposition of smaller particles (dp < 4 μm) increased due to gravitational sedimentation, while that of larger particles decreased due to inertial impaction. The obtained formulas of Stokes number and Re can predict the deposition efficiency due to the combined mechanism in the present model, and the prediction can be used to assess the dose-effect of atmospheric aerosols on the human body. Diseases in deeper generations are mainly attributed to the deposition of smaller particles under lower inhalation rates, while diseases at the proximal generations mainly result from the deposition of larger particles under higher inhalation rates.

Anatomy matters: The role of the subject-specific respiratory tract on aerosol deposition - A CFD study

The COVID-19 pandemic is one of the greatest challenges to humanity nowadays. COVID-19 virus can replicate in the host's larynx region, which is in contrast to other viruses that replicate in lungs only, i.e. SARS. This is conjectured to support a fast spread of COVID-19. However, there is sparse research in this field about quantitative comparison of virus load in the larynx for varying susceptible individuals. In this regard the lung geometry itself could influence the risk of reproducing more pathogens and consequently exhaling more virus. Disadvantageously, there are only sparse lung geometries available. To still be able to investigate realistic geometrical deviations we employ three different digital replicas of human airways up to the 7 th level of bifurcation, representing two realistic lungs (male and female) as well as a more simplified experimental model. Our aim is to investigate the influence of breathing scenarios on aerosol deposition in anatomically different, realistic human airways. In this context, we employ three levels of cardiovascular activity as well as reported experimental particle size distributions by means of Computational Fluid Dynamics (CFD) with special focus on the larynx region to enable new insights into the local virus loads in human respiratory tracts. In addition, the influence of more realistic boundary conditions is investigated by performing transient simulations of a complete respiratory cycle in the upper lung regions of the considered respiratory models, focusing in particular on deposition in the oral cavity, the laryngeal region, and trachea, while simplifying the tracheobronchial tree. The aerosol deposition is modeled via OpenFOAM\protect \relax \special {t4ht=®} by employing an Euler-Lagrangian frame including steady and unsteady Reynolds Averaged Navier-Stokes (RANS) resolved turbulent flow using the k- ω -SST and k- ω -SST DES turbulence models. We observed that the respiratory geometry altered the local deposition patterns, especially in the laryngeal region. Despite the larynx region, the effects of varying flow rate for the airway geometries considered were found to be similar in the majority of respiratory tract regions. For all particle size distributions considered, localized particle accumulation occurred in the larynx of all considered lung models, which were more pronounced for larger particle size distributions. Moreover, it was found, that employing transient simulations instead of steady-state analysis, the overall particle deposition pattern is maintained, however with a stronger intensity in the transient cases.

Effects of hygroscopic growth of ambient urban aerosol particles on their modelled regional and local deposition in healthy and COPD-compromised human respiratory system

Total, regional and local deposition fractions of urban-type aerosol particles with diameters of 50, 75, 110 and 145 nm were modelled and studied in their dry state and after their hygroscopic growth using a Stochastic Lung Model and a Computational Fluid and Particle Dynamics method. Healthy subjects and patients with severe chronic obstructive pulmonary disease (COPD) were considered. The hygroscopic growth factors (HGFs) adopted were determined experimentally and represent a real urban-type environment. The hygroscopic growth of particles resulted in decrease of the deposition fractions in all major parts of the healthy respiratory system and the extent of the deposited fractions was rising monotonically with particle size. In the extrathoracic (ET) region, the relative decrease was between 7% and 13%. In the lungs the deposition decreased by 11-16%. The decrease of deposition fraction due to hygroscopic growth was more accentuated in the conductive airways (up to 25%) and less pronounced towards the terminal airways. The spatial distribution of the deposited particles remained highly inhomogeneous with some areas containing thousands times more particles than the average number of particles per unit surface area. For COPD patients, the hygroscopic growth produced similar deposition alterations in the ET region than for healthy subjects. In the conductive airways, however, the particle growth caused a substantial relative decrease in the deposition fractions. In contrast, the relative depositions of hygroscopic particles increased in the acinar region.

Simulation of the effect of mucociliary clearance on the bronchial distribution of inhaled radon progenies and related cellular damage using a new deposition and clearance model for the lung

Most of the current dosimetry models of inhaled short-lived radon decay products assume uniform activity distributions along the bronchial airways. In reality, however, both deposition and clearance patterns of inhaled radon progenies are highly inhomogeneous. Consequently, a new deposition-clearance model has been developed that accounts for such inhomogeneities and applied together with biophysical models of cell death and cell transformation. The scope of this study was to apply this model which is based on computational fluid and particle dynamics methods, in an effort to reveal the effect of mucociliary clearance on the bronchial distribution of deposited radon progenies. Furthermore, the influence of mucociliary clearance on the spatial distribution of biological damage due to alpha-decay of the deposited radon progenies was also studied. The results obtained demonstrate that both deposition and clearance of inhaled radon progenies are highly non-uniform within a human airway bifurcation unit. Due to the topology of the carinal ridge, a slow clearance zone emerged in this region, which is the location where most of the radio-aerosols deposit. In spite of the slow mucus movement in this zone, the initial degree of inhomogeneity of the activity due to the nonuniform deposition decreased by a factor of about 3 by considering the effect of mucociliary clearance. In the peak of the airway bifurcation, the computed cell death and cell transformation probabilities were lower when considering deposition and clearance simultaneously, compared to the case when only deposition was considered. However, cellular damage remained clustered.

Airflow and Particle Transport Prediction through Stenosis Airways

Airflow and particle transport in the human lung system is influenced by biological and other factors such as breathing pattern, particle properties, and deposition mechanisms. Most of the studies to date have analyzed airflow characterization and aerosol transport in idealized and realistic models. Precise airflow characterization for airway stenosis in a digital reference model is lacking in the literature. This study presents a numerical simulation of airflow and particle transport through a stenosis section of the airway. A realistic CT-scan-based mouth–throat and upper airway model was used for the numerical calculations. Three different models of a healthy lung and of airway stenosis of the left and right lung were used for the calculations. The ANSYS FLUENT solver, based on the finite volume discretization technique, was used as a numerical tool. Proper grid refinement and validation were performed. The numerical results show a complex-velocity flow field for airway stenosis, where airflow velocity magnitude at the stenosis section was found to be higher than that in healthy airways. Pressure drops at the mouth–throat and in the upper airways show a nonlinear trend. Comprehensive pressure analysis of stenosis airways would increase our knowledge of the safe mechanical ventilation of the lung. The turbulence intensities at the stenosis sections of the right and left lung were found to be different. Deposition efficiency (DE) increased with flow rate and particle size. The findings of the present study increase our understanding of airflow patterns in airway stenosis under various disease conditions. More comprehensive stenosis analysis is required to further improve knowledge of the field.

A Review of Respiratory Anatomical Development, Air Flow Characterization and Particle Deposition

The understanding of complex inhalation and transport processes of pollutant particles through the human respiratory system is important for investigations into dosimetry and respiratory health effects in various settings, such as environmental or occupational health. The studies over the last few decades for micro- and nanoparticle transport and deposition have advanced the understanding of drug-aerosol impacts in the mouth-throat and the upper airways. However, most of the Lagrangian and Eulerian studies have utilized the non-realistic symmetric anatomical model for airflow and particle deposition predictions. Recent improvements to visualization techniques using high-resolution computed tomography (CT) data and the resultant development of three dimensional (3-D) anatomical models support the realistic representation of lung geometry. Yet, the selection of different modelling approaches to analyze the transitional flow behavior and the use of different inlet and outlet conditions provide a dissimilar prediction of particle deposition in the human lung. Moreover, incorporation of relevant physical and appropriate boundary conditions are important factors to consider for the more accurate prediction of transitional flow and particle transport in human lung. This review critically appraises currently available literature on airflow and particle transport mechanism in the lungs, as well as numerical simulations with the aim to explore processes involved. Numerical studies found that both the Euler–Lagrange (E-L) and Euler–Euler methods do not influence nanoparticle (particle diameter ≤50 nm) deposition patterns at a flow rate ≤25 L/min. Furthermore, numerical studies demonstrated that turbulence dispersion does not significantly affect nanoparticle deposition patterns. This critical review aims to develop the field and increase the state-of-the-art in human lung modelling.

Euler-Lagrange Prediction of Diesel-Exhaust Polydisperse Particle Transport and Deposition in Lung: Anatomy and Turbulence Effects

In clinical assessments, the correlation between atmospheric air pollution and respiratory damage is highly complicated. Epidemiological studies show that atmospheric air pollution is largely responsible for the global proliferation of pulmonary disease. This is particularly significant, since most Computational Fluid Dynamics (CFD) studies to date have used monodisperse particles, which may not accurately reflect realistic inhalation patterns, since atmospheric aerosols are mostly polydisperse. The aim of this study is to investigate the anatomy and turbulent effects on polydisperse particle transport and deposition (TD) in the upper airways. The Euler-Lagrange approach is used for polydisperse particle TD prediction in both laminar and turbulent conditions. Various anatomical models are adopted to investigate the polydisperse particle TD under different flow conditions. Rossin-Rammler diameter distribution is used for the distribution of the initial particle diameter. The numerical results illustrate that airflow rate distribution at the right lung of a realistic model is higher than a non-realistic model. The CFD study also shows that turbulence effects on deposition are higher for larger diameter particles than with particles of smaller diameter. A significant amount of polydisperse particles are also shown to be deposited at the tracheal wall for CT-based model, whereas particles are mostly deposited at the carinal angle for the non-realistic model. A comprehensive, polydisperse particle TD analysis would enhance understanding of the realistic deposition pattern and decrease unwanted therapeutic aerosol deposition at the extrathoracic airways.

Polydisperse Microparticle Transport and Deposition to the Terminal Bronchioles in a Heterogeneous Vasculature Tree

The atmospheric particles from different sources, and the therapeutic particles from various drug delivery devices, exhibit a complex size distribution, and the particles are mostly polydisperse. The limited available in vitro, and the wide range of in silico models have improved understanding of the relationship between monodisperse particle deposition and therapeutic aerosol transport. However, comprehensive polydisperse transport and deposition (TD) data for the terminal airways is still unavailable. Therefore, to benefit future drug therapeutics, the present numerical model illustrates detailed polydisperse particle TD in the terminal bronchioles for the first time. Euler-Lagrange approach and Rosin-Rammler diameter distribution is used for polydisperse particles. The numerical results show higher deposition efficiency (DE) in the right lung. Specifically, the larger the particle diameter (d_p > 5 μm), the higher the DE at the bifurcation area of the upper airways is, whereas for the smaller particle (d_p < 5 μm), the DE is higher at the bifurcation wall. The overall deposition pattern shows a different deposition hot spot for different diameter particle. These comprehensive lobe-specific polydisperse particle deposition studies will increase understanding of actual inhalation for particle TD, which could potentially increase the efficiency of pharmaceutical aerosol delivery at the targeted position of the terminal airways.

Effect of morphology on nanoparticle transport and deposition in human upper tracheobronchial airways

Aerosol transport and deposition in human lungs has attracted considerable attention in the past few years, as it has significant value to the study of toxicity consequence as well as therapeutic potential in occupational health and medical applications. In reproducing human tracheobronchial airways, two approaches were frequently taken: (1) anatomical realistic reconstruction through image scans (e.g. CT and MRI) or cadaver casts and (2) mathematical description using simplified models. Strengths and limitations are primarily focused on accuracy, resolution, repeatability, and computational\physical expenses. While both approaches were reported in literature, detailed comparison of aerosol transport and deposition in the two representations were scarcely performed, largely due to the challenge to acquire comprehensive data from the irregular structured airway replicas (approach 1). To fill the gap, the current study performed a numerical comparison of nanoparticle transport and deposition in human upper tracheobronchial airways by using an anatomical realistic reconstruction (through CT scans) and a mathematically simplified airway model. As the first step, the current study was focused on the variation in breathing airflow pattern and the effect towards fate of the inhaled nanoparticles in human upper tracheobronchial airways. The study provided important information to understand geometric sensitivity of nanoparticle modeling in the human tracheobronchial tree and is of significant value to predict the whole lung uptake of inhaled nanoparticles in the human respiratory system.

The Creation and Statistical Evaluation of a Deterministic Model of the Human Bronchial Tree from HRCT Images

A quantitative description of the morphology of lung structure is essential prior to any form of predictive modeling of ventilation or aerosol deposition implemented within the lung. The human lung is a very complex organ, with airway structures that span two orders of magnitude and having a multitude of interfaces between air, tissue and blood. As such, current medical imaging protocols cannot provide medical practitioners and researchers with in-vivo knowledge of deeper lung structures. In this work a detailed algorithm for the generation of an individualized 3D deterministic model of the conducting part of the human tracheo-bronchial tree is described. Distinct initial conditions were obtained from the high-resolution computed tomography (HRCT) images of seven healthy volunteers. The algorithm developed is fractal in nature and is implemented as a self-similar space sub-division procedure. The expansion process utilizes physiologically realistic relationships and thresholds to produce an anatomically consistent human airway tree. The model was validated through extensive statistical analysis of the results and comparison of the most common morphological features with previously published morphometric studies and other equivalent models. The resulting trees were shown to be in good agreement with published human lung geometric characteristics and can be used to study, among other things, structure-function relationships in simulation studies.

Does the number of irradiated cells influence the spatial distribution of bystander effects?

There is growing evidence that the radiation effects at low doses are not adequately described by a simple linear extrapolation from high doses, due, among others, to bystander effects. Though several studies have been published on this topic, the explanation of the mechanisms describing the bystander effects remains unclear. This study aims at understanding how the bystander signals are or can be propagated in the cell culture, namely if the number of irradiated cells influences the bystander response. An A549 cell line was exposed to several doses of α-particles, being the bystander response quantified in two non-irradiated areas. The radius of irradiated areas differs by a factor of 2, and the non-irradiated areas were optimally designed to have the same number of cells. Our results show evidence for bystander effects occurring in cells far away from the irradiated ones, meaning that bystander signals can easily spread throughout the cell culture. Additionally, our study highlights that the damage caused by radiation on the surrounding of irradiated areas could be different according to the number of irradiated cells, i.e., for the same dose value; the overall cellular damage could be different.

Superiority of nebulized corticosteroids over dry powder inhalers in certain patients with cough variant asthma or cough-predominant asthma

BACKGROUND

The particle distribution might differ between nebulizer therapy and metered-dose inhaler (MDI) or dry powder inhaler (DPI) therapy because the particles repeatedly enter/re-enter the airways with the nebulizer. Inhaled corticosteroids (ICS) were administered with a nebulizer to assess the benefit of changes in the distribution of particles in patients with cough variant asthma (CVA) and cough-predominant asthma (CPA).

METHODS

Patients whose symptoms were not controlled by their current therapy were enrolled. In patients receiving high-dose ICS by MDI or DPI (ICS-MDI/DPI), steroid therapy was switched to 1,320μg/day of nebulized dexamethasone (1,600μg as dexamethasone sodium phosphate) (chronic steroid-independent group). In patients receiving systemic steroids regardless of their ICS-MDI/DPI therapy, nebulized dexamethasone was added and any concurrent ICS-MDI/DPI therapy was halted to detect a steroid-sparing effect (chronic steroid-dependent group). In patients with acute exacerbation of CVA or CPA and persistent symptoms despite systemic corticosteroids, nebulized dexamethasone was added to assess its effect (acute group).

RESULTS

Superior symptom control was achieved in 10 out of 12 steroid-independent patients, 3 out of 6 steroid-dependent patients, and all 7 acute patients.

CONCLUSIONS

Delivery of ICS via a nebulizer has advantages over ICS-MDI/DPI in some patients with CVA or CPA.

3D-modelling of radon-induced cellular radiobiological effects in bronchial airway bifurcations: direct versus bystander effects

PURPOSE

The primary objective of this paper was to investigate the distribution of radiation doses and the related biological responses in cells of a central airway bifurcation of the human lung of a hypothetical worker of the New Mexico uranium mines during approximately 12 hours of exposure to short-lived radon progenies.

MATERIALS AND METHODS

State-of-the-art computational modelling techniques were applied to simulate the relevant biophysical and biological processes in a central human airway bifurcation.

RESULTS

The non-uniform deposition pattern of inhaled radon daughters caused a non-uniform distribution of energy deposition among cells, and of related cell inactivation and cell transformation probabilities. When damage propagation via bystander signalling was assessed, it produced more cell killing and cell transformation events than did direct effects. If bystander signalling was considered, variations of the average probabilities of cell killing and cell transformation were supra-linear over time.

CONCLUSIONS

Our results are very sensitive to the radiobiological parameters, derived from in vitro experiments (e.g., range of bystander signalling), applied in this work and suggest that these parameters may not be directly applicable to realistic three-dimensional (3D) epithelium models.

Quantification of airway deposition of intact and fragmented pollens

Although pollen is one of the most widespread agents that can cause allergy, its airway transport and deposition is far from being fully explored. The objective of this study was to characterize the airway deposition of pollens and to contribute to the debate related to the increasing number of asthma attacks registered after thunderstorms. For the quantification of the deposition of inhaled pollens in the airways computer simulations were performed. Our results demonstrated that smaller and fragmented pollens may penetrate into the thoracic airways and deposit there, supporting the theory that fragmented pollen particles are responsible for the increasing incidence of asthma attacks following thunderstorms. Pollen deposition results also suggest that children are the most exposed to the allergic effects of pollens. Finally, pollens between 0.5 and 20 μm deposit more efficiently in the lung of asthmatics than in the healthy lung, especially in the bronchial region.

Effect of site-specific bronchial radon progeny deposition on the spatial and temporal distributions of cellular responses

Inhaled short-lived radon progenies may deposit in bronchial airways and interact with the epithelium by the emission of alpha particles. Simulation of the related radiobiological effects requires the knowledge of space and time distributions of alpha particle hits and biological endpoints. Present modelling efforts include simulation of radioaerosol deposition patterns in a central bronchial airway bifurcation, modelling of human bronchial epithelium, generation of alpha particle tracks, and computation of spatio-temporal distributions of cell nucleus hits, cell killing and cell transformation events. Simulation results indicate that the preferential radionuclide deposition at carinal ridges plays an important role in the space and time evolution of the biological events. While multiple hits are generally rare for low cumulative exposures, their probability may be quite high at the carinal ridges of the airway bifurcations. Likewise, cell killing and transformation events also occur with higher probability in this area. In the case of uniform surface activities, successive hits as well as cell killing and transformation events within a restricted area (say 0.5 mm(2)) are well separated in time. However, in the case of realistic inhomogeneous deposition, they occur more frequently within the mean cycle time of cells located at the carinal ridge even at low cumulative doses. The site-specificity of radionuclide deposition impacts not only on direct, but also on non-targeted radiobiological effects due to intercellular communication. Incorporation of present results into mechanistic models of carcinogenesis may provide useful information concerning the dose-effect relationship in the low-dose range.

Cellular burdens and biological effects on tissue level caused by inhaled radon progenies

In the case of radon exposure, the spatial distribution of deposited radioactive particles is highly inhomogeneous in the central airways. The object of this research is to investigate the consequences of this heterogeneity regarding cellular burdens in the bronchial epithelium and to study the possible biological effects at tissue level. Applying computational fluid and particle dynamics techniques, the deposition distribution of inhaled radon daughters has been determined in a bronchial airway model for 23 min of work in the New Mexico uranium mine corresponding to 0.0129 WLM exposure. A numerical epithelium model based on experimental data has been utilised in order to quantify cellular hits and doses. Finally, a carcinogenesis model considering cell death-induced cell-cycle shortening has been applied to assess the biological responses. Present computations reveal that cellular dose may reach 1.5 Gy, which is several orders of magnitude higher than tissue dose. The results are in agreement with the histological finding that the uneven deposition distribution of radon progenies may lead to inhomogeneous spatial distribution of tumours in the bronchial airways. In addition, at the macroscopic level, the relationship between cancer risk and radiation burden seems to be non-linear.

Regional deposition of particles in an image-based airway model: large-eddy simulation and left-right lung ventilation asymmetry

Regional deposition and ventilation of particles by generation, lobe and lung during steady inhalation in a computed tomography (CT) based human airway model are investigated numerically. The airway model consists of a seven-generation human airway tree, with oral cavity, pharynx and larynx. The turbulent flow in the upper respiratory tract is simulated by large-eddy simulation. The flow boundary conditions at the peripheral airways are derived from CT images at two lung volumes to produce physiologically-realistic regional ventilation. Particles with diameter equal to or greater than 2.5 microns are selected for study because smaller particles tend to penetrate to the more distal parts of the lung. The current generational particle deposition efficiencies agree well with existing measurement data. Generational deposition efficiencies exhibit similar dependence on particle Stokes number regardless of generation, whereas deposition and ventilation efficiencies vary by lobe and lung, depending on airway morphology and airflow ventilation. In particular, regardless of particle size, the left lung receives a greater proportion of the particle bolus as compared to the right lung in spite of greater flow ventilation to the right lung. This observation is supported by the left-right lung asymmetry of particle ventilation observed in medical imaging. It is found that the particle-laden turbulent laryngeal jet flow, coupled with the unique geometrical features of the airway, causes a disproportionate amount of particles to enter the left lung.

New developments in aerosol dosimetry

Dosimetry provides information linking environmental exposures to sites of deposition, removal from these sites, and translocation of deposited materials. Dosimetry also aids in extrapolating laboratory animal and in vitro data to humans. Recent progress has shed light on: properties of particles in relation to their fates in the body; influence of age, gender, body size, and lung diseases on inhaled particle doses; particle movement to the brain via the olfactory nerves; and particle deposition hot spots in the respiratory tract. Ultrafine size has emerged as an important dosimetric characteristic. Particle count, composition, and surface properties are recognized as potentially important toxicology-related considerations. Differences in body size influence airway sizes, inhaled particle deposition, specific ventilation, and specific doses (e.g. per unit body mass). Related to body size, age, gender, species, and strain are also dosimetric considerations. Diseases, such as chronic obstructive pulmonary disease (COPD) and bronchitis, produce uneven doses within the respiratory tract. Traditional concepts of the translocation and clearance of deposited particles have been challenged. Ultrafine particles can translocate to the brain via olfactory nerves, and from the lung to other organs. The clearance rates of particles from tracheobronchial airways are slowed by respiratory tract infections, but newer evidence implies that slow particle clearance from this region also exists in healthy lungs. Finally, hot spots of particle deposition are seen in hollow models, lung tissue, and dosimetric simulations. Local doses to groups of epithelial cells can be much greater than those to surrounding cells. The new insights challenge dosimetry scientists.

Modeling the dynamics of airway constriction: effects of agonist transport and binding

Recent advances have revealed that during exogenous airway challenge, airway diameters cannot be adequately predicted by their initial diameters. Furthermore, airway diameters can also vary greatly in time on scales shorter than a breath. To better understand these phenomena, we developed a multiscale model that allowed us to simulate aerosol challenge in the airways during ventilation. The model incorporates agonist-receptor binding kinetics to govern the temporal response of airway smooth muscle contraction on individual airway segments, which, together with airway wall mechanics, determines local airway caliber. Global agonist transport and deposition are coupled with pressure-driven flow, linking local airway constrictions with global flow dynamics. During the course of challenge, airway constriction alters the flow pattern, redistributing the agonist to less constricted regions. This results in a negative feedback that may be a protective property of the normal lung. As a consequence, repetitive challenge can cause spatial constriction patterns to evolve in time, resulting in a loss of predictability of airway diameters. Additionally, the model offers new insights into several phenomena including the intra- and interbreath dynamics of airway constriction throughout the tree structure.