Particle and inhalation exposure in human and monkey computational airway models

Computational fluid-particle dynamics modeling of ultrafine to coarse particles deposition in the human respiratory system, down to the terminal bronchiole

BACKGROUND AND OBJECTIVES

Suspended respirable airborne particles are associated with human health risks and especially particles within the range of ultrafine (< 0.1 μm) or fine (< 2.5 μm) have a high possibility of penetrating the lung region, which is concerned to be closely related to the bronchial or alveoli tissue dosimetry. Nature complex structure of the respiratory system requires much effort to explore and comprehend the flow and the inhaled particle dynamics for precise health risk assessment. Therefore, this study applied the computational fluid-particle dynamics (CFPD) method to elucidate the deposition characteristics of ultrafine-to-coarse particles in the human respiratory tract from nostrils to the 16^th generation of terminal bronchi.

METHODS

The realistic bronchi up to the 8^th generation are precisely and perfectly generated from computed tomography (CT) images, and an artificial model compensates for the 9^th-16^th bronchioles. Herein, the steady airflow is simulated at constant breathing flow rates of 7.5, 15, and 30 L/min, reproducing human resting-intense activity. Then, trajectories of the particle size ranging from 0.002 - 10 μm are tracked using a discrete phase model.

RESULTS

Here, we report reliable results of airflow patterns and particle deposition efficiency in the human respiratory system validated against experimental data. The individual-related focal point of ultrafine and fine particles deposition rates was actualized at the 8^th generation; whilst the hot-spot of the deposited coarse particles was found in the 6^th generation. Lobar deposition characterizes the dominance of coarse particles deposited in the right lower lobe, whereas the left upper-lower and right lower lobes simultaneously occupy high deposition rates for ultrafine particles. Finally, the results indicate a higher deposition in the right lung compared to its counterpart.

CONCLUSIONS

From the results, the developed realistic human respiratory system down to the terminal bronchiole in this study, in coupling with the CFPD method, delivers the accurate prediction of a wide range of particles in terms of particle dosimetry and visualization of site-specific in the consecutive respiratory system. In addition, the series of CFPD analyses and their results are to offer in-depth information on particle behavior in human bronchioles, which may benefit health risk assessment or drug delivery studies.

Large eddy simulation of droplet transport and deposition in the human respiratory tract to evaluate inhalation risk

As evidenced by the worldwide pandemic, respiratory infectious diseases and their airborne transmission must be studied to safeguard public health. This study focuses on the emission and transport of speech-generated droplets, which can pose risk of infection depending on the loudness of the speech, its duration and the initial angle of exhalation. We have numerically investigated the transport of these droplets into the human respiratory tract by way of a natural breathing cycle in order to predict the infection probability of three strains of SARS-CoV-2 on a person who is listening at a one-meter distance. Numerical methods were used to set the boundary conditions of the speaking and breathing models and large eddy simulation (LES) was used for the unsteady simulation of approximately 10 breathing cycles. Four different mouth angles when speaking were contrasted to evaluate real conditions of human communication and the possibility of infection. Breathed virions were counted using two different approaches: the breathing zone of influence and direction deposition on the tissue. Our results show that infection probability drastically changes based on the mouth angle and the breathing zone of influence overpredicts the inhalation risk in all cases. We conclude that to portray real conditions, the probability of infection should be based on direct tissue deposition results to avoid overprediction and that several mouth angles must be considered in future analyses.

Detailed comparative analysis of environmental microparticle deposition characteristics between human and monkey nasal cavities using a surface mapping technique

Inhaled particulate matter is associated with nasal diseases such as allergic rhinitis, rhinosinusitis and neural disorders. Its health risks on humans are usually evaluated by measurements on monkeys as they share close phylogenetic relationship. However, the reliability of cross-species toxicological extrapolation is in doubt due to physiological and anatomical variations, which greatly undermine the reliability of these expensive human surrogate models. This study numerically investigated in-depth microparticle transport and deposition characteristics on human and monkey (Macaca fuscata) nasal cavities that were reconstructed from CT-images. Deposition characteristics of 1-30μm particles were investigated under resting and active breathing conditions. Similar trends were observed for total deposition efficiencies and a single correlation using Stokes Number was fitted for both species and both breathing conditions, which is convenient for monkey-human extrapolation. Regional deposition patterns were carefully compared using the surface mapping technique. Deposition patterns of low, medium and high inertial particles, classified based on their total deposition efficiencies, were further analyzed in the 3D view and the mapped 2D view, which allows locating particle depositions on specific nasal regions. According to the particle intensity contours and regional deposition profiles, the major differences were observed at the vestibule and the floor of the nasal cavity, where higher deposition intensities of medium and high inertial particles were shown in the monkey case than the human case. Comparisons of airflow streamlines indicated that the cross-species variations of microparticle deposition patterns are mainly contributed by two factors. First, the more oblique directions of monkey nostrils result in a sharper airflow turn in the vestibule region. Second, the monkey's relatively narrower nasal valves lead to higher impaction of medium and high inertial particles on the nasal cavity floor. The methods and findings in this study would contribute to an improved cross-species toxicological extrapolation between human and monkey nasal cavities.

Detailed comparison of anatomy and airflow dynamics in human and cynomolgus monkey nasal cavity

Nonhuman primates are occasionally used as laboratory models for sophisticated medical research as they bear the closest resemblance to humans in morphometry and physiological functions. A range of nonhuman primate species have been employed in the inhalation toxicity, nasal drug delivery and respiratory viral infection studies, and they provided valuable insight to disease pathogenesis while other laboratory animals such as rodents cannot recapitulate due to the lesser degree of similarity in metabolism, anatomy and cellular response to that of humans. It is anticipated that nonhuman primate models of respiratory diseases will continue to be instrumental for translating biomedical research for improvement of human health, and the confidence in laboratory data extrapolation between species will play a pivotal role. From the morphometry and flow dynamics point of view, this study performed a detailed comparative analysis between human and a cynomolgus monkey nasal airway, with intention to provide high-fidelity qualitative and quantitative linkage between the two species for more effective laboratory data extrapolation. The study revealed that cynomolgus monkey could be a good human surrogate in nasal inhalation studies; however, care should be given for interspecies data extrapolation as subtle differences in anatomy and airflow dynamics were present between the two species.

Ga-68 EDTA aerosols in evaluation of inhaled-particle deposition and clearance of obstructive pulmonary diseases: A pilot prospective study compared with Galligas

PURPOSE

68-gallium (Ga-68) ethylenediaminetetraacetic acid (EDTA) aerosols and Galligas were compared in evaluation of inhaled-particle deposition and clearance in volunteers with or without obstructive pulmonary diseases.

METHODS

Nonsmoking healthy volunteers, healthy smokers, asthma patients and patients with chronic obstructive pulmonary disease (COPD) were recruited to undergo the dynamic lung ventilation positron emission tomography/computerized tomography (PET/CT) scans within two consecutive days. The inhaled particles were Ga-68-labelled carbon nanoparticles (Galligas, 30-60 nm in size) and Ga-68-labelled EDTA aerosols (1-2 μm in size), respectively. The volunteers' lung function parameters were measured for comparison.

RESULTS

Central deposition and inhomogeneity of both tracers were negatively correlated with lung function parameters, including the ratio of forced expiratory volume at 1 second to forced vital capacity (FEV₁ /FVC). The central or hilum deposition of Galligas, but not 68-gallium (Ga-68) EDTA, was negatively correlated with the maximal expiratory flow at 25%, 50% and 75% of the forced vital capacity. Compared with Galligas, Ga-68 EDTA aerosols were more concentrated in the central region in all groups except for the healthy nonsmokers. Ventilation inhomogeneity was more evident when using Ga-68 EDTA aerosols, especially in patients with COPD and asthma patients. In the healthy smokers, the central region accumulated more Ga-68 EDTA at 30 minutes after inhalation than immediately after inhalation. Ga-68 EDTA cleared faster in lungs than Galligas.

CONCLUSIONS

Both Galligas and Ga-68 EDTA aerosols can be used for PET/CT lung ventilation scan. However, Ga-68 EDTA aerosols showed more advantages in diagnosis and evaluation of obstructive airway diseases by revealing the inhaled-particle deposition and clearance.

TRANSITION AND LAMINAR FLOWS IN A REALISTIC GEOMETRY OF HUMAN UPPER AIRWAY

In this study, a realistic respiratory airway model extending from oral to the end of the trachea including all the key details of the passage was produced. A series of CT scan images were used to generate the topological data of airway cross-sections that were used to generate the computational model, as well as the three-dimensional (3D) printed model of the passage for experimental study. The airflow velocity field and pressure drop in the airway for different breathing rates of 5, 7.5, 10, and 12.5[Formula: see text]L/min were investigated numerically (by laminar and transition models) and experimentally. The velocity distributions, pressure variation, and streamlines along the oral–trachea airway model were studied. The maximum pressure drop was shown to occur in the narrowest part of the larynx region. It was also concluded that the laryngeal jet could significantly influence the airway flow patterns in the trachea. A comparison between the numerical results and experimental data showed that the transition [Formula: see text]–kl–[Formula: see text] model can give better predictions of pressure losses, especially for flow rates higher than 10[Formula: see text]L/min. The simulation results for the velocity profiles in the trachea were also compared with the available particle image velocimetry (PIV) data and earlier simulations. Despite inter-personal variability and difference in the flow regime, the qualitative agreement was found.

Biomedical and biophysical limits to mathematical modeling of pulmonary system mechanics: a scoping review on aerosol and drug delivery

Undoubtedly, the construction of the biomechanical geometry systems with the help of computer tomography (CT) and magnetic resonance imaging (MRI) has made a significant advancement in studying in vitro numerical models as accurately as possible. However, some simplifying assumptions in the computational studies of the respiratory system have caused errors and deviations from the in vivo actual state. The most important of these hypotheses is how to generate volume from the point cloud exported from CT or MRI images, not paying attention to the wall thickness and its effect in computational fluid dynamic method, statistical logic of aerosol trap in software; and most importantly, the viscoelastic effect of respiratory tract wall in living tissue pointed in the fluid–structure interaction method. So that applying the viscoelastic dynamic mesh effect in the form of the moving deforming mesh can be very effective in achieving more appropriate response quality. Also, changing the volume fraction of the pulmonary extracellular matrix constituents leads to changes in elastic modulus (storage modulus) and the viscous modulus (loss modulus) of lung tissue. Therefore, in the biomedical computational methods where the model wall is considered flexible, the viscoelastic properties of the texture must be considered correctly.

The impact of nasal adhesions on airflow and mucosal cooling - A computational fluid dynamics analysis

Nasal adhesions are a known postoperative complication following surgical procedures for nasal airway obstruction (NAO); and are a common cause of surgical failure, with patients often reporting significant NAO, despite relatively minor adhesion size. Division of such nasal adhesions often provides much greater relief than anticipated, based on the minimal reduction in cross-sectional area associated with the adhesion. The available literature regarding nasal adhesions provides little evidence examining their quantitative and qualitative effects on nasal airflow using objective measures. This study examined the impact of nasal adhesions at various anatomical sites on nasal airflow and mucosal cooling using computational fluid dynamics (CFD). A high-resolution CT scan of the paranasal sinuses of a 25-year-old, healthy female patient was segmented to create a three-dimensional nasal airway model. Virtual nasal adhesions of 2.5 mm diameter were added to various locations within the nasal cavity, representing common sites seen following NAO surgery. A series of models with single adhesions were created. CFD analysis was performed on each model and compared with a baseline no-adhesion model, comparing airflow and heat and mass transfer. The nasal adhesions resulted in no significant change in bulk airflow patterns through the nasal cavity. However, significant changes were observed in local airflow and mucosal cooling around and immediately downstream to the nasal adhesions. These were most evident with anterior nasal adhesions at the internal valve and anterior inferior turbinate. Postoperative nasal adhesions create local airflow disruption, resulting in reduced local mucosal cooling on critical surfaces, explaining the exaggerated perception of nasal obstruction. In particular, anteriorly located adhesions created greater disruption to local airflow and mucosal cooling, explaining their associated greater subjective sensation of obstruction.

Computational modelling of nasal respiratory flow

CFD has emerged as a promising diagnostic tool for clinical trials, with tremendous potential. However, for real clinical applications to be useful, overall statistical findings from large population samples (e.g., multiple cases and models) are needed. Fully resolved solutions are not a priority, but rather rapid solutions with fast turn-around times are desired. This leads to the issue of what are the minimum modelling criteria for achieving adequate accuracy in respiratory flows for large-scale clinical applications, with a view to rapid turnaround times. This study simulated a highly-resolved solution using the large eddy simulation (LES) method as a reference case for comparison with lower resolution models that included larger time steps and no turbulence modelling. Differences in solutions were quantified by pressure loss, flow resistance, unsteadiness, turbulence intensity, and hysteresis effects from multiple cycles. The results demonstrated that sufficient accuracy could be achieved with lower resolution models if the mean flow was considered. Furthermore, to achieve an established transient result unaffected by the initial start-up quiescent effects, the results need to be taken from at least the second respiration cycle. It was also found that the exhalation phase exhibited strong turbulence. The results are expected to provide guidance for future modelling efforts for clinical and engineering applications requiring large numbers of cases using simplified modelling approaches.

Investigation of airflow at different activity conditions in a realistic model of human upper respiratory tract

In the present study, the turbulent flows inside a realistic model of the upper respiratory tract were investigated numerically and experimentally. The airway model included the geometrical details of the oral cavity to the end of the trachea that was based on a series of CT-scan images. The topological data of the respiratory tract were used for generating the computational model as well as the 3D-printed model that was used in the experimental pressure drop measurement. Different airflow rates of 30, 45, and 60 L/min, which correspond to the light, semi-light, and heavy activity breathing conditions, were investigated numerically using turbulence and transition models, as well as experimentally. Simulation results for airflow properties, including velocity vectors, pressure drops, streamlines, eddy viscosity, and turbulent kinetic energy contours in the oral-trachea airway model, were presented. The simulated pressure drop was compared with the experimental data, and reasonable agreement was found. The obtained results showed that the maximum pressure drop occurs in the narrowest part of the larynx region. A comparison between the numerical results and experimental data showed that the transition (γ-Reθ) SST model predicts higher pressure losses, especially at higher breathing rates. Formations of the secondary flows in the oropharynx and trachea regions were also observed. In addition, the simulation results showed that in the trachea region, the secondary flow structures dissipated faster for the flow rate of 60 L/min compared to the lower breathing rates of 30 and 45 L/min.

In silico approaches to respiratory nasal flows: A review

The engineering discipline of in silico fluid dynamics delivers quantitative information on airflow behaviour in the nasal regions with unprecedented detail, often beyond the reach of traditional experiments. The ability to provide visualisation and analysis of flow properties such as velocity and pressure fields, as well as wall shear stress, dynamically during the respiratory cycle may give significant insight to clinicians. Yet, there remains ongoing challenges to advance the state-of-the-art further, including for example the lack of comprehensive CFD modelling on varied cohorts of patients. The present article embodies a review of previous and current in silico approaches to simulating nasal airflows. The review discusses specific modelling techniques required to accommodate physiologically- and clinically-relevant findings. It also provides a critical summary of the reported results in the literature followed by an outlook on the challenges and topics anticipated to drive research into the future.