Particle deposition in a CT-scanned human lung airway

The particle deposition in a computerized tomography (CT)-scanned human lung was numerically investigated. The five-generation airway is extracted from the trachea to segmental bronchi of a 60-year-old Chinese male patient. Computations were carried out in the flow rate range of 210–630 ml/s (Reynolds number range of 1000–3000) and particle size of 2–10 μm (Stokes number range of 0.0007–0.049). To count the effect of laryngeal jet on trachea inlet, the trachea was extended and modified to simulate the larynx, consequently the inlet velocity profile is biased towards the rear wall. The laryngeal jet-induced turbulence was simulated using low Reynolds number (LRN) κ–ω turbulent model. Particle deposition patterns, deposition efficiency and deposition factor were studied in detail. The turbulent flow has significant effect on the particle deposition, and the present deposition factor is compared well with the available data.     

Airway branching morphology of mature and immature rabbit lungs.

The scheme of Horsfield et al. for describing the pulmonary airway tree (J Appl Physiol 52: 21-26, 1982) catalogs each airway according to its order and the difference in order of its two daughters (denoted Delta). Although this scheme captures the natural asymmetry in the airway tree, it is still deterministic, because it assumes that all airways of a given order are the same; yet such variability is extremely important in determining the overall behavior of the lungs. We therefore analyzed complete lung lobes from three mature and two immature rabbits and determined the Horsfield order and Delta of every airway down to the terminal bronchioles. We also measured the diameter of each airway. This allowed us to determine the average structure of the rabbit airway tree, the variation about this average, and also how the structures of mature and immature airway trees compare. We found some variation in branching asymmetry and airway diameter at a given order between animals but no evidence of systematic differences in structure between mature and immature lungs. We found evidence of a difference in the branching structure of the peripheral vs. the central part of the airway tree (the break point being around order 20). We also determined the nature of the variation in Delta and diameter as a function of order, which should be valuable for the development of computer models seeking to encapsulate the naturally occurring regional variation in airway geometry in the normal rabbit lung.     

A model of the mechanics of airway narrowing

To examine the interaction between airway smooth muscle shortening and airway wall thickening on changes in pulmonary resistance, we have developed a model of the tracheobronchial tree that allows simulation of the mechanisms involved in airway narrowing. The model is based on the symmetrical dichotomous branching tracheobronchial tree as described by Weibel and uses fluid dynamic equations proposed by Pedley et al. to calculate inspiratory resistance during quiet tidal breathing. To allow for changes in lung volume, we used the airway pressure-area curves developed by Lambert et al. The model is easily implemented with a spreadsheet and personal computer that allows calculation of total and regional pulmonary resistance. At each airway generation in the model, provision is made for airway wall thickness, the maximal airway smooth muscle shortening achievable, and an S-shaped dose-response relationship to describe smooth muscle shortening. To test the validity of the model, we compared pressure-flow curves generated with the model with measurements of pulmonary resistance while normal subjects breathed air and 20% O2-80% He at a variety of lung volumes. By simulating progressive airway smooth muscle shortening, realistic pulmonary resistance vs. dose-response curves were produced. We conclude that this model provides realistic estimates of pulmonary resistance and shows potential for examining the various mechanisms that could produce excessive airway narrowing in disease.     

Human variation in the peripheral air-space deposition of inhaled particles

Intersubject variability in both peripheral air-space dimensions and breathing pattern [tidal volume (VT) and respiratory frequency (f)] may play a role in determining intersubject variation in the fractional deposition of inhaled particles that primarily deposit in the lung periphery (i.e., distal to conducting airways). In healthy subjects breathing spontaneously at rest, we measured the deposition fraction (DF) of a 2.6-microns monodisperse aerosol by Tyndallometry while simultaneous measurement of VT and f were made. Under these conditions particle deposition occurs primarily in the peripheral air spaces of the lung. As an index of peripheral air-space size, we used measurements of aerosol recovery (RC) as a function of breath-hold time (t) (Gebhart et al. J. Appl. Physiol. 51: 465-476, 1981). In each subject, we measured RC (aerosol expired/aerosol inspired) of a 1.0-micron monodisperse aerosol as a function of breath-hold time for inspiratory capacity breaths of aerosol. The half time (t1/2) (the breath-hold time to reach 50% RC with no breath hold) is proportional to a mean diameter (D) of air spaces filled with aerosol. In the 10 subjects studied, we found a variable DF, range 0.04-0.44 [0.25 +/- 0.12 (SD)]. DF correlated most closely with 1/f, or the period of breathing (r = 0.96, P less than 0.01). There was no significant correlation between DF and t1/2 as an index of peripheral air-space size. In fact there was little deviation in t1/2 in these normal subjects [coefficient of variation (CV) = 0.12].(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

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.     

20-hydroxyeicosatetraenoic acid (20-HETE): structural determinants for renal vasoconstriction

The effects of natural and synthetic eicosanoids on the diameter of rat interlobular arteries studied in vitro were compared to that of the potent, endogenous vasoconstrictor 20-HETE. Vasoconstrictor activity was optimum for chain lengths of 20-22 carbons with at least one olefin or epoxide between located between C(13)-C(15) and an oxygen substituent at C(20)-C(22). The presence of delta (Zou et al. Am. J. Physiol. 1996, 270, R228; Gebremedhin, D. et al. Am. J. Physiol. 1998, 507, 771)-, delta (Carroll et al. Am. J. Physiol. 1996, 271, R863; Vazquez et al. Life Sci. 1995, 56, 1455)-, or delta (Imig et al. Hypertension 2000, 35, 307; Lopez et al. Amer. J. Physiol. 2001, 281, F420)-olefins had no influence on the vasoconstrictor response whereas the introduction of a C(7)-thiomethylene enhanced potency. A sulfonamide or alcohol, but not a lactone, could replace the C(1)-carboxylate. These data were used to construct a putative binding domain map of the 20-HETE receptor consisting of: (i) a comparatively open, hydrophilic binding site accommodating the C(1)-functionality; (ii) a hydrophobic trough spanning the olefins; (iii) a shallow pocket containing a critical pi-pi binding site in the vicinity of the pi (Ito et al. Am. J. Physiol. 1998, 274, F395; Quigley, R.; Baum, M.; Reddy, K. M.; Griener, J. C.; Falck, J. R. Am. J. Physiol. 2000, 278, F949)-olefin; and (iv) an oxyphilic binding site proximate to the omega-terminus.     

Geometric hysteresis of alveolated ductal architecture

Low Reynolds number airflow in the pulmonary acinus and aerosol particle kinetics therein are significantly conditioned by the nature of the tidal motion of alveolar duct geometry. At least two components of the ductal structure are known to exhibit stress-strain hysteresis: smooth muscle within the alveolar entrance rings, and surfactant at the air-tissue interface. We hypothesize that the geometric hysteresis of the alveolar duct is largely determined by the interaction of the amount of smooth muscle and connective tissue in ductal rings, septal tissue properties, and surface tension-surface area characteristics of surfactant. To test this hypothesis, we have extended the well-known structural model of the alveolar duct by Wilson and Bachofen (1982, "A Model for Mechanical Structure of the Alveolar Duct," J. Appl. Physiol. 52(4), pp. 1064-1070) by adding realistic elastic and hysteretic properties of (1) the alveolar entrance ring, (2) septal tissue, and (3) surfactant. With realistic values for tissue and surface properties, we conclude that: (1) there is a significant, and underappreciated, amount of geometric hysteresis in alveolar ductal architecture; and (2) the contribution of smooth muscle and surfactant to geometric hysteresis are of opposite senses, tending toward cancellation. Quantitatively, the geometric hysteresis found experimentally by Miki et al. (1993, "Geometric Hysteresis in Pulmonary Surface-to-Volume Ratio during Tidal Breathing," J. Appl. Physiol. 75(4), pp. 1630-1636) is consistent with little or no smooth muscle tone in anesthetized rabbits in control conditions, and with substantial smooth muscle activation following methacholine challenge. The observed local hysteretic boundary motion of the acinar duct would result in irreversible acinar flow fields, which might be important mechanistic contributors to aerosol mixing and deposition deep in the lung.     

Response of normal subjects to inspiratory resistive unloading

The purpose of this study was to examine the role of the normal inspiratory resistive load in the regulation of respiratory motor output in resting conscious humans. We used a recently described device (J. Appl. Physiol. 62: 2491-2499, 1987) to make mouth pressure during inspiration positive and proportional to inspiratory flow, thus causing inspiratory resistive unloading (IRUL); the magnitude of IRUL (delta R = -3.0 cmH2O.1(-1).s) was set so as to unload most (approximately 86% of the normal inspiratory resistance. Six conscious normal humans were studied. Driving pressure (DP) was calculated according to the method of Younes et al. (J. Appl. Physiol. 51: 963-1001, 1981), which provides the equivalent of occlusion pressure at functional residual capacity throughout the breath. IRUL resulted in small but significant changes in minute ventilation (0.6 1/min) and in end-tidal CO2 concentration (-0.11%) with no significant change in tidal volume or respiratory frequency. There was a significant shortening of the duration (neural inspiratory time) of the rising phase of the DP waveform and the shape of the rising phase became more convex to the time axis. There was no change in the average rate of rise of DP or in the duration or shape of the declining phase. We conclude that 1) the normal inspiratory resistance is an important determinant of the duration and shape of the rising phase of DP and 2) the neural responses elicited by the normal inspiratory resistance are similar to those observed with added inspiratory resistive loads.     

Pulmonary mechanics during rapid mechanical ventilation in rabbits with saline-lavaged lungs

Infants with respiratory failure are frequently mechanically ventilated at rates exceeding 60 breaths/min. We analyzed the effect of ventilatory rates of 30, 60, and 90 breaths/min (inspiratory times of 0.6, 0.3, and 0.2 s, respectively) on the pressure-flow relationships of the lungs of anesthetized paralyzed rabbits after saline lavage. Tidal volume and functional residual capacity were maintained constant. We computed effective inspiratory and expiratory resistance and compliance of the lungs by dividing changes in transpulmonary pressure into resistive and elastic components with a multiple linear regression. We found that mean pulmonary resistance was lower at higher ventilatory rates, while pulmonary compliance was independent of ventilatory rate. The transpulmonary pressure developed by the ventilator during inspiration approximated a linear ramp. Gas flow became constant and the pressure-volume relationship linear during the last portion of inspiration. Even at a ventilatory rate of 90 breaths/min, 28-56% of the tidal volume was delivered with a constant inspiratory flow. Our findings are consistent with the model of Bates et al. (J. Appl. Physiol. 58: 1840-1848, 1985), wherein the distribution of gas flow within the lungs depends predominantly on resistive factors while inspiratory flow is increasing, and on elastic factors while inspiratory flow is constant. This dynamic behavior of the surfactant-depleted lungs suggests that, even with very short inspiratory times, distribution of gas flow within the lungs is in large part determined by elastic factors. Unless the inspiratory time is further shortened, gas flow may be directed to areas of increased resistance, resulting in hyperinflation and barotrauma.     

Micron-sized intrapulmonary particle deposition in the developing rat lung

Little is known about the effects of postnatal developmental changes in lung architecture and breathing patterns on intrapulmonary particle deposition. We measured deposition in the developing Wistar-Kyoto rat, whose lung development largely parallels that of humans. Deposition of 2-μm sebacate particles was determined in anesthetized, intubated, spontaneously breathing rats on postnatal days (P) 7 to 90 by aerosol photometry (Karrasch S, Eder G, Bolle I, Tsuda A, Schulz H. J Appl Physiol 107: 1293-1299, 2009). Respiratory parameters were determined by body plethysmography. Tidal volume increased substantially from P7 (0.19 ml) to P90 (2.1 ml) while respiratory rate declined from 182 to 107/min. Breath-specific deposition was lowest (9%) at P7 and P90 and markedly higher at P35 (almost 16%). Structural changes of the alveolar region include a ninefold increase in surface area (Bolle I, Eder G, Takenaka S, Ganguly K, Karrasch S, Zeller C, Neuner M, Kreyling WG, Tsuda A, Schulz H. J Appl Physiol 104: 1167-1176, 2008). Particle deposition per unit of time and surface area peaked at P35 and showed a minimum at P90. At an inhaled particle number concentration of 10(5)/cm(3), there was an estimated 450, 690, and 330 particles/(min × cm(2)) at P7, P35, and P90, respectively. Multiple regression models showed that deposition depends on the mean linear intercept as structural component and the breathing parameters, tidal volume, and respiratory rate (r(2) > 0.9). In conclusion, micron-sized particle deposition was dependent on the stage of postnatal lung development. A maximum was observed during late alveolarization (P35), which corresponds to human lungs of about eight years of age. Children at this age may therefore be more susceptible to micron-sized airborne environmental health hazards.     

Significance of airway resistance for the pattern of breathing and lung volumes in exercising humans

The effects of increased airway resistance on lung volumes and pattern of breathing were studied in eight subjects performing leg exercise on a cycle ergometer. Airway resistance was changed 1) by increasing the density (D) of the respired gas by a factor of 4.2 and changing the inspired gas from O2 at 1.3 bar to air at 6 bar and 2) by increasing airway flow rates by exposing the subjects to incremental work loads of 0-200 W. Increased gas D caused a slower and deeper respiration at rest and during exercise and, at work loads greater than 120 W, depressed the responses of ventilation and mean inspiratory flow. Raised airway resistance induced by increases in D and/or airway flow rates altered respiratory timing by increasing the ratio of inspiratory time (TI) to total breath duration. Furthermore, analyses of the relationships between tidal volume and TI and between end-inspiratory volume and TI revealed elevation of Hering-Breuer inspiratory volume thresholds. We propose that this elevation, and hence exercise-induced increases of tidal volume, can largely be explained by previous observations that the threshold of the inspiratory off-switch mechanisms depends on central inspiratory activity (cf. C. von Euler, J. Appl. Physiol. 55: 1647-1659, 1983), which in turn increases with airway resistance (Acta Physiol. Scand. 120: 557-565, 1984).     

Modeling the bifurcating flow in a CT-scanned human lung airway

The inspiratory flow characteristics in a CT-scanned human lung model were numerically investigated using low Reynolds number (LRN) kappa-omega turbulent model. The five-generation airway is extracted from the trachea to segmental bronchi of a 60-year-old Chinese male patient. Computations were carried out in the Reynolds number range of 900-2100, corresponding to mouth-air breathing rates of 190-440ml/s. Flow patterns on the Re=2100 and flow rate distribution were presented. In this model, the flow pattern is very complex. To count the effect of laryngeal jet on trachea inlet, the trachea was extended and modified to simulate the larynx, consequently the inlet velocity profile is biased towards the rear wall. In the inferior lobar bronchi, there are two stems in which the axial velocity is stronger but secondary velocity is weaker. Secondary flow in the lateral bronchi is stronger than the medial ones. With increasing Re, the air flow increases in the middle, inferior lobes and left main bronchus, i.e., flow biases to left and downward.     

Properties of steady maximal expiratory flow within excised canine central airways

Using our transistor model of the lung during forced expiration (J. Appl. Physiol. 62: 2013-2025, 1987), we recently predicted that 1) axially arranged choke points can exist simultaneously during forced expiration with sufficient effort, and 2) overall maximal expiratory flow may be relatively insensitive to nonuniform airways obstruction because of flow interdependence between parallel upstream branches. We tested these hypotheses in excised central airways obtained from five canine lungs. Steady expiratory flow was induced by supplying constant upstream pressure (Pupstream = 0-16 cmH2O) to the bronchi of both lungs while lowering pressure at the tracheal airway opening (16 to -140 cmH2O). Intra-airway pressure profiles obtained during steady maximal expiratory flow disclosed a single choke point in the midtrachea when Pupstream was high (2-16 cmH2O). However, when Pupstream was low (0 cmH2O), two choke sites were evident: the tracheal site persisted, but another upstream choke point (main carina or both main bronchi) was added. Flow interdependence was studied by comparing maximal expiratory flow through each lung before and after introduction of a unilateral external resistance upstream of the bronchi of one lung. When this unilateral resistance was added, ipsilateral flow always fell, but changes in flow through the contralateral lung depended on the site of the most upstream choke. When a single choke existed in the trachea, addition of the external resistance increased contralateral flow by 38 +/- 28% (SD, P less than 0.003). In contrast, when the most upstream choke existed at the main carina or in the bronchi, addition of the external resistance had no effect on contralateral maximal expiratory flow.(ABSTRACT TRUNCATED AT 250 WORDS)     

Deposition Velocity of PM2.5 in the Winter and Spring above Deciduous and Coniferous Forests in Beijing,China

To estimate the deposition effect of PM2.5 (particle matter with aerodynamic diameter <2.5 µm) in forests in northern China, we used the gradient method to measure the deposition velocity of PM2.5 during the winter and spring above a deciduous forest in Olympic Forest Park and above a coniferous forest in Jiufeng National Forest Park. Six aerosol samplers were placed on two towers at each site at heights of 9, 12 and 15 m above the ground surface. The sample filters were exchanged every four hours at 6∶00 AM, 10∶00 AM, 2∶00 PM, 6∶00 PM, 10∶00 PM, and 2∶00 AM. The daytime and nighttime deposition velocities in Jiufeng Park and Olympic Park were compared in this study. The February deposition velocities in Jiufeng Park were 1.2±1.3 and 0.7±0.7 cm s⁻¹ during the day and night, respectively. The May deposition velocities in Olympic Park were 0.9±0.8 and 0.4±0.5 cm s⁻¹ during the day and night, respectively. The May deposition velocities in Jiufeng Park were 1.1±1.2 and 0.6±0.5 cm s⁻¹ during the day and night, respectively. The deposition velocities above Jiufeng National Forest Park were higher than those above Olympic Forest Park. The measured values were smaller than the simulated values obtained by the Ruijgrok et al. (1997) and Wesely et al. (1985) models. However, the reproducibility of the Ruijgrok et al. (1997) model was better than that of the Wesely et al. (1985) model. The Hicks et al. (1977) model was used to analyze additional forest parameters to calculate the PM2.5 deposition, which could better reflect the role of the forest in PM2.5 deposition.     

Effect of increased static lung recoil on bronchial dimesions of excised lungs.

We measured bronchial diameters and lengths during static deflation and inflation in eight excised dog lobes before and after static lung recoil (Pst(L)) had been significantly increased by cooling the lobe for 48 h at 4 degrees C and ventilating it for 3 h. In control lobes, bronchial diameters were the same at any volume even though Pst(L) was different during inflation and deflation. These results agree with those of Hughes et al. (J. Appl. Physiol. 32: 25-35, 1972). However, when Pst(L) was increased, diameters at a given volume were significantly increased over control values; diameters at a given pressure were nearly the same as the controls. Therefore, under these conditions, bronchial diameter did not conform to lung volume. The ventilation process appeared to alter the circumferential elastic properties of the bronchi because diameters at all pressures were slightly larger after ventilation. Bronchial length-volume relationships were the same in both control and ventilated lobes. Thus, when Pst(L) was markedly increased, diameter corresponded best to lung recoil and length to lung volume.     

Numerical simulation of airflow and microparticle deposition in a synchrotron micro-CT-based pulmonary acinus model

The acinus consists of complex, branched alveolar ducts and numerous surrounding alveoli, and so in this study, we hypothesized that the particle deposition can be much influenced by the complex acinar geometry, and simulated the airflow and particle deposition (density = 1.0 g/cm³, diameter = 1 and 3 μm) numerically in a pulmonary acinar model based on synchrotron micro-CT of the mammalian lung. We assumed that the fluid–structure interaction was neglected and that alveolar flow was induced by the expansion and contraction of the acinar model with the volume changing sinusoidally with time as the moving boundary conditions. The alveolar flow was dominated by radial flows, and a weak recirculating flow was observed at the proximal side of alveoli during the entire respiratory cycle, despite the maximum Reynolds number at the inlet being 0.029. Under zero gravity, the particle deposition rate after single breathing was less than 0.01, although the particles were transported deeply into the acinus after inspiration. Under a gravitational field, the deposition rate and map were influenced strongly by gravity orientation. In the case of a particle diameter of 1 μm, the rate increased dramatically and mostly non-deposited particles remained in the model, indicating that the rate would increase further after repeated breathing. At a particle diameter of 3 μm, the rate was 1.0 and all particles were deposited during single breathing. Our results show that the particle deposition rate in realistic pulmonary acinar model is higher than in an idealized model.     

Increased lung volume limits endurance of inspiratory muscles

We examined the influence of lung volume on the ability of normal subjects to sustain breathing against inspiratory resistive loading. Four normal subjects breathed on a closed circuit in which inspiration was loaded by a flow resistor. Subjects were assigned a series of breathing tasks over a range of pressures and flows. In each task there was a specified resistor and also targets for either mean esophageal or airway opening pressure, respiratory frequency, and duty cycle. Endurance was assessed as the length of time to failure of the assigned task. The prime experimental variable was lung volume, which was increased by approximately 1 liter during some tasks; 8 cmH2O continuous positive airway pressure was applied to increase lung volume without increasing elastic load. As previously shown (McCool et al.J. Appl. Physiol. 60: 299-303, 1986), for tasks that could be sustained for the same time, there was an inverse linear relationship of mean esophageal pressure with inspiratory flow rate. This trade-off of pressure and flow was apparent both with and without the increase of lung volume. Comparable tasks, however, could not be sustained as long at the higher lung volumes. This effect of volume on endurance was greater for tasks characterized by high inspiratory pressures and low flow rates than for tasks that could be sustained for the same time but that had lower inspiratory pressures and higher flow rates. This is probably due to the effects of shortening of the sarcomere on fatiguability. Increased lung volume, per se, may contribute to respiratory failure because of increased inspiratory muscle fatiguability by mechanisms independent of elastic load.