Effects of kinematic hardening of mucus polymers in an airway closure model

Abstract

The formation of a liquid plug inside a human airway, known as airway closure, is computationally studied by considering the elastoviscoplastic (EVP) properties of the pulmonary mucus covering the airway walls for a range of liquid film thicknesses and Laplace numbers. The airway is modeled as a rigid tube lined with a single layer of an EVP liquid. The Saramito–Herschel–Bulkley (Saramito-HB) model is coupled with an Isotropic Kinematic Hardening model (Saramito-HB-IKH) to allow energy dissipation at low strain rates. The rheological model is fitted to the experimental data under healthy and cystic fibrosis (CF) conditions. Yielded/unyielded regions and stresses on the airway wall are examined throughout the closure process. Yielding is found to begin near the closure in the Saramito-HB model, whereas it occurs noticeably earlier in the Saramito-HB-IKH model. The kinematic hardening is seen to have a notable effect on the closure time, especially for the CF case, with the effect being more pronounced at low Laplace numbers and initial film thicknesses. Finally, standalone effects of rheological properties on wall stresses are examined considering their physiological values as baseline.