A two-component theory of aerosol deposition in lung airways

The deposition of aerosol particles in the human lung airways is due to two distinct mechanisms. One is by direct deposition resulting from diffusion, sedimentation and impaction as the aerosol moves in and out of the lung. The other is an indirect mechanism by which particles are transported mechanically from the tidal air to the residential air and eventually captured by the airways due to intrinsic particle motion. This last mechanism is not well understood at present. Using a trumpet airway model constructed from Weibel's data, a two-component theory is developed. In this theory, the particle concentrations in the airways and the alveoli at a given airway depth are considered to be quantitatively different. This difference in concentrations will cause a net mixing between the tidal and residential aerosol as the aerosol is breathed in and out. A distribution parameter is then introduced to account for the distribution of ventilation. The effect of intrinsic particle motion on the aerosol mixing is also included. From this theory, total and regional deposition in the lung at the steady mouth breathing without pause is calculated for several different respiratory cycles. The results agree reasonably well with the experimental data.     

Particle deposition in obstructed airways

One approach to tackle the particle deposition in human lungs in close proximity is to develop an understanding of the particle motion in bifurcation airways. Chronic obstructive pulmonary disease (COPD) is one of the most common diseases in humans. COPD always results in inflammation that leads to narrowing and obstructing of the airways. The obstructive airways can alter the respiratory flow and particle deposition significantly. In order to study the effect of obstruction on particle deposition, four three-dimensional four-generation lung models based on the 23-generation model of Weibel [1963. Morphometry of the Human Lung. New York Academic Press, Springer, Berlin.] have been generated. The fully three-dimensional incompressible laminar Navier-Stokes equations are solved using computational fluid dynamics (CFD) solver on structured hexahedral meshes. Subsequently, a symmetric four-generation airway model serves as the reference and the other three models are considered to be obstructed at each generation, respectively. The calculation results show that the obstructive airway has significant influence on the particle deposition down-stream of the obstruction. The skewed velocity profile in the bifurcation airway is modified by the throat; consequently, more particles impact on the divider which results in higher deposition efficiency.     

Theoretical evaluation of polytripeptide collagen models

The influence of inserting certain residues (X) into a polytripeptide sequence conformed into a poly-L -proline II helix is examined theoretically. It is found that for sequences such as -Gly-Pro-X- and -Gly-X-Pro-, the introduction of glycyl, L -alanyl or L -seryl residues in the X position destabilizes the helix so that it is no longer the most stable intramolecular form. On the other hand, L -prolyl and L -hydroxyprolyl residues cause the PP II helix to be most stable. Of the many stable intramolecular forms, the majority will not pack efficiently to form fiber or solid-state structures. The Rich-Crick and Ramachandran collagen model structures were examined in terms of a Gly-Pro-Ala sequence, the Ramachandran, one-hydrogen-bond structure, being the most stable. However, another triple-strand structure for (Gly-Pro-Ala)_n, is much more energetically favorable. Hence, it may be concluded that none of the aforementioned is an entirely satisfactory collagen model. The new triple helix conformation proposed by Traub, Yonath, and Segal for (Gly-Pro-Pro) is found to give a more favorable intramolecular conformation for (Gly-Pro-Ala)_n than those derived from other collagen models. It is concluded that the collagen molecule derives its stability from interchain interactions in proline-sparse regions and intrachain stability in proline-rich regions.     

Models of human lung airways and their application to inhaled particle deposition

Models of the human respiratory tract were developed based on detailed morphometric measurements of a silicone rubber cast of the human tracheobronchial airways. Emphasis was placed on the “Typical Path Lung Model” which used one typical pathway to represent a portion of the lung, such as a lobe, or to represent the whole lung. The models contain geometrical parameters, including airway segment diameters, lengths, branching angles and angles of inclination to gravity, which are needed for estimating inhaled particle deposition. Aerosol depositions for various breathing patterns and particle sizes were calculated using these lung models and the modified Findeisen-Landahl computational scheme. The results agree reasonably well with recent experimental data. Regional deposition, including lobar deposition fractions, are also calculated and compared with results based on the ICRP lung deposition model.     

Aerosol transport and deposition in sequentially bifurcating airways

Deposition patterns and efficiencies of a dilute suspension of inhaled particles in three-dimensional double bifurcating airway models for both in-plane and 90 deg out-of-plane configurations have been numerically simulated assuming steady, laminar, constant-property air flow with symmetry about the first bifurcation. Particle diameters of 3, 5, and 7 microns were used in the simulation, while the inlet Stokes and Reynolds numbers varied from 0.037 to 0.23 and 500 to 2000, respectively. Comparisons between these results and experimental data based on the same geometric configuration showed good agreement. The overall trend of the particle deposition efficiency, i.e., an exponential increase with Stokes number, was somewhat similar for all bifurcations. However, the deposition efficiency of the first bifurcation was always larger than that of the second bifurcation, while in general the particle efficiency of the out-of-plane configuration was larger than that of the in-plane configuration. The local deposition patterns consistently showed that the majority of the deposition occurred in the carinal region. The distribution pattern in the first bifurcation for both configurations were symmetric about the carina, which was a direct result of the uniaxial flow at the inlet. The deposition patterns about the second carina showed increased asymmetry due to highly nonuniform flow generated by the first bifurcation and were extremely sensitive to bifurcation orientation. Based on the deposition variations between bifurcation levels and orientations, the use of single bifurcation models was determined to be inadequate to resolve the complex fluid-particle interactions that occur in multigenerational airways.     

Targeted drug aerosol deposition analysis for a four-generation lung airway model with hemispherical tumors

One important research area of broad interest is the development of highly efficient drug delivery systems for desired site deposition and uptake. For example, controlled drug aerosol release and targeting to specific regions of the lung is a novel way to combat lung diseases, diabetes, virus infections, cancers, etc. Determination of feasible air-particle streams is a prerequisite for the development of such delivery devices, say, smart inhalers. The concept of "controlled particle release and targeting" is introduced and results are discussed for a representative model of bronchial lung airways afflicted with hemispherical tumors of different sizes and locations. It is shown that under normal particle inlet conditions a particle mass fraction of only up to 11% may deposit on the surface of a specific tumor with critical radius r/R approximately 1.25, while a controlled particle release achieves deposition fractions of 35 to 92% for a realistic combination of inlet Stokes and Reynolds numbers, depending mainly on tumor size. Furthermore, with the controlled release and targeting approach nearby healthy tissue is hardly impacted by the typically aggressive drug aerosols. Assuming laminar, quasi-steady, three-dimensional air flow and spherical non-interacting micron-particles in sequentially bifurcating rigid airways, the results were obtained using a validated commercial finite-volume code with user-enhanced programs on a high-end engineering workstation. The new concept is generic and hence should be applicable to other regions of the respiratory system as well.     

Computational fluid dynamics simulations of particle deposition in large-scale, multigenerational lung models

Computational fluid dynamics (CFD) has emerged as a useful tool for the prediction of airflow and particle transport within the human lung airway. Several published studies have demonstrated the use of Eulerian finite-volume CFD simulations coupled with Lagrangian particle tracking methods to determine local and regional particle deposition rates in small subsections of the bronchopulmonary tree. However, the simulation of particle transport and deposition in large-scale models encompassing more than a few generations is less common, due in part to the sheer size and complexity of the human lung airway. Highly resolved, fully coupled flowfield solution and particle tracking in the entire lung, for example, is currently an intractable problem and will remain so for the foreseeable future. This paper adopts a previously reported methodology for simulating large-scale regions of the lung airway (Walters, D. K., and Luke, W. H., 2010, "A Method for Three-Dimensional Navier-Stokes Simulations of Large-Scale Regions of the Human Lung Airway," ASME J. Fluids Eng., 132(5), p. 051101), which was shown to produce results similar to fully resolved geometries using approximate, reduced geometry models. The methodology is extended here to particle transport and deposition simulations. Lagrangian particle tracking simulations are performed in combination with Eulerian simulations of the airflow in an idealized representation of the human lung airway tree. Results using the reduced models are compared with those using the fully resolved models for an eight-generation region of the conducting zone. The agreement between fully resolved and reduced geometry simulations indicates that the new method can provide an accurate alternative for large-scale CFD simulations while potentially reducing the computational cost of these simulations by several orders of magnitude.     

Capillary-elastic instabilities of liquid-lined lung airways

To model the competition between capillary and elastic forces in controlling the shape of a small lung airway and its interior liquid lining, we compute the equilibrium configurations of a liquid-lined, externally pressurized, buckled elastic tube. We impose axial uniformity and assume that the liquid wets the tube wall with zero contact angle. Non-zero surface tension has a profound effect on the tube's quasi-steady inflation-deflation characteristics. At low liquid volumes, hysteresis arises through two distinct mechanisms, depending on the buckling wavenumber. Sufficient compression always leads to abrupt and irreversible collapse and flooding of the tube; flooding is promoted by increasing liquid volumes or surface tension. The model captures mechanisms whereby capillary-elastic instabilities can lead to airway closure.     

Respiratory deposition model of an inhaled aerosol bolus

The bolus delivery method is designed to deliver a dose to the desired location in the lung, and it has the advantage of fewer side effects and a more efficient way of delivery. Based upon the lung deposition model developed for continuously inhaling aerosols of constant concentration, a mathematical model of aerosol bolus deposition is proposed. The calculated results show that the recovery depends on the bolus penetration depth, flow rate, particle size, breath holding time and bolus volume. Three sets of published experimental data with different controlling factors (particle size, flow rate and breath holding time) are adopted to make the quantitative comparisons with the calculated results. The predictions and data for the low intrinsic motion particles (～1 μm) have good agreement, as do the coarse particles in the shallow airways region. For females, the recovery was found to be consistently lower than that for males.     

Respiratory deposition model of an inhaled aerosol bolus

The bolus delivery method is designed to deliver a dose to the desired location in the lung, and it has the advantage of fewer side effects and a more efficient way of delivery. Based upon the lung deposition model developed for continuously inhaling aerosols of constant concentration, a mathematical model of aerosol bolus deposition is proposed. The calculated results show that the recovery depends on the bolus penetration depth, flow rate, particle size, breath holding time and bolus volume. Three sets of published experimental data with different controlling factors (particle size, flow rate and breath holding time) are adopted to make the quantitative comparisons with the calculated results. The predictions and data for the low intrinsic motion particles (～1 μm) have good agreement, as do the coarse particles in the shallow airways region. For females, the recovery was found to be consistently lower than that for males.     

Theoretical study of nasal deposition

A theory is derived to calculate the regional and total deposition of aerosol particles in the nasal passages during inhalation. The particle size studied range from 0.2 to 10.0 μm diameter. The deposition is calculated in five regions; (I) the region filled with nasal hair, (II) the nasal valve, (III) the expansion region, (IV) the turbinate region and (V) the posterior bend. Equations are derived to determine the deposition caused by direct impaction on the nasal hairs and bends of the passages. The calculations show the deposition due to direct impaction does not account for the amount or location of deposited particles measured in experiments. Secondary flows have been speculated to exist in the expansion region after the nasal valve and an equation is derived to estimate the deposition caused by the secondary flows. The calculated deposition, due to direct impaction and secondary flows, shows general agreement with the experiment as to the predicted amount and location of deposited particles.     

Calculation of regional deposition of aerosol in the respiratory tract

The removal of air-borne particles in the respiratory tract is treated to enable regional deposition to be inferred from measurement of expired aerosol as well as predicted from theory of the primary removal processes. The analysis uses the analogy of a continuous tubular filter-bed and includes consideration of respiratory pauses and the mechanical mixing of gas flow. Derived equations relate regional deposition, distribution of aerosol in the expired air, and efficiency of removal at different depths in the respiratory tract.     

Effects of curved inlet tubes on air flow and particle deposition in bifurcating lung models

In vivo bifurcating airways are complex and the airway segments leading to the bifurcations are not always straight, but curved to various degrees. How do such curved inlet tubes influence the motion as well as local deposition and hence the biological responses of inhaled particulate matter in lung airways? In this paper steady laminar dilute suspension flows of micron-particles are simulated in realistic double bifurcations with curved inlet tubes, i.e., 0 degrees < or =theta< or =90 degrees, using a commercial finite-volume code with user-enhanced programs. The resulting air-flow patterns as well as particle transport and wall depositions were analyzed for different flow inlet conditions, i.e., uniform and parabolic velocity profiles, and geometric configurations. The curved inlet segments have quite pronounced effects on air-flow, particle motion and wall deposition in the downstream bifurcating airways. In contrast to straight double bifurcations, those with bent parent tubes also exhibit irregular variations in particle deposition efficiencies as a function of Stokes number and Reynolds number. There are fewer particles deposited at mildly curved inlet segments, but the particle deposition efficiencies at the downstream sequential bifurcations vary much when compared to those with straight inlets. Under certain flow conditions in sharply curved lung airways, relatively high, localized particle depositions may take place. The findings provide necessary information for toxicologic or therapeutic impact assessments and for global lung dosimetry models of inhaled particulate matter.     

Precipitation of submicron charged particles in human lung airways

Precipitation of charged particles in a tube by their own space charge is investigated theoretically, when the number density of the particles is large enough so that the potential is a smooth function given by Poisson’s equation, and when the number density is small so that only the image force is important. These two approaches have been applied to the data given by Weibel for the human lung, to determine the deposition probabibilities at different generations for submicron particles when the particle density is 1×10⁵ particles/cm³. The results indicate that the electrostatic dispersion can only lead to a small effect on the lung deposition, the predominant effect is due to the image force exerted on the particles.     

A comparative anatomical study of the mammalian uterotubal junction

Among eight species of mammals in this study (cattle, sheep, pig-tail and rhesus monkeys, rabbit, pig, rat, and dog) four basic patterns of anatomical structure at the uterotubal junction are described. The classification of types is based upon the presence or absence of an intramural portion of the oviduct and of isthmal folds or plicae projecting into the lumen of the uterine cornu. Histological variations are reported for three tissues: epithelial and connective of the mucosa and smooth muscle of the tunica muscularis. In the epithelium during the estrous cycle the differences recorded include: (a) absence of ciliated cells in the distal end of the oviduct in rat and dog; (b) variations in ciliated and nonciliated cells in (1) cell height, (2) location, shape and stainability of the nucleus, and (3) in amount and stainability of apical cytoplasm; (c) presence of lymphoblast-like cells which appear to migrate through the epithelium from the lamina propria. The connective tissue of the mucosa, as a circular layer and as cores for the mucosal folds, shows variations in thickness and in relative density of cells and fibers of the matrix. Emphasis is given to the presence of an inner longitudinal layer of smooth muscle in the tunica muscularis of the distal oviduct in six of the eight species.