CFD Simulation and Experimental Validation of Fluid Flow and Particle Transport in a Model of Alveolated Airways

Optimizing cabin air inlet velocities and personal risk assessment: Introducing the Personal Contamination Ratio (PCR) method for enhanced aircraft cabin infection risk evaluation

Recurrent epidemics of respiratory infections have drawn attention from the academic community and the general public in recent years. Aircraft plays a pivotal role in facilitating the cross-regional transmission of pathogens. In this study, we initially utilized an Airbus A320 model for computational fluid dynamics (CFD) simulations, subsequently validating the model's efficacy in characterizing cabin airflow patterns through comparison with empirical data. Building upon this validated framework, we investigate the transport dynamics of droplets of varying sizes under three air supply velocities. The Euler-Lagrangian method is employed to meticulously track key parameters associated with droplet transport, enabling a comprehensive analysis of particle behavior within the cabin environment. This study integrates acquired data into a novel PCR (Personal Contamination Rate) equation to assess individual contamination rates. Numerical simulations demonstrate that increasing air supply velocity leads to enhanced stability in the movement of larger particles compared to smaller ones. Results show that the number of potential infections in the cabin decreases by 51.8% at the highest air supply velocity compared to the base air supply velocity, and the total exposure risk rate reduced by 26.4%. Thus, optimizing air supply velocity within a specific range effectively reduces the potential infection area. In contrast to previous research, this study provides a more comprehensive analysis of droplet movement dynamics across various particle sizes. We introduce an improved method for calculating the breathing zone, thereby enhancing droplet counting accuracy. These findings have significant implications for improving non-pharmacological public health interventions and optimizing cabin ventilation system design.

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.

Bioengineering and Clinical Translation of Human Lung and its Components

Clinical lung transplantation has rapidly established itself as the gold standard of treatment for end-stage lung diseases in a restricted group of patients since the first successful lung transplant occurred. Although significant progress has been made in lung transplantation, there are still numerous obstacles on the path to clinical success. The development of bioartificial lung grafts using patient-derived cells may serve as an alternative treatment modality; however, challenges include developing appropriate scaffold materials, advanced culture strategies for lung-specific multiple cell populations, and fully matured constructs to ensure increased transplant lifetime following implantation. This review highlights the development of tissue-engineered tracheal and lung equivalents over the past two decades, key problems in lung transplantation in a clinical environment, the advancements made in scaffolds, bioprinting technologies, bioreactors, organoids, and organ-on-a-chip technologies. The review aims to fill the lacuna in existing literature toward a holistic bioartificial lung tissue, including trachea, capillaries, airways, bifurcating bronchioles, lung disease models, and their clinical translation. Herein, the efforts are on bridging the application of lung tissue engineering methods in a clinical environment as it is thought that tissue engineering holds enormous promise for overcoming the challenges associated with the clinical translation of bioengineered human lung and its components.

Microparticle Transport and Sedimentation in a Rhythmically Expanding Alveolar Chip

Understanding the mechanism of particle transport and sedimentation in pulmonary alveolus is important for deciphering the causes of respiratory diseases and helping the development of drug delivery. In this study, taking advantage of the microfluidic technique, an experimental platform was developed to study particle behavior in a rhythmically expanding alveolar chip for a sufficient number of cycles. The alveolar flow patterns at different generations were measured for two cases with the gravity direction parallel or vertical to the alveolar duct. Affected by both the vortex flow inside the alveoli and the shear flow in the duct simultaneously, it was observed that particles inside the alveoli either escaped from the inlet of the alveolar duct or stayed in the alveoli, revealing the irreversibility of particle transport in the alveoli. At the earlier acinar generations, particles were inclined to deposit on the distal alveolar wall. The settling rates of particles of different sizes in the alveoli were also compared. This study provides valuable data for understanding particle transport and sedimentation in the alveoli.

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.

Visualizing the Flow Patterns in an Expanding and Contracting Pulmonary Alveolated Duct Based on Microcomputed Tomography Images

We visualized the flow patterns in an alveolated duct model with breathing-like expanding and contracting wall motions using particle image velocimetry, and then, we investigated the effect of acinar deformation on the flow patterns. We reconstructed a compliant, scaled-up model of an alveolated duct from synchrotron microcomputed tomography images of a mammalian lung. The alveolated duct did not include any bifurcation, and its entire surface was covered with alveoli. We embedded the alveolated duct in a sealed container that was filled with fluid. We oscillated the fluid in the duct and container simultaneously and independently to control the flow and duct volume. We examined the flow patterns in alveoli, with the Reynolds number (Re) at 0.03 or 0.22 and the acinar volume change at 0%, 20%, or 80%. At the same Re, the heterogeneous deformation induced different inspiration and expiration flow patterns, and the recirculating regions in alveoli changed during respiratory cycle. During a larger acinar deformation at Re = 0.03, the flow patterns tended to change from recirculating flow to radial flow during inspiration and vice versa during expiration. Additionally, the alveolar geometric characteristics, particularly the angle between the alveolar duct and mouth, affected these differences in flow patterns. At Re = 0.22, recirculating flow patterns tended to form during inspiration and expiration, regardless of the magnitude of the acinar deformation. Our in vitro experiments suggest that the alveolated flows with nonself-similar and heterogeneous wall motions may promote particle mixing and deposition.

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.

Airborne Transmission of COVID-19: Aerosol Dispersion, Lung Deposition, and Virus-Receptor Interactions

Coronavirus disease 2019 (COVID-19), due to infection by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is now causing a global pandemic. Aerosol transmission of COVID-19, although plausible, has not been confirmed by the World Health Organization (WHO) as a general transmission route. Considering the rapid spread of SARS-CoV-2, especially nosocomial outbreaks and other superspreading events, there is an urgent need to study the possibility of airborne transmission and its impact on the lung, the primary body organ attacked by the virus. Here, we review the complete pathway of airborne transmission of SARS-CoV-2 from aerosol dispersion in air to subsequent biological uptake after inhalation. In particular, we first review the aerodynamic and colloidal mechanisms by which aerosols disperse and transmit in air and deposit onto surfaces. We then review the fundamental mechanisms that govern regional deposition of micro- and nanoparticles in the lung. Focus is given to biophysical interactions between particles and the pulmonary surfactant film, the initial alveolar-capillary barrier and first-line host defense system against inhaled particles and pathogens. Finally, we summarize the current understanding about the structural dynamics of the SARS-CoV-2 spike protein and its interactions with receptors at the atomistic and molecular scales, primarily as revealed by molecular dynamics simulations. This review provides urgent and multidisciplinary knowledge toward understanding the airborne transmission of SARS-CoV-2 and its health impact on the respiratory system.

Efficient Low Shear Flow-based Trapping of Biological Entities

Capturing cells or biological entities is an important and challenging step toward in-vitro studies of cells under a precisely controlled microscale environment. In this work, we have developed a compact and efficient microdevice for on-chip trapping of micro-sized particles. This hydrodynamics-based trapping system allows the isolation of polystyrene micro-particles with a shorter time while inducing a less hydrodynamic deformation and stress on the particles or cells both after and before trapping. A numerical simulation was carried out to design a hydrodynamic trapping mechanism and optimize the geometric and fluidic parameters affecting the trapping efficiency of the microfluidic network. By using the finite element analysis, the velocity field, pressure field, and hydrodynamic force on the micro particles were studied. Finally, a PDMS microfluidic device was fabricated to test the device's ability to trap polystyrene microspheres. Computational fluid analysis and experimental testing showed a high trapping efficiency that is more than 90%. This microdevice can be used for single cell studies including their biological, physical and chemical characterization.

A Novel Approach for 3D-Structural Identification through Video Recording: Magnified Tracking

Advancements in optical imaging devices and computer vision algorithms allow the exploration of novel diagnostic techniques for use within engineering systems. A recent field of application lies in the adoption of such devices for non-contact vibrational response recordings of structures, allowing high spatial density measurements without the burden of heavy cabling associated with conventional technologies. This, however, is not a straightforward task due to the typically low-amplitude displacement response of structures under ambient operational conditions. A novel framework, namely Magnified Tracking (MT), is proposed herein to overcome this limitation through the synergistic use of two computer vision techniques. The recently proposed phase-based motion magnification (PBMM) framework, for amplifying motion in a video within a defined frequency band, is coupled with motion tracking by means of particle tracking velocimetry (PTV). An experimental campaign was conducted to validate a proof-of-concept, where the dynamic response of a shear frame was measured both by conventional sensors as well as a video camera setup, and cross-compared to prove the feasibility of the proposed non-contact approach. The methodology was explored both in 2D and 3D configurations, with PTV revealing a powerful tool for the measurement of perceptible motion. When MT is utilized for tracking “imperceptible” structural responses (i.e., below PTV sensitivity), via the use of PBMM around the resonant frequencies of the structure, the amplified motion reveals the operational deflection shapes, which are otherwise intractable. The modal results extracted from the magnified videos, using PTV, demonstrate MT to be a viable non-contact alternative for 3D modal identification with the benefit of a spatially dense measurement grid.

Including auditory tube function on models is relevant to assess water exposure after tympanostomy tubes-Multiphase computerized fluid dynamics model

INTRODUCTION

Myringotomy with tympanostomy tube is the most common otologic surgery and some patients are still advised to avoid water. However, there is no evidence supporting this, with published papers questioning the need for this advice.

METHODS

A Multiphase Computational Fluid Dynamics (CFD) model was created using computerized tomography images of a child's healthy ear. It was then used to study the flow of fluids through the external ear, tympanic cavity, and auditory tube, with and without submersion.

RESULTS

The model accurately described the behavior of the air retained in the patient's nasopharynx and tympanic cavity. A simulated elevation of pressure in the external auditory canal without submersion, without increase of pressure in the nasopharynx, demonstrated that fluids promptly crossed the tympanostomy tube into the middle ear. However, simulated elevation of pressure in the external auditory canal with concurrent elevation of air pressure in the nasopharynx during submersion did not lead to passive tube opening nor to any detectable flow through the tympanostomy tube.

CONCLUSIONS

In the model, submersion increases pressure in the nasopharynx which offsets the pressure in the external auditory canal. So, in the absence of a pressure gradient, no passive tubal opening took place, and no air or fluid flow was detected through the transtympanic tube. This model now includes the exhaust function of the auditory tube in the model and shows its relevance.

Computational Fluid Dynamics Modeling of the Human Pulmonary Arteries with Experimental Validation

Pulmonary hypertension (PH) is a chronic progressive disease characterized by elevated pulmonary arterial pressure, caused by an increase in pulmonary arterial impedance. Computational fluid dynamics (CFD) can be used to identify metrics representative of the stage of PH disease. However, experimental validation of CFD models is often not pursued due to the geometric complexity of the model or uncertainties in the reproduction of the required flow conditions. The goal of this work is to validate experimentally a CFD model of a pulmonary artery phantom using a particle image velocimetry (PIV) technique. Rapid prototyping was used for the construction of the patient-specific pulmonary geometry, derived from chest computed tomography angiography images. CFD simulations were performed with the pulmonary model with a Reynolds number matching those of the experiments. Flow rates, the velocity field, and shear stress distributions obtained with the CFD simulations were compared to their counterparts from the PIV flow visualization experiments. Computationally predicted flow rates were within 1% of the experimental measurements for three of the four branches of the CFD model. The mean velocities in four transversal planes of study were within 5.9 to 13.1% of the experimental mean velocities. Shear stresses were qualitatively similar between the two methods with some discrepancies in the regions of high velocity gradients. The fluid flow differences between the CFD model and the PIV phantom are attributed to experimental inaccuracies and the relative compliance of the phantom. This comparative analysis yielded valuable information on the accuracy of CFD predicted hemodynamics in pulmonary circulation models.

Numerical study of the effects of bronchial structural abnormalities on respiratory flow distribution

BACKGROUND

The anatomical configurations of respiratory tract would be directly associated with their ventilatory function. It is necessary to fully understand the association between airway configurations and their functions as well as the interactions between different airway segments. In this study, we developed a respiratory airway model to investigate the effects of bronchial structural abnormalities on flow distribution in the bronchi and upper airway.

METHODS

Derived from computed tomography (CT) scanner data, three-dimensional (3D) finite element (FE) models of healthy human respiratory tracts were developed with anatomically realistic configurations, including the nasal cavity, oral cavity, pharynx, larynx, trachea, and partial bronchi. Abnormal bronchial configurations were built to correspond to four common bronchial diseases. Through numerical simulation, airway configurations of normal and abnormal bronchi were obtained, and flow patterns were compared between normal and abnormal respiratory tracts, as well as the effects of lower airway changes on flow distribution in the upper airway.

RESULTS

The simulation results showed that during inspiration, abnormal bronchial structures can cause flow redistribution in each generation of bronchi and have significant effects on flow distribution in the daughter bronchi of abnormal segments, but no effect on flow distribution of the upper airway. During expiration, abnormal bronchus structures had a remarkable influence on flow distribution in the trachea, while there was no significant difference in flow distribution when airflow passed from the vocal cords and entered the laryngeal cavity.

CONCLUSIONS

Therefore, abnormal bronchial structures can affect the downstream flow distribution and cause flow redistribution throughout the entire bronchial branches. During expiration, the configurations of the trachea and glottis can diminish the effects of abnormal bronchial structures on flow distribution.

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

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

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.

Systems physiology of the airways in health and obstructive pulmonary disease

Fresh air entering the mouth and nose is brought to the blood-gas barrier in the lungs by a repetitively branching network of airways. Provided the individual airway branches remain patent, this airway tree achieves an enormous amplification in cross-sectional area from the trachea to the terminal bronchioles. Obstructive lung diseases such as asthma occur when airway patency becomes compromised. Understanding the pathophysiology of these obstructive diseases thus begins with a consideration of the factors that determine the caliber of an individual airway, which include the force balance between the inward elastic recoil of the airway wall, the outward tethering forces of its parenchymal attachments, and any additional forces due to contraction of airway smooth muscle. Other factors may also contribute significantly to airway narrowing, such as thickening of the airway wall and accumulation of secretions in the lumen. Airway obstruction becomes particularly severe when these various factors occur in concert. However, the effect of airway abnormalities on lung function cannot be fully understood only in terms of what happens to a single airway because narrowing throughout the airway tree is invariably heterogeneous and interdependent. Obstructive lung pathologies thus manifest as emergent phenomena arising from the way in which the airway tree behaves a system. These emergent phenomena are studied with clinical measurements of lung function made by spirometry and by mechanical impedance measured with the forced oscillation technique. Anatomically based computational models are linking these measurements to underlying anatomic structure in systems physiology terms. WIREs Syst Biol Med 2016, 8:423-437. doi: 10.1002/wsbm.1347 For further resources related to this article, please visit the WIREs website.

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.

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.

Controlled, Parametric, Individualized, 2-D and 3-D Imaging Measurements of Aerosol Deposition in the Respiratory Tract of Asthmatic Human Subjects for Model Validation

BACKGROUND

Computer modeling is used to predict inhaled aerosol deposition in the lung based on definition of the aerosol characteristics and the breathing pattern and airway anatomy of the subject. Validation of the models is limited by the lack of detailed experimental data. Three-dimensional imaging provides an opportunity to address this unmet need.

METHODS

Radioactive aerosol was administered to six male asthmatic subjects on two occasions under carefully monitored input conditions. Input parameters varied in particle size, depth of breathing, and carrier gas. The aerosol distribution was measured by combined single photon emission computed tomography and x-ray computer tomography (SPECT/CT) and airway anatomy by high resolution CT. The deposition distribution was measured by both a 2D and 3D analysis and described in terms of the percentage of inhaled aerosol deposited in sections of the respiratory tract and in both spatial and anatomical subdivisions within each lung. The percentage deposition in the conducting airways was also assessed by 24 h clearance.

RESULTS

A set of imaging data of aerosol deposition has thus been produced in which the input parameters of inhalation are well described. The results in asthmatics were compared to previous measurements in healthy controls using an identical inhalation protocol. The percentages of deposition in extra-thoracic and thoracic compartments of the airways were not significantly affected by disease, but the regional pulmonary deposition pattern was, with asthma leading to increased deposition in the conducting airways.

CONCLUSIONS

The dataset acquired in this study will be useful in validating computer models of aerosol deposition in asthmatic subjects. Asthma did not affect the fraction of inhaled aerosol depositing in the lungs, but gave rise to a more central deposition pattern. The use of 3D SPECT imaging in combination with 24 h clearance measurements enables differentiation of deposition between bronchial and bronchiolar airways.

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.

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.

Biomimetics of fetal alveolar flow phenomena using microfluidics

At the onset of life in utero, the respiratory system begins as a liquid-filled tubular organ and undergoes significant morphological changes during fetal development towards establishing a respiratory organ optimized for gas exchange. As airspace morphology evolves, respiratory alveolar flows have been hypothesized to exhibit evolving flow patterns. In the present study, we have investigated flow topologies during increasing phases of embryonic life within an anatomically inspired microfluidic device, reproducing real-scale features of fetal airways representative of three distinct phases of in utero gestation. Micro-particle image velocimetry measurements, supported by computational fluid dynamics simulations, reveal distinct respiratory alveolar flow patterns throughout different stages of fetal life. While attached, streamlined flows characterize the shallow structures of premature alveoli indicative of the onset of saccular stage, separated recirculating vortex flows become the signature of developed and extruded alveoli characteristic of the advanced stages of fetal development. To further mimic physiological aspects of the cellular environment of developing airways, our biomimetic devices integrate an alveolar epithelium using the A549 cell line, recreating a confluent monolayer that produces pulmonary surfactant. Overall, our in vitro biomimetic fetal airways model delivers a robust and reliable platform combining key features of alveolar morphology, flow patterns, and physiological aspects of fetal lungs developing in utero.

Optimization of flow assisted entrapment of pollen grains in a microfluidic platform for tip growth analysis

A biocompatible polydimethylsiloxane (PDMS) biomicrofluidic platform is designed, fabricated and tested to study protuberance growth of single plant cells in a micro-vitro environment. The design consists of an inlet to introduce the cell suspension into the chip, three outlets to conduct the medium or cells out of the chip, a main distribution chamber and eight microchannels connected to the main chamber to guide the growth of tip growing plant cells. The test cells used here were pollen grains which produce cylindrical protrusions called pollen tubes. The goal was to adjust the design of the microfluidic network with the aim to enhance the uniformly distributed positioning of pollen grains at the entrances of the microchannels and to provide identical fluid flow conditions for growing pollen tubes along each microchannel. Computational fluid analysis and experimental testing were carried out to estimate the trapping efficiencies of the different designs.

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.

Airway gas flow

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

Computational and bioengineered lungs as alternatives to whole animal, isolated organ, and cell-based lung models

Development of lung models for testing a drug substance or delivery system has been an intensive area of research. However, a model that mimics physiological and anatomical features of human lungs is yet to be established. Although in vitro lung models, developed and fine-tuned over the past few decades, were instrumental for the development of many commercially available drugs, they are suboptimal in reproducing the physiological microenvironment and complex anatomy of human lungs. Similarly, intersubject variability and high costs have been major limitations of using animals in the development and discovery of drugs used in the treatment of respiratory disorders. To address the complexity and limitations associated with in vivo and in vitro models, attempts have been made to develop in silico and tissue-engineered lung models that allow incorporation of various mechanical and biological factors that are otherwise difficult to reproduce in conventional cell or organ-based systems. The in silico models utilize the information obtained from in vitro and in vivo models and apply computational algorithms to incorporate multiple physiological parameters that can affect drug deposition, distribution, and disposition upon administration via the lungs. Bioengineered lungs, on the other hand, exhibit significant promise due to recent advances in stem or progenitor cell technologies. However, bioengineered approaches have met with limited success in terms of development of various components of the human respiratory system. In this review, we summarize the approaches used and advancements made toward the development of in silico and tissue-engineered lung models and discuss potential challenges associated with the development and efficacy of these models.

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

The aerosol bolus technique can be used to estimate the degree of convective mixing in the lung; however, contributions of different lung compartments to measured dispersion cannot be differentiated unambiguously. To estimate dispersion in the distal lung, we studied the effect of gravity and airway asymmetry on the dispersion of 1 μm-diameter particle boluses in three-dimensional computational models of the lung periphery, ranging from a single alveolar sac to four-generation (g4) structures of bifurcating airways that deformed homogeneously during breathing. Boluses were introduced at the beginning of a 2-s inhalation, immediately followed by a 3-s exhalation. Dispersion was estimated by the half-width of the exhaled bolus. Dispersion was significantly affected by the spatial orientation of the models in normal gravity and was less in zero gravity than in normal gravity. Dispersion was strongly correlated with model volume in both normal and zero gravity. Predicted pulmonary dispersion based on a symmetric g4 acinar model was 391 ml and 238 ml under normal and zero gravity, respectively. These results accounted for a significant amount of dispersion measured experimentally. In zero gravity, predicted dispersion in a highly asymmetric model accounted for ∼20% of that obtained in a symmetric model with comparable volume and number of alveolated branches, whereas normal gravity dispersions were comparable in both models. These results suggest that gravitational sedimentation and not geometrical asymmetry is the dominant factor in aerosol dispersion in the lung periphery.

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.

Multi-scale lung modeling

Multi-scale modeling of biological systems has recently become fashionable due to the growing power of digital computers as well as to the growing realization that integrative systems behavior is as important to life as is the genome. While it is true that the behavior of a living organism must ultimately be traceable to all its components and their myriad interactions, attempting to codify this in its entirety in a model misses the insights gained from understanding how collections of system components at one level of scale conspire to produce qualitatively different behavior at higher levels. The essence of multi-scale modeling thus lies not in the inclusion of every conceivable biological detail, but rather in the judicious selection of emergent phenomena appropriate to the level of scale being modeled. These principles are exemplified in recent computational models of the lung. Airways responsiveness, for example, is an organ-level manifestation of events that begin at the molecular level within airway smooth muscle cells, yet it is not necessary to invoke all these molecular events to accurately describe the contraction dynamics of a cell, nor is it necessary to invoke all phenomena observable at the level of the cell to account for the changes in overall lung function that occur following methacholine challenge. Similarly, the regulation of pulmonary vascular tone has complex origins within the individual smooth muscle cells that line the blood vessels but, again, many of the fine details of cell behavior average out at the level of the organ to produce an effect on pulmonary vascular pressure that can be described in much simpler terms. The art of multi-scale lung modeling thus reduces not to being limitlessly inclusive, but rather to knowing what biological details to leave out.