A MONTE-CARLO CALCULATION OF THE DEPOSITION EFFICIENCY OF INHALED PARTICLES IN LOWER AIRWAYS

Steady airflow and particle transport in a straight cylindrical tube a nd in a three-dimensional model of a bifurcating airway were calculate d. The model domains were based on human lung data. For simplicity, ai rways in the bifurcating model were assumed to have a rectangular cros s-section with hydraulic diameters similar to those of generations 16 and 17 of the Weibel model. Steady flow with a Reynolds number of unit y within this domain was solved numerically by a finite element method using a commercial software package. Parabolic inlet conditions were assumed. Using the calculated flow profile, individual particle trajec tories within the air stream were simulated by numerically solving the fundamental differential equations of particle motion. The simulation included the effects of diffusion and sedimentation. Deposition effic iencies by individual mechanisms and the total deposition efficiency w ere obtained using a Monte Carlo method. For parabolic flow in a strai ght tube, the results compared very well with existing formulas in the literature. For both geometries, the sum of the efficiencies due to i ndividual mechanisms was greater than the numerical prediction of tota l efficiency. The total deposition efficiency, eta(total), was accurat ely expressed in terms of the efficiencies of sedimentation, eta(s), a nd diffusion, eta(d), using the empirical formula eta(total) = (eta(d) p + eta(s)p)1/p. The best value of p depended on the particular geomet ry.