The dynamics of a thin film of Newtonian fluid coating the inner surface of an elastic circular tube is analysed. This problem is motivated by an interest in the closure of small airways of the lungs either by formation of a liquid bridge, the collapse of the airway wall or a combination of both processes. Liquid bridge formation is due to the destabilization of the liquid film that coats the inner surface of airways, while wall collapse can be due to either the high surface tension of the air–liquid interface or the flexibility of the wall.Nonlinear evolution equations for the film thickness and wall position are derived using lubrication theory, but an accurate representation of the curvatures of both the liquid and wall interfaces is employed which is valid for thick films. These approximations allow closure to be predicted. In addition, these approximations are justified by comparison with rigid-wall results obtained by solving the full Navier–Stokes equations and because fluid inertia only becomes important in the very late stages of closure. The linear stability of these equations is examined using normal-mode analysis for infinitesimal disturbances and the nonlinear stability is investigated by solving the governing equations numerically using the method of lines. Solutions show that there is a critical film thickness, strongly dependent on fluid and wall properties, above which unstable waves grow to form liquid bridges. The critical film thickness decreases with increasing surface tension or wall compliance since waves grow faster. Even for relatively stiff airways, the volume of fluid in the liquid lining required for closure can be approximately 70% of the volume for the rigid-tube case. Wall damping is an important effect only when the airway is sufficiently compliant. Airway closure occurs more rapidly with increasing unperturbed film thickness, surface tension and wall flexibility and decreasing wall damping.

We use lubrication theory and matched asymptotic expansions to model the quasisteady propagation of a liquid plug or bolus through an elastic tube. In the limit of small capillary number, asymptotic expressions are found for the pressure drop across the bolus and the thickness of the liquid film left behind, as functions of the capillary number, the thickness of the liquid lining ahead of the bolus and the elastic characteristics of the tube wall. These results generalize the well-known theory for the low capillary number motion of a bubble through a rigid tube (Bretherton 1961). As in that theory, both the pressure drop across the bolus and the thickness of the film it leaves behind vary like the two-thirds power of the capillary number. In our generalized theory, the coefficients in the power laws depend on the elastic properties of the tube.For a given thickness of the liquid lining ahead of the bolus, we identify a critical imposed pressure drop above which the bolus will eventually rupture, and hence the tube will reopen. We find that generically a tube with smaller hoop tension or smaller longitudinal tension is easier to reopen. This flow regime is fundamental to reopening of pulmonary airways, which may become plugged through disease or by instilled/aspirated fluids.

Many medical therapies require liquid plugs to be instilled into and delivered throughout the pulmonary airways. Improving these treatments requires a better understanding of how liquid distributes throughout these airways. In this study, gravitational and surface mechanisms determining the distribution of instilled liquids are examined experimentally using a bench-top model of a symmetrically bifurcating airway. A liquid plug was instilled into the parent tube and driven through the bifurcation by a syringe pump. The effect of gravity was adjusted by changing the roll angle (phi) and pitch angle (gamma) of the bifurcation (phi = gamma =0 deg was isogravitational). Phi determines the relative gravitational orientation of the two daughter tubes: when phi not equal to 0 deg, one daughter tube was lower (gravitationally favored) compared to the other. Gamma determines the component of gravity acting along the axial direction of the parent tube: when gamma not equal to 0 deg, a nonzero component of gravity acts along the axial direction of the parent tube. A splitting ratio Rs, is defined as the ratio of the liquid volume in the upper daughter to the lower just after plug splitting. We measured the splitting ratio, Rs, as a function of: the parent-tube capillary number (Cap); the Bond number (Bo); phi; gamma; and the presence of pre-existing plugs initially blocking either daughter tube. A critical capillary number (Cac) was found to exist below which no liquid entered the upper daughter (Rs = 0), and above which Rs increased and leveled off with Cap. Cac increased while Rs decreased with increasing phi, gamma, and Bo for blocked and unblocked cases at a given Cap > Ca,. Compared to the nonblockage cases, Rs decreased (increased) at a given Cap while Cac increased (decreased) with an upper (lower) liquid blockage. More liquid entered the unblocked daughter with a blockage in one daughter tube, and this effect was larger with larger gravity effect. A simple theoretical model that predicts Rs and Cac is in qualitative agreement with the experiments over a wide range of parameters.

Liquid plugs may form in pulmonary airways during the process of liquid instillation or removal in many clinical treatments. During inspiration the plug may split at airway bifurcations and lead to a nonuniform final liquid distribution, which can adversely affect treatment outcomes. In this paper, a combination of bench top experimental and theoretical studies is presented to study the effects of inertia and gravity on plug splitting in an airway bifurcation model to simulate the liquid distributions in large airways. The splitting ratio, Rs, is defined as the ratio of the plug volume entering the upper (gravitationally opposed) daughter tube to the lower (gravitationally favored) one. Rs is measured as a function of parent tube Reynolds number, Rep; gravitational orientations for roll angle, phi, and pitch angle, gamma; parent plug length LP; and the presence of pre-existing plug blockages in downstream daughter tubes. Results show that increasing Rep causes more homogeneous splitting. A critical Reynolds number Rec is found to exist so that when Rep < or = Rec, Rs = 0, i.e., no liquid enters the upper daughter tube. Rec increases while Rs decreases with increasing the gravitational effect, i.e., increasing phi and gamma. When a blockage exists in the lower daughter, Rec is only found at phi = 60 deg in the range of Rep studied, and the resulting total mass ratio can be as high as 6, which also asymptotes to a finite value for different phi as Rep increases. Inertia is further demonstrated to cause more homogeneous plug splitting from a comparison study of Rs versus Cap (another characteristic speed) for three liquids: water, glycerin, and LB-400X. A theoretical model based on entrance flow for the plug in the daughters is developed and predicts Rs versus Rep. The frictional pressure drop, as a part of the total pressure drop, is estimated by two fitting parameters and shows a linear relationship with Rep. The theory provides a good prediction on liquid plug splitting and well simulates the liquid distributions in the large airways of human lungs.

Liquid plugs may form in pulmonary airways during the process of liquid instillation or removal in many clinical treatments. Studies have shown that the effectiveness of these treatments may depend on how liquids distribute in the lung. Better understanding of the fundamental fluid mechanics of liquid plug transport will facilitate treatment strategies. In this paper, we develop a numerical model of steady plug propagation driven by gravity and pressure in a two-dimensional liquid-lined channel oriented at an angle α with respect to gravity. We investigate the effects of gravity through the Bond number, Bo, and α; the plug propagation speed through the capillary number, Ca, or the Reynolds number, Re; the plug length LP, and the surfactant concentration C0. Without gravity, i.e., Bo=0, the plug is symmetric, and there are two regimes for the flow: two wall layers and two trapped vortices in the core. There is no flow interaction between the upper and lower half plug domains. When Bo≠0 and α≠0, π, fluid is found to flow from the upper precursor film, through the core and into the lower trailing film. Then the number of vortices can be zero, one, or two, depending on the flow parameters. The vortices have stagnation points on the interface when C0=0, however when the surfactant is present (C0&gt;0), the vortices detach from the interface and create saddle points inside the core. The front meniscus develops a capillary surface wave extending into the precursor film. This is where the film is thinnest and thus the wall shear stress is highest, as high as ∼100dyn∕cm2 in adult airways, which indicates a significant risk of pulmonary airway epithelial cell damage. Adding surfactant can decrease the peak magnitude of the shear stress, thus reducing the risk of cell damage. The prebifurcation asymmetry of the plug is quantified by the volume ratio, Vr, defined as the ratio of the liquid above to that below the center line of the channel. Vr is found to increase with LP, Ca, Re, and C0, while it decreases with Bo. The total mass left behind in the trailing films increases with Bo for any α at α&gt;2π∕5, Ca and α for any value of Bo&gt;0.