Particle Deposition in Large-Scale Human Tracheobronchial Airways Predicted...

[Full title: Particle Deposition in Large-Scale Human Tracheobronchial Airways Predicted by Single-Path Modelling]

Citation: Ou, C.; Hang, J.; Hua, J.; Li, Y.; Deng, Q.; Zhao, B.; Ling, H Particle Deposition in Large-Scale Human Tracheobronchial Airways Predicted by Single-Path Modelling Int. J. Environ. Res. Public Health 2023 , 20 , 4583. https://doi.org/10.3390/ ijerph 20054583 Academic Editors: Lili Wang, Xingru Li and Minjuan Huang Received: 29 January 2023 Revised: 22 February 2023 Accepted: 28 February 2023 Published: 4 March 2023 Copyright: © 2023 by the authors Licensee MDPI, Basel, Switzerland This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/) International Journal of Environmental Research and Public Health Article Particle Deposition in Large-Scale Human Tracheobronchial Airways Predicted by Single-Path Modelling Cuiyun Ou 1 , Jian Hang 1 , Jiajia Hua 2 , Yuguo Li 3 , Qihong Deng 4 , Bo Zhao 2 and Hong Ling 1, * 1 School of Atmospheric Sciences, Sun Yat-sen University, and Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China 2 China Meteorological Administration Xiong’an Atmospheric Boundary Layer Key Laboratory, Baoding 071800, China 3 Department of Mechanical Engineering, The University of Hong Kong, Hong Kong, China 4 College of Public Health, Zhengzhou University, Zhengzhou 450001, China * Correspondence: lingh 23@mail.sysu.edu.cn Abstract: The health effects of particles are directly related to their deposition patterns (deposition site and amount) in human airways. However, estimating the particle trajectory in a large-scale human lung airway model is still a challenge. In this work, a truncated single-path, large-scale human airway model (G 3–G 10) with a stochastically coupled boundary method were employed to investigate the particle trajectory and the roles of their deposition mechanisms. The deposition patterns of particles with diameters ( dp ) of 1–10 µ m are investigated under various inlet Reynolds numbers ( Re = 100–2000). Inertial impaction, gravitational sedimentation, and combined mechanism were considered. With the increasing airway generations, the deposition of smaller particles ( dp < 4 µ m) increased due to gravitational sedimentation, while that of larger particles decreased due to inertial impaction. The obtained formulas of Stokes number and Re can predict the deposition efficiency due to the combined mechanism in the present model, and the prediction can be used to assess the dose-effect of atmospheric aerosols on the human body. Diseases in deeper generations are mainly attributed to the deposition of smaller particles under lower inhalation rates, while diseases at the proximal generations mainly result from the deposition of larger particles under higher inhalation rates Keywords: computational fluid dynamics (CFD); large-scale human airway; single-path modelling; particle deposition mechanism; deposition prediction 1. Introduction The adverse health effects of particulate matters (PM) have been studied in epidemiological and toxicological works [ 1 – 5 ]. Previous studies have found that an increased concentration of PM was associated with an increased risk of all-cause and cardiopulmonary mortality, and a higher risk has been found for smaller-sized particles [ 6 – 11 ]. There are four main exposure routes for the human body to particulate pollutants, namely ingestion, inhalation, injection and dermal contact [ 12 ]. Among them, inhalation is considered as significant because fine airborne particles (particle diameter dp ≤ 2.5 µ m) can penetrate into the human lung. The novel coronavirus disease (COVID-19) also spreads to susceptible individuals by inhalable droplets or droplet nuclei [ 13 ]; therefore, it is essential to know the health effects of particles. Particle deposition patterns, both the rate and the site, in human airways can provide a linkage between PM concentration and health effects [ 14 – 16 ]. Computational fluid dynamic (CFD) modelling has been demonstrated to be an effective approach to predict particle deposition in human airways. Most CFD studies have mainly focused on the small-scale modelling of tracheobronchial airways by using single-, doubleand triple-bifurcations [ 17 – 19 ]. They have found that simple bifurcation models can sufficiently interpret the general airflow structure and particle transport in the certain Int. J. Environ. Res. Public Health 2023 , 20 , 4583. https://doi.org/10.3390/ijerph 20054583 https://www.mdpi.com/journal/ijerph

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 2 of 15 airways; however, Comer et al. [ 20 ] indicated that a simple bifurcation model might not provide a realistic view of particle deposition in the respiratory airways [ 21 ]. Accordingly, several studies have developed new methods to model the entire tracheobronchial (TB) airway by integrating the local bifurcation model in parallel and series [ 22 – 25 ]. However, it is difficult to determine the boundary conditions between local bifurcations, which have significant influences on the prediction of particle transport and the deposition in human airways [ 23 ]. Therefore, airway modelling with more generations is necessary The direct simulation of particle deposition in multiple-generation airways has become the focus of researchers. van Ertbruggen et al. [ 26 ] proposed a three-dimensional airway model in the upper-7-generation airways, and then Gemci et al. [ 27 ] studied the airflow in the upper-17-generation airways. However, these models are difficult for wide application due to their consideration of large numbers of branches. Walters and Luke [ 28 ] suggested a stochastic coupling method to simulate the airflow pattern and particle deposition in the truncated larger-scale airways (G 4 to G 12). They found that branches could be partly removed without affecting the deposition simulation in the remaining branches [ 29 ]. Tian et al. [ 30 ] developed a stochastic individual path (SIP) approach to characterize the transport and deposition of particles throughout the TB region. More applications and validations are warranted for further investigation In the present work, we utilized the aforementioned stochastic coupling method. The main idea of this method is to map the static pressure at the resolved airway flow paths to the outlets of the unresolved flow paths in the randomly truncated single-path airway model. Particles with sizes ranging from 1 to 10 µ m were computed. Two deposition mechanisms, namely inertial impaction and gravitational sedimentation, were taken into account. Meanwhile, the effects of the particle sizes and airflow rates were evaluated 2. Methods 2.1. Truncated Airway Model As obtained in our previous work [ 31 ], we could see that the airflow in lower generations (lower than G 11) was too small, and the dominant particle-deposition mechanism should be Brownian diffusion, which was not the concern of the present work. Hence, in the present study, an eight-generation airway model (G 3–G 10) was established based on the data from Weibel [ 32 ]. Then, we randomly truncated one of the two daughter branches from the middle downstream of each generation. The truncated branch was named an unresolved outlet, while one of the two branches in G 10 remained as the resolved outlet. Hence, there were seven unresolved outlets and one resolved outlet in the model as shown in Figure 1 . Parameters of the model are summed in the accompanying table, including the diameter of generation Z ( D Z ), the length of generation Z ( L Z ) and the angle between 2 daughter branches (2 Θ ). The curvature radius of the ring connecting the daughter branches ( R Z ), which can be calculated by 2 D Z +1 , is also illustrated in Figure 1 . This method is an example of the single flow-path model originally proposed by Tian et al. [ 30 ]. To make the truncated geometry faithfully represent the effect of the removed airways, a special boundary condition must be set to the unresolved outlets, which is described in detail in Section 2.3 . 2.2. Governing Equations The airflow in the airway model is assumed to be Stokes flow (laminar, steady, and incompressible with constant density and viscosity). The non-dimensional governing equations for the continuity and momentum are given as follows: ∇· u = 0 (1) u ·∇ u = −∇ p + 1 Re ∇ 2 u (2) where u and p are the dimensionless velocity vector and the pressure, respectively u = u ∗ / U , p = p ∗ / ρ U 2 , where u ∗ and p ∗ are the local velocity and the pressure, respectively;

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 3 of 15 U and D are the mean velocity at the inlet and the diameter of the inlet, respectively; ρ is the air density (1.225 kg/m 3 ). The air Reynolds number ( Re ) is defined as Re = ρ UD / µ , where µ is the air dynamic viscosity (1.7894 × 10 − 5 kg/(ms)) Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 3 of 15 Figure 1. Truncated eight-generation model based on Weibel lung morphology [32]. 2.2. Governing Equations The airflow in the airway model is assumed to be Stokes flow (laminar, steady, and incompressible with constant density and viscosity). The non-dimensional governing equations for the continuity and momentum are given as follows: ∇ ∙ {?} = 0 (1) {?} ∙ ∇{?} = −∇{?} + 1 {?}{?} ∇ {?} (2) where u and p are the dimensionless velocity vector and the pressure, respectively. {?} = {?} ∗ {?} ⁄ , {?} = {?} ∗ ρ{?} ⁄ , where {?} ∗ and {?} ∗ are the local velocity and the pressure, respectively; U and D are the mean velocity at the inlet and the diameter of the inlet, respectively; ρ is the air density (1.225 kg/m 3 ). The air Reynolds number ( Re ) is defined as {?}{?} = {?}{?}{?} {?} ⁄ , where {?} is the air dynamic viscosity (1.7894 × 10 − 5 kg/(ms)). The particle phase is assumed to be sufficiently diluted, and subsequently, the transport can be modeled by a one-way coupled Lagrangian method in which the particle motion is described as: d {?} ∗ d{?} = {?} {?} ∗ − {?} ∗ + {?} (3) where {?} ∗ and {?} ∗ are the particle displacement and the velocity, respectively; {?} is the gravitational acceleration; {?} represents the drag force coefficient per unit particle mass and is defined as: {?} = {?} {?} {?}{?} {?} (24{?} {?}) ⁄ (4) where {?} and {?} are the Cunningham correction factor and the drag coefficient [33], respectively. {?}{?} is the particle Reynolds number defined as: {?}{?} = {?} {?} ∗ − {?} ∗ {?}{?} {?} ⁄ (5) {?} is the Stokes number and is defined as: {?} = {?} {?}{?} {?}{?} (18{?}{?}) ⁄ (6) Figure 1. Truncated eight-generation model based on Weibel lung morphology [ 32 ]. The particle phase is assumed to be sufficiently diluted, and subsequently, the transport can be modeled by a one-way coupled Lagrangian method in which the particle motion is described as: d 2 x ∗ p d t 2 = f D u ∗ − u ∗ p + g (3) where x ∗ p and u ∗ p are the particle displacement and the velocity, respectively; g is the gravitational acceleration; f D represents the drag force coefficient per unit particle mass and is defined as: f D = C c C D Re p U / ( 24 S tk D ) (4) where C c and C D are the Cunningham correction factor and the drag coefficient [ 33 ], respectively Re p is the particle Reynolds number defined as: Re p = ρ u ∗ − u ∗ p dp / µ (5) S tk is the Stokes number and is defined as: S tk = ρ p dp 2 UC c / ( 18 µ D ) (6) where dp and ρ p are the particle diameter and the density (2000 kg/m 3 ), respectively 2.3. Boundary Conditions Uniform velocity profiles were employed at the inlet. According to the range of oral and/or nasal inhalation airflow rate ( Q = 5 to 60 L/min, corresponding to the human activities from sleep to intense exercise), the inlet velocity U ranged from 0.26 m/s to 5.2 m/s. The corresponding Reynolds numbers ranged from 100 to 2000 in a step of 100. A no-slip velocity boundary condition was applied at all walls. Pressure outlet was employed at the resolved outlets, which was assumed to be equal as atmospheric pressure. We

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 4 of 15 mapped the plane-averaged static pressure at the middle cross-section of the resolved flow-path to the unresolved outlet in the same generation [ 29 ], as the airflow was assumed to be symmetrical at the bifurcation. At the beginning of the calculation for the airflow field, we first set the pressure outlet to all the outlets, including both the resolved and unresolved outlets. After iterating several times, the mapping processes were conducted, and then the calculation was completed until the convergence was obtained After obtaining the flow field, particles were released uniformly at the inlet. In this work, we injected particles with the same size for each injection. The trap condition was applied to all the walls with the assumption that particles were trapped as soon as they touched the wall surfaces. In other words, the deposition occurred when particles reached the airway wall, and the escape condition was employed for the outlets 2.4. Particle Deposition Efficiency The regional deposition of particles in human airways can be quantified in terms of the deposition efficiency ( DE ), which is defined as: DE = No d / No e (7) where No d is the number of particles deposited; No e is the number of particles entered. We calculated DE separately to distinguish the contribution of different deposition mechanisms in the airway model DE i corresponds to the inertial impaction DE g is the gravitational sedimentation DE i+g is contributed by the combined mechanism. Their relationship is described by the following equation [ 34 ]: DE i + g = DE i + DE g − DE i · DE g (8) 2.5. Numerical Methods The governing equations were numerically solved using a finite volume method by a commercial CFD software ANSYS FLUENT (Version 12.0.16; ANSYS Inc., Canonsburg, PA, USA). The pre-processing tool UNIGRAPHICS (Version 5.0; Siemens PLM Software, Berlin, Germany) was used to build the truncated model, and a structured mesh was generated by GAMBIT (Version 2.4.6; ANSYS Inc., Canonsburg, PA, USA). Refined boundary layers were employed near walls and junctions where the velocity gradient may be larger. The number of final cells for the symmetric geometry was 1,570,000. This number was determined by the grid-independent tests for the flow field solution and particle deposition efficiency, according to the grid convergence index applied by Longest and Vinchurkar [ 35 ]. Particles were injected similar to a Dirac function after the steady airflow was achieved, and the disperse phase model was used to calculate the particle trajectories The status of particles (e.g., trap, escape, etc.) was recorded once the track of particles was fulfilled. This process was repeated. Then, the average number of particles trapped in specific regions was obtained, which was used for the statistics of the deposition efficiency. Different numbers of released particles (i.e., 10,000, 20,000, 50,000, and 100,000) were considered for the independence tests. No significant change in deposition was found when the number of released particles was larger than 50,000. This method was validated by comparing it with others’ work in our previous paper [ 35 ]. 3. Results Figure 2 shows the contours of the axial velocity magnitude and the secondary velocity vector at different cross-sections of different generations when Re = 1000 (i.e., an inhalation rate of 30 L/min). Both axial and secondary air velocities were higher near the inner wall of each divider and decreased with each generation. It also shows that the secondary airflow was axisymmetric on the plane of the corresponding divider, and the air moved from the outer side to the inner side in each axis of symmetry, and then produced two symmetric vortices near the axis. The vortices also weakened with each generation.

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 5 of 15 Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 5 of 15 3. Results Figure 2 shows the contours of the axial velocity magnitude and the secondary velocity vector at different cross-sections of different generations when Re = 1000 (i.e., an inhalation rate of 30 L/min). Both axial and secondary air velocities were higher near the inner wall of each divider and decreased with each generation. It also shows that the secondary airflow was axisymmetric on the plane of the corresponding divider, and the air moved from the outer side to the inner side in each axis of symmetry, and then produced two symmetric vortices near the axis. The vortices also weakened with each generation. Figure 2. Axial air velocity and secondary airflow vector at various sections ( Re = 1000). Figure 3 presents the mean axial air velocity at the beginning section of each generation and the corresponding deposition patterns of 10- μ m particles (i.e., dp = 10 μ m. Particles with other sizes are named following the same way) when {?}{?} = 1000 . The axial air flow rate decreased when the generation number increased, changing from 2.6 m/s (G 3) to 0.56 m/s (G 9). The corresponding deposition patterns show that 10- μ m particles were deposited more in the upper parts than that of in the lower parts. Figure 2 B illustrates there are hot spots of deposition for 10- μ m particles, particularly at the carinal ridge of each bifurcation. Figure 4 displays the particle deposition in the whole G 3–G 10 truncated airway model due to different mechanisms for particles of different sizes (1, 3, and 5 μ m). In general, the deposition related to inertial impaction ( DE i ) increases with rising Re , while that ascribed to the gravitational sedimentation ( DE g ) decreases. For 1- μ m particles, the inertial impaction has few influences on the deposition throughout the entire considered Re range. Gravitational sedimentation is dominant when Re ranges from 100 to 2000, and the deposition decreases with the increasing Re . For 3- μ m particles, a U-shape curve of DE i+g can be found. When {?}{?} 800 , the gravitational sedimentation was dominant. Hence, DE i+g decreases rapidly with the increasing Re . When {?}{?} > 1200 , inertial impaction prevails, and the DE i+g curve rapidly rises up with the rising Re . When 800 {?}{?} 1200 , the two mechanisms balance each other, and thus the DE i+g curve is stable at a low level. For 5- μ m particles, the knee point (when gravitational sedimentation balances inertial impaction) of the DE i+g curve moves from higher Re to lower Re, compared to that of 3- μ m particles. Figure 2. Axial air velocity and secondary airflow vector at various sections ( Re = 1000) Figure 3 presents the mean axial air velocity at the beginning section of each generation and the corresponding deposition patterns of 10- µ m particles (i.e., dp = 10 µ m. Particles with other sizes are named following the same way) when Re = 1000. The axial air flow rate decreased when the generation number increased, changing from 2.6 m/s (G 3) to 0.56 m/s (G 9). The corresponding deposition patterns show that 10- µ m particles were deposited more in the upper parts than that of in the lower parts. Figure 2 B illustrates there are hot spots of deposition for 10- µ m particles, particularly at the carinal ridge of each bifurcation Figure 4 displays the particle deposition in the whole G 3–G 10 truncated airway model due to different mechanisms for particles of different sizes (1, 3, and 5 µ m). In general, the deposition related to inertial impaction ( DE i ) increases with rising Re , while that ascribed to the gravitational sedimentation ( DE g ) decreases. For 1- µ m particles, the inertial impaction has few influences on the deposition throughout the entire considered Re range. Gravitational sedimentation is dominant when Re ranges from 100 to 2000, and the deposition decreases with the increasing Re . For 3- µ m particles, a U-shape curve of DE i+g can be found. When Re < 800, the gravitational sedimentation was dominant. Hence, DE i+g decreases rapidly with the increasing Re . When Re > 1200, inertial impaction prevails, and the DE i+g curve rapidly rises up with the rising Re . When 800 < Re < 1200, the two mechanisms balance each other, and thus the DE i+g curve is stable at a low level. For 5- µ m particles, the knee point (when gravitational sedimentation balances inertial impaction) of the DE i+g curve moves from higher Re to lower Re, compared to that of 3- µ m particles Figure 5 shows curves of DE with particle diameter ( dp ) ranging from 1 to 10 µ m in logarithmic coordinates. Generally, depositions dominated by all three mechanisms are enhanced with the increasing dp . For the inertial impaction, DE i increases slowly below a critical dp (e.g., dp = 5 µ m when Re = 100), while it increases rapidly after crossing the knee point. The larger the Re is, the smaller the critical dp is. For the gravitational sedimentation, DE g increases with the rising dp (i.e., a linear increase in the logarithmic coordinates). For the combined mechanism, DE i+g curves display similar trends as DE g curves when Re ≤ 500. When Re > 500, the variation in DE i+g is more similar to that of DE i , only the critical dp is slightly changed for each Re .

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 6 of 15 Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 6 of 15 ( a ) ( b ) Figure 3. Airflow structure and corresponding particle deposition pattern, including the deposition fraction and deposition sites ( Re = 1000, dp = 10 μ m). ( a ) axial velocity and deposition fraction; ( b ) deposition pattern of dp = 10 μ m. ( a ) dp = 1 μ m Figure 3. Airflow structure and corresponding particle deposition pattern, including the deposition fraction and deposition sites ( Re = 1000, dp = 10 µ m). ( a ) axial velocity and deposition fraction; ( b ) deposition pattern of dp = 10 µ m Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 6 of 15 ( a ) ( b ) Figure 3. Airflow structure and corresponding particle deposition pattern, including the deposition fraction and deposition sites ( Re = 1000, dp = 10 μ m). ( a ) axial velocity and deposition fraction; ( b ) deposition pattern of dp = 10 μ m. ( a ) dp = 1 μ m Figure 4. Cont .

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 7 of 15 Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 7 of 15 ( b ) dp = 3 μ m ( c ) dp = 5 μ m Figure 4. Particle deposition efficiency ( DE ) due to different mechanisms for three different particle diameters ( dp ). Figure 5 shows curves of DE with particle diameter ( dp ) ranging from 1 to 10 μ m in logarithmic coordinates. Generally, depositions dominated by all three mechanisms are enhanced with the increasing dp . For the inertial impaction, DE i increases slowly below a critical dp (e.g., {?}{?} = 5 μm when {?}{?} = 100 ), while it increases rapidly after crossing the knee point. The larger the Re is, the smaller the critical dp is. For the gravitational sedimentation, DE g increases with the rising dp (i.e., a linear increase in the logarithmic coordinates). For the combined mechanism, DE i+g curves display similar trends as DE g curves when {?}{?} ≤ 500 . When {?}{?} > 500 , the variation in DE i+g is more similar to that of DE i , only the critical dp is slightly changed for each Re . ( a ) Figure 4. Particle deposition efficiency ( DE ) due to different mechanisms for three different particle diameters ( dp ) Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 8 of 15 ( a ) ( b ) ( c ) Figure 5. Particle deposition efficiency ( DE ) vs. particle diameter ( dp ) due to different mechanisms, including ( a ) inertial impaction; ( b ) gravitational sedimentation; ( c ) combined mechanism Figure 6 depicts the deposition fraction in different generations (G 3–G 10). The deposition fraction due to inertial impaction increases with generations for small particles ( {?}{?} ≤ 4 μm ), while it decreases for larger particles ( {?}{?} ≥ 5 μm ). The deposition fraction due to gravitational sedimentation increases slowly from G 3 to G 9 (except the relatively high deposition fraction at G 4) for particles from 1 to 7 μ m but decreases if dp = 8–10 μ m. Figure 5. Cont .

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 8 of 15 Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 8 of 15 ( a ) ( b ) ( c ) Figure 5. Particle deposition efficiency ( DE ) vs. particle diameter ( dp ) due to different mechanisms, including ( a ) inertial impaction; ( b ) gravitational sedimentation; ( c ) combined mechanism Figure 6 depicts the deposition fraction in different generations (G 3–G 10). The deposition fraction due to inertial impaction increases with generations for small particles ( {?}{?} ≤ 4 μm ), while it decreases for larger particles ( {?}{?} ≥ 5 μm ). The deposition fraction due to gravitational sedimentation increases slowly from G 3 to G 9 (except the relatively high deposition fraction at G 4) for particles from 1 to 7 μ m but decreases if dp = 8–10 μ m. Figure 5. Particle deposition efficiency ( DE ) vs. particle diameter ( dp ) due to different mechanisms, including ( a ) inertial impaction; ( b ) gravitational sedimentation; ( c ) combined mechanism Figure 6 depicts the deposition fraction in different generations (G 3–G 10). The deposition fraction due to inertial impaction increases with generations for small particles ( dp ≤ 4 µ m), while it decreases for larger particles ( dp ≥ 5 µ m). The deposition fraction due to gravitational sedimentation increases slowly from G 3 to G 9 (except the relatively high deposition fraction at G 4) for particles from 1 to 7 µ m but decreases if dp = 8–10 µ m Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 8 of 15 ( b ) ( c ) Figure 5. Particle deposition efficiency ( DE ) vs. particle diameter ( dp ) due to different mechanisms, including ( a ) inertial impaction; ( b ) gravitational sedimentation; ( c ) combined mechanism Figure 6 depicts the deposition fraction in different generations (G 3–G 10). The deposition fraction due to inertial impaction increases with generations for small particles ( {?}{?} ≤ 4 μm ), while it decreases for larger particles ( {?}{?} ≥ 5 μm ). The deposition fraction due to gravitational sedimentation increases slowly from G 3 to G 9 (except the relatively high deposition fraction at G 4) for particles from 1 to 7 μ m but decreases if dp = 8–10 μ m. ( a ) ( b ) Figure 6. Generational deposition fractions for particles in multi-size ( dp = 1 μ m to 10 μ m) when Re = 1000, ( a ) inertial impaction; ( b ) gravitational sedimentation Figure 6. Generational deposition fractions for particles in multi-size ( dp = 1 µ m to 10 µ m) when Re = 1000, ( a ) inertial impaction; ( b ) gravitational sedimentation Figure 7 shows the “S” shape variation curves of DE in terms of non-dimensional parameters. Figure 7 a presents the growth curves of DE i that associated with S tk only. The fitting curve can be represented by a logistic regression as following: DE i ( % ) = 100 1 − 1/ ( 19.65 S tk 2.0996 + 1 ) (9) Figure 7 b illustrates the variation in DE g in terms of S tk / Re 1.6 . It can be fitted to a logistic regression curve such as the following: DE g ( % ) = 100 1 − 1/ ( 3.25 × 10 8 ( S tk / Re 1.6 ) 1.517 + 1 ) (10)

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 9 of 15 Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 9 of 15 Figure 7 shows the “S” shape variation curves of DE in terms of non-dimensional parameters. Figure 7 a presents the growth curves of DE i that associated with S tk only. The fitting curve can be represented by a logistic regression as following: {?}{?} (%) = 100(1 − 1 ( ⁄ 19.65{?} + 1)) (9) Figure 7 b illustrates the variation in DE g in terms of {?} / {?}{?} . It can be fitted to a logistic regression curve such as the following: {?}{?} (%) = 100 1 − 1 (3.25 × 10 ({?} {?}{?} ) ⁄ ⁄ + 1) (10) ( a ) Inertial impaction ( b ) Gravitational sedimentation Figure 7. Variation of particle deposition efficiency ( DE ) in terms of non-dimensional parameters: ( a ) S tk for inertial impaction and ( b ) S tk / Re 1.6 for gravitational sedimentation. 4. Discussion This study predicted the fate of particles and the deposition mechanisms in a largescale TB airway model. We first assessed the influences of the airflow and respiratory parameters on the particle deposition patterns. We then analyzed the factors that affected the DE by using three different mechanisms: inertial impaction, gravitational sedimentation, and a combined mechanism. Finally, we investigated the generational deposition fractions resulting from the mechanisms. We established that the trajectories of particles in the human lung can be influenced by the airflow structure, dp , and deposition mechanism. Furthermore, we obtained a formula in terms of S tk and Re to predict the deposition efficiency ( DE ) due to the combined mechanism in the present model. Finally, we gave a suggestion for assessing the dose-effect of surrounding particulate pollutants on the human body. We used the method of three-dimensional (3-D) truncated single-path airway model. Gemci et al. [27] modelled the whole computer-tomography-based TB (G 0–G 16) airway including 1453 bronchi, although they truncated nearly 99% of the bronchi to simplify the computation. They determined the outflow boundary conditions at outlets by using mass distributions, leading the predicting pressure drop approximately 30% lower than the experimental results obtained by Hyatt and Wilcox [36]. To reduce the computational expenses effectively, numerous researchers have put their efforts into developing new methods to model the airflow and the particle transport patterns in the large-scale human lung airway model. Walters and Luke [28] proposed a method that stochastically mapped the plane-averaged static pressure at the middle cross-section from one of the resolved flow paths to the unresolved outlets in the same generation in a symmetrical G 4–G 12 airway model. The stochastically individual path (SIP) airway modelling approach suggested by Tian et al. [30] considered the influence of different lung lobes on boundary conditions at upper generations of TB airways. Different ventilation modes in the five lung lobes were Figure 7. Variation of particle deposition efficiency ( DE ) in terms of non-dimensional parameters: ( a ) S tk for inertial impaction and ( b ) S tk / Re 1.6 for gravitational sedimentation 4. Discussion This study predicted the fate of particles and the deposition mechanisms in a largescale TB airway model. We first assessed the influences of the airflow and respiratory parameters on the particle deposition patterns. We then analyzed the factors that affected the DE by using three different mechanisms: inertial impaction, gravitational sedimentation, and a combined mechanism. Finally, we investigated the generational deposition fractions resulting from the mechanisms. We established that the trajectories of particles in the human lung can be influenced by the airflow structure, dp , and deposition mechanism Furthermore, we obtained a formula in terms of S tk and Re to predict the deposition efficiency ( DE ) due to the combined mechanism in the present model. Finally, we gave a suggestion for assessing the dose-effect of surrounding particulate pollutants on the human body We used the method of three-dimensional (3-D) truncated single-path airway model Gemci et al. [ 27 ] modelled the whole computer-tomography-based TB (G 0–G 16) airway including 1453 bronchi, although they truncated nearly 99% of the bronchi to simplify the computation. They determined the outflow boundary conditions at outlets by using mass distributions, leading the predicting pressure drop approximately 30% lower than the experimental results obtained by Hyatt and Wilcox [ 36 ]. To reduce the computational expenses effectively, numerous researchers have put their efforts into developing new methods to model the airflow and the particle transport patterns in the large-scale human lung airway model. Walters and Luke [ 28 ] proposed a method that stochastically mapped the plane-averaged static pressure at the middle cross-section from one of the resolved flow paths to the unresolved outlets in the same generation in a symmetrical G 4–G 12 airway model. The stochastically individual path (SIP) airway modelling approach suggested by Tian et al. [ 30 ] considered the influence of different lung lobes on boundary conditions at upper generations of TB airways. Different ventilation modes in the five lung lobes were used to determine the flow rate at the corresponding bronchi for outflow boundary conditions. Tena et al. [ 37 ] proposed a method that applied velocities mapped to the unresolved sections from similar locations in fully resolved airway sections to model the lung airway up to G 16. The approaches proposed by Walters and Luke [ 28 ] and Tian et al. [ 30 ] were efficient in simplifying the airway geometry for modelling particle deposition patterns in large-scale human lung airways; however, there are insufficient applications or validations. Our previous study [ 35 ] adopted the method of Tian et al. [ 28 ] to compare the particle deposition patterns of human and ferret TB airway models (G 0–G 16) We found that this method could effectively complete such a large-scale airway simulation; therefore, the present work adopted this method to simulate the particle deposition in the

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 10 of 15 single-path G 3–G 10 airway model that can be considered as one conducting airway of the five lung lobes We found that the airflow structures can influence the particle deposition, which is consistent with previous studies [ 38 – 40 ]. For example, Comer et al. [ 41 ] found that the three-dimensional airflow influenced the particle deposition in a double-bifurcation model. Zhang and Kleinstreuer [ 42 ] indicated that the particle deposition pattern in a triplebifurcation model was highly influenced by the airflow structure. Rahimi-Gorji et al. [ 43 ] found that the maximum velocity change occurred in the larynx might cause more particle deposition at this location. Likewise, in our G 3–G 10 model we found the particle deposition pattern was greatly affected by the airflow structure. When the airflow velocity decreases with the rising generation, the generational deposition fraction of smaller particles (e.g., dp = 4 µ m) increases, however, that of larger particles (e.g., dp = 10 µ m) decreases We found that inertial impaction and gravitational sedimentation play important roles in the deposition fraction at different sites of the airway. Inertial impaction leads the increased deposition fraction of smaller particles (1 µ m ≤ dp ≤ 4 µ m) and the decreased deposition fraction of larger particles (5 µ m ≤ dp ≤ 10 µ m) with the rising generation Gravitational sedimentation dominates the deposition of smaller particles but not for larger ones. Thus, the gravity only works on smaller particles depositing at lower generations Therefore, larger particles mainly deposit at upper generations due to inertial impaction, while smaller ones are mainly concentrated at lower generations due to gravitational sedimentation. We also noticed an exception to this rule that a relatively high deposition fraction was observed at G 4 due to gravity for particles in the entire size range. This exception may be influenced by the inlet boundary [ 21 ] as well as the morphological configuration (i.e., the bifurcation orientation [ 20 ], the lengths and the diameters of each generation tube [ 44 ], and the geometry parameters of each bifurcation divider [ 45 ]). The model applied in this work can be easily extended to the whole TB airway; hence, it can be used to precisely derive the particle trajectories in the lung and can help better understand the health effect of PM The strength of our study lies in the fact that we separated the roles for different deposition mechanisms in large-scale human lung TB airways. The deposition mechanism is one of the physical factors that determine the particle deposition pattern in the human respiratory tract [ 46 ]. Particle deposition in the lung airways mainly occurs owing to three mechanisms: inertial impaction, gravitational sedimentation, and Brownian diffusion [ 47 ]. Previous studies have focused on either the deposition mechanisms in simple bifurcation airway models [ 48 – 50 ] or the deposition patterns in complicated airway models but without precise descriptions for the mechanisms [ 51 , 52 ], which cannot shed light on the deposition mechanisms in depth. There are studies of gravitational deposition that have concentrated on the small bronchial airways and alveolar region. Sedimentation has been found playing the dominant role in these regions [ 53 , 54 ]. In our work, we investigated the variation in DE i , DE g , and DE i+g in terms of Re and dp . The sedimentation was found dominating the inertial deposition under small- Re conditions, while the inertial impaction become dominant over sedimentation with the increasing Re DE i , DE g , and DE i+g were found increasing with dp The findings are in good agreement with many previous studies, indicating that particle deposition in the TB region increases with dp rising from 1 to 10 µ m [ 47 , 55 , 56 ]. Another merit of this study is that we established the relationships between nondimensional parameters and DE i , DE g , and DE i+g , respectively. First, we found that DE i varied with S tk in a consistent regression curve for the entire range of inlet Re and obtained a statistical regression equation of S tk . Then, we found that DE g varied with S tk / Re 1.6 in a consistent way. Finally, the contribution of the combined mechanism increased with the particle size; however, it decreased with Re when Re ≤ 500 and increased with Re when Re ≥ 1000. Therefore, we failed to find a single parameter to characterize DE i+g throughout the whole range of Re and the particle size. Nevertheless, we can describe it separately using S tk (when Re ≤ 500) or S tk / Re 1.6 (when Re ≥ 1000). Cheng et al. [ 57 ] found that oropharyngeal deposition (equivalent to DE i+g in the present work) had a good

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 11 of 15 correlation with S tk only, while the Reynolds number had no effect on deposition for flow rates ranging from 15 to 60 L/min under which conditions of inertial impaction dominated over gravity. Kleinstreuer et al. [ 48 ] suggested that DE could be correlated only with S tk under normal breathing conditions (i.e., Q ≥ 15 L/min), while under slow breathing conditions (i.e., Q ≤ 7.5 L/min), sedimentation became increasingly dominant over inertial impaction. This finding agrees well with our present results. Chen et al. [ 17 ] suggested that deposition in the airway model of generations G 11–G 14 were determined by S tk and the sedimentation parameter (which is also a function of airflow rate) for particles in the size range of 1–7 µ m. This conclusion also agrees well with our present results The third merit of our study is that we probably have provided a formula in terms of Stokes number ( S tk ) and Reynolds number ( Re ) to predict the deposition efficiency due to the combined mechanism ( DE i+g ). We first obtained the equation of DE i+g by replacing the formula (9) and formula (10) into formula (8) as follows: DE i + g ( % ) = 100 − 100/ (( 19.65 S 2.0996 tk + 1 )( 3.25 × 10 8 ( S tk / Re 1.6 ) 1.517 + 1 )) (11) Then, we calculated the value of DE i+g for particles in the size range of 1–10 µ m under different Re , namely values predicted by equation. Finally, the comparisons of the values between those predicted by equation and those simulated numerically are displayed in Figure 8 . From this figure, we can find that the accuracy of the prediction decreases with Re but increases with particle size, which can be obviously found by calculating the relative error between prediction value and simulation value, as shown in Table 1 . When Re = 500, the relative error ranges from − 71% to 36% with a mean value of − 21%, and the maximum value occurs when dp = 1 µ m, while the minimum occurs when dp = 7 µ m. When Re = 2000, the relative error ranges from − 79% to − 6% with a mean value of − 40%, and the maximum value occurs when dp = 1 µ m, while the minimum occurs when dp = 10 µ m. In other words, the prediction of the deposition efficiency in the large-scale airway model can assess the dose-effect of surrounding particles to human body. The larger the particles are, or the lower the inhalation rates are, the better the accuracy of the prediction can be Int. J. Environ. Res. Public Health 2023 , 20 , x FOR PEER REVIEW 12 of 15 model, connecting the research of inhalable aerosols with in vivo deposition, exposure, and morbidity of respiratory diseases. Using the current work as a crucial tool, a series of the interdisciplinary studies of air pollution and public health are planned. ( a ) ( b ) ( c ) ( d ) Figure 8. Comparisons of the values of DE i+g between the prediction and simulation, including ( a ) Re = 500; ( b ) Re = 1000; ( c ) Re = 1500; ( d ) Re = 2000. Table 1. The relative error between the prediction value and simulation value. dp ( μ m) Re = 500 Re = 1000 Re = 1500 Re = 2000 1 (0.72) (0.72) (0.79) (0.79) 2 (0.36) (0.42) (0.49) (0.43) 3 (0.01) (0.19) (0.41) (0.55) 4 0.22 (0.20) (0.50) (0.61) 5 0.37 (0.26) (0.55) (0.57) 6 0.19 (0.23) (0.45) (0.42) 7 0.02 (0.22) (0.34) (0.28) 8 (0.06) (0.20) (0.23) (0.16) 9 (0.10) (0.14) (0.13) (0.10) 10 (0.10) (0.10) (0.07) (0.06) Mean relative error (0.21) (0.27) (0.40) (0.40) Note: Values in parenthesis are negative values, which means underprediction. 5. Conclusions This work uses a truncated single-path large-scale human airway model (G 3–G 10) and a stochastically coupled boundary method for investigating the particle trajectory and the roles of deposition mechanisms in the human lung airway. Particle size ranging from 1 to 10 μ m, various inlet Reynolds numbers ( Re = 100–2000), and 3 deposition mechanisms (inertial deposition, gravitational sedimentation, and combined mechanism) are considered in this research. We found the truncated single-path lung model efficiently estimates the fate of particles in the lung. Particle size and inhalation rate were found to be key factors that determine the deposition site, amount, and mechanism. A formula is obtained in this work Figure 8. Comparisons of the values of DE i+g between the prediction and simulation, including ( a ) Re = 500; ( b ) Re = 1000; ( c ) Re = 1500; ( d ) Re = 2000.

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 12 of 15 Table 1. The relative error between the prediction value and simulation value dp ( µ m) Re = 500 Re = 1000 Re = 1500 Re = 2000 1 (0.72) (0.72) (0.79) (0.79) 2 (0.36) (0.42) (0.49) (0.43) 3 (0.01) (0.19) (0.41) (0.55) 4 0.22 (0.20) (0.50) (0.61) 5 0.37 (0.26) (0.55) (0.57) 6 0.19 (0.23) (0.45) (0.42) 7 0.02 (0.22) (0.34) (0.28) 8 (0.06) (0.20) (0.23) (0.16) 9 (0.10) (0.14) (0.13) (0.10) 10 (0.10) (0.10) (0.07) (0.06) Mean relative error (0.21) (0.27) (0.40) (0.40) Note: Values in parenthesis are negative values, which means underprediction Our model can be further improved, as it is a symmetrical airway model, while a realistic human airway should be asymmetrical, and asymmetry also affects the particle deposition [ 58 , 59 ]. The asymmetric geometries, such as different lobes, will lead to heterogeneous ventilation in different lung regions due to the nature of each region volume. A realistic, fluctuating inlet velocity profile should be considered in the future work because the inlet airflow patterns also affect the particle deposition [ 60 – 63 ]. The orientation of the gravity may play an important role in the deposition due to sedimentation, while the present work only takes into account the orientation perpendicular to the inlet of the model Moreover, the particle release pattern also affects the particle deposition pattern [ 44 , 64 ], which certainly warrants further study Furthermore, this work is essential for a better understanding of the correlations between human health and characteristics of aerosols (e.g., concentration distribution, dispersion, physical and chemical composition, toxicology, etc.). In our future work, we plan to couple the current work with an atmospheric dispersion model and an exposure model, connecting the research of inhalable aerosols with in vivo deposition, exposure, and morbidity of respiratory diseases. Using the current work as a crucial tool, a series of the interdisciplinary studies of air pollution and public health are planned 5. Conclusions This work uses a truncated single-path large-scale human airway model (G 3–G 10) and a stochastically coupled boundary method for investigating the particle trajectory and the roles of deposition mechanisms in the human lung airway. Particle size ranging from 1 to 10 µ m, various inlet Reynolds numbers ( Re = 100–2000), and 3 deposition mechanisms (inertial deposition, gravitational sedimentation, and combined mechanism) are considered in this research. We found the truncated single-path lung model efficiently estimates the fate of particles in the lung. Particle size and inhalation rate were found to be key factors that determine the deposition site, amount, and mechanism. A formula is obtained in this work for predicting the deposition efficiency due to the combined mechanism ( DE i+g ), using non-dimensional parameters, Stokes number ( S tk ), and Reynolds number ( Re ). Our findings indicate that the prediction of the deposition efficiency in the large-scale airway model can assess the dose-effect of surrounding particulate pollutants to human body. The larger the particles are, and/or the lower the inhalation rates are, the better accuracy the prediction can be. Our results imply that diseases that occur in deeper generations are mainly due to the deposition of smaller particles under lower inhalation rates. Meanwhile, diseases that occur in the proximal generations mainly result from the deposition of larger particles under higher inhalation rates Author Contributions: Conceptualization, C.O. and Q.D.; data duration, C.O.; writing—original draft preparation, C.O. and H.L.; writing—review and editing, C.O., H.L., J.H. (Jian Hang), Y.L., J.H (Jiajia Hua), and B.Z.; All authors have read and agreed to the published version of the manuscript.

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 13 of 15 Funding: The research was funded by the National Natural Science Foundation of China (No 42005069, No. 42175094, No. 42175095) and the R&D Plan Project in Key Fields of Guangdong Province (212020012620600004). The support from the Innovation Group Project of the Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (No. 311020001) is also gratefully acknowledged Data Availability Statement: The data presented in this study are available in the figures of this paper Conflicts of Interest: The authors declare no conflict of interest References 1 Wyzga, R.; Rohr, A.C. Long-term particulate matter exposure: Attributing health effects to individual PM components J. Air Waste Manag. Assoc 2015 , 65 , 523–543. [ CrossRef ] [ PubMed ] 2 Dockery, D. Health effects of particulate air pollution Ann. Epidemiol 2009 , 19 , 257–263. [ CrossRef ] 3 Xue, T.; Zhang, Q. Associating ambient exposure to fine particles and human fertility rates in China Environ. Pollut 2018 , 235 , 497–504. [ CrossRef ] 4 Achilleos, S.; Kioumourtzoglou, M.; Wu, C.; Schwartz, J.; Koutrakis, P.; Papatheodorou, S. Acute effects of fine particulate matter constituents on mortality: A systematic review and meta-regression analysis Environ. Int 2017 , 109 , 89–100. [ CrossRef ] 5 Neuberger, M.; Schimek, M.; Horak, F.; Moshammer, H.; Kundi, M.; Frischer, T.; Gomiscek, B.; Puxbaum, H.; Hauck, H. Acute effects of particulate matter on respiratory diseases, symptoms and functions: Epidemiological results of the Austrian Project on Health Effects of Particulate Matter (AUPHEP) Atmos. Environ 2004 , 38 , 3971–3981. [ CrossRef ] 6 An, Z.; Jin, Y.; Li, J.; Li, W.; Wu, W. Impact of Particulate air pollution on cardiovascular health Curr. Allergy Asthma Rep 2018 , 18 , 15. [ CrossRef ] [ PubMed ] 7 Scungio, M.; Stabile, L.; Rizza, V.; Pacitto, A.; Russi, A.; Buonanno, G. Lung cancer risk assessment due to traffic-generated particles exposure in urban street canyons: A numerical modelling approach Sci. Total Environ 2018 , 631–632 , 1109–1116 [ CrossRef ] 8 Pope, C., III; Burnett, R.T.; Thun, M.J.; Calle, E.E.; Krewski, D.; Ito, K.; Thurston, G.D. Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution JAMA 2002 , 287 , 1132–1141. [ CrossRef ] 9 Gao, P.; Hu, J.; Song, J.; Chen, X.; Ou, C.; Wang, H.; Sha, C.; Hang, J.; Xing, B. Inhalation bioaccessibility of polycyclic aromatic hydrocarbons in heavy PM 2.5 pollution days: Implications for public health risk assessment in northern China Environ. Pollut 2019 , 255 , 113296. [ CrossRef ] [ PubMed ] 10 Al-Hemoud, A.; Gasana, J.; Al-Dabbous, A.; Alajeel, A.; Al-Shatti, A.; Behbehani, W.; Malak, M. Exposure levels of air pollution (PM 2.5 ) and associated health risk in Kuwait Environ. Res 2019 , 179 , 108730. [ CrossRef ] 11 Prata, J.; da Costa, J.; Lopes, I.; Duarte, A.; Rocha-Santos, T. Environmental exposure to microplastics: An overview on possible human health effects Sci. Total Environ 2020 , 702 , 134455. [ CrossRef ] [ PubMed ] 12 Thompson, J. Airborne Particulate Matter: Human Exposure and Health Effects J. Occup. Environ. Med 2018 , 60 , 392–423 [ CrossRef ] [ PubMed ] 13 Ou, C.; Hu, S.; Luo, K.; Yang, H.; Hang, J.; Cheng, P.; Hai, Z.; Xiao, S.; Qian, H.; Xiao, S.; et al. Insufficient ventilation led to a probable long-range airborne transmission of SARS-CoV-2 on two buses Build. Environ 2022 , 207 , 108414 14 Fern á ndez Tena, A.; Casan Clar à , P. Deposition of inhaled particles in the lungs Arch. Bronconeumol 2012 , 48 , 240–246. [ CrossRef ] [ PubMed ] 15 Deng, Q.; Deng, L.; Miao, Y.; Guo, X.; Li, Y. Particle deposition in the human lung: Health implications of particulate matter from different sources Environ. Res 2019 , 169 , 237–245. [ CrossRef ] 16 Ou, C.; Hang, J.; Deng, Q. Particle deposition in human lung airways: Effects of airflow, particle size, and mechanisms Aerosol Air Qual. Res 2020 , 20 , 2846–2858. [ CrossRef ] 17 Chen, X.; Feng, Y.; Zhong, W.; Sun, B.; Tao, F. Numerical investigation of particle deposition in a triple bifurcation airway due to gravitational sedimentation and inertial impaction Powder Technol 2018 , 323 , 284–293. [ CrossRef ] 18 Deng, Q.; Ou, C.; Chen, J.; Xiang, Y. Particle deposition in tracheobronchial airways of an infant, child and adult Sci. Total Environ 2018 , 612 , 339–346. [ CrossRef ] [ PubMed ] 19 Feng, Y.; Kleinstreuer, C. Micron-particle transport, interactions and deposition in triple lung-airway bifurcations using a novel modeling approach J. Aerosol Sci 2014 , 71 , 1–15. [ CrossRef ] 20 Comer, J.; Kleinstreuer, C.; Hyun, S.; Kim, C. Aerosol transport and deposition in sequentially bifurcating airways J. Biomech Eng 2000 , 122 , 152–158. [ CrossRef ] 21 Yang, X.; Liu, Y.; So, R.; Yang, J. The effect of inlet velocity profile on the bifurcation COPD airway flow Comput. Biol. Med 2006 , 36 , 181–194. [ CrossRef ] 22 Kleinstreuer, C.; Zhang, Z. An Adjustable Triple-bifurcation unit model for air-particle flow simulations in human tracheobronchial airways J. Biomech. Eng 2008 , 131 , 021007. [ CrossRef ] 23 Kolanjiyil, A.; Kleinstreuer, C. Nanoparticle mass transfer from lung airways to systemic regions—Part I: Whole-lung aerosol dynamics J. Biomech. Eng 2013 , 135 , 121003. [ CrossRef ]

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 14 of 15 24 Nowak, N.; Kakade, P.; Annapragada, A. Computational fluid dynamics simulation of airflow and aerosol deposition in human lungs Ann. Biomed. Eng 2003 , 31 , 374–390. [ CrossRef ] 25 Zhang, Z.; Kleinstreuer, C.; Kim, C. Airflow and nanoparticle deposition in a 16-generation tracheobronchial airway model Ann Biomed. Eng 2008 , 36 , 2095–2110. [ CrossRef ] [ PubMed ] 26 van Ertbruggen, C.; Hirsch, C.; Paiva, M. Anatomically based three-dimensional model of airways to simulate flow and particle transport using computational fluid dynamics J. Appl. Physiol 2005 , 98 , 970–980. [ CrossRef ] [ PubMed ] 27 Gemci, T.; Ponyavin, V.; Chen, Y.; Chen, H.; Collins, R. Computational model of airflow in upper 17 generations of human respiratory tract J. Biomech 2008 , 41 , 2047–2054. [ CrossRef ] [ PubMed ] 28 Walters, D.; Luke, W. A Method for Three-dimensional Navier–Stokes simulations of large-scale regions of the human lung Airway J. Fluids Eng 2010 , 132 , 051101. [ CrossRef ] 29 Walters, D.; Luke, W. Computational fluid dynamics simulations of particle deposition in large-scale, multigenerational lung models J. Biomech. Eng 2010 , 133 , 011003–011008. [ CrossRef ] 30 Tian, G.; Longest, P.; Su, G.; Walenga, R.; Hindle, M. Development of a stochastic individual path (SIP) model for predicting the tracheobronchial deposition of pharmaceutical aerosols: Effects of transient inhalation and sampling the airways J. Aerosol Sci 2011 , 42 , 781–799. [ CrossRef ] 31 Ou, C.; Li, Y.; Wei, J.; Yen, H.; Deng, Q. Numerical modeling of particle deposition in ferret airways: A comparison with humans Aerosol Sci. Technol 2017 , 51 , 477–487. [ CrossRef ] 32 Weibel, E.; Cournand, A.; Richards, D Morphometry of the Human Lung ; Springer: Berlin/Heidelberg, Germany, 1963 33 Haider, A.; Levenspiel, O. Drag coefficient and terminal velocity of spherical and nonspherical particles Powder Technol 1989 , 58 , 63–70. [ CrossRef ] 34 Yeh, H.; Schum, G. Models of human lung airways and their application to inhaled particle deposition Bull. Math. Biol 1980 , 42 , 461–480. [ CrossRef ] 35 Longest, P.; Vinchurkar, S. Effects of mesh style and grid convergence on particle deposition in bifurcating airway models with comparisons to experimental data Med. Eng. Phys 2007 , 29 , 350–366. [ CrossRef ] 36 Hyatt, R.; Wilcox, R. The pressure-flow relationships of the intrathoracic airway in man J. Clin. Investig 1963 , 42 , 29–39 [ CrossRef ] 37 Tena, A.; Joaqu í n, F.; Eduardo, Á .; Pere, C.; Keith, W. Design of a numerical model of lung by means of a special boundary condition in the truncated branches Int. J. Numer. Method Biomed. Eng 2017 , 33 , e 2830. [ CrossRef ] [ PubMed ] 38 Farkas, Á .; Bal á sh á zy, I. Quantification of particle deposition in asymmetrical tracheobronchial model geometry Comput. Biol Med 2008 , 38 , 508–518. [ CrossRef ] 39 Inthavong, K.; Choi, L.; Tu, J.; Ding, S.; Thien, F. Micron particle deposition in a tracheobronchial airway model under different breathing conditions Med. Eng. Phys 2010 , 32 , 1198–1212. [ CrossRef ] [ PubMed ] 40 Li, Z.; Kleinstreuer, C.; Zhang, Z. Simulation of airflow fields and microparticle deposition in realistic human lung airway models. Part I: Airflow patterns Eur. J. Mech. B Fluids 2007 , 26 , 632–649. [ CrossRef ] 41 Comer, J.; Kleinstreuer, C.; Kim, C. Flow structures and particle deposition patterns in double-bifurcation airway models. Part 2 Aerosol transport and deposition J. Fluid Mech 2001 , 435 , 55–80. [ CrossRef ] 42 Zhang, Z.; Kleinstreuer, C. Transient airflow structures and particle transport in a sequentially branching lung airway model Phys Fluids 2002 , 14 , 862–880. [ CrossRef ] 43 Rahimi-Gorji, M.; Gorji, T.; Gorji-Bandpy, M. Details of regional particle deposition and airflow structures in a realistic model of human tracheobronchial airways: Two-phase flow simulation Comput. Biol. Med 2016 , 74 , 1–17. [ CrossRef ] [ PubMed ] 44 Li, Z.; Kleinstreuer, C.; Zhang, Z. Particle deposition in the human tracheobronchial airways due to transient inspiratory flow patterns J. Aerosol Sci 2007 , 38 , 625–644. [ CrossRef ] 45 Heged ˝us, C.; Bal á sh á zy, I.; Farkas, Á . Detailed mathematical description of the geometry of airway bifurcations Respir. Physiol Neurobiol 2004 , 141 , 99–114. [ CrossRef ] 46 Hofmann, W. Modelling inhaled particle deposition in the human lung—A review J. Aerosol Sci 2011 , 42 , 693–724. [ CrossRef ] 47 Darquenne, C. Aerosol deposition in health and disease J. Aerosol Med. Pulm. Drug Deliv 2012 , 25 , 140–147. [ CrossRef ] 48 Kleinstreuer, C.; Zhang, Z.; Kim, C. Combined inertial and gravitational deposition of microparticles in small model airways of a human respiratory system J. Aerosol Sci 2007 , 38 , 1047–1061. [ CrossRef ] 49 Lee, J.; Goo, J.; Chung, M. Characteristics of inertial deposition in a double bifurcation J. Aerosol Sci 1996 , 27 , 119–138. [ CrossRef ] 50 Longest, P.; Xi, J. Computational investigation of particle inertia effects on submicron aerosol deposition in the respiratory tract J Aerosol Sci 2007 , 38 , 111–130. [ CrossRef ] 51 Heyder, J. Deposition of inhaled particles in the human respiratory tract and consequences for regional targeting in respiratory drug delivery Proc. Am. Thorac Soc 2004 , 1 , 315–320. [ CrossRef ] [ PubMed ] 52 Soni, B.; Aliabadi, S. Large-scale CFD simulations of airflow and particle deposition in lung airway Comput. Fluids 2013 , 88 , 804–812. [ CrossRef ] 53 Fishler, R.; Ostrovski, Y.; Lu, C.; Sznitman, J. Streamline crossing: An essential mechanism for aerosol dispersion in the pulmonary acinus J. Biomech 2017 , 50 , 222–227. [ CrossRef ] [ PubMed ] 54 Harrington, L.; Prisk, G.; Darquenne, C. Importance of the bifurcation zone and branch orientation in simulated aerosol deposition in the alveolar zone of the human lung J. Aerosol Sci 2006 , 37 , 37–62. [ CrossRef ]

Int. J. Environ. Res. Public Health 2023 , 20 , 4583 15 of 15 55 Balashazy, I.; Alfoldy, B.; Molnar, A.; Hofmann, W.; Szoke, I.; Kis, E. Aerosol drug delivery optimization by computational methods for the characterization of total and regional deposition of therapeutic aerosols in the respiratory system Curr. Comput Aided Drug Des 2007 , 3 , 13–32. [ CrossRef ] 56 Martonen, T.; Katz, I. Deposition patterns of aerosolized drugs within human lungs: Effects of ventilatory parameters Pharm. Res 1993 , 10 , 871–878. [ CrossRef ] [ PubMed ] 57 Cheng, Y.; Zhou, Y.; Su, W. Deposition of particles in human mouth-throat replicas and a USP induction port J. Aerosol Med. Pulm Drug Deliv 2015 , 28 , 147. [ CrossRef ] 58 Sul, B.; Wallqvist, A.; Morris, M.; Reifman, J.; Rakesh, V. A computational study of the respiratory airflow characteristics in normal and obstructed human airways Comput. Biol. Med 2014 , 52 , 130–143. [ CrossRef ] 59 Xi, J.; Yuan, J.; Yang, M.; Si, X.; Zhou, Y.; Cheng, Y. Parametric study on mouth–throat geometrical factors on deposition of orally inhaled aerosols J. Aerosol Sci 2016 , 99 , 94–106. [ CrossRef ] 60 Lai, T.; Morsi, S.; Das, S.; Owida, A. Numerical analysis of particle deposition in asymmetrical human upper airways under different inhalation cycles J. Mech. Med. Biol 2013 , 13 , 75–277. [ CrossRef ] 61 Park, S.; Wexler, A. Size-dependent deposition of particles in the human lung at steady-state breathing J. Aerosol Sci 2008 , 39 , 266–276. [ CrossRef ] 62 Phuong, N.; Yamashita, M.; Yoo, S.; Ito, K. Prediction of convective heat transfer coefficient of human upper and lower airway surfaces in steady and unsteady breathing conditions Build. Environ 2016 , 100 , 172–185. [ CrossRef ] 63 Rahimi-Gorji, M.; Pourmehran, O.; Gorji-Bandpy, M.; Gorji, T. CFD simulation of airflow behavior and particle transport and deposition in different breathing conditions through the realistic model of human airways J. Mol. Liq 2015 , 209 , 121–133 [ CrossRef ] 64 Zhang, Z.; Kleinstreuer, C. Effect of particle inlet distributions on deposition in a triple bifurcation lung airway model J. Aerosol Med 2001 , 14 , 13–29. [ CrossRef ] [ PubMed ] Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Other Environmental Sciences Concepts:

[back to top]

Discover the significance of concepts within the article: ‘Particle Deposition in Large-Scale Human Tracheobronchial Airways Predicted...’. Further sources in the context of Environmental Sciences might help you critically compare this page with similair documents:

Editor, Deposition, Responsibility, Air pollution, Particle, Particle size, PubMed, Dose effect, Particulate matter, Reynolds number, Particle deposition, Gravitational Sedimentation, Airflow rate, Numerical method, Computational fluid dynamics, Adverse health effect, Inertial impaction, CrossRef, Airflow, Conducting airway, Individual author(s), Injury to people, Statements, opinions and data, Three dimensional, Lung lobes, MDPI, Statements opinions and data, Ideas methods instructions or products, Static pressure, Deposition pattern, Human airways, Boundary condition, Gravitational acceleration, Governing equations, Finite volume method, Deposition mechanism, Deposition Efficiency, Atmospheric Aerosols, Drag Coefficient, No-Slip Velocity Boundary Condition, Airflow Structure, Grid Convergence Index, Pressure Outlet, Inhalation Rate, Particle deposition efficiency, Particle transport, Brownian diffusion, Atmospheric dispersion model, CFD simulation, Breathing conditions, Realistic model.