The effect of viscoelasticity in an airway closure model

Abstrac

Flow in a slowly-tapering channel with oscillating walls

The flow of a fluid in a channel with walls inclined at an angle to each other is investigated at arbitrary Reynolds number. The flow is driven by an oscillatory motion of the wall incorporating a time-periodic displacement perpendicular to the channel centreline. The gap between the walls varies linearly with distance along the channel and is a prescribed periodic function of time. An approximate solution is constructed assuming that the angle of inclination of the walls is small. At leading order the flow corresponds to that in a channel with parallel, vertically oscillating walls examined by Hall and Papageorgiou \cite{HP}. A careful study of the governing partial differential system for the first order approximation controlling the tapering flow due to the wall inclination is conducted. It is found that as the Reynolds number is increased from zero the tapering flow loses symmetry and undergoes exponential growth in time. The loss of symmetry occurs at a lower Reynolds number than the symmetry-breaking for the parallel-wall flow. A window of asymmetric, time-periodic solutions is found at higher Reynolds number, and these are reached via a quasiperiodic transient from a given set of initial conditions. Beyond this window stability is again lost to exponentially growing solutions as the Reynolds number is increased

Capillary instability of a two-layer annular film: an airway closure model

Capillary instability of a two-layer liquid film lining a rigid tube is studied computationally as a model for liquid plug formation and closure of human airways. The two-layer liquid consists of a serous layer, also called the periciliary liquid layer, at the inner side and a mucus layer at the outer side. Together, they form the airway surface liquid lining the airway wall and surrounding an air core. Liquid plug formation occurs due to Plateau–Rayleigh instability when the liquid film thickness exceeds a critical value. Numerical simulations are performed for the entire closure process, including the pre- and post-coalescence phases. The mechanical stresses and their gradients on the airway wall are investigated for physiologically relevant ranges of the mucus-to-serous thickness ratio, the viscosity ratio, and the air–mucus and serous–mucus surface tensions encompassing healthy and pathological conditions of a typical adult human lung. The growth rate of the two-layer model is found to be higher in comparison with a one-layer equivalent configuration. This leads to a much sooner closure in the two-layer model than that in the corresponding one-layer model. Moreover, it is found that the serous layer generally provides an effective protection to the pulmonary epithelium against high shear stress excursions and their gradients. A linear stability analysis is also performed, and the results are found to be in good qualitative agreement with the simulations. Finally, a secondary coalescence that may occur during the post-closure phase is investigated

A new index for characterizing micro-bead motion in a flow induced by ciliary beating: Part I, experimental analysis.

Mucociliary clearance is one of the major lines of defense of the respiratory system. The mucus layer coating the pulmonary airways is moved along and out of the lung by the activity of motile cilia, thus expelling the particles trapped in it. Here we compare ex vivo measurements of a Newtonian flow induced by cilia beating (using micro-beads as tracers) and a mathematical model of this fluid flow, presented in greater detail in a second companion article. Samples of nasal epithelial cells placed in water are recorded by high-speed video-microscopy and ciliary beat pattern is inferred. Automatic tracking of micro-beads, used as markers of the flow generated by cilia motion, enables us also to assess the velocity profile as a function of the distance above the cilia. This profile is shown to be essentially parabolic. The obtained experimental data are used to feed a 2D mathematical and numerical model of the coupling between cilia, fluid, and micro-bead motion. From the model and the experimental measurements, the shear stress exerted by the cilia is deduced. Finally, this shear stress, which can easily be measured in the clinical setting, is proposed as a new index for characterizing the efficiency of ciliary beating