Monte Carlo simulation of simultaneous gas flow and diffusion in an asymmetric distal pulmonary airway model

Abstract

In order to evaluate the effect of anatomic asymmetries on the gas concentration distribution in the pulmonary airways, a Monte Carlo simulation of combined bulk flow and molecular diffusion was carried out in a realistic distal airway model (Parkeret al., 1971). This airway model, composed of branches distal to the 0.5-ram diameter airways, contained an upper symmetric segment consisting of four generations of conducting airways and a lower asymmetric segment of alveolar ducts and sacs arranged in five transport paths of varying lengths. In accounting for the volume increases of these ducts and sacs occurring during normal respiration, uniform alveolar filling rates and a fixed length-to-diameter ratio of all airways were assumed. For a pulse injection of inert tracer gas, the simulation was employed to determine the longitudinal concentration profiles in the conducting airways. In the alveolated airways, not only were the longitudinal profiles determined along each path, but radial transport from the core to the periphery of the airways was considered.

The results of the simulations indicate that geometric asymmetries alone contribute substantially to regional concentration variations in the distal airways. For example, when a gas bolus is injected at mid*inspiration, there are concentration differences as great as 40% between two points along different transport paths located equi-distant from the proximal end of the model. As viewed from the terminal end of the model (acinus), average concentration differences as large as 6-to-1 exist between the longest and shortest transport paths respectively for gas boli introduced near the end of inspiration. The results further indicate because of large radial diffusion rates, no significant concentration differences exist between the periphery a-ld the central core of alveolated airways.

Simulation of the expired concentration profiles indicate that boll injected very late during inspiration exhibit a sloping tail, unlike the earlier injected boll whose tails are virtually horizontal. Through the use of superposition teehniqnes, it was found that these sloping tails correspond to an alveolar slope of 1.5 vol% between 750 and 1250 ml expired for a continuous washing of tracer. This result is in disagreement with other transport analyses which did not directly account for the effect of geometric asymmetries.