Secondary velocity fields in the conducting airways of the human lung

An understanding of flow and dispersion in the human respiratory airways is necessary to assess the toxicological impact of inhaled particulate matter as well as to optimize the design of inhalable pharmaceutical aerosols and their delivery systems. Secondary flows affect dispersion in the lung by mixing solute in the lumen cross section. The goal of this study is to measure and interpret these secondary velocity fields using in vitro lung models. Particle image velocimetry experiments were conducted in a three-generational, anatomically accurate model of the conducting region of the lung. Inspiration and expiration flows were examined under steady and oscillatory flow conditions. Results illustrate secondary flow fields as a function of flow direction, Reynolds number, axial location with respect to the bifurcation junction, generation, branch, phase in the oscillatory cycle, and Womersley number. Critical Dean number for the formation of secondary vortices in the airways, as well as the strength and development length of secondary flow, is characterized. The normalized secondary velocity magnitude was similar on inspiration and expiration and did not vary appreciably with generation or branch. Oscillatory flow fields were not significantly different from corresponding steady flow fields up to a Womersley number of 1 and no instabilities related to shear were detected on flow reversal. These observations were qualitatively interpreted with respect to the simple streaming, augmented dispersion, and steady streaming convective dispersion mechanisms.     

Effects of gas properties and waveform asymmetry on gas transport in a branching tube network

Local gas transport coefficients, quantifying longitudinal dispersion through a symmetrical constant-diameter tube network, have been measured during oscillation with both symmetrical and nonsymmetrical waveforms. Experiments were carried out over a range of conditions that would prevail in the central to lower airways during high-frequency ventilation at moderate frequency (5 Hz) and tidal volume (15-80 ml). Gas transport coefficients resulting from oscillation of three different resident-trace gas pairs were measured using a new analytic technique. This technique allowed rapid determination of the transport coefficient distribution along the entire network. Results demonstrate a small but significant influence attributable to changes in gas properties that is similar to that found in a straight tube and indicate that augmented dispersion is an important mechanism of axial transport. Gas transport coefficients were found to be unaffected by changes in flow waveform symmetry, suggesting that previously reported improvements in gas exchange associated with decreasing inspiratory to expiratory time ratios are not due to a change in local conditions such as asymmetry in the velocity profile.     

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.     

Aerosol deposition characteristics in distal acinar airways under cyclic breathing conditions

Although the major mechanisms of aerosol deposition in the lung are known, detailed quantitative data in anatomically realistic models are still lacking, especially in the acinar airways. In this study, an algorithm was developed to build multigenerational three-dimensional models of alveolated airways with arbitrary bifurcation angles and spherical alveolar shape. Using computational fluid dynamics, the deposition of 1- and 3-μm aerosol particles was predicted in models of human alveolar sac and terminal acinar bifurcation under rhythmic wall motion for two breathing conditions (functional residual capacity = 3 liter, tidal volume = 0.5 and 0.9 liter, breathing period = 4 s). Particles entering the model during one inspiration period were tracked for multiple breathing cycles until all particles deposited or escaped from the model. Flow recirculation inside alveoli occurred only during transition between inspiration and expiration and accounted for no more than 1% of the whole cycle. Weak flow irreversibility and convective transport were observed in both models. The average deposition efficiency was similar for both breathing conditions and for both models. Under normal gravity, total deposition was ~33 and 75%, of which ~67 and 96% occurred during the first cycle, for 1- and 3-μm particles, respectively. Under zero gravity, total deposition was ~2-5% for both particle sizes. These results support previous findings that gravitational sedimentation is the dominant deposition mechanism for micrometer-sized aerosols in acinar airways. The results also showed that moving walls and multiple breathing cycles are needed for accurate estimation of aerosol deposition in acinar airways.     

Convective flow dominates aerosol delivery to the lung segments

Most previous computational studies on aerosol transport in models of the central airways of the human lung have focused on deposition, rather than transport of particles through these airways to the subtended lung regions. Using a model of the bronchial tree extending from the trachea to the segmental bronchi (J Appl Physiol 98: 970-980, 2005), we predicted aerosol delivery to the lung segments. Transport of 0.5- to 10-μm-diameter particles was computed at various gravity levels (0-1.6 G) during steady inspiration (100-500 ml/s). For each condition, the normalized aerosol distribution among the lung segments was compared with the normalized flow distribution by calculating the ratio (R(i)) of the number of particles exiting each segmental bronchus i to the flow. When R(i) = 1, particle transport was directly proportional to segmental flow. Flow and particle characteristics were represented by the Stokes number (Stk) in the trachea. For Stk < 0.01, R(i) values were close to 1 and were unaffected by gravity. For Stk > 0.01, R(i) varied greatly among the different outlets (R(i) = 0.30-1.93 in normal gravity for 10-μm particles at 500 ml/s) and was affected by gravity and inertia. These data suggest that, for Stk < 0.01, ventilation defines the delivery of aerosol to lung segments and that the use of aerosol tracers is a valid technique to visualize ventilation in different parts of the lung. At higher Stokes numbers, inertia, but not gravitational sedimentation, is the second major factor affecting the transport of large particles in the lung.     

Convective and diffusive gas transport in canine intrapulmonary airways.

The significance of convective and diffusive gas transport in the respiratory system was assessed from the response of combined inert gas and particle boluses inhaled into the conducting airways. Particles, considered as "nondiffusing gas," served as tracers for convection and two inert gases with widely different diffusive characteristics (He and SF6) as tracers for convection and diffusion. Six-milliliter boluses labeled with monodisperse di-2-ethylhexyl sebacate droplets of 0.86-microns aerodynamic diameter, 2% He, and 2% SF6 were inspired by three anesthetized mechanically ventilated beagle dogs to volumetric lung depths up to 170 ml. Mixing between inspired and residual air caused dispersion of the inspired bolus, which was quantified in terms of the bolus half-width. Dispersion of particles increased with increasing lung depth to which the boluses were inhaled. The increase followed a power law with exponents less than 0.5 (mean 0.39), indicating that the effect of convective mixing per unit volume was reduced with depth. Within the pulmonary dead space, the behavior of the inert gases He and SF6 was similar to that of the particles, suggesting that gas transport was almost solely due to convection. Beyond the dead space, dispersion of He and SF6 increased more rapidly than dispersion of particles, indicating that diffusion became significant. The gas and particle bolus technique offers a suitable approach to differential analysis of gas transport in intrapulmonary airways of lungs.     

Dispersion of 0.5- to 2-micron aerosol in microG and hypergravity as a probe of convective inhomogeneity in the lung.

We used aerosol boluses to study convective gas mixing in the lung of four healthy subjects on the ground (1 G) and during short periods of microgravity (microG) and hypergravity ( approximately 1. 6 G). Boluses of 0.5-, 1-, and 2-micron-diameter particles were inhaled at different points in an inspiration from residual volume to 1 liter above functional residual capacity. The volume of air inhaled after the bolus [the penetration volume (Vp)] ranged from 150 to 1,500 ml. Aerosol concentration and flow rate were continuously measured at the mouth. The dispersion, deposition, and position of the bolus in the expired gas were calculated from these data. For each particle size, both bolus dispersion and deposition increased with Vp and were gravity dependent, with the largest dispersion and deposition occurring for the largest G level. Whereas intrinsic particle motions (diffusion, sedimentation, inertia) did not influence dispersion at shallow depths, we found that sedimentation significantly affected dispersion in the distal part of the lung (Vp >500 ml). For 0.5-micron-diameter particles for which sedimentation velocity is low, the differences between dispersion in microG and 1 G likely reflect the differences in gravitational convective inhomogeneity of ventilation between microG and 1 G.     

Ozone exposure alters tracheobronchial mucociliary function in humans

Mucociliary function is a primary defense mechanism of the tracheobronchial airways, and yet the response of this system to an inhalational hazard, such as ozone, is undefined in humans. Utilizing noninvasive techniques to measure deposition and retention of insoluble radiolabeled particles on airway mucous membranes, we studied the effect on mucus transport of 0.2 and 0.4 ppm ozone compared with filtered air (FA) in seven healthy males. During 2-h chamber exposures, subjects alternated between periods of rest and light exercise with hourly spirometric measurement of lung function. Mechanical and mucociliary function responses to ozone by lung airways appeared concentration dependent. Reduction in particle retention was significant (P less than 0.005) (i.e., transport of lung mucus was increased during exposure to 0.4 ppm ozone and was coincident with impaired lung function; e.g., forced vital capacity and midmaximal flow rate fell by 12 and 16%, respectively, and forced expiratory volume at 1 s by 5%, of preexposure values). Regional analysis indicated that mucus flow from distal airways into central bronchi was significantly increased (P less than 0.025) by 0.2 ppm ozone. This peripheral effect, however, was buffered by only a marginal influence of 0.2 ppm ozone on larger bronchi, such that the resultant mucus transport for all airways of the lung in aggregate differed only slightly from FA exposures. These data may reflect differences in regional diffusion of ozone along the respiratory tract, rather than tissue sensitivity. In conclusion, mucociliary function of humans is acutely stimulated by ozone and may result from fluid additions to the mucus layer from mucosal and submucosal secretory cells and/or alteration of epithelial permeability.     

Bolus contaminant dispersion in oscillating flow in curved tubes.

The investigation of longitudinal dispersion of tracer substances in unsteady flows has biomechanical application in the study of heat and mass transport within the bronchial airways during normal, abnormal, and artificial pulmonary ventilation. To model the effects of airway curvature on intrapulmonary gas transport, we have measured local gas dispersion in axially uniform helical tubes of slight pitch during volume-cycled oscillatory flow. Following a small argon bolus injection into the flow field, the time-averaged effective diffusion coefficient (Deff/Dmol) for axial transport of the contaminant was evaluated from the time-dependent local argon concentration measured with a mass spectrometer. The value of (Deff/Dmol) is extracted from the curve of concentration versus time by two techniques yielding identical results. Experiments were conducted in two helical coiled tubes (delta = 0.031, lambda = 0.022 or delta = 0.085, lambda = 0.060) over a range of 2 < alpha < 15, 3 < A < 15, where delta is the ratio of tube radius to radius of curvature, lambda is the ratio of pitch height to radius of curvature, alpha is the Womersley parameter or dimensionless frequency, and A is the stroke amplitude or dimensionless tidal volume. Experimental results show that, when compared to transport in straight tubes, the effective diffusivity markedly increases in the presence of axial curvature. Results also compare favorably to mathematical predictions of bolus dispersion in a curved tube over the ranges of frequency and tidal volume studied.     

Aerosol bolus dispersion and convective mixing in human and dog lungs and physical models.

The dispersion of aerosol boluses in the lung is a probe for convective mixing and has been proposed as a marker for abnormal lung function. To better understand the factors underlying this phenomenon, aerosol dispersion was compared in human subjects, dogs, and various physical models. In all systems, dispersion increased with the volumetric penetration of the aerosol bolus. The rate of this increase was 83% greater in humans compared with dogs. Dispersion in dogs was close to that in a packed bed with beads of 2.5 mm. Aerosol dispersion decreased with increasing flow rate in human subjects. An artificial larynx inserted into the straight tube caused a 33% increase in dispersion. In humans, aerosol dispersion was significantly correlated with forced expired flow between 25 and 75% of vital capacity. A 2-s pause between inspiration and expiration increased dispersion 23-58% in three isolated dog lungs but did not affect dispersion in the packed bed. The data suggest that lung geometry, flow rate, particle mobility, and the larynx all significantly affect aerosol dispersion by influencing the reversibility of aerosol transport between inspiration and expiration.     

Inspiratory flow rate and dead space in dogs

It is generally accepted that a stationary concentration front is established in the tracheobronchial tree during the inspiratory phase of single- and multiple-breath washouts. The anatomic position of this front, which is determined by the balance between diffusive flux toward the airway opening and convective flux toward the periphery, is frequently used to predict the effects of molecular diffusivity and inspiratory flow rate on dead space. Although there is substantial experimental evidence supporting the predictive effect of molecular diffusivity, there is little evidence regarding the effect of convective flow. This study confirmed the predictions for the effects of molecular diffusivity but contradicted those for the effects of inspiratory flow. We measured dead space by multiple- and single-breath inert gas washout techniques and also measured physiological dead space in dogs for inspiratory flow rates of 10-71 ml.kg-1.s-1. None of the three measures of dead space increased over the entire range of flow rate, as predicted by contemporary gas transport models. A possible explanation for these findings is that axial dispersion coefficients in the anatomic region where stationary fronts are believed to develop (respiratory bronchioles and alveolar ducts) significantly increase with convective flow rate rather than remain equal to molecular diffusivity.     

Inertial deposition effects: a study of aerosol mechanics in the trachea using laser Doppler velocimetry and fluorescent dye

This study characterizes the axial velocity and axial turbulence intensity patterns noted in the tracheal portion of a cadaver-based throat model at two different steady flow rates (18.1 and 41.1 LPM.) This characterization was performed using Phase Doppler Interferometry (Laser Doppler Velocimetry). Deposition, as assessed qualitatively using fluorescent dye, is related to the position of the laryngeal jet within the trachea. The position of the jet is dependent on the downstream conditions of the model. It is proposed therefore that lung/airway conditions may have important effects on aerosol deposition within the throat. There is no correspondence noted between regions of high axial turbulence intensity and deposition.     

A model of constant-flow ventilation in a dog lung

A semiempirical model of constant-flow ventilation (CFV) is developed to test the hypothesis that a three-zone serial model with the following characteristics can explain the adequate CO2 transport observed during CFV: 1) a zone of jet recirculation immediately downstream of the catheter in which convection dominates; 2) a zone influenced by turbulence but with little or no bulk flow; and 3) a peripheral zone, free of turbulence, in which transport is governed by molecular and augmented diffusion. Interactions between turbulent eddies and cardiogenic oscillations are included using a modification of Taylor dispersion theory according to the formulation of Kamm et al. Predicted values for arterial PCO2 are reasonably similar to experimental results for He-O2, air, and SF6-O2 mixtures for catheter flow rates from 0.2 to 1.6 l/s. Specific impedance to gas exchange was found to be largest immediately proximal to the end of turbulent mixing zone, where transport is governed by low-level eddy mixing and molecular diffusion. Simulations suggest that, during CFV, cardiogenic oscillations augment gas exchange primarily by promoting turbulent eddy dispersion in the distal airways and by extending the length of the turbulent mixing zone. Even small displacements of the catheter are shown to have a dramatic effect on gas exchange.     

Computational model of particle deposition in the nasal cavity under steady and dynamic flow

A computational model for flow and particle deposition in a three-dimensional representation of the human nasal cavity is developed. Simulations of steady state and dynamic airflow during inhalation are performed at flow rates of 9–60 l/min. Depositions for particles of size 0.5–20 μm are determined and compared with experimental and simulation results from the literature in terms of deposition efficiencies. The nasal model is validated by comparison with experimental and simulation results from the literature for particle deposition under steady-state flow. The distribution of deposited particles in the nasal cavity is presented in terms of an axial deposition distribution as well as a bivariate axial deposition and particle size distribution. Simulations of dynamic airflow and particle deposition during an inhalation cycle are performed for different nasal cavity outlet pressure variations and different particle injections. The total particle deposition efficiency under dynamic flow is found to depend strongly on the dynamics of airflow as well as the type of particle injection.     

Diffusivity dependence of multiple-breath washouts of lung periphery

Subpleural concentrations of He and SF6 were measured during multiple-breath washouts from isolated dog lungs. Tidal volume, inspiratory flow, and frequency were in the normal range of canine ventilation. For each gas, there was a local minimum in concentration during inspiration (Cinsp) and a local maximum in concentration during exhalation (Cexp). SF6 exhibited a deeper inspiratory trough than He for each breath of every washout. For large tidal volumes (10-20 ml/kg), Cexp approximated a single exponential decay and He was cleared more rapidly than SF6. For small tidal volumes (2.5 ml/kg), Cexp was multiexponential and SF6 was cleared more rapidly than He. Cinsp/Cexp (a measure of the depth of the inspiratory trough) and the kinetics of Cexp decay were determined for washouts using a tidal volume of 10 and 20 ml/kg and different inspiratory flows. Under all conditions, an increase of inspiratory flow resulted in a deeper inspiratory trough for both He and SF6. For washouts using 10 ml/kg and 60 breaths/min, an increase of inspiratory flow increased the clearance of both gases. In washouts using lower ventilatory frequencies, gas clearance was independent of inspiratory flow. These findings are contrary to predictions of contemporary models of convection and diffusion in the lung. This study suggests that convective axial mixing and radial diffusion in the airways are important determinants of pulmonary gas transport.     

Pressure change and gas mixing induced by oscillations in a closed system

In an attempt to delineate some mechanical behaviors found in branching airways, pressure transmission, gas motion, and mixing were studied during high-frequency oscillation (HFO) in an idealized system consisting of a large straight tube and a rigid sphere linked together by a small straight tube. Depending on the frequency f, and on the unsteadiness dimensionless parameter alpha, pressure amplitude in the large tube is either strongly attenuated or amplified in the sphere. This finding may provide a theoretical basis for the pressure resonance phenomenon observed in the lung by previous investigators. Gas compression in the closed volume causes convective mixing throughout the system. The measured dispersion was found to be proportional to f(VT/A)2, in agreement with a recent report. However, bulk convective mixing was sufficient to explain the dispersion for oscillatory volumes (VT) as small as 80 percent of the small tube volume, as has been previously suggested.     

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.     

Low reynolds number viscous flow in an alveolated duct

Flow visualization studies and supplementary numerical simulations are carried out on slow flow through a model alveolated duct. The results reveal that the type of flow that develops in the alveoli, or cavities, is controlled by the ratio of the depth to the width of the cavity and by the ratio of cavity volume to duct volume. While weak, the slowly rotating flow in the cavity is thought to be important to the convective transport of heat and mass transfer to, or from, the walls of the cavity. The relevance of these finding to particle transport and deposition deep in the lung is discussed.     

Acinar flow irreversibility caused by perturbations in reversible alveolar wall motion.

Mixing associated with "stretch-and-fold" convective flow patterns has recently been demonstrated to play a potentially important role in aerosol transport and deposition deep in the lung (J. P. Butler and A. Tsuda. J. Appl. Physiol. 83: 800-809, 1997), but the origin of this potent mechanism is not well characterized. In this study we hypothesized that even a small degree of asynchrony in otherwise reversible alveolar wall motion is sufficient to cause flow irreversibility and stretch-and-fold convective mixing. We tested this hypothesis using a large-scale acinar model consisting of a T-shaped junction of three short, straight, square ducts. The model was filled with silicone oil, and alveolar wall motion was simulated by pistons in two of the ducts. The pistons were driven to generate a low-Reynolds-number cyclic flow with a small amount of asynchrony in boundary motion adjusted to match the degree of geometric (as distinguished from pressure-volume) hysteresis found in rabbit lungs (H. Miki, J. P. Butler, R. A. Rogers, and J. Lehr. J. Appl. Physiol. 75: 1630-1636, 1993). Tracer dye was introduced into the system, and its motion was monitored. The results showed that even a slight asynchrony in boundary motion leads to flow irreversibility with complicated swirling tracer patterns. Importantly, the kinematic irreversibility resulted in stretching of the tracer with narrowing of the separation between adjacent tracer lines, and when the cycle-by-cycle narrowing of lateral distance reached the slowly growing diffusion distance of the tracer, mixing abruptly took place. This coupling of evolving convective flow patterns with diffusion is the essence of the stretch-and-fold mechanism. We conclude that even a small degree of boundary asynchrony can give rise to stretch-and-fold convective mixing, thereby leading to transport and deposition of fine and ultrafine aerosol particles deep in the lung.