Flow and Particle Deposition Patterns in a Realistic Human Double Bifurcation Airway Model

Study on the coal dust deposition fraction and site in the upper respiratory tract under different particle sizes and labor intensities

In order to study the dust exposure amount and coal dust deposition rule of coal miners under different labor intensity in coal mine environment, an airflow-particle two-phase coupling calculation model of human upper respiratory tract was established based on Euler-Lagrange framework, and the airflow field in the upper respiratory tract and the characteristics of coal dust deposition were simulated and studied. By comparing the experimental data, the relative error of simulation is in the range of 1.5 %-11.2 %. The results showed that the total deposition fraction of 1 μm dust was the smallest (0.61-1.20 %), and was relatively less affected by respiratory intensity, and the overall distribution was uniform. When the dust particle size increased to 7.07 μm, the total dust deposition fraction in the nasal cavity, pharynx and larynx was in the range of 11.10 %-20.91 %, and increased with the respiratory intensity. When the dust particle size was large, the dust particles of 20 μm and 80 μm were mostly concentrated in the front of the nasal cavity, and the deposition amount of 80 μm dust was about 99.52 %. It was found that with the increase of dust particle size or the increase of labor intensity, the possibility of dust being transported into lungs became smaller. The fitting function of 7.07 μm dust escape rate and labor intensity was obtained, for example, Y_7.07μm = 91.73-0.22n (n is labor intensity), and the escape rate of dust with 7.07 μm particle size was up to 88.90 %. Most of them escape from the upper respiratory tract and enter the lungs, which provides theoretical guidance for quantifying the accumulated dust exposure amount in the lungs and monitoring respiratory dust concentration.

Impact of sleep posture and breathing pattern on soft palate flutter and pharynx vibration in a pediatric airway using fluid-structure interaction

Snoring is a common condition in the general population, and the management of snoring requires a better understanding of its mechanism through a fluid-structure interaction (FSI) perspective. Despite the recent popularity of numerical FSI techniques, outstanding challenges are accurately predicting airway deformation and its vibration during snoring due to complex airway morphology. In addition, there still needs to be more understanding of snoring inhibition when lying on the side, and the possible effect of airflow rates, as well as nose or mouth-nose breathing, on snoring remains to be investigated. In this study, an FSI method verified against in vitro models was introduced to predict upper airway deformation and vibration. The technique was applied to predict airway aerodynamics, soft palate flutter, and airway vibration in four sleep postures (supine, left/right lying, and sitting positions) and four breathing patterns (mouth-nose, nose, mouth, and unilateral nose breathing). It was found that, at given elastic properties of soft tissues, the evaluated flutter frequency of 19.8 Hz in inspiration was in good agreement with the reported frequency of snoring sound in literature. Reduction in flutter and vibrations due to the mouth-nose airflow proportion changes were also noticed when having side-lying and sitting positions. Breathing through the mouth results in larger airway deformation than breathing through the nose or mouth-nose. These results collectively demonstrate the potential of FSI for studying the physics of airway vibration and clarify to some degree the reason for snoring inhibition during sleep postures and breathing patterns.

Particle inhalability of a standing mannequin with large airways in a ventilated room

This study presents a series of numerical simulations for airflow field and particle dispersion and deposition around a mannequin inside a ventilated room. A 3-D airway system of a volunteer subject with a large respiratory system was reconstructed from the nostril inlet to the end of the tracheobronchial tree 4th generation and was integrated into a standing mannequin at the center of a room. The room ventilation system supplied air through a diffuser and expelled air via a damper in three modes. The airflow field was first evaluated by solving the governing equations and the k-ω SST transitional turbulence model using the Ansys-Fluent software. Then spherical particles with various diameters were released into the room, and their trajectories were evaluated using the Lagrangian approach. Aspiration fraction and particle deposition for inhalation flow rates of 15 and 30 L/min were analyzed using a modified discrete random walk (DRW) stochastic model using a user-defined function (UDF) coupled to the Ansys-Fluent discrete phase model. For the first ventilation mode, a recirculation flow region formed behind the mannequin that led the airflow streamlines to the breathing zone. A recirculation flow formed in front of the face for the second ventilation mode that led the airflow streamlines out of the mannequin breathing zone. For the third mode, however, there was no strong recirculation flow zone around the mannequin. Simulation results showed that the aspiration fraction in the first ventilation mode was higher than the other modes. In addition, the regional deposition rates and deposition patterns of particles inside the respiratory system were presented for each region. Accordingly, most large particles were trapped in the nasal passage; however, some large particles penetrated deeper into the airway due to the large airway size. For the higher breathing rate, the percentage of large escaped particles from the lobe branches dropped by a factor of 7 compared to the lower breathing rate.

Structural airway imaging metrics are differentially associated with persistent chronic bronchitis

BACKGROUND

Chronic bronchitis (CB) is strongly associated with cigarette smoking, but not all smokers develop CB. We aimed to evaluate whether measures of structural airway disease on CT are differentially associated with CB.

METHODS

In smokers between ages 45 and 80 years, and with Global Initiative for Obstructive Lung Disease stages 0-4, CB was defined by the classic definition. Airway disease on CT was quantified by (i) wall area percent (WA%) of segmental airways; (ii) Pi10, the square root of the wall area of a hypothetical airway with 10 mm internal perimeter; (iii) total airway count (TAC) and (iv) airway fractal dimension (AFD), a measure of the complex branching pattern and remodelling of airways. CB was also assessed at the 5-year follow-up visit.

MEASUREMENTS AND MAIN RESULTS

Of 8917 participants, 1734 (19.4%) had CB at baseline. Airway measures were significantly worse in those with CB compared with those without CB: WA% 54.5 (8.8) versus 49.8 (8.3); Pi10 2.58 (0.67) versus 2.28 (0.59) mm; TAC 156.7 (81.6) versus 177.8 (91.1); AFD 1.477 (0.091) versus 1.497 (0.092) (all p<0.001). On follow-up of 5517 participants at 5 years, 399 (7.2%) had persistent CB. With adjustment for between-visits changes in smoking status and lung function, greater WA% and Pi10 were associated with significantly associated with persistent CB, adjusted OR per SD change 1.75, 95% CI 1.56 to 1.97; p<0.001 and 1.66, 95% CI 1.42 to 1.86; p<0.001, respectively. Higher AFD and TAC were associated with significantly lower odds of persistent CB, adjusted OR per SD change 0.76, 95% CI 0.67 to 0.86; p<0.001 and 0.69, 95% CI 0.60 to 0.80; p<0.001, respectively.

CONCLUSIONS

Higher baseline AFD and TAC are associated with a lower risk of persistent CB, irrespective of changes in smoking status, suggesting preserved airway structure can confer a reserve against CB.

Simulation of nanoparticle transport and adsorption in a microfluidic lung-on-a-chip device

The effect of air-borne nanoparticles (NPs) on human health is an active area of research, with clinical relevance evidenced by the current COVID-19 pandemic. As in vitro models for such studies, lung-on-a-chip (LOAC) devices can represent key physical and physiological aspects of alveolar tissues. However, widespread adoption of the LOAC device for NP testing has been hampered by low intra-laboratory and inter-laboratory reproducibility. To complement ongoing experimental work, we carried out finite-element simulations of the deposition of NPs on the epithelial layer of a well-established LOAC design. We solved the Navier-Stokes equations for the fluid flow in a three-dimensional domain and studied the particle transport using Eulerian advection-diffusion for fine NPs and Lagrangian particle tracking for coarse NPs. Using Langmuir and Frumkin kinetics for surface adsorption and desorption, we investigated NP adsorption under different exercise and breath-holding patterns. Conditions mimicking physical exercise, through changes in air-flow volume and breathing frequency, enhance particle deposition. Puff profiles typical of smoking, with breath-holding between inhalation and exhalation, also increase particle deposition per breathing cycle. Lagrangian particle tracking shows Brownian motion and gravitational settling to be two key factors, which may cooperate or compete with each other for different particle sizes. Comparisons are made with experimental data where possible and they show qualitative and semi-quantitative agreement. These results suggest that computer simulations can potentially inform and accelerate the design and application of LOAC devices for analyzing particulate- and microbe-alveolar interactions.

Computational modeling of airway and pulmonary vascular structure and function: development of a "lung physiome"

Computational models of lung structure and function necessarily span multiple spatial and temporal scales, i.e., dynamic molecular interactions give rise to whole organ function, and the link between these scales cannot be fully understood if only molecular or organ-level function is considered. Here, we review progress in constructing multiscale finite element models of lung structure and function that are aimed at providing a computational framework for bridging the spatial scales from molecular to whole organ. These include structural models of the intact lung, embedded models of the pulmonary airways that couple to model lung tissue, and models of the pulmonary vasculature that account for distinct structural differences at the extra- and intra-acinar levels. Biophysically based functional models for tissue deformation, pulmonary blood flow, and airway bronchoconstriction are also described. The development of these advanced multiscale models has led to a better understanding of complex physiological mechanisms that govern regional lung perfusion and emergent heterogeneity during bronchoconstriction.

Image-based modeling of lung structure and function

The current state-of-the-art in image-based modeling allows derivation of patient-specific models of the lung, lobes, airways, and pulmonary vascular trees. The application of traditional engineering analyses of fluid and structural mechanics to image-based subject-specific models has the potential to provide new insight into structure-function relationships in the individual via functional interpretation that complements imaging and experimental studies. Three major issues that are encountered in studies of airflow through the bronchial airways are the representation of airway geometry, the imposition of physiological boundary conditions, and the treatment of turbulence. Here we review some efforts to resolve each of these issues, with particular focus on image-based models that have been developed to simulate airflow from the mouth to the terminal bronchiole, and subjected to physiologically meaningful boundary conditions via image registration and soft-tissue mechanics models.

The lung physiome: merging imaging-based measures with predictive computational models

Global measurements of the lung provided by standard pulmonary function tests do not give insight into the regional basis of lung function and lung disease. Advances in imaging methodologies, computer technologies, and subject-specific simulations are creating new opportunities to study structure-function relationships in the lung through multidisciplinary research. The digital Human Lung Atlas is an image-based resource compiled from male and female subjects spanning several decades of age. The Atlas comprises both structural and functional measures, and includes computational models derived to match individual subjects for personalized prediction of function. The computational models in the Atlas form part of the Lung Physiome project, which is an international effort to develop integrative models of lung function at all levels of biological organization. The computational models provide mechanistic interpretation of imaging measures; the Atlas provides structural data on which to base model geometry, and functional data against which to test hypotheses. The example of simulating airflow on a subject-specific basis is considered. Methods for deriving multiscale models of the airway geometry for individual subjects in the Atlas are outlined, and methods for modeling turbulent flows in the airway are reviewed.

Flow analyses in the lower airways: patient-specific model and boundary conditions

Computational fluid dynamics (CFD) is increasingly applied in the respiratory domain. The ability to simulate the flow through a bifurcating tubular system has increased the insight into the internal flow dynamics and the particular characteristics of respiratory flows such as secondary motions and inertial effects. The next step in the evolution is to apply the technique to patient-specific cases, in order to provide more information about pathological airways. This study presents a patient-specific approach where both the geometry and the boundary conditions (BC) are based on individual imaging methods using computed tomography (CT). The internal flow distribution of a 73-year-old female suffering from chronic obstructive pulmonary disease (COPD) is assessed. The validation is performed through the comparison of lung ventilation with gamma scintigraphy. The results show that in order to obtain agreement within the accuracy limits of the gamma scintigraphy scan, both the patient-specific geometry and the BC (driving pressure) play a crucial role. A minimal invasive test (CT scan) supplied enough information to perform an accurate CFD analysis. In the end it was possible to capture the pathological features of the respiratory system using the imaging and computational fluid dynamics techniques. This brings the introduction of this new technique in the clinical practice one step closer.