The understanding and quantitative assessment of air flow fields and local micron-particle wall concentrations in tracheobronchial airways are very important for estimating the health risks of inhaled particulate pollutants, developing algebraic transfer functions of global lung deposition models used in dose-response analyses, and/or determining proper drug-aerosol delivery to target sites in the lung. In this paper (Part 1) the theory, model geometries, and air flow results are provided. In a companion paper (Part 2, Comer et al. 2001), the history of particle deposition patterns and comparisons with measured data sets are reported. Decoupling of the naturally dilute particle suspension makes it feasible to present the results in two parts.Considering a Reynolds number range of 500 [les ] ReD [les ] 2000, it is assumed that the air flow is steady, incompressible and laminar and that the tubular double bifurcations, i.e. Weibel's generations G3–G5, are three-dimensional, rigid, and smooth with rounded as well as sharp carinal ridges for symmetric planar, and just rounded carinas for 90° non-planar configurations. The employed finite-volume code CFX (AEA Technology) and its user-enhanced FORTRAN programs were validated with experimental velocity data points for a single bifurcation. The resulting air flow structures are analysed for relatively low (ReD = 500) and high (ReD = 2000) Reynolds numbers. Sequential pressure drops due to viscous effects were calculated and compared, extending a method proposed by Pedley et al. (1977). Such detailed results for bifurcating lung airways are most useful in the development of global algebraic lung models.

The NASA Goddard Space Flight Center (GSFC) Scanning Raman Lidar (SRL) participated in the International H2O Project (IHOP), which occurred in May and June 2002 in the midwestern part of the United States. The SRL received extensive optical modifications prior to and during the IHOP campaign that added new measurement capabilities and enabled unprecedented daytime water vapor measurements by a Raman lidar system. Improvements were also realized in nighttime upper-tropospheric water vapor measurements. The other new measurements that were added to the SRL for the IHOP deployment included rotational Raman temperature, depolarization, cloud liquid water, and cirrus cloud ice water content. In this first of two parts, the details of the operational configuration of the SRL during IHOP are provided along with a description of the analysis and calibration procedures for water vapor mixing ratio, aerosol depolarization, and cirrus cloud extinction-to-backscatter ratio. For the first time, a Raman water vapor lidar calibration is performed, taking full account of the temperature sensitivity of water vapor and nitrogen Raman scattering. Part II presents case studies that permit the daytime and nighttime error statistics to be quantified.

No abstract

The flow theory and air flow structures in symmetric double-bifurcation airway models assuming steady laminar, incompressible flow, unaffected by the presence of aerosols, has been described in a companion paper (Part 1). The validated computer simulation results showed highly vortical flow fields, especially around the second bifurcations, indicating potentially complex particle distributions and deposition patterns. In this paper (Part 2), assuming spherical non-interacting aerosols that stick to the wall when touching the surface, the history of depositing particles is described. Specifically, the finite-volume code CFX (AEA Technology) with user-enhanced FORTRAN programs were validated with experimental data of particle deposition efficiencies as a function of the Stokes number for planar single and double bifurcations. The resulting deposition patterns, particle distributions, trajectories and time evolution were analysed in the light of the air flow structures for relatively low (ReD1 = 500) and high (ReD1 = 2000) Reynolds numbers and representative Stokes numbers, i.e. StD1 = 0.04 and StD1 = 0.12. Particle deposition patterns and surface concentrations are largely a function of the local Stokes number, but they also depend on the fluid–particle inlet conditions as well as airway geometry factors. While particles introduced at low inlet Reynolds numbers (e.g. ReD1 = 500) follow the axial air flow, secondary and vortical flows become important at higher Reynolds numbers, causing the formation of particle-free zones near the tube centres and subsequently elevated particle concentrations near the walls. Sharp or mildly rounded carinal ridges have little effect on the deposition efficiencies but may influence local deposition patterns. In contrast, more drastic geometric changes to the basic double-bifurcation model, e.g. the 90°-non-planar configuration, alter both the aerosol wall distributions and surface concentrations considerably.

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.