Direct Numerical Simulations of Inhalation in a 23-Generation Lung Model

The air flows in the proximal and distal portions of the human lungs are interconnected: the lower Reynolds number in the deeper generations causes a progressive flow regularization, but mass conservation requires flow rate oscillations to propagate through the airway bifurcations. To explain how these competing effects shape the flow state in the deeper generations, we have performed the first high-fidelity numerical simulation of the air flow in a lung model that includes 23 successive bifurcations of a single planar airway. Turbulence modeling or assumptions on flow regimes are not required. The chosen flow rate is steady on average, and representative of the peak inspiratory flow reached by adult patients breathing through therapeutic inhalers. As expected, advection becomes progressively less important after each bifurcation, until a time-dependent Stokes regime governed solely by viscous diffusion is established in the smallest generations. However, fluctuations in this regime are relatively fast and large with respect to the mean flow, which is in contrast with the commonly agreed picture that only the breathing frequency is relevant at the scale of the alveoli. We demonstrate that the characteristic frequency and amplitude of these fluctuations are linked to the flow in the upper part of the bronchial tree, as they originate from the time-dependent flow splitting in the upper bifurcations. Even though these fluctuations are observed here in an idealized, rigid lung model, our findings suggest that the assumptions usually adopted in many of the current lung models might need to be revised.