Numerical simulation of airflow and microparticle deposition in a synchrotron micro-CT-based pulmonary acinus model

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.

Correlating Local Volumetric Tissue Strains with Global Lung Mechanics Measurements

The mechanics of breathing is a fascinating and vital process. The lung has complexities and subtle heterogeneities in structure across length scales that influence mechanics and function. This study establishes an experimental pipeline for capturing alveolar deformations during a respiratory cycle using synchrotron radiation micro-computed tomography (SR-micro-CT). Rodent lungs were mechanically ventilated and imaged at various time points during the respiratory cycle. Pressure-Volume (P-V) characteristics were recorded to capture any changes in overall lung mechanical behaviour during the experiment. A sequence of tomograms was collected from the lungs within the intact thoracic cavity. Digital volume correlation (DVC) was used to compute the three-dimensional strain field at the alveolar level from the time sequence of reconstructed tomograms. Regional differences in ventilation were highlighted during the respiratory cycle, relating the local strains within the lung tissue to the global ventilation measurements. Strains locally reached approximately 150% compared to the averaged regional deformations of approximately 80-100%. Redistribution of air within the lungs was observed during cycling. Regions which were relatively poorly ventilated (low deformations compared to its neighbouring region) were deforming more uniformly at later stages of the experiment (consistent with its neighbouring region). Such heterogenous phenomena are common in everyday breathing. In pathological lungs, some of these non-uniformities in deformation behaviour can become exaggerated, leading to poor function or further damage. The technique presented can help characterize the multiscale biomechanical nature of a given pathology to improve patient management strategies, considering both the local and global lung mechanics.

Functional lung imaging with synchrotron radiation: Methods and preclinical applications

Many lung disease processes are characterized by structural and functional heterogeneity that is not directly appreciable with traditional physiological measurements. Experimental methods and lung function modeling to study regional lung function are crucial for better understanding of disease mechanisms and for targeting treatment. Synchrotron radiation offers useful properties to this end: coherence, utilized in phase-contrast imaging, and high flux and a wide energy spectrum which allow the selection of very narrow energy bands of radiation, thus allowing imaging at very specific energies. K-edge subtraction imaging (KES) has thus been developed at synchrotrons for both human and small animal imaging. The unique properties of synchrotron radiation extend X-ray computed tomography (CT) capabilities to quantitatively assess lung morphology, and also to map regional lung ventilation, perfusion, inflammation and biomechanical properties, with microscopic spatial resolution. Four-dimensional imaging, allows the investigation of the dynamics of regional lung functional parameters simultaneously with structural deformation of the lung as a function of time. This review summarizes synchrotron radiation imaging methods and overviews examples of its application in the study of disease mechanisms in preclinical animal models, as well as the potential for clinical translation both through the knowledge gained using these techniques and transfer of imaging technology to laboratory X-ray sources.

MRI and CT lung biomarkers: Towards an in vivo understanding of lung biomechanics

BACKGROUND

The biomechanical properties of the lung are necessarily dependent on its structure and function, both of which are complex and change over time and space. This makes in vivo evaluation of lung biomechanics and a deep understanding of lung biomarkers, very challenging. In patients and animal models of lung disease, in vivo evaluations of lung structure and function are typically made at the mouth and include spirometry, multiple-breath gas washout tests and the forced oscillation technique. These techniques, and the biomarkers they provide, incorporate the properties of the whole organ system including the parenchyma, large and small airways, mouth, diaphragm and intercostal muscles. Unfortunately, these well-established measurements mask regional differences, limiting their ability to probe the lung's gross and micro-biomechanical properties which vary widely throughout the organ and its subcompartments. Pulmonary imaging has the advantage in providing regional, non-invasive measurements of healthy and diseased lung, in vivo. Here we summarize well-established and emerging lung imaging tools and biomarkers and how they may be used to generate lung biomechanical measurements.

METHODS

We review well-established and emerging lung anatomical, microstructural and functional imaging biomarkers generated using synchrotron x-ray tomographic-microscopy (SRXTM), micro-x-ray computed-tomography (micro-CT), clinical CT as well as magnetic resonance imaging (MRI).

FINDINGS

Pulmonary imaging provides measurements of lung structure, function and biomechanics with high spatial and temporal resolution. Imaging biomarkers that reflect the biomechanical properties of the lung are now being validated to provide a deeper understanding of the lung that cannot be achieved using measurements made at the mouth.

The Effectiveness of An Averaged Airway Model in Predicting the Airflow and Particle Transport Through the Airway

Background: In this study, we proposed an averaged airway model design based on four healthy subjects and numerically evaluated its effectiveness for predicting the airflow and particle transport through an airway. Methods: Direct-averaged models of the conducting airways of four subjects were restored by averaging the three-dimensional (3D) skeletons of four healthy airways, which were calculated using an inverse 3D thinning algorithm. We simulated the airflow and particle transport in the individual and the averaged airway models using computational fluid dynamics. Results: The bifurcation geometry differs even among healthy subjects, but the averaged model retains the typical geometrical characteristics of the airways. The Reynolds number of the averaged model varied within the range found in the individual subject models, and the averaged model had similar inspiratory flow characteristics as the individual subject models. The deposition fractions at almost all individual lobes ranged within the variation observed in the subjects, however, the deposition fraction was higher in only one lobe. The deposition distribution at the main bifurcation point differed among the healthy subjects, but the characteristics of the averaged model fell within the variation observed in the individual subject models. On the contrary, the deposition fraction of the averaged model was higher than that of the average of the individual subject models and deviated from the range observed in the subject models. Conclusion: These results indicate that the direct-averaged model may be useful for predicting the individual airflow and particle transport on a macroscopic scale.

Challenge in particle delivery to cells in a microfluidic device

Micro and nanotechnology can potentially revolutionize drug delivery systems. Novel microfluidic systems have been employed for the cell culture applications and drug delivery by micro and nanocarriers. Cells in the microchannels are under static and dynamic flow perfusion of culture media that provides nutrition and removes waste from the cells. This exerts hydrostatic and hydrodynamic forces on the cells. These forces can considerably affect the functions of the living cells. In this paper, we simulated the flow of air, culture medium, and the particle transport and deposition in the microchannels under different angles of connection inlet. It was found that the shear stress induced by the medium culture flow is not so high to damage the cells and that it is roughly uniform in the cell culture section (CCS). However, the local shear stresses in the other parts of the microchip differ by changing the angles of the connection inlet. The results showed that the particle deposition was a function of the particle size, the properties of the fluid, and the flow rate. At a lower air flow rate, both small and large particles deposited in the entrance region and none of them reached the CCS. Once the airflow rate increased, the drag of the flow could overcome the diffusion of the small particles and deliver them to the CCS so that more than 88% of the 100 nm and 98% of the 200 nm particles deposited in the CCS. However, larger particles with average diameters in micrometers could not reach the CCS by the airflow even at high flow rate. In contrast, our findings indicated that both small and large particles could be delivered to the CCS by liquid flow. Our experimental data confirm that microparticles (with diameters of 5 and 20 μm) suspended in a liquid can reach the CCS at a well-adjusted flow rate. Consequently, a liquid carrier is suggested to transport large particles through microchannels. As a powerful tool, these numerical simulations provide a nearly complete understanding of the flow field and particle patterns in microchips which can significantly lower the trial and error in the experiment tests and accordingly save researchers considerable cost and time for drug delivery to the cell in the microchip by micro/nanocarriers.

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.

3D patient-specific numerical modeling of the soft palate considering adhesion from the tongue

Collapse of the soft palate in the upper airway contributes to obstructive sleeping apnea (OSA). In this study, we investigate the influence of the adhesion from the tongue on the soft palate global response. This is achieved using a cohesive zone finite element approach. A traction-separation law is determined to describe the adhesion effect from the surface tension of the lining liquid between the soft palate and the tongue. According to pull-off experimental tests of human lining liquid from the oral surface of the soft palate, the corresponding cohesive properties, including the critical normal traction stress and the failure separation displacement, are obtained. The 3D patient-specific soft palate geometry is accounted for, based on one specific patient's computed tomography (CT) images. The calculation results show that influence of the adhesion from the tongue surface on the global response of the soft palate depends on the length ratio between the cohesive length and the soft palate length. When the length of the cohesive zone is smaller than half of the soft palate length, the adhesion's influence is negligible. When the adhesion length is larger than 70 percent of soft palate length, the adhesion force contributes to preventing the soft palate from collapsing towards to the pharynx wall, i.e. the closing pressure is more negative than in the no adhesion case. These results may provide useful information to the clinical treatment of OSA patients.

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.

Palatal implant surgery effectiveness in treatment of obstructive sleep apnea: A numerical method with 3D patient-specific geometries

Obstructive sleep apnea (OSA) affects a large percentage of the population and is increasingly recognized as a major global health problem. One surgical procedure for OSA is to implant polyethylene (PET) material into the soft palate, but its efficacy remains to be discussed. In this study, we provide input to this topic based on numerical simulations. Three 3 dimensional (3D) soft palate finite element models including mouth-close and mouth-open cases were created based on three patient-specific computed tomography (CT) images. A simplified material modeling approach with the Neo-Hookean material model was applied, and nonlinear geometry was accounted for. Young's modulus for the implant material was obtained from uniaxial tests, and the PET implant pillars were inserted to the 3D soft palate model. With the finite element model, we designed different surgical schemes and investigated their efficacy with respect to avoiding the soft palate collapse. Several pillar schemes were tested, including different placement directions, different placement positions, different settings for the radius and the array parameters of the implant pillars, and different Young's moduli for the pillars. Based on our simulation results, the longitudinal-direction implant surgery improved the stiffness of the soft palate to a small degree, and implanting in the transverse direction was evaluated to be a good choice for improving the existing surgical scheme. In addition, the Young's modulus of the polyethylene material implants has an influence on the reinforcement efficacy of the soft palate.

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).

A simplified parametrised model for lung microstructures capable of mimicking realistic geometrical and mechanical properties

The respiratory zone of mammalian lungs contains several millions of so-called alveoli. The geometrical and mechanical properties of this microstructure are crucial for respiration and influence the macroscopic behaviour of the entire organ in health and disease. Hence, if computational models are sought to gain more insight into lung behaviour, predict lung states in certain scenarios or suggest better treatment options in early stages of respiratory dysfunction, an adequate representation of this microstructure is essential. However, investigating the real alveolar architecture requires complex medical-imaging methods and would be computationally extremely expensive. Even worse, there is currently no way of obtaining the real patient-specific microstructure in vivo. Hence, we present a fast and easy to compute parametrised model of lung microstructures based on tetrakaidecahedra which can represent both geometrical and mechanical properties of the parenchyma. We show that gas transport pathways and stress and strain distributions are comparable to real alveolar microstructures and even capable of capturing variations present in biology. The created parametrised lung microstructure models can be utilized in finite element simulations to study, e.g., alveolar flow phenomena, particle deposition, or alveolar stresses and strains during mechanical ventilation. Due to the simpler geometry of the parametrised microgeometries compared to imaging-based microstructures, remarkable savings in CPU time can be achieved. We show that our model requires a minimum of 10% of the computational time for computing the same strain state in structural mechanics simulations compared to imaging-based alveolar microstructures.

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.

The role of anisotropic expansion for pulmonary acinar aerosol deposition

Lung deformations at the local pulmonary acinar scale are intrinsically anisotropic. Despite progress in imaging modalities, the true heterogeneous nature of acinar expansion during breathing remains controversial, where our understanding of inhaled aerosol deposition still widely emanates from studies under self-similar, isotropic wall motions. Building on recent 3D models of multi-generation acinar networks, we explore in numerical simulations how different hypothesized scenarios of anisotropic expansion influence deposition outcomes of inhaled aerosols in the acinar depths. While the broader range of particles acknowledged to reach the acinar region (dp=0.005-5.0μm) are largely unaffected by the details of anisotropic expansion under tidal breathing, our results suggest nevertheless that anisotropy modulates the deposition sites and fractions for a narrow band of sub-micron particles (dp~0.5-0.75μm), where the fate of aerosols is greatly intertwined with local convective flows. Our findings underscore how intrinsic aerosol motion (i.e. diffusion, sedimentation) undermines the role of anisotropic wall expansion that is often attributed in determining aerosol mixing and acinar deposition.

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.