Justus Kavita Mutuku1, Wen-Che Hou This email address is being protected from spambots. You need JavaScript enabled to view it.1, Wei-Hsin Chen  2,3,4

1 Department of Environmental Engineering, National Cheng Kung University, Tainan 70101, Taiwan
2 Department of Aeronautics and Astronautics, National Cheng Kung University, Tainan 70101, Taiwan
3 Department of Chemical and Materials Engineering, College of Engineering, Tunghai University, Taichung 407302, Taiwan
4 Department of Mechanical Engineering, National Chin-Yi University of Technology, Taichung 41170, Taiwan


Received: March 9, 2020
Revised: April 16, 2020
Accepted: April 17, 2020
Download Citation: ||https://doi.org/10.4209/aaqr.2020.03.0107 

  • Download: PDF


Cite this article:

Mutuku, J.K., Hou, W.C. and Chen, W.H. (2020). Two-phase Flow Dynamics and PM2.5 Deposition in Healthy and Obstructed Human Airways during Inhalation. Aerosol Air Qual. Res. 20: 1094-1110. https://doi.org/10.4209/aaqr.2020.03.0107


HIGHLIGHTS


  • This study numerically investigated two-phase flow in the human airways.
  • Inhalation curves for rest, light activity, and moderate exercise were used.
  • DFs for different PM2.5 sizes were 0.12–1.18% in healthy and 0.05–0.49% in COPD.
  • Important factors in COPD; jet flows, skewed mass flow rates and dean vortices.
  • Depositions at the carina regions could be attributed to inertial impaction.
 

ABSTRACT


Lately, there has been an unexpected increase in the global burden of diseases due to inhalation of ambient air contaminated with fine particulate matter (PM2.5). The adverse effects due to inhalation of PM2.5 are more exaggerated for patients with obstructive lung diseases such as Chronic obstructive pulmonary disease (COPD). In this study, numerical investigations were carried out on the effects of particle diameters and different Reynold numbers at the inlet on deposition patterns and fractions in COPD patients and healthy individuals. Three real and transient inhalation curves with average Reynolds numbers of 164.3, 362.4, and 606.4 were used to represent rest, light activity, and moderate exercise inhalation statuses. Four median diameters including 0.075 µm, 0.15 µm, 0.3 µm, and 0.6 µm were injected at the inlet of the control volume (G5–G8) of the Weibel’s lung geometry to represent the particle size distribution for a PM2.5 concentration of 50 µg m–3. Deposition fractions and patterns were obtained from tracking a total of 350031, 692596, and 833553 PM2.5 particles corresponding to rest, light activity and, moderate exercise respectively. Deposition fractions (DFs) for the different sizes of PM2.5 ranged between 0.12% and 1.18% in the healthy airway geometry and between 0.05% and 0.49% in the COPD case. The deposition patterns were skewed for a COPD case due to jet flow phenomena, skewed mass flow rates, and the induced dean vortices. While depositions at the carina regions could be attributed to inertial impaction, those along the bifurcations could be attributed to centrifugal forces, complex secondary flows, and inertial impaction. The effects of the deposition mechanisms varied between the two geometries, among different particle diameters and inhalation statuses.


Keywords: Aerosol; Chronic obstructive pulmonary disease (COPD); Obstructed; Numerical methods; Two-phase flow; Deposition fractions (DFs).


NOMENCLATURE


A             Amplitude of the real breathing curve (Reynolds number)
D             Diameter (mm)
F             Force (N)
G             Generation of Weibel’s airway
C             Carina region of the Weibel’s airway
P             Pressure (Pa)
Q             Flow rate (L s–1)
R             Radius (mm)
Ui            Airflow velocity (m s–1)
Fpf           Viscous forces on the fluid molecules.
X           Characterizes the Cartesian coordinates,
σij            Shear stress
r              Radial coordinate
Re           Reynolds number
t              Time (s)
Umean       Average velocity at the airway’s entrance (m s–1)
Umax        Max velocity at the airway’s entrance (m s–1)
V             Control volume

Greek Letters

ρ             Density (kg m–3)
σ             Stress (N m–2)
µ             Viscosity Coefficient
υn            Averaged inlet velocity of parent branch (m s–1)
θ              Angular coordinate (°)


Subscripts

mean     average value
peak       peak value
I              y-coordinates
j              z coordinates


INTRODUCTION


Scientific evidence shows that obstructive respiratory diseases are mostly caused by exposure of the smooth muscles to irritants such as particulate matter in ambient air which causes contraction of the smooth muscles and subsequent inflammation (Sul et al., 2014). Studies have linked fine particulates (PM2.5) to a wide range of acute and chronic health defects as well as their exacerbation (Brook et al., 2010; Li et al., 2018; Jiang et al., 2019). Potential pathways for the adverse health effects include; inflammation, oxidative stress, and altered immunity (Ni et al., 2015; Chao et al., 2018; Zhang et al., 2018). A case in point is the mounting evidence suggesting that air pollution from the combustion of fossil fuels is responsible for several adverse health effects in urban areas (Valavanidis et al., 2008; Lelieveld et al., 2015). This information is only available for developed and middle-income countries (Cohen et al., 2005). Model estimations depict a worse situation in developing countries although the true situation remains conclusively unknown (Soriano et al., 2017).

The type of resultant adverse impacts due to inhalation of PM2.5 is dependent on their chemical composition and mass concentrations (Valavanidis et al., 2008; Brook et al., 2010). The toxicity of particulate matter is inversely proportional to the particle size. As an implication, fine and ultra-fine particulate matters tend to be most toxic among the total suspended solids (Valavanidis et al., 2008). Some studies have shown a steady and autonomous relationship between the concentration of fine particulate matter in the ambient air and the serious health defects caused by air pollution (Li et al., 2018). Further, a statistically significant association has been established between long term exposure to PM2.5 and reduced life expectancy (Pope III et al., 2009). Specifically, reports on ambient air pollution and mortality show that PM2.5 is responsible for more than 9% of the global annual mortality rates due to cardiopulmonary disease, cancer of the respiratory system, and acute conditions of the respiratory system in children under the age of 5 years (Cohen et al., 2005). In 2012, the residential areas of more than 80% of the global population was believed to exceed annual WHO AQG for fine particulates of 10 µg m3 (Rao et al., 2012).

Environmental protection agencies apply indexes for the enforcement of air quality standards and in efforts to protect the public health from harmful pollutants suspended in the air (Valavanidis et al., 2008). The creation of these indices requires an extensive understanding of the chain of processes involved with air toxics from emission, intake by humans and their physiological effects. Earlier studies have shown a fascinating, but not very conclusive understanding of the possible pathways that associate exposure to aerosols in the ambient air and especially PM2.5 and mortality caused by cardiovascular disease (Brook et al., 2010). As such, a lot of recent studies are focused on aerosols transportation, transformation and deposition inside the human lungs upon inhalation with ambient air.

The bulk of atmospheric aerosols have a smooth and spherical shape and their density range between 1.1–2 g cm–3 (Liu et al., 2015). In other studies, PM2.5 densities in Pittsburg and Beijing were found to be 1.5 ± 0.3 g cm–3 and 1.41 ± 0.40 g cm–3 respectively (He et al., 2001; Rees et al., 2004). The variation of velocity with time at the mouth of the nasal cavity can be satisfactorily be represented by a sinewave. Real inhalation curves vary slightly from sinewaves of equivalent amplitude. Previous numerical investigations in the past have used both uniform velocity and parabolic velocity distribution at the inlet, the latter is a better representative of the real situation in the human airways (Lambert et al., 2011; Augusto et al., 2016; Kolanjiyil and Kleinstreuer, 2017).

COPD is characterized by persistent resistance of air-flow in the human airways. This condition is often progressive and the response to inhalation of toxins by COPD patients increases with time (Viegas et al., 1996). Histological studies have characterized this condition with a constriction along the affected airway tubes. The constriction is irregular in human-specific lung models while in idealized geometric models, the constriction is usually axisymmetric (Luo et al., 2007; Chen et al., 2012; Mutuku and Chen, 2018).

The resistance of airflow which is important in the characterization of obstructions in the human lung is represented by the ratio of total pressure drop (Sul et al., 2014). Previous research has shown a direct proportionality between pressure drop and air-flow rate due to amplified secondary flows. The resistance of flow is thus a very important property, especially when dealing with obstructed human airways (Yanai et al., 1992; Kang et al., 2011). There is a significant scale of research on air dynamics, aerosol motion, particle deposition and particle clearance in healthy and obstructed human lungs using numerical methods. Numerical analysis has been used extensively in the past to characterize experimentally inaccessible measurements (Chen, 2001; Chen et al., 2014; Chen et al., 2018b; Zhang et al., 2019). Most of these studies involve multi-phase processes (Chen, 2002; Morel, 2015). In a study by Chen et al. (2012), the deposition fractions (DFs) of pharmaceutical aerosols in obstructed human airways for unsteady inlet conditions exceeded those of a steady one. Studies conducted on obstructed human airways before have concentrated on air flow dynamics, pressure distribution, etiology, and pathogenesis in exacerbation of COPD (Jain et al., 2013; Sul et al., 2014; Ni et al., 2015). However, some important information is still missing in the literature concerning the deposition of PM2.5 with different particle size distribution in the human airways affected by COPD. This paper is aimed at filling the missing gap by applying numerical analysis tools in the assessment of airflow, aerosol motion, and deposition in healthy and obstructed human airways. The effects of particle size distribution, Reynolds number at the inlet, and geometry on deposition patterns and fractions inside the human lungs are discussed. Exceptional features and essential aspects of results from the numerical analyses are also discussed. Evaluating the deposition patterns and fractions results would give supplementary information needed to understand the exacerbation of COPD as a result of inhalation of PM2.5.


NUMERICAL FORMULATIONS



Solid Geometry and Mesh of the Bifurcation

The real human lung is very complicated and simplification of the geometry for numerical analysis is the first step. There are two well-known and popular lung geometries for the human lungs; Weibel model which is symmetrical and Horsfield model which is asymmetrical (Weibel, 1965; Horsfield et al., 1971). In a study comparing the particle deposition fractions for the Weibel and Horsfield models found higher deposition fractions in Horsfield's model with respect to the whole lung’s model (Kim and Iglesias, 1989). The Weibel model has higher accuracy and reliability in terms of physiological representation of the human lungs and thus it is more popular among researchers (Weibel, 1963). A comprehensive depiction of the Weibel human lung geometry can be found in other literature (Tsuda et al., 2008; Guha et al., 2016).

The oral cavity and upper regions of the respiratory system are rarely affected by chronic inflammation due to the high efficiency of the lung clearance mechanisms in the upper sections. Therefore, COPD in G5–G8 is one of the few possibilities for a patient affected by airway obstruction and for this reason was adopted for the study. The solid geometries for the healthy and obstructed three-dimensional geometry and mesh grids for G5-G8 of the Weibel’s model were created using SOLIDWORKS and are as depicted in Fig. 1. The solid geometry starts at G5 as the first parent generation which branches into two at a bifurcation angle of 35° with regard to the plane z = 0. The demarcation of the computational domain into grids was based on unstructured tetrahedral elements on the ANSYS FLUENT 19.0 platform. Important parameters in the creation of a physiologically accurate lung model including diameter, length, branching angle, and radius of the curvature are provided in previous works (Chen et al., 2012; Mutuku and Chen, 2018). The carina regions were smoothened after lofting mother and daughter branches using the fillet options. In order to differentiate deposition fractions in the two mechanisms under study here; inertial impaction and complex secondary flow phenomenon, the solid geometries were separated into carina regions and curved tube sections as shown in Fig. 2. The geometry of airways deformed due to COPD was characterized by an axisymmetric constriction at the generation affected by the obstruction. For this study, the narrowest part of the constriction had a 50% reduction in diameter as compared to the unobstructed one, as shown in Fig. 1(b). The extend of obstruction was chosen in a way that the volume reduction and airway surface reduction are consistent with histological studies conducted on obstructed airways (Zhang and Papadakis, 2010; Chen et al., 2012).

Fig. 1. Schematics of (a) block structure for healthy geometry (G5–G8) (b) block structure geometry for G5–G8 wherein G6-1 is affected by COPD (c) mesh for the healthy geometry, and (d) mesh for the airways affected by COPD on G6-2. 

Fig. 2. Labeled branches (G5–G8) and carina regions for (a) healthy geometry, and (b) geometry affected by COPD on G6-1.Fig. 2. Labeled branches (G5–G8) and carina regions for (a) healthy geometry, and (b) geometry affected by COPD on G6-1.


Governing Equations

Herein, investigations concerning airflow fields as well as PM2.5’s deposition patterns and fractions for inhalation during three different statuses were investigated. The air flow was assumed to be incompressible, isothermal, and laminar. The density (ρ) and kinematic viscosity (µ) of the air were assumed to remain constant throughout the inhalation cycle. In addition, the force of gravity on the fine particulates was assumed to be insignificant (Johnson and Jackson, 1987; Zhang et al., 2001). Conservation equations for continuity equation and momentum were applied to govern the gas flow and discrete particle motion.

where, xi is the coordinate in the cartesian system, t stands for the time elapsed, and Ui and P represent the local velocity of air flow and pressure, respectively. Air density is represented by ρ while σij represents stress. The y and z coordinates are represented by the subscripts i and j . Lastly, Fpf signifies the viscous forces on the fine particulate matter.

The motion of fine particulates in the control volume was governed by Newton's second law of motion which is expressed as:

Assuming the effects of external forces, temperature changes and magnetic field were insignificant and hence ignored, the trajectories of the fine particulate matter can be expressed as:

where an individual PM2.5’s mass was denoted as mp, while its diameter and velocity were denoted as d and up, respectively. The inertial force exerted by the air-flow to the fine particulate matter was assumed to be the only external force. Further, the fine particulates were assumed to be indestructible and non-distortable under centrifugal forces at the bifurcation as well as the complex secondary flow phenomena in the radial direction of the tube sections.

The equation for drag coefficient (Cd) for fine particulates could be expressed as,

where α1, α2, and α3 are constants obtained through experiments. The PM2.5’s Reynolds number (ReN) can be expressed as.

 

The concept of a fluid volume fraction αf was used to express the ratio of the volume of the aerosol to that of the air in which it was suspended, and it was expressed as

where the particles, their orders in a sequence, and their volume were denoted by p, k, and Vk,p, respectively. The total volume of the computational domain was denoted as V. For a pure gaseous phase flow αf = 1.


Boundary Conditions


Velocity Parameters at the Inlet

The flow curve of average air-flow velocity versus time at the inlet was described using a real inhalation curve (Mutuku and Chen, 2018). Real inhalation curves were extrapolated from a study whereby a spirometry test on the velocity variation at the mouth with time had been carried out (Tena et al., 2015). In most of the studies involving aerosol and air motion in the human lungs, the breathing statuses are delineated based on the level of human activity and consequently Reynolds number of flows. The most commonly studied breathing statuses include rest condition (sedentary), light activity, and moderate exercise. For ordinary breathing, the expansion and contraction of the airways are insignificant as compared to the axial airflow velocity, therefore the walls of the airways can be assumed to be rigid (Hughes et al., 1972). The three real inhalation transient curves applied for these investigations are depicted in Fig. 3 which shows the instantaneous and average Reynolds number (Remean) variations with inhalation time at the inlet (Zhang and Kleinstreuer, 2002). An equation described by Weibel (1965) for obtaining the average inlet velocity in this case G5 (Vn) was applied in this study

where the generation’s number, diameter, and the airflow volumetric at the inlet are denoted by n, D, and Qin, respectively.

The velocity profile at the inlet’s circular cross-section was described as a completely formed parabolic flow using a user-defined function (UDF) (Mutuku and Chen, 2018). The airways were smooth and rigid, and a non-slip boundary condition was implemented on the tube’s walls.

Fig. 3. Real inhalation waveforms for (a) rest (b) light activity and (c) moderate exercise.Fig. 3. Real inhalation waveforms for (a) rest (b) light activity and (c) moderate exercise.



Injection of Particles and their Motion in the Control Volume

The PM2.5 particle size distribution applied for this study was obtained from earlier studies (Chen et al., 1997; Lin et al., 2008) and it is shown in Figs. 4(a)–4(c). The fine particulates were assumed to have an average particle density of 1.5 g cm–3 (Rees et al., 2004). The particles had a deterministic-parabolic distribution at the inlet which was earlier investigated by Zhang and Kleinstreuer (2001) and found to give depositional fractions similar to that of a random distribution. Further information concerning the number of fine particulates inhaled during a single breathing cycle is summarized in Table 1. The ratio of the local concentration of PM2.5 (C(r)) to the average concentration (C) can be represented by the following equation.

where n is the number of particles in a given ring, for instance ra  r  rb. The radius at the inlet and mathematical function for truncating the number of the particles are denoted by R and Int[ ], respectively.

For consistency purposes, the same particle size distribution for a PM2.5 concentration of 50 µg m–3 was applied in this study, similar to the one applied by Chen et al. (2018a). Therefore, a constant number of fine particulates were injected at the inlet of the control volume for every time step that is; 6474, 12810 and 25846 for rest, light activity, and moderate exercise, respectively. During the entire numerical calculation, the momentum of the particles was caused by the forces of the surrounding mass of air. The boundary conditions for the discrete phase at the inlet and outlet of the control volume were set as “escape into the control volume” and “escape out of the control volume”, respectively (Chen et al., 2018a). Since the lung’s walls are moist and covered with mucus, the trapping fraction for all the particles that touched the walls was 100%. The gauge pressure was set as zero at all outlets in both the healthy and diseased human airway geometries (Inthavong et al., 2010). PM2.5 disintegrates only under extreme environmental conditions such as high temperatures and pressures which don’t usually happen in the natural ambient environment or inside the human airways (Zhang et al., 2016). Therefore, in the mathematical formulation, one of the simplifying assumptions is that the particles do not break or undergo any form of deformation during the full inhalation cycle. 


Numerical Analysis

Simulation of two-phase flow in the human airways applies one-way coupling whereby the gas phase affects the solid phase only, but there are no feedback forces from the dispersed phase (Comer et al., 2000; Luo et al., 2007; Rahimi-Gorji et al., 2015). The discrete phase model of CFD is activated by defining several parameters of the dispersed phase, among them their position in the control volume under investigation, velocity, diameter, temperature, mass flow, and time of injection (Rahimi-Gorji et al., 2015). The trajectory calculations use the initial location and parameters of the fluid flow for the calculations.

Conservation equations for both the gas flow and discrete phase motion were solved using ANSYS FLUENT 19.0. Coupled algorithm was adopted because the flow in the region under investigation was laminar and the convergence was determined by pressure-velocity coupling. Previous investigations for instance by Rahimi-Gorji et al. (2015) used the same algorithm. During these investigations, the Navier-Stokes equations and coupled algorithm were solved simultaneously. A three-order MUSCL scheme and Green-Gauss Node were applied as the method of discretization of the fluid’s momentum and base for gradient selection respectively. A standard pressure interpolation scheme and a relaxation coefficient of 0.75 for both velocity and pressure were applied. Before running the numerical simulations, a full breathing cycle was run with the only gaseous flow to get rid of startup effects (Inthavong et al., 2010).

The solutions for the solid phase motion was obtained using a hybrid of Euler and Lagrange approaches provided as the discrete phase model (DPM) in ANSYS FLUENT 19.0 platform. The effects of the inertial force of the continuous phase on the trajectories for the solid phase and their reaction forces were solved for each element in the computational domain. The solution for this couple of conservation equations was attained through solving the equations for the discrete and continuous phase equations alternatively until they both converged. To ensure a smooth breathing curve for all inhalation statuses, a time step of Δt = 0.035 s was selected for solving the air-flow fields and paths of the suspended particles. For a single time-step, the solutions of the flow fields converged when the continuity residual had a relative error of less than 10–5.


Dimensionless Numbers

The inlet’s Reynolds number (Red) is related to airflow velocity at the geometry’s inlet as shown below (Verbanck et al., 2011)

where ρ is the fluid’s density (=1.225 kg m–3), μ represents the viscosity coefficient of air (1.789 ´ 10–5 kg m–1 s–1), and D is the diameter at the inlet of the parent tube (G5) which is 3.5 ´ 10–3 m. The flow velocity at the entrance of G0 (i.e., the mouth’s entrance) is uniform. But as the flow proceeds deeper into the lower generations, the velocity profile becomes parabolic distribution due to the drag force on walls of the airways. The velocity distribution at the inlet is parabolic, and the maximum flow velocity (UMAX) at the middle of the cross-section is twice the Umean. Since the fluid’s density and the viscous coefficient are fixed, the larger the Umean and UMAX, the greater the Reynolds number, hence the greater the inertial force of the suspended particles.

The formula for deposition fraction (DF) represents the number of particles deposited as a percentage of the total number of particles entering a particular section of the control volume (Kim and Iglesias, 1989), and is defined as:

For a fixed number of particles in the tube, the larger the DF, the higher the number of particles deposited.


Model Validation and Grid Independence

The axial flow velocity diagrams, velocity contours, and secondary flow vectors for selected cross-sections as well as deposition patterns during inhalation were validated whereby the air flow rate for the inlet specified case was 250 mL s–1. Schematics for the validation of axial air velocity contours as well as secondary vectors are shown in Fig. 5(a). From the cross-sectional velocity and secondary vector distribution diagrams, it was observed that there were four vortices generated in all the three cross-sections, which was also consistent with the phenomenon observed in an earlier study (Nazridoust and Asgharian, 2008). The overall trends for deposition from this study are consistent with the literature, as shown in Fig. 5(b), and the particles’ deposition patterns were significantly influenced by the different boundary conditions imposed on the breathing cycles. The deposition fractions for 0.1, 1 and 10 µm particles were 7.4, 0.11, and 0.029%, respectively. This closely agreed with earlier findings (Nazridoust and Asgharian, 2008), thereby confirming the accuracy and confidence of the developed model for the prediction of PM2.5 motion and deposition patterns as well as fractions. Through the above verification, it was clear that the results from the coupled algorithm can be applied to predict air and particle motion in the human airways (Chen et al., 2018a). Therefore, the coupled algorithm was used to simulate the instantaneous two-phase flow in the multiple bifurcation respiratory tracts.

Three different grids with 0.44 million, 1.37 million and 1.87 million elements were applied to check the grid’s independence for the deposition fraction of the particles in a COPD geometry during inhalation for light activity as shown in Fig. 6. With 12806 particles inhaled during the light activity which was characterized by 57 time-steps of 0.035 s each, it was discovered that 1.37 million elements in the mesh were enough to capture the particle deposition patterns and fractions within a reasonable time frame and without significantly compromising on accuracy and reliability of the numerical analysis. Similarly, 1.37 million elements were adopted for simulating the healthy case. 

Fig. 5. Validation for (a) axial velocity distribution and (b)Validation for particle deposition with three different diameter sizes 10 µm, 1 µm and 0.1 µm using results from an investigation by Nazridoust and Asgharian (2008). Fig. 5. Validation for (a) axial velocity distribution and (b)Validation for particle deposition with three different diameter sizes 10 µm, 1 µm and 0.1 µm using results from an investigation by Nazridoust and Asgharian (2008). 

Fig. 6. Grid independence test for the COPD mesh.Fig. 6. Grid independence test for the COPD mesh.


RESULTS AND DISCUSSION



Fluid Flow Dynamics


Axial Airflow Fields

Airflow fields are very key in determining the deposition patterns of particulate matter in human airways. Several parameters can be used to explain fluid flow in the human airways such as velocity distribution, stagnation zones, recirculation zones, jet flows, and shear stress on the walls (Zhang et al., 2002; Kleinstreuer and Zhang, 2010; Chen et al., 2012). The airflow velocity distribution affects pressure distribution which was an important parameter in determining hotspots for aerosol deposition (Luo et al., 2007; Kang et al., 2011). The velocity distribution inside the control volume was responsible for the formation of stagnation and recirculation zones as well as jet flow phenomena in the constriction as shown in Fig. 7 (Chen et al., 2012; Chen et al., 2018a). The jets were extracted at the cross-section in the middle of each generation under investigation. These flow phenomena, as a result of velocity distribution, greatly affected particle deposition patterns in both healthy and deformed human airways. Jets tend to form in obstructed sections according to the results from earlier analysis (Kleinstreuer and Zhang, 2010; Sul et al., 2014). Boundary layer separation in the regions near the constrictions leads to the recirculation phenomenon as it can be seen in the vicinity of the obstruction in Fig 7(b) (Sul et al., 2014). Further complex secondary flow phenomena can lead to the formation of Dean vortices. Prior to investigating the two-phase flow, detailed investigations and discussions for the pure gas-phase flow for the geometries discussed here were performed (Mutuku and Chen, 2018). 

Fig. 7. Stream lines of flow and jet flow phenomenon for (a) healthy human airways and (b) COPD geometry during (i) rest, (ii) light activity and, (iii) moderate exercise.Fig. 7. Stream lines of flow and jet flow phenomenon for (a) healthy human airways and (b) COPD geometry during (i) rest, (ii) light activity and, (iii) moderate exercise.


Wall Shear Stress

Wall shear stress is an important factor in the immune response mechanism of the lungs through the release of ATP, which increases the intracellular concentration of calcium for cell engulfment and anti-inflammatory response (Gronski et al., 2009). Fig. 8 depicts that the wall shear stress was maximum in the vicinity of the carina regions as well as around the constriction in the geometry affected by COPD. Jet flow phenomenon in the geometry affected by COPD induced high-velocity gradients in the radial direction. Since wall shear stress is proportional to the velocity gradients perpendicular to the axial airflow direction, the location of the maximum shear stress adhered to the regions where high axial velocities were greatly skewed. Low shear stresses (below 0.3 Pa) have a positive impact on the barrier function of the epithelial lining of the airway walls (Sul et al., 2014). On the other hand, highly exaggerated shear stresses can damage the epithelial cells.

During rest, the healthy generation had shear stress exceeding 0.3 Pa at carina 5 (C5) only, as compared to the COPD geometry whereby exaggerated shear stresses were present at C5, the obstructed generation (G6-1), as well as the carina 6-2 (C6-2). This was because high resistance of flow at the obstructed section led to a higher flow rate in the unobstructed generation. Therefore, high shear stress was generated once the amplified velocity of flow impacted on the carina region. During the light activity, all carina regions except C7-3 and C7-4 had high shear stresses. Additionally, the wall area in the vicinity of the constriction with the highest shear stresses is larger than that observed in the same geometry during rest. The regions of airway walls affected by high shear stress during moderate exercise corresponded to those of light activity but with a larger effective wall area for the hotspots, as shown in Fig. 8. Overall, the geometry affected by COPD had more regions of high shear stress compared to the healthy human airway, indicating a high likelihood of epithelial damage in the obstructed human airways. 

Fig. 8. Wall shear stresses for (a) healthy human airways and (b) COPD geometry during (i) rest, (ii) light activity and, (iii) moderate exercise.Fig. 8. Wall shear stresses for (a) healthy human airways and (b) COPD geometry during (i) rest, (ii) light activity and, (iii) moderate exercise.


Deposition Patterns for PM2.5

Three real and transient inhalation curves with an average Reynolds number of 164.3, 362.4, and 606.4 were applied for these investigations as shown in Fig. 3. A PM2.5 concentration of 50 µg m–3 was applied where particles with four different diameters including 0.075 µm, 0.15 µm, 0.3 µm, and 0.6 µm as presented in Fig. 4 were injected at the inlet of G5. The PM2.5 particles exceeding 1 µm were excluded from these investigations because they are easily deposited in the oral, nasal and upper sections of the human airways. In this simulation, G5 is the inlet of the control volume and only particles with less than 1 µm in diameter reach this section of the respiratory system. The results presented herein depict the deposition fractions and patterns obtained from tracking a total of 350031, 692596 and 833553 PM2.5 particles during rest, light activity, and moderate exercise inhalation statuses, respectively. 

Fig. 4. Distribution of PM2.5 injected at the inlet for each time step during (a) rest and (b) light activity and (c) moderate exercise. Fig. 4. Distribution of PM2.5 injected at the inlet for each time step during (a) rest and (b) light activity and (c) moderate exercise.

There are seven main deposition mechanisms for aerosols in the human airways including (1) turbulent mixing; (2) inertial impaction; (3) impaction due to complex secondary flows; (4) centrifugal forces; (5) gravitational sedimentation, (6) Brownian motion; and (7) electrostatic precipitation (Finlay and Martin, 2008; Darquenne, 2012). The dominance of any of these mechanisms of deposition varies as the particles advance from the oral cavity, through the upper section of the human lungs, central, and the lower sections, and later into the alveolar region. The effect of gravitational sedimentation in the overall deposition is usually more pronounced for aerosols exceeding 10 µm (Hofmann et al., 1995). The present study focused on particles with diameters smaller than 2.5 µm. Therefore, the effect of gravity on particles was ignored. Additionally, diffusional forces are more dominant in the alveolar section of the human lungs only and therefore was ignored for this study. In these investigations, the impaction of the PM2.5 on the airway walls due to inertial forces, centrifugal forces, and complex secondary flow phenomena were investigated. The most important obstacles for inertial impaction in this study are the carina regions. Therefore, to study the hotspots for deposition, both the healthy and obstructed geometries were segmented into; carina region and other sections of the control volumes herein referred to as tube walls as depicted in Fig. 2.

Several factors influencing the deposition fractions and patterns of inhaled particulate matter include (1) geometry of the tracheobronchial model; (2) inhalation status; and (3) chemical and physical properties of the particles (Daigle et al., 2003; Zhang and Papadakis, 2010; Verbanck et al., 2011; Li et al., 2018). The chemical properties of PM2.5 were not considered for this study; therefore, the particles were assumed to be inert. Since they had a constant density, the PM2.5 diameters can be directly linked to their Stokes’ numbers. Similar to previous investigations, the results showed that deposition fractions for aerosols are dependent on their Stokes numbers since inertial impaction is the main deposition mechanism in the human airways (Cheng et al., 1999; Zhang et al., 2002). As it can be seen in Table 2, deposition fractions in the carina regions dominated DFs in the whole geometry. The fine particulates under investigation here were assumed to have a spherical shape and consequently, their Stokes number increased with an increase in diameter. 


Deposition on the Carina Regions

Aerosols deposit on carina regions due to abrupt geometrical changes. Therefore, inertial impaction was the main mechanism credited for deposition (Huang and Zhang, 2011). Herein, the stokes number was the ratio of PM2.5 stopping distance to characteristic distance from the carina regions (Wessel and Righi, 1988). As such fine particulates with higher stokes number tend to deviate from the streamlines of airflow at the carina region, while low stokes number particulates adhere to the streamlines of flow.

As can be seen in the graphic representation of deposition patterns in Figs. 9 and 10, the most important hot spot for PM2.5 deposition from this study was carina 5 (C5). It is the first carina region in this control volume, and it received the highest DF of fine particulates, as seen in Table. 2. The property of deviation of PM2.5 from the streamlines of flow in the two-phase flow under investigation is best described using the stokes number which is expressed as;

where St is the Stokes number, dp and ρp are the diameter and density of the PM2.5, respectively; u and μ are properties of the gas phase that are the average velocity and dynamic viscosity, respectively, and  is the diameter of the airway. Overall, deposition by inertial impaction increases with an increase in the Stokes number at the inlet of the geometry. For all the carina regions (C5 to C7-4) deposition fractions for the healthy geometry exceeded that of COPD geometry for 12 out of 21 instances. Overall, deposition fractions were highest on carina 5 and lowest for carina 7-1 and 7-4 for both geometries. Important factors for the deposition of the PM2.5 on the carina regions include (1) the position in the control volume, (2) airflow velocity, and (3) particle distribution at the inlet which was parabolic deterministic in this case (Zhang and Kleinstreuer, 2001). Distal carina regions had lower deposition fractions as compared to the proximal ones.

Since the effective cross-sectional area increases as air advances from the inlet (G5) to the outlet (G8), the velocity decreases. This consequently led to a drop in the inertial impaction forces as the particles progress from the inlet to the outlet. For carina 5 and 7-2, particle deposition in the healthy human airway exceeded that in the COPD geometry for all inhalation statuses. In carina 6-1 ad 6-2, particle deposition in the COPD geometry exceeded that in the healthy geometry during all inhalation statuses. In carina 6-2, the deposition fraction of PM2.5 for the COPD geometry exceeded that deposited the healthy geometry by 22.33, 231.36, and 248.52% for rest, light activity, and moderate exercise respectively. In carina 6-1 which succeeded the obstructed generation, total deposition fractions in the COPD geometry exceeded the ones in the healthy generation for all inhalation statuses by 137.86, 77.05, and 47.45% for rest, light activity, and moderate exercise, respectively. This can be attributed to jet flow phenomena at the obstructed generation. High axial velocities associated with the jet flow led to amplified DFs at the succeeding carina region (C6-1). Due to the high resistance of flow in the obstructed generation (G6-1), the axial velocity of flow in G6-2 was higher compared to the corresponding generation in a healthy geometry. Consequently, the deposition in C6-2 and C6-1 in COPD human geometry exceeded those in a healthy geometry for all inhalation statuses.

Overall PM2.5 deposition fractions at all the carina regions of the healthy human airway geometry exceeded that of COPD by 9.28, 18.20, and 13.05% for rest, light activity, and moderate exercise respectively. Fig. 11 depicts a scatter of the deposition fractions against the particles’ Stokes numbers at the carina regions, tube walls, and overall deposition at the end of the full inhalation cycle for a COPD and healthy geometry. From the scatter in Fig. 11(a), it is clear that deposition in the carina region was more strongly related to the Stokes number in a healthy geometry than in a COPD geometry. Overall, the deposition of PM2.5 in the carina regions was more important as it was responsible for 64–88% of the overall deposition fraction. 


Deposition along the Tube Walls

Even though deposition along the tube walls was mostly attributed to complex secondary flow phenomenon and centrifugal forces, inertial impaction could not be fully ruled out. Indeed, deposition at the inner walls of the bifurcation can be attributed to inertial impaction. There was a jet flow phenomenon in the generation affected by COPD (G6-1), and hence the deposition fractions at this generation were very low as compared to its corresponding generation in the healthy geometry. The presence of an obstruction in G6-1 leads to amplified PM2.5 deposition in G7-2, as shown in Fig. 9. As a result, there was a significant asymmetry in the deposition patterns in G7-2 and G7-3 unlike in the healthy human airway geometry where deposition in G7-3 slightly exceeded that of G7-2. The second most important obstacle in a human airway bifurcation besides the carina regions were the curved sections of the airway tubes. Centrifugal forces played a key role in the axial velocity distribution at the curvatures. As such, there was a slightly amplified deposition of fine particulates on the inner tube walls immediately past the carina regions.

Fig. 9. PM2.5 deposition in healthy human airway bifurcation (G5– G8) at the (i) carina region, (ii) tube walls and (iii) whole geometry during (a) rest, (b) light activity and, (c) moderate exercise. Fig. 9. PM2.5 deposition in healthy human airway bifurcation (G5– G8) at the (i) carina region, (ii) tube walls and (iii) whole geometry during (a) rest, (b) light activity and, (c) moderate exercise.


Overall Deposition

PM2.5 deposition fractions in the healthy human geometry exceeded deposition in a COPD human airway geometry during all inhalation statuses except during moderate exercise. The total deposition fractions of PM2.5 for the healthy geometry were 1.22, 2.50 and, 1.86%; while that of the COPD geometry were 0.51, 0.91and, 2.11 for rest, light activity and moderate exercise respectively. The relative disparity in the deposition fractions of PM2.5 between the healthy and COPD case for PM2.5 was 59.19 and 33.30% during rest and light activity respectively. On the other hand, during moderate exercise, deposition fractions in the COPD geometry exceeded that in the healthy geometry by 0.37%. From the graphical representation of deposition fractions in Fig. 12, it can be concluded that the dominance of deposition mechanisms under study namely; inertial impaction, centrifugal forces and complex secondary flow phenomena varied differently for a particular diameter size at different average Reynolds numbers. The particle distribution at the inlet also played a role as particles near the airway walls tend to deposit more easily at the slightest deviation from the streamlines in the axial direction of flow. Further, particles at the rings that are closest to the middle of the cross-section deposited easily upon impaction at the carina region. For a healthy human airway geometry, deposition fractions increased as the Reynolds number increased from rest to light activity. However, as the Reynolds number of flows increased from light activity to moderate exercise, the DFs for the 3 and 6 µm particle diameters decreased slightly. This can be explained by an exaggerated jet flow phenomenon throughout the bifurcation which purged the particulate matter to the exit of the control volume at G8. The stokes numbers of the PM2.5 particles at the inlet of G5 ranged from 2.012 × 0–8 to 5.55 × 10–5. A graphical representation of the overall deposition patterns of the particles is shown in Fig. 10, whereby deposition fractions at the carina region together with overall PM2.5 deposition fractions could be closely linked to the stokes numbers of the particulate matters, as shown in Figs 10(a) and 10(c) respectively. This is however not the case for the deposition fraction in the tube section as shown in Fig. 10(b) which are mostly due to complex secondary flows and whose presence is independent of the particles’ Stokes numbers. 

Fig. 10. PM2.5 deposition in an obstructed human airway bifurcation (G5– G8) at the (i) carina region, (ii) tube walls and (iii) whole geometry during (a) rest, (b) light activity and, (c) moderate exercise.Fig. 10. PM2.5 deposition in an obstructed human airway bifurcation (G5– G8) at the (i) carina region, (ii) tube walls and (iii) whole geometry during (a) rest, (b) light activity and, (c) moderate exercise. 

Fig. 11. A scatter for the deposition fractions vs Stokes number for the (a) carina regions, (b) tube walls and, (c) whole geometry.Fig. 11. A scatter for the deposition fractions vs Stokes number for the (a) carina regions, (b) tube walls and, (c) whole geometry. 

Fig. 12. Overall deposition fractions in (a) healthy human airway geometry and (b) human airway geometry affected by COPD.Fig. 12. Overall deposition fractions in (a) healthy human airway geometry and (b) human airway geometry affected by COPD.


CONCLUSIONS


CFD was applied to investigate airflow and localized deposition of fine particulate matters in healthy and obstructed human airways. Herein, complex physical phenomena were broken down and derived from otherwise experimentally inaccessible regions. Wall shear stresses in COPD geometry were found to exceed those in the healthy human airway geometry. The carina regions were responsible for 64–88% of overall PM2.5 deposition. The most important hotspot for the deposition fine particulate matters in this study was found to be carina 5. Jet flow at the constriction affected by COPD led to amplified DFs at the succeeding carina region. The presence of an obstruction in G6-1 leads to a more amplified PM2.5 deposition in G7-2 for the COPD as compared to the corresponding generation in the healthy human airway geometry. The intensity of the stagnation and recirculation zones increased with an increase in Reynolds number from 164.3–606.4. Overall PM2.5 deposition in the carina regions as well as the curved sections was due to inertial impaction, complex secondary flow phenomena and centrifugal forces. PM2.5 deposition in the human airway affected by COPD was skewed. Overall, the DFs of PM2.5 increased with an increase in Stokes number. Total deposition fraction of PM2.5 in a COPD geometry exceeded that of a healthy geometry during the moderate exercise inhalation status only. However, the direct correlation between DFs and stokes number was lacking in a few instances due to very low ranges of Stokes number.


ACKNOWLEDGMENTS


The authors acknowledge financial support from the Ministry of Science and Technology Taiwan, ROC, under the grant numbers MOST 108-2622-E-006-017-CC1 and 109-3116-F-006-016-CC1 for this research.


REFERENCES


  1. Augusto, L., Lopes, G. and Goncalves, J. (2016). A CFD study of deposition of pharmaceutical aerosols under different respiratory conditions. Braz. J. Chem. Eng. 33: 549–558. [Publisher Site]

  2. Brook, R.D., Rajagopalan, S., Pope, C.A., Brook, J.R., Bhatnagar, A., Diez-Roux, A.V., Holguin, F., Hong, Y., Luepker, R.V. and Mittleman, M.A. (2010). Particulate matter air pollution and cardiovascular disease: An update to the scientific statement from the American Heart Association. Circulation 121: 2331–2378. [Publisher Site]

  3. Chao, H.R., Hsu, J.W., Ku, H.Y., Wang, S.L., Huang, H.B., Liou, S.H. and Tsou, T.C. (2018). Inflammatory response and PM2.5 exposure of urban traffic conductors. Aerosol Air Qual. Res. 18: 2633–2642. [Publisher Site]

  4. Chen, S.J., Liao, S.H., Jian, W.J. and Lin, C.C. (1997). Particle size distribution of aerosol carbons in ambient Air. Environ. Int. 23: 475–488. [Publisher Site]

  5. Chen, W.H. (2001). Dynamics of Sulfur dioxide absorption in a raindrop falling at terminal velocity. Atmos. Environ. 35: 4777–4790. [Publisher Site]

  6. Chen, W.H. (2002). An analysis of gas absorption by a liquid aerosol in a stationary environment. Atmos. Environ. 36: 3671–3683. [Publisher Site]

  7. Chen, W.H., Tsai, M.H. and Hung, C.I. (2014). Characterization of transient CO2 transport in two convecting aerosol droplets in tandem. Aerosol Air Qual. Res. 14: 207–219. [Publisher Site]

  8. Chen, W.H., Lee, K.H., Mutuku, J.K. and Hwang, C.J. (2018a). Flow dynamics and PM2.5 deposition in healthy and asthmatic airways at different inhalation statuses. Aerosol Air Qual. Res. 18: 866–883. [Publisher Site]

  9. Chen, X., Zhong, W., Sun, B., Jin, B. and Zhou, X. (2012). Study on gas/solid flow in an obstructed pulmonary airway with transient flow based on CFD–DPM approach. Powder Technol. 217: 252–260. [Publisher Site]

  10. Chen, X., Feng, Y., Zhong, W., Sun, B. and Tao, F. (2018b). Numerical investigation of particle deposition in a triple bifurcation airway due to gravitational sedimentation and inertial impaction. Powder Technol. 323: 284–293. [Publisher Site]

  11. Cheng, Y.S., Zhou, Y. and Chen, B.T. (1999). Particle deposition in a cast of human oral airways. Aerosol Sci. Technol. 31: 286–300. [Publisher Site]

  12. Cohen, A.J., Ross Anderson, H., Ostro, B., Pandey, K.D., Krzyzanowski, M., Künzli, N., Gutschmidt, K., Pope, A., Romieu, I. and Samet, J.M. (2005). The global burden of disease due to outdoor air pollution. J. Toxicol. Environ. Health Part A 68: 1301–1307. [Publisher Site]

  13. Comer, J., Kleinstreuer, C., Hyun, S. and Kim, C. (2000). Aerosol transport and deposition in sequentially bifurcating airways. J. Biomech. Eng. 122: 152–158. [Publisher Site]

  14. Daigle, C.C., Chalupa, D.C., Gibb, F.R., Morrow, P.E., Oberdörster, G., Utell, M.J. and Frampton, M.W. (2003). Ultrafine particle deposition in humans during rest and exercise. Inhalation Toxicol. 15: 539–552. [Publisher Site]

  15. Darquenne, C. (2012). Aerosol deposition in health and disease. J. Aerosol Med. Pulm Drug Delivery 25: 140–147. [Publisher Site]

  16. Finlay, W.H. and Martin, A.R. (2008). Recent Advances in predictive understanding of respiratory tract deposition. J. Aerosol Med. Pulm Drug Delivery 21: 189–206. [Publisher Site]

  17. Gronski, M.A., Kinchen, J.M., Juncadella, I.J., Franc, N.C. and Ravichandran, K.S. (2009). An essential role for calcium flux in phagocytes for apoptotic cell engulfment and the anti-inflammatory response. Cell Death Differ. 16: 1323–1331. [Publisher Site]

  18. Guha, A., Pradhan, K. and Halder, P.K. (2016). Finding order in complexity: A study of the fluid dynamics in a three-dimensional branching network. Phys. Fluids 28: 123602. [Publisher Site]
  19. He, K., Yang, F., Ma, Y., Zhang, Q., Yao, X., Chan, C.K., Cadle, S., Chan, T. and Mulawa, P. (2001). The characteristics of PM2.5 in Beijing, China. Atmos. Environ. 35: 4959–4970. [Publisher Site]

  20. Horsfield, K., Dart, G., Olson, D.E., Filley, G.F. and Cumming, G. (1971). Models of the Human bronchial tree. J. Appl. Phycol. 31: 207–217. [Publisher Site]

  21. Huang, J. and Zhang, L. (2011). Numerical simulation of micro-particle deposition in a realistic human upper respiratory tract model during transient breathing cycle. Particuology 9: 424–431. [Publisher Site]

  22. Hughes, J., Hoppin Jr, F. and Mead, J. (1972). Effect of lung inflation on bronchial length and diameter in excised lungs. J. Appl. Phycol. 32: 25–35. [Publisher Site]

  23. Inthavong, K., Choi, L.T., Tu, J., Ding, S. and Thien, F. (2010). Micron particle deposition in a tracheobronchial airway model under different breathing conditions. Med. Eng. Phys. 32: 1198–1212. [Publisher Site]

  24. Jain, V.V., Abejie, B., Bashir, M.H., Tyner, T. and Vempilly, J. (2013). Lung volume abnormalities and its correlation to spirometric and demographic variables in adult asthma. J. Asthma 50: 600–605. [Publisher Site]
  25. Jiang, N., Liu, X., Wang, S., Yu, X., Yin, S., Duan, S., Wang, S., Zhang, R. and Li, S. (2019). Pollution characterization, source identification, and health risks of atmospheric-particle-bound heavy metals in PM10 and PM2.5 at multiple sites in an emerging megacity in the central region of China. Aerosol Air Qual. Res. 19: 247–271. [Publisher Site]

  26. Johnson, P.C. and Jackson, R. (1987). Frictional–collisional constitutive relations for granular materials, with application to plane shearing. J. Fluid Mech. 176: 67–93. [Publisher Site]

  27. Kang, M.Y., Hwang, J. and Lee, J.W. (2011). Effect of geometric variations on pressure loss for a model bifurcation of the human lung airway. J. Biomech. 44: 1196–1199. [Publisher Site]
  28. Kim, C.S. and Iglesias, A.J. (1989). Deposition of inhaled particles in bifurcating airway Models I. inspiratory deposition. J. Aerosol Med. 2: 1–14. [Publisher Site]

  29. Kleinstreuer, C. and Zhang, Z. (2010). Airflow and particle transport in the human respiratory system. Annu. Rev. Fluid Mech. 42: 301–334. [Publisher Site]

  30. Kolanjiyil, A.V. and Kleinstreuer, C. (2017). Computational analysis of aerosol-dynamics in a human whole-lung airway model. J. Aerosol Sci. 114: 301–316. [Publisher Site]

  31. Lambert, A.R., O'shaughnessy, P.T., Tawhai, M.H., Hoffman, E.A. and Lin, C.L. (2011). Regional deposition of particles in an image-based airway model: Large-eddy simulation and left-right lung ventilation asymmetry. Aerosol Sci. Technol. 45: 11–25. [Publisher Site]

  32. Lelieveld, J., Evans, J.S., Fnais, M., Giannadaki, D. and Pozzer, A. (2015). The Contribution of outdoor air pollution sources to premature mortality on a global scale. Nature 525: 367–371. [Publisher Site]

  33. Li, Y., Wang, P., Hu, C., Wang, K., Chang, Q., Liu, L., Han, Z., Shao, Y., Zhai, Y., Zuo, Z., Gong, Z. and Wu, Y. (2018). Protein corona of airborne nanoscale PM2.5 induces aberrant proliferation of human lung fibroblasts based on a 3D organotypic culture. Sci. Rep. 8: 1939. [Publisher Site]

  34. Lin, C.H., Wu, Y.L., Lai, C.H., Watson, J.G. and Chow, J.C. (2008). Air quality measurements from the southern particulate matter supersite in Taiwan. Aerosol Air Qual. Res. 8: 233–264. [Publisher Site]

  35. Liu, Z., Hu, B., Ji, D., Wang, Y., Wang, M. and Wang, Y. (2015). Diurnal and seasonal variation of the pm2.5 apparent particle density in Beijing, China. Atmos. Environ. 120: 328–338. [Publisher Site]

  36. Luo, H., Liu, Y. and Yang, X. (2007). Particle deposition in obstructed airways. J. Biomech. 40: 3096–3104. [Publisher Site]

  37. Morel, C. (2015). Mathematical Modeling of Disperse Two-Phase Flows. Springer, New York. [Publisher Site]

  38. Mutuku, J.K. and Chen, W.H. (2018). Flow characterization in healthy airways and airways with Chronic Obstructive Pulmonary Disease (COPD) during different inhalation conditions. Aerosol Air Qual. Res. 18: 2680–2694. [Publisher Site]

  39. Nazridoust, K. and Asgharian, B. (2008). Unsteady-state airflow and particle deposition in a three-generation human lung geometry. Inhalation Toxicol. 20: 595–610. [Publisher Site]

  40. Ni, L., Chuang, C.C. and Zuo, L. (2015). Fine particulate matter in acute exacerbation of COPD. Front Physiol. 6: 294. [Publisher Site]

  41. Pope III, C.A., Ezzati, M. and Dockery, D.W. (2009). Fine-particulate air pollution and life expectancy in the United States. N. Engl. J. Med. 360: 376–386. [Publisher Site]

  42. Rahimi-Gorji, M., Pourmehran, O., Gorji-Bandpy, M. and Gorji, T. (2015). CFD simulation of airflow behavior and particle transport and deposition in different breathing conditions through the realistic model of human airways. J. Mol. Liq. 209: 121–133. [Publisher Site]

  43. Rao, S., Chirkov, V., Dentener, F., Van Dingenen, R., Pachauri, S., Purohit, P., Amann, M., Heyes, C., Kinney, P., Kolp, P., Klimont, Z., Riahi, K. and Schoepp, W. (2012). Environmental modeling and methods for estimation of the global health impacts of air pollution. Environ. Model. Assess. 17: 613–622. [Publisher Site]

  44. Rees, S.L., Robinson, A.L., Khlystov, A., Stanier, C.O. and Pandis, S.N. (2004). Mass balance closure and the federal reference method for PM2.5 in Pittsburgh, Pennsylvania. Atmos. Environ. 38: 3305–3318. [Publisher Site]

  45. Soriano, J.B., Abajobir, A.A., Abate, K.H., Abera, S.F., Agrawal, A., Ahmed, M.B., Aichour, A.N., Aichour, I., Aichour, M.T.E. and Alam, K. (2017). Global, regional, and national deaths, prevalence, disability-adjusted life years, and years lived with disability for chronic obstructive pulmonary disease and asthma, 1990-2015: A systematic analysis for the Global Burden of Disease Study 2015. Lancet Respir. Med. 5: 691–706. [Publisher Site]

  46. Sul, B., Wallqvist, A., Morris, M.J., Reifman, J. and Rakesh, V. (2014). A computational study of the respiratory airflow Characteristics in normal and obstructed Human airways. Comput. Biol. Med. 52: 130–143. [Publisher Site]

  47. Tena, A., Francos, J., Alvarez, E. and Casan, P. (2015). A three dimensional in SILICO model for the simulation of inspiratory and expiratory airflow in humans. Eng. Appl. Comput. Fluid Mech. 9: 187–198. [Publisher Site]

  48. Tsuda, A., Henry, F.S. and Butler, J.P. (2008). Gas and aerosol mixing in the acinus. Respir. Physiol. Neurobiol. 163: 139–149. [Publisher Site]

  49. Valavanidis, A., Fiotakis, K. and Vlachogianni, T. (2008). Airborne particulate matter and human health: Toxicological assessment and importance of size and composition of particles for oxidative damage and carcinogenic mechanisms. J. Environ. Sci. Health., Part C 26: 339–362. [Publisher Site]

  50. Verbanck, S., Kalsi, H., Biddiscombe, M., Agnihotri, V., Belkassem, B., Lacor, C. and Usmani, O. (2011). Inspiratory and expiratory aerosol deposition in the upper airway. Inhalation Toxicol. 23: 104–111. [Publisher Site]

  51. Viegas, C.A., Ferrer, A., Montserrat, J.M., Barbera, J.A., Roca, J. and Rodriguez-Roisin, R. (1996). Ventilation-perfusion response after fenoterol in hypoxemic patients with stable COPD. Chest 110: 71–77. [Publisher Site]

  52. Weibel, E.R. (1963). Geometric and Dimensional Airway Models of Conductive, Transitory and Respiratory Zones of the Human Lung, In Morphometry of the Human Lung, Springer Berlin Heidelberg, Berlin, Heidelberg, pp. 136–142. [Publisher Site]

  53. Weibel, E.R. (1965). Morphometry of the human lung. The Anesthesiology 26: 367–367. [Publisher Site]

  54. Wessel, R.A. and Righi, J. (1988). Generalized correlations for inertial impaction of particles on a circular cylinder. Aerosol Sci. Technol. 9: 29–60. [Publisher Site]

  55. Yanai, M., Sekizawa, K., Ohrui, T., Sasaki, H. and Takishima, T. (1992). Site of airway obstruction in pulmonary disease: Direct measurement of intrabronchial pressure. J. Appl. Phycol. 72: 1016–1023. [Publisher Site]

  56. Zhang, H. and Papadakis, G. (2010). Computational analysis of flow structure and particle deposition in a single asthmatic human airway bifurcation. J. Biomech. 43: 2453–2459. [Publisher Site]

  57. Zhang, P., Duan, J., Chen, G. and Wang, W. (2019). Numerical investigation on gas-solid flow in a circumfluent cyclone separator. Aerosol Air Qual. Res. 19: 971–980. [Publisher Site]

  58. Zhang, X., Kang, J., Chen, H., Yao, M. and Wang, J. (2018). PM2.5 meets blood: In vivo damages and immune defense. Aerosol Air Qual. Res. 18: 456–470. [Publisher Site]

  59. Zhang, Z. and Kleinstreuer, C. (2001). Effect of particle inlet distributions on deposition in a triple bifurcation lung airway model. J. Aerosol Med. 14: 13–29. [Publisher Site]

  60. Zhang, Z., Kleinstreuer, C. and Kim, C. (2001). Flow structure and particle transport in a triple bifurcation airway model. J. Fluids Eng. 123: 320–330. [Publisher Site]

  61. Zhang, Z. and Kleinstreuer, C. (2002). Transient Airflow structures and particle transport in a sequentially branching lung airway model. Phys. Fluids 14: 862–880. [Publisher Site]

  62. Zhang, Z., Kleinstreuer, C. and Kim, C. (2002). Gas–solid two-phase flow in a triple bifurcation lung airway model. Int. J. Multiphase Flow 28: 1021–1046. [Publisher Site]

Latest Articles