CFD Simulation of Aerosol Deposition in an Anatomically Based Human Large–Medium Airway Model

Automated bidirectional coupling of multiscale models of aerosol dosimetry: validation with subject-specific deposition data

Assessing the toxicity of airborne particulate matter or the efficacy of inhaled drug depends upon accurate estimates of deposited fraction of inhaled materials. In silico approaches can provide important insights into site- or airway-specific deposition of inhaled aerosols in the respiratory system. In this study, we improved on our recently developed 3D/1D model that simulate aerosol transport and deposition in the whole lung over multiple breath cycles (J. Aerosol Sci 151:105647). A subject-specific multiscale lung model of a healthy male subject using computational fluid-particle dynamics (CFPD) in a 3D model of the oral cavity through the large bronchial airways entering each lobe was bidirectionally coupled with a recently improved Multiple Path Particle Dosimetry (MPPD) model to predict aerosol deposition over the entire respiratory tract over multiple breaths for four conditions matching experimental aerosol exposures in the same subject from which the model was developed. These include two particle sizes (1 and 2.9 μm) and two subject-specific breathing rates of ~300 ml/s (slow breathing) and ~750 ml/s (fast breathing) at a target tidal volume of 1 L. In silico predictions of retained fraction were 0.31 and 0.29 for 1 μm and 0.66 and 0.62 for 2.9 μm during slow and fast breathing, respectively, and compared well with experimental data (1 μm: 0.31±0.01 (slow) and 0.27±0.01 (fast), 2.9 μm: 0.63±0.03 (slow) and 0.68±0.02 (fast)). These results provide a great deal of confidence in the validity and reliability of our approach.

In vitro-in silico correlation of three-dimensional turbulent flows in an idealized mouth-throat model

There exists an ongoing need to improve the validity and accuracy of computational fluid dynamics (CFD) simulations of turbulent airflows in the extra-thoracic and upper airways. Yet, a knowledge gap remains in providing experimentally-resolved 3D flow benchmarks with sufficient data density and completeness for useful comparison with widely-employed numerical schemes. Motivated by such shortcomings, the present work details to the best of our knowledge the first attempt to deliver in vitro-in silico correlations of 3D respiratory airflows in a generalized mouth-throat model and thereby assess the performance of Large Eddy Simulations (LES) and Reynolds-Averaged Numerical Simulations (RANS). Numerical predictions are compared against 3D volumetric flow measurements using Tomographic Particle Image Velocimetry (TPIV) at three steady inhalation flowrates varying from shallow to deep inhalation conditions. We find that a RANS k-ω SST model adequately predicts velocity flow patterns for Reynolds numbers spanning 1'500 to 7'000, supporting results in close proximity to a more computationally-expensive LES model. Yet, RANS significantly underestimates turbulent kinetic energy (TKE), thus underlining the advantages of LES as a higher-order turbulence modeling scheme. In an effort to bridge future endevours across respiratory research disciplines, we provide end users with the present in vitro-in silico correlation data for improved predictive CFD models towards inhalation therapy and therapeutic or toxic dosimetry endpoints.

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.

Fluid dynamic assessment of positive end-expiratory pressure in a tracheostomy tube connector during respiration

High-flow oxygen therapy using a tracheostomy tube is a promising clinical approach to reduce the work of breathing in tracheostomized patients. Positive end-expiratory pressure (PEEP) is usually applied during oxygen inflow to improve oxygenation by preventing end-expiratory lung collapse. However, much is still unknown about the geometrical effects of PEEP, especially regarding tracheostomy tube connectors (or adapters). Quantifying the degree of end-expiratory pressure (EEP) that takes patient-specific spirometry into account would be useful in this regard, but no such framework has been established yet. Thus, a platform to assess PEEP under respiration was developed, wherein three-dimensional simulation of airflow in a tracheostomy tube connector is coupled with a lumped lung model. The numerical model successfully reflected the magnitude of EEP measured experimentally using a lung phantom. Numerical simulations were further performed to quantify the effects of geometrical parameters on PEEP, such as inlet angles and rate of stenosis in the connector. Although sharp inlet angles increased the magnitude of EEP, they cannot be expected to achieve clinically reasonable PEEP. On the other hand, geometrical constriction in the connector can potentially result in PEEP as obtained with conventional nasal cannulae.

Natural Ventilation and Aerosol Particles Dispersion Indoors

Aerosol pollutant particles indoors significantly affect public health. The conventional wisdom is that natural ventilation will alleviate the dispersion of airborne or aerosol particles. However, we show that the problem is far more complex and that natural ventilation should be applied under specific conditions to be effective. We performed several simulations of a simplified (and easily reproducible) room with a window opening and aerosol particles stratified layers. Opening a window can scatter particles present in stratified layers indoors and potentially contribute to the degradation of indoor air quality for a significant period of time. Moreover, we show that thermal instabilities arising from the temperature gradients due to temperature differences between the indoor and outdoor environment spread the particles randomly indoors, adversely affecting air quality and architectural design. Recommendations for more efficient natural ventilation minimizing aerosol pollutant particles dispersed indoors are provided.

Aerosol Transport Modeling: The Key Link Between Lung Infections of Individuals and Populations

The recent COVID-19 pandemic has propelled the field of aerosol science to the forefront, particularly the central role of virus-laden respiratory droplets and aerosols. The pandemic has also highlighted the critical need, and value for, an information bridge between epidemiological models (that inform policymakers to develop public health responses) and within-host models (that inform the public and health care providers how individuals develop respiratory infections). Here, we review existing data and models of generation of respiratory droplets and aerosols, their exhalation and inhalation, and the fate of infectious droplet transport and deposition throughout the respiratory tract. We then articulate how aerosol transport modeling can serve as a bridge between and guide calibration of within-host and epidemiological models, forming a comprehensive tool to formulate and test hypotheses about respiratory tract exposure and infection within and between individuals.

Advancing therapies for viral infections using mechanistic computational models of the dynamic interplay between the virus and host immune response

The COVID-19 pandemic has highlighted a need for improved frameworks for drug discovery, repurposing, clinical trial design and therapy optimization and personalization. Mechanistic computational models can play an important role in developing these frameworks. We discuss how mechanistic models, which consider viral entry, replication in target cells, viral spread in the body, immune response, and the complex factors involved in tissue and organ damage and recovery, can clarify the mechanisms of humoral and cellular immune responses to the virus, viral distribution and replication in tissues, the origins of pathogenesis and patient-to-patient heterogeneity in responses. These models are already improving our understanding of the mechanisms of action of antivirals and immune modulators. We discuss how closer collaboration between the experimentalists, clinicians and modelers could result in more predictive models which may guide therapies for viral infections, improving survival and leading to faster and more complete recovery.

SARS-CoV-2 droplet deposition path and its effects on the human upper airway in the oral inhalation

BACKGROUND AND OBJECTIVE

It is crucial to study the uptake of viral droplets in the human respiratory system to control, prevent, and treat diseases.

METHODS

In this study, a well-verified real anatomical model was used; the passage of air in the human upper respiratory system computed using high-quality Computer Tomography (CT) images. Then, the airflow field, along with the coronavirus micro-droplets injection, was examined in this realistic model using the Fluid-Structure Interaction (FSI) method. The Discrete Phase Model (DPM) was used to solve the field, and with the help of it, the accurate assessment of the temporal and spatial motion of the deposition in the virus-impregnated droplets was obtained in vitro in the upper respiratory system.

RESULTS

The results show that the amount of deposited micro-droplets in the nasal cavity area is meager at the inhalation only through the oral. However, it has the most residence time in this area. The most and least droplet absorption occurred in the oral cavity and larynx-trachea, respectively. Deposition efficiency is about 100% in 30 L/min flow rate and 10 μm diameter; in other words, no droplet enters the lungs. This study's other achievements include the relatively inverse relationship between droplets deposition efficiency in some parts of the upper airway, which have the most deformation in the tract.

CONCLUSIONS

Utilization of a realistic model with accurate and precise computational analysis can end speculation about the deposition zone, accumulation, and the effects of the COVID-19 virus on the upper respiratory tract. On the other hand, recognizing the virus-containing droplet location can ease understanding the areas where the virus can first infect in the upper respiratory tract.

Role of CFD based in silico modelling in establishing an in vitro-in vivo correlation of aerosol deposition in the respiratory tract

Effective evaluation and prediction of aerosol transport deposition in the human respiratory tracts are critical to aerosol drug delivery and evaluation of inhalation products. Establishment of an in vitro-in vivo correlation (IVIVC) requires the understanding of flow and aerosol behaviour and underlying mechanisms at the microscopic scale. The achievement of the aim can be facilitated via computational fluid dynamics (CFD) based in silico modelling which treats the aerosol delivery as a two-phase flow. CFD modelling research, in particular coupling with discrete phase model (DPM) and discrete element method (DEM) approaches, has been rapidly developed in the past two decades. This paper reviews the recent development in this area. The paper covers the following aspects: geometric models of the respiratory tract, CFD turbulence models for gas phase and its coupling with DPM/DEM for aerosols, and CFD investigation of the effects of key factors associated with geometric variations, flow and powder characteristics. The review showed that in silico study based on CFD models can effectively evaluate and predict aerosol deposition pattern in human respiratory tracts. The review concludes with recommendations on future research to improve in silico prediction to achieve better IVIVC.

A dynamical overview of droplets in the transmission of respiratory infectious diseases

The outbreak of the coronavirus disease has drawn public attention to the transmission of infectious pathogens, and as major carriers of those pathogens, respiratory droplets play an important role in the process of transmission. This Review describes respiratory droplets from a physical and mechanical perspective, especially their correlation with the transmission of infectious pathogens. It covers the important aspects of (i) the generation and expulsion of droplets during respiratory activities, (ii) the transport and evolution of respiratory droplets in the ambient environment, and (iii) the inhalation and deposition of droplets in the human respiratory tract. State-of-the-art experimental, computational, and theoretical models and results are presented, and the corresponding knowledge gaps are identified. This Review stresses the multidisciplinary nature of its subject and appeals for collaboration among different fields to fight the present pandemic.

Heterogeneity in lobar and near-acini deposition of inhaled aerosol in the mouse lung

Laboratory animals are often used to derive health risk from environmental exposure or to assess the therapeutic effect of a drug delivered by inhaled therapy. Knowledge of the in-situ distribution of deposited particles on airway and alveolar surfaces is essential in any assessment of these effects. A unique database including both high-resolution lung anatomy and deposition data in four strains of laboratory mice have been recently made publicly available to the research community (https://doi.org/10.25820/9arg-9w56). Using these data, we investigated the effect of particle size on the distribution of deposited particles at the lobar and near-acini level. Analysis was performed on a total of 33 mice where 3, 16 and 14 animals were exposed to 0.5μm, 1μm and 2μm particles, respectively. Ratio of normalized deposition to normalized volume was calculated for each lobe (DV lobe ). At the near-acini level, the skew and standard deviation of the frequency distribution of particle deposition were calculated. Significant deviation above 1 was found for DV ratio in the cranial lobe (DV Cranial ). DV Middle , DV Caudal and DV Accessory were all significantly <1 and lower than DV left (p<0.01). At the near-acini level, skew and standard deviation were positively correlated with particle size and the presence of hot spots (high deposition) were mainly found in the apical region of the lung. These results highlight the uneven distribution of deposited particles in the mouse lung. Thus, depending on the lung sample location, individual analysis to determine overall deposition may either underestimate or overestimate total lung burden, at least for micron-sized particles.

Efficient bi-directional coupling of 3D Computational Fluid-Particle Dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. As a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 μm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFPD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. These methods can be used as the foundation for multi-scale modeling of the full respiratory system.

Innovative preclinical models for pulmonary drug delivery research

Introduction: Pulmonary drug delivery is a complex field of research combining physics which drive aerosol transport and deposition and biology which underpins efficacy and toxicity of inhaled drugs. A myriad of preclinical methods, ranging from in-silico to in-vitro, ex-vivo and in-vivo, can be implemented.Areas covered: The present review covers in-silico mathematical and computational fluid dynamics modelization of aerosol deposition, cascade impactor technology to estimated drug delivery and deposition, advanced in-vitro cell culture methods and associated aerosol exposure, lung-on-chip technology, ex-vivo modeling, in-vivo inhaled drug delivery, lung imaging, and longitudinal pharmacokinetic analysis.Expert opinion: No single preclinical model can be advocated; all methods are fundamentally complementary and should be implemented based on benefits and drawbacks to answer specific scientific questions. The overall best scientific strategy depends, among others, on the product under investigations, inhalation device design, disease of interest, clinical patient population, previous knowledge. Preclinical testing is not to be separated from clinical evaluation, as small proof-of-concept clinical studies or conversely large-scale clinical big data may inform preclinical testing. The extend of expertise required for such translational research is unlikely to be found in one single laboratory calling for the setup of multinational large-scale research consortiums.

In Silico Optimization of Fiber-Shaped Aerosols in Inhalation Therapy for Augmented Targeting and Deposition across the Respiratory Tract

Motivated by a desire to uncover new opportunities for designing the size and shape of fiber-shaped aerosols towards improved pulmonary drug delivery deposition outcomes, we explore the transport and deposition characteristics of fibers under physiologically inspired inhalation conditions in silico, mimicking a dry powder inhaler (DPI) maneuver in adult lung models. Here, using computational fluid dynamics (CFD) simulations, we resolve the transient translational and rotational motion of inhaled micron-sized ellipsoid particles under the influence of aerodynamic (i.e., drag, lift) and gravitational forces in a respiratory tract model spanning the first seven bifurcating generations (i.e., from the mouth to upper airways), coupled to a more distal airway model representing nine generations of the mid-bronchial tree. Aerosol deposition efficiencies are quantified as a function of the equivalent diameter (dp) and geometrical aspect ratio (AR), and these are compared to outcomes with traditional spherical particles of equivalent mass. Our results help elucidate how deposition patterns are intimately coupled to dp and AR, whereby high AR fibers in the narrow range of dp = 6-7 µm yield the highest deposition efficiency for targeting the upper- and mid-bronchi, whereas fibers in the range of dp= 4-6 µm are anticipated to cross through the conducting regions and reach the deeper lung regions. Our efforts underscore previously uncovered opportunities to design the shape and size of fiber-like aerosols towards targeted pulmonary drug delivery with increased deposition efficiencies, in particular by leveraging their large payloads for deep lung deposition.

Numerical modeling of particle deposition in the conducting airways of asthmatic children

Mounting evidence has linked long- and short-term exposure to particulate air pollution with the incidence and exacerbation of asthma in children, but the biological pathogenesis is unclear. We examined the deposition of particles in the airways of asthmatic children. A planar and symmetric model of airways for 4-year-old asthmatic children was considered. Airflow and particle deposition in the upper (G3-G6) and lower (G9-G12) conducting airways were numerically investigated using computation fluid dynamics (CFD) method. We considered the manifestation of moderate (30% reduction in airway diameter) and severe (60% reduction) asthma. Micron particles (1-10 µm) were considered. We found that particle deposition in the asthmatic children was significantly higher than that in healthy children. The deposition efficiency increased slowly with particle size for healthy children, but increased rapidly for asthmatic children, such that smaller particles could be deposited in the conducting airways of asthmatics. For healthy children, particles were deposited by inertial impaction and gravitational sedimentation respectively in the upper and lower airways, but deposited by inertial impaction in asthmatic children. The severity of the asthma increased the particle deposition in the airways. Our study indicated that asthmatic children were more susceptible to the effect of particulate air pollution. The constricted airways increased the particle deposition by inertial impaction, which may be the biological pathogenesis that causes the hospitalization of asthma in children. Avoiding exposure during air pollution events will be an effective measure to reduce the asthma attack.

Assessing Airflow Sensitivity to Healthy and Diseased Lung Conditions in a Computational Fluid Dynamics Model Validated In Vitro

Computational models are useful for understanding respiratory physiology. Crucial to such models are the boundary conditions specifying the flow conditions at truncated airway branches (terminal flow rates). However, most studies make assumptions about these values, which are difficult to obtain in vivo. We developed a computational fluid dynamics (CFD) model of airflows for steady expiration to investigate how terminal flows affect airflow patterns in respiratory airways. First, we measured in vitro airflow patterns in a physical airway model, using particle image velocimetry (PIV). The measured and computed airflow patterns agreed well, validating our CFD model. Next, we used the lobar flow fractions from a healthy or chronic obstructive pulmonary disease (COPD) subject as constraints to derive different terminal flow rates (i.e., three healthy and one COPD) and computed the corresponding airflow patterns in the same geometry. To assess airflow sensitivity to the boundary conditions, we used the correlation coefficient of the shape similarity (R) and the root-mean-square of the velocity magnitude difference (Drms) between two velocity contours. Airflow patterns in the central airways were similar across healthy conditions (minimum R, 0.80) despite variations in terminal flow rates but markedly different for COPD (minimum R, 0.26; maximum Drms, ten times that of healthy cases). In contrast, those in the upper airway were similar for all cases. Our findings quantify how variability in terminal and lobar flows contributes to airflow patterns in respiratory airways. They highlight the importance of using lobar flow fractions to examine physiologically relevant airflow characteristics.

CFD analysis of the flow structure in a monkey upper airway validated by PIV experiments

Inhalation exposure to airborne contaminants has adverse effects on humans; however, related research is typically conducted using in vivo/in vitro tests on animals. Extrapolating the test results is complicated by anatomical and physiological differences between animals and humans and a lack of understanding of the transport mechanism inside their respective respiratory tracts. This study determined the detailed air-flow structure in the upper airway of a monkey. A steady computational fluid dynamics simulation, which was validated by previous particle image velocimetry measurements, was adopted for flow rates of 4 L/min and 10 L/min to analyze the flow structure from the nasal/oral cavities to the trachea region in a monkey airway model. The low Reynolds number type k-ε model provided a reasonably accurate prediction of the airflow in a monkey upper airway. Furthermore, it was confirmed that large velocity gradients were generated in the nasal vestibule and larynx regions, as well as increased turbulent air kinetic energy and wall sheer stress.

Lung Aerosol Dynamics of Airborne Influenza A Virus-Laden Droplets and the Resultant Immune System Responses: An In Silico Study

Influenza A Virus (IAV) replications start from the deposition of inhaled virus-laden droplets on the epithelial cells in the pulmonary tracts. In order to understand the local deposition patterns and within-host dynamics of infectious aerosols, accurate information of high-resolution imaging capabilities, as well as real-time flow cytometry analysis, are required for tracking infected cells, virus agents, and immune system responses. However, clinical and animal studies are in deficit to meet the above-mentioned demands, due to their limited operational flexibility and imaging resolution. Therefore, this study developed an experimentally validated multiscale epidemiological computational model, i.e., the Computational Fluid-Particle Dynamics (CFPD) plus Host Cell Dynamics (HCD) model, to predict the transport and deposition of the low-strain IAV-laden droplets, as well as the resultant regional immune system responses. The hygroscopic growth and shrinkage of IAV-laden droplets were accurately modeled. The subject-specific respiratory system was discretized by generating the new polyhedral-core mesh. By simulating both mouth and nasal breathing scenarios, the inhalations of isotonic IAV-laden droplets with three different compositions were achieved. It is the first time that parametric analysis was performed using the multiscale model on how different exposure conditions can influence the virus aerodynamics in the lung and the subsequent immune system responses. Numerical results show a higher viral accretion followed by a faster immune system response in the supraglottic region when droplets with the higher salt concentration were inhaled. Consequently, more severe symptoms and longer recovery are expected at the pharynx. Furthermore, local deposition maps of IAV-laden droplets and post-deposition infection dynamics provide informative and direct evidence which significantly enhance the fundamental understanding of the underlying mechanisms for upper airway and lower airway infections.

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.

Analysis of the volume of fluid (VOF) method for the simulation of the mucus clearance process with CFD

The clearance of mucus through coughing is a complex, multiphase process, which is affected principally by mucus viscosity and airflow velocity; however, it is also critically affected by the thickness of the two layers of mucus-the serous and gel layers-and oscillation level. The present study examines the effects of the latter parameters more closely. To do so, the mucus clearance process is simulated with a transient 3D volume of fluid (VOF) multiphase model in ANSYS Fluent. The model includes mucus' bilayer properties and a wide range of boundary conditions were tested. The model was analysed in both a straight tube and a realistic trachea. Ultimately, the model was able to both capture air-mucus interface wave evolution and predict the overall behaviour of the clearance process. The results were consistent with experimental clearance data and numerical airflow simulations, which indicates our methodology is appropriate for future studies. Ultimately, the mere presence of the serous layer was found to increase mucus clearance by more than 15 percent. An oscillating flow enhanced clearance by up to 5 percent. Interestingly, interface wave steepness was found to be inversely correlated with mucus thickness, but directly with mucus velocity, which suggests it will be an interesting parameter for further study.

Fluid dynamic assessment of tracheal flow in infants with congenital tracheal stenosis before and after surgery

Tracheal flow in infants with congenital tracheal stenosis (CTS) was numerically investigated using subject-specific airway models before and after reconstructive surgery. We quantified tracheal flow based on airway resistance during inhalation, and compared it between controls and patients before and after surgery. The airway resistance in each subject was assessed using geometrical parameters of the trachea: the minimum cross-sectional area A_min, the minimum cross-sectional area normalized by the standard deviation of the cross-sectional area A_min/σ_A, the area ratio of the minimum and maximum cross-sectional area A_min/A_max, and ratio of the normalized standard deviation of cross-sectional area to the mean cross-sectional area σ_A/A_mean. Our numerical results demonstrated that such geometrical parameters could be used to assess the severity of CTS. Since subjects can be more clearly categorized as controls and most preoperative patients in terms of the airway resistance, a simulation using subject-specific airway models can lead us to a precise understanding of tracheal flow, and also provide knowledge about therapeutic decision. Our numerical results also demonstrated that significant surgical expansion of cross-sectional area did not help recover tracheal flow because of expansion loss. These results will be helpful not only when making therapeutic decisions about surgery but also when assessing quality of life in postoperative patients. Graphical abstract.

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.

A Quasi-3D compartmental multi-scale approach to detect and quantify diseased regional lung constriction using spirometry data

Spirometry is a widely used pulmonary function test to detect the airflow limitations associated with various obstructive lung diseases, such as asthma, chronic obstructive pulmonary disease, and even obesity-related complications. These conditions arise due to the change in the airway resistance, alveolar compliance, and inductance values. Currently, zero-dimensional compartmental models are commonly used for calibrating these resistance, compliance, and inductance values, ie, solving the inverse spirometry problem. However, zero-dimensional compartments cannot capture the flow physics or the spatial geometry effects, thereby generating a low fidelity prediction of the diseased lung. Computational fluid dynamics (CFD) models offer higher fidelity solutions but may be impractical for certain applications due to the duration of these simulations. Recently, a novel, fast-running, and robust Quasi-3D (Q3D) wire model for simulating the airflow in the human lung airway was developed by CFD Research Corporation. This Q3D method preserved the 3D spatial nature of the airways and was favorably validated against CFD solutions. In the present study, the Q3D compartmental multi-scale combination is further improved to predict regional lung constriction of diseased lungs using spirometry data. The Q3D mesh is resolved up to the eighth lung airway generation. The remainder of the airways and the alveoli sections are modeled using a compartmental approach. The Q3D geometry is then split into different spatial sections, and the resistance values in these regions are obtained using parameter inversion. Finally, the airway diameter values are then reduced to create the actual diseased lung model, corresponding to these resistance values. This diseased lung model can be used for patient-specific drug deposition predictions and the subsequent optimization of the orally inhaled drug products.

Airflow Simulations in Infant, Child, and Adult Pulmonary Conducting Airways

The airway structure continuously evolves from birth to adulthood, influencing airflow dynamics and respiratory mechanics. We currently know very little about how airflow patterns change throughout early life and its impact on airway resistance, namely because of experimental limitations. To uncover differences in respiratory dynamics between age groups, we performed subject-specific airflow simulations in an infant, child, and adult conducting airways. Airflow throughout the respiration cycle was calculated by coupling image-based models of the conducting airways to the global respiratory mechanics, where flow was driven by a pressure differential. Trachea diameter was 19, 9, and 4.5 mm for the adult (36 years, female), child (6 years, male), and infant (0.25 years, female), respectively. Mean Reynolds number within the trachea was nearly the same for each subject (1100) and Womersley number was above unity for all three subjects and largest for the adult, highlighting the significance of transient effects. In general, air speeds and airway resistances within the conducting airways were inversely correlated with age; the 3D pressure drop was highest in the infant model. These simulations provide new insight into age-dependent flow dynamics throughout the respiration cycle within subject-specific airways.

Computational fluid dynamics evaluation of excessive dynamic airway collapse

BACKGROUND

Excessive dynamic airway collapse, which is often caused by the collapse of the posterior membrane wall during exhalation, is often misdiagnosed with other diseases; stents can provide support for the collapsing airways. The standard pulmonary function tests do not necessarily show change in functional breathing condition for evaluation of these type of diseases.

METHODS

Flow characteristics through a patient's airways with excessive dynamic airway collapse have been numerically investigated. A stent was placed to support the collapsing airway and to improve breathing conditions. Computed tomography images of the patient's pre- and post-stenting were used for generating 3-Dimensional models of the airways, and were imported into a computational fluid dynamics software for simulation of realistic air flow behavior. Unsteady simulations of the inspiratory phase and expiratory phase were performed with patient-specific boundary conditions for pre- and post-intervention cases to investigate the effect of stent placement on flow characteristic and possible improvements.

FINDINGS

Results of post-stent condition show reduced pressure, velocity magnitude and wall shear stress during expiration. The variation in wall shear stress, velocity magnitude and pressure drop is negligible during inspiration.

INTERPRETATION

Although Spirometry tests do not show significant improvements, computational fluid dynamics results show significant improvements in pre- and post-treatment results, suggesting improvement in breathing condition.

Particulate Deposition in a Patient With Tracheal Stenosis

The presence of obstructions such as tracheal stenosis has important effects on respiratory functions. Tracheal stenosis impacts the therapeutic efficacy of inhaled medications as a result of alterations in particle transport and deposition pattern. This study explores the effects of the presence and absence of stenosis/obstruction in the trachea on air flow characteristics and particle depositions. Computational fluid dynamics (CFD) simulations were performed on three-dimensional (3D) patient-specific models created from computed tomography (CT) images. The analyzed model was generated from a subject with tracheal stenosis and includes the airway tree up to eight generations. CT scans of expiratory and inspiratory phases were used for patient-specific boundary conditions. Pre- and post-intervention CFD simulations' comparison reveals the effect of the stenosis on the characteristics of air flow, transport, and depositions of particles with diameters of 1, 2.5, 4, 6, 8, and 10 μm. Results indicate that the existence of the stenosis inflicts a major pressure force on the flow of inhaled air, leading to an increased deposition of particles both above and below the stenosis. Comparisons of the decrease in pressure in each generation between pre- and post-tracheal stenosis intervention demonstrated a significant reduction in pressure following the stenosis, which was maintained in all downstream generations. Good agreements were found using experimental validation of CFD findings with a model of the control subject up to the third generation, constructed via additive layer manufacturing from CT images.

Morphological and functional properties of the conducting human airways investigated by in vivo computed tomography and in vitro MRI

The accurate representation of the human airway anatomy is crucial for understanding and modeling the structure-function relationship in both healthy and diseased lungs. The present knowledge in this area is based on morphometric studies of excised lung casts, partially complemented by in vivo studies in which computed tomography (CT) was used on a small number of subjects. In the present study, we analyzed CT scans of a cohort of healthy subjects and obtained comprehensive morphometric information down to the seventh generation of bronchial branching, including airway diameter, length, branching angle, and rotation angle. Although some of the geometric parameters (such as the child-to-parent branch diameter ratio) are found to be in line with accepted values, for others (such as the branch length-to-diameter ratio) our findings challenge the common assumptions. We also evaluated several metrics of self-similarity, including the fractal dimension of the airway tree. Additionally, we used phase-contrast magnetic resonance imaging (MRI) to obtain the volumetric flow field in the three-dimensional-printed airway model of one of the subjects during steady inhalation. This is used to relate structural and functional parameters and, in particular, to close the power-law relationship between branch flow rate and diameter. The diameter exponent is found to be significantly lower than in the usually assumed Poiseuille regime, which we attribute to the strong secondary (i.e., transverse) velocity component. The strength of the secondary velocity with respect to the axial component exceeds the levels found in idealized airway models and persists within the first seven generations. NEW & NOTEWORTHY We performed a comprehensive computed tomography-based study of the conductive airway morphology in normal human subjects, including branch diameter, length, and mutual angles. We found significant departure from classic homothetic relationships. We also carried out MRI measurements of the three-dimensional inspiratory flow in an anatomy-based model and directly assessed structure-function relationships that have so far been assumed. We found that strong secondary flows (i.e., transverse velocity components) persist through the first seven generations of bronchial branching.

Regional aerosol deposition in the human airways: The SimInhale benchmark case and a critical assessment of in silico methods

Regional deposition effects are important in the pulmonary delivery of drugs intended for the topical treatment of respiratory ailments. They also play a critical role in the systemic delivery of drugs with limited lung bioavailability. In recent years, significant improvements in the quality of pulmonary imaging have taken place, however the resolution of current imaging modalities remains inadequate for quantifying regional deposition. Computational Fluid-Particle Dynamics (CFPD) can fill this gap by providing detailed information about regional deposition in the extrathoracic and conducting airways. It is therefore not surprising that the last 15years have seen an exponential growth in the application of CFPD methods in this area. Survey of the recent literature however, reveals a wide variability in the range of modelling approaches used and in the assumptions made about important physical processes taking place during aerosol inhalation. The purpose of this work is to provide a concise critical review of the computational approaches used to date, and to present a benchmark case for validation of future studies in the upper airways. In the spirit of providing the wider community with a reference for quality assurance of CFPD studies, in vitro deposition measurements have been conducted in a human-based model of the upper airways, and several groups within MP1404 SimInhale have computed the same case using a variety of simulation and discretization approaches. Here, we report the results of this collaborative effort and provide a critical discussion of the performance of the various simulation methods. The benchmark case, in vitro deposition data and in silico results will be published online and made available to the wider community. Particle image velocimetry measurements of the flow, as well as additional numerical results from the community, will be appended to the online database as they become available in the future.

A quasi-3D wire approach to model pulmonary airflow in human airways

The models used for modeling the airflow in the human airways are either 0-dimensional compartmental or full 3-dimensional (3D) computational fluid dynamics (CFD) models. In the former, airways are treated as compartments, and the computations are performed with several assumptions, thereby generating a low-fidelity solution. The CFD method displays extremely high fidelity since the solution is obtained by solving the conservation equations in a physiologically consistent geometry. However, CFD models (1) require millions of degrees of freedom to accurately describe the geometry and to reduce the discretization errors, (2) have convergence problems, and (3) require several days to simulate a few breathing cycles. In this paper, we present a novel, fast-running, and robust quasi-3D wire model for modeling the airflow in the human lung airway. The wire mesh is obtained by contracting the high-fidelity lung airway surface mesh to a system of connected wires, with well-defined radii. The conservation equations are then solved in each wire. These wire meshes have around O(1000) degrees of freedom and hence are 3000 to 25 000 times faster than their CFD counterparts. The 3D spatial nature is also preserved since these wires are contracted out of the actual lung STL surface. The pressure readings between the 2 approaches showed minor difference (maximum error = 15%). In general, this formulation is fast and robust, allows geometric changes, and delivers high-fidelity solutions. Hence, this approach has great potential for more complicated problems including modeling of constricted/diseased lung sections and for calibrating the lung flow resistances through parameter inversion.

Validating Whole-Airway CFD Predictions of DPI Aerosol Deposition at Multiple Flow Rates

BACKGROUND

The objective of this study was to compare aerosol deposition predictions of a new whole-airway CFD model with available in vivo data for a dry powder inhaler (DPI) considered across multiple inhalation waveforms, which affect both the particle size distribution (PSD) and particle deposition.

METHODS

The Novolizer DPI with a budesonide formulation was selected based on the availability of 2D gamma scintigraphy data in humans for three different well-defined inhalation waveforms. Initial in vitro cascade impaction experiments were conducted at multiple constant (square-wave) particle sizing flow rates to characterize PSDs. The whole-airway CFD modeling approach implemented the experimentally determined PSDs at the point of aerosol formation in the inhaler. Complete characteristic airway geometries for an adult were evaluated through the lobar bronchi, followed by stochastic individual pathway (SIP) approximations through the tracheobronchial region and new acinar moving wall models of the alveolar region.

RESULTS

It was determined that the PSD used for each inhalation waveform should be based on a constant particle sizing flow rate equal to the average of the inhalation waveform's peak inspiratory flow rate (PIFR) and mean flow rate [i.e., AVG(PIFR, Mean)]. Using this technique, agreement with the in vivo data was acceptable with <15% relative differences averaged across the three regions considered for all inhalation waveforms. Defining a peripheral to central deposition ratio (P/C) based on alveolar and tracheobronchial compartments, respectively, large flow-rate-dependent differences were observed, which were not evident in the original 2D in vivo data.

CONCLUSIONS

The agreement between the CFD predictions and in vivo data was dependent on accurate initial estimates of the PSD, emphasizing the need for a combination in vitro-in silico approach. Furthermore, use of the AVG(PIFR, Mean) value was identified as a potentially useful method for characterizing a DPI aerosol at a constant flow rate.

Lactation in the Human Breast From a Fluid Dynamics Point of View

This study is a collaborative effort among lactation specialists and fluid dynamic engineers. The paper presents clinical results for suckling pressure pattern in lactating human breast as well as a 3D computational fluid dynamics (CFD) modeling of milk flow using these clinical inputs. The investigation starts with a careful, statistically representative measurement of suckling vacuum pressure, milk flow rate, and milk intake in a group of infants. The results from clinical data show that suckling action does not occur with constant suckling rate but changes in a rhythmic manner for infants. These pressure profiles are then used as the boundary condition for the CFD study using commercial ansys fluent software. For the geometric model of the ductal system of the human breast, this work takes advantage of a recent advance in the development of a validated phantom that has been produced as a ground truth for the imaging applications for the breast. The geometric model is introduced into CFD simulations with the aforementioned boundary conditions. The results for milk intake from the CFD simulation and clinical data were compared and cross validated. Also, the variation of milk intake versus suckling pressure are presented and analyzed. Both the clinical and CFD simulation show that the maximum milk flow rate is not related to the largest vacuum pressure or longest feeding duration indicating other factors influence the milk intake by infants.

Effect of static vs. dynamic imaging on particle transport in CT-based numerical models of human central airways

Advances in quantitative computed tomography (CT) has provided methods to assess the detailed structure of the pulmonary airways and parenchyma, providing the means of applying computational fluid dynamics-based modeling to better understand subject-specific differences in structure-to-function relationships. Most of the previous numerical studies, seeking to predict patterns of inhaled particle deposition, have considered airway geometry and regional ventilation derived from static images. Because geometric alterations of the airway and parenchyma associated with regional ventilation may greatly affect particle transport, we have sought to investigate the effect of rigid vs. deforming airways, linear vs. nonlinear airway deformations, and step-wise static vs. dynamic imaging on particle deposition with varying numbers of intermediate lung volume increments. Airway geometry and regional ventilation at different time points were defined by four-dimensional (space and time) dynamic or static CT images. Laminar, transitional, and turbulent air flows were reproduced with a three-dimensional eddy-resolving computational fluid dynamics model. Finally, trajectories of particles were computed with the Lagrangian tracking algorithm. The results demonstrated that static-imaging-based models can contribute 7% uncertainty to overall particle distribution and deposition primarily due to regional flow rate (ventilation) differences as opposed to geometric alterations. The effect of rigid vs. deforming airways on serial distribution of particles over generations was significantly smaller than reported in a previous study that used the symmetric Weibel geometric model with smaller flow rate. Rigid vs. deforming airways were also shown to affect parallel particle distribution over lobes by 8% and the differences associated with use of static vs. dynamic imaging was 18%. These differences demonstrate that estimates derived from static vs. dynamic imaging can significantly affect the assessment of particle distribution heterogeneity. The effect of linear vs. nonlinear airway deformations was within the uncertainty due to mesh size.

INHALED AEROSOL DOSIMETRY: SOME CURRENT RESEARCH NEEDS

After the presentation of 60 papers at the conference "Advancing Aerosol Dosimetry Research" (October 24-25, 2014 in Irvine, CA, USA), attendees submitted written descriptions of needed research. About 40 research needs were submitted. The suggestions fell into six broad categories: 1) Access to detailed anatomic data; 2) Access to subject-specific aerosol deposition datasets; 3) Improving current inhaled aerosol deposition models; 4) Some current experimental data needs and hot topics; 5) Linking exposure and deposition modeling to health endpoints; and 6) Developing guidelines for appropriate validation of dosimetry and risk assessment models. Summaries of suggestions are provided here as an update on research needs related to inhaled aerosol dosimetry modeling. Taken together, the recommendations support the overarching need for increased collaborations between dose modelers and those that use the models for risk assessments, aerosol medicine applications, design of toxicology experiments, and extrapolation across species. This paper is only a snapshot in time of perceived research needs from the conference attendees; it does not carry the approval of any agency or other group that plans research priorities or that funds research.

Total and regional deposition of inhaled aerosols in supine healthy subjects and subjects with mild-to-moderate COPD

Despite substantial development of sophisticated subject-specific computational models of aerosol transport and deposition in human lungs, experimental validation of predictions from these new models is sparse. We collected aerosol retention and exhalation profiles in seven healthy volunteers and six subjects with mild-to-moderate COPD (FEV1 = 50-80%predicted) in the supine posture. Total deposition was measured during continuous breathing of 1 and 2.9 μm-diameter particles (tidal volume of 1 L, flow rate of 0.3 L/s and 0.75 L/s). Bolus inhalations of 1 μm particles were performed to penetration volumes of 200, 500 and 800 mL (flow rate of 0.5 L/s). Aerosol bolus dispersion (H), deposition, and mode shift (MS) were calculated from these data. There was no significant difference in total deposition between healthy subjects and those with COPD. Total deposition increased with increasing particle size and also with increasing flow rate. Similarly, there was no significant difference in aerosol bolus deposition between subject groups. Yet, the rate of increase in dispersion and of decrease in MS with increasing penetration volume was higher in subjects with COPD than in healthy volunteers (H: 0.798 ± 0.205 vs. 0.527 ± 0.122 mL/mL, p=0.01; MS: -0.271±0.129 vs. -0.145 ± 0.076 mL/mL, p=0.05) indicating larger ventilation inhomogeneities (based on H) and increased flow sequencing (based on MS) in the COPD than in the healthy group. In conclusion, in the supine posture, deposition appears to lack sensitivity for assessing the effect of lung morphology and/or ventilation distribution alteration induced by mild-to-moderate lung disease on the fate of inhaled aerosols. However, other parameters such as aerosol bolus dispersion and mode shift may be more sensitive parameters for evaluating models of lungs with moderate disease.

Developmental model of an automatic production of the human bronchial tree based on L-system

BACKGROUND AND OBJECTIVE

The human lungs exchange air with the external environment via the conducting airways. The application of an anatomically accurate model of the conducting airways can be helpful for simulating gas exchange and fluid distribution throughout the bronchial tree in the lung.

METHODS

In the current study, Lindenmayer system (L-system) has been formulated to generate the bronchial tree structure in a human lung. It has been considered that the structure of the bronchial tree is divided into two main segments: 1) The central airways (from the trachea to segmental bronchi) and 2) the dichotomous structure (from segmental bronchi to terminal bronchioles). Two sets of parametric rewriting rules which can be used to develop central and peripheral airways have been proposed; the first set used to develop central airways consists of seven rules, while the second rule set contains four rules.

RESULTS

The proposed model is capable of generating bronchial tree inside the volume of the host lung; and comparison of the resulting model with those reported in the literature shows that the morphometric characteristics of L-system structure are in good agreement with their corresponding experimental data.

CONCLUSION

The resulting model can be used to obtain a mathematical model required for the study of transport phenomena occurring in the lung during respiration.

Particle transport in the human respiratory tract: formulation of a nodal inverse distance weighted Eulerian-Lagrangian transport and implementation of the Wind-Kessel algorithm for an oral delivery

This paper is the first in a series wherein efficient computational methods are developed and implemented to accurately quantify the transport, deposition, and clearance of the microsized particles (range of interest: 2 to 10 µm) in the human respiratory tract. In particular, this paper (part I) deals with (i) development of a detailed 3D computational finite volume mesh comprising of the NOPL (nasal, oral, pharyngeal and larynx), trachea and several airway generations; (ii) use of CFD Research Corporation's finite volume Computational Biology (CoBi) flow solver to obtain the flow physics for an oral inhalation simulation; (iii) implement a novel and accurate nodal inverse distance weighted Eulerian-Lagrangian formulation to accurately obtain the deposition, and (iv) development of Wind-Kessel boundary condition algorithm. This new Wind-Kessel boundary condition algorithm allows the 'escaped' particles to reenter the airway through the outlets, thereby to an extent accounting for the drawbacks of having a finite number of lung generations in the computational mesh. The deposition rates in the NOPL, trachea, the first and second bifurcation were computed, and they were in reasonable accord with the Typical Path Length model. The quantitatively validated results indicate that these developments will be useful for (i) obtaining depositions in diseased lungs (because of asthma and COPD), for which there are no empirical models, and (ii) obtaining the secondary clearance (mucociliary clearance) of the deposited particles. Copyright © 2015 John Wiley & Sons, Ltd.

Bridging the Gap Between Science and Clinical Efficacy: Physiology, Imaging, and Modeling of Aerosols in the Lung

Development of a new drug for the treatment of lung disease is a complex and time consuming process involving numerous disciplines of basic and applied sciences. During the 2015 Congress of the International Society for Aerosols in Medicine, a group of experts including aerosol scientists, physiologists, modelers, imagers, and clinicians participated in a workshop aiming at bridging the gap between basic research and clinical efficacy of inhaled drugs. This publication summarizes the current consensus on the topic. It begins with a short description of basic concepts of aerosol transport and a discussion on targeting strategies of inhaled aerosols to the lungs. It is followed by a description of both computational and biological lung models, and the use of imaging techniques to determine aerosol deposition distribution (ADD) in the lung. Finally, the importance of ADD to clinical efficacy is discussed. Several gaps were identified between basic science and clinical efficacy. One gap between scientific research aimed at predicting, controlling, and measuring ADD and the clinical use of inhaled aerosols is the considerable challenge of obtaining, in a single study, accurate information describing the optimal lung regions to be targeted, the effectiveness of targeting determined from ADD, and some measure of the drug's effectiveness. Other identified gaps were the language and methodology barriers that exist among disciplines, along with the significant regulatory hurdles that need to be overcome for novel drugs and/or therapies to reach the marketplace and benefit the patient. Despite these gaps, much progress has been made in recent years to improve clinical efficacy of inhaled drugs. Also, the recent efforts by many funding agencies and industry to support multidisciplinary networks including basic science researchers, R&D scientists, and clinicians will go a long way to further reduce the gap between science and clinical efficacy.

Distribution of aerosolized particles in healthy and emphysematous rat lungs: comparison between experimental and numerical studies

In silico models of airflow and particle deposition in the lungs are increasingly used to determine the therapeutic or toxic effects of inhaled aerosols. While computational methods have advanced significantly, relatively few studies have directly compared model predictions to experimental data. Furthermore, few prior studies have examined the influence of emphysema on particle deposition. In this work we performed airflow and particle simulations to compare numerical predictions to data from our previous aerosol exposure experiments. Employing an image-based 3D rat airway geometry, we first compared steady flow simulations to coupled 3D-0D unsteady simulations in the healthy rat lung. Then, in 3D-0D simulations, the influence of emphysema was investigated by matching disease location to the experimental study. In both the healthy unsteady and steady simulations, good agreement was found between numerical predictions of aerosol delivery and experimental deposition data. However, deposition patterns in the 3D geometry differed between the unsteady and steady cases. On the contrary, satisfactory agreement was not found between the numerical predictions and experimental data for the emphysematous lungs. This indicates that the deposition rate downstream of the 3D geometry is likely proportional to airflow delivery in the healthy lungs, but not in the emphysematous lungs. Including small airway collapse, variations in downstream airway size and tissue properties, and tracking particles throughout expiration may result in a more favorable agreement in future studies.

Nanoparticle mass transfer from lung airways to systemic regions--Part I: Whole-lung aerosol dynamics

This is a two-part paper describing inhaled nanoparticle (NP) transport and deposition in a model of a human respiratory tract (Part I) as well as NP-mass transfer across barriers into systemic regions (Part II). Specifically, combining high-resolution computer simulation results of inhaled NP deposition in the human airways (Part I) with a multicompartmental model for NP-mass transfer (Part II) allows for the prediction of temporal NP accumulation in the blood and lymphatic systems as well as in organs. An understanding of nanoparticle transport and deposition in human respiratory airways is of great importance, as exposure to nanomaterial has been found to cause serious lung diseases, while the use of nanodrugs may have superior therapeutic effects. In Part I, the fluid-particle dynamics of a dilute NP suspension was simulated for the entire respiratory tract, assuming steady inhalation and planar airways. Thus, a realistic airway configuration was considered from nose/mouth to generation 3, and then an idealized triple-bifurcation unit was repeated in series and parallel to cover the remaining generations. Using the current model, the deposition of NPs in distinct regions of the lung, namely extrathoracic, bronchial, bronchiolar, and alveolar, was calculated. The region-specific NP-deposition results for the human lung model were used in Part II to determine the multicompartmental model parameters from experimental retention and clearance data in human lungs. The quantitative, experimentally validated results are useful in diverse fields, such as toxicology for exposure-risk analysis of ubiquitous nanomaterial as well as in pharmacology for nanodrug development and targeting.

Development of characteristic upper tracheobronchial airway models for testing pharmaceutical aerosol delivery

Characteristic models of the upper conducting airways are needed to evaluate the performance of existing pharmaceutical inhalers and to develop new respiratory drug delivery strategies. Previous studies have focused on the development of characteristic mouth-throat (MT) geometries for orally inhaled products; however, characteristic upper tracheobronchial (TB) geometries are currently not available. In this study, a new characteristic model of the upper TB airways for an average adult male was developed based on an analysis of new and existing anatomical data. Validated computational fluid dynamics (CFD) simulations were used to evaluate the deposition of monodisperse and realistic polydisperse aerosols from multiple inhalers. Comparisons of deposition results between the new model and a simpler geometry were used to identify the effects of different anatomical features on aerosol deposition. The CFD simulations demonstrated a good match to regional pharmaceutical aerosol deposition from in vitro experiments in the same geometry. The deposition of both monodisperse and pharmaceutical aerosols was increased in the new TB geometry as a result of additional anatomical detail on a regional and highly localized basis. Tracheal features including an accurate coronal angle, asymmetry, and curvature produced a skewed laryngeal jet and significantly increased regional deposition. Branch curvature and realistic cross-sections increased deposition in the remainder of the TB model. A hexahedral mesh style was utilized to provide the best solution. In conclusion, a number of physiological features in the upper TB region were shown to influence deposition and should be included in a characteristic model of respiratory drug delivery.

Early airway structural changes in cystic fibrosis pigs as a determinant of particle distribution and deposition

The pathogenesis of cystic fibrosis (CF) airway disease is not well understood. A porcine CF model was recently generated, and these animals develop lung disease similar to humans with CF. At birth, before infection and inflammation, CF pigs have airways that are irregularly shaped and have a reduced caliber compared to non-CF pigs. We hypothesized that these airway structural abnormalities affect airflow patterns and particle distribution. To test this hypothesis we used computational fluid dynamics (CFD) on airway geometries obtained by computed tomography of newborn non-CF and CF pigs. For the same flow rate, newborn CF pig airways exhibited higher air velocity and resistance compared to non-CF. Moreover we found that, at the carina bifurcation, particles greater than 5-μm preferably distributed to the right CF lung despite almost equal airflow ventilation in non-CF and CF. CFD modeling also predicted that deposition efficiency was greater in CF compared to non-CF for 5- and 10-μm particles. These differences were most significant in the airways included in the geometry supplying the right caudal, right accessory, left caudal, and left cranial lobes. The irregular particle distribution and increased deposition in newborn CF pig airways suggest that early airway structural abnormalities might contribute to CF disease pathogenesis.