Aerosol bolus dispersion in acinar airways--influence of gravity and airway asymmetry

Exploring the role of electrostatic deposition on inhaled aerosols in alveolated microchannels

Large amounts of net electrical charge are known to accumulate on inhaled aerosols during their generation using commonly-available inhalers. This effect often leads to superfluous deposition in the extra-thoracic airways at the cost of more efficient inhalation therapy. Since the electrostatic force is inversely proportional to the square of the distance between an aerosol and the airway wall, its role has long been recognized as potentially significant in the deep lungs. Yet, with the complexity of exploring such phenomenon directly at the acinar scales, in vitro experiments have been largely limited to upper airways models. Here, we devise a microfluidic alveolated airway channel coated with conductive material to quantify in vitro the significance of electrostatic effects on inhaled aerosol deposition. Specifically, our aerosol exposure assays showcase inhaled spherical particles of 0.2, 0.5, and 1.1 μm that are recognized to reach the acinar regions, whereby deposition is typically attributed to the leading roles of diffusion and sedimentation. In our experiments, electrostatic effects are observed to largely prevent aerosols from depositing inside alveolar cavities. Rather, deposition is overwhelmingly biased along the inter-alveolar septal spaces, even when aerosols are charged with only a few elementary charges. Our observations give new insight into the role of electrostatics at the acinar scales and emphasize how charged particles under 2 µm may rapidly overshadow the traditionally accepted dominance of diffusion or sedimentation when considering aerosol deposition phenomena in the deep lungs.

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.

Fluid-structure interaction analysis of airflow, structural mechanics and aerosol dynamics in a four-generation acinar model

Elucidating the aerosol dynamics in the pulmonary acinar region is imperative for both health risk assessment and inhalation therapy, especially nowadays with the occurrence of the global COVID-19 pandemic. During respiration, the chest's outward elastic recoil and the lungs' inward elastic recoil lead to a change of transmural pressure, which drives the lungs to expand and contract to inhale and expel airflow and aerosol. In contrast to research using predefined wall motion, we developed a four-generation acinar model and applied an oscillatory pressure on the model outface to generate structure deformation and airflow. With such tools at hand, we performed a computational simulation that addressed both the airflow characteristic, structural mechanics, and aerosol dynamics in the human pulmonary acinar region. Our results showed that there is no recirculating flow in the sac. The structural displacement and stress were found to be positively related to the change of model volume and peaked at the end of inspiration. It was noteworthy that the stress distribution on the acinar wall was significantly heterogeneous, and obvious concentrations of stress were found at the junction of the alveoli and the ducts or the junction of the alveoli and alveoli in the sac. Our result demonstrated the effect of breathing cycles and aerosol diameter on deposition fraction and location of aerosols in the size range of 0.1-5 μm. Multiple respiratory cycles were found necessary for adequate deposition or escape of submicron particles while having a negligible influence on the transport of large particles, which were dominated by gravity. Our study can provide new insights into the further investigation of airflow, structural mechanics, and aerosol dynamics in the acinar depth.

Influence of alveolar mixing and multiple breaths of aerosol intake on particle deposition in the human lungs

Predictive dosimetry models play an important role in assessing health effect of inhaled particulate matter and in optimizing delivery of inhaled pharmaceutical aerosols. In this study, the commonly used 1D Multiple-Path Particle Dosimetry model (MPPD) was improved by including a mechanistically based model component for alveolar mixing of particles and by extending the model capabilities to account for multiple breaths of aerosol intake. These modifications increased the retained fraction of particles and consequently particle deposition predictions in the deep lung during tidal breathing. Comparison with an existing dataset (J. Aerosol Sci., 99:27-39, 2016) obtained under two breathing conditions referred to as slow and fast breathing showed significant differences in 1 μm particle deposition between predictions based on subject-specific breathing patterns and lung volume (slow: 30 ± 1%, fast: 21 ± 1%, (average ± standard deviation), N = 7) and measurements (slow: 43 ± 9%, fast: 30 ± 5%) when the prior version of MPPD (single breath and no mixing, J. Aerosol Sci., 151:105647, 2021) was used. Adding a mixing model and multiple breaths moved the predictions (slow: 34 ± 2%, fast:25 ± 2%) closer to the range of deposition measurements. For 2.9 μm particles, predictions from both the original (slow: 70 ± 2%, fast: 57 ± 2%) and the revised MPPD model (slow: 71 ± 2%, fast: 59 ± 3%) compared well with experiments (slow: 67 ± 8%, fast: 58 ± 10%). This was expected as suspended fraction of 2.9 μm particles was small and thus the addition of alveolar mixing and multi breath capability only slightly increased the retained fraction for particles of this size and greater. The revised 1D model improves dose predictions in the deep lung and support human risk assessment from exposure to airborne particles.

Recent advances in the understanding of alveolar flow

Understanding the dynamics of airflow in alveoli and its effect on the behavior of particle transport and deposition is important for understanding lung functions and the cause of many lung diseases. The studies on these areas have drawn substantial attention over the last few decades. This Review discusses the recent progress in the investigation of behavior of airflow in alveoli. The information obtained from studies on the structure of the lung airway tree and alveolar topology is provided first. The current research progress on the modeling of alveoli is then reviewed. The alveolar cell parameters at different generation of branches, issues to model real alveolar flow, and the current numerical and experimental approaches are discussed. The findings on flow behavior, in particular, flow patterns and the mechanism of chaotic flow generation in the alveoli are reviewed next. The different flow patterns under different geometrical and flow conditions are discussed. Finally, developments on microfluidic devices such as lung-on-a-chip devices are reviewed. The issues of current devices are discussed.

[Numerical simulation on the deposition characteristics of inhaled particles in human pulmonary acinus region under the influence of multi-factors]

Research on the deposition of inhaled particles in human pulmonary acinus region is important to the pathogenesis investigation, prevention and treatment of lung diseases. Most of the current research focus on the final deposition fraction of inhaled particles in human acinar region, but little is involved in their dynamic deposition characteristics. In this paper, five multi-alveolar models, G3-G7, were built. The evaluation parameter 1/4 deposition time was introduced to study the particle deposition speed. The deposition characteristics of particles in the diameter ranging 0.1-5 μm were numerically simulated and summarized under the influence of factors such as the generation and structure of model, particle diameter and respiratory mode, shedding some new light on the further research of transport of inhaled particles. The results showed that the generation and structure of model had a significance effect on the deposition of particles. 0.1 μm particles were dominated by Brownian diffusion, which experienced a high deposition fraction, a fast deposition speed and a logarithmic deposition curve, while 5 μm particles were dominated by gravitational sedimentation, with a high deposition fraction, a fast deposition speed and an S-shaped deposition curve. The deposition of 0.3-1 μm particles were influenced greatly by convention and varied with the change of respiratory mode. The research methods and results in this paper can provide theoretical basis and data support for the further exploration of the mechanism, prevention and treatment of lung diseases.

Alveolar size effects on nanoparticle deposition in rhythmically expanding-contracting terminal alveolar models

Significant differences in alveolar size exist in humans of different ages, gender, health, and among different species. The effects of alveolar sizes, as well as the accompanying breathing frequencies, on regional and local dosimetry of inhaled nanoparticles have not been sufficiently studied. Despite a well-accepted qualitative understanding of the advection-diffusion-sedimentation mechanism in the acinar region, a quantitative picture of the interactions among these factors remains inchoate. The objective of this study is to quantify the effects of alveolar size on the regional and local deposition of inhaled nanoparticles in alveolar models of varying complexities and to understand the dynamic interactions among different deposition mechanisms. Three different models were considered that retained 1, 4, and 45 alveoli, respectively. For each model, the baseline geometry was scaled by ¼, ½, 2, 4, and 8 times by volume. Temporal evolution and spatial distribution of particle deposition were tracked using a discrete-phase Lagrangian model. Lower retentions of inhaled nanoparticles were observed in the larger alveoli under the same respiration frequency, while similar retentions were found among different geometrical scales if breathing frequencies allometrically matched the alveolar size. Dimensional analysis reveals a manifold deposition mechanism with tantamount contributions from advection, diffusion, and gravitational sedimentation, each of which can become dominant depending on the location in the alveoli. Results of this study indicate that empirical correlations obtained from one sub-population cannot be directly applied to others, nor can they be simply scaled as a function of the alveolar size or respiration frequency due to the regime-transiting deposition mechanism that is both localized and dynamic.

Airflow and Particle Deposition in Acinar Models with Interalveolar Septal Walls and Different Alveolar Numbers

Unique features exist in acinar units such as multiple alveoli, interalveolar septal walls, and pores of Kohn. However, the effects of such features on airflow and particle deposition remain not well quantified due to their structural complexity. This study aims to numerically investigate particle dynamics in acinar models with interalveolar septal walls and pores of Kohn. A simplified 4-alveoli model with well-defined geometries and a physiologically realistic 45-alveoli model was developed. A well-validated Lagrangian tracking model was used to simulate particle trajectories in the acinar models with rhythmically expanding and contracting wall motions. Both spatial and temporal dosimetries in the acinar models were analyzed. Results show that collateral ventilation exists among alveoli due to pressure imbalance. The size of interalveolar septal aperture significantly alters the spatial deposition pattern, while it has an insignificant effect on the total deposition rate. Surprisingly, the deposition rate in the 45-alveoli model is lower than that in the 4-alveoli model, indicating a stronger particle dispersion in more complex models. The gravity orientation angle has a decreasing effect on acinar deposition rates with an increasing number of alveoli retained in the model; such an effect is nearly negligible in the 45-alveoli model. Breath-holding increased particle deposition in the acinar region, which was most significant in the alveoli proximal to the duct. Increasing inhalation depth only slightly increases the fraction of deposited particles over particles entering the alveolar model but has a large influence on dispensing particles to the peripheral alveoli. Results of this study indicate that an empirical correlation for acinar deposition can be developed based on alveolar models with reduced complexity; however, what level of geometry complexity would be sufficient is yet to be determined.

Deposition of bolus and continuously inhaled aerosols in rhythmically moving terminal alveoli

The particle dynamics in an oscillating alveolus under tidal breathing can be dramatically different from those in a static alveolus. Despite its close relevance to pulmonary drug delivery and health risk from airborne exposure, quantifications of alveolar deposition are scarce due to its inaccessibility to in vivo measurement instruments, tiny size to replicate in vitro, and dynamic wall motions to model. The objective of this study is to introduce a numerical method to quantify alveolar deposition with continuous particle release in a rhythmically oscillating alveolus by integrating the deposition curves for bolus aerosols and use this method to develop correlations applicable in assessing alveolar drug delivery efficiency or dosimetry of inhaled toxicants. An idealized blind-end terminal alveolus model was developed with rhythmically moving alveolar boundary conditions in phase with tidal breathing. The dynamic wall expansion mode and magnitude were based on experimentally measured chest wall motions and tidal volumes. A well-validated Lagrangian tracking model was used to simulate the transport and deposition of inhaled micrometer particles. Large differences were observed between dynamic and static alveoli in particle motion, deposition onset, and final alveolar deposition fraction. Alveolar deposition of bolus aerosols is highly sensitive to breath-holding duration, particle release time, and alveolar dimension. For 1 µm particles, there exists a cut-off release time (zero bolus deposition), which decreases with alveolar size (i.e., 1.0 s in a 0.2-mm-diameter alveolus and 0.56 s in a 0.8-mm-diameter alveolus). The cumulative alveolar deposition was predicted to be 39% for a 0.2-mm-diameter alveolus, 22% for a 0.4-mm-diameter alveolus, and 10% for a 0.8-mm-diameter alveolus. A cumulative alveolar deposition correlation was developed for inhalation delivery with a prescribed period of drug release and the second correlation for the time variation of alveolar deposition of ambient aerosols, both of which captured the relative dependence of the particle release time and alveolar dimension.

Capturing complexity in pulmonary system modelling

Respiratory disease is a significant problem worldwide, and it is a problem with increasing prevalence. Pathology in the upper airways and lung is very difficult to diagnose and treat, as response to disease is often heterogeneous across patients. Computational models have long been used to help understand respiratory function, and these models have evolved alongside increases in the resolution of medical imaging and increased capability of functional imaging, advances in biological knowledge, mathematical techniques and computational power. The benefits of increasingly complex and realistic geometric and biophysical models of the respiratory system are that they are able to capture heterogeneity in patient response to disease and predict emergent function across spatial scales from the delicate alveolar structures to the whole organ level. However, with increasing complexity, models become harder to solve and in some cases harder to validate, which can reduce their impact clinically. Here, we review the evolution of complexity in computational models of the respiratory system, including successes in translation of models into the clinical arena. We also highlight major challenges in modelling the respiratory system, while making use of the evolving functional data that are available for model parameterisation and testing.

One (sub-)acinus for all: Fate of inhaled aerosols in heterogeneous pulmonary acinar structures

Computational Fluid Dynamics (CFD) have offered an attractive gateway to investigate in silico respiratory flows and aerosol transport in the depths of the lungs. Yet, not only do existing models lack sufficient anatomical realism in capturing the heterogeneity and morphometry of the acinar environment, numerical simulations have been widely restricted to domains capturing a mere few percent of a single acinus. Here, we present to the best of our knowledge the most detailed and comprehensive in silico simulations to date on the fate of aerosols in the acinar depths. Our heterogeneous acinar domains represent complete sub-acinar models (i.e. 1/8th of a full acinus) based on the recent algorithm of Koshiyama & Wada (2015), capturing statistics of human acinar morphometry (Ochs et al. 2004). Our simulations deliver high-resolution, 3D spatial-temporal data on aerosol transport and deposition, emphasizing how variances in acinar heterogeneity only play a minor role in determining general deposition outcomes. With such tools at hand, we revisit whole-lung deposition predictions (i.e. ICRP) based on past 1D lung models. While our findings under quiet breathing substantiate general deposition trends obtained with past predictions in the alveolar regions, we underscore how deposition fractions are anticipated to increase, in particular during deep inhalation. For such inhalation maneuver, our simulations support the notion of significantly augmented deposition for all aerosol sizes (0.005-5.0μm). Overall, our efforts not only help consolidate our mechanistic understanding of inhaled aerosol transport in the acinar depths but also continue to bridge the gap between "bottom-up" in silico models and regional deposition predictions from whole-lung models. Such quantifications provide what is deemed more accurate deposition predictions in morphometrically-faithful models and are particularly useful in assessing inhalation strategies for deep airway deposition (e.g. systemic delivery).

Image-based computational fluid dynamics in the lung: virtual reality or new clinical practice?

The development and implementation of personalized medicine is paramount to improving the efficiency and efficacy of patient care. In the respiratory system, function is largely dictated by the choreographed movement of air and blood to the gas exchange surface. The passage of air begins in the upper airways, either via the mouth or nose, and terminates at the alveolar interface, while blood flows from the heart to the alveoli and back again. Computational fluid dynamics (CFD) is a well-established tool for predicting fluid flows and pressure distributions within complex systems. Traditionally CFD has been used to aid in the effective or improved design of a system or device; however, it has become increasingly exploited in biological and medical-based applications further broadening the scope of this computational technique. In this review, we discuss the advancement in application of CFD to the respiratory system and the contributions CFD is currently making toward improving precision medicine. The key areas CFD has been applied to in the pulmonary system are in predicting fluid transport and aerosol distribution within the airways. Here we focus our discussion on fluid flows and in particular on image-based clinically focused CFD in the ventilatory system. We discuss studies spanning from the paranasal sinuses through the conducting airways down to the level of the alveolar airways. The combination of imaging and CFD is enabling improved device design in aerosol transport, improved biomarkers of lung function in clinical trials, and improved predictions and assessment of surgical interventions in the nasal sinuses. WIREs Syst Biol Med 2017, 9:e1392. doi: 10.1002/wsbm.1392 For further resources related to this article, please visit the WIREs website.

Streamline crossing: An essential mechanism for aerosol dispersion in the pulmonary acinus

The dispersion of inhaled microparticles in the pulmonary acinus of the lungs is often attributed to the complex interplay between convective mixing, due to irreversible flows, and intrinsic particle motion (i.e. gravity and diffusion). However, the role of each mechanism, the exact nature of such interplay between them and their relative importance still remain unclear. To gain insight into these dispersive mechanisms, we track liquid-suspended microparticles and extract their effective diffusivities inside an anatomically-inspired microfluidic acinar model. Such results are then compared to experiments and numerical simulations in a straight channel. While alveoli of the proximal acinar generations exhibit convective mixing characteristics that lead to irreversible particle trajectories, this local effect is overshadowed by a more dominant dispersion mechanism across the ductal branching network that arises from small but significant streamline crossing due to intrinsic diffusional motion in the presence of high velocity gradients. We anticipate that for true airborne particles, which exhibit much higher intrinsic motion, streamline crossing would be even more significant.

A Microfluidic Model of Biomimetically Breathing Pulmonary Acinar Airways

Quantifying respiratory flow characteristics in the pulmonary acinar depths and how they influence inhaled aerosol transport is critical towards optimizing drug inhalation techniques as well as predicting deposition patterns of potentially toxic airborne particles in the pulmonary alveoli. Here, soft-lithography techniques are used to fabricate complex acinar-like airway structures at the truthful anatomical length-scales that reproduce physiological acinar flow phenomena in an optically accessible system. The microfluidic device features 5 generations of bifurcating alveolated ducts with periodically expanding and contracting walls. Wall actuation is achieved by altering the pressure inside water-filled chambers surrounding the thin PDMS acinar channel walls both from the sides and the top of the device. In contrast to common multilayer microfluidic devices, where the stacking of several PDMS molds is required, a simple method is presented to fabricate the top chamber by embedding the barrel section of a syringe into the PDMS mold. This novel microfluidic setup delivers physiological breathing motions which in turn give rise to characteristic acinar air-flows. In the current study, micro particle image velocimetry (µPIV) with liquid suspended particles was used to quantify such air flows based on hydrodynamic similarity matching. The good agreement between µPIV results and expected acinar flow phenomena suggest that the microfluidic platform may serve in the near future as an attractive in vitro tool to investigate directly airborne representative particle transport and deposition in the acinar regions of the lungs.

Revisiting pulmonary acinar particle transport: convection, sedimentation, diffusion, and their interplay

It is largely acknowledged that inhaled particles ranging from 0.001 to 10 m are able to reach and deposit in the alveolated regions of the lungs. To date, however, the bulk of numerical studies have focused mainly on micrometer sized particles whose transport kinematics are governed by convection and sedimentation, thereby capturing only a small fraction of the wider range of aerosols leading to acinar deposition. Too little is still known about the local acinar transport dynamics of inhaled (ultra)fine particles affected by diffusion and convection. Our study aims to fill this gap by numerically simulating the transport characteristics of particle sizes spanning three orders of magnitude (0.01-5 m) covering diffusive, convective, and gravitational aerosol motion across a multigenerational acinar network. By characterizing the deposition patterns as a function of particle size, we find that submicrometer particles [formulae see text (0.1 m)] reach deep into the acinar structure and are prone to deposit near alveolar openings; meanwhile, other particle sizes are restricted to accessing alveolar cavities in proximal generations. Our findings underline that a precise understanding of acinar aerosol transport, and ultrafine particles in particular, is contingent upon resolving the complex convective-diffusive interplay in determining their irreversible kinematics and local deposition sites.

Observation of cardiogenic flow oscillations in healthy subjects with hyperpolarized 3He MRI

Recently, dynamic MRI of hyperpolarized (3)He during inhalation revealed an alternation of the image intensity between left and right lungs with a cardiac origin (Sun Y, Butler JP, Ferrigno M, Albert MS, Loring SH. Respir Physiol Neurobiol 185: 468-471, 2013). This effect is investigated further using dynamic and phase-contrast flow MRI with inhaled (3)He during slow inhalations (flow rate ∼100 ml/s) to elucidate airflow dynamics in the main lobes in six healthy subjects. The ventilation MR signal and gas inflow in the left lower lobe (LLL) of the lungs were found to oscillate clearly at the cardiac frequency in all subjects, whereas the MR signals in the other parts of the lungs had a similar oscillatory behavior but were smaller in magnitude and in anti-phase to the signal in the left lower lung. The airflow in the main bronchi showed periodic oscillations at the frequency of the cardiac cycle. In four of the subjects, backflows were observed for a short period of time of the cardiac cycle, demonstrating a pendelluft effect at the carina bifurcation between the left and right lungs. Additional (1)H structural MR images of the lung volume and synchronized ECG recording revealed that maximum inspiratory flow rates in the LLL of the lungs occurred during systole when the corresponding left lung volume increased, whereas the opposite effect was observed during diastole, with gas flow redirected to the other parts of the lung. In conclusion, cardiogenic flow oscillations have a significant effect on regional gas flow and distribution within the lungs.

Particle dynamics and deposition in true-scale pulmonary acinar models

Particle transport phenomena in the deep alveolated airways of the lungs (i.e. pulmonary acinus) govern deposition outcomes following inhalation of hazardous or pharmaceutical aerosols. Yet, there is still a dearth of experimental tools for resolving acinar particle dynamics and validating numerical simulations. Here, we present a true-scale experimental model of acinar structures consisting of bifurcating alveolated ducts that capture breathing-like wall motion and ensuing respiratory acinar flows. We study experimentally captured trajectories of inhaled polydispersed smoke particles (0.2 to 1 μm in diameter), demonstrating how intrinsic particle motion, i.e. gravity and diffusion, is crucial in determining dispersion and deposition of aerosols through a streamline crossing mechanism, a phenomenon paramount during flow reversal and locally within alveolar cavities. A simple conceptual framework is constructed for predicting the fate of inhaled particles near an alveolus by identifying capture and escape zones and considering how streamline crossing may shift particles between them. In addition, we examine the effect of particle size on detailed deposition patterns of monodispersed microspheres between 0.1–2 μm. Our experiments underline local modifications in the deposition patterns due to gravity for particles ≥0.5 μm compared to smaller particles, and show good agreement with corresponding numerical simulations.

Investigating the pharmacokinetics and biological distribution of silver-loaded polyphosphoester-based nanoparticles using (111) Ag as a radiotracer

Purified (111) Ag was used as a radiotracer to investigate silver loading and release, pharmacokinetics, and biodistribution of polyphosphoester-based degradable shell crosslinked knedel-like (SCK) nanoparticles as a comparison to the previously reported small molecule, N-heterocyclic silver carbene complex analog (SCC1) for the delivery of therapeutic silver ions in mouse models. Biodistribution studies were conducted by aerosol administration of (111) Ag acetate, [(111) Ag]SCC1, and [(111) Ag]SCK doses directly into the lungs of C57BL/6 mice. Nebulization of the (111) Ag antimicrobials resulted in an average uptake of 1.07 ± 0.12% of the total aerosolized dose given per mouse. The average dose taken into the lungs of mice was estimated to be 2.6 ± 0.3% of the dose inhaled per mouse for [(111) Ag]SCC1 and twice as much dose was observed for the [(111) Ag]SCKs (5.0 ± 0.3% and 5.9 ± 0.8% for [(111) Ag]aSCK and [(111) Ag]zSCK, respectively) at 1 h post administration (p.a.). [(111) Ag]SCKs also exhibited higher dose retention in the lungs; 62-68% for [(111) Ag]SCKs and 43% for [(111) Ag]SCC1 of the initial 1 h dose were observed in the lungs at 24 h p.a.. This study demonstrates the utility of (111) Ag as a useful tool for monitoring the pharmacokinetics of silver-loaded antimicrobials in vivo.

Deposition of Particles in the Alveolar Airways: Inhalation and Breath-Hold with Pharmaceutical Aerosols

Previous studies have demonstrated that factors such as airway wall motion, inhalation waveform, and geometric complexity influence the deposition of aerosols in the alveolar airways. However, deposition fraction correlations are not available that account for these factors in determining alveolar deposition. The objective of this study was to generate a new space-filling model of the pulmonary acinus region and implement this model to develop correlations of aerosol deposition that can be used to predict the alveolar dose of inhaled pharmaceutical products. A series of acinar models was constructed containing different numbers of alveolar duct generations based on space-filling 14-hedron elements. Selected ventilation waveforms were quick-and-deep and slow-and-deep inhalation consistent with the use of most pharmaceutical aerosol inhalers. Computational fluid dynamics simulations were used to predict aerosol transport and deposition in the series of acinar models across various orientations with gravity where ventilation was driven by wall motion. Primary findings indicated that increasing the number of alveolar duct generations beyond 3 had a negligible impact on total acinar deposition, and total acinar deposition was not affected by gravity orientation angle. A characteristic model containing three alveolar duct generations (D3) was then used to develop correlations of aerosol deposition in the alveolar airways as a function of particle size and particle residence time in the geometry. An alveolar deposition parameter was determined in which deposition correlated with d²t over the first half of inhalation followed by correlation with dt², where d is the aerodynamic diameter of the particles and t is the potential particle residence time in the alveolar model. Optimal breath-hold times to allow 95% deposition of inhaled 1, 2, and 3 μm particles once inside the alveolar region were approximately >10, 2.7, and 1.2 s, respectively. Coupling of the deposition correlations with previous stochastic individual path (SIP) model predictions of tracheobronchial deposition was demonstrated to predict alveolar dose of commercial pharmaceutical products. In conclusion, this study completes an initiative to determine the fate of inhaled pharmaceutical aerosols throughout the respiratory airways using CFD simulations.

Role of Alveolar Topology on Acinar Flows and Convective Mixing

Due to experimental challenges, computational simulations are often sought to quantify inhaled aerosol transport in the pulmonary acinus. Commonly, these are performed using generic alveolar topologies, including spheres, toroids, and polyhedra, to mimic the complex acinar morphology. Yet, local acinar flows and ensuing particle transport are anticipated to be influenced by the specific morphological structures. We have assessed a range of acinar models under self-similar breathing conditions with respect to alveolar flow patterns, convective flow mixing, and deposition of fine particles (1.3 μm diameter). By tracking passive tracers over cumulative breathing cycles, we find that irreversible flow mixing correlates with the location and strength of the recirculating vortex inside the cavity. Such effects are strongest in proximal acinar generations where the ratio of alveolar to ductal flow rates is low and interalveolar disparities are most apparent. Our results for multi-alveolated acinar ducts highlight that fine 1 μm inhaled particles subject to alveolar flows are sensitive to the alveolar topology, underlining interalveolar disparities in particle deposition patterns. Despite the simplicity of the acinar models investigated, our findings suggest that alveolar topologies influence more significantly local flow patterns and deposition sites of fine particles for upper generations emphasizing the importance of the selected acinar model. In distal acinar generations, however, the alveolar geometry primarily needs to mimic the space-filling alveolar arrangement dictated by lung morphology.