Computational analysis of human upper airway aerodynamics

There is a considerable interest in understanding transient human upper airway aerodynamics, especially in view of assessing the effects of various ventilation therapies. Experimental analyses in a patient-specific manner pose challenges as the upper airway consists of a narrow confined region with complex anatomy. Pressure measurements are feasible, but, for example, PIV experiments require special measures to accommodate for the light refraction by the model. Computational fluid dynamics can bridge the gap between limited experimental data and detailed flow features. This work aims to validate the use of combined lattice Boltzmann method and a large eddy scale model for simulating respiration, and to identify clinical features of the flow and show the clinical potential of the method. Airflow was computationally analyzed during a realistic, transient, breathing profile in an upper airway geometry ranging from nose to trachea, and the resulting pressure calculations were compared against in vitro experiments. Simulations were conducted on meshes containing about 1 billion cells to ensure accuracy and to capture intrinsic flow features. Airway pressures obtained from simulations and in vitro experiments are in good agreement both during inhalation and exhalation. High velocity pharyngeal and laryngeal jets and recirculation in the region of the olfactory cleft are observed. Graphical Abstract The Lattice-Boltzmann Method combined with Large Eddy Simulations was used to compute the aerodynamics in a human upper airway geometry. The left side of this graphical abstract shows the velocity and vorticity (middle figure in bottom row, and right figure of the right bottom figure) profiles at peak exhalation. The simulations were validated against experiments on a 3D-print of the geometry (shown in the top figures on the right hand side). The pressure drop (right bottom corner) shows a good agreement between experiments and simulations.