Numerical investigation of the three-dimensional flow in a human lung model

Abstract

The flow field at inspiration and expiration in the upper human airways consisting of the trachea down to the sixth generation of the bronchial tree is numerically simulated. The three-dimensional steady flow at a hydraulic diameter-based Reynolds number Re D=1250 is computed via a lattice-Boltzmann method (LBM). The simulation is validated by the experimental data based on particle-image velocimetry (PIV) measurements. The good agreement between numerical and experimental results is evidenced by comparing velocity contours and distributions in a defined reference plane. The results show the LBM to be an accurate tool to numerically predict flow structures in the human lung. Using an automatic Cartesian grid generator, the overall process time from meshing to a steady-state solution is <12 h. Moreover, the numerical simulation allows a closer analysis of the secondary flow structures than in the experimental investigation. The three-dimensional streamline patterns reveal some insight on the air exchange mechanism at inspiration and expiration. At inspiration, the slower near-wall tracheal flow enters through the right principal bronchus into the right upper lobar bronchus. The bulk mass flux in the trachea is nearly evenly distributed over the left upper, center and lower lobar bronchi and the right center and lower bronchi. At expiration, the air from the right upper lobar bronchus enters the right center of the trachea and displaces the airflow from the lower and center right bronchi such that the tracheal positions of the streamlines at inspiration and expiration are switched. The flow in the left bronchi does not show this kind of switching. The findings emphasize the impact of the asymmetry of the lung geometry on the respiratory air exchange mechanism.