Computational simulation of temperature and velocity distribution in human upper respiratory airway during inhalation of hot air

The present study provides an accurate simulation of velocity and temperature distributions of inhalation thermal injury in a human upper airway, including vestibule, nasal cavity, paranasal sinuses, nasopharynx, oropharynx, larynx, and upper part of main bronchus. To this end, a series of CT scan images, taken from an adult woman, was used to construct a three dimensional model. The airway walls temperature was adjusted according to existing in vivo temperature measurements. Also, in order to cover all breathing activities, five different breathing flow rates (10, 15, 20, 30, and 40 l/min) and different ambient air temperatures (100, 200, 300, 400, and 500 °C) were studied. Different flow regimes, including laminar, transitional, and turbulence were considered and the simulations were validated using reliable experimental data. The results show that nostrils, vestibule, and nasal cavity are damaged more than other part of airway. Finally, In order to obtain the heat flux through the walls, correlations for Nusselt number for each individual parts of airway (vestibule, main upper airway, nasopharynx etc.,) are proposed.