1. Introduction

Rhinitis and pulmonary pathology frequently occur together. The epidemiologic, pathophysiologic, and clinical data are so compelling that the concept of "one airway, one disease" is accepted [1]. Several mechanisms have been held responsible for the interaction between the upper and lower airways (see Review [2]). A good example is the high prevalence of chronic rhinosinusitis in patients with asthma. Recent data provide ample evidence of existence of a systemic pathway between allergic rhinitis and asthma. It has been found that the severity of asthma directly correlates with the severity of sinus disease [3, 4]. Patients with chronic rhinosinusitis and asthma constitute the most severe form of unified respiratory tract disease, which is characterized by older age of the patients, greater duration of nasal symptoms, extent of sinus radiological changes, more prominent systemic inflammation markers, greater bronchial obstruction, and incidence of perennial allergic rhinitis [5].

Treating asthma with inhaled steroid therapy appears to be important in treating refractory chronic rhinosinusitis associated with asthma. In [6], the results of topical inhalation treatment of acute bacterial rhinosinusitis and standard systemic antibiotic therapy are compared. It has been found that topical inhalation therapy of acute bacterial rhinosinusitis may provide better treatment options, because systemic antibiotics can be associated with different adverse effects.

change of air velocity follows from the requirement of conservation airflow before

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance

The trachea (0th generation) has D(0) = 1.44 cm and V(0) = 290 cm/s. The branches of the 18th generation correspond to terminal bronchioles. The radius R(i) = D(i)/2, and length L(i) of each branch is assumed to be equal to three times

In the adopted model, the time tp for the passage of the air channel is the same

We can compare this time with the characteristic inhalation time tin = 1.25 s (we assume that the breathing frequency is 12 1/min, and for each period, there are four identical phase durations: inhale, pause, exhale, and pause). The ratio tin/tp ≈ 100 can be considered as the possible number of TBT generations that air could pass during inhalation. As was to be expected, this number is significantly greater than the total number of generations (≈ 25) of TBT. The passage of the airways to the

An important dimensionless quantity in air mechanics is the Reynolds number (Re). It is used in the scaling of similar but different sized flow situations [11]. The value of Re is the ratio of inertial forces to viscous forces within an air which is subjected to relative internal movement due to different air velocities, which is known as a boundary layer in the case of a bounding surface such as the interior of an airway. In the air channel of TBT, the Reynolds number can be defined as

where ν is the air kinematic viscosity. The Reynolds number is important for different situations where a fluid is in relative motion to a surface. Matching the Reynolds numbers in different airways is sufficient to guarantee similitude of airflows. Table 1 shows the Reynolds numbers characterizing the air flow in the

It can be seen that in almost all branches of the TBT, the flow is laminar,

number (Stk). It is defined as the ratio of the characteristic time of an aerosol particle to a characteristic time of the flow. In the air channel of TBT, the Stokes

<sup>∗</sup> dp

i= 0 2 4 6 8 10 12 14 16 18 20 Re = 1392.0 552.4 219.2 87.0 34.5 13.7 5.4 2.2 0.9 0.3 0.1

Stk ¼ ð Þ 1=18 ∗ ρp=ρ<sup>a</sup>

The behavior of particles suspended in an air flow is characterized [11] by Stokes

2

level of the alveoli takes no more than 10% of the total inhalation time.

tp ¼ L ið Þ=V iðÞ¼ 3 ∗ D ið Þ=V iðÞ¼ 3 ∗ D 0ð Þ=V 0ð Þ ≈0:015 sec : (3)

V iðÞ¼ V 0ð Þ <sup>∗</sup> <sup>2</sup>�i=<sup>3</sup> (2)

Re iðÞ¼ V ið Þ ∗ R ið Þ=ν (4)

<sup>=</sup><sup>ν</sup> <sup>∗</sup> ½ � V ið Þ=R ið Þ (5)

and after dividing branches of TBT:

DOI: http://dx.doi.org/10.5772/intechopen.84254

the diameter.

for all generations:

2.1 Air flow characteristics

respiratory tract.

Table 1.

103

however with Re > 1.

number can be defined as

Reynolds numbers in the branches of the TBT.

Aerosol in lungs are the task the solution of which is necessary in various areas. For example, in medical applications, aerosol therapy allows the medication to be delivered directly to those parts of the lungs where it should act.

In this chapter, we consider both deposition and clearance kinetics throughout the respiratory tract. Aerosol particles can penetrate the human lungs deep into the smallest respiratory tracts and parts of the lungs during inhalation. The effects, caused by inhaled particles, depend on the site at which they deposit within the respiratory system. Information concerning particle deposition in human lungs is of special interest because it is useful for a quantitative risk evaluation from exposure to aerosol particles. It is also vital to the effective administration of pharmaceutical aerosols by inhalation to targeted regions of the respiratory tracts. The deposition of particles in the lungs for given exposure conditions depends on parameters such as the airflow field, aerosol particle properties, breathing pattern, and geometric airway characteristics. There have been many investigations on this subject. Some recent reviews have been carried out and reported in [7–9]. In Section 3 of this chapter, we describe the aerosol inertial capture processes. Particles move along curvilinear trajectories in self-twisted vortex air flows. Deposition occurs over several cycles of inhalation and exhalation. The result is the distribution of deposition of inhaled aerosol particles in the airway generations in the respiratory tract of the lung.

Once deposited, specific clearance mechanisms will become effective, leading to a clearance of the deposited substances or a retention in different compartments of the respiratory tract. The main clearance mechanism for insoluble particles deposited in the conducting airways is via the mucociliary escalator. For the calculation of retention curves, this bronchial clearance model was applied to predicted deposition patterns for different particle sizes and breathing conditions. The importance of taking into account the temporal characteristics of the mucociliary escalation of aerosol particles deposited in the lungs has been demonstrated.

#### 2. Model of human respiratory system

Morphometric models were initiated by E. Weibel [10] who assigned lengths and diameters to a symmetrical tracheobronchial tree (TBT) based on measurements of human airway casts. Weibel's airway tree model assumes dichotomous branching and creates a fractal structure, which is physiologically consistent in the amount of air delivered to the terminal airways and the spatial arrangement of branches within the lung. A generation is defined as all airways occurring at a specific number of bifurcations from the trachea.

The number N (i) of the branches of the ith generation and the diameter D (i) are given by the formulas

$$\mathbf{N(i)} = \mathbf{2^i} \text{; } \mathbf{D(i)} = \mathbf{D(0)} \ast \mathbf{2^{-i/3}} \tag{1}$$

The airway diameter decreases distally, but because of the increasing number of tubes, the total cross section increases and the air velocity decreases. The law of

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance DOI: http://dx.doi.org/10.5772/intechopen.84254

change of air velocity follows from the requirement of conservation airflow before and after dividing branches of TBT:

$$\mathbf{V}(\mathbf{i}) = \mathbf{V}(\mathbf{0}) \* \mathbf{2}^{-i/3} \tag{2}$$

The trachea (0th generation) has D(0) = 1.44 cm and V(0) = 290 cm/s. The branches of the 18th generation correspond to terminal bronchioles. The radius R(i) = D(i)/2, and length L(i) of each branch is assumed to be equal to three times the diameter.

In the adopted model, the time tp for the passage of the air channel is the same for all generations:

$$\mathbf{t\_p} = \mathbf{L(i)}/\mathbf{V(i)} = \mathbf{3} \ast \mathbf{D(i)}/\mathbf{V(i)} = \mathbf{3} \ast \mathbf{D(0)}/\mathbf{V(0)} \approx \mathbf{0.015} \text{ sec.}\tag{3}$$

We can compare this time with the characteristic inhalation time tin = 1.25 s (we assume that the breathing frequency is 12 1/min, and for each period, there are four identical phase durations: inhale, pause, exhale, and pause). The ratio tin/tp ≈ 100 can be considered as the possible number of TBT generations that air could pass during inhalation. As was to be expected, this number is significantly greater than the total number of generations (≈ 25) of TBT. The passage of the airways to the level of the alveoli takes no more than 10% of the total inhalation time.

#### 2.1 Air flow characteristics

Treating asthma with inhaled steroid therapy appears to be important in treating refractory chronic rhinosinusitis associated with asthma. In [6], the results of topical inhalation treatment of acute bacterial rhinosinusitis and standard systemic antibiotic therapy are compared. It has been found that topical inhalation therapy of acute bacterial rhinosinusitis may provide better treatment options, because sys-

Aerosol in lungs are the task the solution of which is necessary in various areas. For example, in medical applications, aerosol therapy allows the medication to be

In this chapter, we consider both deposition and clearance kinetics throughout the respiratory tract. Aerosol particles can penetrate the human lungs deep into the smallest respiratory tracts and parts of the lungs during inhalation. The effects, caused by inhaled particles, depend on the site at which they deposit within the respiratory system. Information concerning particle deposition in human lungs is of special interest because it is useful for a quantitative risk evaluation from exposure to aerosol particles. It is also vital to the effective administration of pharmaceutical aerosols by inhalation to targeted regions of the respiratory tracts. The deposition of particles in the lungs for given exposure conditions depends on parameters such as the airflow field, aerosol particle properties, breathing pattern, and geometric airway characteristics. There have been many investigations on this subject. Some recent reviews have been carried out and reported in [7–9]. In Section 3 of this chapter, we describe the aerosol inertial capture processes. Particles move along curvilinear trajectories in self-twisted vortex air flows. Deposition occurs over several cycles of inhalation and exhalation. The result is the distribution of deposition of inhaled aerosol particles in the airway generations in the respiratory tract of the

Once deposited, specific clearance mechanisms will become effective, leading to a clearance of the deposited substances or a retention in different compartments of the respiratory tract. The main clearance mechanism for insoluble particles deposited in the conducting airways is via the mucociliary escalator. For the calculation of retention curves, this bronchial clearance model was applied to predicted deposition patterns for different particle sizes and breathing conditions. The importance of taking into account the temporal characteristics of the mucociliary escalation of

Morphometric models were initiated by E. Weibel [10] who assigned lengths and diameters to a symmetrical tracheobronchial tree (TBT) based on measurements of human airway casts. Weibel's airway tree model assumes dichotomous branching and creates a fractal structure, which is physiologically consistent in the amount of air delivered to the terminal airways and the spatial arrangement of branches within the lung. A generation is defined as all airways occurring at a

The number N (i) of the branches of the ith generation and the diameter D (i)

The airway diameter decreases distally, but because of the increasing number of tubes, the total cross section increases and the air velocity decreases. The law of

N iðÞ¼ 2i ; D iðÞ¼ D 0ð Þ <sup>∗</sup> <sup>2</sup>�i=<sup>3</sup> (1)

temic antibiotics can be associated with different adverse effects.

delivered directly to those parts of the lungs where it should act.

aerosol particles deposited in the lungs has been demonstrated.

2. Model of human respiratory system

specific number of bifurcations from the trachea.

are given by the formulas

102

lung.

Rhinosinusitis

An important dimensionless quantity in air mechanics is the Reynolds number (Re). It is used in the scaling of similar but different sized flow situations [11]. The value of Re is the ratio of inertial forces to viscous forces within an air which is subjected to relative internal movement due to different air velocities, which is known as a boundary layer in the case of a bounding surface such as the interior of an airway. In the air channel of TBT, the Reynolds number can be defined as

$$\text{Re}(\mathbf{i}) = \mathbf{V}(\mathbf{i}) \* \mathbf{R}(\mathbf{i}) / \nu \tag{4}$$

where ν is the air kinematic viscosity. The Reynolds number is important for different situations where a fluid is in relative motion to a surface. Matching the Reynolds numbers in different airways is sufficient to guarantee similitude of airflows. Table 1 shows the Reynolds numbers characterizing the air flow in the respiratory tract.

It can be seen that in almost all branches of the TBT, the flow is laminar, however with Re > 1.

The behavior of particles suspended in an air flow is characterized [11] by Stokes number (Stk). It is defined as the ratio of the characteristic time of an aerosol particle to a characteristic time of the flow. In the air channel of TBT, the Stokes number can be defined as


Stk ¼ ð Þ 1=18 ∗ ρp=ρ<sup>a</sup> <sup>∗</sup> dp 2 <sup>=</sup><sup>ν</sup> <sup>∗</sup> ½ � V ið Þ=R ið Þ (5)

#### Table 1.

Reynolds numbers in the branches of the TBT.


Table 2.

Stokes numbers for particles of different diameters.

where ρ<sup>p</sup> and ρ<sup>a</sup> are the particle and air densities and dp is the particle diameter. Value R(i)/V(i) is the characteristic time of the flow. A particle with a low Stokes number follows fluid streamlines, while a particle with a large Stokes number is dominated by its inertia and continues along its initial trajectory. For Stk ≫ 1, particles will detach from a flow especially where the flow decelerates abruptly. For Stk ≪ 1, particles follow fluid streamlines closely.

In the accepted TBT model, the ratio of air velocity to diameter is the same for all generations; therefore, the Stokes number is the same in all branches, and it changes depending on the particle diameter. The corresponding results are shown in Table 2.

It can be seen that particles with a diameter of up to 10 μm (Stk ≪ 1) are involved in the movement of air. Only particles with a large diameter can "break away" from the flow and, for example, settle onto the wall directly due to the inertial mechanism.

#### 2.2 Mucociliary escalator

The kinetics of mechanical particle transport from the TBT via the mucociliary escalator to the larynx cannot be measured directly. Average tracheal mucus velocities in healthy nonsmokers, as measured by noninvasive radiological techniques, ranged from 4 to 6 mm/min consequently; both ICRP [12] and NCRP [13] committees adopted a value of 5.5 mm/min. Because of the very limited experimental data on mucus clearance velocities in human bronchial airways, mucus velocity in other bronchial airways has to be calculated, using assumptions about TBT morphology and properties of the mucus flow in these airways. The law of change of mucus velocity U(i) follows from the requirement of conservation flows before and after dividing branches of TBT:

$$\mathbf{2} \ast \pi \ast \mathbf{D}(\mathbf{i}) \ast \mathbf{U}(\mathbf{i}) = \pi \ast \mathbf{D}(\mathbf{i} - \mathbf{1}) \ast \mathbf{U}(\mathbf{i} - \mathbf{1}) \tag{6}$$

The diameter of each airway is determined by (1). So, we have

$$\mathbf{U}(\mathbf{i}) = \mathbf{U}(\mathbf{0}) \* 2 - 2\mathbf{i}/3 \tag{7}$$

wall move more slowly due to the effect of deceleration of the air flow in the boundary layer. During inhalation, they may not reach the end of the channel. They

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance

center. A cylindrical coordinate system (r,z) is used. The airway radius is R; length is L.

Trajectories of particles entering the air flow in a straight pipe with round cross section at distances rс from the

P iðÞ¼ <sup>π</sup>R<sup>2</sup> � <sup>π</sup>r<sup>c</sup>

An important circumstance is the air swirling flow at the point of bifurcation of the branches of TBT [7]. Numerical studies on the bifurcation flow of lung models consistently identified different types of secondary flow pattern such as swirl-flow pattern or double vortex, also known as Dean flow. At inspiration, the rate of swirling flow at the inlet to the branch of TBT is 20–25% of the axial velocity. In such a flow, the centrifugal force shifts the particles to the wall. Here, the particle is additionally inhibited. After several cycles of "inhale-exhale," the particle will deposit on the wall. The capture of particles in some branch of the TBT occurs if the transit time of this branch is longer than the inhalation time tin≈ 1.25 s (see Section 1). The transit time increases with increasing distance from the channel axis. The particles that fly into the airway at a distance r > rc from axis are captured. This condition determines the probability P (i) of particle capture in the ith branch

<sup>2</sup> <sup>=</sup>πR<sup>2</sup> <sup>¼</sup> <sup>1</sup> � rc

This mechanism is the most effective even for small aerosol particles.

There are numerous research studies focusing on airflow in the respiratory system (see Review [14]). The results of these studies provide an almost complete picture of the airflow in the lungs. It is easy to recognize the "classical" flow

The axial velocity distribution at the entrance can be considered uniform. Further downstream, a Prandtl boundary layer forms along the walls. Its thickness is gradually increasing. In sufficiently long pipes at a distance Le from the entrance, the boundary layer fills the entire cross section. It establishes a Poiseuille flow with a parabolic distribution of the axial velocity along the radius r. The steady-state flow does not depend on the velocity distribution at the inlet. The distance at which the

For a small (compared to the channel radius) thickness of the boundary layer,

Poiseuille flow is formed is called the input length Le. For Le we have

the axial air velocity Vz is described by the relation (Blasius solution):

2

Le ≈0:1 ∗ D ∗ Re (9)

=R<sup>2</sup> (8)

turn back during exhalation.

DOI: http://dx.doi.org/10.5772/intechopen.84254

3.1 Air flow characteristics

patterns in smooth cylindrical tubes (see [11]).

of the TBT:

105

Figure 1.

This model uses an empirically derived expression for mucociliary transport rates in TBT.

For the calculation of retention characteristics, this bronchial clearance model was applied to predicted deposition patterns for different particle sizes and breathing conditions.

#### 3. Inertial capture

Inertial capture of aerosol particles from air flow is illustrated by Figure 1. Particles entering the channel near the center fly through it during the time of the tp ≈ 0.015 sec (see above Section 2.1.). Particles entering the channel near the

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance DOI: http://dx.doi.org/10.5772/intechopen.84254

Figure 1.

where ρ<sup>p</sup> and ρ<sup>a</sup> are the particle and air densities and dp is the particle diameter. Value R(i)/V(i) is the characteristic time of the flow. A particle with a low Stokes number follows fluid streamlines, while a particle with a large Stokes number is dominated by its inertia and continues along its initial trajectory. For Stk ≫ 1, particles will detach from a flow especially where the flow decelerates abruptly. For

dp, μm 1 2 4 6 10 20 40 Stk = 0.0015 0.0060 0.024 0.054 0.149 0.60 2.39

In the accepted TBT model, the ratio of air velocity to diameter is the same for all generations; therefore, the Stokes number is the same in all branches, and it changes depending on the particle diameter. The corresponding results are shown in

The kinetics of mechanical particle transport from the TBT via the mucociliary escalator to the larynx cannot be measured directly. Average tracheal mucus velocities in healthy nonsmokers, as measured by noninvasive radiological techniques, ranged from 4 to 6 mm/min consequently; both ICRP [12] and NCRP [13] committees adopted a value of 5.5 mm/min. Because of the very limited experimental data on mucus clearance velocities in human bronchial airways, mucus velocity in other bronchial airways has to be calculated, using assumptions about TBT morphology and properties of the mucus flow in these airways. The law of change of mucus velocity U(i) follows from the requirement of conservation flows before and after

The diameter of each airway is determined by (1). So, we have

This model uses an empirically derived expression for mucociliary transport

For the calculation of retention characteristics, this bronchial clearance model was applied to predicted deposition patterns for different particle sizes and breath-

Inertial capture of aerosol particles from air flow is illustrated by Figure 1. Particles entering the channel near the center fly through it during the time of the tp ≈ 0.015 sec (see above Section 2.1.). Particles entering the channel near the

2 ∗ π ∗ D ið Þ∗ U iðÞ¼ π ∗ D ið Þ � 1 ∗ U ið Þ � 1 (6)

U iðÞ¼ U 0ð Þ ∗ 2 � 2i=3 (7)

It can be seen that particles with a diameter of up to 10 μm (Stk ≪ 1) are involved in the movement of air. Only particles with a large diameter can "break away" from the flow and, for example, settle onto the wall directly due to the

Stk ≪ 1, particles follow fluid streamlines closely.

Stokes numbers for particles of different diameters.

Table 2.

Table 2.

Rhinosinusitis

inertial mechanism.

2.2 Mucociliary escalator

dividing branches of TBT:

rates in TBT.

ing conditions.

104

3. Inertial capture

Trajectories of particles entering the air flow in a straight pipe with round cross section at distances rс from the center. A cylindrical coordinate system (r,z) is used. The airway radius is R; length is L.

wall move more slowly due to the effect of deceleration of the air flow in the boundary layer. During inhalation, they may not reach the end of the channel. They turn back during exhalation.

An important circumstance is the air swirling flow at the point of bifurcation of the branches of TBT [7]. Numerical studies on the bifurcation flow of lung models consistently identified different types of secondary flow pattern such as swirl-flow pattern or double vortex, also known as Dean flow. At inspiration, the rate of swirling flow at the inlet to the branch of TBT is 20–25% of the axial velocity.

In such a flow, the centrifugal force shifts the particles to the wall. Here, the particle is additionally inhibited. After several cycles of "inhale-exhale," the particle will deposit on the wall. The capture of particles in some branch of the TBT occurs if the transit time of this branch is longer than the inhalation time tin≈ 1.25 s (see Section 1). The transit time increases with increasing distance from the channel axis. The particles that fly into the airway at a distance r > rc from axis are captured. This condition determines the probability P (i) of particle capture in the ith branch of the TBT:

$$P(\mathbf{i}) = \left(\pi \mathbf{R}^2 - \pi r\_\mathbf{c}^2\right) / \pi \mathbf{R}^2 = \mathbf{1} - \mathbf{r\_c}^2 / \mathbf{R}^2 \tag{8}$$

This mechanism is the most effective even for small aerosol particles.

#### 3.1 Air flow characteristics

There are numerous research studies focusing on airflow in the respiratory system (see Review [14]). The results of these studies provide an almost complete picture of the airflow in the lungs. It is easy to recognize the "classical" flow patterns in smooth cylindrical tubes (see [11]).

The axial velocity distribution at the entrance can be considered uniform. Further downstream, a Prandtl boundary layer forms along the walls. Its thickness is gradually increasing. In sufficiently long pipes at a distance Le from the entrance, the boundary layer fills the entire cross section. It establishes a Poiseuille flow with a parabolic distribution of the axial velocity along the radius r. The steady-state flow does not depend on the velocity distribution at the inlet. The distance at which the Poiseuille flow is formed is called the input length Le. For Le we have

$$\mathbf{L\_e} \approx \mathbf{0}.\mathbf{1} \ast \mathbf{D} \ast \mathbf{R\_e} \tag{9}$$

For a small (compared to the channel radius) thickness of the boundary layer, the axial air velocity Vz is described by the relation (Blasius solution):

$$\mathbf{V}\_{\mathbf{z}} = \mathbf{V}\_{\text{in}} \ast \mathbf{f} \left( \begin{smallmatrix} \zeta \\ \end{smallmatrix} \right) \tag{10}$$

particles into the zone of slow near-wall flow can significantly increase the time of flight of the particle through the channel. The definition of the transit time can be obtained, for example, by integrating the equation of motion for x (τ = t\*Vin/R). The estimate of the channel transit time τc, taking into account the removal of a

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance

As λ increases to large values, the time span τc grows exponentially and can reach values greater than the inspiratory time. This can be considered the effect of trapping a particle in the corresponding channel. According to Eq. (16), the value of τ<sup>c</sup> is different for particles entering the channel at different distances y1 from the center of the channel. Only particles are captured that enter the channel at distances

The dimensionless inspiratory time τin is defined in Section 2.1; it is ≈ 250. The right-hand side of Eq. (17) varies with the generation number and with the characteristics (mass and size) of aerosol particles. Changes in the entrance radius of the capture of particles at the channel entrance lead to changes (according to

The results of the calculation of the probability of the capture of particles of different diameters in the generation of the respiratory tract TBT are presented in

It can be seen that the air cleaning mechanism due to the capture of aerosols in the airways is quite effective even for particles with a diameter of 1–2 μm. Particles with a diameter of 4 μm and more are almost completely captured in channels with a generation number < 19, that is, to respiratory branches of TBT. Particles of large diameters (>10 μm) are captured in the upper respiratory tract, that is, the

Dependence of the probability P of the capture of a particle (ordinates) with diameter d in the respiratory tracts

Eq. (8)) of the probability of their settling in the ith generation.

=Stk ∗ Г<sup>2</sup> (16)

<sup>¼</sup> <sup>Г</sup><sup>2</sup> <sup>∗</sup> Stk <sup>∗</sup> <sup>τ</sup>in � <sup>λ</sup> =½ � exp ð Þ� <sup>λ</sup> <sup>1</sup> (17)

particle into the zone of a slow near-wall flow, has the form:

greater than that determined from the equation

DOI: http://dx.doi.org/10.5772/intechopen.84254

у1= 1 � y1

3.3 Results

Figure 2.

Figure 2.

107

of the ith generation (abscissa axis) of TBT.

τ<sup>с</sup> ¼ λ þ ½ � exp ð Þ� λ 1 ∗ у1= 1 � y1

where Vin is the velocity far from the wall, which can be equated to the velocity at the entrance to the channel, and ζ is the self-similar variable:

$$\boldsymbol{\zeta} = (\mathbf{R} - \mathbf{r}) \ast \left[\mathbf{V}\_{\text{in}} / \mathbf{z} \,\nu\right]^{1/2} \tag{11}$$

The function f (ζ) is the result of a numerical solution of the continuity and Navier–Stokes equations after the transition to the self-similar variable ζ. The function f (ζ) initially grows linearly with increasing ζ and then approaches to 1. The scale of change is ζ<sup>0</sup> ≈ 3.

The swirling of the air flows was found in the numerical simulation of currents in the TBT channels [15, 16]. The distribution of the azimuthal velocity V<sup>φ</sup> along the radius in the cross section of the channel is well interpolated by the function

$$\mathbf{V\_{op}} = \mathbf{S} \ast \mathbf{V\_{in}} \ast \left(\mathbf{r}/\mathbf{R}\right) \ast \left(\mathbf{1} - \mathbf{r}/\mathbf{R}\right) \tag{12}$$

This function describes the real decrease in azimuthal velocity as it approaches the axis and the channel wall. The intensity of the vortex (swirl number) S is defined as the ratio of the average over the cross-sectional azimuthal velocity V<sup>φ</sup> to axial veloсity Vz. At the entrance to the air channel of the TBD, the value of S = S0 ≈ 0.15–0.25 [15], and then S (z) decreases exponentially. The characteristic decay distance is close to the input length Le [17].

#### 3.2 Model of the capture of aerosol particles

This data is sufficient to make an estimate of the probability of the capture in the ith branch of the TBT. The motion of aerosol particles along the z axis is determined by the speed [Eq. (10)]. In the radial direction, the aerosol particle moves under the action of centrifugal force:

$$\mathbf{a} \ast \pi \ast \mathbf{d\_p} \ast \rho\_\mathbf{a} \ast \nu \ast (\mathbf{dr}/\mathbf{dt}) = \mathbf{m} \ast \mathbf{V\_\phi}^2/\mathbf{r} \tag{13}$$

Time does not explicitly enter into Eqs. (10) and (13). They define (r, z) trajectory. In dimensionless variables x = z/R and y = r/R, this trajectory is described by the formula

$$\mathbf{y} = \mathbf{y}\_1 / \left\{ \mathbf{y}\_1 + (\mathbf{1} - \mathbf{y}\_1) \ast \exp \left[ -(2/3) \ast \Gamma^2 \ast \mathbf{Stk} \ast \zeta\_\mathbf{o} \ast \mathbf{x}^{3/2} / \mathrm{Re}^{1/2} \right] \right\} \tag{14}$$

Here, y1 = r1/R is the dimensionless radius of the particle entrance. The full trajectory will be obtained after twisting (r, z) trajectories around the axis of the channel in accordance with the rotation of the particles with the azimuthal velocity V<sup>φ</sup> along with the flow.

From Eq. (14), it follows that if Re ≫ 1 (initial generation of TBT) and/or Stk ≪ 1 (small particle size), then у ≈ у<sup>1</sup> all the way. In this instance, the particles do not shift to the wall. The shift is noticeable in those air channels where

$$\lambda = (\mathbf{2}/\mathbf{3}) \ast \Gamma^2 \ast \mathbf{Stk} \ast \zeta\_\mathbf{o} \ast \mathbf{x}\_\mathbf{c}^{3/2} / \mathrm{Re}^{1/2} \geq \mathbf{1} \tag{15}$$

Even if λ ≪ 1, always y < 1, that is, direct deposition of particles due to their centrifugal demolition to the channel walls does not occur. However, the removal of Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance DOI: http://dx.doi.org/10.5772/intechopen.84254

particles into the zone of slow near-wall flow can significantly increase the time of flight of the particle through the channel. The definition of the transit time can be obtained, for example, by integrating the equation of motion for x (τ = t\*Vin/R). The estimate of the channel transit time τc, taking into account the removal of a particle into the zone of a slow near-wall flow, has the form:

$$\pi\_{\rm c} = \left\{ \lambda + \left[ \exp \left( \lambda \right) - \mathbf{1} \right] \* \mathbf{y}\_1 / \left( \mathbf{1} - \mathbf{y}\_1 \right) \right\} / \mathbf{Stk} \* \Gamma^2 \tag{16}$$

As λ increases to large values, the time span τc grows exponentially and can reach values greater than the inspiratory time. This can be considered the effect of trapping a particle in the corresponding channel. According to Eq. (16), the value of τ<sup>c</sup> is different for particles entering the channel at different distances y1 from the center of the channel. Only particles are captured that enter the channel at distances greater than that determined from the equation

$$\mathbf{y}\_1/\left(\mathbf{1} - \mathbf{y}\_1\right) = \left[\Gamma^2 \ast \mathbf{Stk} \ast \tau\_{\mathrm{in}} - \lambda\right]/\left[\exp\left(\lambda\right) - \mathbf{1}\right] \tag{17}$$

The dimensionless inspiratory time τin is defined in Section 2.1; it is ≈ 250.

The right-hand side of Eq. (17) varies with the generation number and with the characteristics (mass and size) of aerosol particles. Changes in the entrance radius of the capture of particles at the channel entrance lead to changes (according to Eq. (8)) of the probability of their settling in the ith generation.

#### 3.3 Results

Vz ¼ Vin ∗ f ð Þζ (10)

V<sup>φ</sup> ¼ S ∗ Vin ∗ ð Þ r=R ∗ ð Þ 1 � r=R (12)

2

=r (13)

=Re<sup>1</sup>=<sup>2</sup> . 1 (15)

(14)

<sup>1</sup>=<sup>2</sup> (11)

where Vin is the velocity far from the wall, which can be equated to the velocity

ζ ¼ ð Þ R � r ∗ ½ � Vin=z ν

The function f (ζ) is the result of a numerical solution of the continuity and Navier–Stokes equations after the transition to the self-similar variable ζ. The function f (ζ) initially grows linearly with increasing ζ and then approaches to 1. The

The swirling of the air flows was found in the numerical simulation of currents in the TBT channels [15, 16]. The distribution of the azimuthal velocity V<sup>φ</sup> along the radius in the cross section of the channel is well interpolated by the function

This function describes the real decrease in azimuthal velocity as it approaches

This data is sufficient to make an estimate of the probability of the capture in the ith branch of the TBT. The motion of aerosol particles along the z axis is determined by the speed [Eq. (10)]. In the radial direction, the aerosol particle moves under the

3 ∗ π ∗ dp ∗ ρ<sup>a</sup> ∗ ν ∗ ð Þ¼ dr=dt m ∗ V<sup>φ</sup>

� � <sup>∗</sup> exp : �ð Þ <sup>2</sup>=<sup>3</sup> <sup>∗</sup> <sup>Г</sup><sup>2</sup> <sup>∗</sup> Stk <sup>∗</sup> <sup>ζ</sup><sup>o</sup> <sup>∗</sup> x3<sup>=</sup><sup>2</sup> <sup>=</sup>Re<sup>1</sup>=<sup>2</sup> n o h i

Here, y1 = r1/R is the dimensionless radius of the particle entrance. The full trajectory will be obtained after twisting (r, z) trajectories around the axis of the channel in accordance with the rotation of the particles with the azimuthal velocity

shift to the wall. The shift is noticeable in those air channels where

<sup>λ</sup> <sup>¼</sup> ð Þ <sup>2</sup>=<sup>3</sup> <sup>∗</sup> <sup>Г</sup><sup>2</sup> <sup>∗</sup> Stk <sup>∗</sup> <sup>ζ</sup><sup>o</sup> <sup>∗</sup> <sup>x</sup><sup>с</sup>

From Eq. (14), it follows that if Re ≫ 1 (initial generation of TBT) and/or Stk ≪ 1 (small particle size), then у ≈ у<sup>1</sup> all the way. In this instance, the particles do not

Even if λ ≪ 1, always y < 1, that is, direct deposition of particles due to their centrifugal demolition to the channel walls does not occur. However, the removal of

3=2

Time does not explicitly enter into Eqs. (10) and (13). They define (r, z) trajectory. In dimensionless variables x = z/R and y = r/R, this trajectory is described by

the axis and the channel wall. The intensity of the vortex (swirl number) S is defined as the ratio of the average over the cross-sectional azimuthal velocity V<sup>φ</sup> to axial veloсity Vz. At the entrance to the air channel of the TBD, the value of S = S0 ≈ 0.15–0.25 [15], and then S (z) decreases exponentially. The characteristic

decay distance is close to the input length Le [17].

3.2 Model of the capture of aerosol particles

at the entrance to the channel, and ζ is the self-similar variable:

scale of change is ζ<sup>0</sup> ≈ 3.

Rhinosinusitis

action of centrifugal force:

y ¼ y1= y1 þ 1 � y1

V<sup>φ</sup> along with the flow.

the formula

106

The results of the calculation of the probability of the capture of particles of different diameters in the generation of the respiratory tract TBT are presented in Figure 2.

It can be seen that the air cleaning mechanism due to the capture of aerosols in the airways is quite effective even for particles with a diameter of 1–2 μm. Particles with a diameter of 4 μm and more are almost completely captured in channels with a generation number < 19, that is, to respiratory branches of TBT. Particles of large diameters (>10 μm) are captured in the upper respiratory tract, that is, the

#### Figure 2.

Dependence of the probability P of the capture of a particle (ordinates) with diameter d in the respiratory tracts of the ith generation (abscissa axis) of TBT.

Figure 3. The distribution of trapped particles of different diameters d by the numbers of generation of TBT.

trachea, the zonal extrapulmonary bronchi, and the intrapulmonary subsegmental bronchi.

The probability of capture (formula Eq.(8)) determines the proportion of captured particles from the total number of particles entering the channel. It is possible to determine the probability Pa (i) of the particles reaching the ith generation:

$$Pa(i) = \prod\_{k=0}^{k=i-1} \left[1 - P(k)\right] \tag{18}$$

for the initial estimates, it is advisable to limit ourselves to the simpler initial

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance

We introduce adequate variables to describe the dynamics of mucociliary cleaning of the lungs. We will describe the flow j (z, t) of the mucociliary escalation by the product of the velocity U(z) and the density n (z, t). The latter is the number of aerosol particles per unit length of the corresponding generation of TBT. Нere and later, t is time, and z is the coordinate along the TBT generations. It is counted from the larynx. Here, the advantage of considering generations instead of TBT branches is used explicitly. Generation is the entire set of TBT branches with the same division number. If you follow the flow of particles along the branches of the TBT, it will have gaps (vary approximately twice) at the points of division. If the flow is attributed to generation, it will be continuous. The continuity equation will be

> ∂n ∂t

� <sup>∂</sup>Un

Here, q (z, t) denotes the density of aerosol deposition on the inner surface of the TBT (per unit generation length and per unit time). Since precipitation occurs faster (in a few breath-exhalation cycles) than cleaning the lungs, sedimentation

Here, the function H (t) describes the temporal dependence of the entry of aerosol particles into the lungs, and Q (z) is the distribution of the trapped particles along the TBT generations. The latter function is directly related to the probability of capture Pi found in Section 3. To determine the relationship, let us establish the correspondence of the z coordinate to the generation number i. Choosing the origin (z = 0) at the beginning of the trachea, we get that the end of the ith generation

z ¼ zi � ∑

zi <sup>¼</sup> l0 <sup>∗</sup> <sup>1</sup> � <sup>2</sup>�ð Þ <sup>i</sup>þ<sup>1</sup> <sup>=</sup><sup>3</sup>

rate in the lungs W and the concentration of particles R in inhaled air:

k¼i k¼0

<sup>1</sup> � <sup>2</sup>�1=<sup>3</sup> <sup>¼</sup> <sup>Z</sup> <sup>∗</sup> <sup>1</sup> � <sup>2</sup>�iþ<sup>1</sup>

Here, lk is the length of the kth generation. In the model of E. Weibel, it changes with the number k as well as the diameter: lk = l0\*2�k/3. The geometric progression

Here, Z = l0/ (1–2�1/3) ≈ 29.1 cm is the total length of the TBT generations. The distribution of trapped particles Q (z) can be determined through the air exchange

We can express the mucociliary velocity U in terms of the z coordinate along the

<sup>∂</sup><sup>z</sup> <sup>¼</sup> q xð Þ ;<sup>t</sup> (20)

q zð Þ¼ ;t Q zð Þ ∗ H tð Þ (21)

lk (22)

(23)

3

Q zð Þ¼ W∗ R ∗ Pi=li (24)

U zð Þ¼ U0 <sup>∗</sup> ð Þ <sup>1</sup> � <sup>z</sup>=<sup>Z</sup> <sup>2</sup> (25)

morphometric model of E. Weibel.

DOI: http://dx.doi.org/10.5772/intechopen.84254

4.1 Formulation of the problem

density can be factorized:

corresponds to the coordinate

(Eq. (22)) is summed. For zi we have

TBT generations:

109

The probability of reaching the deep parts of the lungs is almost zero for particles with sizes larger than 10 microns. Particles with sizes of 4–6 microns can penetrate into the lungs only up to the terminal bronchioles. Only particles with a size of less than 2 microns can penetrate into the respiratory branches of the TBT, but the cleaning coefficient for them is close to 1. These data also indicate a high efficiency of the capture of particles as a mechanism for cleaning air from aerosols.

Finally, we present data on the probability of trapping particles in the airway of the TBT. Denoting this value as Рс (i), we can calculate it using the obvious formula

$$\mathbf{P\_{c}(i)} = \mathbf{P\_{a}(i)} \ast \mathbf{P(i)}\tag{19}$$

The results of the calculation of Pc (i) are shown in Figure 3.

The latter result gives a clear picture of the areas of the lungs in which the capture of aerosol particles of different diameters occurs. It is seen that even the smallest particles are mainly trapped in the bronchioles to the respiratory part. Large particles (with a diameter of 6 microns and more) do not penetrate into the terminal bronchioles at all.

#### 4. Model of mucociliary clearance of the lungs

As a rule, the modification of the model of E. Weibel leads to more realistic conclusions about the work of the lungs. In addition, it allows us to describe a number of effects that are impossible in the "classical" TBT model. There are, however, negative aspects of such a modification. Asymmetric statistical models are analyzed only by numerical methods. At the same time, however, the visibility of the results and the possibility of generalization the conclusions are lost. Therefore,

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance DOI: http://dx.doi.org/10.5772/intechopen.84254

for the initial estimates, it is advisable to limit ourselves to the simpler initial morphometric model of E. Weibel.

#### 4.1 Formulation of the problem

trachea, the zonal extrapulmonary bronchi, and the intrapulmonary subsegmental

The distribution of trapped particles of different diameters d by the numbers of generation of TBT.

РаðÞ¼ i

The results of the calculation of Pc (i) are shown in Figure 3.

4. Model of mucociliary clearance of the lungs

The latter result gives a clear picture of the areas of the lungs in which the capture of aerosol particles of different diameters occurs. It is seen that even the smallest particles are mainly trapped in the bronchioles to the respiratory part. Large particles (with a diameter of 6 microns and more) do not penetrate into the

As a rule, the modification of the model of E. Weibel leads to more realistic conclusions about the work of the lungs. In addition, it allows us to describe a number of effects that are impossible in the "classical" TBT model. There are, however, negative aspects of such a modification. Asymmetric statistical models are analyzed only by numerical methods. At the same time, however, the visibility of the results and the possibility of generalization the conclusions are lost. Therefore,

k <sup>Y</sup><sup>¼</sup>i�<sup>1</sup> k¼0

cles with sizes larger than 10 microns. Particles with sizes of 4–6 microns can penetrate into the lungs only up to the terminal bronchioles. Only particles with a size of less than 2 microns can penetrate into the respiratory branches of the TBT, but the cleaning coefficient for them is close to 1. These data also indicate a high efficiency of the capture of particles as a mechanism for cleaning air from aerosols. Finally, we present data on the probability of trapping particles in the airway of the TBT. Denoting this value as Рс (i), we can calculate it using the obvious

The probability of capture (formula Eq.(8)) determines the proportion of captured particles from the total number of particles entering the channel. It is possible to determine the probability Pa (i) of the particles reaching the ith generation:

The probability of reaching the deep parts of the lungs is almost zero for parti-

½ � 1 � P kð Þ (18)

PcðÞ¼ i Pað Þi ∗ P ið Þ (19)

bronchi.

Figure 3.

Rhinosinusitis

formula

108

terminal bronchioles at all.

We introduce adequate variables to describe the dynamics of mucociliary cleaning of the lungs. We will describe the flow j (z, t) of the mucociliary escalation by the product of the velocity U(z) and the density n (z, t). The latter is the number of aerosol particles per unit length of the corresponding generation of TBT. Нere and later, t is time, and z is the coordinate along the TBT generations. It is counted from the larynx. Here, the advantage of considering generations instead of TBT branches is used explicitly. Generation is the entire set of TBT branches with the same division number. If you follow the flow of particles along the branches of the TBT, it will have gaps (vary approximately twice) at the points of division. If the flow is attributed to generation, it will be continuous. The continuity equation will be

$$\frac{\partial n}{\partial t} - \frac{\partial Un}{\partial \mathbf{z}} = \mathbf{q}(\mathbf{x}, \mathbf{t}) \tag{20}$$

Here, q (z, t) denotes the density of aerosol deposition on the inner surface of the TBT (per unit generation length and per unit time). Since precipitation occurs faster (in a few breath-exhalation cycles) than cleaning the lungs, sedimentation density can be factorized:

$$\mathbf{q(z,t) = Q(z) \* H(t)}\tag{21}$$

Here, the function H (t) describes the temporal dependence of the entry of aerosol particles into the lungs, and Q (z) is the distribution of the trapped particles along the TBT generations. The latter function is directly related to the probability of capture Pi found in Section 3. To determine the relationship, let us establish the correspondence of the z coordinate to the generation number i. Choosing the origin (z = 0) at the beginning of the trachea, we get that the end of the ith generation corresponds to the coordinate

$$\mathbf{z} = \mathbf{z}\_{\mathbf{i}} \equiv \sum\_{\mathbf{k}=\mathbf{0}}^{\mathbf{k}=\mathbf{i}} \mathbf{l}\_{\mathbf{k}} \tag{22}$$

Here, lk is the length of the kth generation. In the model of E. Weibel, it changes with the number k as well as the diameter: lk = l0\*2�k/3. The geometric progression (Eq. (22)) is summed. For zi we have

$$\mathbf{z}\_{i} = \mathbf{l}\_{0} \* \frac{\mathbf{1} - \mathbf{2}^{-(i+1)/3}}{\mathbf{1} - \mathbf{2}^{-1/3}} = \mathbf{Z} \* \left(\mathbf{1} - \mathbf{2}^{-\frac{i+1}{3}}\right) \tag{23}$$

Here, Z = l0/ (1–2�1/3) ≈ 29.1 cm is the total length of the TBT generations. The distribution of trapped particles Q (z) can be determined through the air exchange rate in the lungs W and the concentration of particles R in inhaled air:

$$\mathbf{Q}(\mathbf{z}) = \mathbf{W} \* \mathbf{R} \* \mathbf{P}\_{\mathbf{i}}/\mathbf{l}\_{\mathbf{i}} \tag{24}$$

We can express the mucociliary velocity U in terms of the z coordinate along the TBT generations:

$$\mathbf{U}(\mathbf{z}) = \mathbf{U}\_0 \* (\mathbf{1} - \mathbf{z}/\mathbf{Z})^2 \tag{25}$$

#### Rhinosinusitis

where U0 ≈ 2.4 \* 10�<sup>2</sup> cm/s is the velocity at the entrance to the TBT (near the larynx). With such values of the scale of length and velocity, the time scale is T = Z/U0 ≈ 1.2\*10<sup>3</sup> sec = 0.337 hour. The relation (Eq. (25)) is valid only at the bifurcation points, but in the future, it is also advisable to use it at intermediate values of z.

Instead of the density n (z, t), we introduce a new variable—the stream j (z, t) = n (z, t) \* U (z), and instead of the z coordinate, a new variable

$$\mathfrak{G}(\mathbf{x}) = \int\_0^\mathbf{x} \mathbf{dx} / \mathbf{U}(\mathbf{x}) \tag{26}$$

The meaning of this variable is the time required for the escalation of particles from depth z to the TBT. For the dependence U (z) given by formula (Eq. (25)), we obtain

$$\frac{8}{\text{T}} = \frac{\text{z}/\text{Z}}{\text{1} - \text{z}/\text{Z}} ; \frac{\text{z}}{\text{Z}} = \frac{8/\text{T}}{\text{1} + 8/\text{T}} \tag{27}$$

From Eq. (20) we obtain the equation for the flow:

$$\frac{\partial \dot{\jmath}}{\partial t} - \frac{\partial \dot{\jmath}}{\partial \theta} = \mathbf{q}(\theta, \mathbf{t}) \* \mathbf{U}(\theta) \tag{28}$$

Integration over time of the effluent (30) gives us

Pi=l ð Þ<sup>i</sup> dz ¼ . dz ¼ l ð Þ¼ idi . T0WR∑Pi ¼ T0W (31)

H tðÞ¼ ηð Þt ∗ ηð Þ t0 � t (32)

9 >=

>;

Q z<sup>0</sup> ð Þdz<sup>0</sup> (34)

(33)

0 if t , 0

<sup>ϑ</sup> Q ϑ<sup>0</sup> ð Þv ϑ<sup>0</sup> ð Þdϑ<sup>0</sup> if t0 . t . 0

<sup>ϑ</sup>þt�t0 <sup>Q</sup> <sup>ϑ</sup><sup>0</sup> ð Þ<sup>v</sup> <sup>ϑ</sup><sup>0</sup> ð Þdϑ<sup>0</sup> if t . t0

Interesting is the observed value of the flow of precipitated particles at the exit from the TBT, with ϑ = 0. If we return from the variable ϑ to the coordinate z, we

> ðz2 z1

z2 ðÞ¼ t U0 ∗ t=½ � 1 þ U0 ∗ t=Z ; z1 ðÞ¼ t 0 for t , t0 and z1 ðÞ¼ t z2 ð Þ t � t0 for t . t0:

The last equality indicates the preservation of the number of aerosol particles: the total number of particles released (left side of Eq. (31)) is equal to the number of

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance

We get other results in the case of prolonged (many hours) inhalation of dusty air. This situation is realized in the case of industrial dust exposure, which continues during the work shift. We assume that such an effect begins at t = 0 and stops at t = t0 = 8 hours. In this case, the function H (t) in the relation (Eq. (21)) can be

Here, η (t) is the unit Dirac function, which is zero for negative values of the

In this case, the solution (Eq. (29)) of Eq. (28) can be written as

j 0ð Þ¼ ;t . 0

Ð <sup>ϑ</sup>þ<sup>t</sup>

Ð <sup>ϑ</sup>þ<sup>t</sup>

8 ><

>:

0

Dependence on time t (hour) of the particle flux at the exit of TBT.

DOI: http://dx.doi.org/10.5772/intechopen.84254

inhalation particles (right side of Eq. (31)).

jð Þ¼ ϑ;t

ð<sup>∞</sup> 0

Figure 4.

get

111

where

j oð Þ ;<sup>t</sup> dt <sup>¼</sup> T0WR <sup>ð</sup><sup>Z</sup>

given by the formula

argument and one for positive.

Its solution for the natural boundary conditions (j ! 0 as ϑ ! ∞ or t ! - ∞) has the form:

$$\mathbf{j}(\theta, \mathbf{t}) = \int\_0^\infty \mathbf{q}(\theta + \mathbf{u}, \mathbf{t} - \mathbf{u}) \* \mathbf{U}(\theta + \mathbf{u}) \* \mathbf{du} \tag{29}$$

Solution (Eq. (29)) of the problem of the dynamics of mucociliary clearance of the lungs allows us to estimate the parameters of escalation in several practically important cases.

#### 4.2. Examples of solutions

Let us consider, for example, the dynamics of excretion of particles inhaled in a short time T0. In this case, the time dependence of the settling of particles can be chosen in the form H (t) = T0\*δ (t), where δ (t) is the Dirac impulse delta function.

The presence of the delta function allows us to determine the integral (29) and record the outflow of TBT (with z = 0) of the flow of settled particles in the form

$$\mathbf{j(z=0,t)} = \mathbf{T\_0}\mathbf{W}\mathbf{R}[\mathbf{U(z)} \* \mathbf{P\_i/l\_i}]\_{z=z(t)}\tag{30}$$

The content of square brackets in this ratio is taken for z, which varies with time according to Eq. (27). The results of calculations of the time dependence of the removal of particles of different diameters are shown in Figure 4.

As was to be expected, particles of large diameters (10–20 μm), which are captured in proximal sections of the TBT (generation with number i = 6–11), are removed in the first hours (5–9 hours). This can be explained by the high rate of mucociliary escalation (U = 10�<sup>3</sup> –10�<sup>4</sup> cm/s). The maximum is clearly visible in the time dependence of the flow. Small particles (d = 1–4 μm) are removed much more slowly (1 day or more), and the time dependence of the flux decreases monotonously.

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance DOI: http://dx.doi.org/10.5772/intechopen.84254

#### Figure 4.

where U0 ≈ 2.4 \* 10�<sup>2</sup> cm/s is the velocity at the entrance to the TBT (near the

larynx). With such values of the scale of length and velocity, the time scale is T = Z/U0 ≈ 1.2\*10<sup>3</sup> sec = 0.337 hour. The relation (Eq. (25)) is valid only at the bifurcation points, but in the future, it is also advisable to use it at intermediate

Instead of the density n (z, t), we introduce a new variable—the stream j (z, t) = n (z, t) \* U (z), and instead of the z coordinate, a new variable

> ðx 0

The meaning of this variable is the time required for the escalation of particles from depth z to the TBT. For the dependence U (z) given by formula (Eq. (25)), we

z

Its solution for the natural boundary conditions (j ! 0 as ϑ ! ∞ or t ! - ∞) has

Solution (Eq. (29)) of the problem of the dynamics of mucociliary clearance of the lungs allows us to estimate the parameters of escalation in several practically

Let us consider, for example, the dynamics of excretion of particles inhaled in a short time T0. In this case, the time dependence of the settling of particles can be chosen in the form H (t) = T0\*δ (t), where δ (t) is the Dirac impulse delta function. The presence of the delta function allows us to determine the integral (29) and record the outflow of TBT (with z = 0) of the flow of settled particles in the form

The content of square brackets in this ratio is taken for z, which varies with time

according to Eq. (27). The results of calculations of the time dependence of the

As was to be expected, particles of large diameters (10–20 μm), which are captured in proximal sections of the TBT (generation with number i = 6–11), are removed in the first hours (5–9 hours). This can be explained by the high rate of

the time dependence of the flow. Small particles (d = 1–4 μm) are removed much more slowly (1 day or more), and the time dependence of the flux decreases

removal of particles of different diameters are shown in Figure 4.

j zð Þ¼ <sup>¼</sup> <sup>0</sup>;<sup>t</sup> T0WR U zð Þ <sup>∗</sup> Pi=li ½ �<sup>z</sup> <sup>¼</sup> z tð Þ (30)

–10�<sup>4</sup> cm/s). The maximum is clearly visible in

<sup>Z</sup> <sup>¼</sup> <sup>ϑ</sup>=<sup>T</sup>

dx=U xð Þ (26)

<sup>1</sup> <sup>þ</sup> <sup>ϑ</sup>=<sup>T</sup> (27)

<sup>∂</sup><sup>ϑ</sup> <sup>¼</sup> <sup>q</sup>ð Þ <sup>ϑ</sup>;<sup>t</sup> <sup>∗</sup> <sup>U</sup>ð Þ <sup>ϑ</sup> (28)

qð Þ ϑ þ u;t � u ∗ Uð Þ ϑ þ u ∗ du (29)

ϑð Þ¼ x

ϑ <sup>T</sup> <sup>¼</sup> <sup>z</sup>=<sup>Z</sup> <sup>1</sup> � <sup>z</sup>=<sup>Z</sup> ;

From Eq. (20) we obtain the equation for the flow:

jð Þ¼ ϑ;t

∂j ∂t � ∂j

ð<sup>∞</sup> 0

values of z.

Rhinosinusitis

obtain

the form:

important cases.

4.2. Examples of solutions

mucociliary escalation (U = 10�<sup>3</sup>

monotonously.

110

Dependence on time t (hour) of the particle flux at the exit of TBT.

Integration over time of the effluent (30) gives us

$$\int\_{0}^{\infty} \mathbf{j(o,t)dt} = \mathbf{T\_0} \mathbf{W} \mathbf{R} \int\_{0}^{Z} (\mathbf{P\_i/l\_i}) d\mathbf{z} = \mathbf{\hat{z}} \,(d\mathbf{z} = \mathbf{l\_i} \mathbf{di}) = \mathbf{\hat{z}} \,\mathbf{W} \mathbf{R} \sum \mathbf{P\_i} = \mathbf{T\_0} \mathbf{W} \tag{31}$$

The last equality indicates the preservation of the number of aerosol particles: the total number of particles released (left side of Eq. (31)) is equal to the number of inhalation particles (right side of Eq. (31)).

We get other results in the case of prolonged (many hours) inhalation of dusty air. This situation is realized in the case of industrial dust exposure, which continues during the work shift. We assume that such an effect begins at t = 0 and stops at t = t0 = 8 hours. In this case, the function H (t) in the relation (Eq. (21)) can be given by the formula

$$\mathbf{H}(\mathbf{t}) = \boldsymbol{\eta}(\mathbf{t}) \* \boldsymbol{\eta}(\mathbf{t}\_0 - \mathbf{t}) \tag{32}$$

Here, η (t) is the unit Dirac function, which is zero for negative values of the argument and one for positive.

In this case, the solution (Eq. (29)) of Eq. (28) can be written as

$$\mathbf{j}(\boldsymbol{\theta}, \mathbf{t}) = \left\{ \begin{array}{c} \mathbf{0} \\ \int\_{\boldsymbol{\theta}}^{\boldsymbol{\theta} + \mathbf{t}} \mathbf{Q}(\boldsymbol{\theta}') \mathbf{v}(\boldsymbol{\theta}') d\boldsymbol{\theta}' \\ \int\_{\boldsymbol{\theta} + \mathbf{t} - \mathbf{t}\_0}^{\boldsymbol{\theta} + \mathbf{t}} \mathbf{Q}(\boldsymbol{\theta}') \mathbf{v}(\boldsymbol{\theta}') d\boldsymbol{\theta}' \end{array} \text{ if } \mathbf{t} > \mathbf{t} > \mathbf{0} \right\} \tag{33}$$

Interesting is the observed value of the flow of precipitated particles at the exit from the TBT, with ϑ = 0. If we return from the variable ϑ to the coordinate z, we get

$$\mathbf{j}(\mathbf{0}, \mathbf{t} \ge \mathbf{0}) = \int\_{\mathbf{z}\_1}^{\mathbf{z}\_2} \mathbf{Q}(\mathbf{z}') \mathbf{d} \mathbf{z}' \tag{34}$$

where

$$\mathbf{z}\_2\ (\mathbf{t}) = \mathbf{U}\_0 \ast \mathbf{t} / [\mathbf{1} + \mathbf{U}\_0 \ast \mathbf{t} / \mathbf{Z}] \\ \vdots \\ \mathbf{z}\_1\ (\mathbf{t}) = \mathbf{0} \text{ for } \mathbf{t} \le \mathbf{t}\_0 \text{ and } \mathbf{z}\_1\ (\mathbf{t}) = \mathbf{z}\_2 \ (\mathbf{t} - \mathbf{t}\_0) \text{ for } \mathbf{t} \ge \mathbf{t}\_0$$

Formula (34) clearly demonstrates where the deposited particles are drawn from at a certain point in time.

Dependences on time t are depicted graphically in Figure 5. Along the abscissa, the time (in hours) after the start of inhalation of aerosol particles (10 μm in diameter) is plotted. Inhalation is assumed to occur over a period of time from t = 0 to t0 = 8 hours. The boundaries z1 and z2 of the integration range are represented by the corresponding curves. The sedimentation density distribution of particles along the z coordinate (the integrand function Q (z) in the formula (Eq. (34)) is shown in the same graph. It can be seen that in the first hours after the start of inhalation, particles are removed from the TBT sections that are at a distance from the beginning to ≈ 0.8 \* Z. A small part of the aerosol particles is captured here. Over time, the length of the cleaning section increases. After 5 hours, the removal of particles from a depth of ≈ 0.94 \* Z begins. From here, a significant part of the trapped aerosol particles is removed. The maximum elimination occurs at t ≈ 10 hours after the start. After this, the boundaries of the dust removal area are shifted further into TBT, where relatively little dust accumulates. By 15–20 hours after the start of inhalation (7–12 hours after the end), almost all captured particles are removed from TBT.

Thus, breaks lasting 7-12 hours eliminate the risk of accumulation of aerosol particles with a diameter of 10 mkm in the lungs. It is assumed, of course, that the mass of the trapped aerosol particles is not too large, so that its removal does not exceed the mass transfer capabilities during mucociliary escalation.

complete removal of such particles, it is necessary not less than 2 days. More time is

The fact of accumulation of small particles during periodic inhalation of dusty air (during working hours), alternating with periods of rest when mucociliary clearance of the lungs occurs, raises the question of the average number of particles

Eq. (20), which describes the dynamics of mucociliary cleaning of the lungs, we can determine the density n (z, t) of the captured particles n(z, t) = j (z, t) / U (z). And after integrating by TBT generation, we have the total amount of N (t) particles in

> ð<sup>∞</sup> 0 dϑ ð<sup>∞</sup> 0

dϑUð Þ ϑ Qð Þ ϑ

Now, the time dependence of the capture of particles H (t) is a sequence of rectangular pulses of the form (Eq. (32)), duration t0 (working time), and periodicity t1 (1 day). A rigorous computation of the double integral (Eq. (36)) is possible, but not interesting. It is clear that this is a certain function of time, which varies with the period of the day around a certain mean value <N>. It can be estimated quite accurately. Indeed, the calculation of the average < N > in time from N (t) reduces to replacing in the integral H (t-τ) by <H (t-τ)>. The last value does not

dϑvð Þ ϑ Qð Þ ϑ ϑ ¼ , H . ∗

As in formulas (Eq. (26–27)), the meaning of the variable ϑ (z) is the time required for the escalation of particles from depth z to the beginning of the TBT.

ðX 0

Obvious substitutions of order and integration variables simplify this expression

ðϑ 0

dτqð Þ ϑ þ τ;t � τ vð Þ ϑ þ τ (35)

dτHð Þ t � τ (36)

dxQ xð Þϑð Þ x (37)

characterizing the result of such accumulation. Using solution (Eq. (33)) of

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance

needed to remove smaller particles.

DOI: http://dx.doi.org/10.5772/intechopen.84254

the lungs:

Figure 6.

to the form:

113

N tðÞ¼ <sup>ð</sup><sup>Z</sup>

0

dz j zð Þ ;<sup>t</sup> U zð Þ <sup>¼</sup> ð<sup>∞</sup> 0

N tðÞ¼ <sup>ð</sup><sup>∞</sup>

depend on τ; therefore, from Eq. (36) we can get

ð<sup>∞</sup> 0

, N . ¼ , H . ∗

0

dϑjð Þ¼ ϑ;t

The time dependence of the elimination of aerosol particles trapped in 8 hours of work.

Cleaning the lungs from smaller dust particles occurs differently. The most probable capture depth is ≈ 0.98 \* Z for particles with g = 4 μm. In these generations of TBT, the rate of mucociliary escalation is 10 times less than in the trachea. Cleaning time is much longer.

The results of direct calculations of the flow at the TBT outlet confirm these conclusions (see Figure 6).

It is seen that particles with d = 20 μm practically do not stay in the lungs; they come out of TBT as they arrive, of course, if the mass of trapped particles does not exceed the mass transfer capacity during mucociliary escalation. Particles with d < 4 μm are deposited in the deep sections of TBT. They do not have time to get out of the lungs during the break. Their share exceeds 30%. This indicates the possibility of accumulation in the lungs of trapped particles of small diameter. For almost

Figure 5.

The dependence on time t (hour) of the boundaries of removal from the TBT trapped particles with a diameter of 10 mkm.

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance DOI: http://dx.doi.org/10.5772/intechopen.84254

#### Figure 6.

Formula (34) clearly demonstrates where the deposited particles are drawn from

Dependences on time t are depicted graphically in Figure 5. Along the abscissa,

Thus, breaks lasting 7-12 hours eliminate the risk of accumulation of aerosol particles with a diameter of 10 mkm in the lungs. It is assumed, of course, that the mass of the trapped aerosol particles is not too large, so that its removal does not

Cleaning the lungs from smaller dust particles occurs differently. The most probable capture depth is ≈ 0.98 \* Z for particles with g = 4 μm. In these generations of TBT, the rate of mucociliary escalation is 10 times less than in the trachea.

The results of direct calculations of the flow at the TBT outlet confirm these

It is seen that particles with d = 20 μm practically do not stay in the lungs; they come out of TBT as they arrive, of course, if the mass of trapped particles does not exceed the mass transfer capacity during mucociliary escalation. Particles with d < 4 μm are deposited in the deep sections of TBT. They do not have time to get out of the lungs during the break. Their share exceeds 30%. This indicates the possibility of accumulation in the lungs of trapped particles of small diameter. For almost

The dependence on time t (hour) of the boundaries of removal from the TBT trapped particles with a diameter

exceed the mass transfer capabilities during mucociliary escalation.

the time (in hours) after the start of inhalation of aerosol particles (10 μm in diameter) is plotted. Inhalation is assumed to occur over a period of time from t = 0 to t0 = 8 hours. The boundaries z1 and z2 of the integration range are represented by the corresponding curves. The sedimentation density distribution of particles along the z coordinate (the integrand function Q (z) in the formula (Eq. (34)) is shown in the same graph. It can be seen that in the first hours after the start of inhalation, particles are removed from the TBT sections that are at a distance from the beginning to ≈ 0.8 \* Z. A small part of the aerosol particles is captured here. Over time, the length of the cleaning section increases. After 5 hours, the removal of particles from a depth of ≈ 0.94 \* Z begins. From here, a significant part of the trapped aerosol particles is removed. The maximum elimination occurs at t ≈ 10 hours after the start. After this, the boundaries of the dust removal area are shifted further into TBT, where relatively little dust accumulates. By 15–20 hours after the start of inhalation (7–12 hours after the end), almost all captured particles are removed

at a certain point in time.

Rhinosinusitis

from TBT.

Figure 5.

of 10 mkm.

112

Cleaning time is much longer.

conclusions (see Figure 6).

The time dependence of the elimination of aerosol particles trapped in 8 hours of work.

complete removal of such particles, it is necessary not less than 2 days. More time is needed to remove smaller particles.

The fact of accumulation of small particles during periodic inhalation of dusty air (during working hours), alternating with periods of rest when mucociliary clearance of the lungs occurs, raises the question of the average number of particles characterizing the result of such accumulation. Using solution (Eq. (33)) of Eq. (20), which describes the dynamics of mucociliary cleaning of the lungs, we can determine the density n (z, t) of the captured particles n(z, t) = j (z, t) / U (z). And after integrating by TBT generation, we have the total amount of N (t) particles in the lungs:

$$\mathbf{N(t)} = \int\_{0}^{Z} d\mathbf{z} \frac{\mathbf{j(z,t)}}{\mathbf{U(z)}} = \int\_{0}^{\infty} d\theta \mathbf{j(\theta, t)} = \int\_{0}^{\infty} d\theta \int\_{0}^{\infty} d\tau q(\theta + \tau, \mathbf{t} - \tau) \mathbf{v}(\theta + \tau) \tag{35}$$

Obvious substitutions of order and integration variables simplify this expression to the form:

$$\mathbf{N(t)} = \int\_0^\infty d\theta U(\theta) \mathbf{Q(\theta)} \int\_0^\theta d\tau H(\mathbf{t} - \tau) \tag{36}$$

Now, the time dependence of the capture of particles H (t) is a sequence of rectangular pulses of the form (Eq. (32)), duration t0 (working time), and periodicity t1 (1 day). A rigorous computation of the double integral (Eq. (36)) is possible, but not interesting. It is clear that this is a certain function of time, which varies with the period of the day around a certain mean value <N>. It can be estimated quite accurately. Indeed, the calculation of the average < N > in time from N (t) reduces to replacing in the integral H (t-τ) by <H (t-τ)>. The last value does not depend on τ; therefore, from Eq. (36) we can get

$$\mathbf{x} \le \mathbf{N} \ge \mathbf{w} \le \mathbf{w} \int\_0^\infty d\theta v(\theta) \mathbf{Q}(\theta) \theta = \le H \ge \ast \int\_0^\mathbf{X} \mathbf{dx} \mathbf{Q}(\mathbf{x}) \theta(\mathbf{x}) \tag{37}$$

As in formulas (Eq. (26–27)), the meaning of the variable ϑ (z) is the time required for the escalation of particles from depth z to the beginning of the TBT.


Table 3.

The average exit time for particles of different diameters.

The same technique that was used in the transition in formula (Eq. (21)) from the integral to the sum over the generation numbers leads to the expression

$$ =  \* W\*R\* < \emptyset> \tag{38}$$

Here, the average exit time < ϑ > is determined by the probability distribution of the capture of particles of the corresponding diameter:

$$<\emptyset> = \sum \mathbf{P}\_{\mathbf{i}} \ast \theta\_{\mathbf{i}} \tag{39}$$

Author details

115

Gennady Fedorovitch

NTM-Defense, Moscow, Russia

provided the original work is properly cited.

\*Address all correspondence to: fedorgv@gmail.com

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance

DOI: http://dx.doi.org/10.5772/intechopen.84254

The relation (Eq. (38)) is quite clearly interpreted: the average number of aerosol particles in the lungs is equal to their number in the volume of air, which is determined by pulmonary air exchange over the average time of exiting <ϑ > .

The obtained values of the average time of exit are in complete agreement with the results on the time dependence of the output stream for various inhalation regimes given above. The average exit time for particles of different diameters is given in Table 3.

#### 5. Conclusions

We investigated the deposition, retention, and clearance of inhaled aerosol particles in the lungs. Each of these processes is intensively investigated in numerous papers. As a rule, for their analysis, complex computer simulation of air and hydrodynamic flows in the branches of TBT is used.

The phenomenological model we have constructed is a simplified picture of real processes. It is used to study the key aspects of aerosol behavior in the lungs. This is a self-consistent and closed set of assumptions about the process, reflecting the basic, though not all, properties of a real phenomenon. Simulation is a research method with several goals:


The model reflects the unity of the main functional structures and, at the same time, the specifics of the processes occurring in them. Despite the descriptive nature of the proposed modeling, its meaning is that it gives the plot accuracy and compositional clarity, which are not inherent in real objects of modeling.

#### Conflict of interest

The author does not have a conflict of interest.

Aerosol Particles in Lungs: Theoretical Modeling of Deposition and Mucociliary Clearance DOI: http://dx.doi.org/10.5772/intechopen.84254
