Flow and particle dispersion in a pulmonary alveolus--part I: velocity measurements and convective particle transport

Design of a multilayer lung chip with multigenerational alveolar ducts to investigate the inhaled particle deposition

We present the development and application of a multilayer microfluidic lung chip designed to accurately replicate the human respiratory bronchi, providing an innovative platform for controlled particle deposition in the lung. By employing a quantitative control method of fluid velocity through the deformation of an elastic PDMS membrane, this platform mimics the passive breathing process in humans and allows for precise simulation of the respiration cycle. We utilized time-lapse photography of fluorescent particles in a water/glycerol solution to qualitatively observe fluid morphology in the channel, while a chip-aerosol exposure device combined with microscopy imaging was employed to visualise aerosol deposition. Both experimental and numerical simulation results showed that particle concentration decreased towards the distal generations of the lung, and that changes in breathing pattern significantly affected particle deposition trends. Furthermore, we found that increasing the residence time of particles in the channel facilitated deeper particle deposition, achievable by adjusting parameters such as breath-hold time, exhalation time, respiration cycle length, and tidal volume. The proposed microfluidic lung chip device has significant potential for future research in respiratory health and inhaled drug delivery, providing an efficient, cost-effective, and ethical alternative to traditional in vivo studies.

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.

Revisiting Airflow and Aerosol Transport Phenomena in the Deep Lungs with Microfluidics

The dynamics of respiratory airflows and the associated transport mechanisms of inhaled aerosols characteristic of the deep regions of the lungs are of broad interest in assessing both respiratory health risks and inhalation therapy outcomes. In the present review, we present a comprehensive discussion of our current understanding of airflow and aerosol transport phenomena that take place within the unique and complex anatomical environment of the deep lungs, characterized by submillimeter 3D alveolated airspaces and nominally slow resident airflows, known as low-Reynolds-number flows. We exemplify the advances brought forward by experimental efforts, in conjunction with numerical simulations, to revisit past mechanistic theories of respiratory airflow and particle transport in the distal acinar regions. Most significantly, we highlight how microfluidic-based platforms spanning the past decade have accelerated opportunities to deliver anatomically inspired in vitro solutions that capture with sufficient realism and accuracy the leading mechanisms governing both respiratory airflow and aerosol transport at true scale. Despite ongoing challenges and limitations with microfabrication techniques, the efforts witnessed in recent years have provided previously unattainable in vitro quantifications on the local transport properties in the deep pulmonary acinar airways. These may ultimately provide new opportunities to explore improved strategies of inhaled drug delivery to the deep acinar regions by investigating further the mechanistic interactions between airborne particulate carriers and respiratory airflows at the pulmonary microscales.

Study on Deposition Characteristics of Microparticles in Terminal Pulmonary Acini by IB-LBM

As an indicator of health risk, the deposition of microparticles in terminal pulmonary acini is of great significance in the medical field. To control particulate pollution and optimize aerosol delivery, it is necessary to perform an in-depth study of the microparticle deposition in terminal pulmonary acini; however, little research has been done on this topic. This paper proposes a respiratory movement model of terminal pulmonary acini using an immersed boundary–lattice Boltzmann method. In addition, we explored the effect of gravity direction, respiratory rate, microparticle diameter, and other parameters on the microparticles deposition process and distribution, under the airflow in the acinar wall. It was found that the deposition of microparticles is sensitive to gravity direction, and the growth of the respiratory rate increases the rate of microparticle migration and deposition. It was observed that the gravity effect is enhanced by increasing the diameter of microparticles, causing a high deposition and dispersion rate. The study reveals the dynamic correlation between the respiration process and the movement of microparticles, which is of reference value to figure out the pathogenicity mechanism of inhalable particles and to optimize the aerosol delivery.

Investigation on Microparticle Transport and Deposition Mechanics in Rhythmically Expanding Alveolar Chip

The transport and deposition of micro/nanoparticles in the lungs under respiration has an important impact on human health. Here, we presented a real-scale alveolar chip with movable alveolar walls based on the microfluidics to experimentally study particle transport in human lung alveoli under rhythmical respiratory. A new method of mixing particles in aqueous solution, instead of air, was proposed for visualization of particle transport in the alveoli. Our novel design can track the particle trajectories under different force conditions for multiple periods. The method proposed in this study gives us better resolution and clearer images without losing any details when mapping the particle velocities. More detailed particle trajectories under multiple forces with different directions in an alveolus are presented. The effects of flow patterns, drag force, gravity and gravity directions are evaluated. By tracing the particle trajectories in the alveoli, we find that the drag force contributes to the reversible motion of particles. However, compared to drag force, the gravity is the decisive factor for particle deposition in the alveoli.

Septal destruction enhances chaotic mixing and increases cellular doses of nanoparticles in emphysematous acinus

One hallmark of emphysema is the breakdown of inter-alveolar septal walls in pulmonary acini. How the acinar dosimetry of environmental aerosols varies at different stages of emphysema remains unclear; this is specifically pertinent to users of tobacco products, which is the leading cause of emphysema. The objective of this study is to systematically assess the impacts of septal destruction on the behavior and fate of nanoparticles (1–800 nm) in a pyramid-shaped sub-acinar model consisting of 496 alveoli. Four diseased geometry variants were created by gradually removing the septal walls from the base model. Particle motions within the acinar region were tracked for particles raging 1–800 nm at four emphysema stages using a well-tested Lagrangian tracking model. Both spatial profile and temporal variation of particle deposition were predicted in healthy and diseased sub-acinar geometries on both a total and regional basis. Results show large differences in airflow and particle dynamics among different emphysema stages. Large differences in particle dynamics are also observed among different particle sizes, with one order of magnitude’s variation in the speeds of particles of 1, 10, and 200 nm. The destruction of septal walls also changed the deposition mechanisms, shifting from connective diffusion to chaotic mixing with emphysema progression. The sub-acinar dosimetry became less sensitive to particle size variation with more septal destructions. The lowest retention rate was found at 200–500 nm in the healthy sub-acinar geometry, but at 800 nm in all emphysematous models considered. The acinus-averaged dose for nanoparticles (1–800 nm) increases with aggravating septal destructions, indicating an even higher risk to the acinus at later emphysema stages.

Nanoparticle Deposition in Rhythmically Moving Acinar Models with Interalveolar Septal Apertures

Pulmonary delivery of nanomedicines has been extensively studied in recent years because of their enhanced biocompatibility, sustained-release properties, and surface modification capability. The lung as a target also offers many advantages over other routers, such as large surface area, noninvasive, quick therapeutic onset, and avoiding first-pass metabolism. However, nanoparticles smaller than 0.26 µm typically escape phagocytosis and remain in the alveoli for a long time, leading to particle accumulation and invoking tissue responses. It is imperative to understand the behavior and fates of inhaled nanoparticles in the alveoli to reliably assess therapeutic outcomes of nanomedicines or health risk of environmental toxins. The objective of this study is to numerically investigate nanoparticle deposition in a duct-alveolar model with varying sizes of inter-alveolar septal apertures (pores). A discrete phase Lagrangian model was implemented to track nanoparticle trajectories under the influence of rhythmic wall expansion and contraction. Both temporal and spatial dosimetry in the alveoli were computed. Wall motions are essential for nanoparticles to penetrate the acinar region and deposit in the alveoli. The level of aerosol irreversibility (i.e., mixing of inhaled nanoparticles with residual air in the alveolar airspace) is determined by the particle diffusivity, which in turn, dictates the fraction of particles being exhaled out. When deposition in the upper airways was not considered, high alveolar deposition rates (74–95%) were predicted for all nanoparticles considered (1–1000 nm), which were released into the alveoli at the beginning of the inhalation. The pore size notably affects the deposition pattern of inhaled nanoparticles but exerts a low impact upon the total deposition fractions. This finding indicates that consistent pulmonary doses of nanomedicine are possible in emphysema patients if breathing maneuver with the same tidal volume can be performed.

Toxicant Deposition and Transport in Alveolus: A Classical Density Functional Prediction

The deposition and transport of toxicants on pulmonary surfactant are important processes in human health and medical care. We have introduced classical density functional theory (CDFT) to provide insight into this process. Nine typical toxicants in PM2.5 were considered, and their free energy and structural information have been examined. The free energy profile indicates that PbO, As₂O₃, and CdO are the three toxicants most easily deposited in the pulmonary alveolus, which is consistent with survey data. CuO appears to be the easiest toxicant to transport through the surfactant. Structural analysis indicates that the toxicants tend to pass through the surfactant with rotation. The configuration of the pulmonary surfactant was examined by extending our previous work to polymer systems, and it appears that both the configurational entropy and the direct interaction between the surfactant and the toxicant dominate the configuration of the pulmonary surfactant.

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.

Onset of alveolar recirculation in the developing lungs and its consequence on nanoparticle deposition in the pulmonary acinus

The structure of the gas exchange region of the human lung (the pulmonary acinus) undergoes profound change in the first few years of life. In this paper, we investigate numerically how the change in alveolar shape with time affects the rate of nanoparticle deposition deep in the lung during postnatal development. As human infant data is unavailable, we use a rat model of lung development. The process of postnatal lung development in the rat is remarkably similar to that of the human, and the structure of the rat acinus is indistinguishable from that of the human acinus. The current numerical predictions support our group's recent in vivo findings, which were also obtained by using growing rat lung models, that nanoparticle deposition in infants is strongly affected by the change in the structure of the pulmonary acinus. In humans, this major structural change occurs over the first 2 yr of life. Our current predictions would suggest that human infants at the age of ∼ 2 yr might be most at risk to the harmful effects of air pollution. Our results also suggest that dose estimates for inhalation therapies using nanoparticles, based on fully developed adult lungs with simple body weight scaling, are likely to overestimate deposition by up to 55% for newborns and underestimate deposition by up to 17% for 2-yr-old infants.

The simultaneous role of an alveolus as flow mixer and flow feeder for the deposition of inhaled submicron particles

In an effort to understand the fate of inhaled submicron particles in the small sacs, or alveoli, comprising the gas-exchange region of the lung, we calculated the flow in three-dimensional (3D) rhythmically expanding models of alveolated ducts. Since convection toward the alveolar walls is a precursor to particle deposition, it was the goal of this paper to investigate the streamline maps' dependence upon alveoli location along the acinar tree. On the alveolar midplane, the recirculating flow pattern exhibited closed streamlines with a stagnation saddle point. Off the midplane we found no closed streamlines but nested, funnel-like, spiral, structures (reminiscent of Russian nesting dolls) that were directed towards the expanding walls in inspiration, and away from the contracting walls in expiration. These nested, funnel-like, structures were surrounded by air that flowed into the cavity from the central channel over inspiration and flowed from the cavity to the central channel over expiration. We also found that fluid particle tracks exhibited similar nested funnel-like spiral structures. We conclude that these unique alveolar flow structures may be of importance in enhancing deposition. In addition, due to inertia, the nested, funnel-like, structures change shape and position slightly during a breathing cycle, resulting in flow mixing. Also, each inspiration feeds a fresh supply of particle-laden air from the central channel to the region surrounding the mixing region. Thus, this combination of flow mixer and flow feeder makes each individual alveolus an effective mixing unit, which is likely to play an important role in determining the overall efficiency of convective mixing in the acinus.

Multiscale imaging and registration-driven model for pulmonary acinar mechanics in the mouse

A registration-based multiscale method to obtain a deforming geometric model of mouse acinus is presented. An intact mouse lung was fixed by means of vascular perfusion at a hydrostatic inflation pressure of 20 cmH(2)O. Microcomputed tomography (μCT) scans were obtained at multiple resolutions. Substructural morphometric analysis of a complete acinus was performed by computing a surface-to-volume (S/V) ratio directly from the 3D reconstruction of the acinar geometry. A geometric similarity is observed to exist in the acinus where S/V is approximately preserved anywhere in the model. Using multiscale registration, the shape of the acinus at an elevated inflation pressure of 25 cmH(2)O is estimated. Changes in the alveolar geometry suggest that the deformation within the acinus is not isotropic. In particular, the expansion of the acinus (from 20 to 25 cmH(2)O) is accompanied by an increase in both surface area and volume in such a way that the S/V ratio is not significantly altered. The developed method forms a useful tool in registration-driven fluid and solid mechanics studies as displacement of the alveolar wall becomes available in a discrete sense.

Flow and Particle Dispersion in Lung Acini: Effect of Geometric and Dynamic Parameters During Synchronous Ventilation

The human lung comprises about 300 million alveoli which are located on bronchioles between the 17th to 24th generations of the acinar tree, with a progressively higher population density in the deeper branches (lower acini). The alveolar size and aspect ratio change with generation number. Due to successive bifurcation, the flow velocity magnitude also decreases as the bronchiole diameter decreases from the upper to lower acini. As a result, fluid dynamic parameters such as Reynolds (Re) and Womersley (α) numbers progressively decrease with increasing generation number. In order to characterize alveolar flow patterns and inhaled particle transport during synchronous ventilation, we have conducted measurements for a range of dimensionless parameters physiologically relevant to the upper acini. Acinar airflow patterns were measured using a simplified in vitro alveolar model consisting of a single transparent elastic truncated sphere (representing the alveolus) mounted over a circular hole on the side of a rigid circular tube (representing the bronchiole). The model alveolus was capable of expanding and contracting in-phase with the oscillatory flow through the bronchiole thereby simulating synchronous ventilation. Realistic breathing conditions were achieved by exercising the model over a range of progressively varying geometric and dynamic parameters to simulate the environment within several generations of the acinar tree. Particle image velocimetry was used to measure the resulting flow patterns. Next, we used the measured flow fields to calculate particle trajectories to obtain particle transport and deposition statistics for massless and finite-size particles under the influence of flow advection and gravity. Our study shows that the geometric parameters (β and ΔV/V) primarily affect the velocity magnitudes, whereas the dynamic parameters (Re and α) distort the flow symmetry while also altering the velocity magnitudes. Consequently, the dynamic parameters have a greater influence on the particle trajectories and deposition statistics compared to the geometric parameters. The results from this study can benefit pulmonary research into the risk assessment of toxicological inhaled aerosols, and the pharmaceutical industry by providing better insight into the flow patterns and particle transport of inhalable therapeutics in the acini.

Aerosol deposition characteristics in distal acinar airways under cyclic breathing conditions

Although the major mechanisms of aerosol deposition in the lung are known, detailed quantitative data in anatomically realistic models are still lacking, especially in the acinar airways. In this study, an algorithm was developed to build multigenerational three-dimensional models of alveolated airways with arbitrary bifurcation angles and spherical alveolar shape. Using computational fluid dynamics, the deposition of 1- and 3-μm aerosol particles was predicted in models of human alveolar sac and terminal acinar bifurcation under rhythmic wall motion for two breathing conditions (functional residual capacity = 3 liter, tidal volume = 0.5 and 0.9 liter, breathing period = 4 s). Particles entering the model during one inspiration period were tracked for multiple breathing cycles until all particles deposited or escaped from the model. Flow recirculation inside alveoli occurred only during transition between inspiration and expiration and accounted for no more than 1% of the whole cycle. Weak flow irreversibility and convective transport were observed in both models. The average deposition efficiency was similar for both breathing conditions and for both models. Under normal gravity, total deposition was ~33 and 75%, of which ~67 and 96% occurred during the first cycle, for 1- and 3-μm particles, respectively. Under zero gravity, total deposition was ~2-5% for both particle sizes. These results support previous findings that gravitational sedimentation is the dominant deposition mechanism for micrometer-sized aerosols in acinar airways. The results also showed that moving walls and multiple breathing cycles are needed for accurate estimation of aerosol deposition in acinar airways.

Flow and particle dispersion in a pulmonary alveolus--part II: effect of gravity on particle transport

The acinar region of the human lung comprises about 300x10(6) alveoli, which are responsible for gas exchange between the lung and the blood. As discussed in Part I (Chhabra and Prasad, "Flow and Particle Dispersion in a Pulmonary Alveolus-Part I: Velocity Measurements and Convective Particle Transport," ASME J. Biomech. Eng., 132, p. 051009), the deposition of aerosols in the acinar region can either be detrimental to gas exchange (as in the case of harmful particulate matter) or beneficial (as in the case of inhalable pharmaceuticals). We measured the flow field inside an in-vitro model of a single alveolus mounted on a bronchiole and calculated the transport and deposition of massless particles in Part I. This paper focuses on the transport and deposition of finite-sized particles ranging from 0.25 microm to 4 microm under the combined influence of flow-induced advection (computed from velocity maps obtained by particle image velocimetry) and gravitational settling. Particles were introduced during the first inhalation cycle and their trajectories and deposition statistics were calculated for subsequent cycles for three different particle sizes (0.25 microm, 1 microm, and 4 microm) and three alveolar orientations. The key outcome of the study is that particles <or=0.25 microm follow the fluid streamlines quite closely, whereas midsize particles (d(p)=1 microm) deviate to some extent from streamlines and exhibit complex trajectories. The motion of large particles >or=4 microm is dominated by gravitational settling and shows little effect of fluid advection. Additionally, small and midsize particles deposit at about two-thirds height in the alveolus irrespective of the gravitational orientation whereas the deposition of large particles is governed primarily by the orientation of the gravity vector.