Numerical investigation of mucociliary clearance using power law and thixotropic mucus layers under discrete and continuous cilia motion

Mucus layer movement inside the airway system is an important phenomenon as the first defensive mechanism against pathogens. This research deals with the mucus velocity variations inside the nasal cavity using two different power law and thixotropic mucus layers. The cilia movement is replaced with four cyclic velocity profiles at the lower boundary of the mucus layer, while the upper boundary is exposed to the free-slip condition. The effects of boundary conditions and different fluid parameters are evaluated on the mucus flow. Furthermore, the replacement of power law and thixotropic mucus layers with a high viscous Newtonian mucus is examined under the free-slip condition at the mucus upper boundary. The adaptation rate is used as the criteria for replacing fluids instead of each other. The results show the mucus flow has enough time to adjust the changes from the lower boundary and the recovery stroke does not affect the mucus velocity in the effective stroke. Moreover, it is observed that the mucus flow variations are the same under the influence of recovery, breakdown, and breakdown exponent parameters. However, the effects of the exponent parameter on the mucus flow are more than the other two parameters in the recovery stroke. It is concluded that the assumption for replacing the power law mucus with a high viscous Newtonian one is acceptable. However, this assumption leads to the maximum error of 98.5% for thixotropic mucus in the recovery stroke.