Coupled discrete phase model and Eulerian wall film model for numerical simulation of respiratory droplet generation during coughing

Computational fluid dynamics is widely used to simulate droplet-spreading behavior due to respiratory events. However, droplet generation inside the body, such as the number, mass, and particle size distribution, has not been quantitatively analyzed. The aim of this study was to identify quantitative characteristics of droplet generation during coughing. Airflow simulations were performed by coupling the discrete phase model and Eulerian wall film model to reproduce shear-induced stripping of airway mucosa. An ideal airway model with symmetric bifurcations was constructed, and the wall domain was covered by a mucous liquid film. The results of the transient airflow simulation indicated that the droplets had a wide particle size distribution of 0.1–400 µm, and smaller droplets were generated in larger numbers. In addition, the total mass and number of droplets generated increased with an increasing airflow. The total mass of the droplets also increased with an increasing mucous viscosity, and the largest number and size of droplets were obtained at a viscosity of 8 mPa s. The simulation methods used in this study can be used to quantify the particle size distribution and maximum particle diameter under various conditions.