Simulation of ozone uptake distribution in the human airways by orthogonal collocation on finite elements

Ozone transport in a rigid single-pathway anatomic model of the lung was analyzed by a stable convergent numerical algorithm, the method of orthogonal collocation on finite elements. The simulations predicted the dynamic behavior of gas phase concentration profiles for both ozone and an insoluble inert gas. An internal quasi-stationary diffusion front was observed during early inspiration for both gases. In addition, the absorptive distribution of ozone in lung airways was computed as a total dose as well as a tissue dose. The total dose of ozone decreased along the airway path from the mouth. However, the tissue dose of ozone increased along the conducting airways, reached a maximal dose in the terminal bronchioles, and decreased sharply in the respiratory airways.